第二章原子构与性质§21氢原子和类氢原子的薛定谔方程及其

第二章原子构与性质§21氢原子和类氢原子的薛定谔方程及其
第二章原子构与性质§21氢原子和类氢原子的薛定谔方程及其

第二章 原子结构与性质

§2.1.氢原子和类氢原子的薛定谔方程及其解 2.1.1.单电子原子的薛定谔方程

H 原子和He +、Li 2+ 等类氢离子是单原子,它们的核电荷数为Z ,若把原子的质量中心放在坐标原点上,绕核运动的电子离核的距离为r ,电子的电荷为-e ,其静电作用势能为:

r Ze V 0

2

4πε-=

将势能代入薛定谔方程:

得 0)(2

2282

=ψ+

+

ψ?r

Ze h m

E π

或ψ=ψ-

?-

E r

Ze m

h ][2

2

2

2

为了解题方便,将x 、y 、z 变量换成极坐标变量r 、θ、φ。 其关系:φθcos sin r x = φθsin sin r y =

φcos r z =

2222

z y x r

++=

2

1

)

/(cos 2

22z y x Z ++=θ

x y tg /=φ

})(sin )({2

222

sin 1sin 121

2

φθθθ

θθ

??????????++=

?

r r

r r 代入薛定谔方程:

)()(sin )(2

2

22

222228sin 1

1sin 1121=ψ+

+

++???ψ???????

r

Ze h m

r r r r

r E r πφθθθ

θθ

2.1.2.分离变量§法:

上述的方程是含三个度量的偏微分方程,要解这个方程可用度数分离法将其化为三个分别只含一个度量的常微分方程求解。

含:)()()(),,(φθθΦΘ=Φψr R r 代入方程:并乘以ΘΦR r θ

22sin 移项

得:

)

(sin )(sin )(228sin 2sin 12222

2V E r r h

u d d d d

dr dR dr

d

R d d ----=

ΘΘΦΦθθ

πθθ

θθφ左边不含r 、θ,右边不含φ,欲左右两边相等必等于同一个常数(-m 2 )

Φ-=Φ

222m d d φ

, 而右边可为:(除以sin θ)

)(sin )()(sin

1sin 821

22

2

θ

θθ

πθ

d d d d m h

ur dr dR dr

d

R V E r ΘΘ-=

-+ 则有:

K d d d d m =-ΘΘ)(sin sin

1sin 22

θ

θ

θ

θθ

K E r r

Ze h

ur dr dR dr

d

R =++)()(2

2

2

2821

π

2.1.

3.方程解的结果 2.1.3.1.Φ(φ)方程的解

022

2=Φ+Φ

m d d φ

这是一个常系数二阶齐次线性方程,有两个复函数的独立解。 |)|(]exp[m m im A m

±==Φφ

Φ符合波函数品优条件:连续、单值、电子边界条件(归一) 1]exp[]exp[202

20

*=-A =ΦΦ??

φφφφπ

πd im im d m

m π21=A ]exp[][21φπ

im m

α、φ周期变化,Φm 值不变 )2()(πφ

φ+Φ=Φm m

]2exp[]exp[)]2(exp[]exp[πφπφφim im im im =+= 得: 1]2exp[=πim

根据Euler 公式 φφφm i m im sin cos ]exp[+= 1)2sin()2cos(=+ππm i m 故m 取值必须为:)(,2,1,0量子化ΛΛ±±=m

2.1.

3.2.○H (θ)方程的解。

Θ=-ΘΘ两边乘k d d d d m )(sin sin

1sin 2

θ

θ

θ

θ

0)(sin 22

sin sin 1

=Θ+Θ-Θk m d d d d

θ

θθ

θθ

要使方程得到收敛解,并有确定值,k 必须限制(解过程很复杂)。 ΛΛ,2,1,0)1(=+=l l l k

l d d m l m l l l m

l m l m l m l )1(cos )cos 1(}{2cos 2|)!|(2|)!|)(12(!

21,|

||

|2

|

|2

1-??

-=

Θ+++-+θθθ

(缔合勒让德函数) 2.1.3.3.R (r )方程的解。 同理对R (r )方程,将)1(+=l l k

代入

0])([)(2

02

2

2)1(4221=-+

++R E r r l l r Ze m

dr dR dr d

r

πεη

解得:

)]

([}){(1

22

2130

])!([(2)!1(32,l n d d d d l l n n l n na Z l n e e

e R l

n l n l n l +-+--++++-?

-=

ρρ

ρρρρρ

其中:0

2na Zr

=

ρ,并有:Λ

3,2,12

22422=-=

n E h

n Z me n π

最后得到完整的单电子波函数:

φ

θφθφθψim m l i l l n i a Zr i m m l l n m l n e P e c c r R r na Zr

)(cos })({)

()()(),,(||11

,,,,00

--+-=∑

=Φ?Θ?=

),()(,,,,φθψm l l n m l n Y r R =

上述各波函数都归一化,即

??

?

?===ΘΘ=ΦΦ∞π

ππ

φθθθθφ0*2

*20

0**1sin ,11

sin ,

1d d Y Y dr Rr R d d

φθθτφθθππ

d drd r d d drd d r sin ,1sin 2220

*00

==ψψ?

??∞

并有量子数n 、l 、m 。

2.1.4.结论

1量子数的物理意义。 量子数以整数跳跃取值,不连续,量子化的体现。 A 、主量子数 l n n >=Λ

3,2,1

主量子数决定了电子状态的能量,基态时n=1。

Λ

3,2,1)(595.13)

(595.1310178.22

2

22

421812=?

-=-=?-=-=

-n eV E eV J E E n n h

n Z me n π

B 、角量子数 m l n l

≥-=1,3,2,1,0Λ

将角动量平方算符2

∧M 作用于单电子波函数,可得:

ψ+=ψM ∧2

22

))(1(πh l l

1,2,1,0)1(||))(1(2222-=?+=+=M n l l l M l l h

h Λπ

π或

由此可见,上式中量子数l 决定了电子的角动量大小,故称角量子数。原子的角动量和原子的磁矩有关,只要有角动量也就有磁矩。其关系为:

e c

m e

z c

m e

e

e

M -=-=22μμ

βμππ)1()()1(||422+=+=+=l l l l l l l c

m eh

h c

m e e

e

βμπ

m m h

c m e z

e

-==-22

其中,21

41027.9-?==c

m eh e

πβ

J /K 称波尔磁子。 C 、磁量子数 m

角动量度Z 方向的分量与m 有关,即: l m m m h z h z

±±±==M ψ=ψM ∧

,,2,1,022Λπ

π

这也说明,Z 方向是磁场方向,m 决定了角动量在这方向上的分量大小,是量子化 的,称磁量子数。

电磁学指出,磁铁在磁场中取向不同时,能量会不同,其表示: H =H -=H -=E βθ

μm u z cos

由此可见,原子像小磁铁在无外场时,n ,l 相同,m 不同时,能量本来相同,但当处于外加磁场中,m 不同时,能量就不同了,称在磁场中能级分裂,如: 1=l

此现象称为塞曼(Zelman )效应,证实了角动量在磁场方向分量的量子化。

2),()()()()(),,(,φθφθφθY r R r R r m l =ΦΘ=ψ 称为单电子波函

数,也称为原子轨道或原子轨迹。 ○

3量子数取值是有限制的。 )

(2,1,0:)1(,3,2,1,0:,3,2,1:l m n l n ±±±-ΛΛ

Λ

同一个n 之下有

∑-===+1

2

)12(n l l n l 个不同波函数,则说明同一个En 之下有n 个独立状态,称简并度(能量相同,运动方式不同)。

2.1.5.角函数(球谐函数)

2.1.5.1.几种常用的球谐函数Y (θ,φ) ),(Y ,φθm l ——球谐函数

φ

π

φ

π

π

φ

π

φ

π

π

πθθθθθθθ

i i i i e e e e ±-±-±-?=

Y ?=

Y -=

Y ?=Y ?=

Y =Y =Y 23215221651221652083118311431041

00sin cos sin )

1cos 3(sin sin cos

2.1.5.2.原子轨道角函数

对同一个能量下的简并态,n ψψψΛ,,21

221

1

Eψ=ψH Eψ=ψH ∧

如这些波函数组成新的线性波函数,仍然是该能量的本征函数即:是一个新的简并态。 )()(22112211ψ+ψE =ψ+ψH ∧

c c c c 应用光谱学的习惯表示,将 Λ3,2,1,0=l

等状态化为:Λ

f d p s ,,,,而在

每一个l 之下有12+l 个独立状态用: 222,,,,,,,y x yz xz xy z z y x -等加

以区别。

φ

θφ

θφ

θφθθθφ

θφ

θθπ

π

π

ππ

π

π

π

φφπ

π

2sin sin )(2cos sin )(sin 2sin )(cos 2sin )()

1cos 3(cos cos sin )()(cos sin )()(sin )()(2161522222

12161522222

12241512212141512212

1

216520243102232111112122321

2232111112

14100=

Y -Y =

=

Y +Y =-=Y -Y ==

Y +Y =

-=Y ==Y =P =Y -Y =P =+=

Y +Y =

P =Y =-i i z i i y i i x dxy y dx dyz dxz dz e e S

2.1.6.波函数和电子云的图形

波函数ψ是三维空间坐标的函数,可用图形表示出来,从而使抽象的数学表达式成为具体的图形。 2.1.6.1.角度分布图

○a Y 角度分布图(Y 的图形) )(41

00常数π

=Y =Y s

球面(任意方向都一样) π

41=r θπ

cos 4310=

Y =Y pz

与φ无关

XY 平面几率为零)0,0()()90(===Y X 节面οθ

在Z 轴正负方向上几率最大)180,0(οο=θ

两个在XY 平面上下相切的球

)

(180,90,0max

)({cos 0)

1cos 3(''

24454161253

12161520极大值节面οοοοο===?±==Y

-=

Y =Y θθ

θθπ

Y dz

同理: 0Y 0,X ,YZ :==平面为节面x P

0Y 0,y ,XZ :==平面为节面y

P

00

{0{0:====?=Y =平面平面YZ ZX x y xy dxy

00

0{0{0:====?=Y =平面平面YZ XY x z xz dxz 000

{0{0:====?=Y =平面平面XZ XY y z yz dyz

y x y x y x y x y x y dx

-==?=-+?=--{0))((0:2222

Y 正负交替排布,交界的地方为节面,与n 无关。 ○

b Y 2 的角度分布图 Y 2

(电子云),表示电子出现在径向θ、φ方向上单位立体角的几率,还表示同一球面上各点几率密度|ψ|2 的相对大小。 电子云的角度分布与原子轨道角度分布是相似的。

区别:Ⅰ、电子云的角度分布比原子轨道角度分布要“瘦”一点;

Ⅱ、原子轨道分布有正、负之分,而电子云均为正值。

注:角度分布正负只与Y有关,整个ψ的正负还要考虑R的符号。化学反应与电子运动有关,因此波函数的性质尤其是Y的图形和正负号有重要作用。 (2)径向分布函数和径向分布图

对于ψ=RY,其中,R~r 的函数在球坐标表达式中,空间某点(r ,θ,φ)附近的体积元在其中找到的电子的几率为:

2

2222

200

2

2

2

2202

02)()(sin ||sin ||||R r r D dr

r D dr R r d d dr R r d drd r R d ===Y =Y =ψ?????φ

θθφ

θθτππππ

)(D r 称为径向分布函数,其物理意义是:dr r )(D 代表在半径为r 和

dr r +的两个球壳夹层内找到的电子几率,电子云随半径变化情况。

)(D r ~r 作用称为径向分布图,有如下结论: ○1n 决定电子分布离核远近。

○2核上无电子出现,但n 大的电子也有在核近的地方活动。 ○3有)(l n -峰;)1(--l n 个节面。

b. 2Y 的角度分布图

2Y

(电子云),表示电子出现在径向θ、φ方向上单位立体角内的几率,还表示同一球面上各点几率密度2||ψ的相对大小。

电子云的角度2Y 分布与原子轨道角度Y 分布是相似的。

区别:

1、电子云的角度分布比原子轨道角度分布要“瘦”一点。

2、原子轨道分布有正负之分,而电子云均为正值。

注:角度分布正负只与Y 有关,整个Ψ的正负还要考虑它的符号,化学反应与电子运动有关,因此波函数的性质,尤其是Y 的图形和正负号有重要作用。 2.1.6.2.径向分布函数和径向分布图。

对于的函数其中r R RY

~=ψ 在球坐标表达式中,空间某点(r, θ,φ)

附近的体积元φθθτd drd r d sin 2= 在其中找到的电子的几率为:

dr

r D R r r D dr R r d d Y dr R r d drd r Y R d )()(sin ||sin ||||2

22220

20

2

2

222

202

=====ψ?

??

?

?令φθθφ

θθτπ

ππ

π

D(r)称为径向分布函数,其物理意义是:dr r D )(代表在半径为r 和r+dr 的

两个球壳夹层内找到的电子几率,电子元随半径变化情况。

D(r)~r 作用称为径向分布图(P56)。有如下结论: ○

1n 决定电子分布离核远近。 ○

2核上无电子出现,但n 大的电子也有在核近的地方活动。 ○3有(n-l )峰;(n-l-1)个节面。

2.1.6.

3.等密度图和界面图

1等密度图:通常把在纸面上不易表达的三维数值,通过原子核及某些坐标轴的截面上,把面上各点的r, θ,φ值代入Ψ中,然后要据Ψ值的正负和大小画出等密(值)线,称为等密度图。

2界面图:为了了解电子分布的几率,可取一个等密度面,使在面内出现的几率达到总几率的一定百分数,如:

τ

τd d ??????ψ=ψ=||5.0||9.0

注:电子云是指电子长时间出现的结果(统计结果),并非无限多个电子构成“云”。 §2.2.多电子原子的结构 2.2.1.氦原子 2.2.1.1.原子单位

把质量用电子质量e m ,电荷为质子电荷e ,长度为?a ,能量为电子伏特(1哈密顿)称为原子单位,则波动方程可简化为:

Eψ=ψ-ψ?-Z

r 221

2.2.1.2.轨道近似法

氦原子有两个电子,在原子单位下,波动方程为: Eψ=ψ

+--

?+?-Z Z ])([122

1

1

2

22

121r r r

),,,(222;1,11φθφθr r ψ=ψ

近似法,忽略两电子的排斥势

Eψ=ψ-?-+ψ-?-)()(2

1

2

2

22122121r r

分离变量:)2()1()2,1(21ψψ=ψ 简化为两个氢原子方程:

11122121)(1

ψE =ψ-?-r 22222

221)(2

ψE =ψ-?-r

21E +E =E 波函数223

1

23

21

221

1)2()2()2()1(r r e e --=

ψ=ψπ

π

类氢原子:ev ev n n 4.546.1346.13212

2

-=?-=E =E ?-=E Z

)0.79:(8.108ev ev --=E 实验值 偏差较大。

2.2.2.多电子原子结构

多电子原子与单电子原子最主要的区别是含有两个以上的电子,这样在电子间就存在静电库仑作用,其排斥能为:

ij

ij ij r e r e E 22

4==?πε

对含N 个电子的原子,其哈密顿算符应为:

∑∑∑∑≠==Λ

+Z -?-

=H i i j ij

N

i i i N

i r e r e m

h 2122

1

22

218π

第一项为动能,第二项为单电子位能,第三项是电子间排斥能;由于第三项涉及到两个电子间的坐标:222)()()(j i j i j i ij z z y y x x r -+-+-= 不可能分离,求解困难,要采用近似法处理。

2.2.2.1.中心力场近似

这种中心力场模型主要认为除考虑的电子外其它电子所产生的一个有效平均场是一种球对称场。这有效平均场是指每一个电子在这种场中其它电子坐标都在对电子排斥能求平均的过程中被去除掉了,唯独只剩下各该电子自己坐标作为变量,能量算符在形式上变成和其它电子的相对位置无关。

由于是球形对称,其排斥能可写)(i i r u 只与i r 有关,与θ、φ无关。在这种

近似下,其中某个单电子薛定谔方程:

)()()](48[)(22

22i i i i i i i i i i i r r r u r ze m h r ψE =ψ+-?-=ψH Λ

οπεπ

那么整个体系的波函数:(有n 个电子)

n

n E ++E +E +E =E ψψ?ψ?ψ=ψΛΛ321321

总结起来,在这种近似下,多电子原子可以看作是Z 个电子在各自的原子轨道上运动,这些轨道像氢原子轨道一样,用n 、l 、m 表示,角度部分与氢原子轨道相同,只是r 部分不同,能量不仅与n 、l 有关,而且与其它电子有关。整个体系能量为各个电子能量之和,整个波函数是所有电子的波函数的重组。

2.2.2.2.屏蔽效应(Slater 规则)

在中心力场近似下,我们可以看出:电子的排斥能一定会抵消核吸引能的作

用。因此可假设: i

i ij i i r e r e r u 2

2)(σ== i σ—电子i 的屏蔽常数,相当于抵消了i σ个原子核正电荷的作用。这样i 电子

就好像处在一个以原子核为中心的单中心平均有效场中,其位能形式:

i

i i i i i i r e Z r e Z r e r Ze r V 2

*222)()(-=--=+-=σσ

i Z Z σ-=*称为有效核电荷,其哈密顿算符与单电子原子算符相似:

i

i i r e Z m

h 2

*2

228-?-

=H Λ

π

其能量: )(6

.13)(6

.132

2*2

2

ev n

Z

n

Z i i -=--=E σ

1屏蔽常数的影响因素 ○

a 外层电子对内层电子屏蔽作用小(有一定屏蔽作用)。 ○

b 内层电子对外层电子屏蔽作用大。(0.85~1.0) ○

c 同层电子l 大的屏蔽作用小。(0.2~0.45)

考虑电子i 本身 ↑↑↑↑i i l n σσ 考虑其余电子时 ↓↑

↑i i l n σσ

2多电子原子能级计算: )()(6

.132

2

ev n Z i i σ--=E

a n 相同,l 不同时,*

3*333p d p

d l E E E ↑

○b n , l 都不同时,有效主量子数 n’ 不同

)

()

()(4.0')(7.0'4343离子原子离子原子s d s d l n n l n n E E +=+=

这表明考虑屏蔽效应后,电子能量与 n, l, 有关,而且与其它电子数目和状态有关。

例如:6

.131

85

.035

.01

70.53.0613.01

)

(1

,01*-=E ===>=-=E ==++==s n c b a n s a n cz by ax l 碳原子电子光谱数据σ x — 同n , l 的电子数 y — (n-1)的电子数 z=(n-2)或更

内层的电子数

Co :27626224333221s d p s p s s

对于4s 电子,1.231011585.0135.0=?+?+?=σ )(93.124

)1.2327(6.132

2

4ev s -=--=E

K :162622433221s p s p s s

对于4s 电子,8.16101885.0=?+?=σ

)

(81.47

.3')

(11.44

)8.1619(6

.1342

2

4ev n ev s s

-=E =-=--=E

2.2.2.

3.钻穿效应

我们已经知道,对于n 相同 l 不同的轨道能量是:l 越大 ,E 越高。

nf nd np ns <<<

这就是钻穿效应,也就是屏蔽多少的问题,屏蔽越多,能量越高,反之越低。一般讲,在原子核附近出现几率较大的电子屏蔽较少,直接接受核电荷的吸引,能量较低;相反,在原子核附近出现几率较少,屏蔽多,能量则高。

根据量子力学求平均值的原理:多电子原子的位能平均值 V E E

V 2

12==

???∞Λ

=Z -=Z -=

E ψZ -ψ=E Z -

=022

*2

2

**2**2

*))(()(21)(21)(21)(r D Rr Ddr r

e dr Rr r e R d r

e r

e r V τ

由此可见径向分布函数D 对原子轨道的能量有重要影响。在原子核附近,r 值小、Z*值大、V 值低。对于n 相同l 不同的轨道在核附近出现的几率大小可以比较它们的径向分布函数的大小。l 越小,离核越近。 其在核附近电子出现的几率相对大小为:Λnf nd np ns >>> 而轨道能量次序则相反为:Λnf nd np ns <<<

这种由于l 不同轨道电子云径向分布不同、电子云钻到核附近的几率不同,在核附近电子出现几率较大的轨道称钻得较深的轨道,引起能量不同的现象称为钻穿效应。也反映了外层电子对内层电子也有一定的屏蔽作用。因为电子活动不局限于主峰上,有一部分钻到离核很近的内层。

2.2.2.4.电子的核外排布规则

原子处在基态,其核外电子排布遵守三个规则:

1、Pauli 不相容原理:在一个原子中,没有两个电子有完全相同的四个量子数,即一个原子轨道最多只能排两个电子,而且这两个电子自旋方向必须相反。

2、能量最低原理:在不违背Pauli 原理的条件下,电子优先占据能级较低的原子轨道,使得整个原子体系能量处于最低。

3、Hund 规则:在能级高低都相等的轨道上,自旋平行的电子数目最多时,原子的能量最低。因此,全充满,半充满或全空状态较稳定,此时电子云分布近于球形。

原子轨道顺序:1s 2s 2p 3s 3p 4s 3d 4p 5s 4d 5p 6s 4f 5d 6p ……次序“颠倒”的主要原因是钻空效应。原子电子的核外排布称组态。

§2.3.原子状态和原子光谱项

对于多电子原子,可以近似地认为原子中的电子在各自的轨道上运动,其运动状态可由量子数n, l, m, m s 规定。再根据电子核外排布(称为组态)说明了原子中各个单电子所处的状态。但这不足以表达原子整体运动状态。因为各电子间相互作用,不是简单的加和,而且电子轨道运动和自旋运动产生的磁矩也要相互作用,情况也是复杂的。无法与光谱实验联系。

实验表明光谱实验可以表达原子的整体状态,原子状态变化了,可从原子光谱中反映出来。所以对一个原子状态,可由原子的轨道角动量、自旋角动量和总角动量规定,这样就有另一套量子数 n 、L 、S 、J 、M L 、M S 、M J 与光谱实验相联系来表达原子的状态。

2.3.1.角动量偶合

1原子的总轨道角动量M L 为各电子轨道角动量之和。 π

2)

1(h

L L M L += L —原子总轨道角动量量子数;令各个电子的角动量量子数为li ;两个角动量的加和称为角动量偶合。即l —l 偶合,偶合规则为:

||,,1,212121l l l l l l L --++=Λ

原子轨道的角动量在Z 轴方向的分量为:

π

2)(h m M L

Z L = L m —原子磁量子数,L m L ±±±=,,2,1,0Λ

∑=i L m m ,组态中各电子磁量子数m 加和得到各可能的L m 。

2原子的总自旋角动量:

π

2)

1(h S S M S +=

S —原子的自旋量子数,是由两个电子自旋角动量偶合,即S —S 偶合,是整数或半整数。

原子总自旋角动量与磁场方向的分量: π

2)(h m M S

Z S = S m —原子的自旋磁量子数

)(,,2,1,0整数时S

m S ±±±=Λ

)(,,2

1

,21半整数时S

m S ±±±=Λ

∑=S S m M ,组态中各电子可能自旋磁量子数S m 加和得到各种可能的原

子自旋磁量子数。最小为零或1/2。

3 轨道和自旋的作用,将总轨道角动量和总自旋角动量相加得到总角动量的组合,称为L —S 偶合。因此,原子的总角动量为:

π

2)

1(h

J J M J +=

J —原子的总量子数或内量子数,根据角动量偶合规则,J 有如下数值: ||,,1,S L S L S L J --++=Λ

有时当个有时当)12(,)12(,+<+≥L J S L S J S L

总角动量在Z 轴方向的分量: π

2)(h m M J

Z J = )(,,2,1,0为整数J J

m J ±±±=Λ )(,,2,1,0为半整数J J

m J ±±±=Λ

2.3.2. 原子光谱项推求

用原子的量子数(S 、L 、T )来描述原子光谱项,规定:

光谱支项

光谱项J

S S L L

1

21

2++

Λ

ΛI H G F D P S L

6543210 12+S 称光谱的多重性

1等价电子(具有完全相同的主量子数n 和角量子数l ) 如 22p 组态: 2,1,0121=?==L l l

31121,02

1

21或=+=?==S S S S

可能有:D D P

P

S

S 3

1

3

1

3

1 (必须服从保里原理)

m ∑=i L m M

∑=S

S m M

光谱项

+1 0 -1

1

↑↓

2

1

D

2 ↑ ↑ 1 1 3

P 3 ↑ ↓ 1 0 1

D 4 ↓ ↑ 1 0 3

P 5 ↓ ↓ 1 -1 3

P 6 ↑ ↑ 0 1 3

P 7 ↑ ↓ 0 0 1

D 8 ↓ ↑ 0 0 3

P 9 ↓ ↓ 0 -1 3

P 10 ↑↓ 0 0 1

S 11 ↑ ↑ -1 1 3

P 12 ↑ ↓ -1 0 1

D 13 ↓ ↑ -1 0 3

P 14

-1

-1

3

P

15 ↑↓ -2 0

1

D

m

∑=i L m M ∑=S S m M

光谱项

+1

0 -1

1 ↑↓

2 0 1

D 2 ↑ ↑ 1 1 3

P 3 ↑ ↓ 1 0 1

D 4 ↓ ↑ 1 0 3

P 5 ↑ ↑ 0 1 3

P 6 ↑ ↓ 0 0 3

P 7 ↓ ↑ 0 0 1

D 8 ↑↓ 0 0 3

P 9 ↓ ↓ 1 -1 3

P 10 ↓ ↓ 0 -1 1

S 11 ↑ ↑ -1 1 3

P 12 ↑ ↓ -1 0 1

D 13 ↓ ↑ -1 0 3

P 14 ↓ ↓ -1 -1 3

P 15

↑↓

-2

1

D

种种10

)

12)(12(:

091,01,0)12)(12(:152,1,00)

12)(12(:21

1

1==++=±=±=++=±±==++=L S L S L S M M S L S

L M M S L P L M M S L D L )

(!!

v u v u C v u -=

u =每个电子可能存在的状态

15)!

26(!2!

62

62=-=

=C np

受保里原理限制,3D 、1P 、3S 不出现。

当 2,1,0,

2±±==L M L

a 、最大的2±=L M 时,两个电子已有三个量子数相同,则要求S m 不同。两个

电子须反平行,即:0=S m 表明:0=S 不应有 )(13D S =

2

1

1

211

2,

2

D L D L S L J J S S ??=+=++

在0=S m 的情况下,21D 包含 2,1,0±±=L 的5种状态(磁场中) (2J+1)=5 D 3 不出现的原因 b 、余下最大的 1,1,1====S L M M S L 属 03

1323

3,,0

,1,2P P P J P

=

)(9)102()112()122(状态=+?++?++?

c 、S M M S L S L 1

,

0,0====

01

S J =

d 、电子与空位无关, n 个电子与 n 个空位的组态光谱项相同。

8

2915

14

2nd nd nd nd np np np np ====

2不等价电子( n 和 l 至少有一个量子数不同) 如1,31

,232221111====l n l n p p 组态

S

P D S

P

D

S L 1

11

3

3

3

0,10,1,2==

○3光谱项的能量高低顺序 Hund 规则

ⅰ.原子在同一组态时,S 值最大者最稳定。 ⅱ.S 值相同时,L 值最大者最稳定。

ⅲ.一般地,L 和S 相同时,电子少于和等于半充满时,J 值小,能量低;电子多于半充满时,J 值大,能量低。

最新薛定谔方程及其解法

关于薛定谔方程 一.定义及重要性 薛定谔方程(Schrdinger equation)是由奥地利物理学家薛定谔提 出的量子力学中的一个基本方程,也是量子力学的一个基本假定, 其正确性只能靠实验来检验。是将物质波的概念和波动方程相结合 建立的二阶偏微分方程,可描述微观粒子的运动,每个微观系统都 有一个相应的薛定谔方程式,通过解方程可得到波函数的具体形式 以及对应的能量,从而了解微观系统的性质。 薛定谔方程是量子力学最基本的方程,亦是量子力学的一个基 本假定,它的正确性只能靠实验来检验。 二.表达式 三.定态方程 ()() 2 2 2 V r E r m η ψψ + ?? -?= ?? ?? 所谓势场,就是粒子在其中会有势能的场,比如电场就是一个带电粒子的势场;所谓定态,就是假设波函数不随时间变化。 其中,E是粒子本身的能量;v(x,y,z)是描述势场的函数,假设不随时间变化。

2 2 22222 z y x ??????++=? 可化为 d 0)(222 =-+ψψv E h m dx 薛定谔方程的解法 一. 初值解法;欧拉法,龙格库塔法 二. 边值解法;差分法,打靶法,有限元法 龙格库塔法(对欧拉法的完善) 给定初值问题 ). ()()((3) ) ,(),()( ,,(2) )(),( 311212 2111021h O t y t y hk y h t f k y t f k k c k c h y y y c c a y b t a y t f dt dy i i i i i i i i =-???????++==++==?????=≤≤=++的局部截断误差使以下数值解法的值及确定常数ββα βα

实验三 定态薛定谔方程的矩阵解法

实验三 定态薛定谔方程的矩阵解法 一.实验目的 1.掌握定态薛定谔方程的矩阵解法。 2.掌握几种矩阵特征值问题数值解法的原理,会调用相应的子程序求解具体问题。 二.实验内容 1.问题描述 以/2ω/()m ω为长度单位,一维谐振子的哈密顿量为 2 202d H x dx =-+, 其本征值为21n E n =+,本证波函数为 2 /2)()n n x H x ?=-, 其中()n H x 为厄米多项式,满足递推关系 11()2()2()n n n H x xH x nH x +-=-。 用矩阵方法求 2 22d H x x dx =-++ 的本证能量和相应的波函数。 2.问题分析 H E ψψ= 0()|j j j t c ψ?∞ ==>∑ 0||i i j i j i j c E c x Ec ??∞ =+<>=∑ 11|j j j x ???-+>=>>

11||||j j j j x x ????-+<>= <>= 0010010 112111,211,11,1 n n n n n n n n n n n n E x c c x E x c c E x E x c c x E c c -------?????????????????????????=??????????????????????? ? 3.程序编写 子程序及调用方法见《FORTRAN 常用算法程序集(第二版)》第三章 徐士良,P97 4.实验要求 ◆用恰当的算法求解以上实对称三对角矩阵的特征值问题。 ◆取n=8,给出H 的全部特征值和相应的特征向量。 5.实验步骤 ● 启动软件开发环境Microsoft Developer Studio 。 ● 创建新工作区shiyan03。 ● 创建新项目xm3。 ● 创建源程序文件xm3.f90,编辑输入源程序文本。 ● 编译、构建、运行、调试程序。 6.实验结果 程序设计:

大学物理-一维定态薛定谔方程的应用

一维定态薛定谔方程 的应用 授课人: 物理科学与技术学院

势 阱 日常生活中的各种井(阱) 物理学中研究微观粒子运动状态时常用的模型,因其势能函数曲线的形状如同井而得名 水井 窨井 陷阱 U x O a U

() U x x O a ∞ ∞00()0 , x a U x x x a ≤≤?=?∞<>? 这是一个理想化的物理模型, 应用定态薛定谔方程求解波函数, 有利于进一步理解在微观系统中 能量量子化和概率密度等概念 这样的势能函数称为 一维无限深势阱

建立定态薛定谔方程并求解 假设微观粒子质量为 ,由 m 22 2d ()()()2d U x x E x m x ψψ??-+=???? x a U x 0()0≤≤=阱内( ) : 22 2d ()()2d x E x m x ψψ-= x x a U x 0 , ()<>→∞ 阱外( ): 令: 2 22mE k =得通解: ()sin() x A kx ψ?=+ 微观粒子的能量不可能达到 无穷大,所以粒子不可能在阱外出现,或者说粒子在阱外出现的概率为零。 ()0 x ψ≡222 d 0d k x ψψ+=

利用标准条件确定 和 k ?因 在整个 轴上必须连续 x ()x ψsin() 0()0 0 0 A kx x a x x x ?ψ+≤≤?=? <>?,(0)sin 0 A ψ?== a A ka ()sin()0 ψ?=+=求归一化的波函数 一维无限深势阱中 微观粒子的波函数 2220π()d sin d a n x x A x x a ψ+∞-∞=??221 A a =?= 2A a = n a x x a x a x x a π2sin 0()00 , ψ? ≤≤?=??<>?() π ()sin 1,2,3n x A x n a ψ==??, 0?=π n k a =()1,2,3n =???,

第二章原子结构与性质§氢原子和类氢原子的薛定谔方程及其

第二章 原子结构 ?性质 §2.1.氢原子和 ?氢原子的 ?定谔方程 ?其解 2.1.1.单电子原子?的 定谔方?程 H 原子和H ?e +、Li 2+ 等 氢离子?是单原子,它们的核电?荷数为Z ,若把原子的?质量中心 ?在坐标原 ?上,绕核运动的?电子离核的?距离为r ,电子的电荷?为-e ,其静电作 ?势能为: r Ze V 0 2 4πε-= 将势能代 ? 定谔方程?: 得 0)(2 2282 =ψ+ + ψ?r Ze h m E π 或ψ=ψ- ?- E r Ze m h ][2 2 2 2 8π 为了解题方?便,将x 、y 、z 变量换 ?极坐标变量?r 、θ、φ。 其关系:φθcos sin r x = φθsin sin r y = φcos r z = 2222 z y x r ++= 2 1 ) /(cos 2 22z y x Z ++=θ x y tg /=φ })(sin )({2 222 sin 1sin 121 2 φθθθ θθ ??????????++= ? r r r r 代 定谔?方程: )()(sin )(2 2 22 222228sin 1 1sin 1121=ψ+ + ++???ψ??????? r Ze h m r r r r r E r πφθθθ θθ 2.1.2.分离变量§法:

上述的方程?是含三个 ?量的偏微分?方程,要解这个方?程可 数?分离法将其?化为三个分?别只含一个? 量的常微?分方程求解?。 含:)()()(),,(φθθΦΘ=Φψr R r 代 方程:并乘以ΘΦR r θ 22sin 移项 可 得: ) (s )(s )(228s i 2s i n 1222 2 2V E r r h u d d d d dr dR dr d R d d ---- =ΘΘΦΦθθ πθθ θθφ左边不含r ?、θ,右边不含 ?,欲左右两边?相等必等 ?同一个常数?(-m 2 ) Φ-=Φ 222m d d φ , 而右边可为?:(除以sin ?θ) )(sin )()(sin 1sin 821 22 2 2θ θ θθ πθ d d d d m h ur dr dR dr d R V E r ΘΘ-= -+ 则有: K d d d d m =-ΘΘ)(sin sin 1sin 22 θ θ θ θθ K E r r Ze h ur dr dR dr d R =++)()(2 2 2 2821 π 2.1. 3.方程解的 ?果 2.1.3.1.Φ(φ)方程的解 022 2=Φ+Φ m d d φ 这是一个常?系数二阶 ?次线性方程?,有两个复 ?数的独立解?。 |)|(]exp[m m im A m ±==Φφ Φ符合波 ?数品优条 ?:连续、单值、电子边界条? (归一) 1]exp[]exp[202 20 *=-A =ΦΦ?? φφφφπ πd im im d m m π21=A ]exp[][21φπ im m =Φ α、φ周期变化?,Φm 值不变 )2()(πφ φ+Φ=Φm m

薛定谔方程与提出背景

薛定谔方程 在一维空间里,一个单独粒子运动于位势中的含时薛定谔方程为 ;(1) 其中,是质量,是位置,是相依于时间的波函数,是约化普朗克常数,是位势。类似地,在三维空间里,一个单独粒子运动于位势中的含时薛定谔方程为 。(2) 假若,系统有个粒子,则波函数是定义于 -位形空间,所有可能的粒子位置空间。用方程表达, 。 其中,波函数的第个参数是第个粒子的位置。所以,第个粒子的位置是。 不含时薛定谔方程 不含时薛定谔方程不相依于时间,又称为本征能量薛定谔方程,或定态薛定谔方程。顾名思义,本征能量薛定谔方程,可以用来计算粒子的本征能量与其它相关的量子性质。 应用分离变量法,猜想的函数形式为 ; 其中,是分离常数,是对应于的函数.稍回儿,我们会察觉就是能量. 代入这猜想解,经过一番运算,含时薛定谔方程 (1) 会变为不含时薛定谔方程: 。 类似地,方程 (2) 变为

。 历史背景与发展 爱因斯坦诠释普朗克的量子为光子,光波的粒子;也就是说,光波具有粒子的性质,一种很奇奥的波粒二象性。他建议光子的能量与频率成正比。在相对论里,能量与动量之间的关系跟频率与波数之间的关系相同,所以,连带地,光子的动量与波数成正比。 1924年,路易·德布罗意提出一个惊人的假设,每一种粒子都具有波粒二象性。电子也有这种性质。电子是一种波动,是电子波。电子的能量与动量决定了它的物质波的频率与波数。1927年,克林顿·戴维和雷斯特·革末将缓慢移动的电子射击于镍晶体标靶。然后,测量反射的强度,侦测结果与X射线根据布拉格定律 (Bragg's law) 计算的衍射图案相同。戴维森-革末实验彻底的证明了德布罗意假说。 薛定谔夜以继日地思考这些先进理论,既然粒子具有波粒二象性,应该会有一个反应这特性的波动方程,能够正确地描述粒子的量子行为。于是,薛定谔试着寻找一个波动方程。哈密顿先前的研究引导著薛定谔的思路,在牛顿力学与光学之间,有一种类比,隐蔽地暗藏于一个察觉里。这察觉就是,在零波长极限,实际光学系统趋向几何光学系统;也就是说,光射线的轨道会变成明确的路径,遵守最小作用量原理。哈密顿相信,在零波长极限,波传播会变为明确的运动。可是,他并没有设计出一个方程来描述这波行为。这也是薛定谔所成就的。他很清楚,经典力学的哈密顿原理,广为学术界所知地,对应于光学的费马原理。借着哈密顿-雅可比方程,他成功地创建了薛定谔方程。薛定谔用自己设计的方程来计算氢原子的谱线,得到了与用玻尔模型计算出的能级相同的答案。 但是,薛定谔对这结果并不满足,因为,索末菲似乎已经正确地计算出氢原子光谱线精细结构常数的相对论性的修正。薛定谔试着用相对论的能量动量关系式,来寻找一个相对论性方程(现今称为克莱因-高登方程),可以描述电子在库仑位势的量子行为。薛定谔计算出这方程的定态波函数。可是,相对论性的修正与索末菲的公式有分歧。虽然如此,他认为先前非相对论性的部分,仍旧含有足够的新结果。因此,决定暂时不发表相对论性的修正,只把他的波动方程与氢原子光谱分析结果,写为一篇论文。1926年,正式发表于物理学界[2]。从此,给予了量子力学一个新的发展平台。 薛定谔方程漂亮地解释了的行为,但并没有解释的意义。薛定谔曾尝试解释代表电荷的密度,但却失败了。1926年,就在薛定谔第四篇的论文发表之后几天,马克斯·玻恩提出概率幅的概念,成功地解释了的物理意义[3]。可是,薛定谔本人一直不承认这种统计或概率的表示方法,和它所伴随的非连续性波函数坍缩。就像爱因斯坦的认为量子力学是基本为确定性理论的统计近似,薛定谔永远无法接受哥本哈根诠释。在他有生最后一年,他写给马克斯·玻恩的一封信,薛定谔清楚地表明了这看法。 含时薛定谔方程导引

第二章 薛定谔方程 习题

第二章 薛定谔方程 习题 (课本44页) 证明在定态中,概率流密度与时间无关。 证明:当一个系统处于定态时,其波函数),(t r ?可以写作, ?? ? ??-=Et i r t r ex p )(),(φ? 于是便有, ?? ? ??=Et i r t r ex p )(),(**φ? 根据概率流密度的定义式有, ?????????? ??-??? ??-???? ????? ??-= ??????????????? ??-???? ??-????????? ?????? ??-=?-?≡????????ψψψψt iE t iE t iE t iE m i t iE t iE t iE t iE m i m i J exp exp exp exp 2exp exp exp exp 2) (2* ***** 即有,)(2)(2****φφφφ?????-?=?-?=m i m i J 显然,在定态中概率流密度与时间无关。从某种意义上说明上述波函数称为定态波函数是名副其实的。 — 由下列两定态波函数计算概率流密度:⑴ )exp(1 1ikr r =?,⑵ )exp(12ikr r -=?。 从所得结果说明1?表示向外传播的球面波,2?表示向内(即向原点)传播的球面波。 解:在解本题之前,首先给出一个函数f 的梯度在球坐标系下的表达式,即 ? θθ? θ??+??+??=?f r e f r e r f e f r sin 1?1?? ⑴ 首先求解函数1?的概率流密度

r ikr ikr r ikr ikr ikr r ikr e mr k r ike r e e r e r ike r e e r e m i r ikr r ikr r ikr r ikr m i m i J ???2)exp()exp()exp()exp(2) (22221*1*111 =?????????? ??+--???? ??-+-=?? ? ???---?=?-?=---???? 可见,概率流密度1J 与r 同号,这便意味着1J 的指向是向外的,即1?表示向外传播的球面 波。 ⑵ 同理,可以得到2?的概率流密度 r ikr ikr r ikr ikr ikr r ikr e mr k r ike r e e r e r ike r e e r e m i r ikr r ikr r ikr r ikr m i m i J ???2)exp()exp()exp()exp(2) (22222*2*222 -=? ????????? ??---???? ??+-=?? ? ??-?-?-=?-?=---???? 这里的负号,即为概率流密度2J 与r 的符号相反,意味着概率流密度2J 的指向是向内 的,即波函数2?表示向内传播的球面波。 < 一粒子在一维势场 ?? ? ??>∞≤≤<∞=a x a x x x U ,0,00, )( 中运动,求粒子的能级和对应的波函数。 解:在量子力学中,一维薛定谔方程扮演着非常重要的角色。 其一,一维问题是微分方程中最简单、最基础的问题,通过解一维薛定谔方程,不但可以了解到量子力学中不同于经典力学的结果,如能量的量子化和势垒的贯穿等,还可以解更高维薛定谔方程的基础,如经典的氢原子的结构问题和现代的黑洞的结构问题,这些问题通过分离变量,最终化成求解一维薛定谔方程问题。 其二,随着现代科学技术的发展,在实验室中已经制成了一维的或准一维的系统,这样,求解一维薛定谔方程对于理解这些系统的性质起着至关重要的作用。 一维薛定谔方程的求解一般有两大类:一类是束缚态的求解,即求解束缚态的能级及相 x

薛定谔方程

第一章 薛定谔方程 §1.1.波函数及其物理意义 1. 波函数: 用波函数描述微观客体的运动状态。 例:一维自由粒子的波函数 推广 :三维自由粒子波函数 2. 波函数的强度——模的平方 3. 波函数的统计解释 用光栅衍射与电子衍射对比的方式理解波函数的统计解释。 t 时刻,出现在空间(x,y,z )点附近单位体积内的粒子数与总粒子数之比。 t 时刻,粒子出现在空间(x,y,z )点附近单位体积内的概率。 t 时刻,粒子在空间分布的概率密度 4、 波函数的归一化条件和标准条件 归一化条件 粒子在整个空间出现的概率为1 标准条件:一般情况下, 有关特殊情况波函数所满足的条件参看曾谨言教程。 对微观客体的数学描述: 脱离日常生活经验,避免借用经典语言引起的表观矛盾 §1.2. 薛定谔方程 是量子力学的基本假设之一,只能建立,不能推导,其正确性由实验检验。 1. 建立 (简单→复杂, 特殊→一般) 一维自由粒子的振幅方程 非相对论考虑 2. 一维定态薛定谔方程 2 |),,,(|t z y x ψ1d d d d d ||2===?=ψ???N N N N V V N N V V V . 是单值、有限、连续的ψ0)(2d )(d 222=ψ+ψx mE x x 0)()(2d )(d 222=ψ-+ψx U E m x x

3. 三维定态薛定谔方程 4. 一般形式薛定谔方程 5. 多粒子体系的薛定谔方程 讨论: 1、薛定谔方程也称波动方程,描述在势场U 中粒子状态随时间的变化规律。 2 、建立方程而不是推导方程,正确性由实验验证。薛定谔方程实质上是一种基本假设,不能从其他更基本原理或方程推导出来,它的正确性由它解出的结果是否符合实验来检验。 3、薛定谔方程是线性方程。是微观粒子的基本方程,相当于牛顿方程。 4、自由粒子波函数必须是复数形式,否则不满足自由粒子薛定谔方程。 5、薛定谔方程是非相对论的方程。 量子力学的中心任务就是求解薛定谔方程。 求解问题的思路: 1. 写出具体问题中势函数U (r )的形式代入方程 2. 用分离变量法求解 3. 用归一化条件和标准条件确定积分常数 4. 讨论解的物理意义, 薛定谔的另一伟大科学贡献 《What is life ?》 薛定谔(Schroding,1897-1961)奥地利人,因发现原子理论的有效的新形式一波动力学与狄拉克(Dirac,1902-1984)因创立相对论性的波动方程一狄拉克方程,共同分享了1933年度诺贝尔物理学奖 定态薛定谔方程 一.定态薛定谔方程条件:V (r,t )=V(r), 与t 无关。用分离变量法, 令Ψ=φ(r)f(t),代入薛定谔方程,得两个方程: 此称定态薛定谔方程 整个定态波函数形式: ),,,(),,,()],,,(2[),,,(2121212221t r r t r r V t r r m t r r t i i i i ψ+ψ?-=ψ??∑)t (Ef t )t (f i =?? Et i ce )t (f -=)r (E )r ()r (V )r (m ?=?+??-222Et i e )r ( -?=ψ

固体物理学 1-5-薛定谔方程应用举例II

薛定谔方程应用举例II---原子系统
? 氢原子 ? 电子自旋 ? 多电子原子
1

氢原子的定态薛定谔方程
?原子由一个原子核和核外电子构成,属于多粒子体系。多粒 子体系的总能量等于每一个粒子的能量与粒子间相互作用能量 之和。
?氢原子包括一个原子核和电子,库仑场是各向同性的,哈密 顿量可记作(绝热近似):
H?
=
?
h2 2me
?2
+
qeU(r)
me为电子质量,qe是电子电荷。U(r)为原子核静电场中的库 仑势,记作:
U(r) = ? Zqe = ? Z h2
4πε0r a1meqer
Z为核的电荷数,a1 = 4πε0?2/(meqe2) = 0.529?,为氢原子的第
一波尔轨道半径。
2

??? ?
h2 2me
?2
?
Zh 2 a1meqer
??ψ
?
(r)
=
E

(r)
中心力场问题,采用球坐标,薛定谔方程为:
? ?? ??
h2 2me
?
????
1 r2
? ?r
r2
? ?r
?
L?2 r2
???? ?
Zh2
?
?ψ (r,?,θ ) =
a1mer ??
E ?ψ (r,?,θ )
用分离变量法求解,令:
ψ (r,θ ,φ) = R(r) ?Y (?,θ )
分别求解径向波函数R(r)和角向波函数Y(?,θ)。
3

非线性薛定谔方程数值解的MATLAB仿真

[键入作者姓名] [键入文档标题] ——利用分步快速傅里叶变换对光纤中光信号的传输方程进行数值求解

1、非线性薛定谔方程 非线性薛定谔方程(nonlinear Schrodinger equation ,NLSE)是奥地利物理学家薛定谔于1926 年提出的,应用在量子力学系统中。由于量子力学主要研究粒子的动力学运动状态,所以不能运用牛顿力学公式来表示。通常在量子力学中,研究系统的状态一般通过波函数(x ,t)来表示。而对波函数的研究主要是求解非线性薛定谔方程。本文主要研究光脉冲在光纤中传输状态下的演变。 一般情况下,光脉冲信号在光纤中传输时,同时受到光纤的色散和非线性效应的影响。通过Maxwell 方程,考虑到光纤的色散和非线性效应,可以推导出光信号在光纤中的传输方程,即非线性薛定谔方程。NLSE 是非线性偏微分方程,一般很难直接求出解析解,于是通过数值方法进行求解。具体分为两大类:(1)分布有限差分法(split-step finite differencemethod ,SSFD);(2)分步傅里叶变换法(split-step Fourier transform method ,SSFT)。一般情况,在达到相同精度,由于分步傅里叶变换法采用运算速度快的快速傅里叶变换,所以相比较有限差分法运算速度快一到两个数量级。于是本文介绍分步傅里叶变换法来对光纤中光信号的传输方程,即非线性薛定谔方程进行数值求解。并通过MATLAB 软件对结果数值仿真。 非线性薛定谔方程的基本形式为: 22||t xx iu u u u =+ 其中u 是未知的复值函数. 目前,采用分步傅立叶算法(Split step Fourier Method)求解非线性薛定谔方程的数值解应用比较多。分步傅立叶方法最早是在1937年开始应用的,这种方法己经被证明是相同精度下数值求解非线性薛定愕方程最快的方法,部分原因是它采用了快速傅立叶变换算法(Fast Fourier Transform Algorithm)。基于MATLAB 科学计算软件以及MATLAB 强大的符号计算功能,完全可以实现分步傅立叶数值算法来对脉冲形状和频谱进行仿真。 一般情况下,光脉冲沿光纤传播时受到色散和非线性效应的共同作用,假设当传输距离 很小的时候,两者相互独立作用,那么,根据这种思想可建立如下分步傅立叶数值算法的数 学模型: 把待求解的非线性薛定谔方程写成以下形式: ??()U D N U z ?=+? (I ) (II )

第二章 薛定谔方程

第二章 薛定谔方程 本章介绍:本章将系统介绍波动力学。波函数统计解释和态叠加原理是量子力学的两个基本假设。薛定谔方程是波动力学的核心。在一定的边界条件和初始条件下求解薛定谔方程,可以给出许多能与实验直接比较的结果。 §2.1 波函数的统计解释 §2.1.1 波动—粒子两重性矛盾的分析按照德布罗意的观点,和每个粒子相联系的都有一个波。怎样理解粒子性和波动性之间的联系,这是量子力学首先遇到的根本问题。 2.1.1 波动—粒子两重性矛盾的分析能否认为波是由粒子组成? 粒子的单缝和双缝实验表明,如减小入射粒子强度,让粒子近似的一个一个从粒子源射出,实验发现,虽然开始时底片上的感光点是无规则的,但只要时间足够长,感光点足够多,底片上仍然会出现衍射条纹。如果波是由粒子做成,那末,波的干涉、衍射必然依赖于粒子间的相互作用。这和上述实验结果相矛盾,实际上,单个粒子也具有波动性的。 能否认为粒子是由波组成? 比如说,电子是三维空间的物质波包,波包的大小即电子的大小,波包的速度即电子的速度,但物质波包是色散的,即使原来的物质波包很小,但经过一段时间后,也会扩散到很大的空间去,或者形象地说,随着时间的推移,粒子将越来越“胖”,这与实验相矛盾 经典物理对自然界所形成的基本物理图像中有两类物理体系: ◆一类是实物粒子 ◆另一类是相互作用场(波)经典粒子是以同时确定的坐标和动量来描述其运动状态,粒子的运动遵从经典力学规律,在运动过程中具有确定严格的轨道。粒子的能量,动量在粒子限度的空间小区域集中;当其与其它物理体系作用时,只与粒子所在处附近的粒子相互作用,并遵从能量、动量的单个交换传递过程,其经典物理过程是粒子的碰撞;“定域”是粒子运动的特征。经典波动则是以场量(振幅、相位等)来描述其运动状态,遵从经典波动方程,波的能量和动量周期性分布于波所传播的空间而不是集中在空间一点,即波的能量、动量是空间广延的。波与其他物质体系相互作用时,可同时与波所在广延空间内的所有物理体系相互作用,其能量可连续变化,波满足叠加原理,“非定域”是波动性运动的特性。◆◆在经典物理中,粒子和波各为一类宏观体系的呈现,反映着两类对象,两种物质形态,其运动特点是不相容的,即具有粒子性运动的物质不会具有波动性;反之具有波动性运动的物质不会具有粒子性。综上所述,微观粒子既不是经典的粒子又不是经典的波,或者说它既是量子概念的粒子又是量子概念的波。其量子概念中的粒子性表示他们是具有一定的能量、动量和质量等粒子的属性,但不具有确定的运动轨道,运动规律不遵从牛顿定律;其量子概念中的波动性并不是指某个实在物理量在空间的波动,而是指用波函数的模的平方表示在空间某处粒子被发现的概率。◆现在被物理学家们普遍接受的波函数解释是玻恩提出的统计解释。他认为,粒子在衍射或干涉实验中所揭示的波动性质,既可以看成是大量粒子在同一实验中的统计结果,也可以认为是单个粒子在多次相同实验中显示的统计结果。 ◆玻恩的统计解释:波函数在某一时刻在空间的强度,即其振幅绝对值的平方与在这一点找到粒子的几率成正比,和粒子联系的波是概率波 §2.1.2 波函数统计解释 波函数的的特点:1.由于 2 |),(|t r ψ给出在 t 时刻,粒子在 r 处出现的几率密度,因此原 则上可由统计平均公式:? ?>= <)(r f 。在这种意义下,波函数),(t r ψ描述了微观粒子的运

非线性薛定谔方程的孤子解和怪波解

非线性薛定谔方程的孤子解和怪波解 摘要:光纤中光波的传输模型一直是当前研究的热点理论模型之一,从非线性薛定谔方程到金格堡-朗道方程,都试图对其进行更好的阐释,其次对于非线性动力学系统中,非线性薛定谔方程的解有呈现出非常多有趣的特征,对于其中特定解的研究能够让我们了解脉冲演化的本质,所以本文主要从孤子解的传输入手,并且简单介绍了怪波解的解形式。 薛定谔方程又称薛定谔波动方程,是量子力学的一个基本方程,同时又是量子力学的基本假设之一,由奥地利物理学家薛定谔1926年在《量子化就是本征值问题》中提出的,它在量子力学中的地位非常重要,相当于牛顿定律对于经典力学一样。 随着人们对世界的不断探索,非线性现象逐渐走进人们的视野,这种现象一般大都用非线性偏微分方程的数学模型来描述,显然线性方程已经不能满足人们的需求。 1973年,Hasegawa从含有非线性项的色散方程中推导出了非线性薛定谔方程。非线性薛定谔方程(NLS)是普适性很强的一个基本方程,最简单的形式是: 其中为常数。因为这个方程在几乎所有的物理分支及其他科学领域得到了广泛的应用,如超导,光孤子在光纤中传播,光波导,等离子体中的Langnui波等,所以许多学者对此方程的研究投入了很大的热情,至今还在生机勃勃的向前发展着。 1 分步傅里叶法计算演化过程 对于处理非线性性薛定谔方程,常用的数值仿真方式为分步傅里叶方法,为了简单起见,只考虑二阶色散和自相位调制,不考虑高阶色散、自陡以及四波混频等高阶非线性效应。上述方程中做 2 β为二阶色散,γ表示Kerr效应系数,g和α分别代表光纤中的增益和损耗。对上述方程转化到频域,先不考虑增益和损耗。可以得到 2 k k k k k dA i A i a a dz βγ =?+F. 其中2 2 2 k i β β ?=Ω 令() exp k k A B i z β =?可以得到 () 2exp k k k k dB i a a i z dz γβ =-? F 以上方程可以用四阶龙格库塔直接求解,但是速度较慢,所以我们需要做差分处理。 ()() ()()() 2 exp k k k k k B z z B z i a z a z i z z γβ +?- =-? ? F 再利用() exp k k A B i z β =?可以得到 ()()()() ()()() 2 2 exp exp exp k k k k k k k k A z z A i a z a z z i z a z i a z z i z γβ γβ ?? +?=+??? ?? ?? ?? ≈????? ?? F F 然后做傅里叶反变换就可以得到最终的结果 ()()()() 2 1exp exp - k k k k a z z a z i a z z i z γβ ?? +?=????? ?? F F

薛定谔方程及其解法

关于薛定谔方程 一. 定义及重要性 薛定谔方程(Schrdinger equation )是由奥地利物理 学家薛定谔提出的量子力学中的一个基本方程,也是量子力学的一个基本假定,其正确性只能靠实验来检验。是将物质波的概念和波动方程相结合建立的二阶偏微分方程,可描述微观粒子的运动,每个微观系统都有一个相应的薛定谔方程式,通过解方程可得到波函数的具体形式以及对应的能量,从而了解微观系统的性质。 薛定谔方程是量子力学最基本的方程,亦是量子力学的一个基本假定,它的正确性只能靠实验来检验。 二. 表达式 三. 定态方程 ()()2 22V r E r m ηψψ+??-?=???? 所谓势场,就是粒子在其中会有势能的场,比如电场就是一个带电粒子的势场;所谓定态,就是假设波函数不随时间变化。 其中,E 是粒子本身的能量;v(x ,y ,z)是描述势场的函数,假设不随时间变化。

2 2 22222z y x ?? ????++=? 可化为d 0)(222=-+ψψ v E h m dx 薛定谔方程的解法 一. 初值解法;欧拉法,龙格库塔法 二. 边值解法;差分法,打靶法,有限元法 龙格库塔法(对欧拉法的完善) 给定初值问题 ).()()((3) ) ,() ,() ( ,,(2) )() ,( 3112122111021h O t y t y hk y h t f k y t f k k c k c h y y y c c a y b t a y t f dt dy i i i i i i i i =-???????++==++==????? =≤≤=++的局部截断误差使以下数值解法 的值及确定常数ββα βα

大学物理练习题 氢原子理论 薛定谔方程

练习二十三 氢原子理论 薛定谔方程 一、选择题 1. 已知氢原子从基态激发到某一定态所需能量为10.19eV ,若氢原子从能量为?0.85eV 的状态跃迁到上述定态时,所发射的光子的能量为 (A ) 2.56eV 。 (B ) 3.41eV 。 (C ) 4.25eV 。 (D ) 9.95eV 。 2. 氢原子光谱的巴耳末系中波长最长的谱线用λ1表示,其次波长用λ2表示,则它们的比值λ1/λ2为 (A ) 9/8。 (B ) 19/9。 (C ) 27/20。 (D ) 20/27。 3. 根据氢原子理论,氢原子在n =5的轨道上的动量矩与在第一激发态的轨道动量矩之比为: (A ) 5/2。 (B ) 5/3。 (C ) 5/4。 (D ) 5。 4. 将波函数在空间各点的振幅同时增大D 倍,则粒子在空间的分布几率将 (A ) 增大D 2倍。 (B ) 增大2D 倍。 (C ) 增大D 倍。 (D ) 不变。 5. 一维无限深势阱中,已知势阱宽度为a 。 应用不确定关系估计势阱中质量为m 的粒子的零点能量为: (A ) ?/(ma 2)。 (B ) ?2/(2ma 2)。 (C ) ?2/(2ma )。 (D ) ?/(2ma 2)。 6. 由于微观粒子具有波粒二象性,在量子力学中用波函数Ψ(x ,y ,z ,t )来表示粒子的状态,波函数Ψ (A ) 只需满足归一化条件。 (B ) 只需满足单值、有界、连续的条件。 (C ) 只需满足连续与归一化条件。 (D ) 必须满足单值、有界、连续及归一化条件。 7. 反映微观粒子运动的基本方程是 (A ) 牛顿定律方程。 (B ) 麦克斯韦电磁场方程。 (C ) 薛丁格方程。 (D ) 以上均不是。 8. 已知一维运动粒子的波函数为 ()()?? ???==?0e x cx x kx ψψ00<≥x x 则粒子出现概率最大的位置是x =

薛定谔方程对氢原子的应用

(16.4.4) (16.4.5) (图16.4a )球极坐标 薛定谔方程对氢原子的应用 (一)氢原子的薛定谔方程 前一节讨论一维运动自由粒子的薛定谔方程及 其定态解.本节要讨论氢原子中电子的运动,这与 前一节有两点不同: (1)氢原子电子作三维空间运动,因此,薛定 谔方程(16.3.3)中的波函数ψ(x,t )应换成ψ(x,y,z,t ) 或ψ(r ,t ),而22x ??应换成=??+??+??222222z y x ▽2.此▽2称为拉普拉斯算符或拉氏算符. ??????<<的薛定谔方程三维运动自由粒子)c (v 222222222z y x )m 2/(t i ??+??+??=?=?ψ?-=?ψ? (16.4.1) (2)氢原子的电子不是自由粒子,它受到氢核的库仑力,此力的作用可用它们的电势能E p 表示.因此,氢原子电子的薛定谔方程可表示如下??,见〔附录16D 〕. ??????<<的薛定谔方程氢原子电子)c (v p 2p k p 22E )m 2/p (E E E E )m 2/(t i +=+=ψ+ψ?-=?ψ? (16.4.2) *(二)氢原子的定态薛定谔方程 定态解是解决氢原子各种问题的基础.参照(16.3.4)至(16.3.6)式,可把(16.4.2)式中的波函数ψ(r ,t )分离为空间部分u (r )和时间部分f (t ),并参照(16.3.10)式写出氢原子的定态薛定谔方程,见〔附录16E 〕. ψ(r ,t )=u (r )f (t ), f (t )=C /iEt e - (16.4.3) ??????<<的定态薛定谔方程氢原子电子)c (v r 4e E 0u )E E )(/m 2(u 02p p 22πε-==-+? 氢核的质量比电子的大得多,可认为氢核不动,电子绕核转动.其电势能可表成E p =-e 2/4πε0r .此势能E p 只与电子至氢核的距离r 有关,而与方向无关,即具有球对称性,应用球极坐标较为方便.如(图16.4a ),O 表氢核,e 表电子,r 为e 至O 的距离.θ为r 与z 轴的夹角,θ称天顶角或极角.?为r 在xOy 平面的投影与x 轴的夹角.故有 x=rsin θcos ?; y=rsin θsin ?; z=rcos θ (16.4.6) 拉氏算符 2222222z y x ??+??+??=?改用球坐标(r,θ,?)表示如下:?? ()() 22222222sin r 1sin sin r 1r r r r 1???θ+θ??θθ ??θ+????=?(16.4.7) 将此▽2算符代入(16.4.4)式,便得到以球坐标表示的氢原子定态薛定谔方程. ? 郭敦仁《量子力学初步》18—19,34—35页,1978年版. ? 程守洙、江之永编,王志符、朱讠永春等修订《普通物理学》第3册177—180页,1982年修订本. ? 郭敦仁《量子力学初步》35—45页,1978年版. ? 周世勋编《量子力学》59—72页,1961年版.

§16.3 一维定态薛定谔方程的建立和求解举例

§16.3 一维定态薛定谔方程的建立和求解举例 (一)一维运动自由粒子的薛定谔方程 波函数随时间和空间而变化的基本方程,是薛定谔于1926年提出的,称为薛定谔波动方程,简称波动方程或薛定谔方程,它成为量子力学的基本方程. 将(16.2.14)式分别对t 和x 求导,然后从这两式消去E 、p 、和ψ,便可得到一维运动自由粒子的薛定谔方程: ψ-=?ψ?)/iE (t 即ψ=?ψ?E t i (16.3.1) ψ=?ψ ?22)/ip (x 2 ψ=ψ ?-2222p ????? ?????<<的薛定谔方程自由粒子轴运动的沿)c x (v 方程(16.3.3)中不含有能量E 和动量p ,表明此方程是不受E 和p 的数值限制的普遍方程. 请同学们自己试一试,如果上述波函数不用复数表式(16.2.14),改用类似于(16.2.1)式的余弦函数或正弦函数表式,就不会得到合乎要求的薛定谔方程(16.3.3)式?. 这薛定谔方程不是根据直接实验结果归纳而得,也不是由经典波动理论或其他理论推导出来的,它是在物质波假设的基础上,参照经典波动方程而建立起来的.薛定谔方程在微观领域中得到广泛的应用,它推导出来的结果,都与相关实验结果符合得很好,这才是薛定谔方程正确反映微观领域客观规律的最有力的证明. (二)一维运动自由粒子的定态薛定谔方程?? 上述薛定谔方程(16.3.3)是偏微分方程,从此方程可解出波函数ψ(x ,t ).在量子力学中最重要的解,是可把波函数ψ(x,t )分离成空间部分u (x )和时间部分f (t )两函数的乘积的特解,即 〔一维运动自由粒子的定态波函数〕 ψ(x,t )=u (x )f (t )(16.3.4) 将此式代入(16.3.3)式得: 22 2dx u d )t (f )m 2/(dt df )x (u i -= 两边除以ψ=uf 得: 22 2dx u d u 1)m 2/(dt df f 1i -= 此式左边是时间t 的函数,右边是坐标x 的函数.已知t 与x 是互相独立的自变量,左右两边相等,必须是两边都等于同一常量E ,即 ? 郭敦仁《量子力学初步》16—17页,人民教育出版社1978年版. ? 郭敦仁《量子力学初步》21—22页,人民教育出版社1978年版. ? 周世勋编《量子力学》32—33页,上海科学技术出版社1961年版.

求解非线性薛定谔方程的一类数值解法

求解非线性薛定谔方程的一类数值解法 张艳敏,刘明鼎 (青岛理工大学琴岛学院,山东青岛266106) 摘要:利用非标准有限差分方法构造了求解非线性薛定谔方程的两个非标准有限差分格式。对于离散后的差分格式,把关于时间和空间的步长函数作为分母逼近导数项。对于非线性项,通过非局部的离散方法计算了这两个非标准有限差分格式的局部截断误差。数值实验结果验证了非标准有限差分格式的有效性。关键词:非线性薛定谔方程;局部截断误差;数值解法中图分类号:O241.82 文献标识码:A 文章编号:2095-7726(2019)03-0008-03 薛定谔方程是物理学中量子力学的一个重要方程,可以用于研究深水波浪理论。柱(球)非线性薛定谔方程常用于描述单色波的一维自调适、光学的自陷现象、固体中的热脉冲传播和等离子体中的Langnui 波[1–5],因此对于此类方程的研究具有非常重要的意义。 薛定谔方程有线性和非线性两种,在本文中,我们研究的是非线性薛定谔方程。非线性薛定谔方程解的解析表达式是很难得到的,因此求解此类方程最常用的就是数值解法。求非线性薛定谔方程数值解的方法主要有差分方法、配置谱方法[6]、有限元方法[7]和平均离散梯度方法[8]等。在本文中,我们利用非标准有限差分方法研究了非线性薛定谔方程的数值解,这种方法已在求解偏微分方程中得到了广泛的应用[9],其优点是对非线性项作非局部离散,对导数项作离散后用步长函数作分母,这样不仅能保持差分方程的数值解与原方程的解析解具有相同的正性,而且能保持较好的数值稳定性。 1非标准有限差分格式的构造 现在我们利用文献[10-12]给出的方法构造非线 性薛定谔方程的两种非标准有限差分格式,要考虑的非线性薛定谔方程为 (1)相应的初边值条件为 其中:为虚数单位;、、和均为连续函数;和均为正数。 为了得到非线性薛定谔方程的差分格式,需要对式 (1)进行离散。首先,需要利用网格对区域进行分割,取空间步长时间步长其 次,在网格点处,定义数值解其中,且下面将分别构造式(1)的两种非标准有限差分格式。 1.1第一种非标准有限差分格式的构造 为了构造式(1)的第一种非标准有限差分格式,我们利用R.E.Mickens 提出的构造非标准有限差分格式的原理[10]和文献[13-14]中提到的方法,并利用给定的记号,对式(1)进行离散。离散后的差分方程为 其中,和为分母函数,且,且分母是通过步长函数逼近得到的。 从式(4)可以看出,和分别取代了和分母函数的选择依据了薛定谔方程解的性质[4]。 记对式(4)进行整理,可得第36卷第3期Vol.36No.3 新乡学院学报 Journal of Xinxiang University 2019年3月Mar.2019 收稿日期:2018-12-21 基金项目:山东省高校科技计划项目(J17KB053);青岛理工大学琴岛学院教育教学研究重点项目(2018003A)作者简介:张艳敏(1981—),女,山东东营人,副教授,硕士,研究方向:偏微分方程数值分析。通信作者:刘明鼎(1982—),男,辽宁大连人,副教授,硕士,研究方向:偏微分方程数值分析。 222 (,)(,) (,)(,)(,),i u x t u x t u x t u x t g x t t x ??=++??(,0)(), u x f x =(2) 01(0,)(), (,)()u t p t u L t p t =ìí =?。 (3)0,0;x H t T <£<£i (,)g x t ()f x 0()p x 1()p x H T [0,][0,]H T ′,h H M =Δt T N =。(,)m n x t (,),n m m n u u x t =(0,1,2,,),m x mh m M ==L Δ(0,1,2,,),n t n t n N ==L ,M N ++??Z Z 。111212 2(),i n n n n n m m m m m n n n m m m u u u u u u u g D D ++---+=++(4) 1D 2D 12exp(Δ)1,D t D =-=24sin ()2 h 1D 2D Δt 2,h 11122 ,,D R D R D ==

相关文档
最新文档