Mecánica cuántica/Matriz densidad y evolución temporal

Ecuación de Schrödinger: evolución temporal

¿Cuál es la dinámica? ¿Cómo cambia mi sistema con el tiempo?

Postulado 5: La evolución temporal del vector $|\alpha (t)⟩$ viene dada por la ecuación de Schrödinger

donde $H$ es el operador Hamiltoniano asociado a la energía del sistema (es hermítico, $H=H^{†}$).

• La ecuación de Schrödinger es determinista. Si se conoce el estado en un instante $|\alpha (t_0)⟩$, el estado queda completamente determinado en cualquier instante posterior $|\alpha (t)⟩$ (la indeterminación cuántica aparece al medir, no en la dinámica).
• La ecuación de Schrödinger para los bras es

$-i\hslash \frac{d}{dt}⟨\alpha (t)|=⟨\alpha (t)|H,$

por lo que, considerando que el operador A puede depender del tiempo,

$\frac{\partial }{\partial t}A\ne 0,$

obtendremos una bellísima ecuación

$\begin{array}{ccc}i\hslash \frac{d}{dt}⟨A{⟩}_{\alpha }& =& i\hslash \frac{d}{dt}[⟨\alpha |A|\alpha ⟩]\\ & =& i\hslash \frac{d}{dt}[⟨\alpha |]A|\alpha ⟩+i\hslash ⟨\alpha |A\frac{d}{dt}[|\alpha ⟩]+i\hslash ⟨\alpha |\frac{d}{dlt}\left[A\right]|\alpha ⟩\\ & =& i\hslash [\frac{-1}{i\hslash }⟨\alpha |H]A|\alpha ⟩+i\hslash ⟨\alpha |A[\frac{1}{i\hslash }H|\alpha ⟩]+i\hslash {⟨\frac{dA}{dt}⟩}_{\alpha }\\ & =& -⟨\alpha |HA|\alpha ⟩+⟨\alpha |AH|\alpha ⟩+i\hslash {⟨\frac{dA}{dt}⟩}_{\alpha }\\ & =& ⟨\alpha |\underset{[A,H]}{\underbrace{(AH-HA)}}|\alpha ⟩+i\hslash {⟨\frac{dA}{dt}⟩}_{\alpha }.\end{array}$

En el último paso solo hemos usado la definición de suma (o resta) de operadores. Tenemos por tanto que

Si el operador A no depende explícitamente del tiempo $A\ne A(t)$,

$\frac{d}{dt}⟨A{⟩}_{\alpha }=\frac{1}{i\hslash }⟨[A,H]{⟩}_{\alpha }.$

Es decir, la derivada temporal del valor esperado del operador A será, salvo constante, el conmutador $[A,H]$. Si A conmuta con H, $[A,H]=0$, entonces, necesariamente

$⟨A{⟩}_{\alpha }=\text{cte}$.

Esto es formalmente idéntico a lo que ocurre en Mecánica Clásica, donde la derivada temporal de cualquier función dinámica A (los observables en dicha formulación, donde no hay valores esperados ya que su valor está siempre deterministamente concretado) viene dada por los corchetes de Poisson de A con la función Hamiltoniana

$\dot{A}=\frac{d}{dt}A=\{A,H\}+\frac{\partial A}{\partial t}\text{(en Mecánica Clásica)}.$

• $i\hslash \frac{d\rho }{dt}=[H,\rho (t)]$
• $⟨H{⟩}_{\alpha }=\text{cte}\Rightarrow (\frac{\partial H}{\partial t}=0\iff \text{Sistema conservativo})$.

El operador evolución temporal

¿Cómo evolucionan los estados propios de $H$?

$H|E⟩=E|E⟩$

¿Cómo es la ecuación de Schödinger para un estado propio?

$i\hslash \frac{d}{dt}|E⟩=H|E⟩=E|E⟩$

$\Rightarrow |E(t)⟩=|E(0)⟩\underset{fasetemporal}{\underbrace{exp(-\frac{i}{\hslash }E(t-t_0))}}$

• La energía se repite.
• La ecuación de Schrödinguer es lineal: sus soluciones pueden componerse linealmente (si $|{\alpha }_1⟩$ y $|{\alpha }_2⟩$ son soluciones $\Rightarrow$ $a_1|{\alpha }_1⟩+a_2|{\alpha }_2⟩$ también es solución)
• La evolución entre un sistema inicial y uno final $t$ se puede describir mediante un operador lineal, el operador de evolución temporal $U(t,t_0)$ tal que

$|\alpha (t)⟩=U(t,t_0)|\alpha (t_0)⟩$

$U(t,t_0)$ tiene las siguientes propiedades:

• $U(t_0,t_0)=ℐ\text{(operador identidad)}$

Podemos obtener la forma explícita del operador sustituyendo $|\alpha (t)⟩=U(t,t_0)|\alpha (t_0)⟩$ en la ecuación de Schrödinguer

$i\hslash \frac{d}{dt}U(t,t_0)|\alpha (t_0)=HU(t,t_0)|\alpha (t_0)⟩$

de donde obtenemos la igualdad de operadores

$\frac{d}{dt}U(t,t_0)=\frac{1}{i\hslash }HU(t,t_0)=-\frac{i}{\hslash }HU(t,t_0).$

Si $H$ no depende de $t$ (sistema conservativo), al igual que con las ecuaciones diferenciales de funciones, se tiene que

• $U(t,t_0)=e^{-\frac{i}{\hslash }H(t-t_0)}$

En general, el desarrollo en serie de Taylor de $e^A$ es $e^A=ℐ+A+\frac{1}{2!}{A}^2+\cdots$ pero, en la base en la que la representación matricial del operador $A$ es diagonal

$A=\left(\begin{array}{cc}{a}_{11}& 0\\ 0& a_{22}\end{array}\right),A^n\equiv AA\cdots A=\left(\begin{array}{cc}{a}_{11}^n& 0\\ 0& a_{22}^n\end{array}\right)$

se tiene que

$e^A=\left(\begin{array}{cc}1& 0\\ 0& 1\end{array}\right)+\left(\begin{array}{cc}{a}_{11}& 0\\ 0& a_{22}\end{array}\right)+\left(\begin{array}{cc}\frac{a_{11}}{2!}& 0\\ 0& \frac{a_{11}}{2!}\end{array}\right)+\cdots =\left(\begin{array}{cc}{e}^{a_{11}}& 0\\ 0& e^{a_{22}}\end{array}\right).$

De manera más más compacta

${\left(e^A\right)}_{ij}=e^{A_{ij}}=\left(\begin{array}{cc}{e}^{a_{11}}& 0\\ 0& e^{a_{22}}\end{array}\right).$

• $U^{-1}(t,t^{\prime })=U(t^{\prime },t)$ ya que

$U(t,t^{\prime })U(t^{\prime },t)|\alpha (t)⟩=|\alpha (t)⟩$

Ahora, veremos como a primer orden de su desarrollo de Taylor, el operador evolución temporal es unitario ($U=U^{-1}$)

$U(t,t_0)=e^{-\frac{i}{\hslash }H(t-t_0)}\Rightarrow U(t+dt,t)\approx ℐ-\frac{i}{\hslash }Hdt$

$U^{-1}(t,t_0)=e^{\frac{i}{\hslash }H(t-t_0)}\Rightarrow U^{-1}(t+dt,t)\approx ℐ+\frac{i}{\hslash }Hdt=ℐ+\frac{i}{\hslash }{H}^{†}dt$

$U{U}^{-1}=ℐ+𝒪\left(d{t}^2\right).$

Seguimos tras esta última nota

$U(t+dt,t)={ℐ}^{†}+\frac{i}{\hslash }{H}^{-1}dt=ℐ+\frac{i}{\hslash }Hdt$

$U(t+dt,t)=U^{-1}(t+dt,t)$

$\Rightarrow U^{†}=U^{-1}(esunitario).$

• Cualquier tranformación de simetría viene descrita por un operador unitario

$U|\alpha ⟩\to ⟨\alpha |\underset{ℐ}{\underbrace{U{U}^{†}}}|\alpha ⟩.$

Estados estacionarios y constantes del movimiento

Un estado conservativo es aquel que no depende del tiempo $H\ne H(t).$ En un estado conservativo los estados propios de $H$ se denominan estados estacionarios.

$H|\alpha (t_0)⟩=E|\alpha (t_0)⟩$

$|\alpha (t)⟩=e^{-\frac{i}{\hslash }E(t-t_0)}|\alpha (t_0)⟩$

$H|\alpha (t)⟩=E|\alpha (t)⟩.$

Si un observable $A$ conmuta con $H$ (se sobreentiende que nos referimos al operador hermítico asociado al observable $A$) y como habitualmente el observable no depende explícitamente del tiempo,

$\frac{d}{dt}⟨A{⟩}_{\alpha }=0\Rightarrow ⟨A{⟩}_{\alpha }=\text{cte}.$

$A$ se denomina entonces una constante del movimiento.

Ejemplo

Supongamos que nuestro sistema en un tiempo t=0 se encuentra en el estado

$ℍ\ni |\alpha ⟩=\frac{1}{\sqrt{2}}(|+⟩+i|-⟩)$

y que el operador Hamiltoniano se encuentra definido por la siguiente combinación lineal de operadores proyección

$H=\hslash \omega (|+⟩⟨+|-|-⟩⟨+|).$

¿Cuál es $|\alpha (t)⟩$?

Usando la base $\{|+⟩,|-⟩\}$, estos dos vectores son obviamente,

$|+⟩=\left(\begin{array}{c}1\\ 0\end{array}\right)$

$|-⟩=\left(\begin{array}{c}0\\ 1\end{array}\right)$

Expresado en dicha base, el Hamiltoniano toma la forma matricial

$\begin{array}{cc}H& \doteq \hslash \omega [\left(\begin{array}{c}1\\ 0\end{array}\right)\left(\begin{array}{cc}1& 0\end{array}\right)-\left(\begin{array}{c}0\\ 1\end{array}\right)\left(\begin{array}{cc}0& 1\end{array}\right)]\\ & =\hslash \omega [\left(\begin{array}{cc}1& 0\\ 0& 0\end{array}\right)-\left(\begin{array}{cc}0& 0\\ 0& 1\end{array}\right)]\\ & =\hslash \omega \left(\begin{array}{cc}1& 0\\ 0& -1\end{array}\right)\end{array}$

que es diagonal. Como vimos anteriormente, ya que $H\ne H(t)$

$U(t,0)=e^{-\frac{i}{\hslash }Ht}$

y como demostramos, al ser la matriz que representa $H$ diagonal en dicha base, la exponencial de $H$ también es diagonal y tiene la forma

$U(t,0)=e^{-\frac{i}{\hslash }Ht}\doteq \left(\begin{array}{cc}{e}^{-i\omega t}& 0\\ 0& e^{+i\omega t}\end{array}\right)$

$\begin{array}{cc}|\alpha (t)⟩& =U(t,0)|\alpha (t_0)⟩\\ & \doteq \left(\begin{array}{cc}{e}^{-i\omega t}& 0\\ 0& e^{+i\omega t}\end{array}\right)\frac{1}{\sqrt{2}}(\left(\begin{array}{c}1\\ 0\end{array}\right)+i\left(\begin{array}{c}0\\ 1\end{array}\right))\\ & =\frac{1}{\sqrt{2}}\left(\begin{array}{cc}{e}^{-i\omega t}& 0\\ 0& e^{+i\omega t}\end{array}\right)\left(\begin{array}{c}1\\ i\end{array}\right)\\ & =\frac{1}{\sqrt{2}}\left(\begin{array}{c}{e}^{-i\omega t}\\ i{e}^{+i\omega t}\end{array}\right).\end{array}$

Reglas de cuantización

¿Cómo construir los operadores que representan a los observables físicos?

El sistema físico más simple es una partícula escalar sometida a un potencial central.

Los observable físicos más evidentes son la posición $\overrightarrow{r}=(x,y,z)$ y su momento $\overrightarrow{p}=(p_x,p_y,p_z)$.

Postulado 4: Si en un sistema físico las coordenadas son $x_1,\dots ,x_n$ y sus momentos conjugados son $p_1,\dots p_n$, entonces los operadores $X_i$ y $P_i$ que representan a dichos observables deben satisfacer las reglas de conmutación

$[X_i,X_j]=0;[P_i,P_j]=0;[X_i,P_j]=i\hslash {\delta }_{ij}ℐ$

¿Qué forma tiene un operador arbitrario?

Si el sistema tiene un observable cuya expresión clásica es $A(x_i,p_i,t)$ entonces el operador correspondiente se obtiene sustituyendo de manera adecuada las variables $x_i$ y $p_i$ por los operadores $X_i$ y $P_i$ asociados.

El operador Hamiltoniano

¿Qué es el (operador) Hamiltoniano clásico? Realmente una función, la función Hamiltoniana.

Formulación Lagrangiana

Para cualquier sistema físico se puede definir una función que depende de las posiciones y de las velocidades tal que toda la física del sistema dependa de dicha función. Dicha función es la llamada función Lagrangiana y se define como $L=T-V$.

En cada estado la lagrangiana toma un valor.

Se define la acción como la siguiente integral de línea entre $x_1$ y $x_2$

$S=\underset{x_1}{\overset{x_2}{\int }}dtL(x_i,{\dot{x}}_i).$

La naturaleza elige el camino con acción más pequeña (principio de mínima acción).

$\Rightarrow$ Ecuaciones de Euler-Lagrange

$\frac{d}{dt}\left(\frac{\partial L}{\partial {\dot{x}}_i}\right)-\frac{\partial L}{\partial x_i}=0$

Es equivalente a la formulación Newtoniana y no menos extraña.

Se define el momento conjugado a la variable $x_i$ como

$p_i\equiv \frac{\partial L}{\partial {\dot{x}}_j}$

Formulación de Hamilton

$H=H(x_i,p_i,t)={\dot{x}}_i{p}_i-L(x_i,{\dot{x}}_i,t)$

Principio de mínima acción $\Rightarrow$

${\dot{x}}_i=\frac{\partial H}{\partial p_i}$

${\dot{p}}_i=-\frac{\partial H}{\partial x_i}$

Normalmente se cumple que $H=E_T=T+V=\frac{p^2}{2m}+V(\overrightarrow{r})$

Cuánticamente

$H=\frac{P^2}{2m}+V(X_i)$ donde $P=(P_1,P_2,P_3)$ es un operador vectorial (3 operadores).

Se tiene que

$P^2=\sum _i{P}_i^2.$

Los paréntesis de Poisson se definen como

$\{A(x_i,p_j),B(x_i,p_j)\}\stackrel{\mathrm{d}\mathrm{e}\mathrm{f}}{=}\sum [\frac{\partial A}{\partial x_i}\frac{\partial B}{\partial p_i}-\frac{\partial A}{\partial p_i}\frac{\partial B}{\partial x_i}]$

$\{x,p_x\}=1$

$[X,P]=i\hslash ℐ$

A los franceses les gusta mucho la asociación

$\{,\}\to \frac{[,]}{i\hslash }.$

Reglas de superselección

Supongamos que existe un observable $Q$ que conmuta con todos los observables $A_i$

$[Q,A_i]=0$

$\Rightarrow$ Si simultáneamente tengo un estado propio de $Q$ tras cualquier medida de otro observable, el estado continuará siendo propio de $Q$ con el mismo valor propio.

Demostración (suponiendo que no son degenerados)

$|\alpha ⟩\text{ propio de Q}$

$Q|\alpha ⟩=q|\alpha ⟩$

se mide $A$,

$|\alpha ⟩\stackrel{a}{\to }|a⟩$

con probabilidad $|{\alpha }_a{|}^2$

$Q|a⟩=Q\frac{P_a|\alpha ⟩}{|P_a|\alpha ⟩|}$

en la base de autoestados de $A$, $Q$ es diagonal y $P_a$ también

$\frac{P_a}{|P_a|}=Q|\alpha ⟩=q|a⟩(?)\bullet$

Se dice que $Q$ es un operador de super-selección, y los estados "físicos" tendrán siempre un valor de $Q$ bien definido.

Es algo que no me permite mezcla cuántica, combinación lineal, con vectores con distintos valores de q. (Viven en espacios de Hilbert distintos (?)).

Imagen de Heisenberg

Lo que hemos visto hasta ahora se llama imagen de Schrödinger.

$⟨A{⟩}_{\alpha }=⟨\alpha |A|\alpha ⟩=⟨{\alpha }_0|U^{†}A|{\alpha }_0⟩$

(Uno construye una metodología distinta pero de forma que los autovalores evolucionen de la misma forma (?)).