Linear Harmonic Oscillator
Linear Harmonic Oscillator
Linear Harmonic Oscillator
−ħ2 𝑑2𝜓
+ 𝑉(𝑥)𝜓 = 𝐸𝜓
2𝑚 𝑑𝑥 2
2
𝑑 𝜓 2𝑚
𝑂𝑟, 2 + (𝐸 − 𝑉)𝜓 = 0 . . . . . . . . . . . . . . . (1)
𝑑𝑥 ħ2
Where E is the total energy, V is the potential energy and 𝜓 is the wave function for
the particle.
For a linear S.H.O. along the x- axis with angular frequency 𝜔 under a restoring
force proportional to the displacement 𝑥, the potential energy,
1
𝑉 = 𝑚 𝜔2 𝑥 2 ................... (2)
2
𝑑 2 𝜓 2𝑚 1
+ (𝐸 − 𝑚 𝜔2 𝑥 2 ) 𝜓 = 0
𝑑𝑥 2 ħ2 2
𝑑2𝜓 2𝑚 𝐸 𝑚2 𝜔2 2
𝑜𝑟, 2 + ( 2 − 𝑥 )𝜓
𝑑𝑥 ħ ħ2
= 0. . . . . . . . . . . . . . . . . . . . . . . . . . . . . (3)
This is the Schrodinger wave equation for the oscillator.
To simplify the above equation, we introduce a dimension less independent variable
y,
𝑚𝜔
𝑦 = 𝑎𝑥 ; 𝑊ℎ𝑒𝑟𝑒 𝑎 = √
ħ
𝑦 𝑑𝑦
∴𝑥= ∴ =𝑎
𝑎 𝑑𝑥
𝑑𝜓 𝑑𝜓 𝑑𝑦 𝑑𝜓
Now, = =𝑎
𝑑𝑥 𝑑𝑦 𝑑𝑥 𝑑𝑦
𝑑2𝜓 𝑑 2 𝜓 𝑑𝑦 2
𝑑2𝜓
𝐴𝑛𝑑 =𝑎 = 𝑎
𝑑𝑥 2 𝑑𝑦 2 𝑑𝑥 𝑑𝑦 2
𝑑2𝜓
Substituting the values of 𝑥 2 𝑎𝑛𝑑 in equation (3);
𝑑𝑦 2
𝑑2𝜓
2
2𝑚 𝐸
𝑎 + ( − 𝑎2 𝑦 2 ) 𝜓 = 0
𝑑𝑦 2 ħ2
𝑑2𝜓 2𝑚 𝐸
𝑂𝑟, 2 + ( 2 2 − 𝑦 2 ) 𝜓 = 0
𝑑𝑦 𝑎 ħ
2
𝑑 𝜓 2𝐸
𝑂𝑟, 2 + ( − 𝑦2) 𝜓 = 0
𝑑𝑦 ħ𝜔
2
𝑑 𝜓
𝑂𝑟, 2 + ( 𝜆 − 𝑦 2 ) 𝜓 = 0. . . . . . . . . . . . . . . . . . . . . . . . . (4)
𝑑𝑦
2𝐸
𝑊ℎ𝑒𝑟𝑒 𝜆 = . . . . . . . . . . . . . . . . . . . . (5)
ħ𝜔
Solution:- To solve the above equation let us try first an asymptotic solution i.e.
when
𝑦 >> 𝜆.
∴Equations (4)reduces to,
𝑑2𝜓
2
− 𝑦 2 𝜓 = 0 . . . . . . . . . . (6)
𝑑𝑦
Which has approximate solution,
𝑦2
−
𝜓= 𝑒 2 … … … … (7)
𝑦2
[since 𝜓 → 0 𝑤ℎ𝑒𝑛 𝑦 → 0 , 𝑤𝑒 𝑟𝑒𝑗𝑒𝑐𝑡 𝑡ℎ𝑒 𝑓𝑜𝑟𝑚 𝑒 2 ]
𝑑𝜓 𝑦2
−
∴ = −𝑦 𝑒 2 = −𝑦𝜓
𝑑𝑦
&,
𝑑2𝜓 𝑑𝜓
= − (𝜓 + 𝑦 )
𝑑𝑦 2 𝑑𝑦
= −(𝜓 − 𝑦 2 𝜓)
= (𝑦 2 − 1)𝜓
𝑑2𝜓
𝑜𝑟, 2 − (𝑦 2 − 1) 𝜓 = 0. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (8)
𝑑𝑦
For large values of 𝑦 equation (8) is reduced to equation (6). This suggests an
accurate solution of the equation (6) must be in the form,
𝑦2
−
𝜓 = 𝑓(𝑦)𝑒 ..................... (9)
2
𝑑𝜓 𝑦 2 𝑑𝑓 𝑦2
− −
= 𝑒 2 + 𝑓 𝑒 2 (−𝑦)
𝑑𝑦 𝑑𝑦
𝑑𝑓 𝑦2
−
= ( − 𝑓𝑦) 𝑒 2
𝑑𝑦
2
𝑑 𝜓 𝑑2𝑓 𝑑𝑓 −
𝑦2 𝑑𝑓 −
𝑦2
𝑜𝑟, 2 = ( 2 − 𝑦 − 𝑓 ) 𝑒 2 + ( − 𝑓𝑦) 𝑒 2 (−𝑦)
𝑑𝑦 𝑑𝑦 𝑑𝑦 𝑑𝑦
2 𝑦2
𝑑 𝑓 𝑑𝑓 2 −
= ( 2 − 2𝑦 − 𝑓 + 𝑓𝑦 ) 𝑒 2
𝑑𝑦 𝑑𝑦
∴ 𝐸𝑞 (4) 𝑟𝑒𝑑𝑢𝑐𝑒𝑠 𝑡𝑜,
𝑑2𝑓 𝑑𝑓 2 −
𝑦2
2 −
𝑦2
( 2 − 2𝑦 − 𝑓 + 𝑓𝑦 ) 𝑒 2 + ( 𝜆 − 𝑦 )𝑓 𝑒 2 = 0
𝑑𝑦 𝑑𝑦
𝑑2𝑓 𝑑𝑓
𝑂𝑟 , − 2𝑦 + ( 𝜆 − 1) 𝑓
𝑑𝑦 2 𝑑𝑦
= 0. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (10)
Replacing 𝜆 = 2𝑛 + 1 , (𝑛 = 0,1,2 , . . . . . . . . . . . . . . . ) Equation (10) becomes
𝑑2𝑓 𝑑𝑓
− 2𝑦 + 2𝑛 𝑓 = 0
𝑑𝑦 2 𝑑𝑦
This is well known Hermite’s differential equation. The solutions of the equation are
called Hermite’s Polynomials and are given by,
2 𝑑 𝑛 −𝑦2
𝐻𝑛 (𝑦) = (−1) 𝑛 𝑒 𝑦 𝑒 . . . . . . . . . . . . (11)
𝑑𝑦 𝑛
Eigen values :- For physical significance only those solutions of the Hermite eq (10)
are acceptable , for all values of 𝑦 for which ,
𝜆 = 2𝑛 + 1 [ n = 0,1,2,........called quantum
numbers]
2𝐸
Substituting , 𝜆 = ,
ħ𝜔
ħ𝜔
𝐸 = (2𝑛 + 1)
2
1
𝑂𝑟, 𝐸 = (𝑛 + )ħ𝜔
2
Or more generally,
1
𝐸𝑛 = (𝑛 + ) ħ𝜔 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (12)
2
From equation (12) we have the following conclusions,
1. The wave equations for the oscillator is satisfied only for discrete values of
total energies,
2. The lowest energy of the oscillator is obtained by putting n = 0 and it is,
1
𝐸0 = ħ𝜔
2
This is called the ground state energy or the zero point energy of the harmonic
oscillator.
(iii) The eigen values of the total energy depends only on one quantum number ‘n’ ;
all energy levels are hence non degenerate .
(iv) The successive energy levels are equally spaced; separation between two
adjacent energy level being ħ 𝜔0 .
Wave functions :- For each value of ‘n’ there is a different wave function 𝜓𝑛 which
have the general form ,
𝑦2
−
𝜓𝑛 (𝑦) = 𝑁𝑛 𝐻𝑛 (𝑦)𝑒 2
𝑎2 𝑥2
− 𝑚𝜔
Or, 𝜓𝑛 (𝑥) = 𝑁𝑛 𝐻𝑛 (𝑎𝑥)𝑒 2 where a = √
ħ
1 ∞ 2
Or ∫−∞ 𝑁𝑛2 𝑒 −𝑦 𝐻𝑛2 (𝑦)𝑑𝑦 = 1
𝑎
Using the relation of Hermit’s polynomial ,
∞ 2
∫−∞ 𝑒 −𝑦 𝐻𝑛 (𝑦)𝐻𝑚 (𝑦)𝑑𝑦 = 2𝑛 𝑛! √𝜋𝛿𝑚𝑛
𝑁𝑛2
2𝑛 𝑛! √𝜋 = 1
𝑎
1
𝑎 1
Or , 𝑁𝑛 = [ ] 2
2𝑛 𝑛! √𝜋
1 𝑦2
𝑎 1 −
∴ 𝜓𝑛 (𝑦) = [ ] 𝑒
2 2 𝐻𝑛 (𝑦).
2𝑛 𝑛! √𝜋
λ=2n+1 𝑬𝒏 𝑯𝒏 (𝒚) 𝟏 𝒚𝟐
𝒂 𝟏
𝝍𝒏 (𝒚) = [𝟐𝒏 𝒏! ] 𝟐 𝒆− 𝟐 𝑯𝒏 (𝒚).
√ 𝝅
1 1 𝐻0 (𝑦) = 1 1 𝑦2
𝑎
ℏ𝜔 Ψ0 (𝑦) = [ ] 2 𝑒− 2
2 √ 𝜋
3 1 𝐻1 (𝑦)=2y 𝑎 𝑦2
−
ℏ𝜔 Ψ1 (𝑦) = [ ] 2y𝑒 2
2 2√𝜋
5 1 𝐻2 (𝑦) = 4𝑦 2 − 1 𝑎 1
−
𝑦2
ℏ𝜔 Ψ2 (𝑦) = [ ] 2 (4𝑦 2 − 1)𝑒 2
2 8√𝜋
7 1 𝐻3 (𝑦) = 8𝑦 3 − 12𝑦 𝑎 1
ℏ𝜔 Ψ3 (𝑦) = [ ] 2 (8𝑦 3
2 48√𝜋
𝑦2
− 12𝑦)𝑒 − 2
9 1 𝐻4 (𝑦) 𝑎 1
ℏ𝜔 Ψ4 (𝑦) = [ ] 2 (16𝑦 4 − 48𝑦 2
2 = 16𝑦 4 − 48𝑦 2 384√𝜋
+ 12 𝑦2
−
+ 12) 𝑒 2
QuantizedSolutionsOfThe1DSchroedingerEquationForAHarmonicOsc.cdf
Q:Calculate probability density for the 1d LHO classically. Compare with quantum
result.
∴ 𝑑𝑥 = 𝐴𝜔 cos( 𝜔𝑡 + 𝜃) 𝑑𝑡
𝑑𝑥
∴ 𝑑𝑡 =
A𝜔 cos ( 𝜔𝑡 + 𝜃)
𝑑𝑥
=
𝐴𝜔√1 − 𝑠𝑖𝑛2 (𝜔𝑡 + 𝜃
𝑑𝑥
=
𝑥2
𝐴𝜔√1 −
𝐴2
𝑑𝑥
=
𝜔√𝐴2 − 𝑥 2
𝑑𝑥
∴ 𝑃(𝑥)𝑑𝑥 =
𝜋√𝐴2 − 𝑥 2
1
𝑖. 𝑒. 𝑃 (𝑥) = . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (1)
𝜋√𝐴2 − 𝑥 2
According to the quantum mechanics the probability – density,
−𝑎2 𝑥2
2 1
∴ 𝑃0 (𝑥) = |𝜓0 | = ( ) 𝑒 2
𝑥 √
QuantumClassicalCorrespondenceForTheHarmonicOscillator.cdf
𝑚𝜔 𝑖𝑃̂
𝑎̂ = √ 𝑥̂ +
2ħ √2𝑚𝜔ħ
𝑚𝜔 𝑖𝑃̂
𝑎̂+ = √ 𝑥̂ −
2ħ √2𝑚𝜔ħ
̂ 𝑎̂+ ] = 1
[𝑎,
Proof iii): The product gives ,
𝑚𝜔 𝑖𝑃̂ 𝑚𝜔 𝑖𝑃̂
̂ 𝑎̂+ ] = [(√
[𝑎, 𝑥̂ + ) , (√ 𝑥̂ − )]
2ħ √2𝑚𝜔ħ 2ħ √2𝑚𝜔ħ
𝑚𝜔 1
=√ × [(𝑥̂ + 𝑖𝑃̂), (𝑥̂ − 𝑖𝑃̂)]
2ħ √2𝑚𝜔ħ
1
= {[𝑥̂, 𝑥̂] + 𝑖[𝑃̂, 𝑥̂] − 𝑖[𝑥̂, 𝑃̂] + [𝑃̂, 𝑃̂]}
2ℏ
1
= {−2𝑖[𝑥̂, 𝑃̂]}
2ℏ
= {−2𝑖 × 𝑖ℏ}
=1
̂ can be represented by 𝒂
𝑖𝑣)Hamiltonian operator 𝑯 ̂ + as,
̂, 𝒂
̂ = ħ𝜔 (𝑎̂+ 𝑎̂ + 1 )..................................(3)
𝐻
2
1
̂ = ħ𝜔 (𝑎̂𝑎̂+ −
𝐻 )..................................(4)
2
𝑚𝜔 𝑖𝑃̂ 𝑚𝜔 𝑖𝑃̂
ħ𝜔 𝑎̂+ 𝑎̂ = ħ𝜔 (√ 𝑥̂ − ) (√ 𝑥̂ + )
2ħ √2𝑚𝜔ħ 2ħ √2𝑚𝜔ħ
1 𝑃̂ 2 𝑖𝜔
2 2
= 𝑚𝜔 𝑥̂ + + (𝑥̂𝑃̂ − 𝑃̂𝑥̂)
2 2𝑚 2
1
̂−
=𝐻 ħ𝜔
2
̂ = (𝑎̂+ 𝑎̂ + 1) ħ𝜔
∴𝐻
2
1
̂ = (𝑎̂𝑎̂+ −
Similarly, 𝐻 ) ħ𝜔
2
̂𝒂
a) Commutation relation of 𝒂, ̂,
̂ + with 𝑯
1
̂ , 𝑎̂] = ħ𝜔 [𝑎̂+ 𝑎̂ +
[𝐻 , 𝑎̂] = ħ𝜔 [ 𝑎̂+ 𝑎̂, 𝑎̂]
2
= − ħ𝜔 𝑎̂ ......................................... (6)
̂ , 𝑎̂+ ] = ħ𝜔𝑎̂+
Similarly [𝐻
̂ (𝑎̂𝑢) − 𝑎̂(𝐻
Or, 𝐻 ̂ 𝑢) = − ħ𝜔 𝑎̂𝑢
Repeated application of 𝑎̂ to the eigen functions will finally lead to the state of
lowest energy , the ground state say 𝑢0
i.e. 𝑎̂ 𝑢0 = 0
̂ 𝑢0 = ħ𝜔(𝑎̂+ 𝑎̂, + 1) 𝑢0 =
𝐻
1
ħ𝜔𝑢0 ............... (3)
2 2
1
This ground state energy 𝐸0 = ħ𝜔 is called zero point
2
energy.
̂ (𝑎̂+ 𝑢) − 𝑎̂+ (𝐻
𝑂𝑟, 𝐻 ̂ 𝑢) = ħ𝜔𝑎̂+ 𝑢
3
= ħ𝜔 𝑎̂+ 𝑢0
2
̂ (𝑎̂+ 𝑢0 ) = 𝐸1 (𝑎̂+ 𝑢0 )
𝐻
̂ (𝑢1 ) = 𝐸1 (𝑢1 )
𝐻
3
𝑢1 = 𝑎̂+ 𝑢0 is called 1st excited state. 𝐸1 = ħ𝜔 is the corresponding energy eigen
2
value.
Thus proceeding in the similar way, operating on 𝑢0 , n times with 𝑎̂+ , we obtain,
̂ {(𝑎̂+ )𝑛 𝑢0 } = (1 + 𝑛) ħ𝜔 (𝑎̂+ )𝑛 𝑢0
𝐻
2
̂ 𝑢𝑛 = 𝐸𝑛 𝑢𝑛
𝐻
1
Here 𝐸𝑛 = (𝑛 + ) ħ𝜔 is the energy eigen value in the nth state wave
2
function.
Q: Explain number operator
The operator 𝑎̂+ 𝑎̂ is called the number operator as, it determines the quantum
number of the eigen state.
𝐻̂ 1
𝑎̂+ 𝑎̂ = − ,
ℏ𝜔 2
We have,
1 1
𝑎̂+ 𝑎̂ = ̂ −
𝐻
ħ𝜔 2
Thus ,
1 1
𝑎̂+ 𝑎̂ 𝜓𝑛 (𝑥) = [ 𝐻̂ − ] 𝜓𝑛 (𝑥)
ħ𝜔 2
1 1
= [(𝑛 + ) − ] 𝜓𝑛 (𝑥)
2 2
= 𝑛 𝜓𝑛 (𝑥)
Q: Using the raising and lowering operator Eigen function for 1d LHO.
Let 𝑢0 be the ground state wave function. Hence the lowering operator,
𝑚𝜔 𝑖𝑃̂
𝑎̂ = √ 𝑥̂ +
2ħ √2𝑚𝜔ħ
operating on it will produce zero.; i.e.
𝑎̂𝑢0 = 0
𝑚𝜔 𝑖 𝜕
𝑜𝑟, (√ 𝑥̂ + (−𝑖 ℏ )) 𝑢0 = 0
ħ √𝑚𝜔ħ 𝜕𝑥
𝑚𝜔 ℏ 𝜕
𝑜𝑟, (√ 𝑥̂ + ( )) 𝑢0 = 0
ħ √𝑚𝜔ħ 𝜕𝑥
𝑚𝜔
𝐿𝑒𝑡 𝑦 = √ 𝑥 ;
ħ
ℎ𝑒𝑛𝑐𝑒,
𝑑
(𝑦 +)𝑢 = 0
𝑑𝑦 0
𝑑𝑢0
𝑜𝑟, 𝑦𝑢0 + =0
𝑑𝑦
𝑑𝑢0
= −𝑦𝑑𝑦
𝑢0
Solving,
𝑦2
−
𝑢0 = 𝐴0 𝑒 2
𝑚𝜔 2
= 𝐴0 𝑒 − 2ħ
𝑥
𝐴0 2 1 1
𝑜𝑟, . 2. . Γ ( ) = 1
𝛼 2 2
𝛼
ℎ𝑒𝑛𝑐𝑒, 𝐴0 2 =
√𝜋
Therefore normalize wave function,
1 1
𝛼 2 𝑦2 𝑚𝜔 4 − 𝑚𝜔𝑥 2
−
𝑢0 = ( ) 𝑒 2 =( ) 𝑒 2ħ
√𝜋 𝜋ℏ
Eigen value
𝐸𝑛 = [𝑛𝑥 + 𝑛𝑦 + 1]ℏ𝜔 = [𝑛 + 1]ℏ𝜔
It is degenerate, and few energy levels are given below.
ground 0 0 0 ℏ𝜔 𝜓0,0 1
Ext. 0 1 𝜓0,1 2
2nd 2 0 𝜓2,0
1 1 𝜓1,0
0 3 𝜓3,0
Ext 1 2 𝜓2,1
2 1 𝜓1,2
ground 0 0 0 0 3 𝜓0,0,0 1
ℏ𝜔
2
1st 1 1 0 0 𝜓1,0,0
Ext. 0 1 0 5 𝜓0,1,0 2
𝜔
2
0 0 1 𝜓0,0,1
2 0 0 𝜓2,0,0
0 2 0 𝜓0,2,0
2nd 2 0 0 2 7 𝜓0,0,2
ℏ𝜔
2
Ext 1 1 0 𝜓1,1,0 6
1 0 1 𝜓1,0,1
0 1 1 𝜓0,1,1
(𝒏+𝟏)(𝒏+𝟐)
NOTE: In general degeneracy for the nth energy level is 𝟐
Q: A pendulum with 1 gm bob has massless string 250mm. long the period of
pendulum is 1 sec. What is its zero-point energy? Would you think that zero-point
oscillations be detectable? What is the separation between the adjacent energy
levels? 2+1+1
(2011)
3
Explain i) parity, < 𝑥 >, < 𝑝𝑥 > for this state. (2016)
i) 𝜓(𝑥, 𝑡)
ii) <E>
iii) < 𝑃𝑥 >, < 𝑥 > as function of time. (2015)
Q: Calculate [𝑃̂, 𝐻
̂ ] for 1D lho hence comment on the result.
Q: Show that the probability density of 1d LHO at the classical limit of motion to the
1
probability density at the origin is for the ground state.
𝑒
Q:For a 1d LHO in its ground state show that probability of finding it beyond the
classical limit is 0.16
1
Q: i) A 1D harmonic oscillator moves in a potential 𝑉(𝑥) = 𝑘𝑥 2 + 𝐶𝑥. C is a
2
constant. Find the energy eigen value.
1
iii) An electron confined in a potential well 𝑘𝑥 2 is subjected to an electric
2
field 𝜖. Find the shift in the energy level.
1
Q: An electron confined in a potential well 𝑉(𝑥) = 𝑘𝑥 2 , Where K is a constant, is
2
subjected to an electric field e along the x-axis. Find the energy eigen value. CU
(2017,03)