¿Por qué ambas formas equivalentes de esta función delta no dan la respuesta correcta?

Estoy un poco confundido sobre un problema básico que involucra una función delta de Dirac que se integra en una integral múltiple.

El problema original es encontrar la distribución de probabilidad en posición-momento ( z , pag ) espacio de una pelota que rebota hacia arriba y hacia abajo, dado que la pelota alcanza una altura máxima de z = h . Ahora, un argumento básico de conservación de la energía nos da que pag ( z ) = ± metro 2 gramo ( h z ) dónde metro es masa y gramo es la aceleración gravitacional. También sabemos que la densidad de probabilidad de encontrar la pelota en z es inversamente proporcional a la velocidad de la pelota en z , a saber pag ( z ) / metro . Entonces, la distribución de probabilidad es

PAG ( z , pag ) = C 2 gramo ( h z ) [ d ( pag metro 2 gramo ( h z ) ) + d ( pag + metro 2 gramo ( h z ) ) ]

dónde C es un factor de normalización. Encontrar C , podemos simplemente integrar PAG encima z y pag y establecer el resultado igual a 1 , usando las funciones delta para hacer el pag integral e integrador z de 0 a h , y obtenemos la respuesta correcta. Mi pregunta es, ¿por qué no funciona esta estrategia si escribimos

PAG ( z , pag ) = D | pag | d ( z h + pag 2 / ( 2 metro 2 gramo ) ) ?
En este caso, cuando intentamos encontrar la normalización D , si intentamos deshacernos de la función delta integrando sobre z y luego integramos sobre el rango apropiado de pag , fallamos porque el 1 / | pag | la integral diverge, mientras que la 1 / 2 gramo ( h z ) integral converge. ¿Que me estoy perdiendo aqui?

Respuestas (2)

Para el segundo caso, desea dividir por la fuerza, no por la velocidad.

Básicamente, está calculando qué fracción del tiempo pasa en un punto particular en el espacio de fase.

Sin embargo, lo que tienes es una densidad de probabilidad. Entonces PAG ( z , pag ) es algo que multiplicas por un volumen en el espacio de fase para obtener una probabilidad.

Una forma correcta de obtenerlo sería considerar un período completo (por ejemplo, el tiempo de arriba hacia arriba o de abajo hacia abajo o desde cualquier punto del espacio de fase hacia sí mismo).

Un delta de Dirac tiene unidades. Tiene unidades de uno sobre el argumento, así que desde d ( X ) d X = 1 entonces d ( X ) tiene unidades de distancia inversa y d ( pag ) d pag = 1 entonces d ( pag ) tiene unidades de impulso inverso.

Cuando integras sobre el impulso, obtienes una densidad de probabilidad para la posición. Esto tiene unidades que multiplicas por una distancia para obtener una probabilidad de estar en ese rango de distancia.

Estas densidades de probabilidad sobre posiciones son aquellas donde es inversamente proporcional a la velocidad. Entonces, de hecho, después de integrar el impulso, tiene una probabilidad por unidad de longitud y las probabilidades por unidad de longitud son inversamente proporcionales a la velocidad. Así que debería funcionar.

Sin embargo, si integra sobre la posición, entonces tiene una densidad de probabilidad para el impulso. Esto es algo que se multiplica por un rango de impulsos para obtener una probabilidad de tener un impulso en ese rango. No existe una regla de que este tipo de densidad sea inversamente proporcional a la velocidad. Y no lo es, es inversamente proporcional a la fuerza.

Esto es similar a cómo la densidad de energía es diferente para la frecuencia y la longitud de onda porque cada uno es algo que se multiplica por un rango de cosas para obtener la cantidad real de cosas. Vea mi respuesta en https://physics.stackexchange.com/a/165256 si necesita más detalles sobre las densidades en diferentes dominios unidos.

Repasemos de dónde viene la regla de la inversamente proporcional a la velocidad. Como una simplificación durante un intervalo de tiempo Δ t podrías viajar una distancia Δ z = v Δ t Y realmente la probabilidad que estás calculando es Δ t / T dónde T es ese período que mencioné.

Así que la probabilidad de que estés en una región de ancho Δ z es proporcional a la cantidad de tiempo Δ t pasas ahí. Por lo tanto, inversamente proporcional a la velocidad, de hecho es Δ z / ( v T ) . Entonces la densidad espacial es 1 / ( v T ) . Y esta regla no se aplica si pregunta cuánto tiempo pasa algo en una región de tamaño Δ pag . Entonces, no divida por velocidad si desea una densidad de probabilidad de estar dentro de un rango de impulsos de ancho Δ pag .

Para hacer la densidad de probabilidad para los intervalos de momento, use Δ pag = F Δ t entonces la probabilidad es Δ t / T = Δ pag / ( F T ) por lo que la densidad de probabilidad para los intervalos de cantidad de movimiento es 1 / ( F T ) .

Entonces, una densidad es como el delta de Dirac, es algo que integras para obtener la probabilidad, pero debe ser por intervalo de distancia o por intervalo de impulso para que sea algo que puedas integrar.

Y la versión de impulso debería darle una pausa, ya que generalmente necesita tener en cuenta correctamente la fuerza cuando gira en el suelo. Así que el exceso de conteo es un tema a considerar.

TL; DR: la sustitución dentro de la función delta produce un factor jacobiano

(1) d ( F ( v ) )   =   v ( 0 ) , F ( v ( 0 ) ) = 0 1 | F ( v ( 0 ) ) | d ( v v ( 0 ) ) .

Aquí la suma es sobre los ceros v ( 0 ) de la función F ( v ) .

Consideremos por simplicidad la velocidad v en lugar de impulso pag = metro v . Entonces la conservación de la energía

(2) 1 2 v 2 + gramo z   =   gramo h , 0 z h , | v |     2 gramo h ,

produce una parábola en el ( v , z ) avión. Si definimos una función

(3) F ( v )   :=   z h + v 2 2 gramo ,

entonces ec. (1) se convierte

(4) d ( z h + v 2 2 gramo )   =   gramo | v | ± d ( v ± 2 gramo ( h z ) ) .

Nótese en particular que el factor 1 | v | en el lado derecho de la ec. (4) no aparece en el lado izquierdo.