Partícula cuántica neutra en campo magnético no homogéneo

Estoy tratando de entender el experimento de Stern-Gerlach a nivel computacional. Supongamos que tenemos una partícula neutra con momento magnético (por ejemplo, un neutrón) y le aplicamos un campo magnético no homogéneo (dejemos que cambie linealmente con la coordenada). Según tengo entendido, su hamiltoniano se vería así:

H ^ = 2 2 metro 2 + ( mi metro C ) s ^ B

Ahora el operador de giro es

s ^ i = 2 σ i ,

dónde σ i es i ª matriz de Pauli .

Entonces, para el campo magnético B = mi X B 0 X tendríamos la ecuación de Schrödinger 1D (las direcciones Y y Z se pueden separar debido a la simetría de traducción):

2 2 metro 2 ψ X 2 + ( mi 2 metro C ) σ X B 0 X ψ = i ψ t .

Ahora trato de resolver esta ecuación numéricamente, tomando la función de onda inicial de la siguiente forma:

ψ ( X , t = 0 ) = ( ψ 0 ( X ) ψ 0 ( X ) ) ,

dónde ψ 0 ( X ) es un paquete de ondas gaussianas con impulso promedio cero.

Los problemas comienzan cuando selecciono σ X como suele darse:

σ X = ( 0 1 1 0 ) .

La solución parece verse como se muestra a continuación. Es decir, ¡ambos componentes de la función de onda aceleran hacia la izquierda!

ingrese la descripción de la imagen aquí

Pensé, ¿y si elijo otro eje como X , así que intenté hacer lo mismo con σ y :

σ y = ( 0 i i 0 ) .

El resultado en la animación de abajo. Ahora es un poco mejor: la función de onda al menos se divide en dos partes, una hacia la izquierda y otra hacia la derecha. Pero aún así, ambas partes se componen de una mezcla de estados de giro hacia arriba y hacia abajo, por lo que no es realmente lo que uno esperaría del experimento de Stern-Gerlach.

ingrese la descripción de la imagen aquí

Finalmente, probé la última opción: usar σ z :

σ z = ( 1 0 0 1 ) .

El resultado se muestra nuevamente a continuación. Finalmente, obtengo la división en partes de giro "independientes", es decir, una parte de giro va a la izquierda, otra va a la derecha.

ingrese la descripción de la imagen aquí

Ahora, la pregunta : ¿cómo interpretar estos resultados? ¿Por qué la elección del eje activo da como resultado diferencias tan drásticas en los resultados? ¿Cómo debería haberlo hecho para obtener resultados significativos? ¿No debería la permutación de las matrices de Pauli afectar los resultados?

Respuestas (1)

Voy a hacer algunos cambios cosméticos a tu ecuación:

  1. Su campo magnético debe ser B = mi X B 0 X / X 0 , de modo que B y B 0 tienen las mismas dimensiones. Este campo no obedece B = 0 , pero dejaremos eso solo por ahora.

  2. Escribiré m σ B por la energía magnética. Observe que el momento magnético del neutrón es bastante pequeño, alrededor de 50 neV/T. Los llamados neutrones "ultrafríos", con energías por debajo de unos 100 neV, pueden ser confinados por un mínimo de campo magnético, pero para los neutrones fríos o térmicos ~meV, la dirección de Stern-Gerlach es generalmente insignificante. (Si no fuera insignificante, podría usar la separación de Stern-Gerlach para convertir un haz de neutrones en dos haces polarizados que van en diferentes direcciones; por desgracia, todos los polarizadores de neutrones reales absorben un estado de espín).

  3. Tu ecuación de movimiento también es separable en el tiempo, así que solo consideraré la parte espacial; podemos agregar un factor de mi i ω t más tarde.

  4. voy a usar ϕ y x para dar nombres separados a los dos componentes de su espinor.

En notación matricial, su ecuación de Schrödinger independiente del tiempo es

( 2 2 metro X 2 m B 0 X X 0 m B 0 X X 0 2 2 metro X 2 ) ( ϕ x ) = mi ( ϕ x ) .

Esto hace que sea un poco más obvio lo que está sucediendo. Cuando definió su conjunto inicial para tener ϕ = x = ψ 0 ( X ) , establece su sistema en un estado propio de σ X ! Entonces, por supuesto, todo el conjunto se mueve junto: ¡está polarizado a lo largo de la dirección del campo! Si su condición inicial es ϕ = x , debería obtener un paquete moviéndose en la otra dirección.

cuando reemplazas σ X con σ y , z , está cambiando efectivamente la dirección del campo. Recuerda que la energía es m σ B ; si el único término distinto de cero en este producto es σ y B 0 X / X 0 , le está diciendo a su modelo que el campo varía en fuerza con X pero puntos a lo largo del y dirección. (Dado que esos campos obedecen B = 0 , deberías preferirlas de todos modos). Así que en tu segunda y tercera figura estás poniendo una X muestra polarizada en un campo a lo largo y o a lo largo z , y se separa. La separación es a lo largo de la dirección del gradiente de intensidad de campo, en lugar de a lo largo de la dirección del campo; los tramperos de la UCN hablan de "buscadores de campo fuerte" y "buscadores de campo débil".

En el campo apuntando a lo largo z , su muestra se separa en su base de espinor, lo que le da los bonitos colores.

En el campo apuntando a lo largo y , su muestra aún se separa. Pero en lugar de separarse en el ϕ , x base que está utilizando para los colores, los estados diagonales son ϕ ± i x . Esto genera interferencia cuando los dos paquetes de ondas se superponen.

Vaya, esto es interesante. ¿Podría explicar por qué la configuración ϕ = x hace que la función de onda sea un estado propio de σ X ? Pensé que sería un estado propio si configuro, por ejemplo x = 0 , ϕ 0 .
Tarea: demostrar que los vectores propios de σ X son ( 1 , ± 1 ) , y que los vectores propios de σ y son ( 1 , ± i ) . Encuentre los vectores propios de σ z . :-)
OK, parece que empiezo a entenderlo. Pero, ¿por qué se separa en base de espinor para z -campo dirigido, pero no para y -dirigido? ¿No debería ser lo mismo, porque podríamos rotar el marco de referencia alrededor mi X mientras cambia y -campo a z -campo y no han cambiado nada?
Ah, parece que entiendo: esta diferencia no es medible: la densidad de probabilidad de los electrones no cambia, por lo que este no es un problema real.
Creo que estoy de acuerdo contigo: la única razón por la que la respuesta al campo a lo largo z se ve especial es que por casualidad elegiste el z eje para su representación.