examples/scuff-em-benchmarks/materials/matprop.dat

   1 # Drude-Lorentz models, with the same constants as in qpms.
   2 # Compared to qpms, sign of the imaginary part is swapped here due
   3 # to different time-frequency transform convention.
   4
   5 MATERIAL LDSilver
   6 # Electronvolt over reduced Planck's constant; the omegas and gammas
   7 # are defined in these units.
   8    eh = 1519267460583196.5;
   9    omegap = 9.01
  10    f0 = 0.84
  11    f1 = 0.065
  12    f2 = 0.124
  13    f3 = 0.111
  14    f4 = 0.84
  15    f5 = 5.646
  16    omega0 = 0.0
  17    omega1 = 0.816
  18    omega2 = 4.481
  19    omega3 = 8.185
  20    omega4 = 9.083
  21    omega5 = 20.29
  22    gamma0 = 0.053
  23    gamma1 = 3.886
  24    gamma2 = 0.452
  25    gamma3 = 0.065
  26    gamma4 = 0.916
  27    gamma5 = 2.419
  28    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);
  29 ENDMATERIAL
  30
  31 MATERIAL LDGold
  32    eh = 1519267460583196.5;
  33    omegap = 9.03
  34    f0 = 0.76
  35    f1 = 0.024
  36    f2 = 0.010
  37    f3 = 0.071
  38    f4 = 0.601
  39    f5 = 4.384
  40    omega0 = 0.0
  41    omega1 = 0.415
  42    omega2 = 0.83
  43    omega3 = 2.969
  44    omega4 = 4.304
  45    omega5 = 13.32
  46    gamma0 = 0.053
  47    gamma1 = 0.241
  48    gamma2 = 0.345
  49    gamma3 = 0.87
  50    gamma4 = 2.494
  51    gamma5 = 2.214
  52    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);
  53 ENDMATERIAL
  54