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

Hi ballkikhaal - I edited out the part of your question asking about a book, because we limit the number of book recommendation questions on the site. See if anything in this question helps you.
–
David Z♦Oct 10 '12 at 16:23

5 Answers
5

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

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).

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.

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.

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.

thank you for the answer but what i am asking is how to derive it within the regime of effective mass approximation and that too when the mass has a spatial profile... for example in the case of Al/GaAs high electron mobility transistor we have a position dependent mass.
–
baalkikhaalOct 10 '12 at 16:55

Are you looking for a derivation of the Hamiltonian $\hat H$ or the Schrödinger equation? I am positive that neither probability conservation nor the continuity equation have anything to do with the earlier.
–
ClaudiusOct 10 '12 at 16:57

Where do you think Schrodinger equations come from, if not by some sort of a derivation?
–
Ron MaimonOct 10 '12 at 19:31

The answer is not by guesswork, it is from the tight-binding approximation with CP invariance to guarantee that the hopping parameter is real, and then Hermiticity guarantees the hopping is symmetric and equal to the given Hamiltonian. If the hopping is slowly locally varying, then you get the Hamiltonian they say. The Schrodinger equation which is axiomatic is not as specific as the Schrodinger equation in space, which has a specific ansatz for the kinetic term which can be justified from tight binding, as Feynman does in his lectures.
–
Ron MaimonOct 10 '12 at 19:47