[양자역학] 2.2-(1) Example: 파동 패킷 Wave Packet

Free particle의 plane wave ψ(x)=eiℏpxψ(x)=eiℏpx는 Hamiltonian operator의 eigenvector이지만, probability density fucntion이

ψ∗(x)ψ(x)=e−iℏpxeiℏpx=1ψ∗(x)ψ(x)=e−iℏpxeiℏpx=1

즉, 모든 공간에 걸쳐 확률이 균일하게 퍼져있기 때문에, 일정한 공간에서 찾을 확률이 높은 입자의 특성을 표현하기에는 적합하지 않다. 따라서 입자의 위치를 표현하기 위해서는 plane wave를 중첩하여 확률이 작은 영역에 몰려있는 wave function을 사용하는데 이 wave function을 wave packet이라고 한다. 이번 페이지에서는 wave packet의 한 예를 살펴보고 몇가지 특성을 살펴본다.

 

#Gaussian Wave Packet

wave function이 다음과 같이 Gaussian distribution 형태라고 가정하자.

ψ(x)=(12πϵ)14e−x24ϵψ(x)=(12πϵ)14e−x24ϵ

일단, 이 wave function이 normalized 되어있음을 확인하자.

∫∞−∞ψ∗(x)ψ(x) dx=∫∞−∞(12πϵ)12e−x22ϵ dx=1∫−∞∞ψ∗(x)ψ(x) dx=∫−∞∞(12πϵ)12e−x22ϵ dx=1

 

이 wave function의 probability density function ψ∗(x)ψ(x)=(12πϵ)12e−x22ϵψ∗(x)ψ(x)=(12πϵ)12e−x22ϵ의 그래프는 아래 그림과 같다.

 

image by Wolfram Mathematica

 

위의 그래프에서 볼 수 있듯이 ϵϵ이 작을수록 x=0x=0 중심으로 발견할 확률이 모이고, ϵϵ이 커질수록 확률이 흩어진다는 것을 확인할 있다. 이제 이 wave packet의 momentum space wave function을 구해보자.

^ψ(p)=∫∞−∞ψ(x)1√2πℏe−iℏpx dx=∫∞−∞(18π3ϵℏ2)14e−(x24ϵ+iℏpx) dxψ^(p)=∫−∞∞ψ(x)12πℏe−iℏpx dx=∫−∞∞(18π3ϵℏ2)14e−(x24ϵ+iℏpx) dx

이를 풀기 위해, exponential의 지수를 정리하면,

x24ϵ+iℏpx=(x+2iϵℏp)24ϵ+ϵℏ2p2x24ϵ+iℏpx=(x+2iϵℏp)24ϵ+ϵℏ2p2

따라서

^ψ(p)=(18π3ϵℏ2)14e−ϵℏ2p2∫∞−∞e−(x+2iϵℏp)2/4ϵ dxψ^(p)=(18π3ϵℏ2)14e−ϵℏ2p2∫−∞∞e−(x+2iϵℏp)2/4ϵ dx

이제 Cauchy integral theorem을 이용하면(복소함수 적분에 익숙하지 않다면 그냥 Gaussian distribution 적분으로 생각해도 된다)

^ψ(p)=(2ϵπℏ2)14e−ϵℏ2p2ψ^(p)=(2ϵπℏ2)14e−ϵℏ2p2

Momentum space wave function에서 probability density function의 그래프는 다음과 같다.

image by Wolfram Mathematica

 

momentum space와 position space의 그래프를 비교해보면, ϵϵ이 작아질수록 position space의 확률은 집중되는데 반하여 momentum space의 확률은 퍼지게 된다. 반대로 ϵϵ이 커질수록 position space의 확률은 퍼지는데 반해 momentum space의 확률은 집중된다. 이러한 경향의 극단적인 형태가 plane wave이다. plane wave의 position space의 확률은 전 영역으로 균일하게 퍼져있는데 반해 momentum space의 확률은 딱 한 개의 값으로 집중되어 있는 Dirac delta function이 된다. 반대로 position space의 확률이 Dirac delta function 형태라면, momentum space의 확률은 inverse Fourier transform으로부터 전 영역에 균일하게 퍼지게 됨을 알 수 있다.

 

#Uncertainty Principle

이러한 경향은 Heisenberg uncertainty principle(불확정성 원리)으로 이해할 수 있다. 위치를 정확히 하려고 wave packet을 집중시키면 운동량의 wave packet이 넓어지고, 반대로 운동량의 wave packet을 집중시키면 위치의 wave packet이 넓어진다고 해석할 수 있다. 그렇다면 Gaussian wave packet은 얼마나 좋은 wave packet일까? 만약 위치를 정확히 하려고 wave packet을 집중시켰는데, 운동량이 분산되는 것이 너무 심하다면 좋지 않은 wave packet이라고 할 수 있을 것이다.

 

XX의 표준편차를 구하기 위해, XX의 기대값과 X2X2의 기대값을 구하면,

E[X]=∫∞−∞ψ∗(x)xψ(x) dx=∫∞−∞x(12πϵ)12e−x22ϵ dx=0E[X]=∫−∞∞ψ∗(x)xψ(x) dx=∫−∞∞x(12πϵ)12e−x22ϵ dx=0

E[X2]=∫∞−∞ψ∗(x)x2ψ(x) dx=∫∞−∞x2(12πϵ)12e−x22ϵ dx=ϵE[X2]=∫−∞∞ψ∗(x)x2ψ(x) dx=∫−∞∞x2(12πϵ)12e−x22ϵ dx=ϵ

똑같이 PP와 P2P2의 기대값을 momentum space에서 구하면,

E[P]=∫∞−∞^ψ∗(p)p^ψ(p) dp=∫∞−∞p(2ϵπℏ2)12e−2ϵℏ2p2 dp=0E[P]=∫−∞∞ψ^∗(p)pψ^(p) dp=∫−∞∞p(2ϵπℏ2)12e−2ϵℏ2p2 dp=0

E[P2]=∫∞−∞^ψ∗(p)p2^ψ(p) dp=∫∞−∞p2(2ϵπℏ2)12e−2ϵℏ2p2 dp=ℏ24ϵE[P2]=∫−∞∞ψ^∗(p)p2ψ^(p) dp=∫−∞∞p2(2ϵπℏ2)12e−2ϵℏ2p2 dp=ℏ24ϵ

이들로부터 표준편차를 구하면,

(ΔX)=√E[X2]−(E[X])2=√ϵ(ΔX)=E[X2]−(E[X])2=ϵ

(ΔP)=√E[P2]−(E[P])2=ℏ2√ϵ(ΔP)=E[P2]−(E[P])2=ℏ2ϵ

따라서 표준편차들의 곱은

(ΔX)(ΔP)=ℏ2(ΔX)(ΔP)=ℏ2

Gaussian wave packet은 uncertainty principle에서 표준편차들의 최소값에 해당한다는 것을 알 수 있다.1 즉, Gaussian wave packet이 입자의 위치와 운동량을 한정하는 한계라고도 할 수 있다.

 

#Dispersion

Momentum space에서 free particle Hamiltonian operator는 단순히 p22mp22m을 곱해주는 역할이므로 wave packet이 시간에 대하여 어떻게 변하는지 확인하는데 momentum space가 더 편리하다. Time dependent Schrodinger equation

iℏddt^ψ(p,t)=Hψ(p,t)=p22m^ψ(p,t)iℏddtψ^(p,t)=Hψ(p,t)=p22mψ^(p,t)

로부터

^ψ(p,t)=^ψ(p,0)e−iℏp22mtψ^(p,t)=ψ^(p,0)e−iℏp22mt

를 얻는다. 만약 t=0t=0일 때, wave function이 Gaussian wave packet일 때,

^ψ(p,t)=(2ϵπℏ2)14e−(ϵℏ2+it2mℏ)p2ψ^(p,t)=(2ϵπℏ2)14e−(ϵℏ2+it2mℏ)p2

이 wave function의 probability density function은

^ψ∗(p,t)^ψ(p,t)=(2ϵπℏ2)12e−2ϵℏ2p2ψ^∗(p,t)ψ^(p,t)=(2ϵπℏ2)12e−2ϵℏ2p2

로 momentum space에서 probability density function은 변하지 않는다. 이제 이 함수를 다시 inverse Fourier transform을 통해 position space로 바꿔보자.

ψ(x,t)=∫∞−∞^ψ(p,t)1√2πℏeiℏpx dp=∫∞−∞(2ϵπℏ2)141√2πℏe−(ϵℏ2+it2mℏ)p2+iℏpx dpψ(x,t)=∫−∞∞ψ^(p,t)12πℏeiℏpx dp=∫−∞∞(2ϵπℏ2)1412πℏe−(ϵℏ2+it2mℏ)p2+iℏpx dp

지수를 정리하면,

(ϵℏ2+it2mℏ)p2−iℏpx=[(2mϵ+iℏt2mℏ2)12p−i(m4mϵ+2iℏt)12x]2+m4mϵ+2iℏtx2(ϵℏ2+it2mℏ)p2−iℏpx=[(2mϵ+iℏt2mℏ2)12p−i(m4mϵ+2iℏt)12x]2+m4mϵ+2iℏtx2

위에서 한 것과 같은 방식으로 적분하면,

ψ(x,t)=(m2mϵ+iℏt)12(2ϵπ)14e−m4mϵ+2iℏtx2ψ(x,t)=(m2mϵ+iℏt)12(2ϵπ)14e−m4mϵ+2iℏtx2

t=0t=0일 때 처음에 본 Gaussian distribution이 된다는 것을 확인할 수 있다. 이제 probability density function을 구하면,

ψ∗(x,t)ψ(x,t)=√m24m2ϵ2+ℏ2t2√2ϵπe−8m2ϵ16m2ϵ2+4ℏ2t2x2ψ∗(x,t)ψ(x,t)=m24m2ϵ2+ℏ2t22ϵπe−8m2ϵ16m2ϵ2+4ℏ2t2x2

Momentum space와는 다르게, position space에서 probability density function은 시간에 대하여 변한다는 것을 볼 수 있다. Probability density function의 그래프는 다음과 같다.

 

image by Wolfram Mathematica

 

그래프에서 확인할 수 있는 것처럼, 처음에는 x=0x=0에 집중되어 있는 확률이 시간이 지나면서 점점 퍼져가는 것을 볼 수 있다. 만약 t=0t=0에서 position space wave function이 Gaussian wave packet이 아니라 Dirac delta function처럼 정확히 한 곳에 집중되어 있다고 하더라도 마찬가지로 시간이 지나면서 확률이 점점 퍼져간다. 이렇게 시간이 시간이 지나면서 점점 퍼져나가는 wave packet의 특성을 dispersion(분산)이라고 한다.2

 

#Moving Wave Packet

만약 wave packet이 평균 운동량 p0p0의 값으로 움직인다고 한다면, 간단히 momentum space wave function을

^ψ(p)=(2ϵπℏ2)14e−ϵℏ2(p−p0)2ψ^(p)=(2ϵπℏ2)14e−ϵℏ2(p−p0)2

라고 하면 된다. 이를 inverse Fourier transform을 하면,

ψ(x)=∫∞−∞^ψ(p)1√2πℏeiℏpx dp=∫∞−∞(2ϵπℏ2)14e−ϵℏ2(p−p0)21√2πℏeiℏ(p−p0)xeiℏp0x d(p−p0)ψ(x)=∫−∞∞ψ^(p)12πℏeiℏpx dp=∫−∞∞(2ϵπℏ2)14e−ϵℏ2(p−p0)212πℏeiℏ(p−p0)xeiℏp0x d(p−p0)

이므로

ψ(x)=(12πϵ)14e−(x24ϵ−iℏp0x)ψ(x)=(12πϵ)14e−(x24ϵ−iℏp0x)

를 구할 수 있다.

 

 

1. 1.4-(3) Example: 하이젠베르크의 불확정성 원리 Heisenberg's Uncertainty Principle 참고. [본문으로]
2. Dispersion에 대해서는 --EM,group velocity-- 참고. [본문으로]

 

 

 

 

https://archive.ph/Ywrdv

https://elementary-physics.tistory.com/107

[양자역학] 2.2-(1) Example: 파동 패킷 Wave Packet