Schrödinger equation

src: i.ytimg.com

In quantum mechanics, the Schrödinger equation is a mathematical equation that describes the changes over time of a physical system in which quantum effects, such as wave-particle duality, are significant. The equation is a mathematical formulation for studying quantum mechanical systems. It is considered a central result in the study of quantum systems and its derivation was a significant landmark in developing the theory of quantum mechanics. It was named after Erwin Schrödinger, who derived the equation in 1925 and published it in 1926, forming the basis for his work that resulted in Schrödinger being awarded the Nobel Prize in Physics in 1933. The equation is a type of differential equation known as a wave-equation, which serves as a mathematical model of the movement of waves.

In classical mechanics, Newton's second law (F = ma) is used to make a mathematical prediction as to what path a given system will take following a set of known initial conditions. In quantum mechanics, the analogue of Newton's law is Schrödinger's equation for a quantum system (usually atoms, molecules, and subatomic particles whether free, bound, or localised). It is not an algebraic equation, but in general a linear partial differential equation, describing the time-evolution of the system's wave function (also called a "state function").

The concept of a wavefunction is a fundamental postulate of quantum mechanics. Using these postulates, Schrödinger's equation can be derived from the fact that the time-evolution operator must be unitary and must therefore be generated by the exponential of a self-adjoint operator, which is the quantum Hamiltonian. This derivation is explained below.

In the Copenhagen interpretation of quantum mechanics, the wave function is the most complete description that can be given of a physical system. Solutions to Schrödinger's equation describe not only molecular, atomic, and subatomic systems, but also macroscopic systems, possibly even the whole universe. Schrödinger's equation is central to all applications of quantum mechanics including quantum field theory which combines special relativity with quantum mechanics. Theories of quantum gravity, such as string theory, also do not modify Schrödinger's equation.

The Schrödinger equation is not the only way to study quantum mechanical systems and make predictions, as there are other quantum mechanical formulations such as matrix mechanics, introduced by Werner Heisenberg, and path integral formulation, developed chiefly by Richard Feynman. Paul Dirac incorporated matrix mechanics and the Schrödinger equation into a single formulation.

Video Schrödinger equation

Equation

Time-dependent equation

The form of the Schrödinger equation depends on the physical situation (see below for special cases). The most general form is the time-dependent Schrödinger equation (TDSE), which gives a description of a system evolving with time:

where i is the imaginary unit, ? is the reduced Planck constant which is :${\displaystyle \hbar ={\frac {h}{2\pi }}}$, the symbol ?/?t indicates a partial derivative with respect to time t, ? (the Greek letter psi) is the wave function of the quantum system, r and t are the position vector and time respectively, and ? is the Hamiltonian operator (which characterises the total energy of the system under consideration).

The most famous example is the nonrelativistic Schrödinger equation for a single particle moving in an electric field (but not a magnetic field; see the Pauli equation):

where ? is the particle's "reduced mass", V is its potential energy, ?² is the Laplacian (a differential operator), and ? is the wave function (more precisely, in this context, it is called the "position-space wave function"). In plain language, it means "total energy equals kinetic energy plus potential energy", but the terms take unfamiliar forms for reasons explained below.

Given the particular differential operators involved, this is a linear partial differential equation. It is also a diffusion equation, but unlike the heat equation, this one is also a wave equation given the imaginary unit present in the transient term.

The term "Schrödinger equation" can refer to both the general equation, or the specific nonrelativistic version. The general equation is indeed quite general, used throughout quantum mechanics, for everything from the Dirac equation to quantum field theory, by plugging in diverse expressions for the Hamiltonian. The specific nonrelativistic version is a strictly classical approximation to reality and yields accurate results in many situations, but only to a certain extent (see relativistic quantum mechanics and relativistic quantum field theory).

To apply the Schrödinger equation, the Hamiltonian operator is set up for the system, accounting for the kinetic and potential energy of the particles constituting the system, then inserted into the Schrödinger equation. The resulting partial differential equation is solved for the wave function, which contains information about the system.

Time-independent equation

The time-dependent Schrödinger equation described above predicts that wave functions can form standing waves, called stationary states (also called "orbitals", as in atomic orbitals or molecular orbitals). These states are particularly important as their individual study later simplifies the task of solving the time-dependent Schrödinger equation for any state. Stationary states can also be described by a simpler form of the Schrödinger equation, the time-independent Schrödinger equation (TISE).

This is only used when the Hamiltonian itself is not dependent on time explicitly. However, even in this case the total wave function still has a time dependency.

In words, the equation states:

When the Hamiltonian operator acts on a certain wave function ?, and the result is proportional to the same wave function ?, then ? is a stationary state, and the proportionality constant, E, is the energy of the state ?.

In linear algebra terminology, this equation is an eigenvalue equation and in this sense the wave function is an eigenfunction of the Hamiltonian operator.

As before, the most common manifestation is the nonrelativistic Schrödinger equation for a single particle moving in an electric field (but not a magnetic field):

with definitions as above.

The time-independent Schrödinger equation is discussed further below.

Maps Schrödinger equation

Derivation

In the modern understanding of quantum mechanics, Schrödinger's equation may be derived as follows. If the wave-function at time t is given by ${\displaystyle \Psi (t)}$, then by the linearity of quantum mechanics the wave-function at time t' must be given by ${\displaystyle \Psi (t')=U(t',t)\Psi (t)}$, where ${\displaystyle U(t',t)}$ is a linear operator. Since time-evolution must preserve the norm of the wave-function, it follows that ${\displaystyle U(t',t)}$ must be a member of the unitary group of operators acting on wave-functions. We also know that when ${\displaystyle t'=t}$, we must have ${\displaystyle U(t,t)=1}$. Therefore, expanding the operator ${\displaystyle U(t',t)}$ for t' close to t, we can write ${\displaystyle U(t',t)=1-iH(t'-t)}$ where H is a Hermitian operator. This follows from the fact that the Lie algebra corresponding to the unitary group comprises Hermitian operators. Taking the limit as the time-difference ${\displaystyle t'-t}$ becomes very small, we obtain Schrödinger's equation.

So far, H is only an abstract Hermitian operator. However using the correspondence principle it is possible to show that, in the classical limit, the expectation value of H is indeed the classical energy. The correspondence principle does not completely fix the form of the quantum Hamiltonian due to the uncertainty principle and therefore the precise form of the quantum Hamiltonian must be fixed empirically.

src: i.ytimg.com

Implications

The Schrödinger equation and its solutions introduced a breakthrough in thinking about physics. Schrödinger's equation was the first of its type, and solutions led to consequences that were very unusual and unexpected for the time.

Total, kinetic, and potential energy

The overall form of the equation is not unusual or unexpected, as it uses the principle of the conservation of energy. The terms of the nonrelativistic Schrödinger equation can be interpreted as total energy of the system, equal to the system kinetic energy plus the system potential energy. In this respect, it is just the same as in classical physics.

Quantization

The Schrödinger equation predicts that if certain properties of a system are measured, the result may be quantized, meaning that only specific discrete values can occur. One example is energy quantization: the energy of an electron in an atom is always one of the quantized energy levels, a fact discovered via atomic spectroscopy. (Energy quantization is discussed below.) Another example is quantization of angular momentum. This was an assumption in the earlier Bohr model of the atom, but it is a prediction of the Schrödinger equation.

Another result of the Schrödinger equation is that not every measurement gives a quantized result in quantum mechanics. For example, position, momentum, time, and (in some situations) energy can have any value across a continuous range.

Measurement and uncertainty

In classical mechanics, a particle has, at every moment, an exact position and an exact momentum. These values change deterministically as the particle moves according to Newton's laws. Under the Copenhagen interpretation of quantum mechanics, particles do not have exactly determined properties, and when they are measured, the result is randomly drawn from a probability distribution. The Schrödinger equation predicts what the probability distributions are, but fundamentally cannot predict the exact result of each measurement.

The Heisenberg uncertainty principle is the statement of the inherent measurement uncertainty in quantum mechanics. It states that the more precisely a particle's position is known, the less precisely its momentum is known, and vice versa.

The Schrödinger equation describes the (deterministic) evolution of the wave function of a particle. However, even if the wave function is known exactly, the result of a specific measurement on the wave function is uncertain.

Quantum tunneling

In classical physics, when a ball is rolled slowly up a large hill, it will come to a stop and roll back, because it doesn't have enough energy to get over the top of the hill to the other side. However, the Schrödinger equation predicts that there is a small probability that the ball will get to the other side of the hill, even if it has too little energy to reach the top. This is called quantum tunneling. It is related to the distribution of energy: although the ball's assumed position seems to be on one side of the hill, there is a chance of finding it on the other side.

Particles as waves

The nonrelativistic Schrödinger equation is a type of partial differential equation called a wave equation. Therefore, it is often said particles can exhibit behavior usually attributed to waves. In some modern interpretations this description is reversed - the quantum state, i.e. wave, is the only genuine physical reality, and under the appropriate conditions it can show features of particle-like behavior. However, Ballentine shows that such an interpretation has problems. Ballentine points out that whilst it is arguable to associate a physical wave with a single particle, there is still only one Schrödinger wave equation for many particles. He points out:

"If a physical wave field were associated with a particle, or if a particle were identified with a wave packet, then corresponding to N interacting particles there should be N interacting waves in ordinary three-dimensional space. But according to (4.6) that is not the case; instead there is one "wave" function in an abstract 3N-dimensional configuration space. The misinterpretation of psi as a physical wave in ordinary space is possible only because the most common applications of quantum mechanics are to one-particle states, for which configuration space and ordinary space are isomorphic."

Two-slit diffraction is a famous example of the strange behaviors that waves regularly display, that are not intuitively associated with particles. The overlapping waves from the two slits cancel each other out in some locations, and reinforce each other in other locations, causing a complex pattern to emerge. Intuitively, one would not expect this pattern from firing a single particle at the slits, because the particle should pass through one slit or the other, not a complex overlap of both.

However, since the Schrödinger equation is a wave equation, a single particle fired through a double-slit does show this same pattern (figure on right). Note: The experiment must be repeated many times for the complex pattern to emerge. Although this is counterintuitive, the prediction is correct; in particular, electron diffraction and neutron diffraction are well understood and widely used in science and engineering.

Related to diffraction, particles also display superposition and interference.

The superposition property allows the particle to be in a quantum superposition of two or more quantum states at the same time. However, it is noted that a "quantum state" in quantum mechanics means the probability that a system will be, for example at a position x, not that the system will actually be at position x. It does not imply that the particle itself may be in two classical states at once. Indeed, quantum mechanics is generally unable to assign values for properties prior to measurement at all.

Many-Worlds interpretation

In Dublin in 1952 Erwin Schrödinger gave a lecture in which at one point he jocularly warned his audience that what he was about to say might "seem lunatic". It was that, when his Nobel equations seem to be describing several different histories, they are "not alternatives but all really happen simultaneously". This is the earliest known reference to the Many-worlds interpretation of quantum mechanics.

src: skyranch.org

Interpretation of the wave function

The Schrödinger equation provides a way to calculate the wave function of a system and how it changes dynamically in time. However, the Schrödinger equation does not directly say what, exactly, the wave function is. Interpretations of quantum mechanics address questions such as what the relation is between the wave function, the underlying reality, and the results of experimental measurements.

An important aspect is the relationship between the Schrödinger equation and wavefunction collapse. In the oldest Copenhagen interpretation, particles follow the Schrödinger equation except during wavefunction collapse, during which they behave entirely differently. The advent of quantum decoherence theory allowed alternative approaches (such as the Everett many-worlds interpretation and consistent histories), wherein the Schrödinger equation is always satisfied, and wavefunction collapse should be explained as a consequence of the Schrödinger equation.

src: www.mathworks.com

Historical background and development

Following Max Planck's quantization of light (see black body radiation), Albert Einstein interpreted Planck's quanta to be photons, particles of light, and proposed that the energy of a photon is proportional to its frequency, one of the first signs of wave-particle duality. Since energy and momentum are related in the same way as frequency and wavenumber in special relativity, it followed that the momentum p of a photon is inversely proportional to its wavelength ?, or proportional to its wavenumber k:

${\displaystyle p={\frac {h}{\lambda }}=\hbar k}$

where h is Planck's constant. Louis de Broglie hypothesized that this is true for all particles, even particles which have mass such as electrons. He showed that, assuming that the matter waves propagate along with their particle counterparts, electrons form standing waves, meaning that only certain discrete rotational frequencies about the nucleus of an atom are allowed. These quantized orbits correspond to discrete energy levels, and de Broglie reproduced the Bohr model formula for the energy levels. The Bohr model was based on the assumed quantization of angular momentum L according to:

${\displaystyle L=n{h \over 2\pi }=n\hbar .}$

According to de Broglie the electron is described by a wave and a whole number of wavelengths must fit along the circumference of the electron's orbit:

${\displaystyle n\lambda =2\pi r.\,}$

This approach essentially confined the electron wave in one dimension, along a circular orbit of radius r.

In 1921, prior to de Broglie, Arthur C. Lunn at the University of Chicago had used the same argument based on the completion of the relativistic energy-momentum 4-vector to derive what we now call the de Broglie relation. Unlike de Broglie, Lunn went on to formulate the differential equation now known as the Schrödinger equation, and solve for its energy eigenvalues for the hydrogen atom. Unfortunately the paper was rejected by the Physical Review, as recounted by Kamen.

Following up on de Broglie's ideas, physicist Peter Debye made an offhand comment that if particles behaved as waves, they should satisfy some sort of wave equation. Inspired by Debye's remark, Schrödinger decided to find a proper 3-dimensional wave equation for the electron. He was guided by William R. Hamilton's analogy between mechanics and optics, encoded in the observation that the zero-wavelength limit of optics resembles a mechanical system--the trajectories of light rays become sharp tracks that obey Fermat's principle, an analog of the principle of least action. A modern version of his reasoning is reproduced below. The equation he found is:

${\displaystyle i\hbar {\frac {\partial }{\partial t}}\Psi (\mathbf {r} ,\,t)=-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\Psi (\mathbf {r} ,\,t)+V(\mathbf {r} )\Psi (\mathbf {r} ,\,t).}$

However, by that time, Arnold Sommerfeld had refined the Bohr model with relativistic corrections. Schrödinger used the relativistic energy momentum relation to find what is now known as the Klein-Gordon equation in a Coulomb potential (in natural units):

${\displaystyle \left(E+{e^{2} \over r}\right)^{2}\psi (x)=-\nabla ^{2}\psi (x)+m^{2}\psi (x).}$

He found the standing waves of this relativistic equation, but the relativistic corrections disagreed with Sommerfeld's formula. Discouraged, he put away his calculations and secluded himself in an isolated mountain cabin in December 1925.

While at the cabin, Schrödinger decided that his earlier nonrelativistic calculations were novel enough to publish, and decided to leave off the problem of relativistic corrections for the future. Despite the difficulties in solving the differential equation for hydrogen (he had sought help from his friend the mathematician Hermann Weyl) Schrödinger showed that his nonrelativistic version of the wave equation produced the correct spectral energies of hydrogen in a paper published in 1926. In the equation, Schrödinger computed the hydrogen spectral series by treating a hydrogen atom's electron as a wave ?(x, t), moving in a potential well V, created by the proton. This computation accurately reproduced the energy levels of the Bohr model. In a paper, Schrödinger himself explained this equation as follows:

This 1926 paper was enthusiastically endorsed by Einstein, who saw the matter-waves as an intuitive depiction of nature, as opposed to Heisenberg's matrix mechanics, which he considered overly formal.

The Schrödinger equation details the behavior of ? but says nothing of its nature. Schrödinger tried to interpret it as a charge density in his fourth paper, but he was unsuccessful. In 1926, just a few days after Schrödinger's fourth and final paper was published, Max Born successfully interpreted ? as the probability amplitude, whose absolute square is equal to probability density. Schrödinger, though, always opposed a statistical or probabilistic approach, with its associated discontinuities--much like Einstein, who believed that quantum mechanics was a statistical approximation to an underlying deterministic theory--and never reconciled with the Copenhagen interpretation.

Louis de Broglie in his later years proposed a real valued wave function connected to the complex wave function by a proportionality constant and developed the De Broglie-Bohm theory.

src: i.ytimg.com

The wave equation for particles

The Schrödinger equation is a diffusion equation, the solutions are functions which describe wave-like motions. Wave equations in physics can normally be derived from other physical laws - the wave equation for mechanical vibrations on strings and in matter can be derived from Newton's laws, where the wave function represents the displacement of matter, and electromagnetic waves from Maxwell's equations, where the wave functions are electric and magnetic fields. The basis for Schrödinger's equation, on the other hand, is the energy of the system and a separate postulate of quantum mechanics: the wave function is a description of the system. The Schrödinger equation is therefore a new concept in itself; as Feynman put it:

The foundation of the equation is structured to be a linear differential equation based on classical energy conservation, and consistent with the De Broglie relations. The solution is the wave function ?, which contains all the information that can be known about the system. In the Copenhagen interpretation, the modulus of ? is related to the probability the particles are in some spatial configuration at some instant of time. Solving the equation for ? can be used to predict how the particles will behave under the influence of the specified potential and with each other.

The Schrödinger equation was developed principally from the De Broglie hypothesis, a wave equation that would describe particles, and can be constructed as shown informally in the following sections. For a more rigorous description of Schrödinger's equation, see also Resnick et al.

Consistency with energy conservation

The total energy E of a particle is the sum of kinetic energy T and potential energy V, this sum is also the frequent expression for the Hamiltonian H in classical mechanics:

${\displaystyle E=T+V=H\,\!}$

Explicitly, for a particle in one dimension with position x, mass m and momentum p, and potential energy V which generally varies with position and time t:

${\displaystyle E={\frac {p^{2}}{2m}}+V(x,t)=H.}$

For three dimensions, the position vector r and momentum vector p must be used:

${\displaystyle E={\frac {\mathbf {p} \cdot \mathbf {p} }{2m}}+V(\mathbf {r} ,t)=H}$

This formalism can be extended to any fixed number of particles: the total energy of the system is then the total kinetic energies of the particles, plus the total potential energy, again the Hamiltonian. However, there can be interactions between the particles (an N-body problem), so the potential energy V can change as the spatial configuration of particles changes, and possibly with time. The potential energy, in general, is not the sum of the separate potential energies for each particle, it is a function of all the spatial positions of the particles. Explicitly:

${\displaystyle E=\sum _{n=1}^{N}{\frac {\mathbf {p} _{n}\cdot \mathbf {p} _{n}}{2m_{n}}}+V(\mathbf {r} _{1},\mathbf {r} _{2}\cdots \mathbf {r} _{N},t)=H\,\!}$

Linearity

The simplest wavefunction is a plane wave of the form:

${\displaystyle \Psi (\mathbf {r} ,t)=Ae^{i(\mathbf {k} \cdot \mathbf {r} -\omega t)}\,\!}$

where the A is the amplitude, k the wavevector, and ? the angular frequency, of the plane wave. In general, physical situations are not purely described by plane waves, so for generality the superposition principle is required; any wave can be made by superposition of sinusoidal plane waves. So if the equation is linear, a linear combination of plane waves is also an allowed solution. Hence a necessary and separate requirement is that the Schrödinger equation is a linear differential equation.

For discrete k the sum is a superposition of plane waves:

${\displaystyle \Psi (\mathbf {r} ,t)=\sum _{n=1}^{\infty }A_{n}e^{i(\mathbf {k} _{n}\cdot \mathbf {r} -\omega _{n}t)}\,\!}$

for some real amplitude coefficients A_n, and for continuous k the sum becomes an integral, the Fourier transform of a momentum space wavefunction:

${\displaystyle \Psi (\mathbf {r} ,t)={\frac {1}{({\sqrt {2\pi }})^{3}}}\int \Phi (\mathbf {k} )e^{i(\mathbf {k} \cdot \mathbf {r} -\omega t)}d^{3}\mathbf {k} \,\!}$

where d³k = dk_xdk_ydk_z is the differential volume element in k-space, and the integrals are taken over all k-space. The momentum wavefunction ?(k) arises in the integrand since the position and momentum space wavefunctions are Fourier transforms of each other.

Consistency with the De Broglie relations

Einstein's light quanta hypothesis (1905) states that the energy E of a photon is proportional to the frequency ? (or angular frequency, ? = 2??) of the corresponding quantum wavepacket of light:

${\displaystyle E=h\nu =\hbar \omega \,\!}$

Likewise De Broglie's hypothesis (1924) states that any particle can be associated with a wave, and that the momentum p of the particle is inversely proportional to the wavelength ? of such a wave (or proportional to the wavenumber, k = 2?/?), in one dimension, by:

${\displaystyle p={\frac {h}{\lambda }}=\hbar k\;,}$

while in three dimensions, wavelength ? is related to the magnitude of the wavevector k:

${\displaystyle \mathbf {p} =\hbar \mathbf {k} \,,\quad |\mathbf {k} |={\frac {2\pi }{\lambda }}\,.}$

The Planck-Einstein and de Broglie relations illuminate the deep connections between energy with time, and space with momentum, and express wave-particle duality. In practice, natural units comprising ? = 1 are used, as the De Broglie equations reduce to identities: allowing momentum, wavenumber, energy and frequency to be used interchangeably, to prevent duplication of quantities, and reduce the number of dimensions of related quantities. For familiarity SI units are still used in this article.

Schrödinger's insight, late in 1925, was to express the phase of a plane wave as a complex phase factor using these relations:

${\displaystyle \Psi =Ae^{i(\mathbf {k} \cdot \mathbf {r} -\omega t)}=Ae^{i(\mathbf {p} \cdot \mathbf {r} -Et)/\hbar }}$

and to realize that the first order partial derivatives were:

with respect to space:

${\displaystyle \nabla \Psi ={\dfrac {i}{\hbar }}\mathbf {p} Ae^{i(\mathbf {p} \cdot \mathbf {r} -Et)/\hbar }={\dfrac {i}{\hbar }}\mathbf {p} \Psi }$

with respect to time:

${\displaystyle {\dfrac {\partial \Psi }{\partial t}}=-{\dfrac {iE}{\hbar }}Ae^{i(\mathbf {p} \cdot \mathbf {r} -Et)/\hbar }=-{\dfrac {iE}{\hbar }}\Psi }$

Another postulate of quantum mechanics is that all observables are represented by linear Hermitian operators which act on the wavefunction, and the eigenvalues of the operator are the values the observable takes. The previous derivatives are consistent with the energy operator, corresponding to the time derivative,

${\displaystyle {\hat {E}}\Psi =i\hbar {\dfrac {\partial }{\partial t}}\Psi =E\Psi }$

where E are the energy eigenvalues, and the momentum operator, corresponding to the spatial derivatives (the gradient ?),

${\displaystyle {\hat {\mathbf {p} }}\Psi =-i\hbar \nabla \Psi =\mathbf {p} \Psi }$

where p is a vector of the momentum eigenvalues. In the above, the "hats" ( ^ ) indicate these observables are operators, not simply ordinary numbers or vectors. The energy and momentum operators are differential operators, while the potential energy function V is just a multiplicative factor.

Substituting the energy and momentum operators into the classical energy conservation equation obtains the operator:

${\displaystyle E={\dfrac {\mathbf {p} \cdot \mathbf {p} }{2m}}+V\quad \rightarrow \quad {\hat {E}}={\dfrac {{\hat {\mathbf {p} }}\cdot {\hat {\mathbf {p} }}}{2m}}+V}$

so in terms of derivatives with respect to time and space, acting this operator on the wavefunction ? immediately led Schrödinger to his equation:

${\displaystyle i\hbar {\dfrac {\partial \Psi }{\partial t}}=-{\dfrac {\hbar ^{2}}{2m}}\nabla ^{2}\Psi +V\Psi }$

Wave-particle duality can be assessed from these equations as follows. The kinetic energy T is related to the square of momentum p. As the particle's momentum increases, the kinetic energy increases more rapidly, but since the wavenumber |k| increases the wavelength ? decreases. In terms of ordinary scalar and vector quantities (not operators):

${\displaystyle \mathbf {p} \cdot \mathbf {p} \propto \mathbf {k} \cdot \mathbf {k} \propto T\propto {\dfrac {1}{\lambda ^{2}}}}$

The kinetic energy is also proportional to the second spatial derivatives, so it is also proportional to the magnitude of the curvature of the wave, in terms of operators:

${\displaystyle {\hat {T}}\Psi ={\frac {-\hbar ^{2}}{2m}}\nabla \cdot \nabla \Psi \,\propto \,\nabla ^{2}\Psi \,.}$

As the curvature increases, the amplitude of the wave alternates between positive and negative more rapidly, and also shortens the wavelength. So the inverse relation between momentum and wavelength is consistent with the energy the particle has, and so the energy of the particle has a connection to a wave, all in the same mathematical formulation.

Wave and particle motion

Schrödinger required that a wave packet solution near position r with wavevector near k will move along the trajectory determined by classical mechanics for times short enough for the spread in k (and hence in velocity) not to substantially increase the spread in r. Since, for a given spread in k, the spread in velocity is proportional to Planck's constant ?, it is sometimes said that in the limit as ? approaches zero, the equations of classical mechanics are restored from quantum mechanics. Great care is required in how that limit is taken, and in what cases.

The limiting short-wavelength is equivalent to ? tending to zero because this is limiting case of increasing the wave packet localization to the definite position of the particle (see images right). Using the Heisenberg uncertainty principle for position and momentum, the products of uncertainty in position and momentum become zero as ? -> 0:

${\displaystyle \sigma (x)\sigma (p_{x})\geqslant {\frac {\hbar }{2}}\quad \rightarrow \quad \sigma (x)\sigma (p_{x})\geqslant 0\,\!}$

where ? denotes the (root mean square) measurement uncertainty in x and p_x (and similarly for the y and z directions) which implies the position and momentum can only be known to arbitrary precision in this limit.

The Schrödinger equation in its general form

${\displaystyle i\hbar {\frac {\partial }{\partial t}}\Psi \left(\mathbf {r} ,t\right)={\hat {H}}\Psi \left(\mathbf {r} ,t\right)\,\!}$

is closely related to the Hamilton-Jacobi equation (HJE)

${\displaystyle -{\frac {\partial }{\partial t}}S(q_{i},t)=H\left(q_{i},{\frac {\partial S}{\partial q_{i}}},t\right)\,\!}$

where S is action and H is the Hamiltonian function (not operator). Here the generalized coordinates q_i for i = 1, 2, 3 (used in the context of the HJE) can be set to the position in Cartesian coordinates as r = (q₁, q₂, q₃) = (x, y, z).

Substituting

${\displaystyle \Psi ={\sqrt {\rho (\mathbf {r} ,t)}}e^{iS(\mathbf {r} ,t)/\hbar }\,\!}$

where ? is the probability density, into the Schrödinger equation and then taking the limit ? -> 0 in the resulting equation, yields the Hamilton-Jacobi equation.

The implications are as follows:

• The motion of a particle, described by a (short-wavelength) wave packet solution to the Schrödinger equation, is also described by the Hamilton-Jacobi equation of motion.
• The Schrödinger equation includes the wavefunction, so its wave packet solution implies the position of a (quantum) particle is fuzzily spread out in wave fronts. On the contrary, the Hamilton-Jacobi equation applies to a (classical) particle of definite position and momentum, instead the position and momentum at all times (the trajectory) are deterministic and can be simultaneously known.

src: c8.alamy.com

Nonrelativistic quantum mechanics

The quantum mechanics of particles without accounting for the effects of special relativity, for example particles propagating at speeds much less than light, is known as nonrelativistic quantum mechanics. Following are several forms of Schrödinger's equation in this context for different situations: time independence and dependence, one and three spatial dimensions, and one and N particles.

In actuality, the particles constituting the system do not have the numerical labels used in theory. The language of mathematics forces us to label the positions of particles one way or another, otherwise there would be confusion between symbols representing which variables are for which particle.

Time independent

If the Hamiltonian is not an explicit function of time, the equation is separable into a product of spatial and temporal parts. In general, the wavefunction takes the form:

${\displaystyle \Psi ({\text{space coords}},t)=\psi ({\text{space coords}})\tau (t)\,.}$

where ?(space coords) is a function of all the spatial coordinate(s) of the particle(s) constituting the system only, and ?(t) is a function of time only.

Substituting for ? into the Schrödinger equation for the relevant number of particles in the relevant number of dimensions, solving by separation of variables implies the general solution of the time-dependent equation has the form:

${\displaystyle \Psi ({\text{space coords}},t)=\psi ({\text{space coords}})e^{-i{Et/\hbar }}\,.}$

Since the time dependent phase factor is always the same, only the spatial part needs to be solved for in time independent problems. Additionally, the energy operator Ê = i??/?t can always be replaced by the energy eigenvalue E, thus the time independent Schrödinger equation is an eigenvalue equation for the Hamiltonian operator:

${\displaystyle {\hat {H}}\psi =E\psi }$

This is true for any number of particles in any number of dimensions (in a time independent potential). This case describes the standing wave solutions of the time-dependent equation, which are the states with definite energy (instead of a probability distribution of different energies). In physics, these standing waves are called "stationary states" or "energy eigenstates"; in chemistry they are called "atomic orbitals" or "molecular orbitals". Superpositions of energy eigenstates change their properties according to the relative phases between the energy levels.

The energy eigenvalues from this equation form a discrete spectrum of values, so mathematically energy must be quantized. More specifically, the energy eigenstates form a basis - any wavefunction may be written as a sum over the discrete energy states or an integral over continuous energy states, or more generally as an integral over a measure. This is the spectral theorem in mathematics, and in a finite state space it is just a statement of the completeness of the eigenvectors of a Hermitian matrix.

One-dimensional examples

For a particle in one dimension, the Hamiltonian is:

${\displaystyle {\hat {H}}={\frac {{\hat {p}}^{2}}{2m}}+V(x)\,,\quad {\hat {p}}=-i\hbar {\frac {d}{dx}}}$

and substituting this into the general Schrödinger equation gives:

${\displaystyle -{\frac {\hbar ^{2}}{2m}}{\frac {d^{2}}{dx^{2}}}\psi (x)+V(x)\psi (x)=E\psi (x)}$

This is the only case the Schrödinger equation is an ordinary differential equation, rather than a partial differential equation. The general solutions are always of the form:

${\displaystyle \Psi (x,t)=\psi (x)e^{-iEt/\hbar }\,.}$

For N particles in one dimension, the Hamiltonian is:

${\displaystyle {\hat {H}}=\sum _{n=1}^{N}{\frac {{\hat {p}}_{n}^{2}}{2m_{n}}}+V(x_{1},x_{2},\cdots x_{N})\,,\quad {\hat {p}}_{n}=-i\hbar {\frac {\partial }{\partial x_{n}}}}$

where the position of particle n is x_n. The corresponding Schrödinger equation is:

${\displaystyle -{\frac {\hbar ^{2}}{2}}\sum _{n=1}^{N}{\frac {1}{m_{n}}}{\frac {\partial ^{2}}{\partial x_{n}^{2}}}\psi (x_{1},x_{2},\cdots x_{N})+V(x_{1},x_{2},\cdots x_{N})\psi (x_{1},x_{2},\cdots x_{N})=E\psi (x_{1},x_{2},\cdots x_{N})\,.}$

so the general solutions have the form:

${\displaystyle \Psi (x_{1},x_{2},\cdots x_{N},t)=e^{-iEt/\hbar }\psi (x_{1},x_{2}\cdots x_{N})}$

For non-interacting distinguishable particles, the potential of the system only influences each particle separately, so the total potential energy is the sum of potential energies for each particle:

${\displaystyle V(x_{1},x_{2},\cdots x_{N})=\sum _{n=1}^{N}V(x_{n})\,.}$

and the wavefunction can be written as a product of the wavefunctions for each particle:

${\displaystyle \Psi (x_{1},x_{2},\cdots x_{N},t)=e^{-i{Et/\hbar }}\prod _{n=1}^{N}\psi (x_{n})\,,}$

For non-interacting identical particles, the potential is still a sum, but wavefunction is a bit more complicated - it is a sum over the permutations of products of the separate wavefunctions to account for particle exchange. In general for interacting particles, the above decompositions are not possible.

Free particle

For no potential, V = 0, so the particle is free and the equation reads:

${\displaystyle -E\psi ={\frac {\hbar ^{2}}{2m}}{d^{2}\psi \over dx^{2}}\,}$

which has oscillatory solutions for E > 0 (the C_n are arbitrary constants):

${\displaystyle \psi _{E}(x)=C_{1}e^{i{\sqrt {2mE/\hbar ^{2}}}\,x}+C_{2}e^{-i{\sqrt {2mE/\hbar ^{2}}}\,x}\,}$

and exponential solutions for E < 0

${\displaystyle \psi _{-|E|}(x)=C_{1}e^{{\sqrt {2m|E|/\hbar ^{2}}}\,x}+C_{2}e^{-{\sqrt {2m|E|/\hbar ^{2}}}\,x}.\,}$

The exponentially growing solutions have an infinite norm, and are not physical. They are not allowed in a finite volume with periodic or fixed boundary conditions.

See also free particle and wavepacket for more discussion on the free particle.

Constant potential

For a constant potential, V = V₀, the solution is oscillatory for E > V₀ and exponential for E < V₀, corresponding to energies that are allowed or disallowed in classical mechanics. Oscillatory solutions have a classically allowed energy and correspond to actual classical motions, while the exponential solutions have a disallowed energy and describe a small amount of quantum bleeding into the classically disallowed region, due to quantum tunneling. If the potential V₀ grows to infinity, the motion is classically confined to a finite region. Viewed far enough away, every solution is reduced to an exponential; the condition that the exponential is decreasing restricts the energy levels to a discrete set, called the allowed energies.

Harmonic oscillator

The Schrödinger equation for this situation is

${\displaystyle E\psi =-{\frac {\hbar ^{2}}{2m}}{\frac {d^{2}}{dx^{2}}}\psi +{\frac {1}{2}}m\omega ^{2}x^{2}\psi }$

It is a notable quantum system to solve for; since the solutions are exact (but complicated - in terms of Hermite polynomials), and it can describe or at least approximate a wide variety of other systems, including vibrating atoms, molecules, and atoms or ions in lattices, and approximating other potentials near equilibrium points. It is also the basis of perturbation methods in quantum mechanics.

There is a family of solutions - in the position basis they are

${\displaystyle \psi _{n}(x)={\sqrt {\frac {1}{2^{n}\,n!}}}\cdot \left({\frac {m\omega }{\pi \hbar }}\right)^{1/4}\cdot e^{-{\frac {m\omega x^{2}}{2\hbar }}}\cdot H_{n}\left({\sqrt {\frac {m\omega }{\hbar }}}x\right)}$

where n = 0,1,2,..., and the functions H_n are the Hermite polynomials.

Three-dimensional examples

The extension from one dimension to three dimensions is straightforward, all position and momentum operators are replaced by their three-dimensional expressions and the partial derivative with respect to space is replaced by the gradient operator.

The Hamiltonian for one particle in three dimensions is:

${\displaystyle {\hat {H}}={\frac {{\hat {\mathbf {p} }}\cdot {\hat {\mathbf {p} }}}{2m}}+V(\mathbf {r} )\,,\quad {\hat {\mathbf {p} }}=-i\hbar \nabla }$

generating the equation:

${\displaystyle -{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi (\mathbf {r} )+V(\mathbf {r} )\psi (\mathbf {r} )=E\psi (\mathbf {r} )}$

with stationary state solutions of the form:

${\displaystyle \Psi (\mathbf {r} ,t)=\psi (\mathbf {r} )e^{-iEt/\hbar }}$

where the position of the particle is r. Two useful coordinate systems for solving the Schrödinger equation are Cartesian coordinates so that r = (x, y, z) and spherical polar coordinates so that r = (r, ?, ?), although other orthogonal coordinates are useful for solving the equation for systems with certain geometric symmetries.

For N particles in three dimensions, the Hamiltonian is:

${\displaystyle {\hat {H}}=\sum _{n=1}^{N}{\frac {{\hat {\mathbf {p} }}_{n}\cdot {\hat {\mathbf {p} }}_{n}}{2m_{n}}}+V(\mathbf {r} _{1},\mathbf {r} _{2},\cdots \mathbf {r} _{N})\,,\quad {\hat {\mathbf {p} }}_{n}=-i\hbar \nabla _{n}}$

where the position of particle n is r_n and the gradient operators are partial derivatives with respect to the particle's position coordinates. In Cartesian coordinates, for particle n, the position vector is r_n = (x_n, y_n, z_n) while the gradient and Laplacian operator are respectively:

${\displaystyle \nabla _{n}=\mathbf {e} _{x}{\frac {\partial }{\partial x_{n}}}+\mathbf {e} _{y}{\frac {\partial }{\partial y_{n}}}+\mathbf {e} _{z}{\frac {\partial }{\partial z_{n}}}\,,\quad \nabla _{n}^{2}=\nabla _{n}\cdot \nabla _{n}={\frac {\partial ^{2}}{{\partial x_{n}}^{2}}}+{\frac {\partial ^{2}}{{\partial y_{n}}^{2}}}+{\frac {\partial ^{2}}{{\partial z_{n}}^{2}}}}$

The Schrödinger equation is:

${\displaystyle -{\frac {\hbar ^{2}}{2}}\sum _{n=1}^{N}{\frac {1}{m_{n}}}\nabla _{n}^{2}\Psi (\mathbf {r} _{1},\mathbf {r} _{2},\cdots \mathbf {r} _{N})+V(\mathbf {r} _{1},\mathbf {r} _{2},\cdots \mathbf {r} _{N})\Psi (\mathbf {r} _{1},\mathbf {r} _{2},\cdots \mathbf {r} _{N})=E\Psi (\mathbf {r} _{1},\mathbf {r} _{2},\cdots \mathbf {r} _{N})}$

with stationary state solutions:

${\displaystyle \Psi (\mathbf {r} _{1},\mathbf {r} _{2}\cdots \mathbf {r} _{N},t)=e^{-iEt/\hbar }\psi (\mathbf {r} _{1},\mathbf {r} _{2}\cdots \mathbf {r} _{N})}$

Again, for non-interacting distinguishable particles the potential is the sum of particle potentials

${\displaystyle V(\mathbf {r} _{1},\mathbf {r} _{2},\cdots \mathbf {r} _{N})=\sum _{n=1}^{N}V(\mathbf {r} _{n})}$

and the wavefunction is a product of the particle wavefunctions

${\displaystyle \Psi (\mathbf {r} _{1},\mathbf {r} _{2}\cdots \mathbf {r} _{N},t)=e^{-i{Et/\hbar }}\prod _{n=1}^{N}\psi (\mathbf {r} _{n})\,.}$

For non-interacting identical particles, the potential is a sum but the wavefunction is a sum over permutations of products. The previous two equations do not apply to interacting particles.

Following are examples where exact solutions are known. See the main articles for further details.

Hydrogen atom

This form of the Schrödinger equation can be applied to the hydrogen atom:

${\displaystyle E\psi =-{\frac {\hbar ^{2}}{2\mu }}\nabla ^{2}\psi -{\frac {e^{2}}{4\pi \varepsilon _{0}r}}\psi }$

where e is the electron charge, r is the position of the electron (r = |r| is the magnitude of the position), the potential term is due to the Coulomb interaction, wherein ?₀ is the electric constant (permittivity of free space) and

${\displaystyle \mu ={\frac {m_{e}m_{p}}{m_{e}+m_{p}}}}$

is the 2-body reduced mass of the hydrogen nucleus (just a proton) of mass m_p and the electron of mass m_e. The negative sign arises in the potential term since the proton and electron are oppositely charged. The reduced mass in place of the electron mass is used since the electron and proton together orbit each other about a common centre of mass, and constitute a two-body problem to solve. The motion of the electron is of principle interest here, so the equivalent one-body problem is the motion of the electron using the reduced mass.

The wavefunction for hydrogen is a function of the electron's coordinates, and in fact can be separated into functions of each coordinate. Usually this is done in spherical polar coordinates:

${\displaystyle \psi (r,\theta ,\phi )=R(r)Y_{\ell }^{m}(\theta ,\phi )=R(r)\Theta (\theta )\Phi (\phi )}$

where R are radial functions and Ym
l(?, ?) are spherical harmonics of degree l and order m. This is the only atom for which the Schrödinger equation has been solved for exactly. Multi-electron atoms require approximative methods. The family of solutions are:

${\displaystyle \psi _{n\ell m}(r,\theta ,\phi )={\sqrt {{\left({\frac {2}{na_{0}}}\right)}^{3}{\frac {(n-\ell -1)!}{2n[(n+\ell )!]}}}}e^{-r/na_{0}}\left({\frac {2r}{na_{0}}}\right)^{\ell }L_{n-\ell -1}^{2\ell +1}\left({\frac {2r}{na_{0}}}\right)\cdot Y_{\ell }^{m}(\theta ,\phi )}$

where:

• ${\displaystyle a_{0}={\frac {4\pi \varepsilon _{0}\hbar ^{2}}{m_{e}e^{2}}}}$ is the Bohr radius,
• ${\displaystyle L_{n-\ell -1}^{2\ell +1}(\cdots )}$ are the generalized Laguerre polynomials of degree n - l - 1.
• n, l, m are the principal, azimuthal, and magnetic quantum numbers respectively: which take the values:

${\displaystyle {\begin{aligned}n&=1,2,3,\dots \\\ell &=0,1,2,\dots ,n-1\\m&=-\ell ,\dots ,\ell \\\end{aligned}}}$

NB: generalized Laguerre polynomials are defined differently by different authors--see main article on them and the hydrogen atom.

Two-electron atoms or ions

The equation for any two-electron system, such as the neutral helium atom (He, Z = 2), the negative hydrogen ion (H^-, Z = 1), or the positive lithium ion (Li⁺, Z = 3) is:

${\displaystyle E\psi =-\hbar ^{2}\left[{\frac {1}{2\mu }}\left(\nabla _{1}^{2}+\nabla _{2}^{2}\right)+{\frac {1}{M}}\nabla _{1}\cdot \nabla _{2}\right]\psi +{\frac {e^{2}}{4\pi \varepsilon _{0}}}\left[{\frac {1}{r_{12}}}-Z\left({\frac {1}{r_{1}}}+{\frac {1}{r_{2}}}\right)\right]\psi }$

where r₁ is the position of one electron (r₁ = |r₁| is its magnitude), r₂ is the position of the other electron (r₂ = |r₂| is the magnitude), r₁₂ = |r₁₂| is the magnitude of the separation between them given by

${\displaystyle |\mathbf {r} _{12}|=|\mathbf {r} _{2}-\mathbf {r} _{1}|\,\!}$

? is again the two-body reduced mass of an electron with respect to the nucleus of mass M, so this time

${\displaystyle \mu ={\frac {m_{e}M}{m_{e}+M}}\,\!}$

and Z is the atomic number for the element (not a quantum number).

The cross-term of two laplacians

${\displaystyle {\frac {1}{M}}\nabla _{1}\cdot \nabla _{2}\,\!}$

is known as the mass polarization term, which arises due to the motion of atomic nuclei. The wavefunction is a function of the two electron's positions:

${\displaystyle \psi =\psi (\mathbf {r} _{1},\mathbf {r} _{2}).}$

There is no closed form solution for this equation.

Time dependent

This is the equation of motion for the quantum state. In the most general form, it is written:

${\displaystyle i\hbar {\frac {\partial }{\partial t}}\Psi ={\hat {H}}\Psi .}$

and the solution, the wavefunction, is a function of all the particle coordinates of the system and time. Following are specific cases.

For one particle in one dimension, the Hamiltonian

${\displaystyle {\hat {H}}={\frac {{\hat {p}}^{2}}{2m}}+V(x,t)\,,\quad {\hat {p}}=-i\hbar {\frac {\partial }{\partial x}}}$

generates the equation:

${\displaystyle i\hbar {\frac {\partial }{\partial t}}\Psi (x,t)=-{\frac {\hbar ^{2}}{2m}}{\frac {\partial ^{2}}{\partial x^{2}}}\Psi (x,t)+V(x,t)\Psi (x,t)}$

For N particles in one dimension, the Hamiltonian is:

${\displaystyle {\hat {H}}=\sum _{n=1}^{N}{\frac {{\hat {p}}_{n}^{2}}{2m_{n}}}+V(x_{1},x_{2},\cdots x_{N},t)\,,\quad {\hat {p}}_{n}=-i\hbar {\frac {\partial }{\partial x_{n}}}}$

where the position of particle n is x_n, generating the equation:

${\displaystyle i\hbar {\frac {\partial }{\partial t}}\Psi (x_{1},x_{2}\cdots x_{N},t)=-{\frac {\hbar ^{2}}{2}}\sum _{n=1}^{N}{\frac {1}{m_{n}}}{\frac {\partial ^{2}}{\partial x_{n}^{2}}}\Psi (x_{1},x_{2}\cdots x_{N},t)+V(x_{1},x_{2}\cdots x_{N},t)\Psi (x_{1},x_{2}\cdots x_{N},t)\,.}$

For one particle in three dimensions, the Hamiltonian is:

${\displaystyle {\hat {H}}={\frac {{\hat {\mathbf {p} }}\cdot {\hat {\mathbf {p} }}}{2m}}+V(\mathbf {r} ,t)\,,\quad {\hat {\mathbf {p} }}=-i\hbar \nabla }$

generating the equation:

${\displaystyle i\hbar {\frac {\partial }{\partial t}}\Psi (\mathbf {r} ,t)=-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\Psi (\mathbf {r} ,t)+V(\mathbf {r} ,t)\Psi (\mathbf {r} ,t)}$

For N particles in three dimensions, the Hamiltonian is:

${\displaystyle {\hat {H}}=\sum _{n=1}^{N}{\frac {{\hat {\mathbf {p} }}_{n}\cdot {\hat {\mathbf {p} }}_{n}}{2m_{n}}}+V(\mathbf {r} _{1},\mathbf {r} _{2},\cdots \mathbf {r} _{N},t)\,,\quad {\hat {\mathbf {p} }}_{n}=-i\hbar \nabla _{n}}$

where the position of particle n is r_n, generating the equation:

${\displaystyle i\hbar {\frac {\partial }{\partial t}}\Psi (\mathbf {r} _{1},\mathbf {r} _{2},\cdots \mathbf {r} _{N},t)=-{\frac {\hbar ^{2}}{2}}\sum _{n=1}^{N}{\frac {1}{m_{n}}}\nabla _{n}^{2}\Psi (\mathbf {r} _{1},\mathbf {r} _{2},\cdots \mathbf {r} _{N},t)+V(\mathbf {r} _{1},\mathbf {r} _{2},\cdots \mathbf {r} _{N},t)\Psi (\mathbf {r} _{1},\mathbf {r} _{2},\cdots \mathbf {r} _{N},t)}$

This last equation is in a very high dimension, so the solutions are not easy to visualize.

src: i.ytimg.com

Solution methods

src: i.stack.imgur.com

Properties

The Schrödinger equation has the following properties: some are useful, but there are shortcomings. Ultimately, these properties arise from the Hamiltonian used, and the solutions to the equation.

Linearity

In the development above, the Schrödinger equation was made to be linear for generality, though this has other implications. If two wave functions ?₁ and ?₂ are solutions, then so is any linear combination of the two:

${\displaystyle \displaystyle \psi =a\psi _{1}+b\psi _{2}}$

where a and b are any complex numbers (the sum can be extended for any number of wavefunctions). This property allows superpositions of quantum states to be solutions of the Schrödinger equation. Even more generally, it holds that a general solution to the Schrödinger equation can be found by taking a weighted sum over all single state solutions achievable. For example, consider a wave function ?(x, t) such that the wave function is a product of two functions: one time independent, and one time dependent. If states of definite energy found using the time independent Schrödinger equation are given by ?_E(x) with amplitude A_n and time dependent phase factor is given by

${\displaystyle e^{{-iE_{n}t}/\hbar },}$

then a valid general solution is

${\displaystyle \displaystyle \Psi (x,t)=\sum \limits _{n}A_{n}\psi _{E_{n}}(x)e^{{-iE_{n}t}/\hbar }.}$

Additionally, the ability to scale solutions allows one to solve for a wave function without normalizing it first. If one has a set of normalized solutions ?_n, then

${\displaystyle \displaystyle \Psi =\sum \limits _{n}A_{n}\psi _{n}}$

can be normalized by ensuring that

${\displaystyle \displaystyle \sum \limits _{n}|A_{n}|^{2}=1.}$

This is much more convenient than having to verify that

${\displaystyle \displaystyle \int \limits _{-\infty }^{\infty }|\Psi (x)|^{2}\,dx=\int \limits _{-\infty }^{\infty }\Psi (x)\Psi ^{*}(x)\,dx=1.}$

Momentum space Schrödinger equation

The Schrödinger equation ${\displaystyle i\hbar \partial _{t}|\psi \rangle ={\hat {H}}|\psi \rangle }$ is often presented in the position basis form ${\displaystyle i\hbar \partial _{t}\psi (x)=-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi (x)+V(x)\psi (x)}$ (with ${\displaystyle \langle x|\psi \rangle \equiv \psi (x)}$). But as a vector operator equation it has a valid representation in any arbitrary complete basis of kets in Hilbert space. For example, in the momentum space basis the equation reads

${\displaystyle \displaystyle i\hbar \partial _{t}f(p)={\frac {p^{2}}{2m}}f(p)+({\tilde {V}}*f)(p)}$

where ${\displaystyle |p\rangle }$ is the plane wave state of definite momentum ${\displaystyle p}$, ${\displaystyle \langle p|\psi \rangle \equiv f(p)=\int \psi (x)e^{-ipx}dx}$, ${\displaystyle {\tilde {V}}}$ is the Fourier transform of ${\displaystyle V}$, and ${\displaystyle *}$ denotes Fourier convolution.

In the 1D example with absence of a potential, ${\displaystyle {\tilde {V}}=0}$ (or similarly ${\displaystyle {\tilde {V}}(k)\propto \delta (k)}$ in the case of a background potential constant throughout space), each stationary state of energy ${\displaystyle \hbar \omega =q^{2}/2m}$ is of the form

${\displaystyle f_{q}(p)=(c_{+}\delta (p-q)+c_{-}\delta (p+q))e^{-i\omega t}}$

for arbitrary complex coefficients ${\displaystyle c_{\pm }}$. Such a wave function, as expected in free space, is a superposition of plane waves moving right and left with momenta ${\displaystyle \pm q}$; upon momentum measurement the state would collapse to one of definite momentum ${\displaystyle \pm q}$ with probability ${\displaystyle \propto |c_{\pm }|^{2}}$.

A version of the momentum space Schrödinger equation is often used in solid-state physics, as Bloch's theorem ensures the periodic crystal lattice potential couples ${\displaystyle f(p)}$ with ${\displaystyle f(p+K)}$ for only discrete reciprocal lattice vectors ${\displaystyle K}$. This makes it convenient to solve the momentum space Schrödinger equation at each point in the Brillouin zone independently of the other points in the Brillouin zone.

Real energy eigenstates

For the time-independent equation, an additional feature of linearity follows: if two wave functions ?₁ and ?₂ are solutions to the time-independent equation with the same energy E, then so is any linear combination:

${\displaystyle {\hat {H}}(a\psi _{1}+b\psi _{2})=a{\hat {H}}\psi _{1}+b{\hat {H}}\psi _{2}=E(a\psi _{1}+b\psi _{2}).}$

Two different solutions with the same energy are called degenerate.

In an arbitrary potential, if a wave function ? solves the time-independent equation, so does its complex conjugate, denoted ?*. By taking linear combinations, the real and imaginary parts of ? are each solutions. If there is no degeneracy they can only differ by a factor.

In the time-dependent equation, complex conjugate waves move in opposite directions. If ?(x, t) is one solution, then so is ?*(x, -t). The symmetry of complex conjugation is called time-reversal symmetry.

Space and time derivatives

The Schrödinger equation is first order in time and second in space, which describes the time evolution of a quantum state (meaning it determines the future amplitude from the present).

Explicitly for one particle in 3-dimensional Cartesian coordinates - the equation is

${\displaystyle i\hbar {\partial \Psi \over \partial t}=-{\hbar ^{2} \over 2m}\left({\partial ^{2}\Psi \over \partial x^{2}}+{\partial ^{2}\Psi \over \partial y^{2}}+{\partial ^{2}\Psi \over \partial z^{2}}\right)+V(x,y,z,t)\Psi .\,\!}$

The first time partial derivative implies the initial value (at t = 0) of the wavefunction

${\displaystyle \Psi (x,y,z,0)\,\!}$

is an arbitrary constant. Likewise - the second order derivatives with respect to space implies the wavefunction and its first order spatial derivatives

${\displaystyle {\begin{aligned}&\Psi (x_{b},y_{b},z_{b},t)\\&{\frac {\partial }{\partial x}}\Psi (x_{b},y_{b},z_{b},t)\quad {\frac {\partial }{\partial y}}\Psi (x_{b},y_{b},z_{b},t)\quad {\frac {\partial }{\partial z}}\Psi (x_{b},y_{b},z_{b},t)\end{aligned}}\,\!}$

are all arbitrary constants at a given set of points, where x_b, y_b, z_b are a set of points describing boundary b (derivatives are evaluated at the boundaries). Typically there are one or two boundaries, such as the step potential and particle in a box respectively.

As the first order derivatives are arbitrary, the wavefunction can be a continuously differentiable function of space, since at any boundary the gradient of the wavefunction can be matched.

On the contrary, wave equations in physics are usually second order in time, notable are the family of classical wave equations and the quantum Klein-Gordon equation.

Local conservation of probability

The Schrödinger equation is consistent with probability conservation. Multiplying the Schrödinger equation on the right by the complex conjugate wavefunction, and multiplying the wavefunction to the left of the complex conjugate of the Schrödinger equation, and subtracting, gives the continuity equation for probability:

${\displaystyle {\partial \over \partial t}\rho \left(\mathbf {r} ,t\right)+\nabla \cdot \mathbf {j} =0,}$

where

${\displaystyle \rho =|\Psi |^{2}=\Psi ^{*}(\mathbf {r} ,t)\Psi (\mathbf {r} ,t)\,\!}$

is the probability density (probability per unit volume, * denotes complex conjugate), and

${\displaystyle \mathbf {j} ={1 \over 2m}\left(\Psi ^{*}{\hat {\mathbf {p} }}\Psi -\Psi {\hat {\mathbf {p} }}\Psi ^{*}\right)\,\!}$

is the probability current (flow per unit area).

Hence predictions from the Schrödinger equation do not violate probability conservation.

Positive energy

If the potential is bounded from below, meaning there is a minimum value of potential energy, the eigenfunctions of the Schrödinger equation have energy which is also bounded from below. This can be seen most easily by using the variational principle, as follows. (See also below).

For any linear operator  bounded from below, the eigenvector with the smallest eigenvalue is the vector ? that minimizes the quantity

${\displaystyle \langle \psi |{\hat {A}}|\psi \rangle }$

over all ? which are normalized. In this way, the smallest eigenvalue is expressed through the variational principle. For the Schrödinger Hamiltonian ? bounded from below, the smallest eigenvalue is called the ground state energy. That energy is the minimum value of

${\displaystyle \langle \psi |{\hat {H}}|\psi \rangle =\int \psi ^{*}(\mathbf {r} )\left[-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi (\mathbf {r} )+V(\mathbf {r} )\psi (\mathbf {r} )\right]d^{3}\mathbf {r} =\int \left[{\frac {\hbar ^{2}}{2m}}|\nabla \psi |^{2}+V(\mathbf {r} )|\psi |^{2}\right]d^{3}\mathbf {r} =\langle {\hat {H}}\rangle }$

(using integration by parts). Due to the complex modulus of ?² (which is positive definite), the right hand side is always greater than the lowest value of V(x). In particular, the ground state energy is positive when V(x) is everywhere positive.

For potentials which are bounded below and are not infinite over a region, there is a ground state which minimizes the integral above. This lowest energy wavefunction is real and positive definite - meaning the wavefunction can increase and decrease, but is positive for all positions. It physically cannot be negative: if it were, smoothing out the bends at the sign change (to minimize the wavefunction) rapidly reduces the gradient contribution to the integral and hence the kinetic energy, while the potential energy changes linearly and less quickly. The kinetic and potential energy are both changing at different rates, so the total energy is not constant, which can't happen (conservation). The solutions are consistent with Schrödinger equation if this wavefunction is positive definite.

The lack of sign changes also shows that the ground state is nondegenerate, since if there were two ground states with common energy E, not proportional to each other, there would be a linear combination of the two that would also be a ground state resulting in a zero solution.

Analytic continuation to diffusion

The above properties (positive definiteness of energy) allow the analytic continuation of the Schrödinger equation to be identified as a stochastic process. This can be interpreted as the Huygens-Fresnel principle applied to De Broglie waves; the spreading wavefronts are diffusive probability amplitudes. For a free particle (not subject to a potential) in a random walk, substituting ? = it into the time-dependent Schrödinger equation gives:

${\displaystyle {\partial \over \partial \tau }X(\mathbf {r} ,\tau )={\frac {\hbar }{2m}}\nabla ^{2}X(\mathbf {r} ,\tau )\,,\quad X(\mathbf {r} ,\tau )=\Psi (\mathbf {r} ,\tau /i)}$

which has the same form as the diffusion equation, with diffusion coefficient ?/2m. In that case, the diffusivity yields the De Broglie relation in accordance with the Markov process.

Regularity

On the space ${\displaystyle L^{2}}$ of square-integrable densities, the Schrödinger semigroup ${\displaystyle e^{it{\hat {H}}}}$ is a unitary evolution, and therefore surjective. The flows satisfy the Schrödinger equation ${\displaystyle i\partial _{t}u={\widehat {H}}u}$, where the derivative is taken in the distribution sense. However, since ${\displaystyle {\widehat {H}}}$ for most physically reasonable Hamiltonians (e.g., the Laplace operator, possibly modified by a potential) is unbounded in ${\displaystyle L^{2}}$, this shows that the semigroup flows lack Sobolev regularity in general. Instead, solutions of the Schrödinger equation satisfy a Strichartz estimate.

src: i.ytimg.com

Relativistic quantum mechanics

Relativistic quantum mechanics is obtained where quantum mechanics and special relativity simultaneously apply. In general, one wishes to build relativistic wave equations from the relativistic energy-momentum relation

${\displaystyle E^{2}=(pc)^{2}+(m_{0}c^{2})^{2}\,,}$

instead of classical energy equations. The Klein-Gordon equation and the Dirac equation are two such equations. The Klein-Gordon equation,

${\displaystyle {\frac {1}{c^{2}}}{\frac {\partial ^{2}}{\partial t^{2}}}\psi -\nabla ^{2}\psi +{\frac {m^{2}c^{2}}{\hbar ^{2}}}\psi =0.}$,

was the first such equation to be obtained, even before the nonrelativistic one, and applies to massive spinless particles. The Dirac equation arose from taking the "square root" of the Klein-Gordon equation by factorizing the entire relativistic wave operator into a product of two operators - one of these is the operator for the entire Dirac equation. Entire Dirac equation:

${\displaystyle \left(\beta mc^{2}+c\left(\sum _{n\mathop {=} 1}^{3}\alpha _{n}p_{n}\right)\right)\psi =i\hbar {\frac {\partial \psi }{\partial t}}}$

The general form of the Schrödinger equation remains true in relativity, but the Hamiltonian is less obvious. For example, the Dirac Hamiltonian for a particle of mass m and electric charge q in an electromagnetic field (described by the electromagnetic potentials ? and A) is:

${\displaystyle {\hat {H}}_{\text{Dirac}}=\gamma ^{0}\left[c{\boldsymbol {\gamma }}\cdot \left({\hat {\mathbf {p} }}-q\mathbf {A} \right)+mc^{2}+\gamma ^{0}q\phi \right]\,,}$

in which the ? = (?¹, ?², ?³) and ?⁰ are the Dirac gamma matrices related to the spin of the particle. The Dirac equation is true for all spin- ¹/₂ particles, and the solutions to the equation are 4-component spinor fields with two components corresponding to the particle and the other two for the antiparticle.

For the Klein-Gordon equation, the general form of the Schrödinger equation is inconvenient to use, and in practice the Hamiltonian is not expressed in an analogous way to the Dirac Hamiltonian. The equations for relativistic quantum fields can be obtained in other ways, such as starting from a Lagrangian density and using the Euler-Lagrange equations for fields, or use the representation theory of the Lorentz group in which certain representations can be used to fix the equation for a free particle of given spin (and mass).

In general, the Hamiltonian to be substituted in the general Schrödinger equation is not just a function of the position and momentum operators (and possibly time), but also of spin matrices. Also, the solutions to a relativistic wave equation, for a massive particle of spin s, are complex-valued 2(2s + 1)-component spinor fields.

src: nanohub.org

Quantum field theory

The general equation is also valid and used in quantum field theory, both in relativistic and nonrelativistic situations. However, the solution ? is no longer interpreted as a "wave", but should be interpreted as an operator acting on states existing in a Fock space.

src: i.ytimg.com

First order form

The Schrödinger equation can also be derived from a first order form similar to the manner in which the Klein-Gordon equation can be derived from the Dirac equation. In 1D the first order equation is given by

This equation allows for the inclusion of spin in nonrelativistic quantum mechanics. Squaring the above equation yields the Schrödinger equation in 1D. The matrices ${\displaystyle \eta }$ obey the following properties

The 3 dimensional version of the equation is given by

Here ${\displaystyle \eta =(\gamma _{0}+i\gamma _{5})/{\sqrt {2}}}$ is a ${\displaystyle 4\times 4}$ nilpotent matrix and ${\displaystyle \gamma _{i}}$ are the Dirac gamma matrices (${\displaystyle i=1,2,3}$). The Schrödinger equation in 3D can be obtained by squaring the above equation. In the nonrelativistic limit ${\displaystyle E-m\simeq E'}$ and ${\displaystyle E+m\simeq 2m}$, the above equation can be derived from the Dirac equation.

src: nanohub.org

See also

src: i.ytimg.com

Notes

src: nanohub.org

References

• P. A. M. Dirac (1958). The Principles of Quantum Mechanics (4th ed.). Oxford University Press. 
• B.H. Bransden & C.J. Joachain (2000). Quantum Mechanics (2nd ed.). Prentice Hall PTR. ISBN 0-582-35691-1. 
• David J. Griffiths (2004). Introduction to Quantum Mechanics (2nd ed.). Benjamin Cummings. ISBN 0-13-124405-1. 
• Richard Liboff (2002). Introductory Quantum Mechanics (4th ed.). Addison Wesley. ISBN 0-8053-8714-5. 
• David Halliday (2007). Fundamentals of Physics (8th ed.). Wiley. ISBN 0-471-15950-6. 
• Serway, Moses, and Moyer (2004). Modern Physics (3rd ed.). Brooks Cole. ISBN 0-534-49340-8. CS1 maint: Multiple names: authors list (link)
• Schrödinger, Erwin (December 1926). "An Undulatory Theory of the Mechanics of Atoms and Molecules". Phys. Rev. 28 (6): 1049-1070. Bibcode:1926PhRv...28.1049S. doi:10.1103/PhysRev.28.1049. 
• Teschl, Gerald (2009). Mathematical Methods in Quantum Mechanics; With Applications to Schrödinger Operators. Providence: American Mathematical Society. ISBN 978-0-8218-4660-5. 

src: i.ytimg.com

External links

• Hazewinkel, Michiel, ed. (2001) [1994], "Schrödinger equation", Encyclopedia of Mathematics, Springer Science+Business Media B.V. / Kluwer Academic Publishers, ISBN 978-1-55608-010-4 
• Quantum Physics - textbook by Benjamin Crowell with a treatment of the time-independent Schrödinger equation
• Linear Schrödinger Equation at EqWorld: The World of Mathematical Equations.
• Nonlinear Schrödinger Equation at EqWorld: The World of Mathematical Equations.
• The Schrödinger Equation in One Dimension as well as the directory of the book.
• All about 3D Schrödinger Equation
• Mathematical aspects of Schrödinger equations are discussed on the Dispersive PDE Wiki.
• Web-Schrödinger: Interactive solution of the 2D time-dependent and stationary Schrödinger equation
• An alternate reasoning behind the Schrödinger Equation
• Online software-Periodic Potential Lab Solves the time-independent Schrödinger equation for arbitrary periodic potentials.
• What Do You Do With a Wavefunction?
• The Young Double-Slit Experiment

Source of the article : Wikipedia