Energía de pozo potencial infinito para la función de onda constante por partes

Estoy tratando de calcular el valor esperado de la energía para un cierto estado en un pozo de potencial infinito , pero obtengo respuestas contradictorias.

El pozo tiene potencial

V ( X ) = { 0 : 0 < X < L :  en otra parte
que tiene estados propios ϕ norte ( X ) = 2 L pecado ( norte π X L ) con las energías correspondientes mi norte = norte 2 π 2 2 2 metro L 2 .

Ahora, considere el estado

ψ ( X ) = { 0 : 0 < X < L / 2 2 L : L / 2 X < L 0 :  en otra parte
queremos calcular H para este estado. Una forma de hacer esto es simplemente usando la definición: H = 0 L ψ ( X ) H ψ ( X ) d X . Aunque el problema con esto es que ψ ( X ) es constante por partes, y por lo tanto esto le dará 0 . La otra opción es ampliar ψ en términos de la base propia de energía calculando los coeficientes C norte = ϕ norte | ψ , y obteniendo H = norte = 1 | C norte | 2 mi norte . Como mi norte > 0 para cada norte , esta cantidad será estrictamente mayor que 0 y por lo tanto diferirá de la respuesta anterior.

¿Cuál es la discrepancia? Seguramente tiene que ver con la discontinuidad en X = L / 2 , pero no puedo entender cómo lidiar con eso.

Tenga en cuenta que su constante por partes ψ ( X ) no es continua y, por lo tanto, no es una función de onda físicamente disponible. Intenta repetir el cálculo usando C 1 funciones que convergen a ψ y obtendrá contribuciones ilimitadamente crecientes del "pegamento" en 0 y L.
De manera más general, representar funciones discontinuas utilizando series de Fourier no es un asunto trivial y conduce al fenómeno de Gibbs .

Respuestas (1)

El valor esperado

(1) ψ | V | ψ   =   0

del operador de energía potencial V es de hecho cero, pero el valor esperado

(2) ψ | k | ψ   = 2 2 metro R d X   | ψ ( X ) | 2   =   +

del operador de energía cinética k es en realidad infinito para la función de onda

(3) ψ ( X )   =   2 L ( θ ( X L 2 ) θ ( X L ) ) , X R .

Aquí θ es la función escalón de Heaviside .

El operador de energía cinética k := 2 2 metro d 2 d X 2 es un ejemplo de un operador ilimitado , que solo tiene sentido en su dominio D k H dentro del espacio de Hilbert H := L 2 ( R ) de funciones integrables de Lebesgue cuadradas. En particular, es un problema matemático no trivial cómo aplicar el operador diferencial k a la función de onda no diferenciable (3).

El resultado infinito (2) se puede ver (en el nivel físico de rigor) en al menos tres formas (ordenadas con el cálculo computacionalmente más simple primero):

  1. Conecte la función de paso de Heaviside en la ec. (2) para obtener una integral sobre el cuadrado de un par de funciones delta de Dirac situadas en X = L 2 y X = L . Esto es estrictamente hablando matemáticamente mal definido. Físicamente, tiene sentido asignar a la integral el valor infinito, cf. esta publicación Phys.SE.

  2. Calcular las superposiciones C norte = ϕ norte | ψ , y mostrar que la suma

    (4) ψ | H | ψ = norte = 1 | C norte | 2 mi norte   =   +
    diverge Esta conclusión infinita parece físicamente sólida, ya que todos los términos de la serie (4) son no negativos.

  3. Por regularización , como sugiere Emilio Pisanty en un comentario. Definir una función de onda regularizada ψ ε C 1 ( R ) de tal manera que (i) converge ψ ε ψ para ε 0 + , (ii) el valor esperado ψ ε | k | ψ ε es fácil de calcular y finito para ε > 0 . Muestra esa ψ ε | k | ψ ε + diverge para ε 0 + .

Gracias, mi problema con el enfoque (1) es, como dices, matemáticamente no tiene sentido para mí. Estaba pensando en términos de integrales de Lebesgue, donde ψ es 0 casi en todas partes y por lo tanto su integral (o la integral de su cuadrado) es por lo tanto 0 (¿hay algo malo con este proceso de pensamiento?). El enfoque (2) es el más sencillo (en mi opinión) ya que se puede hacer con fuerza bruta.
Este requisito de que las funciones de ondas físicas sean continuamente diferenciables me resulta confuso. Podría ser solo mi falta de comprensión de QM de dimensión infinita, pero tome cualquier función de onda 1D ψ (hazlo tan bonito como quieras) y elige X de modo que ψ ( X ) 0 . Ahora mida si la partícula está a la izquierda o a la derecha de X . Supongamos que encuentra la partícula a la derecha de X , ¿cuál es el estado después de la medición? Al estar acostumbrado a dim. finitos, asumo que 'proyecto ψ ' (colocar pag s i ( X ) = 0 para X < X ) y normalizar. Esta función no será continua, pero no puedo ver el problema.
El problema es que la función de onda ψ en general tiene que pertenecer al dominio D A del operador autoadjunto A si se quiere medir el valor esperado ψ | A | ψ .