<?xml version="1.0"?>

<data>

<geometry>
 
 <grid> <!-- Grid --> <!-- Warning! <grid> should be before the <structure> -->
  <var name = "gpn_x" value = "81" units = "n/a" />
  <var name = "gpn_y" value = "41" units = "n/a" />
  <var name = "gpn_z" value = "41" units = "n/a" />
  
  <var name = "h_common" value = "0.4" units = "nm" />
  <var name = "h_x" value = "1.0" units = "n/a" />
  <var name = "h_y" value = "1.0" units = "n/a" />
  <var name = "h_z" value = "1.0" units = "n/a" />
 </grid>
 
 <structure> <!-- Warning! <structure> should be after the <grid> -->
  
  <!-- Electrolyte -->
  <region id = "4">
   <plus>
    <cuboid>
     <xlim start = "0.0" end = "1.0" units = "fraction" />
    </cuboid>
   </plus>
  </region>
  
  <!-- SiO2 layer -->
  <region id = "3">
   <plus>
    <cuboid>
     <xlim start = "0.21875" end = "0.78125" units = "fraction" />
    </cuboid>
   </plus>

   <minus>
     <!-- pore -->
    <cuboid>
     <xlim start = "0.21875" end = "0.78125" units = "fraction" />
     <ylim start = "0.4375" end = "0.5625" units = "fraction" />
     <zlim start = "0.4375" end = "0.5625" units = "fraction" />
    </cuboid>
   </minus>

  </region>
 
  
 </structure>

</geometry>

<materials>
 <mat id = "1" type = "Si-n">
  <var name = "eps" value = "11.70" units = "n/a" />
  <var name = "rho_D_p" value = "2.0e16" units = "1/cm^3" /> <!-- Si-n, Donors -->
  <var name = "rho_A_n" value = "0.0" units = "1/cm^3" /> <!-- Si-p, Acceptors -->
  <var name = "E_gap" value = "1.124" units = "eV" /> <!-- in Si 1.124 eV  -->
 </mat>
 
 <mat id = "2" type = "Si-p">
  <var name = "eps" value = "11.70" units = "n/a" />
  <var name = "rho_D_p" value = "0.0" units = "1/cm^3" /> <!-- Si-n, Donors -->
  <var name = "rho_A_n" value = "2.0e16" units = "1/cm^3" /> <!-- Si-p, Acceptors -->
  <var name = "E_gap" value = "1.124" units = "eV" /> <!-- in Si 1.124 eV  -->
 </mat>
 
 <mat id = "3" type = "SiO2">
  <var name = "eps" value = "3.9" units = "n/a" />
  <var name = "rho_D_p" value = "0.0" units = "1/cm^3" /> <!-- Si-n, Donors -->
  <var name = "rho_A_n" value = "0.0" units = "1/cm^3" /> <!-- Si-p, Acceptors -->
 </mat>
 
 <mat id = "4" type = "H2O">
  <var name = "eps" value = "78.0" units = "n/a" />
  <var name = "rho_D_p" value = "0.0" units = "1/cm^3" /> <!-- Si-n, Donors -->
  <var name = "rho_A_n" value = "0.0" units = "1/cm^3" /> <!-- Si-p, Acceptors -->
 </mat>
</materials>

<constants> 
 <!-- General constants and variables -->
 <var name = "eps_0" value = "5.52635e-2" units = "e^2/eVnm" /> <!-- 5.52635e-2 e^2/(eV nm) = 5.52635e7 e^2/(eV m) =  8.85e-12 C^2/(J m) -->
 <var name = "Temperature" value = "300" units = "K" />
 <var name = "k_B" value = "8.617343e-5" units = "eV/K" /> <!-- 8.617343(15)e-5 eV/K  =  1.3806504e-23 J/K -->
 <var name = "q_elem" value = "1.0" units = "e" /> <!-- 1e = 1.602e-19 C; -->
 <var name = "h_Planck" value = "4.1356673310e-15" units = "eVs" /> <!-- eVs  6.626e-34 J*s; also? h = 4.13566733(10)e-15 eV*s -->
 
 <!-- Mobility -->
 <!-- mu_n = mu_n_0 * (Temperature/300K)^(-mu_n_alpha), Temperature = 300K by default  -->
 <var name = "mu_n_0" value = "1448.0" units = "cm^2/Vs" /> <!-- cm^2/Vs -->
 <var name = "mu_n_alpha" value = "2.33" units = "n/a" /> <!-- n/a -->
 <var name = "mu_p_0" value = "473.0" units = "cm^2/Vs" /> <!-- cm^2/Vs -->
 <var name = "mu_p_alpha" value = "2.23" units = "n/a" /> <!-- n/a -->
 
 <!-- Semiconductor -->
 <var name = "m_electron_0" value = "9.109e-31" units = "kg" /> <!-- 9.109e-31kg -->
 <var name = "m_electron_eff_relative" value = "1.08" units = "n/a" /> <!-- 1.08 n/a -->
 <var name = "m_hole_eff_relative" value = "0.57" units = "n/a" /> <!-- 0.57 n/a --> 
 <var name = "rho_A_n_max" value = "2.0e16" units = "cm^-3" /> <!-- 5.0e16 cm^-3 --> 
 <var name = "rho_D_p_max" value = "2.0e16" units = "cm^-3" /> <!-- 5.0e16 cm^-3 --> 
 <var name = "E_gap" value = "1.124" units = "eV" /> <!-- 1.124 eV --> 
 
 <!-- Recombination -->
 <!-- R_SRH_tau_n = R_SRH_tau_n_0 / (1 + (rho_A_n_max/R_SRH_N_ref)); -->
 <!-- R_SRH_tau_n, R_SRH_tau_p are characteristic times -->
 <var name = "R_SRH_N_ref" value = "5.0e16" units = "cm^-3" /> <!-- 5.0e16 cm^-3 --> 
 <var name = "R_SRH_tau_n_0" value = "5.0e-5" units = "s" /> <!-- 5.0e-5 s --> 
 <var name = "R_SRH_tau_p_0" value = "5.0e-5" units = "s" /> <!-- 5.0e-5 s -->
 
</constants>

</data>