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

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
1 reply · 2K views
Foracle
Messages
29
Reaction score
8
TL;DR
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: 177
Last edited:
Reply
  • Like
Likes   Reactions: BvU
Physics news on Phys.org
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.