Ả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 Các chƣơng trình chạy 2.1 Sự phụ thuộc hệ số hấp thụ vào nhiệt độ % CAC CHUONG TRINH CHAY % 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 hệ số hấp thụ vào cƣờng độ sóng điện từ mạnh E01 % CAC CHUONG TRINH CHAY % 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 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 hệ số hấp thụ vào lƣợng sóng điện từ mạnh (h*Omeg1) % CAC CHUONG TRINH CHAY % 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 hệ số hấp thụ vào lƣợng sóng điện từ yếu (h*Omeg2) % CAC CHUONG TRINH CHAY % 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 % 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 ... 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 hệ số hấp thụ vào cƣờng độ sóng điện từ mạnh E01 % CAC... thi alpha-E01'); 70 xlabel(' Bien 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 hệ số hấp thụ vào lƣợng sóng điện từ mạnh. .. axis([0.0125 0.0301 -0.5e-8 2e-8]); 2.4 Sự phụ thuộc hệ số hấp thụ vào lƣợng sóng điện từ yếu (h*Omeg2) % CAC CHUONG TRINH CHAY % Alpha phu thuoc vao h*Omeg2: clc;clear all;close all; T1=50; T2=150;