Schrödinger equation

From Wikiversity
Jump to navigation Jump to search

The Schrödinger equation is an equation that is fundamental to quantum theory. The time independent Schrödinger equation looks like

The is called the Hamiltonian Hamiltonian mechanics. denote the wavefunctions.

How to construct the Schrödinger equation for a system[edit | edit source]

The operator is called the Hamiltonian. It contains two parts:

  1. The kinetic energy, denoted by the operator
  2. The potential energy, denoted by the operator

These two operators represent the principal types of energies in any physical system. Putting these two parts together, we can write this as . We can think of the Hamiltonian as a mathematical object that encodes how the energies in a system can be distributed.

Kinetic energy operator[edit | edit source]

In classical mechanics, we know that the kinetic energy of a moving particle is given by

in which corresponds to the mass of the particle, and is its velocity. Noting that classically the linear momentum is given by , so the kinetic energy can also be written as


In order to get the quantum mechanical version of this, we substitute all occurrences of the momentum with


where is the vector analogue of the partial differential operator .

Hence, the kinetic energy operator is given by


Total energy operator and the time dependent Schrödinger equation[edit | edit source]

The total energy operator introduces the time propagation to the wavefunctions so the time dependent Schrödinger equation tells how the quantum system developes is time while ist was found is some state at the original time. It is basically constructed from the condition that its operation on time dependent oscillatory factors characteristic to the plane waves obtained from the classical wave equation describing the mechanical waves like sound or the waves on the water surface should produce the energy quanta as the eigenvalues from the eigen-equation namely

One therefore may readily guess it as

While te total energy in the classical mechanics is normally equal to the Hamiltonian the time dependent Schrödinger equation is obtained by equalizing the two:

How to solve the Schrödinger equation[edit | edit source]

Example 1: The ground state of the hydrogen atom[edit | edit source]

The time independent Schrödinger equation for the electron in the Coulomb potential in the spherical coordinates takes the form:

For the ground state which is the state with the lowest possible we seek the solution in the simplest possible form not dependent on the spherical angles and . The Schrödinger equation greatly simplifies to the second order ordinary partial differential equation of the one radial variable:

The part of it contains derivatives multiplied by the various powers of the variable and we guess the solution in the exponential form since the exponential function is the eigenfunction of the differentiation operator while the derivative of the exponential function is proportional to it:

Substituting it to the Schrödinger equation we obtain the following equation:

The only way this equation is valid for all values of is that the coefficient multiplying the is equal to :

that gives the value of the Bohr radius (the range how the ground state mostly extents radially):

and the ground state energy


Example 2: Time dependent perturbation theory of weak field ionization[edit | edit source]

The first thing one can think to do above finding the ground state or the other allowed energies of the system is to solve the time dependent Schrödinger equation approximately i.e. in the perturbative manner. Let us apply it to the small destruction effects on the ground state of the quantum system and for the simplicity in one spatial dimension which reduces the kinetic energy operator to the second derivative with the proper coefficients. Let the ground state of the atom, for example the atom of Hydrogen will be such that is the smallest of possible values for the eigen-equation:

Let us add the external electric field effect to the Hamiltonian in the simplest linear manner so the Schrödinger equation is


We will assume for the initial condition and will check what happens to that state under the influence of the perturbing term, more specifically where it will go in terms of the expansion in the basis of the unbound states of the system which in very good approximation are simply plane waves. For the convenience of the calculations we added the envelope function normalized to in time so the electric field is turned on and off in the smooth and easy to Fourier transform Gaussian manner:


it is now convenient to absorb the term in the equation in the so-called interaction picture introducing the transformed solution:

The original time-dependent equation greatly simplifies to


where the position operator becomes time dependent namely

With the zero iteration condition we get in the first approximation

we are now interested in the unbound continuum state component in the final solution which tells us about the energy spectrum of the ionized electrons. We can get that by projecting the solution onto the continuum state, simply:

where we extended the integration in time from the minus infinity to infinity since the electric field was localized in time with the Gaussian envelope. The only significant term which remains after the integration is:

Where the is the Fourier transform of the field envelope function and while is broadly localized the is sharply loalized around the point with the condition

expressing the energy conservation condition

stating that while the ground state energy is negative there is no ionization below the threshold (for the photon or field frequencies that low that the kinetic energy would be impossibly negative) and the ejected electron kinetic energy is equal the energy of the absorbed photon lowered by the ionization or the binding energy.

Example 3: Theory of Above Threshold Ionization (ATI)[edit | edit source]

Above Threshold Ionization may be explained by solving the time dependent Schrödinger equation in the approximate manner. The Schrödinger equation for the free electron in the field of the electromagnetic wave in one dimension and in the radiation gauge is given by



so the electric field is given by


we are obtaining the equation for

with the solution




The Schrödinger equation for the electron in the wave field and in the atomic potential will be given by


where, is the Hamiltonian of the free electron in the field. Adding and subtracting the energy of the ground state from which the electron is to be ionized we obtain

Because in the ground state the electron kinetic energy is equall the total energy but with the opposite sign (virial theorem) and only this energy will be left after the fast removal of the electron, we neglect in this equation the sum for all and obtain the approximate equation


where the only remaining of the atomic potential is the constant.

This equation may now be solved using the previous result for the free electron in the field and expanding the ground state into its Fourier components



This equation has therefore the solution

We obtain the ionization spectrum from the formula


telling how much of the plane wave of the free electron with the given kinetic energy is at the and of the ionization process, where is the averaging function of the measuring detector for example


Expanding the factor

with the generalized Bessel functions defined by the inverse transform, we obtain

( is the Fourier transform of the detector function) therefore the sum of sharp or fuzzy maxima localized around the energy condition of emitted electrons depending on the speed i.e the averaging parameter of the detector .

Example 4: Basic theory of the Trojan Wave Packet[edit | edit source]

Because the Trojan wave packet is an approximate eigenstate of the Hamiltonian its search starts from the Schrödinger´s equation for the Hydrogen atom in the field of the circularly polarized wave in the coordinate system rotating with the frequency of the field




and is the operator of z-component of the angular momentum.

For the simplification one can restrict himself to only two spacial dimensions. This equation after neglecting the curvature terms of the Laplacian and assumption that some terms are slowly varying or even constant takes in the polar coordinate system the simplified and separated form


where , and is the point of the potential of the part of the Hamiltonian dependent only on the radial coordinate and therefore also the radius of the classical orbit. Because the radial potential has the minimum the Trojan wavepacket is in this theory the product of a well localized radial Gaussian function and a well localized excited state of the quantum mathematical pendulum, corresponding to the pendulum totally reversed upside down which is also the Mathieu function with the characteristic value (and the energy) equal with a good approximation to . The more precise theory shows that it is the pendulum of a negative and fractional electron mass equal -1/3 (of the electron) and that is why it is on the contrary the ground state with approximately the same energy but in the reversed spectrum.

To notice that let us more keep the correction from the expansion of the centrifugal term

and let us consider its action on the Fourier component of , e.g. .

The additional centrifugal potential will be

It will cause the shift of the minimum and the energy lowering of the radial harmonic oscillator equal to

Replacing back with the operator of the z component of the angular momentum we get

Adding this correction in the Schródinger equation we get the equation with the negative and fractional electron mass -1/3 ():


which allows the Gaussian localization of the electron in the maximum but not in the minimum of the potential in the angular direction.

References[edit | edit source]

[1] Qi-Chang Su and J. H. Eberly, Numerical simulations of multiphoton ionization and above-threshold electron spectra, Phys. Rev. A 38, 3430 (1988)

[2] H. R. Reiss and V. P. Krainov, Approximation for a Coulomb-Volkov solution in strong fields, Phys. Rev. A 50, R910 (1994)

Links[edit | edit source]

Simulations[edit | edit source]