Ảnh hưởng của sóng điện từ mạnh lên hấp thụ sóng điện từ yếu bởi điện tử giam cầm trong siêu mạng pha tạp có kể đến hiệu ứng giam cầm của phonon (trường hợp tán xạ điện tử phonon âm)

e0=1;m0=9.1093892-31;h=1.05459*10^(-34);
m1=.067*m0;e=2.07*e0;n0=1e23;
k=1.3807e-23;
a11=2*pi^2*e^2*k.*T*V^2./((c*sqrt(Xinf)*r0*vs^2*m1^2.*omeg2.^3).*h.^2);
y=0;
for n=1:3
for t=1:mt
for n1=1:3
if (Imnn5(t,n,n1)~=0)&(Imnn5(t,n,n1)~=inf);
al=a11.*(D(0,1,omeg1,omeg2,T,d,E02,n,n1)...
-D(0,-2,omeg1,omeg2,T,d,E02,n,n1)...
-1/2.*(H(0,1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
-H(0,-1,omeg1,omeg2,T,d,E01,E02,gama,n,n1))...
+3/32.*(G(0,1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
-G(0,-1,omeg1,omeg2,T,d,E01,E02,gama,n,n1))...
+1/4.*(H(-1,1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
-H(-1,-1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
+H(1,1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
-H(1,-1,omeg1,omeg2,T,d,E01,E02,gama,n,n1))...
-1/16.*(G(-1,1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
-G(-1,-1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
+G(1,1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
-G(1,-1,omeg1,omeg2,T,d,E01,E02,gama,n,n1))...

68


+1/64.*(G(-2,1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
-G(-2,-1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
+G(2,1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...

-G(2,-1,omeg1,omeg2,T,d,E01,E02,gama,n,n1))).*Imnn5(t,n,n1);
y=y+al;
end
end
end
end
end
2. Các chƣơng trình chạy
2.1. Sự phụ thuộc của hệ số hấp thụ vào nhiệt độ
% CAC CHUONG TRINH CHAY 1
% Bien doi theo T
clc;clear all;close all;
T=linspace(100,300,20);
E011=3.6e9;
E012=3.8e9;
E013=4e9;
E02=8e6;
d=85e-10;omeg1=2.1e12;
omeg2=13.8e12;gama=0;k=1.3807e-23;
V=13.5*1.6021e-19;c=1;Xinf=10.9;
r0=5320;vs=5370;Xif=10.9;
e0=1;m0=9.109389e-31;h=1;
m1=0.067*m0;e=2.07*e0;n0=1e23;
y11=y(omeg1,omeg2,T,d,E011,E02,gama);
y12=y(omeg1,omeg2,T,d,E012,E02,gama);
y13=y(omeg1,omeg2,T,d,E013,E02,gama);
figure(1)
plot(T,y11,'r');grid on;hold on;

69

plot(T,y12,'b');grid on;hold on;
plot(T,y13,'g');grid on;hold on;
title(' Do thi \alpha-T');
xlabel('Nhiet do(K)');
ylabel('he so hap thu \alpha');
legend('E011=3.6e9','E012=3.8e9','E013=4e9');
2.2. Sự phụ thuộc của hệ số hấp thụ vào cƣờng độ sóng điện từ mạnh E01
% CAC CHUONG TRINH CHAY 2
% Alpha phu thuoc vao E01:
clc;clear all;close all;
E01=linspace(1e9,8e9,80);
E02=8e6;
T1=150;
T2=170;
T3=200;
T4=250;
d=85e-10;omeg1=2.1e12;
omeg2=13.8e12;gama=0;k=1.3807e-23;
V=13.5*1.6021e-19;c=1;Xinf=10.9;
r0=5320;vs=5370;Xif=10.9;
e0=1;m0=9.109389e-31;h=1;
m1=0.067*m0;e=2.07*e0;n0=1e23;
y1=y(omeg1,omeg2,T1,d,E01,E02,gama);
plot(E01,y1,'.b');grid on;hold on;
y2=y(omeg1,omeg2,T2,d,E01,E02,gama);
plot(E01,y2,'-r');grid on;hold on;
y3=y(omeg1,omeg2,T3,d,E01,E02,gama);
plot(E01,y3,'.g');grid on;hold on;

y4=y(omeg1,omeg2,T4,d,E01,E02,gama);
plot(E01,y4,'.g');grid on;hold on;
title('Do thi \alpha-E01');

70


xlabel('

Bien do song E01(V/m)');

ylabel('He so hap thu \alpha');
legend('T1=150K','T2=170K','T3=200K', 'T4=250K');
2.3. Sự phụ thuộc của hệ số hấp thụ vào năng lƣợng sóng điện từ mạnh
(h*Omeg1)
% CAC CHUONG TRINH CHAY 3
% Alpha phu thuoc vao h*Omeg1 (truong hop phi tuyen):
clc;clear all;close all;
T1=50;
T2=150;
T3=200;
E01=8.1e7;
E02=8e5
omeg1=linspace(20e12,45e12,100);
N1=omeg1;
omeg2=2e11;
d=85e-10;gama=10;
k=1.3807e-23;V=13.5*1.6021e-19;c=1;Xinf=10.9;
r0=5320;vs=5370;
e0=1.6e-19;m0=9.109389e-31;h=1;

m1=0.067*m0;e=2.07*e0;n0=1e20;
N1=omeg1*1.0544e-34/1.6*1e19;
y1=y(omeg1,omeg2,T1,d,E01,E02,gama);
y2=y(omeg1,omeg2,T2,d,E01,E02,gama);
y3=y(omeg1,omeg2,T3,d,E01,E02,gama);
plot(N1,y1,'.b');grid on;hold on;
plot(N1,y2,'-r');grid on;hold on;
plot(N1,y3,'-g');grid on;hold on;
title('Do thi \alpha-h*omeg1');
xlabel('

Nang luong song h*omeg1(eV)');

ylabel('he so hap thu \alpha');

71


legend('T1=50K','T2=150K','T3=200K');
axis([0.0125 0.0301 -0.5e-8 2e-8]);
2.4. Sự phụ thuộc của hệ số hấp thụ vào năng lƣợng sóng điện từ yếu
(h*Omeg2)
% CAC CHUONG TRINH CHAY 4
% Alpha phu thuoc vao h*Omeg2:
clc;clear all;close all;
T1=50;
T2=150;
T3=250;
E01=8.1e9;
E02=8e7;

omeg2=linspace(30e12,90e12,150);
d=850e-10;omeg1=20e12;gama=0;k=1.3807e-23;
V=13.5*1.6021e-19;c=1;Xinf=10.9;
r0=5320;vs=5370;Xif=10.9;
e0=1;m0=9.109389e-31;h=1;
m1=0.067*m0;e=2.07*e0;n0=1e23;
N1=omeg2*1.0544e-34/1.6*1e19;
y1=y(omeg1,omeg2,T1,d,E01,E02,gama);
y2=y(omeg1,omeg2,T2,d,E01,E02,gama);
y3=y(omeg1,omeg2,T3,d,E01,E02,gama);
plot(N1,y1,'.b');grid on;hold on;
plot(N1,y2,'.r');grid on;hold on;
plot(N1,y3,'.g');grid on;hold on;
title('Do thi \alpha-h*omeg2');
xlabel('

Nang luong song h*omeg2(eV)');

ylabel('he so hap thu \alpha');
legend('T1=50K','T2=150K','T2=250K');
axis([0.019 0.06 -5e-7 0.6e-7]);
2.5. Sự phụ thuộc hệ số hấp thụ vào chu kỳ siêu mạng

72


% CAC CHUONG TRINH CHAY 5
% Alpha phu thuoc vao chu ky sieu mang pha tap d:
clc;clear all;close all;
T1=120;

T2=150;
T3=200;
omeg2=2e12;
omeg1=13.8e12;gama=0;k=1.3807e-23;V=13.5*1.6021e-19;c=1;Xinf=10.9;
r0=5320;vs=5370;Xif=10.9;
e0=1;m0=9.109389e-31;h=1;
m1=0.067*m0;e=2.07*e0;n0=1e23;
al=0;
E01=1e6;
E02=3e4;
d=linspace(50e-10,999.1e-10,200);
y1=y(omeg1,omeg2,T1,d,E01,E02,gama);
plot(d,y1,'.r');grid on;hold on;
y2=y(omeg1,omeg2,T2,d,E01,E02,gama);
plot(d,y2,'.b');grid on;hold on;
y3=y(omeg1,omeg2,T3,d,E01,E02,gama);
plot(d,y3,'.g');grid on;hold on;
title('Do thi \alpha-d');
xlabel('

Chu ky sieu mang d');

ylabel('he so hap thu \alpha');
legend('T1=120K','T2=150K', 'T2=200K');

73