(Numerical) Boundary Value Problem for Schrodinger's Equation

cc94
Messages
19
Reaction score
2

Homework Statement


Suppose we have the standard rectangular potential barrier in 1D, with
$$
V =
\left\{
\!
\begin{aligned}
0 & \,\text{ if } x<0, x>d\\
V_0 & \,\text{ if } x>0,x<d\\
\end{aligned}
\right.
$$

The standard approach to solve for tunneling through the barrier is to match the wavefunctions at the boundaries of the barrier, and solve for the ratio of amplitudes of the outgoing and incoming waves. My issue is that this problem can't be solved numerically, since it is open on the ends, so the spectrum is continuous. Or so I thought, until I found a paper that says that one can turn this problem into a boundary value problem. If you pick some arbitrary $$x_a, x_b$$ to the left and right of the barrier, respectively, then one can use the following equations as boundary conditions to numerically solve the ODE. My first question is, how were these equations derived? And then if it's possible to answer, does this method work in 2D? From simulating some 2D potentials with these boundary conditions, it looks like it works, but I need to know more specifics.

Homework Equations


$$
\left\{
\!
\begin{aligned}
\frac{d\psi}{dx}\Big|_{x=x_a} &+ ik\psi(x_a) = 2ikA\text{exp}(ikx_a),\\
\frac{d\psi}{dx}\Big|_{x=x_b} &- ik\psi(x_b) = 0
\end{aligned}
\right.
$$
where A is amplitude of inc. wave.

The Attempt at a Solution


[/B]
To derive those equations, I'm assuming you follow the same piecewise procedure to match wavefunctions, with the potential being V=0 on both sides. But that doesn't give an answer. Also I don't get why, for example, the first condition only depends on the amplitude of the incoming wave and not the amplitude of the reflected wave, R.
 
Physics news on Phys.org
cc94 said:
Or so I thought, until I found a paper that says that one can turn this problem into a boundary value problem.
Reference, please.
 
DrClaude said:
Reference, please.

Sure, google "Quantum scattering theory and stealth finite element analysis", and it's the second result for me.
 
I didn't solve my first question, but I found an approach for the 2D case in the Quantum Transmitting Boundary Model. I'll mark this solved.
 
Hello everyone, I’m considering a point charge q that oscillates harmonically about the origin along the z-axis, e.g. $$z_{q}(t)= A\sin(wt)$$ In a strongly simplified / quasi-instantaneous approximation I ignore retardation and take the electric field at the position ##r=(x,y,z)## simply to be the “Coulomb field at the charge’s instantaneous position”: $$E(r,t)=\frac{q}{4\pi\varepsilon_{0}}\frac{r-r_{q}(t)}{||r-r_{q}(t)||^{3}}$$ with $$r_{q}(t)=(0,0,z_{q}(t))$$ (I’m aware this isn’t...
Hi, I had an exam and I completely messed up a problem. Especially one part which was necessary for the rest of the problem. Basically, I have a wormhole metric: $$(ds)^2 = -(dt)^2 + (dr)^2 + (r^2 + b^2)( (d\theta)^2 + sin^2 \theta (d\phi)^2 )$$ Where ##b=1## with an orbit only in the equatorial plane. We also know from the question that the orbit must satisfy this relationship: $$\varepsilon = \frac{1}{2} (\frac{dr}{d\tau})^2 + V_{eff}(r)$$ Ultimately, I was tasked to find the initial...
Back
Top