Ejemplo 6.2: Evolución reversible en fibra óptica¶
$\newcommand{\ket}[1]{|{#1}\rangle}\newcommand{\bra}[1]{\langle{#1}|}$Considere el ejemplo anterior de la transmisión de un fotón por una línea óptica:
Ésta se puede modelar como $U_{\alpha}$ donde $\alpha$ se corresponde con el ángulo de rotación de polarización: \begin{align*} U_{\alpha} = \left[\begin{matrix} \cos(\alpha)& -\sin(\alpha)\\\sin(\alpha)& \cos(\alpha) \end{matrix}\right]. \end{align*}
In [1]:
function [U] = fibra(alpha) % Transformación unitaria de una fibra optica
U = [[cos(alpha), -sin(alpha)]; [sin(alpha), cos(alpha)]];
end
- Si el ángulo $\alpha$ es conocido, basta aplicar una rotación con un ángulo $-\alpha$ para recuperar el fotón original. De hecho, en este caso $U_{\alpha}^H = U_{\alpha}^{-1} = U_{-\alpha}$: \begin{align*} U_{\alpha} U_{-\alpha} &= \left[\begin{matrix} \cos(\alpha)& -\sin(\alpha)\\\sin(\alpha)& \cos(\alpha) \end{matrix}\right] \left[\begin{matrix} \cos(\alpha)& \sin(\alpha)\\-\sin(\alpha)& \cos(\alpha) \end{matrix}\right]\\ &= \left[\begin{matrix} \cos(\alpha)^2 + \sin(\alpha)^2 & \cos(\alpha)\sin(\alpha)-\sin(\alpha)\cos(\alpha)\\ \sin(\alpha)\cos(\alpha)-\cos(\alpha)\sin(\alpha)& \sin(\alpha)^2 + \cos(\alpha)^2 \end{matrix}\right]= \boldsymbol{I}, \end{align*} donde hemos utilizado que $\cos(-\alpha)=\cos(\alpha)$ y que $\sin(-\alpha) = -\sin(\alpha)$ al ser funciones par e impar, respectivamente; y la conocida identidad trigonométrica $\cos(\alpha)^2 + \sin(\alpha)^2=1$.
In [3]:
fprintf('\nComprobamos que U_{-alpha} es de hecho la inversa de U_{alpha}:\n\n')
alpha = 0.2;
U = fibra(alpha)
invU = fibra(-alpha)
U*invU
invU*U
Comprobamos que U_{-alpha} es de hecho la inversa de U_{alpha}:
U =
0.9801 -0.1987
0.1987 0.9801
invU =
0.9801 0.1987
-0.1987 0.9801
ans =
1.0000e+00 -1.2308e-17
-1.2308e-17 1.0000e+00
ans =
1.0000e+00 1.2308e-17
1.2308e-17 1.0000e+00
- Si el ángulo $\alpha$ es desconocido, en cambio, no es posible recuperar el fotón original al no poder construir la transformación inversa. En este caso sería necesario estimar de alguna forma la rotación que introduce la fibra para así poder compensarla.
In [5]:
function [rho] = foton(theta) % Matriz de densidad de probabilidad de un foton polarizado
ket_theta = [cos(theta); sin(theta)];
rho = ket_theta * ket_theta';
end
fprintf('\nFibra que introduce una rotación aleatoria:\n\n')
rho = foton(0)
U = fibra(2*pi*rand)
sigma = U*rho*U'
Fibra que introduce una rotación aleatoria: rho = 1 0 0 0 U = -0.3834 -0.9236 0.9236 -0.3834 sigma = 0.1470 -0.3541 -0.3541 0.8530