Infinite Square Well with an Oscillating Wall (Klein-Gordon Equation)

In summary, the conversation discusses the process of numerically solving a relativistic version of an infinite square well with an oscillating wall using the Klein-Gordon equation. The speaker shares that they transformed the spatial coordinate and used NDSolve to solve the system. However, when they plotted the result, the wavefunction blew up when the frequency of the oscillation was high. The speaker then questions the reason for this and whether it has any connection to particle production. After locating their error, it is revealed that the velocity of the wall was not allowed in the units being used. The solution to this was to adjust the value of L(t) to avoid this issue.
  • #1
Foracle
30
8
TL;DR Summary
Numerical solution of the wavefunction blows up when the frequency of oscillation of the wall is high. What does this mean?
I am trying to numerically solve (with Mathematica) a relativistic version of infinite square well with an oscillating wall using Klein-Gordon equation. Firstly, I transform my spatial coordinate ## x \to y = \frac{x}{L[t]} ## to make the wall look static (this transformation is used a lot in solving non-static boundary condition in the non-relativistic case), which brings Klein-Gordon equation to :
Input :
1644917229990.png

Output :
1644917258275.png

All constants have been set to 1
1644917266102.png

I tried to solve this system where ##L(t)=2+sin(1000 t)## using NDSolve :
1644917349008.png

Then I plot my result as a function of ##(y,t)## :
1644917373977.png

The wavefunction ##\psi## blows up. This doesn't happen when I tune the frequency down to 1, ##L(t)=2+sin(t)##
1644917399498.png

My question is why does ##\psi## blow up when the frequency of the oscillation is high and what does it mean? Does it have anything to do with particle production? Or did I just mess up my code?
 

Attachments

  • 1644917332177.png
    1644917332177.png
    7.1 KB · Views: 118
Last edited:
  • Like
Likes BvU
Physics news on Phys.org
  • #2
I have located my error. It was the fact that when I set ##\omega## to 1000, the velocity of the wall, ##\dot{L} = \omega cos(\omega t) \ge 1##, which is not allowed in the units that I am working with ##(c=1)##.
The solution to this is just to set ##L(t) = \frac{1}{\omega} (2+sin(\omega t))##, and this disastrous result doesn't happen anymore.
 

Similar threads

Back
Top