Closed-form solutions to the wave equation

In summary, the solution for the massless wave equation corresponding to Gaussian initial conditions is given by: \psi(t,x) = \frac{1}{a^{3/2} \pi^{3/4} }\frac{1}{2r} \left( F(t+r) - F(t-r) \right)
  • #1
schieghoven
85
1
Hi all,

I'm interested to find a solution to the wave equation corresponding to
Gaussian initial conditions
[tex] \psi(0,x) = e^{-x^2/2} [/tex]
A solution which satisfies these initial conditions is (up to some constant factor)
[tex] \psi(t,x) = \int \frac{d^3k}{(2\pi)^3} e^{-k^2/2 + i(k \cdot x - \omega t)} [/tex]
where \omega = |k|. If we use spherical coordinates ( so that k.x = |k| r cos \theta )
then the angular dependence can be integrated out ... I get
[tex] \psi(t,x) = \int_{0}^{\inf} \frac{dk}{(2\pi)^2} \frac{2k \sin(kr)}{r} e^{-k^2/2 - \omega t} [/tex]
but can't get any further. Any ideas? Maybe this isn't the right direction to go anyway.

The reason I'm looking at this is because I'm interested to see the viability of
Hermite functions as a basis for state space in QM: the above function e^-x^2/2 is the
zeroth Hermite function. Hermite functions form a countable orthonormal basis, they
can be defined to have unit normalization, and all of their moments are finite. None
of these are true for the usual plane wave basis... so I feel the Hermite functions provide
a more concrete realisation of the Hilbert space of states. Of course, they won't be much
use unless I get them in closed form for t>0.

Thanks,
Dave
 
Physics news on Phys.org
  • #2
I have an idea, but I'm not sure whether it will be helpful.. I mean, if you rewrite sin(kr) as (exp(ikr)-exp(-ikr))/2i, then you have something on the form k*exp(f(k)). Then maybe you can use integration by parts to integrate? Maybe it won't work...
 
  • #3
The exponent of the initial solution is a quadratic equation. Complete the square, do a change of variables, done.
 
  • #4
schieghoven said:
Hi all,

I'm interested to find a solution to the wave equation corresponding to
Gaussian initial conditions
[tex] \psi(0,x) = e^{-x^2/2} [/tex]

Do a plane wave decomposition from the x to the p domain. Then presume positive
energy and thus :

[tex] \omega = \sqrt{k^2+m^2}[/tex]

This gives you the time-evolution of the independent plane waves. The last step is
the 3d inverse Fourier from the (t,p) domain to the (t,x) domain. The result comes
down to a convolution between the Gaussian and a second order Bessel K function
with imaginary argument. (A Hankel function)


Regards, Hans.
 
  • #5
Thanks for the three responses, they were a lot of help, and in the end helped me piece together an answer for the massless case... it's pretty epic.

The problem:
Massless propagation of free scalar field with Gaussian initial conditions
[tex]
\Box^2 \psi(t,x) = 0
[/tex]
[tex]
\psi(0,x) = \frac{1}{ a^{3/2} \pi^{3/4} }\exp \left( -\frac{x^2}{2a^2} \right)
[/tex]​
where x is shorthand for the three space dimensions (x,y,z). 'a' is a constant.

The solution is
[tex]
\psi(t,x) = \frac{a^{1/2}}{ \pi^{3/4} 2r } \left( F(t+r) - F(t-r) \right)
[/tex]​
where r = \sqrt(x^2 + y^2 + z^2) and F is given by
[tex]
F(s) = s \exp \left( -\frac{s^2}{2a^2} \right) \left( 1- i \text{Erfi} \frac{s}{\sqrt{2}a} \right)
[/tex]​
The error functions arise because, as pointed out in the responses, the second integral in the original post is of the form k*exp(quadratic), and this can be done by substitution and it gives rise to an error function. Erfi is notation used by Mathematica; it's basically just Erf rotated by 90 degrees in the complex plane.

I've attached two plots which show respectively the real and imaginary parts of the solution \psi(t,r), plotted on the (t,r) plane. Gaussian initial conditions can be seen along the axis t=0. It's kind of nifty, the light-like propagation is clearly visible in each plot, as the packets propagate along the lightcone t=+/-r.

If there's any interest to see the details, I guess I can try write them up somehow.

Thanks again for the tips,

Dave
 

Attachments

  • GaussianWavepacketRealPart.pdf
    716.9 KB · Views: 336
  • GaussianWavepacketImagPart.pdf
    736.8 KB · Views: 222
  • #6
Sorry, that should be
[tex]
\psi(t,x) = \frac{1}{a^{3/2} \pi^{3/4} } \frac{1}{2r} \left( F(t+r) - F(t-r) \right)
[/tex]​
... the answer in the previous post is out by a factor of a^2. For which I earn the 'wrong dimension' hat.
 
  • #7
schieghoven said:
that should be
[tex]
\psi(t,x) = \frac{1}{a^{3/2} \pi^{3/4} } \frac{1}{2r} \left( F(t+r) - F(t-r) \right)
[/tex]​
.

At first glance that looks to be what it should be. With the two terms coming
from the sine used in the radial Fourier transform. Nice images BTW.Regards, Hans
 

1. What is a closed-form solution to the wave equation?

A closed-form solution to the wave equation is an analytical solution that expresses the displacement of a wave at any point in space and time, without the need for numerical calculations. It is obtained by solving the wave equation using mathematical techniques such as separation of variables or Fourier transforms.

2. Why are closed-form solutions important in studying waves?

Closed-form solutions provide a complete and exact description of the behavior of a wave, which can be useful in understanding its properties and predicting its behavior in different scenarios. They also allow for easier and more efficient calculations compared to numerical methods.

3. What are the limitations of closed-form solutions to the wave equation?

Closed-form solutions are only applicable to simple wave systems and cannot account for complicated factors such as nonlinearities, damping, and boundary conditions. They also assume idealized conditions and may not accurately represent real-world phenomena.

4. How do closed-form solutions differ from numerical solutions?

Closed-form solutions provide an exact mathematical expression for the wave behavior, while numerical solutions involve approximating the solution through a series of calculations. Closed-form solutions are generally faster and more accurate, but may not be feasible for complex systems.

5. Can closed-form solutions be applied to all types of waves?

No, closed-form solutions are only applicable to linear waves, which are waves that follow the principle of superposition. Nonlinear waves, such as shock waves and solitons, require different mathematical techniques to find solutions.

Similar threads

Replies
9
Views
721
Replies
17
Views
1K
  • Quantum Physics
Replies
21
Views
1K
Replies
1
Views
584
  • Quantum Physics
Replies
24
Views
590
  • Quantum Physics
Replies
19
Views
1K
Replies
29
Views
3K
  • Quantum Physics
Replies
31
Views
3K
  • Quantum Physics
Replies
10
Views
1K
Replies
9
Views
1K
Back
Top