Sistemas abiertos: Derivación de la ecuación maestra para un

sistema de dos niveles en interacción con un baño térmico

Ramsés Sánchez
2 de Febrero de 2011

En mecánica cuántica la evolución de los sistemas aislados está completamente determinada por
transformaciones unitárias. En la práctica, sin embargo, la consideración de sistemas totalmente aislados
es irreal debido a la inevitable existencia de interacción entre el sistema bajo estudio y el medio ambiente
en que se encuentra inmerso. Esta interacción con el medio genera correlaciones no locales entre ambos
subsistemas, ocasiona que la evolución del subsistema al cual se tiene acceso sea no unitaria y produce
la supresión local de la capacidad de este subsistema para mostrar interferencia. Este último fenómeno
es conocido como decoherencia y es uno de los argumentos utilizados para explicar, de cierta forma,
el origen de las propiedades clásicas de los sistemas a partir del formalismo de la mecánica cuántica.
Algunas referencias significativas en este respecto son las siguientes: [3],[6],[4],[7].
En este trabajo vamos a estudiar el formalismo de la ecuación maestra, el cual nos permite obtener
una ecuación diferencial que describe la dinámica del sistema bajo estudio en interacción con un baño
térmico, de manera que las propiedades del baño aparecen como parámetros en la ecuación diferencial.
Esta ecuación se obtiene al realizar un conjunto de aproximaciones sobre la interacción entre el sistema
y el baño térmico.
Presentamos a continuación la derivación de la ecuación maestra para un sistema de dos niveles
inmerso en un baño térmico descrito por un conjunto de osciladores armónicos acoplados al sistema.

1 Ecuación maestra de un sistema de dos niveles inmerso en un

baño térmico
Consideremos un sistema de dos niveles descrito por la matriz de densidad ρs . Los niveles de este
sistema son etiquetados como |1i = (0 1)T y |2i = (1 0)T , cuyas energı́as respectivas E1 y E2 satisfacen
E2 > E1 . El Hamiltoniano de este sistema puede ser escrito a través de la descomposición espectral
como:

Hs = E1 |1ih1| + E2 |2ih2| (1)

Es conveniente reescribir este operador en términos de la matriz de Pauli σ3 , la cual puede ser expresar
en términos de los {|1i, |2i} como σ3 = |2ih2| − |1ih1|:
1 1
Hs = (E1 + E2 )I + (E2 − E1 )σ3 (2)
2 2
El primer término de la ecuación (2) es una constante que podemos despreciar debido a que ex-
perimentalmente medimos diferencias de energı́a, y por tanto no importa donde tomemos el cero. El
Hamiltoniano del sistema se escribe entonces como:

h̄ω
Hs = 2 σ3 (3)

con ω = (E2 −E
h̄
1)
. Dicho sistema interactúa con un baño térmico descrito por la matriz de densidad ρb
y formado por un conjunto de osciladores armónicos acoplados al sistema de dos niveles. El Hamiltoniano
del baño tiene la forma:

1
 X
Hb = h̄ωi b†i bi ; [bi , b†j ] = δij (4)
i

La interacción entre ambos subsistemas es descrita, utilizando la aproximación dipolar, a través del
siguiente Hamiltoniano:
Hi = −DE (5)
Donde D es el operador de momento dipolar definido como D = eR, R es el operador de posición
y e la carga del electrón. E es el operador asociado al campo eléctrico del baño. Esta aproximación
consiste básicamente en considerar la cantidad k̄. r̄ ∼ aλ0 ≪ 1, con k̄ el número de onda y a0 el radio de
Bohr, debido a que las dimensiones del vector posición r̄ están en el orden de las dimensiones at́omicas.
Esta consideración ocasiona que el operador asociado al vector potencial del baño sea igual a cero en el
calibre apropiado. Los detalles de esta aproximación pueden ser encontrados en [1].
Es conocido que en la representación ({|1i, |2i}) el operador de posición se escribe como:

R = (h1|R|2i|1ih2| + h2|R|1i|2ih1|) (6)

Sustituyendo en la expresión del operador de momento dipolar se tiene que:

D = e (h1|R|2i|1ih2| + h2|R|1i|2ih1|) (7)

Podemos tomar D real: d ≡ eh1|R|2i = eh2|R|1i y definir σ = |1ih2| y σ † = |2ih1| como los operadores
de subida y bajada, cuya representación matricial viene dada por:
   
† 0 1 0 0
σ = ; σ= (8)
0 0 1 0
Entonces:
D = d(σ + σ † ) (9)
Por otro lado podemos considerar el campo eléctrico del baño térmico como la suma de i campos
eléctricos estacionarios: X
E= αi (bi + b†i )sen(k̄.r̄) (10)
i

donde αi es el campo eléctrico por fotón. Algunos detalles al respecto pueden ser encontrados en [5]. El
Hamiltoniano de interacción más general tiene entonces la forma:
X
Hi = h̄gi (σ † + σ)(bi + b†i ) (11)
i

donde gi = −(αi d/h̄)sen(k̄.r̄). Los cuatro términos que aparecen en esta última ecuación pueden ser
interpretados como:
• σ † bi : El sistema absorbe un fotón y es excitado.
• σb†i : El sistema emite un fotón y se des-excita

• σ † b†i : El sistema emite un fotón y es excitado.

• σbi : El sistema absorbe un fotón y se des-excita.
Note que sólo los dos primeros procesos son conservativos. El Hamiltoniano que describe el sistema
más el baño viene dado entonces por:
h̄ω X X
H = H0 + H i = [ σ3 + h̄ωi b†i bi ] + h̄gi (σ † + σ)(bi + b†i ) (12)
2 i i

Queremos determinar una ecuación para ρs (t) en la que las propiedades del baño aparezcan como
parámetros. Para esto, supongamos primero que inicialmente los estados del sistema y del baño no están
correlacionados:
ρsb (0) = ρs (0) ⊗ ρb (0) (13)

2
 Donde ρsb es la matriz de densidad que describe el sistema completo. Supongamos además que el
baño es suficientemente grande y en equilibrio térmico como para verse afectado por el sistema:

ρsb (t) = ρs (t) ⊗ ρb (0) ∀t (14)

El sistema completo se mantiene aislado y evoluciona unitariamente a través de la ecuación de Liou-

ville:
dρsb i
(t) = − [H, ρsb (t)] (15)
dt h̄
Al evolucionar aparecen correlaciones no locales entre el subsistema de dos niveles y el baño térmico.
Como queremos estudiar la evolución del subsistema ρs sin tener que considerar explı́citamente los grados
de libertad del baño, tendremos que tomar una traza sobre los estados del subsistema ρb . Sin embargo,
antes de hacer esto conviene colocarse en el cuadro de interacción. Considere la siguiente transformación
unitaria:
i i
ρsb (t) 7→ ρ¯sb (t) = e h̄ H0 t ρsb (t)e− h̄ H0 t (16)
Derivando respecto al tiempo la ecuación anterior hallamos que:
i
dρ¯ − i H0 t i i i i
dt (t)
sb
= i
h̄ H0 e
h̄ H0 t ρ
sb e h̄ + e h̄ H0 t ∂ρ∂tsb e− h̄ H0 t − h̄i e h̄ H0 t ρsb H0 e− h̄ H0 t
i i
= i
h̄ [H0 , ρ¯
sb ] + e h̄ H0 t ∂ρ∂tsb e− h̄ H0 t
(17)
i
= h̄ [H0 − H̄, ρ¯sb ]

= − h̄i [H̄i , ρ¯sb ]

Necesitamos entonces Hi en el cuadro de interacción:

H̄i = U HiU †

h̄gi (σ † + σ)(bi + b†i )e−( h̄ )Ho t

i i
= e h̄ Ho t
P
i (18)

h̄gi (σ † bi + σb†i )e−( h̄ )Ho t + e h̄ Ho t h̄gi (σ † b†i + σbi )e−( h̄ )Ho t

i i i i
= e h̄ Ho t
P P
i i

Trabajemos el primer término del lado derecho de la ecuación (18). Dado que eαA Be−αA = B +
2
α[A, B] + α2! [A, [A, B]] + ..., tal que A y B son operadores con dominios bien definidos en sus espacios
correspondientes; escribimos:
i i
e h̄ Ho t h̄gi (σ † bi + σb†i )e−( h̄ )Ho t
P
i

h̄gi (σ † bi + σb†i ) + it[(ω σ23 + ωk b†k bk ), h̄gi (σ † bi + σb†i )]

P P P
= i k i (19)
2
− t2 [(ω σ23 + ωk b†k bk ), [(ω σ23 + ωj b†j bj ), h̄gi (σ † bi + σb†i )]] + ...
P P P
k j i

Veamos el primer conmutador de la ecuación (19):

[(ω σ23 + ωk b†k bk ), h̄gi (σ † bi + σb†i )]

P P
k i

σ3 †
+ [ σ23 , σb†i ]) + ωk h̄gi ([b†k bk , σ † bi ] + [b†k bk , σb†i ])
P P
= ωh̄ i gi ([ 2 , σ bi ] i,k

†
bi − σb†i ) + ωk h̄gi (−σ † bk δik + σb†k δik ) (20)
P P
= ωh̄ i gi (σ i,k

h̄gi (σ † bi − σb†i )(ω − ωi )

P
= i

h̄gi (σ † bi )(ω − ωi ) − h̄gi (σb†i )(ω − ωi )

P P
= i i

Donde se utilizó la propiedad [σ, σ † ] = −σ3 . Luego, el segundo conmutador de (19) queda como:

3
 [(ω σ23 + ωk b†k bk ), [(ω σ23 + ωj b†j bj ), h̄gi (σ † bi + σb†i )]]
P P P
k j i

= [(ω σ23 + ωk b†k bk ), h̄gi (σ † bi − σb†i )(ω − ωi )]

P P
k i

gi ([ σ23 , σ † bi ] − [ σ23 , σb†i ]) + ωk h̄gi ([b†k bk , σ † bi ]

P P
= {ωh̄ i i,k

−[b†k bk , σb†i ]}(ω − ωi )

h̄gi (σ † bi )(ω − ωi )2 + h̄gi (σb†i )(ω − ωi )2

P P
= i i

A partir
P de estos dos resultados podemos suponer que, en general, tomando para abreviar A =
(ω σ23 + k h̄ωk b†k bk ) y B = i h̄gi (σ † bi + σb†i ), se cumple que:
P

X X
[A, [A, [A, ... [A[A, B]]... ] = h̄gi (σ † bi )(ω − ωi )n + (−1)n h̄gi (σb†i )(ω − ωi )n (21)
i i

Es sencillo probar la ecuación (21) por inducción. De hecho, supongamos que la expresión anterior
es válida para n conmutadores. Para n + 1 conmutadores tendrı́amos que:

[A, [A, [A, ... [A[A, B]]... ]

= [(ω σ23 + ωk b†k bk ), h̄gi (σ † bi )(ω − ωi )n + (−1)n h̄gi (σb†i )(ω − ωi )n ]

P P P
k i i

h̄gi ω(σ † bi )(ω − ωi )n + h̄gi ωi σ † (−bi )(ω − ωi )n + (−1n ) h̄gi ωσ(−b†i )(ω − ωi )n
P P P
= i i i

+(−1n ) h̄gi ωi (σb†i )(ω − ωi )n

P
i

h̄gi (σ † bi )(ω − ωi )n+1 + (−1)n+1 h̄gi (σb†i )(ω − ωi )n+1

P P
= i i

Este resultado, junto a la expresión para sólo un conmutador (Ec.20) que corresponde al caso inicial,
prueban por inducción la ecuación (21).
Luego, el primer término de la ecuación para el Hamiltoniano de interacción en el cuadro de interacción
(Ec.18) tiene la siguiente forma:
i i
e h̄ Ho t h̄gi (σ † bi + σb†i )e−( h̄ )Ho t
P
i

(it)n (−it)n
h̄gi (σ † bi ) − ωi )n + h̄gi (σb†i ) − ωi )n
P P P P
= i n n! (ω i n n! (ω

h̄gi (σ † bi eit(ω−ωi ) + σb†i e−it(ω−ωi ) )

P
= i

Procediendo de manera análoga con el segundo término de la ecuación (18), escribimos el operador
H̄i como:

H̄i = h̄(σ † G(t) + σG† (t) + σF (t) + σ † F (t)† ) (22)

donde,

= Pi gi bi eit(ω−ωi )
P
G(t)
it(ω+ωi ) (23)
F (t) = i gi bi e

Ahora, recordemos que al escribir el Hamiltoniano de interacción consideramos la aproximación dipo-

lar; para que dicha aproximación sea válida el baño térmico debe estar cerca de la resonancia con la
transición atómica del sistema de dos niveles [2]. En otras palabras, ω ∼ ωi lo cual es lo mismo a decir
que ω − ωi << ω + ωi . Debido a esto los exponentes complejos eit(ω+ωi ) , e−it(ω+ωi ) pueden ser consid-
erados como términos que oscilan rápidamente. Luego en cualquier escala de tiempo apreciable dichas
oscilaciones tienden a cero rápidamente.

4
 Entonces podemos despreciar los términos de rápida oscilación F (t) y F (t)† y tomar el Hamiltoniano
de interacción en la cuadro de interacción como:

H̄i = h̄(σ † G(t) + σG† (t)) (24)

Esta aproximación es conocida como aproximación de onda rotante. La ecuación de evolución del
sistema más el baño se escribe entonces de la forma siguiente:
dρ̄sb
dt (t) = − h̄i [H̄i (t), ρ̄sb (t)]
(25)
= − h̄i [σ † G(t) + σG† (t), ρ̄sb (t)]
Podemos integrar la relación anterior y sustituir ρ̄sb (t) dentro de su solución formal para obtener:
t t t′
i 1
Z Z Z
ρ̄sb (t) = ρ̄sb (0) − dt′ [H̄i (t′ ), ρ̄sb (0)] − dt′ dt′′ [H̄i (t′ ), [H̄i (t′′ ), ρ̄sb (t′′ )]
h̄ 0 h̄2 0 0

Esta solución a segundo orden es conocida como aproximación de Born. Derivando este último
resultado se llega a que:
Z t
dρ̄sb i i
(t) = − [H̄i (t), ρ̄sb (0)] − 2 dt′ [H̄i (t), [H̄i (t′ ), ρ̄sb (t′ )] (26)
dt h̄ h̄ 0
donde ρ̄sb (t′ ) = ρ̄s (t′ ) ⊗ ρ̄b (0).
La ecuación (26) es aún una ecuación complicada debido a que la evolución de ρ̄sb (t) depende de la
historia del sistema pues hay que integrar sobre ρ̄sb (t′ ).
Consideremos ahora la hipótesis de Markoff la cual consiste en suponer que el comportamiento futuro
de un sistema dado depende solamente del estado actual de dicho sistema. Matemáticamente escribimos
Markoff reemplazando ρ̄sb (t′ ) por ρ̄sb (t).
Fı́sicamente la hipótesis de Markoff es bastante razonable. Es natural pensar que el sistema ρ̄s puede
depender de sus estados anteriores pues estos quedan grabados en la memoria del sistema completo a
través de la interacción H̄i ; dicha información se almacena en el baño térmico. Ahora, si el baño es
suficientemente grande y en equilibrio térmico, no esperamos que se vea afectado a largo plazo por la
interacción con ρ̄s . Esto es, el baño no se verá afectado el tiempo suficienta como para influir en el
estado actual del sistema de dos niveles. Entonces la validez de la hipótesis de Markoff depende del
tiempo de correlación del baño contra la escala de tiempo de que un cambio sea significante en ρ̄s .
Algunos argumentos al respecto de las escalas de los tiempos de correlación del baño y el sistema son
presentados en [2].
Con esto en mente, reescribimos (26) como:
Z t
dρ̄sb i 1
(t) = − [H̄i (t), ρ̄sb (0)] − 2 dt′ [H̄i (t), [H̄i (t′ ), ρ̄sb (t)] (27)
dt h̄ h̄ 0
Ahora si estamos en condiciones de continuar con el cálculo de la ecuación maestra. Aplicando una
traza parcial sobre los estados del baño térmico escribimos:
dρ̄s dρ̄sb
(t) = T rb ( (t)) (28)
dt dt
Para calcular la traza parcial es necesario conocer el operador que describe el sistema de osciladores
en equilibrio térmico del baño ρ̄b (0). Para hallarlo consideremos lo siguiente:
Sea |ni i el estado asociado al modo ωi del oscilador.

b†i bi |ni i = ni |ni i (29)

Tal que ni es el número de fotones en el modo ωi .

El estado del baño es descrito por la matriz de densidad:
X
ρ̄bi = Pni |ni ihni | (30)
ni

5
 donde la probabilidad Pni de que el modo i del campo esté excitado con ni fotones viene dado por el
factor de Boltzmann:
−nh̄ωi
e KT
Pni = P −nh̄ωi (31)
n e KT

−nh̄ωi −11
)
donde e KT ≈ e−(10 ≤ 1, por lo que escribimos, usando la serie geométrica:
−nh̄ωi −nh̄ωi
Pni = e KT (1 − e KT ) (32)
Entonces, utilizando las ecuaciones (29), (30) y (32) escribimos el estado del baño como:
P
ρ̄bi = ni Pni |ni ihni |

P −nih̄ωi −h̄ωi
= ni e KT (1 − e KT )|ni ihni | (33)
†
−(b bi )h̄ωi −h̄ωi
P i
= ni e
KT (1 − e KT )|ni ihni |

La ecuación anterior corresponde al estado del baño en la base de autoestados de b†i bi , que ahora
escribimos independientemente de la base como:
†
−(b bi )h̄ωi −h̄ωi
i
ρ̄b = Πi e KT (1 − e KT ) (34)

Con este resultado, volvamos a la ecuación de evolución (Ec.28). Veamos primeramente que la traza
del primer término en la ecuación (27) es cero:

T rb ([H̄i (t), ρ̄sb (0)]) = T rb ([σ † G(t) + σG† (t), ρ̄s (0) ⊗ ρ̄b (0)])

= T rb ([σ † G(t), ρ̄s (0) ⊗ ρ̄b (0)]) + C.C

= [σ † , ρ̄s ]T r(Gρ̄b ) + C.C

donde,
†
−(b bk )h̄ωk −h̄ωk
k
T r( j gj bj e(ω−ωj ) Πk e KT
P
T r(Gρ̄b ) = (1 − e KT ))
†
−ωk −(b bk )ωk
P (ω−ωj ) k
= j gj e (1 − e KT )T r(bj Πk e KT )
−h̄ωk −nkh̄ωk
gj e(ω−ωj ) (1 − e
P
= j
KT )Πk e KT T r(bj )

= 0
ya que, P
T r(bj ) = l hn1 ...ni ...|bj |n1 ...ni ...i =0

Pues |n1 ...ni ...i es la base de b†i bi y en esa base todos los elementos de la diagonal de b†i y bi son cero.
La contribución de los términos que involucran a T rb (σG† , ρs ⊗ ρb ) son cero por un argumento similar
al anterior. Luego T rb ([H̄i (t), ρ̄sb (0)]) = 0 y la ecuación de evolución queda de la siguiente forma:
Z t
dρ̄s 1
(t) = − 2 T rb ( dt′ [H̄i (t), [H̄i (t′ ), ρ̄sb (t)]]) (35)
dt h̄ 0
Dejemos de lado, por ahora, la integral. Trabajando con la traza se llega a que:

6
 T rb ( [H̄i (t), [H̄i (t′ ), ρ̄sb (t)] ])

= T rb ([σ † G(t) + σG† (t), [σ † G(t′ ) + σG† (t′ ), ρ̄s (t) ⊗ ρ̄b (0)]])

= −T r(G(t′ )G† (t)ρ¯b )[σ, ρ̄s σ † ] − T r(G(t)G† (t′ )ρ̄b )[σ ρ̄s , σ † ] − T r(G† (t′ )G(t)ρ̄b )[σ † , ρ̄s σ]

−T r(G† (t)G(t′ )ρ̄b )[σ † ρ̄s , σ]

Las propiedades del baño térmico entran en la ecuación maestra a través de las funciones de correlación
del baño: hG† (t)G(t)ib = T r(G† (t)G(t)ρ̄b ) y hG(t)G† (t)ib = T r(G(t)G† (t)ρ̄b ), las cuales se relacionan
de la forma siguiente: consideremos la cantidad G(t)G† (t). Recordando la relación de conmutación
[bi , b†j ] = δij , se tiene que:
′
G(t)G† (t) gi gj∗ bi b†j eit(ω−ωi ) eit (ω−ωj )
P
= ij

′
gi gj∗ (δij + b†j bi )eit(ω−ωi ) eit (ω−ωj )
P
= ij

′ ′
|gi |2 ei(ω−ωi )(t −t) + gi gj∗ b†j bi eit(ω−ωi ) eit (ω−ωj )
P P
= i ij

′
|gi |2 ei(ω−ωi )(t −t) + G† (t)G(t)
P
= i

Utilizando este último resultado y continuando el cálculo de la traza, escribimos:

T rb ( [H̄i (t), [H̄i (t′ ), ρ̄sb (t)]] )

= −T r(αρ̄b + G† (t′ )G(t)ρ̄b )[σ, ρ̄s σ † ] − T r(αρ̄b + G† (t)G(t′ )ρ̄b )[σ ρ̄s , σ † ]

−T r(G† (t′ )G(t)ρ̄b )[σ † , ρ̄s σ] − T r(G† (t)G(t′ )ρ̄b )[σ † ρ̄s , σ]
′
Con α = i |gi |2 ei(ω−ωi )(t −t) .
P

Consideremos ahora la siguiente cantidad, recordando que T r(Aρ) = hAi y que < b†i bi >=< ni >
′
T r[G† (t′ )G(t)ρ̄b ] = gi∗ gj b†i bj e−i(ω−ωi )t ei(ω−ωj )t ρ̄b ]
P
T r[ ij

′
|gi |2 ei(ω−ωi )(t −t) hni i
P
= i

Debido a que el caso i 6= j arroja un cero.

Entonces, utilizando este resultado y sustituyendo en la integral de la ecuación (35), teniendo en
cuenta que al integrar se tienen que considerar todas las frecuencias posibles para el conjunto de os-
ciladores acoplados del baño; pasamos al lı́mite continuo introduciendo ası́ la función de densidad D(ω)
tal que:

D(ω ′ )dω ′
P R
i 7→

|gi |2 7→ |g(ω ′ )|2 (36)

< ni >7→ n(ω ′ )

Escribimos a partir de (35):
 Rt R∞ ′
dρ̄s
dt (t) = − i
h̄2 [σ, ρ̄s σ † ] ( 0 dτ 0 dω ′ |g(ω ′ )|2 ei(ω−ω )τ D(ω ′ ))

Rt R∞ ′
+[σ ρ̄s , σ † ]( 0 dτ 0 dω ′ |g(ω ′ )|2 e−i(ω−ω )τ D(ω ′ ))
(37)
Rt R∞ ′
+([σ, ρ̄s σ ] + [σ ρ̄s , σ])( 0 dτ 0 dω ′ |g(ω ′ )|2 ei(ω−ω )τ n(ω ′ )D(ω ′ ))
† †

Rt R∞ ′
+([σ ρ̄s , σ † ] + [σ † , ρ̄s σ])( 0 dτ 0 dω ′ |g(ω ′ )|2 e−i(ω−ω )τ n(ω ′ )D(ω ′ )) )

7
 Donde se realizó el cambio de variable τ = t′ − t, dτ = dt′ . En la ecuación (37) tenemos dos tipos de
integrales:
Z t Z ∞
′
A= dτ dω ′ |g(ω ′ )|2 e±i(ω−ω )τ D(ω ′ ) (38)
0 0
Z t Z ∞
′
B= dτ dω ′ |g(ω ′ )|2 e±i(ω−ω )τ n(ω ′ )D(ω ′ ) (39)
0 0
Note que para valores grandes de τ la integral toma valores muy pequeños, mientras que para τ
pequeño la integral tiene valores considerables.
Fı́sicamente podemos pensar que el tiempo de correlación del baño (en equilibrio térmico) es prácticamente
cero respecto al tiempo de correlación del sistema [2], lo cual corrobora la validez de la hipótesis de Markoff
hecha anteriormente. Escribimos entonces:
Z ∞ Z t
′
′ 2
A = lim ′
dω |g(ω )| D(ω ) ′
dτ e±i(ω−ω )τ (40)
t→∞ 0 0
Como ambas integrales convergen, es posible invertir el orden de integración, donde el lı́mite de la
integral viene dado por el valor principal de Cauchy:
Z t
′ 1
lim dτ e±i(ω−ω )τ = πδ(ω − ω ′ ) ± P (41)
t→∞ 0 ω − ω′

Procediendo de la misma forma con la integral B, se llega a lo siguiente:

A = π|g(ω)|2 D(ω) ± Λ (42)

B = π|g(ω)|2 D(ω)n̄(ω) ± Γ (43)

donde,
∞
|g(ω ′ )|2 D(ω ′ )
Z
Λ=P dω ′ (44)
0 ω − ω′
∞
|g(ω ′ )|2 D(ω ′ )n(ω ′ )
Z
Γ=P dω ′ (45)
0 ω − ω′
Sustituyendo en la ecuación (37) se tiene que:
dρ̄s
dt (t) = − h̄i2 ([σ, ρ̄s σ † ] (π|g(ω)|2 D(ω) ± Λ)

+[σ ρ̄s , σ † ](π|g(ω)|2 D(ω) − Λ)

(46)
+([σ, ρ̄s σ † ] + [σ † ρ̄s , σ])(π|g(ω)|2 D(ω)n̄(ω) + Γ)

+([σ ρ̄s , σ † ] + [σ † , ρ̄s σ])(π|g(ω)|2 D(ω)n̄(ω) − Γ))

Reescribiendo (46) se llega a que la ecuación maestra para un sistema de dos niveles inmerso en un
baño térmico tiene la forma siguiente:
γ
ρ̄˙ s (t) = †
2 n̄(ω)([σ ρ̄s , σ] + [σ † , ρ̄s σ]) + Λ[σ † σ, ρ̄s ]
(47)
+ γ2 (1 + n̄(ω))([σ, ρ̄s σ † ] + [σ ρ̄s , σ † ])
con,
γ = 2π|g(ω)|2 D(ω) (48)

8
2 Problemas
1. Utilizando el software Mathematica, resuelva la ecuación maestra para un sistema de dos niveles
en interacción con un baño térmico (Ec.47) considerando lo siguiente:
 
a(t) c(t)
• El sistema es descrito por la matriz de densidad general: ρs (t) =
d(t) b(t)
• Λ = 0 y γ = 1.
A partir de las soluciones generales obtenidas, estudie y utilice los comandos Plot[.] y Manipu-
late[.] para analizar la evolución de las funciones a(t),b(t),c(t) y d(t), variando las condiciones
iniciales del sistema y tomando en cuenta que estas condiciones iniciales deben satisfacer ciertas
restricciones para que ρs (t) pueda ser una matriz de densidad. Considere 0 ≤ n ≤ 1. En partic-
ular, escoja una condición inicial y un valor fijo para n y grafique la evolución de las funciones
a(t),b(t),c(t) y d(t).
¿Cómo evolucionan los términos fuera de la diagonal? Note que estos términos caracterizan la
habilidad del sistema para mostrar interferencia.
2. Suponga ahora que tiene dos sistemas de dos niveles, ambos inmersos en un baño térmico. El
espacio de Hilbert del sistema completo viene dado en este caso por:

H = (HA ⊗ HB ) ⊗ Hbaño (49)

Donde HA y HB son los espacios de Hilbert de los subsistemas de dos niveles. Digamos que
definimos el espacio de Hilbert de ambos subsistemas sin el baño como HT = HA ⊗ HB . Entonces
el espacio de Hilbert del sistema completo se escribe como H = HT ⊗ Hbaño . Este espacio de
Hilbert es equivalente al espacio correspondiente a un único sistema en interacción con el medio
ambiente por lo que la ecuación maestra del par de sistemas 2 × 2 es similar a la ecuación (47):

γ
ρ̄˙ s (t) = †
2 n̄(ω)([σ ρ̄s , σ] + [σ † , ρ̄s σ]) + Λ[σ † σ, ρ̄s ]
(50)
+ γ2 (1 + n̄(ω))([σ, ρ̄s σ † ] + [σ ρ̄s , σ † ])

de manera que ρs es la matriz de densidad 4 × 4 que describe el par de sistemas de dos niveles y
σ ≡ σa + σb ; con σa , σa† y σb , σb† los operadores de subida y bajada para cada subsistema.
Resuelva usando Mathematica la ecuación maestra para dos sistemas 2 × 2 en un baño común
(Ec.50) considerando lo siguiente:

• Los dos subsistemas de dos niveles son descritos por:

 
1 + c(t) + 2d(t) 0 0 a(t) − b(t)
1 0 1 − c(t) a(t) + b(t) 0 
ρs (t) =   (51)
4 0 a(t) + b(t) 1 − c(t) 0 
a(t) − b(t) 0 0 1 + c(t) − 2d(t)
• Λ = 0, n = 0 y γ = 1.
• d(0) = 0.

A partir de las soluciones generales obtenidas para las funciones a(t),b(t),c(t) y d(t), estudie la
evolución de la siguiente función:
p p
C(t) = 2(|a(t) + b(t)| − 1 + c(t) + 2d(t) 1 + c(t) − 2d(t)) (52)

Tome como condiciones iniciales el siguiente subconjunto de estados: a(0) = b(0) = c(0) = 13 − 4F
3
tal que 0 ≤ F ≤ 1. Grafique la evolución de C(t) para tres valores particulares de F : tome F = 1
y otros dos casos de su preferencia.

9
 La función C(t) es conocida como concurrencia del sistema y nos permite conocer el grado de
entrelazamiento cuántico que éste posee. A partir de sus resultados, ¿qué puede decir de la evolución
del entrelazamiento del sistema al interactuar con un baño térmico común?; ¿qué sucede con el
entrelazamiento para la condición inicial con F = 1?
Lo que usted ha encontrado es un subespacio libre de decoherencia.

References
[1] Cohen-Tannoudji C., Dui B., Laloë F. Quantum Mechanics Volume II. Hermann
(1977).
[2] Carmichael H.J., ULB Lectures in Nonlinear Optics Fall 1991: An Open Systems
Approach to Quantum Optics. University of Oregon, Eugene, Oregon 97403.
[3] Joos, E., H. D. Zeh, Z. The emergence of classical properties through interaction
with the environment Phys. B: Condens. Matter 59, 223. (1985)
[4] Joos, E., H. D. Zeh, C. Kiefer, D. Giulini, J. Kupsch, and I.O. Stamatescu,
Decoherence and the Appearance of a Classical World in Quantum Theory, 2nd
ed. Springer, New York. (2003)
[5] Orszag Miguel. Quantum Optics. Edit Springer-Verlag Berlin Heidelberg 2000
[6] Zurek, W. H. Phys. Today 44 10, 36. Decoherence and transition from quantum
to classical. quant-ph/0306072. (1991)
[7] Zurek, W. H. Rev. Mod. Phys. 75, 715. (2003)

10