qpms/drudeparam_data.c

   1 /*! \file drudeparam_data.c
   2  * \brief Lorentz-Drude model parameters for various materials permittivity.
   3  *
   4  * Data taken from drudelorentz.com, originally apparently from
   5  * Rakić & al., Appl. Opt. 22 (1198), pp. 5271–5283,
   6  * Optical Properties of Metallic Films for Vertical-Cavity Optoelectronic Device
   7  * https://dx.doi.org/10.1364/AO.37.005271
   8  */
   9
  10 #include "materials.h"
  11 #include <gsl/gsl_const_mksa.h>
  12
  13 #define EH (GSL_CONST_MKSA_ELECTRON_VOLT / GSL_CONST_MKSA_PLANCKS_CONSTANT_HBAR)
  14
  15 static const qpms_ldparams_t LDPARAMS_AU = {
  16   1, // eps_inf
  17   9.03*EH, // omega_p
  18   6, // n
  19   {
  20     {0.76 ,  0.   *EH, 0.053*EH},
  21     {0.024,  0.415*EH, 0.241*EH},
  22     {0.01 ,  0.83 *EH, 0.345*EH},
  23     {0.071,  2.969*EH, 0.87 *EH},
  24     {0.601,  4.304*EH, 2.494*EH},
  25     {4.384, 13.32 *EH, 2.214*EH}
  26   }
  27 };
  28 const qpms_ldparams_t *const QPMS_LDPARAMS_AU = &LDPARAMS_AU;
  29
  30 static const qpms_ldparams_t LDPARAMS_AG = {
  31   1, // eps_inf
  32   9.01*EH, // omega_p
  33   6, // n
  34   {
  35     {0.845 ,  0.  *EH, 0.048*EH},
  36     {0.065,  0.816*EH, 3.886*EH},
  37     {0.124,  4.481*EH, 0.452*EH},
  38     {0.011,  8.185*EH, 0.065*EH},
  39     {0.84 ,  9.083*EH, 0.916*EH},
  40     {5.646, 20.29 *EH, 2.419*EH}
  41   }
  42 };
  43 const qpms_ldparams_t *const QPMS_LDPARAMS_AG = &LDPARAMS_AG;
  44
  45 static const qpms_ldparams_t LDPARAMS_CU = {
  46   1, // eps_inf
  47   10.83*EH, // omega_p
  48   5, // n
  49   {
  50     {0.575,  0.   *EH, 0.03 *EH},
  51     {0.061,  0.291*EH, 0.378*EH},
  52     {0.104,  2.957*EH, 1.056*EH},
  53     {0.723,  5.3  *EH, 3.213*EH},
  54     {0.638, 11.18 *EH, 4.305*EH}
  55   }
  56 };
  57 const qpms_ldparams_t *const QPMS_LDPARAMS_CU = &LDPARAMS_CU;
  58
  59 static const qpms_ldparams_t LDPARAMS_AL = {
  60   1, // eps_inf
  61   14.98*EH, // omega_p
  62   5, // n
  63   {
  64     {0.523, 0.   *EH, 0.047*EH},
  65     {0.227, 0.162*EH, 0.333*EH},
  66     {0.05 , 1.544*EH, 0.312*EH},
  67     {0.166, 1.808*EH, 1.351*EH},
  68     {0.03 , 3.473*EH, 3.382*EH},
  69   }
  70 };
  71 const qpms_ldparams_t *const QPMS_LDPARAMS_AL = &LDPARAMS_AL;
  72
  73 static const qpms_ldparams_t LDPARAMS_CR = {
  74   1, // eps_inf
  75   10.75*EH, // omega_p
  76   5, // n
  77   {
  78     {0.168, 0.   *EH, 0.047*EH},
  79     {0.151, 0.121*EH, 3.175*EH},
  80     {0.15 , 0.543*EH, 1.305*EH},
  81     {1.149, 1.97 *EH, 2.676*EH},
  82     {0.825, 8.775*EH, 1.335*EH},
  83   }
  84 };
  85 const qpms_ldparams_t *const QPMS_LDPARAMS_CR = &LDPARAMS_CR;
  86
  87 static const qpms_ldparams_t LDPARAMS_TI = {
  88   1, // eps_inf
  89   7.29*EH, // omega_p
  90   5, // n
  91   {
  92     {0.148,  0.   *EH, 0.082*EH},
  93     {0.899,  0.777*EH, 2.276*EH},
  94     {0.393,  1.545*EH, 2.518*EH},
  95     {0.187,  2.509*EH, 1.663*EH},
  96     {0.001, 19.43 *EH, 1.762*EH},
  97   }
  98 };
  99 const qpms_ldparams_t *const QPMS_LDPARAMS_TI = &LDPARAMS_TI;
 100
 101 static const qpms_ldparams_t LDPARAMS_BE = {
 102   1, // eps_inf
 103   18.51*EH, // omega_p
 104   5, // n
 105   {
 106     {0.084, 0.   *EH,   0.035*EH},
 107     {0.031, 0.1  *EH,   1.664*EH},
 108     {0.14 , 1.032*EH,   3.395*EH},
 109     {0.53 , 3.183*EH,   4.454*EH},
 110     {0.13 , 4.604*EH,   1.802*EH},
 111   }
 112 };
 113 const qpms_ldparams_t *const QPMS_LDPARAMS_BE = &LDPARAMS_BE;
 114
 115 static const qpms_ldparams_t LDPARAMS_NI = {
 116   1, // eps_inf
 117   15.92*EH, // omega_p
 118   5, // n
 119   {
 120     {0.096, 0.   *EH, 0.048*EH},
 121     {0.100, 0.174*EH, 4.511*EH},
 122     {0.135, 0.582*EH, 1.334*EH},
 123     {0.106, 1.597*EH, 2.178*EH},
 124     {0.729, 6.089*EH, 6.292*EH},
 125   }
 126 };
 127 const qpms_ldparams_t *const QPMS_LDPARAMS_NI = &LDPARAMS_NI;
 128
 129 static const qpms_ldparams_t LDPARAMS_PD = {
 130   1, // eps_inf
 131   9.72*EH, // omega_p
 132   5, // n
 133   {
 134     {0.33 , 0.   *EH, 0.008*EH},
 135     {0.649, 0.336*EH, 2.95 *EH},
 136     {0.121, 0.501*EH, 0.555*EH},
 137     {0.638, 1.659*EH, 4.621*EH},
 138     {0.453, 5.715*EH, 3.236*EH},
 139   }
 140 };
 141 const qpms_ldparams_t *const QPMS_LDPARAMS_PD = &LDPARAMS_PD;
 142
 143 static const qpms_ldparams_t LDPARAMS_PT = {
 144   1, // eps_inf
 145   9.59*EH, // omega_p
 146   5, // n
 147   {
 148     {0.333, 0.   *EH, 0.080*EH},
 149     {0.191, 0.78 *EH, 0.517*EH},
 150     {0.659, 1.314*EH, 1.838*EH},
 151     {0.547, 3.141*EH, 3.668*EH},
 152     {3.576, 9.249*EH, 8.517*EH},
 153   }
 154 };
 155 const qpms_ldparams_t *const QPMS_LDPARAMS_PT = &LDPARAMS_PT;
 156
 157 static const qpms_ldparams_t LDPARAMS_W = {
 158   1, // eps_inf
 159   13.22*EH, // omega_p
 160   5, // n
 161   {
 162     {0.206, 0.   *EH, 0.064 *EH},
 163     {0.054, 1.004*EH, 0.53  *EH},
 164     {0.166, 1.917*EH, 1.281 *EH},
 165     {0.706, 3.58 *EH, 3.332 *EH},
 166     {2.59 , 7.498*EH, 5.836 *EH},
 167   }
 168 };
 169 const qpms_ldparams_t *const QPMS_LDPARAMS_W = &LDPARAMS_W;