[r]
(1)ClearAT, h, m, kb, chivc, chi0, epxi, Ωl , qd, Lx, Ly , Lz, dt, E0, Ef, kΑ, kΒ, EnΑ, EnΒ, EnΗ, nxΑ, nyΑ, nxΒ, nyΒ, nxΗ, nyΗ, i, j, l, n, EΑ, EΒ, fΑ, fΒ, T, ND, HS, Nq, Ω, q1, q2, q3, q4, d1, d2, d3, d4 , H1, H2, H3, H4, B, P, p1, p2E; T= 200;
h=1.054 10-34; m=6.097 10-32; kb=1.38066 10-23; chivc=10.9; chi0=13.1; epxi=12.5;
Ωl =5.556 1012; qd= 108;
Ly = 200 10-10; Lz = 100 000 10-10; dt=1.6 10-19; E0=10-5; Ef=0.5 10-18; kΑ =109.25; kΒ =109.3556;
nxΗ =1;
nyΗ =1;
nxΑ :=i;
nyΑ := j;
nxΒ := n; nyΒ :=l; EnΑ@Lx_D= h
2Pi2
2 m
nxΑ2
Lx2
+nyΑ 2
(2)EnΗ@Lx_D=
h2Pi2 2 m
nxΗ2
Lx2
+nyΗ 2
Ly2 ;
EnΒ@Lx_D=
h2Pi2
2 m
nxΒ2
Lx2
+ nyΒ 2
Ly2 ;
EΑ@Lx_D= h 2kΑ2
2 m
+EnΑ@LxD;
EΒ@Lx_D= h 2kΒ2
2 m
+EnΒ@LxD;
fΑ@Lx_D= 1 ã
HEΑ@LxD-EfL
kb T +1
;
fΒ@Lx_D= 1 ã
HEΒ@LxD-EfL
kb T +1
;
ND@Lx_D=fΑ@LxD-fΒ@LxD;
HS@Lx_D= 12
2dt2 Ωl
ND@LxDepxiHLx LyL2Lz Pi2
1 chivc
- 1
chi0 ;
Nq= 1
ã
hΩl kb T -1
;
q1@Lx_, Ω_D= kΑ2 + 2 m
h2 H
hΩ -hΩl +EnΗ@LxD-EnΑ@LxDL ;
q2@Lx_, Ω_D= kΑ2 + 2 m
h2 H
hΩ +hΩl -EnΗ@LxD+EnΑ@LxDL ;
q3@Lx_, Ω_D= kΒ2 - 2 m
h2 H
hΩ -hΩl +EnΗ@LxD-EnΒ@LxDL ;
q4@Lx_, Ω_D= kΒ2 - 2 m
h2 H
(3)d1@Lx_,Ω_D=
4 Pi2nxΑ2
Lx2 +
4 Pi2nyΑ2
Ly2 +q1@Lx, ΩD
4 Pi2nxΑ2
Lx2 +
4 Pi2nyΑ2
Ly2 +q1@Lx, ΩD
2 +qd2
+ 2 q1@Lx, ΩD
2
Iq1@Lx, ΩD2+qd2M2
;
d2@Lx_,Ω_D=
4 Pi2nxΑ2
Lx2 +
4 Pi2nyΑ2
Ly2 +q2@Lx, ΩD
4 Pi2nxΑ2
Lx2 +
4 Pi2nyΑ2
Ly2 +q2@Lx, ΩD
2 +qd2
+ 2 q2@Lx, ΩD
2
Iq2@Lx, ΩD2+qd2M2
;
d3@Lx_,Ω_D=
4 Pi2 nxΒ2
Lx2 +
4 Pi2nyΒ2
Ly2 +q3@Lx, ΩD
4 Pi2 nxΒ2
Lx2 +
4 Pi2nyΒ2
Ly2 +q3@Lx, ΩD
2+qd2
+ 2 q3@Lx, ΩD
2
Iq3@Lx, ΩD2 +qd2M2
;
d4@Lx_,Ω_D=
4 Pi2 nxΒ2
Lx2 +
4 Pi2nyΒ2
Ly2 +q4@Lx, ΩD
4 Pi2 nxΒ2
Lx2 +
4 Pi2nyΒ2
Ly2 +q4@Lx, ΩD
2+qd2
+ 2 q4@Lx, ΩD
2
Iq4@Lx, ΩD2+qd2M2 ;
H1@Lx_, Ω_D=fΑ@LxD 1- 1 ã
h2 q1@Lx,ΩD2
2 m +EnΗ@LxD-Ef
kb T +1
+Nq -Nq
1
ã
h2 q1@Lx,ΩD2
2 m +EnΗ@LxD-Ef
kb T +1
;
H2@Lx_, Ω_D= 1- 1 ã
h2 q2@Lx,ΩD2
2 m +EnΗ@LxD-Ef
kb T +1
H1-fΑ@LxD+NqL-Nq fΑ@LxD;
H3@Lx_, Ω_D= 1
ã
h2 q3@Lx,ΩD2
2 m +EnΗ@LxD-Ef
kb T +1
(4)H4@Lx_, Ω_D=fΒ@LxD 1- 1 ã
h2 q4@Lx,ΩD2
2 m +EnΗ@LxD-Ef
kb T +1
+Nq -Nq
1
ã
h2 q4@Lx,ΩD2
2 m +EnΗ@LxD-Ef
kb T +1
;
B@Lx_, Ω_D= HS@LxD Hd1@Lx,ΩDH1@Lx,ΩD
-d2@Lx,ΩDH2@Lx,ΩD+d3@Lx,ΩDH3@Lx,ΩD-d4@Lx,ΩDH4@Lx,ΩDL; P@Lx_, Ω_D= E0
2
2Ω â i=1
â
j=1
â
n=1
â
l=1
2 dt2h kΑkΒ
m2Lx Ly ND@LxDB@Lx, ΩD
HhΩ -HEΒ@LxD-EΑ@LxDLL2+B@Lx, ΩD2
;
p1=PlotA9PA162 10-10, Ω1015 E1012, PA160 10-10, Ω1015E1012, PA158 10-10, Ω1015E1012, PA156 10-10, Ω1015E1012,
PA154 10-10, Ω1015E1012, PA152 10-10, Ω1015E1012, PA150 10-10, Ω1015E1012, PA148 10-10, Ω1015E1012, PA146 10-10, Ω1015E1012, PA144 10-10, Ω1015E1012, PA142 10-10, Ω1015E1012, PA140 10-10, Ω1015E1012,
PA138 10-10, Ω1015E1012, PA136 10-10, Ω1015E1012=, 8Ω, 1.7 , 2.2<E p2=PlotA9PA162 10-10, Ω1015 E10128.10, PA160 10-10, Ω1015E1012 6.90,
PA158 10-10, Ω1015E10126.00, PA156 10-10, Ω1015 E1012 5.00, PA154 10-10, Ω1015E10124.30, PA152 10-10, Ω1015 E1012 3.70, PA150 10-10, Ω1015E10123.20, PA148 10-10, Ω1015 E1012 2.70, PA146 10-10, Ω1015E10122.30, PA144 10-10, Ω1015 E1012 1.95, PA142 10-10, Ω1015E10121.66, PA140 10-10, Ω1015 E1012 1.40,
PA138 10-10, Ω1015E10121.19, PA136 10-10, Ω1015 E1012=, 8Ω, 1.7 , 2.2<E
(5)FindRootAPA160 10 ,Ω10 E10 6.90 3.2,8Ω, 1.876<E; FindRootAPA160 10-10,Ω1015E10126.903.2,8Ω, 2.020<E; FindRootAPA158 10-10,Ω1015E10126.00 3.2,8Ω, 1.869<E; FindRootAPA158 10-10,Ω1015E10126.003.2,8Ω, 2.022<E; FindRootAPA156 10-10,Ω1015E10125.00 3.2,8Ω, 1.863<E; FindRootAPA156 10-10,Ω1015E10125.003.2,8Ω, 2.023<E; FindRootAPA154 10-10,Ω1015
E10124.30 3.2,8Ω, 1.860<E; FindRootAPA154 10-10,Ω1015
E10124.303.2,8Ω, 2.024<E; FindRootAPA152 10-10,Ω1015
E10123.70 3.2,8Ω, 1.852<E; FindRootAPA152 10-10,Ω1015
E10123.703.2,8Ω, 2.031<E; FindRootAPA150 10-10,Ω1015
E10123.20 3.2,8Ω, 1.842<E; FindRootAPA150 10-10,Ω1015
E10123.203.2,8Ω, 2.037<E; FindRootAPA148 10-10,Ω1015
E10122.70 3.2,8Ω, 1.834<E; FindRootAPA148 10-10,Ω1015
E10122.703.2,8Ω, 2.039<E; FindRootAPA146 10-10,Ω1015
E10122.30 3.2,8Ω, 1.822<E; FindRootAPA146 10-10,Ω1015
E10122.303.2,8Ω, 2.046<E; FindRootAPA144 10-10,Ω1015
E10121.95 3.2,8Ω, 1.812<E; FindRootAPA144 10-10,Ω1015
E10121.953.2,8Ω, 2.050<E; FindRootAPA142 10-10,Ω1015
E10121.66 3.2,8Ω, 1.798<E; FindRootAPA142 10-10,Ω1015
E10121.663.2,8Ω, 2.057<E; FindRootAPA140 10-10,Ω1015
E10121.40 3.2,8Ω, 1.782<E; FindRootAPA140 10-10,Ω1015
E10121.403.2,8Ω, 2.062<E; xData=9140 10-10, 142 10-10, 144 10-10, 146 10-10, 148 10-10,
150 10-10, 152 10-10, 154 10-10, 156 10-10, 158 10-10, 160 10-10, 162 10-10 =
yData=90.37706 1015, 0.33666 1015, 0.30678 1015, 0.28405 1015, 0.26399 1015, 0.24743 1015, 0.23128 1015, 0.21698 1015, 0.20471 1015, 0.19420 1015, 0.18703 1015, 0.181264 1015= dataToPlot =TransposeB:xData 1010,
h dt
yData>F
H*ListPlot@dataToPlot ,PlotStyle -> PointSize@0.03D, AxesLabel®8"TEMPERATURE HKL ","LINEWIDTHHeVL"<D*L
ListPlotAdataToPlot , Frame®True, FrameLabel®8"Lx HmL ", "Half-Linewidth HeVL"<, PlotStyle®8PointSize@0.03D<, PlotRange®99136 10-10, 162 10-10=,80, 0.4<=E ListPlotAdataToPlot , Frame®True, FrameTicks®88True, None<,8True, None<<,
FrameLabel ®8"Kich thuoc cua day HAL", "Nua rong pho HeVL"<, PlotStyle®8PointSize@0.03D, Red<, PlotRange ®Automatic
(6)PlotStyle®8PointSize@0.03D, Red<, PlotRange ®Automatic
H*99136 10-10,162 10-10=,80,0.4<=*L, Mesh®Full, InterpolationOrder ®2, Joined®TrueE
1.8 1.9 2.0 2.1 2.2
2
1.8 1.9 2.0 2.1 2.2
2
:
500 000 000 ,
71 000 000 000
, 625 000 000
,
73 000 000 000
,
37 500 000 000
, 200 000 000
,
19 250 000 000
,
77 000 000 000
,
39 500 000 000
,
79 000 000 000
, 62 500 000
,
81 000 000 000>
(7)88140, 0.248388<,8142, 0.221775<,8144, 0.202091<,8146, 0.187118<,
8148, 0.173903<,8150, 0.162995<,8152, 0.152356<,8154, 0.142936<,
8156, 0.134853<,8158, 0.127929<,8160, 0.123206<,8162, 0.119408<<
140 145 150 155 160 0.12
0.14 0.16 0.18 0.20 0.22 0.24
Kich thuoc cua dayHAL
Nua
do
rong
pho
H
eV