Función de Green para la ecuación no homogénea de Klein-Gordon

Estoy tratando de resolver la enorme ecuación de Klein-Gordon en el viejo espacio-tiempo de Minkowski:

( + metro 2 ) ϕ = ρ ( t , X )
dónde = m m = t 2 2 . Entonces uno puede usar el enfoque de la función de Green para encontrar la solución fundamental de la forma
( + metro 2 ) GRAMO metro = d ( X m X m )
dónde GRAMO metro es el conocido propagador de Klein-Gordon. Entonces se obtiene la solución ϕ en el espacio de posición, dado como la solución familiar
ϕ ( X m ) = d 4 X GRAMO metro ( X m , X m ) ρ ( t , X ) ( )
Estaba perfectamente feliz con esto hasta que necesité implementar un ρ y realiza las integrales. Mi mejor apuesta hasta ahora ha sido usar la representación de la función de Bessel que encontré (aquí: http://functions.wolfram.com/Bessel-TypeFunctions/BesselJ/31/02/ ) (que supuse que tiene una generalización) de la forma:
GRAMO metro ( t , t , X , X ) = θ ( t t ) 2 π d ( ( t t ) 2 | X X | 2 ) metro 2 π θ ( t t | X X | ) j 1 ( metro ( t t ) 2 | X X | 2 ) metro ( t t ) 2 | X X | 2 )
Si bien esta es una buena representación de forma cerrada, todavía tengo dificultades para evaluar la integral. ( ) . He buscado durante bastante tiempo en varios lugares ejemplos explícitos de cómo calcular la integral, y hasta ahora he encontrado muy pocos. Mathematica (mi programa computacional de elección) realmente desdeña estas funciones de Heaviside en las integrales y ofrece poca orientación. El único caso que puedo hacer hasta ahora es metro 0 .

Pregunta: Usando la representación de GRAMO metro dado (u otro mejor), ¿cómo se puede hacer realmente para calcular ( ) ? ¿Alguien tiene una referencia en la que se calcule algún ejemplo explícito donde ρ va más allá de un simple d -¿función? Incluso algo como ρ = ρ ( r , θ ) o ρ = ρ ( r ) seria de gran ayuda

¡Gracias!

Por favor, publique el específico ρ ( X , t ) , hay que tratar caso por caso.
Estoy buscando una estrategia genérica aquí, lo que puede ser una ilusión, supongo. Los detalles asociados con ρ no son tan importantes (ya que tengo muchos casos con los que buscaría jugar eventualmente). Quiero saber si la representación que tengo allí es buena para los cálculos y cómo se haría. Idealmente, esperaba que alguien supiera una referencia en la que se calculan algunos ejemplos explícitos. Incluso algo como ρ = a r 2 + b r + C Sería útil, ya que puedo ver esto y espero generalizar. Tal como está, realmente no sé cómo abordarlo.
Si las funciones de Heaviside son el problema, ¿puede intentar resolver el problema sin ellas e implementar una distinción de casos a mano? Los Heaviside parecen implementar la causalidad...
Gracias, @Neuneck, definitivamente vale la pena intentarlo. Sin embargo, con respecto a las integraciones, estoy luchando para evaluar incluso
d 4 X j 1 ( metro ( t t ) 2 | X X | 2 metro ( t t ) 2 | X X | 2
(Eso es cuando ρ es constante). Encontré esta ( fh-jena.de/~rsh/Forschung/Stoer/besint.pdf ) tabla de integrales de funciones de Bessel pero esta no aparece allí. ¿Hay alguna razón para suponer que la integral existe en forma elemental? Supongo que no.
¿Tiene una referencia de cómo se deriva la función de green? Me refiero a la fórmula de wolframio

Respuestas (1)

Soy consciente de que esta es una pregunta antigua y puede considerarse algo obsoleta, pero permítanme responderla de todos modos, aunque solo sea por el bien de la exhaustividad.

La representación del espacio de posición de la función de Klein-Gordon Green (propagador) claramente parece intimidante. El truco consiste en hacer el cálculo en el espacio de cantidad de movimiento, donde el propagador es solo una función racional. Antes de hacer esto, permítanme señalar que en el caso sin masa, metro = 0 , y para una fuente estática, ρ = ρ ( X ) , en realidad se está resolviendo la ecuación de Poisson. Si la fuente es radialmente simétrica, ρ = ρ ( r ) , como se sugiere en la pregunta, la solución es el potencial de Coulomb, ϕ = ϕ ( r ) 1 / r . Teniendo en cuenta la masa que no desaparece, se obtiene un potencial de Yukawa, ϕ Exp ( metro r ) / r .

Esto se puede mostrar explícitamente en términos de integrales de Fourier. Primero, transforma el campo,

ϕ ( k ) = d 4 X mi i k X ϕ ( X ) ,

y de manera similar la fuente, ρ ρ ( k ) . La solución del espacio de cantidad de movimiento de la ecuación de KG es entonces

ϕ ( k ) = ρ ( k ) k 2 metro 2 + i ϵ ,

con k 2 = k 0 2 k 2 y la causa i ϵ -prescripción. (Cambie apropiadamente para soluciones retrasadas/avanzadas). Luego transforme de nuevo al espacio de posición,

ϕ ( X ) = d 4 k ( 2 π ) 4 mi i k X ρ ( k ) k 2 metro 2 + i ϵ . ( )

Como ejemplo, tome una fuente gaussiana, ρ ( r ) = ρ 0 Exp ( α r 2 ) , con 'normalización' ρ 0 = ( α / π ) 3 / 2 . Su transformada de Fourier es ρ ( k ) = 2 π i d ( k 0 ) Exp ( k 2 / 4 α ) . los k 0 -integral en ( ) es por lo tanto trivial, y el d 3 k la integración se puede hacer con la técnica habitual de residuos escribiendo k 2 + metro 2 ( k + i metro ) ( k i metro ) . El resultado es

ϕ ( r ) = mi metro 2 / 4 α mi metro r 4 π r .

En el límite de la fuente puntual ( α ) volvemos a obtener el potencial estándar de Yukawa.

Para fuentes dependientes del tiempo habrá transferencia de energía (no d ( k 0 ) ), y la integral ( ) normalmente será más difícil. Por lo general, uno puede hacer el k 0 -integración vía residuos y el(los) restante(s) usando fase estacionaria como por ejemplo en el Cap. 2.1 de Peskin y Schroeder.

Gracias por tomarse el tiempo para escribir esto. He visto estos enfoques antes, pero de hecho me preocupaba que no funcionaran con términos dependientes del tiempo y el ángulo. Cogeré a Peskin ya Schroeder de la biblioteca y lo intentaré; gracias por la referencia del capitulo. Seguramente, aunque las soluciones en Fourier y el espacio de posición deben ser equivalentes bajo una transformación, ¿implica esto que tenemos algunas representaciones de forma cerrada para las integrales de Bessel? ¡Interesante! Puede ser una vía de investigación interesante por derecho propio =)