Parte de una serie de artículos sobre |
Mecánica cuántica |
---|
This article includes a list of references, related reading, or external links, but its sources remain unclear because it lacks inline citations. (November 2023) |
El oscilador armónico cuántico es el análogo mecánico cuántico del oscilador armónico clásico . Debido a que un potencial suave arbitrario puede ser aproximado generalmente como un potencial armónico en la proximidad de un punto de equilibrio estable , es uno de los sistemas modelo más importantes en mecánica cuántica. Además, es uno de los pocos sistemas mecánicos cuánticos para los que se conoce una solución analítica exacta. [1] [2] [3]
El hamiltoniano de la partícula es: donde m es la masa de la partícula, k es la constante de fuerza, es la frecuencia angular del oscilador, es el operador de posición (dado por x en la base de coordenadas), y es el operador de momento (dado por en la base de coordenadas). El primer término del hamiltoniano representa la energía cinética de la partícula, y el segundo término representa su energía potencial, como en la ley de Hooke . [4]
La ecuación de Schrödinger independiente del tiempo (TISE) es, donde denota un número real (que necesita ser determinado) que especificará un nivel de energía independiente del tiempo , o valor propio , y la solución denota el estado propio de energía de ese nivel . [5]
Luego, resuelva la ecuación diferencial que representa este problema de valores propios en la base de coordenadas, para la función de onda , utilizando un método espectral . Resulta que existe una familia de soluciones. En esta base, ascienden a funciones de Hermite , [6] [7]
Las funciones H n son los polinomios de Hermite de los físicos ,
Los niveles de energía correspondientes son [8] Los valores esperados de posición y momento combinados con la varianza de cada variable se pueden derivar de la función de onda para comprender el comportamiento de los eigenkets de energía. Se muestra que son y debido a la simetría del problema, mientras que:
Se observa que la variación tanto en la posición como en el momento aumenta para los niveles de energía más altos. El nivel de energía más bajo tiene un valor que es su valor mínimo debido a la relación de incertidumbre y también corresponde a una función de onda gaussiana. [9]
Este espectro de energía es notable por tres razones. Primero, las energías están cuantizadas, lo que significa que solo son posibles valores de energía discretos (múltiplos enteros más la mitad de ħω ); esta es una característica general de los sistemas mecánicos cuánticos cuando una partícula está confinada. Segundo, estos niveles de energía discretos están espaciados de manera uniforme, a diferencia del modelo de Bohr del átomo o la partícula en una caja . Tercero, la energía más baja alcanzable (la energía del estado n = 0 , llamado estado fundamental ) no es igual al mínimo del pozo de potencial, sino ħω /2 por encima de él; esto se llama energía de punto cero . Debido a la energía de punto cero, la posición y el momento del oscilador en el estado fundamental no son fijos (como lo serían en un oscilador clásico), sino que tienen un pequeño rango de varianza, de acuerdo con el principio de incertidumbre de Heisenberg .
La densidad de probabilidad del estado fundamental se concentra en el origen, lo que significa que la partícula pasa la mayor parte de su tiempo en el fondo del pozo de potencial, como se esperaría de un estado con poca energía. A medida que aumenta la energía, la densidad de probabilidad alcanza su punto máximo en los "puntos de inflexión" clásicos, donde la energía del estado coincide con la energía potencial (véase la discusión a continuación de los estados altamente excitados). Esto es coherente con el oscilador armónico clásico, en el que la partícula pasa la mayor parte de su tiempo (y, por lo tanto, es más probable que se la encuentre) cerca de los puntos de inflexión, donde se mueve más lentamente. De este modo, se satisface el principio de correspondencia . Además, los paquetes de ondas no dispersivos especiales , con una incertidumbre mínima, llamados estados coherentes , oscilan de forma muy similar a los objetos clásicos, como se ilustra en la figura; no son estados propios del hamiltoniano.
El método del " operador de escalera ", desarrollado por Paul Dirac , permite la extracción de los valores propios de energía sin resolver directamente la ecuación diferencial. [10] Es generalizable a problemas más complicados, en particular en la teoría cuántica de campos . Siguiendo este enfoque, definimos los operadores a y su adjunto a † . Nótese que estos operadores clásicamente son exactamente los generadores de rotación normalizada en el espacio de fase de y , es decir , describen la evolución hacia adelante y hacia atrás en el tiempo de un oscilador armónico clásico. [ aclaración necesaria ]
Estos operadores conducen a la siguiente representación de y ,
El operador a no es hermítico , ya que él mismo y su adjunto a † no son iguales. Los estados propios de energía | n ⟩ , cuando se operan con estos operadores de escalera, dan
A partir de las relaciones anteriores, también podemos definir un operador numérico N , que tiene la siguiente propiedad:
Los siguientes conmutadores se pueden obtener fácilmente sustituyendo la relación de conmutación canónica ,
y el operador de Hamilton se puede expresar como
Por lo tanto, los estados propios de N son también los estados propios de la energía. Para comprobarlo, podemos aplicar a un estado numérico :
Usando la propiedad del operador numérico :
obtenemos:
Por lo tanto, dado que resuelve el TISE para el operador hamiltoniano , es también uno de sus estados propios con el valor propio correspondiente:
QED.
La propiedad de conmutación produce
y de manera similar,
Esto significa que a actúa sobre | n ⟩ para producir, hasta una constante multiplicativa, | n –1⟩ , y a † actúa sobre | n ⟩ para producir | n +1⟩ . Por esta razón, a se denomina operador de aniquilación ("operador de reducción"), y a † operador de creación ("operador de elevación"). Los dos operadores juntos se denominan operadores de escalera .
Dado cualquier estado propio de energía, podemos actuar sobre él con el operador de reducción, a , para producir otro estado propio con ħω menos energía. Mediante la aplicación repetida del operador de reducción, parece que podemos producir estados propios de energía hasta E = −∞ . Sin embargo, dado que
El valor propio más pequeño del operador numérico es 0, y
En este caso, las aplicaciones posteriores del operador de reducción producirán simplemente cero, en lugar de estados propios de energía adicionales. Además, hemos demostrado anteriormente que
Finalmente, actuando sobre |0⟩ con el operador de elevación y multiplicando por factores de normalización adecuados , podemos producir un conjunto infinito de estados propios de energía.
de manera que coincida con el espectro de energía dado en la sección anterior.
Los estados propios arbitrarios se pueden expresar en términos de |0⟩, [11]
El análisis precedente es algebraico, utilizando únicamente las relaciones de conmutación entre los operadores de elevación y descenso. Una vez que se completa el análisis algebraico, se deben abordar cuestiones analíticas. En primer lugar, se debe encontrar el estado fundamental, es decir, la solución de la ecuación . En la representación de posición, esta es la ecuación diferencial de primer orden cuya solución se encuentra fácilmente como la gaussiana [nb 1] Conceptualmente, es importante que solo haya una solución de esta ecuación; si hubiera, digamos, dos estados fundamentales linealmente independientes, obtendríamos dos cadenas independientes de vectores propios para el oscilador armónico. Una vez que se calcula el estado fundamental, se puede demostrar inductivamente que los estados excitados son polinomios de Hermite multiplicados por el estado fundamental gaussiano, utilizando la forma explícita del operador de elevación en la representación de posición. También se puede demostrar que, como se esperaba a partir de la unicidad del estado fundamental, los estados propios de energía de las funciones de Hermite construidos por el método de escalera forman un conjunto ortonormal completo de funciones. [12]
Conectando explícitamente con la sección anterior, el estado fundamental |0⟩ en la representación de la posición está determinado por , por lo tanto , de modo que , y así sucesivamente.
El oscilador armónico cuántico posee escalas naturales de longitud y energía que pueden utilizarse para simplificar el problema. Estas pueden encontrarse mediante la adimensionalización .
El resultado es que, si la energía se mide en unidades de ħω y la distancia en unidades de √ ħ /( mω ) , entonces el hamiltoniano se simplifica a mientras que las funciones propias de energía y los valores propios se simplifican a funciones de Hermite y números enteros desplazados por un medio, donde H n ( x ) son los polinomios de Hermite .
Para evitar confusiones, en este artículo no se adoptarán en su mayoría estas "unidades naturales". Sin embargo, suelen resultar útiles al realizar cálculos, ya que evitan el desorden.
Por ejemplo, la solución fundamental ( propagador ) de H − i∂ t , el operador de Schrödinger dependiente del tiempo para este oscilador, simplemente se reduce al núcleo de Mehler , [13] [14] donde K ( x , y ;0) = δ ( x − y ) . La solución más general para una configuración inicial dada ψ ( x ,0) es entonces simplemente
Los estados coherentes (también conocidos como estados de Glauber) del oscilador armónico son paquetes de ondas no dispersivos especiales , con una incertidumbre mínima σ x σ p = ℏ ⁄ 2 , cuyos valores esperados de observables evolucionan como un sistema clásico. Son vectores propios del operador de aniquilación, no del hamiltoniano, y forman una base sobrecompleta que, en consecuencia, carece de ortogonalidad. [15]
Los estados coherentes están indexados y expresados en la base | n ⟩ como
Dado que los estados coherentes no son estados propios de energía, su evolución temporal no es un simple cambio de fase de la función de onda. Sin embargo, los estados evolucionados en el tiempo también son estados coherentes, pero con un parámetro de cambio de fase α : .
Debido a que y a través de la identidad de Kermack-McCrae, la última forma es equivalente a un operador de desplazamiento unitario que actúa sobre el estado fundamental: . Cálculo de los valores esperados:
donde es la fase aportada por el complejo α . Estas ecuaciones confirman el comportamiento oscilatorio de la partícula.
Las incertidumbres calculadas mediante el método numérico son:
lo que da . Dado que la única función de onda que puede tener la incertidumbre de posición-momento más baja, , es una función de onda gaussiana, y dado que la función de onda de estado coherente tiene una incertidumbre de posición-momento mínima, observamos que la función de onda gaussiana general en mecánica cuántica tiene la forma: Sustituyendo los valores esperados como una función del tiempo, da la función de onda variable en el tiempo requerida.
La probabilidad de cada estado propio de energía se puede calcular para encontrar la distribución de energía de la función de onda:
que corresponde a la distribución de Poisson .
Cuando n es grande, los estados propios se localizan en la región clásica permitida, es decir, la región en la que una partícula clásica con energía E n puede moverse. Los estados propios alcanzan su máximo cerca de los puntos de inflexión: los puntos en los extremos de la región clásicamente permitida donde la partícula clásica cambia de dirección. Este fenómeno se puede verificar a través de la asintótica de los polinomios de Hermite y también a través de la aproximación WKB .
La frecuencia de oscilación en x es proporcional al momento p ( x ) de una partícula clásica de energía E n y posición x . Además, el cuadrado de la amplitud (que determina la densidad de probabilidad) es inversamente proporcional a p ( x ) , lo que refleja el tiempo que la partícula clásica pasa cerca de x . El comportamiento del sistema en una pequeña vecindad del punto de inflexión no tiene una explicación clásica simple, pero se puede modelar utilizando una función de Airy . Utilizando las propiedades de la función de Airy, se puede estimar que la probabilidad de encontrar la partícula fuera de la región clásicamente permitida es aproximadamente Esto también se da, asintóticamente, por la integral
En la formulación del espacio de fases de la mecánica cuántica, los estados propios del oscilador armónico cuántico en varias representaciones diferentes de la distribución de cuasiprobabilidad se pueden escribir en forma cerrada. La más utilizada de ellas es la distribución de cuasiprobabilidad de Wigner .
La distribución de cuasiprobabilidad de Wigner para el estado propio de energía | n ⟩ es, en las unidades naturales descritas anteriormente, [ cita requerida ] donde L n son los polinomios de Laguerre . Este ejemplo ilustra cómo los polinomios de Hermite y Laguerre están vinculados a través de la función de Wigner .
Mientras tanto, la función Q de Husimi de los estados propios del oscilador armónico tiene una forma aún más simple. Si trabajamos en las unidades naturales descritas anteriormente, tenemos Esta afirmación se puede verificar utilizando la transformada de Segal-Bargmann . Específicamente, dado que el operador de elevación en la representación de Segal-Bargmann es simplemente la multiplicación por y el estado fundamental es la función constante 1, los estados del oscilador armónico normalizado en esta representación son simplemente . En este punto, podemos apelar a la fórmula para la función Q de Husimi en términos de la transformada de Segal-Bargmann.
El oscilador armónico unidimensional se puede generalizar fácilmente a N dimensiones, donde N = 1, 2, 3, … . En una dimensión, la posición de la partícula se especificó mediante una única coordenada , x . En N dimensiones, esto se reemplaza por N coordenadas de posición, que etiquetamos x 1 , …, x N . A cada coordenada de posición corresponde un momento; los etiquetamos p 1 , …, p N . Las relaciones de conmutación canónicas entre estos operadores son
El hamiltoniano para este sistema es
Como lo deja claro la forma de este hamiltoniano, el oscilador armónico N -dimensional es exactamente análogo a N osciladores armónicos unidimensionales independientes con la misma masa y constante de resorte. En este caso, las cantidades x 1 , ..., x N se referirían a las posiciones de cada una de las N partículas. Esta es una propiedad conveniente del potencial r 2 , que permite separar la energía potencial en términos que dependen de una coordenada cada uno.
Esta observación hace que la solución sea sencilla. Para un conjunto particular de números cuánticos, las funciones propias de energía para el oscilador N -dimensional se expresan en términos de las funciones propias unidimensionales como:
En el método del operador de escalera, definimos N conjuntos de operadores de escalera,
Mediante un procedimiento análogo al caso unidimensional, podemos demostrar que cada uno de los operadores a i y a † i disminuyen y aumentan la energía en ℏω respectivamente. El hamiltoniano es Este hamiltoniano es invariante bajo el grupo de simetría dinámica U ( N ) (el grupo unitario en N dimensiones), definido por donde es un elemento en la representación matricial definitoria de U ( N ) .
Los niveles de energía del sistema son
Al igual que en el caso unidimensional, la energía está cuantizada. La energía del estado fundamental es N veces la energía fundamental unidimensional, como cabría esperar utilizando la analogía con N osciladores unidimensionales independientes. Hay una diferencia más: en el caso unidimensional, cada nivel de energía corresponde a un estado cuántico único. En N dimensiones, a excepción del estado fundamental, los niveles de energía están degenerados , lo que significa que hay varios estados con la misma energía.
La degeneración se puede calcular con relativa facilidad. Como ejemplo, considere el caso tridimensional: Defina n = n 1 + n 2 + n 3 . Todos los estados con el mismo n tendrán la misma energía. Para un n dado , elegimos un n 1 particular . Entonces n 2 + n 3 = n − n 1 . Hay n − n 1 + 1 pares posibles { n 2 , n 3 } . n 2 puede tomar los valores de 0 a n − n 1 , y para cada n 2 el valor de n 3 es fijo. Por lo tanto, el grado de degeneración es: Fórmula para N y n generales [ siendo g n la dimensión de la representación de potencia n -ésima irreducible simétrica del grupo unitario U ( N ) ]: El caso especial N = 3, dado anteriormente, se sigue directamente de esta ecuación general. Sin embargo, esto solo es cierto para partículas distinguibles, o una partícula en N dimensiones (ya que las dimensiones son distinguibles). Para el caso de N bosones en una trampa armónica unidimensional, la degeneración se escala como el número de formas de particionar un entero n utilizando números enteros menores o iguales a N.
Esto surge debido a la restricción de poner N cuantos en un estado ket donde y , que son las mismas restricciones que en la partición de enteros.
La ecuación de Schrödinger para una partícula en un oscilador armónico tridimensional con simetría esférica se puede resolver explícitamente mediante la separación de variables. Este procedimiento es análogo a la separación realizada en el problema del átomo similar al hidrógeno , pero con un potencial con simetría esférica diferente, donde μ es la masa de la partícula. Como m se utilizará a continuación para el número cuántico magnético, la masa se indica mediante μ , en lugar de m , como se indicó anteriormente en este artículo.
La solución de la ecuación es: [16] donde
son polinomios de Laguerre generalizados ; El orden k del polinomio es un entero no negativo;
El valor propio de la energía es La energía generalmente se describe mediante un único número cuántico
Como k es un entero no negativo, para cada n par tenemos ℓ = 0, 2, …, n − 2, n y para cada n impar tenemos ℓ = 1, 3, …, n − 2, n . El número cuántico magnético m es un entero que satisface − ℓ ≤ m ≤ ℓ , por lo que para cada n y ℓ hay 2 ℓ + 1 estados cuánticos diferentes , etiquetados por m . Por lo tanto, la degeneración en el nivel n es donde la suma comienza desde 0 o 1, según si n es par o impar. Este resultado está de acuerdo con la fórmula de dimensión anterior, y equivale a la dimensionalidad de una representación simétrica de SU(3) , [17] el grupo de degeneración relevante.
La notación de un oscilador armónico se puede extender a una red unidimensional de muchas partículas. Consideremos una cadena armónica mecanocuántica unidimensional de N átomos idénticos. Este es el modelo mecanocuántico más simple de una red, y veremos cómo surgen los fonones a partir de él. El formalismo que desarrollaremos para este modelo es fácilmente generalizable a dos y tres dimensiones.
Al igual que en la sección anterior, denotamos las posiciones de las masas por x 1 , x 2 , … , medidas desde sus posiciones de equilibrio (es decir, x i = 0 si la partícula i está en su posición de equilibrio). En dos o más dimensiones, las x i son cantidades vectoriales. El hamiltoniano para este sistema es
donde m es la masa (supuestamente uniforme) de cada átomo, y x i y p i son los operadores de posición y momento para el i ésimo átomo y la suma se realiza sobre los vecinos más cercanos (nn). Sin embargo, es habitual reescribir el hamiltoniano en términos de los modos normales del vector de onda en lugar de en términos de las coordenadas de la partícula para poder trabajar en el espacio de Fourier más conveniente .
Introducimos, entonces, un conjunto de N "coordenadas normales" Q k , definidas como las transformadas de Fourier discretas de las x s, y N "momentos conjugados" Π definidos como las transformadas de Fourier de las p s,
La cantidad k n será el número de onda del fonón, es decir 2 π dividido por la longitud de onda . Toma valores cuantizados, porque el número de átomos es finito.
Esto conserva las relaciones de conmutación deseadas tanto en el espacio real como en el espacio vectorial de onda.
A partir del resultado general es fácil demostrar, mediante trigonometría elemental, que el término de energía potencial es donde
El hamiltoniano puede escribirse en el espacio vectorial de onda como
Obsérvese que los acoplamientos entre las variables de posición se han transformado; si los Q y los Π fueran hermíticos (lo cual no son), el hamiltoniano transformado describiría N osciladores armónicos desacoplados .
La forma de cuantificación depende de la elección de las condiciones de contorno; para simplificar, imponemos condiciones de contorno periódicas , definiendo el átomo ( N +1) -ésimo como equivalente al primer átomo. Físicamente, esto corresponde a unir la cadena en sus extremos. La cuantificación resultante es
El límite superior de n proviene de la longitud de onda mínima, que es el doble del espaciado reticular a , como se explicó anteriormente.
Los valores propios del oscilador armónico o niveles de energía para el modo ω k son
Si ignoramos la energía del punto cero , entonces los niveles están espaciados uniformemente en
Por lo tanto, se debe suministrar una cantidad exacta de energía ħω a la red del oscilador armónico para impulsarlo al siguiente nivel de energía. En analogía con el caso del fotón , cuando el campo electromagnético está cuantizado, el cuanto de energía vibracional se denomina fonón .
Todos los sistemas cuánticos presentan propiedades ondulatorias y corpusculares. Las propiedades corpusculares del fonón se entienden mejor utilizando los métodos de segunda cuantificación y las técnicas de operadores descritas en otros trabajos. [18]
En el límite continuo , a → 0 , N → ∞ , mientras que Na se mantiene fijo. Las coordenadas canónicas Q k se transforman en los modos de momento desacoplados de un campo escalar, , mientras que el índice de ubicación i ( no la variable dinámica de desplazamiento ) se convierte en el argumento del parámetro x del campo escalar, .