#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 1.18e-4 0.1]}}}
array set $Base {Cstar {Double {[Arrhenius [expr 5.0e22*exp(9.0)] 3.97]}}}

#assume diffusion limited recombination rates...
array set $Base {RecombDefect {String Int}}
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.145]}}}
array set $Base {positive {Double {[Arrhenius 5.68 0.455]}}}
array set $Base {dnegative {Double {[Arrhenius 32.47 0.62]}}}

#volume - collapsed around the vacancy
#volume is negative and equal to 2 * pi * r2 * 3 * (1-nu)/1+nu) * Gam / mu
array set $Base {Volume {Double {-2*3.14159*2.47e-8*2.47e-8*(1.615)*705/7.9e11}}}

#Loop Parameters
array set $Base {gb {Double 1.4285}}
array set $Base {K {Double 0.0}}

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