2.

3: The One-Dimensional Particle in a Box

Imagine a particle of mass m constrained to travel back and forth in a one dimensional box of length a . For convenience, we define
the endpoints of the box to be located at x = 0 and x = a . The derivation of wavefunctions and energy levels and the properties of
the system using the tools of quantum mechanics will be instructive as we move forward in our studies of quantum mechanics.

The Hamiltonian
Whenever we begin a new quantum mechanical problem, the first challenge is to write the Hamiltonian that describes the system.
This always has two parts - a Kinetic Energy term (which is always the same for each particle) and a Potential Energy term (that
is different for each new system.)
The kinetic energy term in one dimension for a single particle is always given by
2 2
ℏ d
^ =−
T
2
2m dx

This operator can be derived from the momentum operator based on the relationship between momentum and kinetic energy that
comes from classical physics. Namely
2
p
T =
2m

As such,
2
p
^
T =
2m

1 d
= (−iℏ )
2m dx

2 2
(−iℏ) d
=
2
2m dx
2 2
ℏ d
=−
2
2m dx

The potential energy function is also fairly simple for this problem. The potential energy is infinite outside of the box (x < 0 and
x > a) and zero every place else. This forces the particle to be in the box at all times. It also limits the relevant space of the

problem to lie between x = 0 and x = a since the infinite potential energy precludes the particle from ever existing outside of the
limits of x = 0 and x = a .

∞ if x <0

U (x) =∣ 0 if 0 ≤x ≤a

∞ if x >a

Figure 2.3.1
So for the problem, limited to the space inside the box, the Hamiltonian can be written
2 2
ℏ d
^
H =−
2m dx2

And the Schrödinger equation can be written as

2 2
ℏ d
− ψ(x) = Eψ(x)
2
2m dx

2.3.1 https://chem.libretexts.org/@go/page/420475
where ψ(x) is the wavefunction describing the state of the particle. There are a number of approaches that can be used to solve this
equation to find the wavefunctions ψ(x) which satisfy the differential equation.

The Solution
We will solve this problem two different ways. First, we will solve it using the de Broglie wavelength (an algebraic solution) and
then using the Schrödinger equation (an eigenvalue/eigenfunction approach.)

The de Broglie Approach

Before trying to solve the problem using Schrödinger’s equation, let’s use the de Broglie condition to solve the problem
algebraically. Recall that de Broglie suggested that a particle can be treated as a wave, the wavelength of which is given by
λ = h/p , where h is Planck’s constant, and p is the momentum of the particle.

The necessary conditions on the de Broglie wave are that the wave itself must vanish at the ends of the box (in order to satisfy the
first postulate, since the particle can never escape the box.) This will happen for very specific wavelengths which are dependent on
the length of the box itself. This is very common in physics for any system with a wave nature. When the wave is constrained to a
specific geometry, the system will "ring" with frequencies (and thus wavelengths) characteristic of the medium and the geometry.
Quantum mechanical systems are no different in that regard.

Figure 2.3.2
What will be required in order to create a standing wave is that the length of the box (a) must be an integral multiple of half de
Broglie wavelengths (λ/2).
λ
a =n
2

Given that the de Broglie wavelength is related to momentum, it is simple to derive the following relationship, indicating the
possible values for momentum.
λ
a =n
2
nh
=
2p

nh
p =
2a

2.3.2 https://chem.libretexts.org/@go/page/420475
Given the relationship between momentum and kinetic energy, the expected expression for energy levels can be derived.
2 2 2 2
p 1 nh n h
E = = ( ) =
2m 2m 2a 8ma2

And since the energy depends on n , the spacings between successive energy levels increases as the energy increases.
2

Figure 2.3.3
Now let’s see if we can derive this expression based on the Schrödinger equation.

The Schrödinger equation: the wavefunctions

The time-independent Schrödinger equation can be written
^
H ψ = Eψ

Where H is the Hamiltonian operator that was derived in section B.2, ψ is the wavefunction describing the system, and E , the
eigenvalue of the Hamiltonian, gives the energy. The wavefunctions are derived so that they are eigenfunctions of the Hamiltonian
operator. Substituting the specific statement of the Hamiltonian
2 2
ℏ d
− ψ = Eψ
2
2m dx

For convenience, we can gather all of the constants in one place by making a substitution

2
2mE
−k =−
2
ℏ

The particular choice if the form of this substitution is made to simplify the solutions by avoiding (for now) imaginary functions.
With the substitution, the Schrödinger equation can be rewritten as
2
d
2
ψ = −k ψ
2
dx

As was the case for the classical wave-on-a-string problem, this is a second order ordinary differential equation, and this has two
linearly independent solutions. A general solution is given by a linear combination of two linearly independent solutions, so one
way to write a solution is

ψ = A sin(kx) + B cos(kx)

Now we can focus on evaluating A, B and k based on the boundary conditions. The boundary conditions are that the wavefunction
must go to 0 at the ends of the box, in accordance with the first postulate.
The first boundary condition, ψ(0) = 0 , yields the following result:

2.3.3 https://chem.libretexts.org/@go/page/420475
 ψ(0) = A sin(k ⋅ 0) + B cos(k ⋅ 0)

= 0 +B = 0

So B = 0 and the cosine term must vanish. Focusing only on what has not vanished from the solutions, the second boundary
condition, ψ(a) = 0 , can be applied.

ψ(a) = A sin(k ⋅ a) = 0

There are two trivial ways to make this true. One is to make A = 0 and the other is to make k = 0 . Both are trivial solutions and
unimportant (but fun to mention in class!) The other way to force the function to 0 at x = a is to insure that the sine function is
zero by forcing

k ⋅ a = nπ

where n is an integer (n = 1, 2, 3 …), since the sine function crosses zero every nπ radians. This is an important point: the
application of a boundary condition leads to the introduction of a quantum number and fixed the results to only functions where
that number has a value taken from a very specific list. In fact, the origin of quantum numbers in all problems is the result of the
application of boundary conditions.
Solving for k and substituting yields
nπx
ψ(x) = A sin( )
a

This is as far as the boundary conditions can get us. The value of A is determined based on the first postulate of quantum
mechanics, which says that the square of the wavefunction must give a probability distribution as to where the particle can be
measured to be. Since all measurements must place the particle in the box, the sum of probabilities at all of the possible locations in
the box must equal unity. This implies the condition that
a
2
∫ (ψ(x)) dx = 1
0

Solving for A yields

a a
nπx
2 2 2
∫ (ψ(x)) dx = A ∫ sin ( )dx
0 0
a
a
2nπx
⎡ ⎤
sin( )
⎢x a ⎥
2
=A ⎢ − ⎥
⎢ 2 4nπ ⎥
( )
⎣ ⎦
a
0

2
a
=A ( − 0 − 0 + 0)
2
−−
2
A =√
a

Notice that the value of A did not depend on the quantum number n. Normalization constants usually do have some dependence on
the quantum numbers that arise from the application of boundary conditions, but this is one of the rare problems in which the
normalization constant does not.

The Schrödinger Equation: the energy levels

Whenever we solve a quantum mechanical problem, there are two important things at which we must look: the energy levels and
the wavefunctions. To chemists, the energy levels are the most important part, as the energy levels govern the chemistry the system
can do. To a physicist, it is the wavefunctions that are important as they contain all of the information about the physical nature of
the system.
The energy levels can be derived using the normalized wavefunctions and the Schrödinger equation.

2.3.4 https://chem.libretexts.org/@go/page/420475
 ^
H ψ = Eψ

2 2
−
− 2
−
−
2
ℏ d 2 nπx ℏ nπ 2 nπx
− √ sin( ) = ( ) √ sin( )
2m dx2 a a 2m a a a
   
^ ψ E ψ
H

Comparison (or solving for E ) yields the following

2 2 2
n π ℏ
E =
2
2ma

which looks similar to, but not exactly like the result produced using the de Broglie relationship. In fact, it is the identical result!
Making the substitution ℏ = h/2π , it is easy to show that
2 2
n h
E =
8ma2

These energy levels depend on n , and so doubling the quantum number n quadruples the energy. Another way of saying this is
2

that the energy level spacings (the difference in energy between two successive levels) increase with increasing n or energy.
It is also interesting to note that the energy levels are given by a real (non-imaginary) expression. This is to be expected since the
energy is the eigenvalue of a Hermitian operator, the Hamiltonian, and thus must be a real value.

Properties of the Wavefunctions

The wavefunctions for the one-dimensional particle in a box problem are given by
−
−
2 nπx
ψn (x) = √ sin( )
a a

These wavefunctions have many important properties.

Orthogonality
Similar to the relationship of Hermitian operators having real eigenvalues, the eigenfunctions of Hermitian operators must be
orthogonal. Our wavefunctions are actually an infinite set of function, any pair of which must cause the inner product integral to
vanish. Mathematically, this looks like
a

∫ ψn (x) ⋅ ψm (x)dx = 0 n ≠m
0

This relationship is easy to verify. To do so, we will make use of the following result taken from a standard table of integrals.
sin[(α − β)x] sin[(α + β)x]
∫ sin(αx) sin(βx)dx = − α ≠β
2(α − β) 2(α + β)

nπ mπ
Noting that α = and β = , substitution into the above relationship yields
a a

π π a

a ⎡ sin[ (n − m)x] sin[ (n + m)x] ⎤

a a
∫ ψn (x) ⋅ ψm (x)dx = ⎢ − ⎥
π π
0
⎣ 2( (n − m)) 2( (n + m)) ⎦
a a 0

⎡ ⎤
sin[π(n − m)] sin[π(n + m)]
=⎢ − − 0 + 0⎥
π π
⎣ 2( (n − m)) 2( (n + m)) ⎦
a a

And since n and m are integer, n − m and n + m must also be integers. And the sine of an integral multiple of π is always zero, it
is easy to show that this function vanishes for any n ≠ m .

2.3.5 https://chem.libretexts.org/@go/page/420475
Normalization
When n = m the integral becomes
a a
2
2 nπx
2
∫ [ ψn (x)] dx = ∫ sin ( )dx
0
a 0
a

which can be evaluated using the result from a table of integrals

x sin(2αx)
2
∫ sin (αx)dx = −
2 4α

nπ
So making the substitution α =
a

a
2nπx
⎡ ⎤
sin( )
a
2 nπx 2 ⎢x a ⎥
2
∫ sin ( )dx = ⎢ − ⎥
a 0 a a ⎢ 2 4(nπ/a) ⎥

⎣ ⎦
0

2 a
= [ − 0 − 0 + 0]
a 2

=1

−
−
2
This result shouldn’t be surprising since the value A =√ was chosen to ensure the result! Specifically, it was chosen so as to
a

normalize the wave functions.

 Example 2.3.1

Show that the wavefunction

−−
−
30
Ψ(x) = √ ⋅ x(a − x)
5
a

is normalized for a particle in a box of length a .

Solution
The wavefunction is normalized if
a

∫ Ψ(x)Ψ(x)dx = 1
0

This can be demonstrated by plugging the wavefunction into the relationship and testing to see if it is true:
a −−
− −−
− a
30 30 30
2 2 2
∫ √ ⋅ x(a − x) √ ⋅ x(a − x)dx = ∫ x (a − 2ax + x ) dx
5 5 5
0 a a a 0
a
30
2 2 3 4
= ∫ (a x − 2ax + x ) dx
5
a 0
a
2 3 4 5
30 a x 2ax x
= [ − + ]
5
a 3 4 5
0

5 5 5
30 a a a
= ( − + − 0 + 0 − 0)
5
a 3 2 5
5 5 5
30 10a 15a 6a
= ( − + )
5
a 30 30 30

5
30 a
= ( )
5
a 30

=1

Therefore the wavefunction is normalized!

2.3.6 https://chem.libretexts.org/@go/page/420475
This page titled 2.3: The One-Dimensional Particle in a Box is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or
curated by Patrick Fleming.

2.3.7 https://chem.libretexts.org/@go/page/420475