#set the interstitial model information
array set $Base {DiffModel {Switch 0 {Equilibrium Numeric}}}
array set $Base {Discrete {Switch 2 {linear log s2 s4}}}
array set $Base {Abs.Error {Double 1.0e7}}

#Values from Tomas at LLNL
array set $Base {D0 {Double {[Arrhenius 0.138 1.37]}}}
array set $Base {Cstar {Double {[Arrhenius [expr 5.0e22*exp(11.2)] 3.7]}}}

#assume diffusion limited recombination rates...
array set $Base {RecombDefect {String Vac}}
array set $Base {Kbulk {Double {4.0 * 3.14159 * ([pdbGetDouble Si Int D0] + [pdbGetDouble Si Vac D0]) * [pdbGetDouble Si LatticeSpacing]}}}

#charge states
array set $Base {neutral {Double 1.0}}
array set $Base {negative {Double {[Arrhenius 5.68 0.48]}}}
array set $Base {positive {Double {[Arrhenius 5.68 0.42]}}}

#volume - inclusion model, so volume is positive
#dilation * 4 pi r3
array set $Base {Volume {Double {0.5*4*3.14159*1.11e-8*1.11e-8*1.11e-8}}}

#Loop parameters
array set $Base {gb {Double 0.7}}
array set $Base {K {Double 5e6}}

return { DiffModel Discrete Abs.Error D0 Capture Kbulk neutral negative positive Volume gb K}
