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

Click For Summary
SUMMARY

This discussion focuses on solving a relativistic infinite square well problem with an oscillating wall using the Klein-Gordon equation in Mathematica. The transformation of spatial coordinates from ##x## to ##y = \frac{x}{L[t]}## is employed to simplify the boundary conditions. The issue of the wavefunction ##\psi## blowing up at high frequencies (##\omega = 1000##) is attributed to the wall's velocity exceeding the speed of light, necessitating a modification of the wall function to ##L(t) = \frac{1}{\omega} (2+\sin(\omega t))## to avoid this problem.

PREREQUISITES
  • Understanding of the Klein-Gordon equation
  • Familiarity with numerical methods in Mathematica
  • Knowledge of relativistic boundary conditions
  • Concept of wavefunction behavior in quantum mechanics
NEXT STEPS
  • Explore advanced numerical methods in Mathematica for solving differential equations
  • Study the implications of relativistic effects on wavefunctions in quantum mechanics
  • Investigate particle production phenomena in oscillating potentials
  • Learn about non-static boundary conditions in quantum systems
USEFUL FOR

Physicists, mathematicians, and computational scientists interested in quantum mechanics, particularly those working with relativistic equations and numerical simulations in Mathematica.

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: 150
Last edited:
  • 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.
 

Similar threads

  • · Replies 3 ·
Replies
3
Views
599
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 9 ·
Replies
9
Views
2K
  • · Replies 45 ·
2
Replies
45
Views
4K
  • · Replies 18 ·
Replies
18
Views
3K
  • · Replies 28 ·
Replies
28
Views
3K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 10 ·
Replies
10
Views
4K
  • · Replies 31 ·
2
Replies
31
Views
4K