Barrera triangular en pozo de potencial infinito

Supongamos que estoy buscando resolver la función de onda para el siguiente potencial 1D:

(1) tu ( X ) = { V 0 a | X | a por | X | < a por | X | > a
Dado que nuestro potencial es simétrico, tenemos soluciones pares e impares, por lo que podemos resolver el sistema para X 0 y luego construir la solución completa para cualquiera ψ mi v mi norte o ψ o d d . Para la región dentro del pozo, podemos proyectar la ecuación de Schrödinger independiente del tiempo. en una forma que es la ecuación diferencial de Airy (tenga en cuenta que es un potencial lineal), a saber
(2) d 2 ψ d z 2 z ψ = 0 ,
dónde
(3) z q ( 1 η ) ,

dónde

(4) η = X a 0 ,       q = 2 metro V 0 a 0 2 a ,       a 0 a a V 0 mi .

Por lo tanto podemos expresar ψ como

(5) ψ ( z ) = C 1 A i ( z ) + C 2 B i ( z ) .

Nuestras condiciones de contorno nos dicen:

(6) ψ ( ζ a ) = 0
(7) ψ ( ζ 0 ) = 0 por ψ o d d
(8) ψ ( ζ 0 ) = 0 por ψ mi v mi norte .

donde estoy denotando ζ 0 = z | X = 0 y ζ a = z | X = a .

Pregunta : ¿Es este sistema exactamente solucionable?

Por lo general, solo tenemos una condición límite cero en este tipo de problemas que nos permite cuantificar los niveles de energía a través de los ceros de la función de Airy. Pero aquí debemos satisfacer simultáneamente dos condiciones de frontera cero. Para aclarar, aquí está nuestro sistema de ecuaciones lineales.

ψ o d d :

(9) C 1 A i ( ζ 0 ) + C 2 B i ( ζ 0 ) = 0
(10) C 1 A i ( ζ a ) + C 2 B i ( ζ a ) = 0

ψ mi v mi norte :

(11) C 1 A i ( ζ 0 ) + C 2 B i ( ζ 0 ) = 0 ,
(12) C 1 A i ( ζ a ) + C 2 B i ( ζ a ) = 0.

Mi idea propuesta sería resolver por el determinante cero de las ecuaciones lineales homogéneas, pero esto me obligaría a una solución numérica. Me gustaría encontrar las energías en términos de los ceros de las funciones de Airy. ¿Ideas?

Confío en que conoces el artículo Ecuación de Airy en Wikipedia. El punto X = a Donde quieras ψ = 0 , cae en z = q ( 1 a / a 0 ) . te dara la razon C 2 / C 1 . Pero, atención que como X disminuye, z aumenta Puedes ver cómo se ve la región. z > 0 en los gráficos del artículo. Te interesará requerir también eso ψ o su derivada, sea cero en X = 0 es decir z = q , que está a la izquierda de z = q ( 1 a / a 0 ) . No parece que en esta región obtendrá otro cero o un máximo (o mínimo) para ψ ( z ) .(Continúo)
En general, un sistema de 2 ecuaciones homogéneas. con 2 variables, puede ser incompatible. Lo que pudo salvar la situación fue el valor de mi , que determina el valor de a 0 . Pero, de nuevo veo un problema en el hecho de que los valores de z caer en la región z > 0 . Me parece que podrías haber obtenido soluciones de tu ecuación incluyendo la condición de simetría o antisimetría, si z tomó valores en la región z < 0 donde las funciones A i y B i oscilar. Pero con z > 0 Soy bastante escéptico.

Respuestas (2)

De hecho, no hay forma de evitar tener que resolver este problema numéricamente. Partiendo del método determinante publicado originalmente, quizás un enfoque alternativo aquí sería resolver el sistema usando solo la energía. Comenzamos con la forma completa de nuestros puntos de interés tomando como ejemplo las funciones de onda impares:

z ( 0 ) = X 0 a ( mi V 0 + 1 )
z ( a ) = X 0 a ( mi V 0 ) .
yo denoto α i y β i como el i los ceros de A i ( z ) y B i ( z ) respectivamente. Igualmente α i y β i los ceros de las derivadas. También X 0 [ 2 metro V 0 2 a ] 1 / 3 - Se editó la notación de la publicación original :(

En primer lugar, sabemos que mi > V metro i norte y por lo tanto mi > 0 . Por lo tanto, vemos que z ( 0 ) > z ( a ) . Ahora, el primer cero de cualquiera A i ( z ) o B i ( z ) debe ocurrir en z < 0 . Para un valor apropiado de mi , podemos ver que si los ceros de las funciones de Airy difieren en X 0 a , entonces nuestra condición se cumpliría para una sola energía. Trabajando con las funciones de Airy del primer tipo, podemos imponer la condición

α i = X 0 a ( mi V 0 )
o
mi = V 0 X 0 a α i .
El criterio sobre el segundo cero es entonces
α j = X 0 a + α i .
Podríamos recorrer numéricamente todos los ceros de la función de Airy y para aquellos con una separación de
α j α i = X 0 a
entonces tendríamos una energía correspondiente de
mi = V 0 X 0 a α i .
Aplicando el mismo concepto a las funciones de onda pares e impares tenemos:
ψ mi v mi norte ( z ) = C 1 A i ( z ) ; mi = V 0 X 0 a α i ; α j α i = X 0 a
ψ o d d ( z ) = C 2 A i ( z ) ; mi = V 0 X 0 a α i ; α j α i = X 0 a
ψ mi v mi norte ( z ) = C 3 B i ( z ) ; mi = V 0 X 0 a β i ; β j β i = X 0 a
ψ o d d ( z ) = C 4 B i ( z ) ; mi = V 0 X 0 a β i ; β j β i = X 0 a
donde la tercera columna denota la condición que deben cumplir nuestros ceros, C i es la constante de normalización apropiada para cada función de onda, y
z [ 2 metro V 0 2 a ] 1 / 3 ( a 0 X ) = X 0 ( a 0 X ) ; a 0 [ a a V 0 mi ] .
Al observar un gráfico de la función de Airy, podemos ver que en z 0 la ``longitud de onda'' es muy pequeña, por lo que debería haber un gran número de posibles ceros que satisfagan nuestros criterios. De la misma manera, en este rango de z , nuestras energías serán muy grandes y tenemos la apariencia del pozo cuadrado infinito. Gracias Sofia y Ali Moh por los útiles comentarios.

Lo que pides es imposible. Los valores propios de energía están dados por los ceros del determinante (como supuso) y no pueden relacionarse con los ceros de la función de aire.

A i ( ζ a ) B i ( ζ 0 ) A i ( ζ 0 ) B i ( ζ a ) = 0
Para la solución impar (y con derivadas en A i para la solución uniforme).

También tenga cuidado, su z no es adimensional como debería ser, por lo que su cambio de variables no es del todo correcto...

¡Vea los comentarios para ver los argumentos físicos que garantizan infinitas soluciones a esta condición!

Lamento que no hayas leído mis comentarios. A través del determinante el OP intenta encontrar las constantes C que permiten cumplir las condiciones de contorno. No es seguro si el determinante se puede hacer cero . Sin embargo, si el determinante puede hacerse cero, a partir del valor a 0 para el cual el determinante se vuelve cero, uno puede encontrar la energía mi .
Leí tu comentario pero no estoy de acuerdo con él. El determinante siempre se puede hacer cero, la energía se puede cambiar continuamente para satisfacer esta condición. Una variable con una condición. Pero mi garantía proviene principalmente de la física, ¡es imposible que no haya muchos valores propios de energía! porque este es un pozo de potencial con una altura infinita, por lo que para energías lo suficientemente altas, la función de onda será en realidad ciega a la protuberancia triangular y sabemos en ese caso que hay un número infinito de valores propios discretos de energía
Así que no puedo decir con certeza si estos valores propios serán comparables a V 0 (la altura de la protuberancia), pero puedo garantizar a partir del argumento físico anterior que existen al menos estados propios con mi >> V 0 .. Más aún, me parece muy poco probable que los valores propios de energía bajos no existan
Una forma equivalente es notar el comportamiento trigonométrico asintótico de las funciones aéreas, lo que valida la intuición física.
No hay necesidad de mi >> V 0 , me parece suficiente mi > V 0 . el problema es por mi < V 0 . Estás mirando el determinante y yo estoy mirando cómo se ven las funciones de onda en la región. z > 0 .
cual asintotica? El soporte de la función de onda es finito.
Entiendes que V ( X ) por | X | > a , ¿Correcto? así que cualquiera mi > 0 representa un estado ligado expresable en términos de la superposición lineal de funciones aéreas. Así que todo lo que digo es que EL DETERMINANTE ESTÁ GARANTIZADO PARA TENER SOLUCIONES, al menos para mi >> V 0 ... No entiendo qué tiene de especial mi < V 0 que dices que hay un problema? olvidarse de z , dibuje el potencial como una función de X y pensar en funciones de onda de estado ligado
Estás tratando de convencerme de que puede haber una solución para mi > V 0 allí donde es fácil. que tiene de especial mi < mi 0 es que en este caso es difícil. Y es difícil precisamente porque el comportamiento en la región z > 0 . Creo que me expliqué. De todos modos, dejo aquí la argumentación.
En cuanto al cambio de variables que se editaron en mi post, no creo que sean correctos. Basado en lo que tenía antes: d 2 d X 2 = d 2 d z 2 [ 2 metro V 0 2 a ] 2 / 3 , dónde z [ 2 metro V 0 2 a ] 1 / 3 ( a 0 X )
@Sofia lo siento en este momento, no puedo entender lo que dices. También en lo que respecta al comportamiento asintótico, está claro que quise decir a 0 , entonces, ¿qué significa "el soporte de la función de onda es finito"?