55 lines
1.8 KiB
Plaintext
55 lines
1.8 KiB
Plaintext
|
# Drude-Lorentz models, with the same constants as in qpms.
|
||
|
|
# Compared to qpms, sign of the imaginary part is swapped here due
|
||
|
|
# to different time-frequency transform convention.
|
||
|
|
|
||
|
|
MATERIAL LDSilver
|
||
|
|
# Electronvolt over reduced Planck's constant; the omegas and gammas
|
||
|
|
# are defined in these units.
|
||
|
|
eh = 1519267460583196.5;
|
||
|
|
omegap = 9.01
|
||
|
|
f0 = 0.84
|
||
|
|
f1 = 0.065
|
||
|
|
f2 = 0.124
|
||
|
|
f3 = 0.111
|
||
|
|
f4 = 0.84
|
||
|
|
f5 = 5.646
|
||
|
|
omega0 = 0.0
|
||
|
|
omega1 = 0.816
|
||
|
|
omega2 = 4.481
|
||
|
|
omega3 = 8.185
|
||
|
|
omega4 = 9.083
|
||
|
|
omega5 = 20.29
|
||
|
|
gamma0 = 0.053
|
||
|
|
gamma1 = 3.886
|
||
|
|
gamma2 = 0.452
|
||
|
|
gamma3 = 0.065
|
||
|
|
gamma4 = 0.916
|
||
|
|
gamma5 = 2.419
|
||
|
|
Eps(w) = f0 * (omegap * eh)^2 / ((omega0 * eh)^2 - w^2 + I * w * gamma0) + f1 * (omegap * eh)^2 / ((omega1 * eh)^2 - w^2 + I * w * gamma1) + f2 * (omegap * eh)^2 / ((omega2 * eh)^2 - w^2 + I * w * gamma2) + f3 * (omegap * eh)^2 / ((omega3 * eh)^2 - w^2 + I * w * gamma3) + f4 * (omegap * eh)^2 / ((omega4 * eh)^2 - w^2 + I * w * gamma4) + f5 * (omegap * eh)^2 / ((omega5 * eh)^2 - w^2 + I * w * gamma5);
|
||
|
|
ENDMATERIAL
|
||
|
|
|
||
|
|
MATERIAL LDGold
|
||
|
|
eh = 1519267460583196.5;
|
||
|
|
omegap = 9.03
|
||
|
|
f0 = 0.76
|
||
|
|
f1 = 0.024
|
||
|
|
f2 = 0.010
|
||
|
|
f3 = 0.071
|
||
|
|
f4 = 0.601
|
||
|
|
f5 = 4.384
|
||
|
|
omega0 = 0.0
|
||
|
|
omega1 = 0.415
|
||
|
|
omega2 = 0.83
|
||
|
|
omega3 = 2.969
|
||
|
|
omega4 = 4.304
|
||
|
|
omega5 = 13.32
|
||
|
|
gamma0 = 0.053
|
||
|
|
gamma1 = 0.241
|
||
|
|
gamma2 = 0.345
|
||
|
|
gamma3 = 0.87
|
||
|
|
gamma4 = 2.494
|
||
|
|
gamma5 = 2.214
|
||
|
|
Eps(w) = f0 * (omegap * eh)^2 / ((omega0 * eh)^2 - w^2 + I * w * gamma0) + f1 * (omegap * eh)^2 / ((omega1 * eh)^2 - w^2 + I * w * gamma1) + f2 * (omegap * eh)^2 / ((omega2 * eh)^2 - w^2 + I * w * gamma2) + f3 * (omegap * eh)^2 / ((omega3 * eh)^2 - w^2 + I * w * gamma3) + f4 * (omegap * eh)^2 / ((omega4 * eh)^2 - w^2 + I * w * gamma4) + f5 * (omegap * eh)^2 / ((omega5 * eh)^2 - w^2 + I * w * gamma5);
|
||
|
|
ENDMATERIAL
|
||
|
|
|