In 1926, the Austrian physicist Erwin Schrödinger (1887-1961) made a fundamental mathematical discovery that had a profound impact on the study of the molecular world (in 1933, Schrödinger was awarded with the Nobel prize in Physics just 7 years later his breakthrough discovery). He discovered that a state of a quantum system composed by particles (such as electrons and nucleons) can be described by postulating the existence of a function of the particle coordinates and time, called state function or wave function (, psi function). This function are solution of a wave equation: the so-called the Schrödinger equation (SE). Although the SE equation can be solved analytically only for relatively simple cases, the development of computer and numerical methods has made possible the application of SE to study complex molecular.

Nowadays, computational quantum chemistry has surged to a level of a fundamental research discipline in chemistry and material science. In recognition of the contribution to the advancement of science, the pioneers of this discipline, Prof. Walter Kohn, and Prof. John A. Pople have been awarded the 2008 Nobel prize in chemistry.

In this primer, we are going to study analytical and numerical solutions of SE for simple quantum systems in stationary conditions, e.g. whose state does not change with the time. For this purpose, we are going to use the time-independent form of SE. This condition is an excellent approximation to describe the state of molecular systems.

The time-independent SE is given in his compact form by the equation.

where the (H hat) indicates the time-independent Hamiltonian operator used to calculate the total energy of the system. I just recall you that an operator is a mathematical transformation that applies to what follows the operator symbol, in this case, the wave function (). We will see in a moment how this looks. The application of the operator produce as results the same wave function multiplied for a number $E.$ This number represented the total energy state of the system (or \em eigenstate\em). This equation is also called *eigenvalue *equation, and, in this case, is also said to be an eigenfunction of with the associated \em eigenvalue,\em

Shall we look more in detail to the Hamiltonian operator, and, for the sake of simplicity, shall we limit ourselves to consider only a mono-dimensional physical system, namely a system defined only along the reference axis *x*. Therefore, the wave function depends on only by *x* or The Hamiltonian operator represents the total energy of the system, and, therefore, is composed by the summation of a kinetic and potential energy operator: The kinetic energy is expressed as

where is the classic momentum that is associated with the partial differential operator The potential energy is usually a function of the coordinates, and it is associated (in our case) with operator Therefore we can write the resulting SE as:

This differential equation can be solved analytically for different types of potentials. A simple classical example is a particle in a one-dimension box within infinitely repulsive barriers.

In this case, the value of is equal zero for and infinite in the walls. The wave function beyond the walls of the box is null. The wave function and energies of the particle inside the box are obtained by the solution of SE with reported in the central part of the Figure 1.

**Derivation of the wavefunction **

The one-dimensionSE is a linear second order differential equation with solutions

using the Euler relation

we can write

or called and we can write

by substituting in SE with (e.g. particle inside the box)

we obtain

that gives

the quantised total energy of the system.

Now we estimate the values of constant *k*, *C*, and *D* of the solution.

By applying the boundary conditions

we obtain for *x=0* that

It is more convenient to represent the parameter *k* as a circular frequency with the De Broglie wavelength

of the particle

For *x=L* the solution gives

This is accomplished for or

and

We still need to calculate the value of the constant *C*. As the square of the wave function gives the probability density and its integral the probability to find an particle in a region of the box then if we integrate over all the *x* axis the probability should be equal to one:

(1)

Using this integral we can find the value of C that *normalize* the function. The integral with can be calculated by part

as with in this case therefore

so we have that

taking the last integral in the second member of the equation to the first member we have

but the first term in the second member is equal to zero and the second equal to L, hence

.

Therefore from Equ. (1), we have

and .

Therefore the wavefunctions of a particle in a box are given by

.

and the associated energies are

.

The quantum number *n* is an integer associated with the different energy levels of the particle. Larger *n* values correspond to higher energy levels of the system. Also note in the graphs, the increase of the number of nodes in the wave function with the increase of the value of *n*.

In Figure 2, the wavefunctions for the first three energy level of the particle in a box are shown.

The separation between the adjacent energy level with quantum number $n$ and $n+1$ is

The probability density for a particle in a box is given by the square of the wavefunction

and the probability to find the particle in a region of the box can be calculated using the integral

.

**READINGS**

- Atkins, P. and Paula, J. (2010) Physical Chemistry. 9th Edition, W. H. Freeman Co., New York.
- I. Levine. Quantum Chemistry. VI edition. Pearson International edition.