U-Fhyfig4.ppi

SPECIATION
  jobTitle                             "Uranium adsorption by ferrihydrite in the presence of carbonate"
  calculationType                      custom
  calculationMethod                    1
  database                             050000c0.tdb   # NEA database from JAEA
  xmin                                 2.0
  xmax                                 10.0
  resolution                           71
  loopmin                              58
  loopmax                              580
  loopint                              522
  numericTags                          <pO2> = "-0.7" \
                                       <pHinit> = 2 \
                                       <mass> = <loop>/650

PLOT
  plotTitle                            "Carbonate adsorption by ferrihydrite<br> \
                                       Hiemstra et al. (2009), Fig. 4 after Zachara et al. (1987)"
  customXcolumn                        pH
  lines                                %sorbed
  labels                               "58 m<sup>2</sup>/L" "580 m<sup>2</sup>/L"
  xtitle                               pH
  ytitle                               "% CO<sub>3</sub> sorbed"
  legendtextsize                       0
  extraText                            extratextfig4.dat 

CHEMISTRY

include cdmusic_hiemstra.dat 

PRINT
 -reset false

PHASES
Fix_H+
   H+ = H+
   log_k 0.0

UO2(OH)2(s)
   UO2(OH)2 + 2H+ = UO2+2 + 2H2O
   log_k 5.3

SELECTED_OUTPUT
    -high_precision       true
    -reset                false

USER_PUNCH
-headings pH logU %sorbed
10 PUNCH -la("H+"),  log10(TOT("C")), 100*SURF("C","Fhy")/SYS("C")
SOLUTION 1
   temp      25
   pH        <pHinit>
   pe        14
   units     mol/kgw
   water     1 # kg
   N         0.1
   Na        0.1 charge
   C         4.6e-6

SURFACE 1
  Fhy_unicOHH0.5 3.5  650 <mass>     # sites/nm2  m2/g  g
  Fhy_unieOHH0.5 2.5
  Fhy_triOH0.5   1.2
    -cap   1.15  0.9              # C1  C2  (in F/m2)
    -cd_music
    -sites_units density

EQUILIBRIUM_PHASES
   O2(g)     <pO2>
   Fix_H+  -<x_axis>  NaOH
      -force_equality true
   UO2(OH)2(s) 0 0

END