Potencial del oscilador armónico, ¿prueba de que las gaussianas siguen siendo gaussianas?

Leí en varios artículos que para un oscilador armónico hamiltoniano en la ecuación de Schrödinger dependiente del tiempo, un paquete de ondas gaussianas sigue siendo gaussiana.

Desafortunadamente, no pude encontrar ninguna prueba para esta afirmación y al tratar de verificarla yo mismo no tuve éxito. Si hago un ansatz general con un paquete de ondas gaussianas esféricamente simétrico con ancho dependiente del tiempo y fase dependiente del tiempo y el espacio

ψ ( t , X ) = ( π a ( t ) 2 ) 3 / 4 Exp ( X 2 2 a ( t ) 2 + i ϕ ( t , X ) )
e insertarlo en la ecuación de Schrödinger
i ψ ˙ ( t , X ) = 1 2 metro Δ ψ ( t , X ) + k 2 X 2 ψ ( t , X )
Obtengo ecuaciones diferenciales relativamente complicadas que involucran las primeras derivadas de a y ϕ así como derivadas espaciales de primer y segundo orden de ϕ . No pude resolver esas ecuaciones o incluso mostrar que existe una solución.

¿Hay una manera fácil de mostrar esto?

¿Hay alguna referencia donde se muestre esto?

¿Cuáles son las soluciones para a ( t ) y ϕ ( t , X ) para una Gaussiana (dadas las condiciones iniciales a ( 0 ) = σ y ϕ ( 0 , X ) = 0 )?

¿Quizás podría intentar ir al espacio de Fourier, utilizando la descomposición espacial de Fourier?
@Danu: ¿Eso no me daría exactamente las mismas ecuaciones? Dado que la transformada de Fourier de una gaussiana es nuevamente una gaussiana y en la hamiltoniana X 2 y pag 2 aparecer simétricamente?
He visto una transformada de Fourier simplificar algunos cálculos que involucran gaussianos antes, pero depende de usted si quiere intentarlo.
Más información sobre paquetes de ondas gaussianas: physics.stackexchange.com/search?q=Gaussian+wave+packet

Respuestas (1)

Hay un poco más de estructura en el problema, y ​​le recomendaría que lo aproveche. En particular, sabes que ϕ es realmente una función lineal de X , o de lo contrario no es un gaussiano; y que las partes real e imaginaria de su coeficiente tienen un significado físico distinto, que puede obtener tomando valores esperados apropiados. Si aplica esto, entonces, puede hacer que su Ansatz sea

X | ψ ( t ) = ψ ( t , X ) = ( π a ( t ) 2 ) 3 / 4 Exp ( ( X X 0 ( t ) ) 2 2 a ( t ) 2 + i pag 0 ( t ) · X / ħ ) .
Aún mejor, puede calcular fácilmente los valores esperados de posición y momento para ser
ψ ( t ) | X | ψ ( t ) = X 0 ( t )
y
ψ ( t ) | pag | ψ ( t ) = pag 0 ( t ) ,
y puede aplicar el teorema de Ehrenfest para obtener ecuaciones de movimiento significativas y fáciles de resolver para estas cantidades. Esto restringe completamente su solución.

Por supuesto, esto debería dejarlo un poco incómodo porque no ha probado que ψ ( t , X ) es una solución del TDSE... pero descubrirlo está a unas pocas diferenciaciones de distancia. Como sabe que es una solución de todos modos, está en un terreno bastante seguro.