黄昆固体物理习题-第三章 晶体的热性质
《固体物理·黄昆》第三章

氢键结合的情况可写成通式:
X-H…Y。 式中 X 、 Y 代表 F 、 O 、 N 等电负 性大而原子半径较小的非金属原 子, X 和 Y 可以是两种相同的元 素,也可以是两种不同的元素。 d F l H F H F
归纳起来,氢键形成的条件是:
A)有与电负性大(X)的原子相结合的氢原子;
B) 有一个电负性也很大,含有孤对电子并带有部分负 电荷的原子(Y); C)X与Y的原子半径都要较小。
氯化钠型 —— NaCl、KCl、AgBr、PbS、MgO (配位数6) 氯化铯型 —— CsCl、 TlBr、 TlI(配位数8)
离子结合成分较大的半导体材料ZnS等(配位数4)
2. 离子晶体结合的性质
1) 系统内能的计算 晶体内能 : 1)所有离子库仑相互作用能(吸引作用)
2) 和重叠排斥能之和(排斥作用)
具体晶体的内聚能(晶格能)参见周期表,有一定的规律性: 惰性气体晶体<碱金属<过渡族金属(共价晶体)
两粒子间的相互作用 相互作用能.
f(r) 和u(r)分别表示相互 作用力和相互作用势 则:
u (r ) f (r ) r
U 排斥 r
f (r )
B rn
u (r )
pij A12= j'
12
12.13188
pij A6= j'
6
14.45392
物理意义:
晶体总的势能:
—— 非极性分子晶体的晶格常数、结合能和体变模量 晶格常数
平衡状态体变模量
晶体的结合能
分子晶体: 常温下是气态的物质如:Cl2,SO2,HCl, H2, O2, He, Ne, Ar, Xe等在低温下依靠范德瓦耳斯力结合成的晶体.
固体物理习题与答案

《固体物理学》习题解答黄昆 原著 韩汝琦改编 (志远解答,仅供参考)第一章 晶体结构1.1、解:实验表明,很多元素的原子或离子都具有或接近于球形对称结构。
因此,可以把这些原子或离子构成的晶体看作是很多刚性球紧密堆积而成。
这样,一个单原子的晶体原胞就可以看作是相同的小球按点阵排列堆积起来的。
它的空间利用率就是这个晶体原胞所包含的点的数目n 和小球体积V 所得到的小球总体积nV 与晶体原胞体积Vc 之比,即:晶体原胞的空间利用率, VcnVx = (1)对于简立方结构:(见教材P2图1-1)a=2r , V=3r 34π,Vc=a 3,n=1 ∴52.06r 8r34a r 34x 3333=π=π=π= (2)对于体心立方:晶胞的体对角线BG=x 334a r 4a 3=⇒= n=2, Vc=a 3∴68.083)r 334(r 342a r 342x 3333≈π=π⨯=π⨯= (3)对于面心立方:晶胞面对角线BC=r 22a ,r 4a 2=⇒= n=4,Vc=a 374.062)r 22(r 344a r 344x 3333≈π=π⨯=π⨯= (4)对于六角密排:a=2r 晶胞面积:S=6260sin a a 6S ABO ⨯⨯=⨯∆=2a 233 晶胞的体积:V=332r 224a 23a 38a 233C S ==⨯=⨯ n=1232126112+⨯+⨯=6个 74.062r224r 346x 33≈π=π⨯= (5)对于金刚石结构,晶胞的体对角线BG=3r 8a r 24a 3=⇒⨯= n=8, Vc=a334.063r 338r 348a r 348x 33333≈π=π⨯=π⨯=1.2、试证:六方密排堆积结构中633.1)38(a c 2/1≈= 证明:在六角密堆积结构中,第一层硬球A 、B 、O 的中心联线形成一个边长a=2r 的正三角形,第二层硬球N 位于球ABO 所围间隙的正上方并与这三个球相切,于是: NA=NB=NO=a=2R.即图中NABO 构成一个正四面体。
固体物理第三章

固体物理第三章班级成绩学号Chapter 3 晶格振动与晶体的热学性质姓名(lattice vibration and its heat characteristics)⼀、简要回答下列问题(answer the following questions):1、在晶格常数为a 的⼀维单原⼦晶格中,波长λ=8a 和波长λ=8a/5的格波所对应的原⼦振动状态有⽆不同? 试画图加以说明。
[答]对于⼀维单原⼦链,由q=2π/λ知,λ=8a 时,q =π/4a ,λ=8a /5时,q =5π/4a ,⼆者的aq 相差π,不是2π的整数倍,因此,两个格波所对应的原⼦振动状态不同。
如上图,当两个格波的位相差为2π的整数倍时,则它们所对应的原⼦的振动状态相同。
2、什么叫简正振动模式?简正振动数⽬、格波数⽬或格波振动模式数⽬是否是⼀回事?[答]在简谐振动下,由N 个原⼦构成的晶体的晶格振动,可等效成3N 个独⽴的谐振⼦的振动,每⼀个谐振⼦的振动模式称为简正振动模式。
格波振动通常是这3N 个简正振动模式的线性叠加。
简正振动数⽬、格波数⽬或格波振动模式数⽬是是⼀回事,其数⽬等于晶体中所有原⼦的⾃由度之和,即等于3N 。
3、晶体中声⼦数⽬是否守恒?在极低温下,晶体中的声⼦数与温度T 之间有什么样的关系?[答]频率为ωi 的格波的平均声⼦数为: 11)(/-=Tk i B en ωω即每⼀个格波的声⼦数都与温度有关,因此晶体中的声⼦数⽬不守恒,它随温度的改变⽽改变。
以德拜模型为例。
晶体中的声⼦数⽬为ωωωωd g n N D)()('0=其中令 T k x B ω= 则 123'2/033233-=x TB e dxx C T k V N D θπ在极低温度下,θD /T →∞,于是 331332332033233)2(23123'T nC T Vk e dx x C T k V N n B x B ∑∞=∞=-=ππ即在温度极低时,晶体中的声⼦数⽬与T 3成正⽐。
黄昆固体物理习题-第三章_晶体的热性质

第三章习题参考解答3.1已知一维单原子链,其中第j 个格波,在第n 个格点引起的位移μnj 为:δj 为任意位相因子。
并已知在较高温度下每个格波的平均能量为kT ,具体计算每个原子的平方平均位移。
)sin(j j j j nj naq t δωαμ++=21)(sin 102=++⎰dt q n t T j j j Tδαω根据=2nj μ22221)(sin j j j j j q n t αδαωα=++解:其中T =2π/ωj 为振动周期,所以:格波的平均动能:∑∙=n njm E 221μN m j j 2241ωα=一维单原子链可以认为是经典的简谐运动,因此有:)(cos 21222j j j j n j q n t m δαωωα++=∑平均动能=平均势能= 格波平均能量=kT 2121其中:M =ρL其中振幅222j j Nm kT ωα=得:kT N m E j j 214122==ωα所以有:22221jj nj Nm kT ωαμ==所以,每个原子的平方平均位移:∑∑∑===222121j j njn Nm kT ωαμμ其中:M =ρL3.2 讨论N个原胞的一维双原子链(相邻原子间距为a),其2N个格波解,当M=m时与一维单原子链结果一一对应。
解:质量为M的原子位于2n-1,2n+1,2n+3……。
质量为m的原子位于2n,2n+2,2n+4 ……。
牛顿运动方程体系有N个原胞,有2N个独立的方程方程的解:A,B有非零解可以得到:两种不同的格波的色散关系为:对应一个q有两支格波:一支声学波和一支光学波,总的格波数目为2N。
当M=m时,色散关系简化为:长波极限情况下与一维单原子晶格格波的色散关系一致3.3质量相同两种原子形成一维双原子链,最近邻原子间的力常数交错等于c和10c,令两种原子的质量相等,并且最近邻间距是a/2,试求在k=0和k=π/a处的ω(k)。
并粗略。
画出色散关系。
固体物理基础第3章-晶格振动与晶体的热学性质

3-2 一维单原子链模型
格波的色散关系 4 2 2 aq sin ( )
m 2 • ω取正值,则有 (3)
(q)
aq 2 sin( ) m 2 • 频率是波数的偶函数
• 色散关系曲线具有周期性, 仅取简约布里渊区的结果即可 • 由正弦函数的性质可知,只有满足 0 2 / m 的格波 才能在一维单原子链晶体中传播,其它频率的格波将被强
原子n和原子n+1间的距离
非平衡位置
原子n和原子n+1间相对位移
a n1 n
n1 n
3-2 一维单原子链模型
• 忽略高阶项,简谐近似考虑原子 振动,相邻原子间相互作用势能 1 d 2v v(a ) ( 2 ) a 2 2 dr • 相邻原子间作用力 dv d 2v f , ( 2 )a d dr • 只考虑相邻原子的作用,第n个原 子受到的作用力
• 连续介质中的波(如声波)可表示为 Ae ,则可看出 • 格波和连续介质波具有完全类似的形式 • 一个格波表示的是所有原子同时做频率为ω的振动 • 格波与连续介质波的主要区别在于(2)式中,aq取值任意加减 2π的整数倍对所有原子的振动没有影响,所以可将波数q取值 限制为 q a a
V
O
a
r
• 第n个原子的运动方程
(n1 n ) (n n1 ) (n1 n1 2n )
(1)
平衡位置
d 2 n m 2 ( n1 n 1 2n ) dt
非平衡位置
——牛顿第二定律F=ma
3-2 一维单原子链模型
• 上述(1)式的解(原子振动位移)具有平面波的形式
a
)
黄昆版固体物理学课后答案解析答案 (2)

《固体物理学》习题解答黄昆 原着 韩汝琦改编 (陈志远解答,仅供参考)第一章 晶体结构、解:实验表明,很多元素的原子或离子都具有或接近于球形对称结构。
因此,可以把这些原子或离子构成的晶体看作是很多刚性球紧密堆积而成。
这样,一个单原子的晶体原胞就可以看作是相同的小球按点阵排列堆积起来的。
它的空间利用率就是这个晶体原胞所包含的点的数目n 和小球体积V 所得到的小球总体积nV 与晶体原胞体积Vc 之比,即:晶体原胞的空间利用率, VcnVx = (1)对于简立方结构:(见教材P2图1-1) a=2r , V=3r 34π,Vc=a 3,n=1∴52.06r8r34a r 34x 3333=π=π=π=(2)对于体心立方:晶胞的体对角线BG=x 334a r 4a 3=⇒= n=2, Vc=a 3∴68.083)r 334(r 342a r 342x 3333≈π=π⨯=π⨯= (3)对于面心立方:晶胞面对角线BC=r 22a ,r 4a 2=⇒= n=4,Vc=a 3(4)对于六角密排:a=2r 晶胞面积:S=6260sin a a 6S ABO ⨯⨯=⨯∆=2a 233 晶胞的体积:V=332r 224a 23a 38a 233C S ==⨯=⨯ n=1232126112+⨯+⨯=6个(5)对于金刚石结构,晶胞的体对角线BG=3r 8a r 24a 3=⇒⨯= n=8, Vc=a 3、试证:六方密排堆积结构中633.1)38(ac 2/1≈=证明:在六角密堆积结构中,第一层硬球A 、B 、O 的中心联线形成一个边长a=2r 的正三角形,第二层硬球N 位于球ABO 所围间隙的正上方并与这三个球相切,于是:NA=NB=NO=a=2R.即图中NABO 构成一个正四面体。
…、证明:面心立方的倒格子是体心立方;体心立方的倒格子是面心立方。
证明:(1)面心立方的正格子基矢(固体物理学原胞基矢):123()2()2()2a a j k a a i k a a i j ⎧=+⎪⎪⎪=+⎨⎪⎪=+⎪⎩r r r r r rr r r由倒格子基矢的定义:1232()b a a π=⨯Ωr r r31230,,22(),0,224,,022a a a a a a a a a a Ω=⋅⨯==r r rQ ,223,,,0,()224,,022i j ka a a a a i j k a a ⨯==-++r rr r r r r r同理可得:232()2()b i j k ab i j k aππ=-+=+-r rr r r r r r 即面心立方的倒格子基矢与体心立方的正格基矢相同。
固体物理课后习题解答(黄昆版)第三章

固体物理课后习题解答(黄昆版)第三章黄昆固体物理习题解答第三章晶格振动与晶体的热学性质3.1 已知⼀维单原⼦链,其中第j个格波,在第个格点引起的位移为,µ= anj j sin(ωj_j+ σj) ,σj为任意个相位因⼦,并已知在较⾼温度下每个格波的平均能量为,具体计算每个原⼦的平⽅平均位移。
解:任意⼀个原⼦的位移是所有格波引起的位移的叠加,即µn= ∑ µnj=∑ a j sin(ωj t naq j+σj)j j(1)µ2 n =∑µjnj∑µj*nj=µj2nj+ µ µnj*nj′j j′由于µ µnj?nj数⽬⾮常⼤的数量级,⽽且取正或取负⼏率相等,因此上式得第2 项与第⼀项µ相⽐是⼀⼩量,可以忽略不计。
所以2= ∑ µ 2njn j由于µnj是时间的周期性函数,其长时间平均等于⼀个周期内的时间平均值为µ 2 = 1 T∫0 2 ω+σ 1 2 j aj sin( t naqjj j)dt a=j(2)T0 2已知较⾼温度下的每个格波的能量为KT,µnj的动能时间平均值为1 L T ?1 ?dµ 2 ?ρw a2 T 1= ∫∫dx0?ρnj?= j j∫0 2 ω+ σ= ρ 2 2 T??dt L a sin( t naq)dt w Lanj T0 0 0 ? 2 ?dt??2T0 j j j j 4 j j其中L 是原⼦链的长度,ρ使质量密度,T0为周期。
1221所以Tnj= ρ w La j j=KT(3)4 2µKT因此将此式代⼊(2)式有nj2 = ρωL 2 jµ所以每个原⼦的平均位移为2== ∑ µ 2= ∑KT= KT∑1n njρωL2ρLω2j j j j j3.2 讨论N 个原胞的⼀维双原⼦链(相邻原⼦间距为a),其2N 格波解,当M=m 时与⼀维单原⼦链的结果⼀⼀对应.解答(初稿)作者季正华- 1 -黄昆固体物理习题解答解:如上图所⽰,质量为M 的原⼦位于2n-1,2n+1,2n+3 ……质量为m 的原⼦位于2n,2n+2,2n+4 ……⽜顿运动⽅程:..mµ2n= ?βµ(22n?µ2n+1 ?µ2n?1)..Mµ2n+1 = ?βµ(22n+1 ?µ2n+2 ?µ2n)体系为N 个原胞,则有2N 个独⽴的⽅程i na q⽅程解的形式:iµ2n=Ae[ωt?(2 ) ] µ2n+1=Be[ω?(2n+1)aq]na qµ=将µ2n=Ae[ωt?(2 ) ]2n+1 Be i[ωt?(2n+1) aq]代回到运动⽅程得到若A、B 有⾮零的解,系数⾏列式满⾜:两种不同的格波的⾊散关系:——第⼀布⾥渊区解答(初稿)作者季正华- 2 -第⼀布⾥渊区允许 q 的数⽬黄昆固体物理习题解答对应⼀个 q 有两⽀格波:⼀⽀声学波和⼀⽀光学波。
黄昆版固体物理习题课_03

01/34
设试探解: 将试探解代入方程得到: 由线性齐次方程组有非零解的条件得到:
习题问题讨论三 ——晶格振动与晶体的热力学性质
02/34
当:
代入原方程组得到: 光学支:
声学支:
当:
光学支:B=0 声学支:A=0
习题问题讨论三 ——晶格振动与晶体的热力学性质
03/34
补充题二、设有一纵波:
沿着一维单原子链传播,原子间距为a,最近邻忽作用的恢复 力常数为 试证明:每个原子对时间平均的总能量为:
考虑到晶体中存在有许多不同频率、不同模式的格波,因此总 的均方位移应对所有不同格波进行求和。又由于各振动模式间 是相互独立的,因此有:
当N足够大时,振动频率趋于连续,求和可以用积分代替
习题问题讨论三 ——晶格振动与晶体的热力学性质
11/34
将德拜模型的频率分布函数及最大频率代入得:
(2)非零温度下相应于某频率的格波的平均能量应为格波能 量和该温度下该格波的平均声子数之积,即:
红色标记原子位于2n, 2n+2, 2n+4 ……
习题问题讨论三 ——晶格振动与晶体的热力学性质
19/34
—— 第2n个原子和第2n+1个原子的运动方程
—— 体系N个原胞,有2N个独立的方程 —— 方程的解
令
习题问题讨论三 ——晶格振动与晶体的热力学性质
20/34
—— A、B有非零的解,系数行列式满足
补充题一、证明在由两种不同质量M,m(M>m)的原子所组成的
一维复式格子中,如果波矢q取边界值
(a为相邻原子
间距),则在声学支上质量为m的轻原子全部保持不动;在光
学支上质量为M的重原子保持不动。
解:如图所示
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
(10 eiaq ) 2M 11
解之得:
2
M
11
121
20(1
c
os
aq)
1 2
当 q 时,cos 1
a
可得
2
M
1
2
1
20
M
2
当q=0时,则有:
0
1
22
M
2
色散关系如图所示:
22
M
20
M
2
M
3.4题目:
略。
系。
3.4 解:(1)设μl,m代表第(l,m)个原子,即第l
要激发的声子所用的电磁波波长在近红外线波段
2k mM
mM
其中,mNa=23 kg
代入数计算,可得: 11.11013 Hz
计算波长,可得:
2
q
2
c
2c
代入数据,计算波长,可得:
17m
与实测同数量级。因为得到的色散关系存在近似。
3.6 计算一维单原子链的频率分布函数 解: 设单原子链长度
2
2
所以有:
E
1 4
m
2j
2 j
N
1 2
kT
其中振幅
2 j
2kT
Nm
2 j
得:
2 nj
1 2
2 j
kT
Nm
2 j
所以,每个原子的平方平均位移:
n
2 nj
1
2
2 j
kT Nm
1
2 j
其中:M=rL
3.2 讨论N个原胞的一维双原子链(相邻原子间距为a),其 2N个格波解,当M=m时与一维单原子链结果一一对应。
j
)dt
1 2
其中T=2π/ωj为振动周期,所以:
2 nj
2 j
sin
2
(
jt
nq
j
j
)
1 2
2 j
格波的平均动能:
E
n
1 2
m
•
2 nj
1m 2
n
2j
2 j
cos
2
(
jt
nq
j
j)
1 4
m
2j
2 j
N
其中:M=rL
一维单原子链可以认为是经典的简谐运动,因此有:
平均动能=平均势能= 1 格波平均能量= 1 kT
cos aky )
这就是色散关系。
(c)k=kx,ky=0, kx=ky的-k图像为:
由(c)的结果代入可得到
3.5题:已知NaCl晶体平均每对离子的相互作用能为:
U (r) q2 / r / r n
其中马德龙常数α=1.75, n=9,平均离子间距为
r0=2.82Å。
(1)试求离子在平衡位置附近的振动频率; (2)计算与频率相当的电磁波的波长,并与NaCl红 外吸收频率的测量值61μm进行比较。
系是:ω=vq,在二维波矢空间内,格波的等频面是一个园, 如图所示,在q—q+dq区间内波速为v的格波数目为:
qy dq
dz
S
(2 )2
,
2qdq
Sd 2v 2
q
o
qx 式中S为二维晶格的总面积,由此可以得
到波速为v的格波的模式密度为:
g ( )
dz
d
S 2v2
由此可以得E为:
3.9题: 写出量子谐振子系统的自由能,证明在经典 极限,自由能为:
解:绿色标记的原子位于2n-1, 2n+1, 2n+3 …… 红色标记原子位于 2n, 2n+2, 2n+4 ……
第2n个原子和第2n+1个原子的运动方程可以写为: 体系N个原胞,有2N个独立的方程 方程的解:
令
A、B有非零的解,系数行列式满足
两种色散关系
—— 两种色散关系
色散关系(2--k)图如 右,这是一个双原子 (例如H2)晶体。
解: 1)光学波的最大频率
光学波的最小频率 声学波的最大频率
2)相应声子的能量
3) 某一特定谐振子具有激发能
的几率
归一化常数
—— 归一化条件
频率为谐振子的平均能量 频率为谐振子的能量
第i个q态的平均数声子
声子数目为:
4)如果用电磁波激发光学波 要激发 的声子所用的电磁波波长在什么波段? 对应电磁波的能量和波长
将上式代入F式中,便得到
F U0 kBT
q
ln( q )
kBT
3.10题:设晶体中每个振子的零点振动能为:1 ,
试用德拜模型求晶体的零点振动能。
2
解:根据量子力学,零点能是谐振子所固有的,与温
度无关,故T =0 K时的振动能E0就是各振动模零点能
之和:
E0
D 0
E0
()g()d
将
E0
1 ,
解:质量为M的原子位于 2n-1, 2n+1, 2n+3 ……。 质量为m的原子位于 2n, 2n+2, 2n+4 ……。
牛顿运动方程 体系有N个原胞,有2N个独立的方程 方程的解:
A , B有 非零解
可以得到: 两种不同的格波的色散关系为:
对应一个q有两支格波:一支声学波和一支光学波, 总的格波数目为2N 。
(,m1 ,m1 2,m )]
设试探解的形式为:
(0)ei[t( akx maky )] ,m
代入运动方程,消去公因子后,得:
M 2 c(eiakx eiakx eiaky eiaky 4)
2c(2 cos akx cos aky )
据此得色散关系 :
2
2c M
(2
cos akx
2
g()
3V 2 2 2C 3
代入上式,
E0
3V
2 2 C 3
. 2
m 3d
0
3Vm4 16 2 C 3
其中V是晶体的体积
由
CV
(T )
d E(T ) dT
3Vk
2 2 C 3
.
m 0
k
T
2
e
k
T
e
k
T
2
1
2d
式给出ωm
1
m
C 6 2
N V
3
把 m 代入上式,便可得到晶体的零点能:
行 m列的原子垂直于晶格平面的位移,当只考虑最 近邻原子间的相互作用时,由于(l +1,m)原子对 它的作用力:
f1 c( 1,m ,m )
第l-1,m原子对它的作用力:
f2 c( ,m ) 1,m
而f1于f2 方向是相反的。同样处理( l,m+1)原 子和( l,m-1)原子对(l,m)原子的作用力f3,f4, 于是得到第(l,m)个原子所受的力:
3.5解:由 U (r) q2 / r / r n
dU 0 dr
可以得到β:
1 n
q
2
r0n1
利用题目所给条件(1),可知:
f
k
dU dr
dU
d
d 2U dr 2
计算、简化,代入β,可得回复力系数k:
k
(n
1)q2
r03
=143N/m
对于NaCl晶体,可以认为是一维双原子晶体,由 黄昆教材p96-97分析可以知道,对于离子晶体只有长 光学波可以和电磁波发生相互作用。所以有
E0
3
16 2
V
3
C
(6 2
N
)
4 3
C
4
V
9 8
NkB (6
2
N V
1
)3
C kB
9 8
R D
3.11题: 一维复式格子中,如果
计算
1) 光学波频率的最大值
的最大值
;
和最小值 ,声学波频率
2) 相应声子的能量 , 和
;
3) 在
下,三种声子数目各为多少?
4) 如果用电磁波激发光学波
要激发的声子所用的电磁波波长在什么波段?
F U0 kBT
q
ln( q )
kBT
解:已知晶体自由能可以表示为
F U0
q
1 2
q
kBT
ln(1
q
e
kBT
)
量子谐振子的自由能为:
F U kBT
q
1 2
q
kBT
ln(1 eq kBT )
对于经典极限,kBT q ,因而有
1 q 0
2 kBT
eq kBT 1 q
kBT
4
F fi c ( 1,m ,m ) ( ,m 1,m ) i 1
c ( ,m1 ,m ) ( ,m ,m1)
c[( 1,m 1,m 2 ,m )
(,m1 ,m1 2,m )]
根据运动方程形式:
M, m
M
d
2l,
dt 2
m
F
c[( 1,m 1,m 2 ,m )
两边微分
将dq和
代入
得到 时
为虚数,有
解方法 2:
振动模式密度函数 已知三维色散关系
q空间的等频率面是球面,q为常数
对于光学波,在
处振动频率具有最大值
频率分布函数
3.8 有N各相同原子组成面积为S的二维晶格,在德拜 近似下计算比热,并论述在低温极限比热正比于T2。
解:德拜模型考虑的格波是弹性波,波速为v的格波的色散关
ei(t naq) n
ei(t naq) n
代入到运动方程,得:
2Mn (10 n neiaq 11n )