固体物理学第三章资料
合集下载
固体物理-第三章 金属自由电子论讲解

N=I0G(EF)+ I1G’(EF)+ I2G’’(EF)+….. 其中, I0=- (-f/E) dE, I1=-(E-EF)(-f/E)dE,
3.1.量子自由电子理论
I2=(1/2!)-(E-EF)2(-f/E) dE 不难算出, I0=1(d-函数积分), I1=0 (根据d-函数的性质) 为了计算I2, 而令h=(E-EF)/kBT,于是, I2=[(kBT)2/2]-{h2/[(eh+1)(e-h+1)] }dh=(pkBT)2/6
波长),可见k为电子的波矢, 是3 维空间矢量. r:电 子的位置矢量。
由波函数的归一化性质:vy*(r) y(r)d(r)=1, v:金属体积, 假设为立方体,边长为L,把3.1.1.3式 代入归一化式子, 得: A=L-3/2=V-1/2, 所以
y(r)= V-1/2eik•r 3.1.1.4, 此即自由电子的本征态。 由周期性边界条件, y(x,y,z)= y(x+L,y,z) = y(x,y+L,z) = y(x,y,z+L)
一状态的电子具有确定的动量ħk和能量ħ2k2/(2m),因而 具有确定的速度,v=ħk/m,故一个k全面反映了自由电子 的一个状态,简称态。
2. k-空间
以kx, ky , kz 为坐标轴建立的 波矢空间叫k-空间。电子的 本征态可以用该空间的一点
来代表。点的坐标由3.1.1.5 式确定。
3.1.量子自由电子理论
T>0K的费米能EF 把3.1.2.2和3.1.3.1代入3.1.3.2, 分步积分, 得:
N= (-2C/3) 0 E3/2(f/E) dE 3.1.3.3 令G(E)= 2C E3/2/3, 3.1.3.3.式化简为 N= 0G(E) (-f/E) dE 3.1.3.4 (-f/E)函数具有类似d函数的特性,仅仅在EF附近kBT范 围内才有显著的值,且为E-EF偶函数. 由于(-f/E)函数 具有这些性质,把G(E)在EF附近展开为泰勒级数, 且积分 下限写成 -,不会影响积分值. 3.1.3.4化为:
3.1.量子自由电子理论
I2=(1/2!)-(E-EF)2(-f/E) dE 不难算出, I0=1(d-函数积分), I1=0 (根据d-函数的性质) 为了计算I2, 而令h=(E-EF)/kBT,于是, I2=[(kBT)2/2]-{h2/[(eh+1)(e-h+1)] }dh=(pkBT)2/6
波长),可见k为电子的波矢, 是3 维空间矢量. r:电 子的位置矢量。
由波函数的归一化性质:vy*(r) y(r)d(r)=1, v:金属体积, 假设为立方体,边长为L,把3.1.1.3式 代入归一化式子, 得: A=L-3/2=V-1/2, 所以
y(r)= V-1/2eik•r 3.1.1.4, 此即自由电子的本征态。 由周期性边界条件, y(x,y,z)= y(x+L,y,z) = y(x,y+L,z) = y(x,y,z+L)
一状态的电子具有确定的动量ħk和能量ħ2k2/(2m),因而 具有确定的速度,v=ħk/m,故一个k全面反映了自由电子 的一个状态,简称态。
2. k-空间
以kx, ky , kz 为坐标轴建立的 波矢空间叫k-空间。电子的 本征态可以用该空间的一点
来代表。点的坐标由3.1.1.5 式确定。
3.1.量子自由电子理论
T>0K的费米能EF 把3.1.2.2和3.1.3.1代入3.1.3.2, 分步积分, 得:
N= (-2C/3) 0 E3/2(f/E) dE 3.1.3.3 令G(E)= 2C E3/2/3, 3.1.3.3.式化简为 N= 0G(E) (-f/E) dE 3.1.3.4 (-f/E)函数具有类似d函数的特性,仅仅在EF附近kBT范 围内才有显著的值,且为E-EF偶函数. 由于(-f/E)函数 具有这些性质,把G(E)在EF附近展开为泰勒级数, 且积分 下限写成 -,不会影响积分值. 3.1.3.4化为:
《固体物理·黄昆》第三章

氢键结合的情况可写成通式:
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等在低温下依靠范德瓦耳斯力结合成的晶体.
固体物理学03_04

K
K K i [ ωt − R l ⋅ q ] k
()
中,波矢的作用体现在不同原胞之间的位相联系: e
K K − iR ( l )⋅q
REVISED TIME: 05-4-9
-2-
CREATED BY XCH
固体物理学_黄昆_第三章 晶格振动与晶体的热学性质_20050406
波矢改变一个倒格矢: Gn = n1b1 + n2b2 + n3b3
V* =
K K K v0 * , v0 * = b1 ⋅ (b2 × b3 ) —— 倒格子原胞体积 N
N N Nv0 V = K K K = = 3 v0 * b1 ⋅ (b2 × b3 ) (2π ) (2π )3
状态密度:
波矢 q 的取值
K
在原子振动波函数 µα ⎜ ⎟ = Ak e
⎛l ⎞ ⎝k ⎠
⎧ b1 h1b1 b1 N ⎧ N1 − < h1 ≤ 1 ⎪− 2 < N ≤ 2 ⎪ 1 2 2 ⎪ ⎪ K K K ⎪ h h h b hb b N K ⎪ N q = 1 b1 + 2 b2 + 3 b3 —— ⎨ − 2 < 2 2 ≤ 2 Ì ⎨ − 2 < h2 ≤ 2 N1 N2 N3 N2 2 2 ⎪ 2 ⎪ 2 N ⎪ N3 ⎪ b3 h3b3 b3 − < h3 ≤ 3 ≤ ⎪ ⎪− < 2 ⎩ 2 N3 2 ⎩ 2
—— q 和 q + Gn 产生的位相一样
K
K
K
b1 ⎧ b1 ⎪ − 2 < qx ≤ 2 ⎪ b K ⎪ b —— q 的取值限制在一个倒格子原胞中,即第一布里渊区: ⎨ − 2 < q y ≤ 2 2 ⎪ 2 b ⎪ b3 − < qz ≤ 3 ⎪ 2 ⎩ 2
固体物理学:第3章 晶格振动

2 2
21 2
cos
qa
1 2
光学支
2 o
1
m
2 1 m
1
2 1
2 2
21
2
cos
qa
2
声学支
2A
1
m
2 1 m
12 22 21 2 cos qa
1 2
三、色散关系
UESTC
ω
当 q=0
ωO
ωA = 0 ωo = 21 2
m
ωA
当
q=
a
a
o
q
a
A
21
m
o
2 2
m
四、格波数
q 2 m
Na
2
Na
m 0 , 1, 2
q
o
波矢q 的取值是分立的,相邻q的“距离”N2a
五、格波数
UESTC
此前研究的晶格原子集体的波动运动就是格波。
晶体中所有原子以相同的频率和振幅在 平衡位置附近作简谐振动,原子的运动状 态在晶体中以波的形式传播,这种简谐波 称为格波。
五、格波数
UESTC
3.1 一维单原子链的振动
一. 物理模型 二. 运动方程 三. 色散关系 四. 波恩-卡曼周期性边界条件 五. 格波数 六. 小结
UESTC
一、物理模型
UESTC
一维简单晶格的振动
平衡位置 振动时偏离 平衡位置
un :第n个原子偏离平衡位置的位移 m :原子质量
一、物理模型
UESTC
V (r) V (0) dV (r) r 1 d 2V (r) r2
UESTC
❖ 对于一维原子链,简约区中波数q的取值总
固体物理第三章总结

时以比T3更快的速度趋于零。 温度越低,与实验吻合的越好。
kBE
局限性
E
kB
D
D
kB
晶体的非简谐效应
1.非简谐效应:
U(
R0
)
U(
R0
)
1 2!
2U R2
R0
2
1 3!
3U R3
R0
3
c 2 g 3
im jm
b1
b2
1010 i 1010 j
m 1 m 1
3.14 1010 i m 1 3.14 1010 j m1
a3 21010 km b3 1010 k m1 3.141010 k m1
S
TO 0,
3.极化声子和电磁声子
0
因为长光学波是极化波,且只有长光学纵波才伴随着宏观
的极化电场,所以长光学纵波声子称为极化声子。 长光学横波与电磁场相耦合,它具有电磁性质,称长光学
横波声子为电磁声子。
1.已知模式密度 ( ) 求:
(1)~+d间隔内的振动模式数;
(2) ~+d间隔内的声子数及晶体中总的声子数;
2
2
2
中的 振~ 动模d式数目:2Lc
2 d ,
v
Sc
2
v2
d ,
Vc
2 2
2
v3
d
一维有一支纵波,二维有一支纵波一支横波,三维有
一支纵波两支横波,纵波与横波速度相等
:
Lc 2 d , 2 v
固体物理第三章

19
格波 —— 短波极限情况 ( q →
πa)源自aq ω = 2 β / m sin( ) 2
ωmax = 2 β / m
长波极限下 ( q → 0) ,相邻两个原子之间的位相差
q(n + 1)a − qna = qa ⇒ 0
—— 一个波长内包含许多原子,晶格看作是连续介质 短波极限下 q ⇒
π
a
2π λ= = 2a q
2
17
格波 —— 长波极限情况
4β 2 aq ω = sin ( ) m 2
2
aq ω=2 sin( ) m 2
当 q→0
β
qa qa sin( ) ≈ 2 2
ω = a β /m q
ω =VElasticq
—— 一维单原子格波的色散关系与连续 介质中弹性波的色散关系一致
18
相邻原子之间的作用力 f = βδ 长波极限情况
o xij = x o − xio j
(3.1.2)
u ij = u j − u i
xn −1
•0
un −1
•0
u
n
xn xn
•0
un +1
xn +1
x
4
a
5
设两原子间的相互作用势能为 ϕ ( xij ) ,且只考虑二 体相互作用,则总的相互作用能为
1 N U = ∑ ϕ ( xij ) 2 i≠ j
4β 2 aq ω = sin ( ) m 2
2
相邻原子位相差 aq ⇒ 2π + aq
π
4a 2a 相邻原子位相差 aq1 = π / 2 2π 5π 两种波矢的格波中,原子 两种波矢的格波中, = 格波2(Green)波矢 q2 = 的振动完全相同, 4a / 5 2a 的振动完全相同,相邻原 相邻原子的位相差 aq2 = 2π + π / 2 子的位相差 − π < aq ≤ π
格波 —— 短波极限情况 ( q →
πa)源自aq ω = 2 β / m sin( ) 2
ωmax = 2 β / m
长波极限下 ( q → 0) ,相邻两个原子之间的位相差
q(n + 1)a − qna = qa ⇒ 0
—— 一个波长内包含许多原子,晶格看作是连续介质 短波极限下 q ⇒
π
a
2π λ= = 2a q
2
17
格波 —— 长波极限情况
4β 2 aq ω = sin ( ) m 2
2
aq ω=2 sin( ) m 2
当 q→0
β
qa qa sin( ) ≈ 2 2
ω = a β /m q
ω =VElasticq
—— 一维单原子格波的色散关系与连续 介质中弹性波的色散关系一致
18
相邻原子之间的作用力 f = βδ 长波极限情况
o xij = x o − xio j
(3.1.2)
u ij = u j − u i
xn −1
•0
un −1
•0
u
n
xn xn
•0
un +1
xn +1
x
4
a
5
设两原子间的相互作用势能为 ϕ ( xij ) ,且只考虑二 体相互作用,则总的相互作用能为
1 N U = ∑ ϕ ( xij ) 2 i≠ j
4β 2 aq ω = sin ( ) m 2
2
相邻原子位相差 aq ⇒ 2π + aq
π
4a 2a 相邻原子位相差 aq1 = π / 2 2π 5π 两种波矢的格波中,原子 两种波矢的格波中, = 格波2(Green)波矢 q2 = 的振动完全相同, 4a / 5 2a 的振动完全相同,相邻原 相邻原子的位相差 aq2 = 2π + π / 2 子的位相差 − π < aq ≤ π
固体物理-第三章

l 1
原 子
上式说明每个坐标gk的振动,都可以分解成3N个简正振动的线 性迭加,Ql新坐标称为简正坐标,所以,我们可以得出结论:N个
的
原子组成晶体的任何一种微振动,可看成3N个简正振动的迭加。
运
动
★简正坐标与原子位移坐标之间的正交变换,
实际上是按付氏展开式把坐标系由位置坐标转
换到状态空间(正格子——倒格子)。
单
体原子集体运动状态的激发单元,它不能脱离固体而单
原
独存在,它并不是一种真实的粒子, 只是一种准粒子;
子
➢声子的作用过程遵从能量守恒和准动量守恒。
晶
➢一种格波即一种振动模式称为一种声子,对于由N个原子
格
组成的一维单原子链,有N个格波,即有N种声子,
3.1 晶体中原子的微振动及其量子化
声子
采用“声子”概念不仅表达简洁、处理问题方便(例晶格与微观粒
3N
动
2 Ak bik Ai 0 k 1, 2,L 3N (9) i 1
方程组(9)又可改写成:
3N
bik 2ik Ai 0 k 1, 2,L 3N (10)
i 1
3.1 晶体中原子的微振动及其量子化
原子的运动方程
3N
bik 2ik Ai 0 k 1, 2,L 3N (10)
3.1 晶体中原子的微振动及其量子化
原子的运动方程
••
gk bik gi 0 k 1, 2,L 3N
(7)
原
gk Ak sin t k 1, 2,L 3N
(8)
子
的 运
(8)式所给出的特解应能够满足方程(7),则将(8)式 代入(7)式,得确定ω与bik之间关系的方程组:
固体物理学第三章

因为光学波的频率比声学波的频率高, eO/kBT 1 大于 eA/kB,T所以1在温度一定情况下, 一个光
学波的声子数目少于一个声学波的声子数目.
计算题:试求由5个原子组成的一维单原子晶格的 格波频率,设原子质量m=8.35×10-27kg,恢复力 常数β=15N·m-1
解:一维单原子链的解为 xn Aei(tqn)a
m2
代入,β,m及q值 则得到五个频率依次为(以rad/sec为单位) 8.06×1013,4.99×1013,0,4.99×1013,8.06×1013
3.5 晶格热容
爱因斯坦模型的计算值 比实际值下降得更陡。
考虑了格波的频率分布,低温时只有长光 学波对比热有重要贡献。
德拜截止频率——有限晶体的振动频率不能 无限增大,而是有一上限。
据周期边界条件 ,x1 xN1
此处N=5,代入上式即得:
ei(5a)q 1 5aqn2(n为整数)
由于格波波矢取值范围:
q
a
a
则:5 n 5
22
故n可取-2,-1,0,1,2这五个值
相应波矢:4,2,0,2,4
5a 5a 5a 5a
由于, 2 sin qa
n1,2. ,N
F=ma
一个质点的振动会影 响到其他质点振动
设方程组的解是一振幅为A,角频率为的简谐 振动 :
xnAeitqna
把它代入到运动方程组中,可得振动频率和 波矢q之间的关系式:
1
2
2
sin
q
a
m 2
即与n无关,表明N 个联立方程都归结为同一个方程。
因为:
即qa改变2π的整数倍,原子的振动是一样的。 这样,q的取值范围只需要控制在±π/a之间即可。 这个区间称为第一布里渊区。
学波的声子数目少于一个声学波的声子数目.
计算题:试求由5个原子组成的一维单原子晶格的 格波频率,设原子质量m=8.35×10-27kg,恢复力 常数β=15N·m-1
解:一维单原子链的解为 xn Aei(tqn)a
m2
代入,β,m及q值 则得到五个频率依次为(以rad/sec为单位) 8.06×1013,4.99×1013,0,4.99×1013,8.06×1013
3.5 晶格热容
爱因斯坦模型的计算值 比实际值下降得更陡。
考虑了格波的频率分布,低温时只有长光 学波对比热有重要贡献。
德拜截止频率——有限晶体的振动频率不能 无限增大,而是有一上限。
据周期边界条件 ,x1 xN1
此处N=5,代入上式即得:
ei(5a)q 1 5aqn2(n为整数)
由于格波波矢取值范围:
q
a
a
则:5 n 5
22
故n可取-2,-1,0,1,2这五个值
相应波矢:4,2,0,2,4
5a 5a 5a 5a
由于, 2 sin qa
n1,2. ,N
F=ma
一个质点的振动会影 响到其他质点振动
设方程组的解是一振幅为A,角频率为的简谐 振动 :
xnAeitqna
把它代入到运动方程组中,可得振动频率和 波矢q之间的关系式:
1
2
2
sin
q
a
m 2
即与n无关,表明N 个联立方程都归结为同一个方程。
因为:
即qa改变2π的整数倍,原子的振动是一样的。 这样,q的取值范围只需要控制在±π/a之间即可。 这个区间称为第一布里渊区。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
dk dt
和F的分量相等;
当F与速度v垂直时,可由冲量定理证明在垂直于v
的方向上,
dk dt
和外力F的分量也相等。
F dk dt
这是电子在外场作用下运动状态变化的基本公式,具 有与经典力学中牛顿定律相似的形式。
一、电子的准动量
由经典功能关系 F dp / dt 得:
k —— 电子的准动量
❖ Bloch电子
❖ 有效质量不仅可以取正,也可以取负,在能带底附近 (E(k)极小),有效质量总是正的;而在能带顶附近 ( E(k)极大), 有效质量总是负的。
有效质量的物理解释
1
m m0x 0
0 mx 0
0
0
mx
2
2a 2 J1
0
❖
在能带顶R点:
k
a
,
,
a
a
2
mx my mz m 2a2J1 0
2
2
mx 2E kx2
2a2J1
cos kxa
1
2
my 2E
k2y
2
2a2J1
cos k y a
1
2
2
mz
2E kz2
2a 2 J1
在能带中 的某处,
d 2E dk 2
0
电子速度的数值最大
与自由电子的速度总是随 能量的增加而单调上升是 完全不同的。
❖ 电子速度的方向为k空间中能量梯度的方向,即沿等 能面的法线方向。
电子的运动方向决定于等能面的形状: 在一般情况下,在k空间中,等能面并不是球面,因
此,v的方向一般并不是k的方向;
0
m*>0; m*<0
2. 三维情况
a
dv dt
d dt
1
k
E
1
dk dt
E kk
分量形式:
a
dv dt
d 1
dt
E k
1
3 dk 1 dt
k
E
k
1
2
3
F
1
2E k k
=1, 2, 3,或
x,y,z
原因:在三维情形,沿k空间的不同方向一般有不同的色散关系,
电子的有效质量比较复杂,表现为一个二级张量。
cos kza
1
在能带底和能带顶电子的有效质量是各向同性的,
退化为一标量,这是立方对称的结果。
❖ 在X点:
k
a
,
0,
0
2
mx 2a2J1 0,
2
my mz 2a2J1 0
有效质量是一个很重要的概念,它把晶体中电子准 经典运动的加速度与外力联系起来。
❖ 有效质量中包含了周期场对电子的作用。在一般情况下, 有效质量是一个张量,在特殊情况下可以退化为标量。
2E
kk
2a2J1 cos ka 0
即kx , ky, kz为张量的主轴方向
, 1, 2, 3
2
2
mx
2E kx2
2a 2 J1
cos kxa
1
2
2
my
2E
k
2 y
2a 2 J1
cos kya
1
2
2
mz 2E kz2
2a2J1
cos kza
1
有效质量的三个主分量均与J1成反比,若原子间距越
当等能面为球面时,或沿某些特殊方向,v才与k的
方向相同。
ky v k
kx
3.2 电子在外电场作用下的加速度,有效 质量,等能面
在外场中,电子所受的力为F,在dt时间内,外场
对电子所做的功为F vdt
功能原理: F vdt dE E d k k
v 1 E k
F
dk dt
v
0
在平行于v的方向上,
kx2
0
0
1
mx
0
0
1 m
1
2
0
2E
k
2 y
0 0
1 my
0
0
0
2
E
0
kz2
0
1 mz
dvx dt
1 mx
Fx ,
dvy dt
1 my
Fy ,
dvz dt
1 mz
Fz
例:求简单立方晶体中,紧束缚近似下s带电子的有效质量, 并讨论其特点。
E k s J0 2J1 cos kxa cos kya coskza
a
dv dt
d dt
1
dE dk
பைடு நூலகம்
1
dk dt
d 2E dk 2
2
F
d2E
dk 2
引入电子的有效质量:
F m* dv dt
2
m* d 2E
dk 2
在周期场中电子的有效质量m*与k有关
❖ 在能带底: E(k)取极小值,
❖ 在能带顶: E(k)取极大值,
d 2E dk 2
0
d 2E dk 2
k
r
eikru k
r
的行为类似于波长
为 2 k 的平面波,再由de Broglie关系得其
具有 k 的动量。
❖ 晶体中的电子在碰撞过程中所贡献的动量为 k 。
二、电子的加速度和有效质量
晶体中电子准经典运动的基本关系式:
v 1 E
{
k
dk
F
dt
由以上两式可直接导出在外力作用下电子的加速度。
1. 一维情况
计算结果:
电子的速度:
vk0
dx dt
1
dE dk
k0
三维情况: 电子速度为 v 1 E k
理解:电子的平均速度与其能量和状态有关。
由于E(k)=E(-k), 所以
v(k) v(k)
❖ 电子运动速度的大小与k的关系
以一维为例:
在能带底和能带顶,E(k)取极值, dE 0
dk
在能带底和能带顶,电子速度v=0
2E
kx2
1 m
1
2E
2 kykx
2E
kzkx
2E
k xk y
2E
k
2 y
2E
k z k y
2E
kxkz
2E
kykz
2E
k
2 z
电子的加速度方向并不一定与外力的方向一致。
倒有效质量张量是对称张量,如将kx、ky、kz取为 张量的主轴方向,可将其对角化。
在主轴坐标系中:
2E
大,J1越小,则有效质量就越大。
❖ 在能带底点: k 0, 0, 0
2
mx my mz m 2a2J1 0
有效质量张量退化为一个标量
2
2
mx 2E kx2
2a2J1
cos kxa
1
2
my 2E
k2y
2
2a2J1
cos k y a
1
2
2
mz
2E kz2
2a 2 J1
cos kza
第三章 外场作用下晶体电子的 运动
固体材料
器件
器件工作需要外加电磁场驱动其中的电子
3.1 晶体中电子的速度
处理晶体中电子在外场中的运动所采用的方法:
❖ 解含外场的波动方程
2
2m
2
U
r
V
r
E
❖ 在一定条件下,把晶体中电子在外场中的运动当作 准经典粒子来处理。
条件:外场较弱、恒定,不考虑电子在不同能带间的 跃迁,不涉及电子的衍射和干涉等。
矩阵形式:
2E
ax
ay
az
kx2
1
2E
2 ky kx 2E
kz
kx
2E
k x k y
2E
k
2 y
2E
kz ky
牛顿定律: a 1 F m
2E
k x k z 2E
ky kz 2E
Fx Fy Fz
k
2 z
这里用二阶张量
1 m*
代替了
1 m
倒有效质量张量: