exact double well schrodinger

7/31/2019 Exact Double Well Schrodinger

1/20

The Chemical Educator

Exact solutions of the quantum double square well potential

Enrique Peacock-Lopez1,

1Department of ChemistryWilliams College

Williamstown, MA 01267(Dated: June 13, 2006)

For a symmetrical quantum double square well potential, we find analytical expressions satisfied bythe quantized energies. Graphical or numerical solutions of the former relations allow us to calculatenormalization constants and construct the first eight solutions of the Schrodinger equation. Withthese exact solutions, we analyze quantum tunneling across a potential barrier and compare ourresults with the experimental data for ammonia.

INTRODUCTION

Most chemistry textbooks [1]-[15] in their quantummechanics sections discuss, to different levels of analysis,the one-dimensional particle in a box (PIB) as a relevantone-dimensional quantum system. A natural applicationof the PIB considers delocalized electrons in hydro-

carbons or other organic compounds. In contrast withphysicists [16]-[22], who usually consider other relevantone-dimensional potential related to scattering, chemistsrarely consider other simple but chemically relevant one-dimensional potentials. In some textbooks [5]-[15] wefind a qualitative discussion of quantum tunneling andbarrier penetration but rarely a quantitative analysis .

Some advanced texts [5]-[15] may include a qualitativediscussion of barrier penetration in two minima or doublewell potentials, but again no quantitative analysis. In achemical system, the two minima potential corresponds,usually, to two equilibrium positions or molecular con-

formations. In the case of ammonia its so-called periodicinversion is a well documented example of quantum tun-neling.

The energy level splitting resulting from barrier pene-tration is an important quantum mechanical effect thatoccurs whether the double well potential is symmetricor not. The internal rotation in CH3CH3 from onestaggered configuration to another through an eclipsedpositions at the top of the barrier is a good exampleof a symmetric potential. In contrast, in the cases ofthe asymmetric H2O2 hindered rotor or beryllium di-cyclopentadienyl we encounter asymmetric double wellpotentials,. Ammonia (N H3), cyanamide (N H2CN),

P H3, and AsH3 are chemical examples of symmetricaldouble well potentials where atoms can tunnel throughthe barrier. In the case of ammonia the double minimarepresents the possible nitrogen positions with respect toplane defined by the three hydrogens.

In the case of ammonia and similar molecules the vi-brational level ting due to barrier penetration has beenobserved through infrared (IR) spectroscopy. The am-monia double well potential can be approximated usingthe Manning potential [23] where the the depth or dis-

sociation energy is estimated to be about 5 eV and theheight of the central barrier to be 0.25 eV. The quanti-tative analysis of the Manning [23] or similar potential[24] is quite difficult and only numerically tractable. Asan alternative, some authors [7, 15, 18, 19, 21] consider aone-dimensional discontinuous infinite depth double wellpotential to emphasize the symmetry of the wave func-

tion and quantum tunneling, but none consider a quan-titative analysis or a finite depth well..

In this paper we consider a double well finite depthpotential as an approximation of the ammonia poten-tial. In the second section we consider the discontinu-ous potential and its more general analytical solutions,as well as the continuity and smoothness conditions im-pose on the solutions of the Schrodinger equation by thequantum mechanical postulates. In the third section, weconsider the boundary conditions that yield the allowedquantized energies, and we calculate step by step the re-sultant transcendental equation for the case of infinitedepth and energy greater that the barrier hight, All of

the algebraic manipulations are straight forward and ac-cessible to chemistry juniors with the typical mathemat-ical background. In section four and five we calculatenumerically the allowed energies and construct the firsteight normalized wave function. In section six we dis-cuss quantum tunneling. Direct comparison between ourresults and experimental data from ammonia is considerin section seven. Finally we summarize and offer somesuggestion of how to use this paper in the classroom.

GENERAL SOLUTION

By dividing into five spatial regions, we set general dif-ferential equations derived from the Schrodinger equation(SE) and construct general solutions. The Double SquareWell Potential (DSWP) is a simplified model of poten-tials found in molecular chemical systems where possibledifferent conformations are separated by an energy bar-rier. Although the DSWP is a discontinuous function ofposition, it gives an appropriate and tractable descrip-tion of the continuous double minima potentials. TheDSWP is defined by the following piecewise function:


2/20

2

V(x) =

VD (L +a2 ) x <

0 a2 < x < (L +a2 )

Vo a2 x

a2

0 (L + a2 ) < x < a2

VD < x (L +a2 )

(1)

where VD is the depth or dissociation energy and Vo theheight of the central barrier separating two equivalentspatial regions. Also notice that the potential is an evenfunction of position, i.e., V(x) = V(x). Due to thissymmetry, the solutions of the SE are either odd or evenfunctions of position.

3 2 1 0 1 2 3

1

2

3

4

5Eepb vs xL

L RIII II I

FIG. 1: The Double Well Square Potential diagram, whereregions L, R, I, II, III are defined

Since the potential rises to the dissociation energy, VD,for |x| L + a/2 the SE,

2

2m

d2 RLd x2

+ VD RL = E RL (2a)

for E < VD can be reduced to

d2 RLd x2

= 2 RL (2b)

with

=

8m2(VD E)

h2(2c)

Now we consider the three internal spatial regions.First, we set the SE in each of the three spatial regions.

In region I, defined by a/2 < x < L + a/2, we have

2

2m

d2 Id x2

= E I (3a)

In region II, |x| a/2, we need to include the constantpotential Vo,

2

2m d

2

IId x2 + Vo II = E II (3b)

Finally, in region III, (L + a/2) < x < a/2,

2

2m

d2 IIId x2

= E III (3c)

For regions I and III, we can reduce the SE to a simplergeneral equation

d2

d x2= 2 (4a)

where we have defined

=

8mE2

h2(4b)

In region II we have two cases. In the first case the SEyields the following general equation:

d2 IId x2

= > 2 >II (6a)

with

>

8m2(E Vo)

h2(6b)

Since , > and < are positive real numbers, the solu-tions to eqs 4a, 6a are the simple and well known sine

and cosine functions. In the case of eq 5a, the solutionsare the hyperbolic sine and cosine, i.e., sinh and cosh.

The most general solution of the SE in each regionis a linear combination of sine and cosine or sinh andcosh functions. Consequently we have six unknown co-efficients that must be determined. The main problemin many of the one-dimensional quantum potentials isnot the solution of the differential equations but find-ing coefficients that satisfy the required continuity andsmoothness conditions at the regions boundaries. In


3/20


4/20

4

{En}. These En are consistent with the the continuityand smoothness condition at the potential boundaries.

The most familiar boundary condition requires, as inthe case of the particle in a box, that the wave func-tion vanishes at the edges of the box, and, given thesymmetry of the potential, we only have to consider oneedge,

I(L + a/2 ) = R (L + a/2) (14)

Using eq 10a, 12a we get

A cos((L + a/2)) + B sin((L + a/2)) =

DR exp ( (L + a/2))(15a)

The slopes of the s yield a second condition

A sin((L + a/2)) B cos((L + a/2)) =

DR exp ( (L + a/2))

(15b)

Now we use eq 15a in eq 15b and rearrange the

A

1 +

B

A

cos ((L + a/2)) =

B

1

A

B

sin ((L + a/2))

(16)

Equation 14 is satisfied in all possible cases, and wecan recast it as

tan ( (L + a/2)) = A

BF

A,B,

. (17)

where we define

F

A,B,

1 +

B

A

1 AB

(18)Equation 17 is a transcendental equation that we need

to solve for the quantized energies, which are a conse-quence of the spatial constraint posed on the system.The exact solutions of eq 17 can only be obtained nu-merically, and we can do so in many different ways. Inour approach we first rewrite eq 17 as

(L + a/2) = n + arctan

A

BF

A,B,

(19)

As we will show, once we fix the values of the mass,

the length of the box, the thickness and the height ofthe barrier, and the depth of the potential, , , A andB depend on the energy, E, and only a finite number ofenergies satisfy the relation expressed by equation eq 19.

Before proceeding with algebraic manipulations andfinding expressions for the coefficients A and B, we rescalethe spatial length and energy as follows:

x =

8mL2 Vo

h2

E

Vo

x

L c X (20a)

x

D2 c2 X (20b)

D

VDVo

(20c)

< x =

1 c2 X (20d)

> x =

c2 1 X (20e)

r a

2 L(20f)

In eq 20a, 2 represent the height of the barrier in unitsof epb defined as

epb h2

8mL2, (21)

which is determined by the the mass of the particle andthe length of the box. Also c2 represents the ratio of the

systems energy to the height of the barrier, and finallyX is the distance in units of L.

Now we are ready to consider the even and odd solu-tions for E > Vo and E < Vo and to analyze the boundaryconditions expressed by eqs 7. First we consider E > Voor c > 1, and we find that the even solution at x = a/2satisfies the following continuity condition:

AS cos ( c r) + BS sin ( c r)

= CS cos

c2 1 r

(22a)

In the case of the smoothness condition we get

AS sin ( c r) BS cos ( c r)

= CS

1 1

c2sin

c2 1 r

(22b)

We can solve for A if we first multiply eq 22a by cos(cr)and eq 22b by sin(cr). Second we add the resultingequations, and since sin2 + cos2 = 1 we find that As isgiven by the following expression:

AS = CS

cos(cr) cos

c2 1 r

+1

1

c2sin(cr) sin

c2 1 r

(23a)

Now we can multiply Eq.(22a) by sin and we substracteq 22b multiplied by cos. So we get

BS = CS

sin(cr) cos

c2 1 r

1

1

c2cos(cr) sin

c2 1 r (23b)

Notice that both A and B are expressed as function ofCs that can be calculated numerically using the required


5/20


6/20


7/20


8/20


9/20


10/20


11/20


12/20


13/20


14/20


15/20


16/20


17/20


18/20


19/20


20/20

exact double well schrodinger

Documents