Derivation of Schrodinger equation for a system with position dependent effective mass How to derive the Schrodinger equation for a system with position dependent effective mass? For example, I encountered this equation when I first studied semiconductor hetero-structures. All the books that I referred take the equation as the starting point without at least giving an insight into the form of the Schrodinger equation which reads as
$$\big[-\frac{\hbar^2}{2}\nabla \frac{1}{m^*}\nabla  + U  \big]\Psi ~=~ E \Psi. $$
I feel that it has something to do with conserving the probability current and should use the continuity equation, but I am not sure.
 A: For derivation of the PDM Schrodinger equation see
K. Young, Phys. Rev. B 39, 13434–13441 (1989)
"Position-dependent effective mass for inhomogeneous semiconductors".
Abstract.:A systematic approach is adopted to extract an effective low-energy Hamiltonian for crystals with a slowly varying inhomogeneity, resolving several controversies. It is shown that the effective mass $m_R$ is, in general, position dependent and enters the kinetic energy operator as $ -\nabla({m_R-1})\nabla/2$. The advantage of using a basis set that exactly diagonalizes the Hamiltonian in the homogeneous limit is emphasized.
A: A Hamiltonian must be self-adjoint. The equation must also reduce to the familiar equation in the case of a constant mass. Now the form of the operator is already determined as the only simple self-adjoint generalization of the position-independent Schroedinger equation to the position dependent case.
If you specialize to 1 dimension, you get the Sturm-Liouville equation. At 
http://en.wikipedia.org/wiki/Sturm-Liouville_theory
you can find a discussion of its self-adjointness. Everything genaralizes to the PDE case.
A: The derivation is straightforward if you consider the source of the effective mass is a slowly varying hopping parameter on a tight-binding (lattice particle) model. Here you have a particle on a square lattice with a probability amplitude to go left, right, up and down, forward, backwards. The main physical requirement is Hermiticity, which in 1d can be used (with a phase choice on the wavefunction) to turn the phase everywhere real.
Once you do this, there is a real amplitude at site n to hop one square to the right r(n) and an amplitude to hop one square to the left, which by hermiticity and reality, must be r(n-1)--- it must be the complex conjugate of the amplitude to hop right from position n-1. So the amplitude equation is
$$ i{dC\over dt} = r(n-1) C_{n-1} - (r(n-1)+r(n))C_n + r(n) C_{n+1} $$
This is, when r is slowly varying, equivalent to the continuum equation found by Taylor expanding and keeping only the most relevant terms:
$$ i {d\psi \over dt} = {1\over 2} {\partial \over \partial x} (r(x) {\partial\over \partial x} \psi(x)) $$
As Feynman noted but never published (Dyson published this comment posthumously, in a paper in American Journal of Physics titles something like "Feynman's derivation of the Maxwell equations from Schrodinger equation"), Dirac's phase trick doesn't work in higher dimensions, because you can't fix all the phases. Then the commutators have a magnetic field addition, and to make it consistent, the magnetic field has to end up obeying Maxwell's equation, since the phase rotation gives a U(1) symmetry. This is not a true derivation of Maxwell's physics from quantum mechanics, it is just a way of showing you need the extra assumption of CP invariance to make the hopping hamiltonian real (which is true).
Then with the extra assumption, you just get
$$ i {d\psi\over dt} = {1\over 2} \nabla \cdot (t(x) \nabla \psi) + V(x) \psi $$
Where I have added back the potential. This is the continuum limit of a tight binding model with spatially slowly varying hopping, or inverse effective mass.
A: In addition to Claudius' and Ron Maimon's answers, I would like to make three comments:


*

*Classically, the Hamiltonian function for the effective mass approximation reads
$$\tag{1} H({\bf r}, {\bf p})~:=~\frac{{\bf p}^2}{2m^*({\bf r})}+V({\bf r}).$$

*Quantum mechanically, when one quantizes the classical model (1), one should pick a self-consistent choice for the Hamiltonian operator $\hat{H}$. It is natural to replace the classical variable ${\bf r}$ and ${\bf p}$ in the Hamiltonian (1) with the operators 
$$\tag{2}\hat{\bf r}~=~{\bf r} \qquad\text{and}\qquad \hat{\bf p}~=~\frac{\hbar}{i}\nabla $$
(in the Schrodinger representation). But which operator ordering prescription should one choose? One natural choice, which (under appropriate boundary conditions) makes the Hamiltonian Hermitian, is
$$\tag{3} \hat{H}~:=~\hat{\bf p}\cdot \frac{1}{2m^*(\hat{\bf r})}\hat{\bf p}+V(\hat{\bf r})~=~-\frac{\hbar^2}{2}\nabla\cdot \frac{1}{m^*({\bf r})}\nabla+V({\bf r}).$$

*Finally, let us mention a somewhat related/generalized Hermitian Hamiltonian operator
$$\tag{4} \hat{H}~=~-\frac{\hbar^2}{2}\Delta_g +V({\bf r}), $$
which may give another useful (anisotropic) effective mass model. Here $\Delta_g$ is the Laplace-Beltrami operator for a Riemannian $3\times 3$ metric $g_{ij}=g_{ij}({\bf r})$, which, roughly speaking, may be viewed as an (anisotropic) effective mass tensor.
A: I would be very surprised if you managed to find a strict mathematical derivation of the Schrödinger equation anywhere – at least I have not encountered one until now. However, it might be worth pointing out that the ‘general’ time-dependent Schrödinger equation, which is often taken as an axiom of quantum mechanics, is usually
$$i \hbar \partial_t \Psi = \hat H \Psi \quad .$$
In the case of a stationary Hamiltonian (usually $U(x,t) \equiv U(x)$), this equation separates and you get the stationary Schrödinger equation, namely
$$ \hat H \Psi = E \Psi \quad ,$$
that is, an eigenvalue equation for the Hamiltonian.
Given this equation, it is then relatively simple to work out the form of the Hamiltonian (in your case, $-\frac{\hbar^2}{2} \nabla \frac{1}{m^\star} \nabla + U$) and plug it into the equation. The exact form of the Hamiltonian is usually guesswork based on observation and analogies to classical mechanics. In general, we have
$$ \hat H = \hat T + \hat U $$
where $\hat T$ and $\hat U$ denote the operators for kinetic and potential energy, correspondingly.
It is worth noting that you can derive the continuity equation (which is identical to probability conservation in this case) from the Schrödinger equation by adding the complex conjugate of the Schrödinger equation to itself.
