level 3
forest353
楼主
(*Parameters, J. Bonse, JAP 2005*)
ClearAll["Global`*"];
s = 10;(*Shape factor, which could be any value above 0 when \
\[Theta]=0, referred to J. E. Sipe - PRB 1983*)
f = 0.1;(*Filling factor, which could be any value above 0 when \
\[Theta]=0, referred to J. E. Sipe - PRB 1983*)
(*\[Epsilon]=(3.69+I*0.0065)^2;(*Dielectric permittivity, unexcited \
state of Si*)*)
\[Epsilon] = -110 +
I*8;(*Dielectric permittivity, obtained from Drude model*)
\[Theta] = 0;(*Incident angle*)
(*Equations referred to J. Bonse - JAP 2005, which is concluded by \
Bonse from J. E. Sipe - PRB 1983*)
rc = (\[Epsilon] - 1)/(\[Epsilon] + 1);
fc = Sqrt[(s)^2 + 1] - s;
gc = (Sqrt[(s)^2 + 4] + s)/2 - Sqrt[(s)^2 + 1];
\[Eta] = 2*Pi*Abs[\[Nu]1 + Conjugate[\[Nu]2]];
(*s-polarization*)
\[Nu]1 = (hss1*(k1y)^2 + hkk1*(k1x)^2)*\[Gamma]t*(Abs[ts])^2;
\[Nu]2 = (hss2*(k2y)^2 + hkk2*(k2x)^2)*\[Gamma]t*(Abs[ts])^2;
(k1y) = (Sin[\[Theta]] + \[Kappa]y)/\[Kappa]1;
(k2y) = (Sin[\[Theta]] - \[Kappa]y)/\[Kappa]2;
(k1x) = \[Kappa]x/\[Kappa]1;
(k2x) = \[Kappa]x/\[Kappa]2;
\[Kappa]1 = Sqrt[(\[Kappa]x)^2 + (Sin[\[Theta]] + \[Kappa]y)^2];
\[Kappa]2 = Sqrt[(\[Kappa]x)^2 + (Sin[\[Theta]] - \[Kappa]y)^2];
hss1 = 2*I/(Sqrt[1 - (\[Kappa]1)^2] +
Sqrt[\[Epsilon] - (\[Kappa]1)^2]);
hss2 = 2*I/(Sqrt[1 - (\[Kappa]2)^2] +
Sqrt[\[Epsilon] - (\[Kappa]2)^2]);
hkk1 = 2*I*
Sqrt[(\[Epsilon] - (\[Kappa]1)^2)*(1 - (\[Kappa]1)^2)]/(\[Epsilon]*
Sqrt[1 - (\[Kappa]1)^2] + Sqrt[\[Epsilon] - (\[Kappa]1)^2]);
hkk2 = 2*I*
Sqrt[(\[Epsilon] - (\[Kappa]2)^2)*(1 - (\[Kappa]2)^2)]/(\[Epsilon]*
Sqrt[1 - (\[Kappa]2)^2] + Sqrt[\[Epsilon] - (\[Kappa]2)^2]);
ts = 2*Abs[
Cos[\[Theta]]]/(Abs[Cos[\[Theta]]] +
Sqrt[\[Epsilon] - (Sin[\[Theta]])^2]);
\[Gamma]t = (\[Epsilon] - 1)/(4*
Pi*(1 + (1 - f)*(\[Epsilon] - 1)*(fc - rc*gc)/2));
(*Equations finished*)
(*Plotting*)
DensityPlot[\[Eta], {\[Kappa]x, -4, 4}, {\[Kappa]y, -4, 4},
PlotPoints -> 100]
Export["03.xls",
Table[\[Eta], {\[Kappa]x, -4, 4, 0.1}, {\[Kappa]y, -4, 4, 0.1}]]
(*
Plot3D[\[Eta],{\[Kappa]x,-4,4},{\[Kappa]y,-4,4},PlotPoints\[Rule]100]
*)
2017年09月11日 14点09分
1
ClearAll["Global`*"];
s = 10;(*Shape factor, which could be any value above 0 when \
\[Theta]=0, referred to J. E. Sipe - PRB 1983*)
f = 0.1;(*Filling factor, which could be any value above 0 when \
\[Theta]=0, referred to J. E. Sipe - PRB 1983*)
(*\[Epsilon]=(3.69+I*0.0065)^2;(*Dielectric permittivity, unexcited \
state of Si*)*)
\[Epsilon] = -110 +
I*8;(*Dielectric permittivity, obtained from Drude model*)
\[Theta] = 0;(*Incident angle*)
(*Equations referred to J. Bonse - JAP 2005, which is concluded by \
Bonse from J. E. Sipe - PRB 1983*)
rc = (\[Epsilon] - 1)/(\[Epsilon] + 1);
fc = Sqrt[(s)^2 + 1] - s;
gc = (Sqrt[(s)^2 + 4] + s)/2 - Sqrt[(s)^2 + 1];
\[Eta] = 2*Pi*Abs[\[Nu]1 + Conjugate[\[Nu]2]];
(*s-polarization*)
\[Nu]1 = (hss1*(k1y)^2 + hkk1*(k1x)^2)*\[Gamma]t*(Abs[ts])^2;
\[Nu]2 = (hss2*(k2y)^2 + hkk2*(k2x)^2)*\[Gamma]t*(Abs[ts])^2;
(k1y) = (Sin[\[Theta]] + \[Kappa]y)/\[Kappa]1;
(k2y) = (Sin[\[Theta]] - \[Kappa]y)/\[Kappa]2;
(k1x) = \[Kappa]x/\[Kappa]1;
(k2x) = \[Kappa]x/\[Kappa]2;
\[Kappa]1 = Sqrt[(\[Kappa]x)^2 + (Sin[\[Theta]] + \[Kappa]y)^2];
\[Kappa]2 = Sqrt[(\[Kappa]x)^2 + (Sin[\[Theta]] - \[Kappa]y)^2];
hss1 = 2*I/(Sqrt[1 - (\[Kappa]1)^2] +
Sqrt[\[Epsilon] - (\[Kappa]1)^2]);
hss2 = 2*I/(Sqrt[1 - (\[Kappa]2)^2] +
Sqrt[\[Epsilon] - (\[Kappa]2)^2]);
hkk1 = 2*I*
Sqrt[(\[Epsilon] - (\[Kappa]1)^2)*(1 - (\[Kappa]1)^2)]/(\[Epsilon]*
Sqrt[1 - (\[Kappa]1)^2] + Sqrt[\[Epsilon] - (\[Kappa]1)^2]);
hkk2 = 2*I*
Sqrt[(\[Epsilon] - (\[Kappa]2)^2)*(1 - (\[Kappa]2)^2)]/(\[Epsilon]*
Sqrt[1 - (\[Kappa]2)^2] + Sqrt[\[Epsilon] - (\[Kappa]2)^2]);
ts = 2*Abs[
Cos[\[Theta]]]/(Abs[Cos[\[Theta]]] +
Sqrt[\[Epsilon] - (Sin[\[Theta]])^2]);
\[Gamma]t = (\[Epsilon] - 1)/(4*
Pi*(1 + (1 - f)*(\[Epsilon] - 1)*(fc - rc*gc)/2));
(*Equations finished*)
(*Plotting*)
DensityPlot[\[Eta], {\[Kappa]x, -4, 4}, {\[Kappa]y, -4, 4},
PlotPoints -> 100]
Export["03.xls",
Table[\[Eta], {\[Kappa]x, -4, 4, 0.1}, {\[Kappa]y, -4, 4, 0.1}]]
(*
Plot3D[\[Eta],{\[Kappa]x,-4,4},{\[Kappa]y,-4,4},PlotPoints\[Rule]100]
*)

