(Numerical) Boundary Value Problem for Schrodinger's Equation

Click For Summary
The discussion focuses on solving a one-dimensional Schrödinger equation for a rectangular potential barrier using boundary value problem techniques. The standard method involves matching wavefunctions at the barrier boundaries, but the continuous spectrum complicates numerical solutions. A paper suggests that by selecting arbitrary points outside the barrier, specific boundary conditions can be established for numerical analysis. Questions arise regarding the derivation of these boundary conditions and their applicability in two dimensions, with some success reported in simulations. The conversation concludes with a reference to relevant literature for further exploration of the topic.
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.
 

Similar threads

Replies
5
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
1K
Replies
7
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
29
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 8 ·
Replies
8
Views
2K