Un sistema de dos niveles bajo Ito-calculus

Estoy considerando la evolución de un sistema cuántico de dos niveles dado por

i ( C ˙ 1 C ˙ 2 ) = ( 1 2 mi Δ ( t ) Δ ( t ) 1 2 mi ) ( C 1 C 2 ) .

Si Δ ( t ) es periódica, el modelo da oscilaciones de Rabi. Pero ahora estoy considerando la función. Δ ( t ) para ser un verdadero ruido blanco gaussiano satisfactorio

Δ ( t ) = 0 , Δ ( t ) Δ ( t ) = Γ 2 d ( t t ) .

En términos del proceso estándar de Wiener W t , la ecuación de Schrödinger se convierte en

i d C 1 = 1 2 mi C 1 d t + Γ C 2 d W t , i d C 2 = Γ C 1 d W t 1 2 mi C 2 d t ,

donde ambos C 1 d W t y C 2 d W t están en el sentido Ito, lo que significa que C 1 y C 2 están en t y d W t es durante [ t , t + d t ] . Luego voy a la imagen de interacción.

C ~ 1 = C 1 mi i 2 mi t , C ~ 2 = C 2 mi i 2 mi t ,

de modo que C ~ 1 y C ~ 2 no evoluciones cuando Γ = 0 . En general, satisfacen

i d C ~ 1 = Γ C ~ 2 mi i mi t d W t , i d C ~ 2 = Γ C ~ 1 mi i mi t d W t .

Ahora, ¿cómo procedo desde aquí para expresar C ~ 1 y C ~ 2 en términos de alguna integral de W t ? También me pregunto si obtengo decoherencia o solo oscilación Rabi porque el sistema absorbe la frecuencia de forma selectiva. mi del ruido W t .

no es tarea...
Bien. Haz que parezca más fácil para que la gente pueda resolverlo.

Respuestas (1)

Para resolver estas ecuaciones

i d C ~ 1 = Γ C ~ 2 mi i mi t d W t , i d C ~ 2 = Γ C ~ 1 mi i mi t d W t .

veamos si podemos encontrar un par de funciones F 1 ( t , X ) y F 2 ( t , X ) tal que C ~ 1 = F 1 ( t , W t ) y C ~ 2 = F 2 ( t , W t ) . Conectando estos en las fórmulas que tenemos

i d F 1 ( t , W t ) = Γ F 2 ( t , W t ) mi i mi t d W t , i d F 2 ( t , W t ) = Γ F 1 ( t , W t ) mi i mi t d W t .

Ahora aplicamos el cálculo de Itô a las diferenciales

i [ F 1 t ( t , W t ) d t + F 1 X ( t , W t ) d W t + 1 / 2 F 1 X X ( t , W t ) d t ] = Γ F 2 ( t , W t ) mi i mi t d W t , i [ F 2 t ( t , W t ) d t + F 2 X ( t , W t ) d W t + 1 / 2 F 2 X X ( t , W t ) d t ] = Γ F 1 ( t , W t ) mi i mi t d W t .

En estos, las letras en la notación de índice son abreviaturas de derivadas parciales con las variables correspondientes. Por lo tanto, obtenemos las siguientes ecuaciones diferenciales parciales

F 1 t ( t , X ) + 1 / 2 F 1 X X ( t , X ) = 0 , F 2 t ( t , X ) + 1 / 2 F 2 X X ( t , X ) = 0 , i F 1 X ( t , X ) = Γ F 2 ( t , X ) mi i mi t , i F 2 X ( t , X ) = Γ F 1 ( t , X ) mi i mi t .

Tomando la derivada de X de la tercera ecuación y multiplicando por i obtenemos

F 1 X X ( t , X ) = Γ i F 2 X ( t , X ) mi i mi t = Γ 2 F 1 ( t , X )

donde en el último paso sustituyo en la cuarta ecuación. Finalmente reemplazando esto en la primera ecuación obtenemos

F 1 t ( t , X ) = Γ 2 2 F 1 ( t , X )

que es una fácil ecuación diferencial parcial lineal de primer orden en t . La solucion es

F 1 ( t , X ) = A ( X ) mi Γ 2 t / 2 .

Con A ( X ) alguna función para determinar a partir de las otras ecuaciones. Por ejemplo, llenando el resultado para F 1 en la primera fórmula, obtenemos la siguiente ecuación diferencial para A ( X ) :

A X X ( X ) + Γ 2 A ( X ) = 0

cuyas soluciones son

A ( X ) = B mi i Γ X + D mi i Γ X

y por lo tanto

F 1 ( t , X ) = ( B mi i Γ X + D mi i Γ X ) mi Γ 2 t / 2

y de la misma manera

F 2 ( t , X ) = ( B mi i Γ X + D mi i Γ X ) mi ( Γ 2 / 2 i mi ) t .

De este modo

C ~ 1 = ( B mi i Γ W t + D mi i Γ W t ) mi Γ 2 t / 2 , C ~ 2 = ( B mi i Γ W t + D mi i Γ W t ) mi ( Γ 2 / 2 i mi ) t .

Creo que uno puede mostrar | C ~ 1 | 2 + | C ~ 2 | 2 = C o norte s t , por lo que las amplitudes de probabilidad no deberían crecer exponencialmente. Pero probaré la idea de C ~ i = F i ( t , W t ) , i = 1 , 2 .
También me pareció extraño, pero podría haber un error de señal en alguna parte. Además, acabo de notar una discrepancia en su formulación donde el cambio en la imagen de interacción tiene un 1 / 2 en los exponentes, pero esos no aparecen en sus ecuaciones a continuación.
Noté un error en mi cálculo. olvidé un 1 / 2 factor yo mismo. voy a arreglarlo
De hecho C ~ i = F i ( t , W t ) es una suposición muy fuerte. asume C ~ i no depende de antes W t con t < t . Si conecta su solución a la ecuación y encuentra que realmente funciona, entonces está bien. De lo contrario, la solución debería ser de alguna forma integral.
Ese es un buen punto, pero no creo que sea el problema aquí.
Descubrí cuál es el problema. Ambos estábamos tan seguros de que | C ~ 1 | 2 + | C ~ 2 | 2 = C o norte s t , pero en realidad no es cierto. Con el lema de Itô, puedes demostrar que d | C ~ 1 | 2 = C ~ 1 d C ~ 1 + C ~ 1 d C ~ 1 + d C ~ 1 d C ~ 1 = C ~ 1 d C ~ 1 + C ~ 1 d C ~ 1 + Γ 2 | C ~ 2 | 2 d t y del mismo modo para d | C ~ 2 | 2 . esto le dará d ( | C ~ 1 | 2 + | C ~ 2 | 2 ) = Γ 2 ( | C ~ 1 | 2 + | C ~ 2 | 2 ) d t lo cual es consistente con mi respuesta.
Gracias. También me di cuenta de esto. A la ecuación de Schrödinger le falta el d t 2 término de Exp ( i H ( Δ ) d t ) porque Δ O ( d t 1 / 2 ) . Entonces el estado permanece normalizado.
La decoherencia ocurre en el sentido de conjunto. Para cada trayectoria aleatoria individual, el spin- 1 / 2 el estado puro sigue siendo un estado puro bajo evoluciones unitarias. Pero como el ruido es desconocido, el estado final es una mezcla aleatoria de muchas posibilidades, perdiendo así coherencia. El mejor formalismo para manejar la mezcla aleatoria es usar operadores de Lindblad. Eventualmente, la matriz de densidad evoluciona al estado incoherente 1 2 ( | | + | | ) .