Thermodynamic integration: redox potential: Difference between revisions
Move the TPT note under Step 4, where the GGA integration it follows on from is described |
|||
| (30 intermediate revisions by the same user not shown) | |||
| Line 1: | Line 1: | ||
Redox potentials of an electrochemical half-cell can be computed by thermodynamic integration, following the procedure of Jinnouchi, Karsai, and Kresse {{Cite|jinnouchi:ncm:2024}}. The Fe<sup>3+</sup>/Fe<sup>2+</sup> half-cell serves as the example throughout (''oxidized state/reduced state''): | |||
<math>\mathrm{Fe^{3+} + e^- \;\rightleftharpoons\; Fe^{2+}}</math> | <math>\mathrm{Fe^{3+} + e^- \;\rightleftharpoons\; Fe^{2+}}</math> | ||
| Line 13: | Line 13: | ||
where U<sub>0</sub> and U<sub>1</sub> are the potential energies of the non-interacting and interacting states, respectively, λ is a coupling parameter, μ is the chemical potential (relative to vacuum), and ''n'' is the number of electrons in the half-cell reaction. | where U<sub>0</sub> and U<sub>1</sub> are the potential energies of the non-interacting and interacting states, respectively, λ is a coupling parameter, μ is the chemical potential (relative to vacuum), and ''n'' is the number of electrons in the half-cell reaction. | ||
Three thermodynamic integrations with machine-learned force fields (MLFFs) are described from ''non-interacting'' (λ=0) to ''interacting'' (λ=1) systems: | |||
{|cellpadding="5" cellspacing="0" border="1" | {|cellpadding="5" cellspacing="0" border="1" | ||
| '''TI''' | | '''TI''' | ||
| Line 23: | Line 20: | ||
| '''Term''' | | '''Term''' | ||
|- | |- | ||
| [[ | | [[Thermodynamic integration between machine-learned force fields|MLFF:MLFF]] | ||
| ''Oxidized'' (Fe<sup>3+</sup>) | | ''Oxidized'' (Fe<sup>3+</sup>) | ||
| ''Reduced'' (Fe<sup>2+</sup>) | | ''Reduced'' (Fe<sup>2+</sup>) | ||
| <math>\Delta A_{ML}</math> | | <math>\Delta A_{\mathrm{ML}}</math> | ||
|- | |- | ||
| [[ | | [[Thermodynamic integration between a machine-learned force field and a density functional|MLFF:GGA]] | ||
| MLFF (Fe<sup>3+</sup>) | | MLFF (Fe<sup>3+</sup>) | ||
| RPBE+D3 (Fe<sup>3+</sup>) | | RPBE+D3 (Fe<sup>3+</sup>) | ||
| <math>\Delta A^{ | | <math>\Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_0</math> | ||
|- | |- | ||
| [[ | | [[Thermodynamic integration between a machine-learned force field and a density functional|MLFF:GGA]] | ||
| MLFF (Fe<sup>2+</sup>) | | MLFF (Fe<sup>2+</sup>) | ||
| RPBE+D3 (Fe<sup>2+</sup>) | | RPBE+D3 (Fe<sup>2+</sup>) | ||
| <math>\Delta A^{ | | <math>\Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_1</math> | ||
|- | |- | ||
|} | |} | ||
Finally, the absolute redox potential on a GGA-level <math>\Delta A^{\mathrm{FP}_\mathrm{GGA}}</math> is calculated: | |||
{{ | |||
:<math>\Delta A^{\mathrm{FP}_\mathrm{GGA}} = \Delta A ^{\mathrm{ML}} + \Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_1 - \Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_0</math> | |||
== Input files == | == Input files == | ||
=== {{FILE|POSCAR}}s === | === {{FILE|POSCAR}}s === | ||
The chemical potential requires four different systems: bulk H<sub>2</sub>O with 64 molecules (<code>64H2O_bulk</code>), Fe<sup>3+</sup> in 64 H<sub>2</sub>O (<code>Fe3P_64H2O</code>), Fe<sup>2+</sup> in 64 H<sub>2</sub>O (<code>Fe2P_64H2O</code>), H<sub>2</sub>O slab with 128 molecules (<code>128H2O_slab</code>) | The chemical potential requires four different systems: bulk H<sub>2</sub>O with 64 molecules (<code>64H2O_bulk</code>), Fe<sup>3+</sup> in 64 H<sub>2</sub>O (<code>Fe3P_64H2O</code>), Fe<sup>2+</sup> in 64 H<sub>2</sub>O (<code>Fe2P_64H2O</code>), H<sub>2</sub>O slab with 128 molecules (<code>128H2O_slab</code>), that is, a bulk water reference, <math>[\mathrm{Fe}(\mathrm{H}_2\mathrm{O})_n]^{3+}</math>, <math>[\mathrm{Fe}(\mathrm{H}_2\mathrm{O})_n]^{2+}</math>, and the water slab. | ||
[[File:Redox_structures.png|800px|thumb|center|'''Figure 1'''. Top left: Fe<sup>3+</sup> in 64 H<sub>2</sub>O (<code>Fe3P_64H2O</code>), top center: Fe<sup>2+</sup> in 64 H<sub>2</sub>O (<code>Fe2P_64H2O</code>), top right: bulk H<sub>2</sub>O with 64 molecules (<code>64H2O_bulk</code>), and bottom: H<sub>2</sub>O slab with 128 molecules (<code>128H2O_slab</code>).]] | [[File:Redox_structures.png|800px|thumb|center|'''Figure 1'''. Top left: Fe<sup>3+</sup> in 64 H<sub>2</sub>O (<code>Fe3P_64H2O</code>), top center: Fe<sup>2+</sup> in 64 H<sub>2</sub>O (<code>Fe2P_64H2O</code>), top right: bulk H<sub>2</sub>O with 64 molecules (<code>64H2O_bulk</code>), and bottom: H<sub>2</sub>O slab with 128 molecules (<code>128H2O_slab</code>).]] | ||
| Line 262: | Line 252: | ||
</div> | </div> | ||
<div class="toccolours mw-customtoggle-poscar-fe3p-64h2o">'''Click to reveal the < | <div class="toccolours mw-customtoggle-poscar-fe3p-64h2o">'''Click to reveal the <math>[\mathrm{Fe}(\mathrm{H}_2\mathrm{O})_n]^{3+}</math> POSCAR'''</div> | ||
<div class="mw-collapsible mw-collapsed" id="mw-customcollapsible-poscar-fe3p-64h2o"> | <div class="mw-collapsible mw-collapsed" id="mw-customcollapsible-poscar-fe3p-64h2o"> | ||
| Line 468: | Line 458: | ||
</div> | </div> | ||
<div class="toccolours mw-customtoggle-poscar-fe2p-64h2o">'''Click to reveal the < | <div class="toccolours mw-customtoggle-poscar-fe2p-64h2o">'''Click to reveal the <math>[\mathrm{Fe}(\mathrm{H}_2\mathrm{O})_n]^{2+}</math> POSCAR'''</div> | ||
<div class="mw-collapsible mw-collapsed" id="mw-customcollapsible-poscar-fe2p-64h2o"> | <div class="mw-collapsible mw-collapsed" id="mw-customcollapsible-poscar-fe2p-64h2o"> | ||
| Line 1,072: | Line 1,062: | ||
=== {{FILE|INCAR}} === | === {{FILE|INCAR}} === | ||
The {{FILE|INCAR}} | The {{FILE|INCAR}} file differs from step to step. Each is given on the how-to page that the corresponding step links to. | ||
=== {{FILE|KPOINTS}} === | === {{FILE|KPOINTS}} === | ||
The {{FILE|KPOINTS}} file contains only the Gamma | The {{FILE|KPOINTS}} file contains only the Γ point: | ||
Gamma-point only | Gamma-point only | ||
| Line 1,089: | Line 1,079: | ||
*<code>PAW_PBE Fe_sv 23Jul2007</code> | *<code>PAW_PBE Fe_sv 23Jul2007</code> | ||
This is different to in Ref. {{Cite|jinnouchi: | This is different to in Ref. {{Cite|jinnouchi:ncm:2024}}, where GW {{FILE|POTCAR}}s are used. | ||
== Step-by-step instructions == | == Step-by-step instructions == | ||
=== 1 | === Step 1: Preparing the MLFFs === | ||
First, [[ | First, [[Training machine-learned force fields for redox reactions|train the MLFFs]] for the water slab (<code>128H2O_slab</code>) and the redox pairs (<math>[\mathrm{Fe}(\mathrm{H}_2\mathrm{O})_n]^{3+}</math> and <math>[\mathrm{Fe}(\mathrm{H}_2\mathrm{O})_n]^{2+}</math>). | ||
=== 2 | === Step 2: Calculating the vacuum potential === | ||
Run a molecular dynamics (MD) simulation, then calculate the local potential <math>\phi</math> of the water slab. Using this, the difference between the Oxygen ''1s'' level between the water slab and the water molecules far away from the ions in the redox cells can be calculated. This provides the [[ | Run a molecular dynamics (MD) simulation, then calculate the local potential <math>\phi</math> of the water slab. Using this, the difference between the Oxygen ''1s'' level between the water slab and the water molecules far away from the ions in the redox cells can be calculated. This provides the [[Vacuum reference|chemical potential relative to the vacuum]] <math>e \Delta \bar{\phi}</math>. | ||
{{NB|important|This step requires using GGAs on many structures, so will be more time-consuming.}} | {{NB|important|This step requires using GGAs on many structures, so will be more time-consuming.}} | ||
=== 3 | === Step 3: TI from oxidized to reduced state === | ||
Perform a thermodynamic integration from the [[ | Perform a thermodynamic integration from the [[Thermodynamic integration between machine-learned force fields|oxidized state (λ=0; Fe<sup>3+</sup>) to the reduced state (λ=1; Fe<sup>2+</sup>)]]. | ||
{{NB|mind|This is a thermodynamic integration between two MLFFs, so is very quick.}} | {{NB|mind|This is a thermodynamic integration between two MLFFs, so is very quick.}} | ||
=== 4 | === Step 4: TI from MLFF to GGA === | ||
Perform a thermodynamic integration from [[ | Perform a thermodynamic integration from [[Thermodynamic integration between a machine-learned force field and a density functional|MLFF (λ=0) to RPBE+D3 (λ=1)]]. This should be done separately for Fe<sup>3+</sup> and Fe<sup>2+</sup>. | ||
{{NB|important|This step requires using GGAs on many structures so will be more time consuming.}} | {{NB|important|This step requires using GGAs on many structures so will be more time consuming.}} | ||
{{NB|mind|Thermodynamic perturbation theory (TPT) from RPBE+D3 to PBE0, for both Fe<sup>3+</sup> and Fe<sup>2+</sup>, is not covered here because of its cost. It can be replicated using Δ-learning from RPBE to PBE0, using the same procedure as from MLFF to RPBE.}} | |||
=== 5 | === Step 5: Analysis === | ||
Once each step is completed, you will have the following values: | Once each step is completed, you will have the following values: | ||
* the chemical potential - <math>e \Delta \bar{\phi}</math> | * the chemical potential - <math>e \Delta \bar{\phi}</math> | ||
* the free energy difference from MLFF to MLFF (with chemical potential reference) | * the free energy difference from MLFF to MLFF (with chemical potential reference) | ||
** from Fe<sup>3+</sup> to Fe<sup>2+</sup> - <math>\Delta A_{ML}</math> | ** from Fe<sup>3+</sup> to Fe<sup>2+</sup> - <math>\Delta A_{\mathrm{ML}}</math> | ||
* the free energy difference from MLFF to GGA | * the free energy difference from MLFF to GGA | ||
** Reduced: Fe<sup>2+</sup> - <math>\Delta A^{ | ** Reduced: Fe<sup>2+</sup> - <math>\Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_1</math> | ||
** Oxidized: Fe<sup>3+</sup> - <math>\Delta A^{ | ** Oxidized: Fe<sup>3+</sup> - <math>\Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_0</math> | ||
:<math>\Delta A^{\mathrm{FP}_\mathrm{GGA}} = \Delta A ^{\mathrm{ML}} + \Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_1 - \Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_0</math> | |||
The vacuum alignment is already contained in <math>\Delta A^{\mathrm{ML}}</math>. Combining the remaining terms as above gives the free energy difference from Fe<sup>3+</sup> to Fe<sup>2+</sup> on the GGA level <math>\Delta A^{\mathrm{FP}_\mathrm{GGA}}</math>. Switch the sign to get the absolute redox potential <math>U^{\mathrm{FP}_\mathrm{GGA}}_{\mathrm{redox}}</math>. Finally, take the standard hydrogen electrode (SHE) (4.44 V {{Cite|trasatti:pac:1986}}) away to compare to experiment <math>U^{\mathrm{FP}_\mathrm{GGA}}_{\mathrm{redox,SHE}}</math>. The results below match the literature to within 30 meV. | |||
The results are summarized in the following table: | |||
{|cellpadding="5" cellspacing="0" border="1" | {|cellpadding="5" cellspacing="0" border="1" | ||
| '''Step''' | | '''Step''' | ||
| '''Term''' | | '''Term''' | ||
| '''Ours (eV)''' | | '''Ours (eV)''' | ||
| '''Lit. {{Cite|jinnouchi: | | '''Lit. {{Cite|jinnouchi:ncm:2024}} (eV)''' | ||
|- | |- | ||
| [[ | | [[Vacuum reference|'''2 - Vacuum reference''']] | ||
| <math>e \Delta \bar{\phi}</math> | | <math>e \Delta \bar{\phi}</math> | ||
| -3.71 | | -3.71 | ||
| -3.74 | | -3.74 | ||
|- | |- | ||
| [[ | | [[Thermodynamic integration between machine-learned force fields|'''3 - TI (MLFF)''']] | ||
| <math>\Delta A_{ML}</math> | | <math>\Delta A_{\mathrm{ML}}</math> | ||
| -4.98 | | -4.98 | ||
| -4.95 | | -4.95 | ||
|- | |- | ||
| [[ | | [[Thermodynamic integration between a machine-learned force field and a density functional|'''4 - TI (GGA)''']] | ||
| <math>\Delta A^{ | | <math>\Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_1 - \Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_0</math> | ||
| -0.055 | | -0.055 | ||
| -0.049 | | -0.049 | ||
|- | |- | ||
| Absolute redox potential | | Absolute redox potential | ||
| <math>U^{ | | <math>U^{\mathrm{FP}_\mathrm{GGA}}_{\mathrm{redox}} </math> | ||
| | | 5.03 | ||
| | | 5.00 | ||
|- | |- | ||
| Redox potential | | Redox potential | ||
| <math>U^{ | | <math>U^{\mathrm{FP}_\mathrm{GGA}}_{\mathrm{redox,SHE}}</math> | ||
| 0.59 | | 0.59 | ||
| 0.56 | | 0.56 | ||
| Line 1,159: | Line 1,146: | ||
== Recommendations and advice == | == Recommendations and advice == | ||
*The | *The vacuum-reference step and the MLFF:GGA integrations are by far the most expensive parts of the workflow, since both evaluate a GGA on many structures. Training the MLFFs in the first place can also be time consuming. The MLFF:MLFF integration, by contrast, is very quick. | ||
<!--*TPT can be performed using Δ-ML. This is not yet available, so it is hidden for now.--> | |||
== Related tags and articles == | == Related tags and articles == | ||
;How-tos | ;How-tos | ||
* [[ | * [[Vacuum reference]] | ||
* [[ | * [[Thermodynamic integration between machine-learned force fields]] | ||
* [[ | * [[Thermodynamic integration between a machine-learned force field and a density functional|Thermodynamic integration between a machine-learned force field and a density functional]] | ||
;Theory | ;Theory | ||
* [[Thermodynamic integration]] | * [[Thermodynamic integration]] | ||
;Tags | ;Tags | ||
* {{TAG|VCAIMAGES}}, {{TAG|NCORE_IN_IMAGE1}}, {{TAG|MDALGO}}, {{TAG|ML_LMLFF}}, {{TAG|NELECT | * {{TAG|VCAIMAGES}}, {{TAG|NCORE_IN_IMAGE1}}, {{TAG|MDALGO}}, {{TAG|ML_LMLFF}}, {{TAG|NELECT}}, {{TAG|ICORELEVEL}}, {{TAG|WRT_POTENTIAL}}, {{TAG|LVHAR}} | ||
== References == | == References == | ||
<references/> | <references/> | ||
[[Category:Howto]][[Category:Advanced molecular-dynamics sampling]][[Category:Machine-learned force fields]] | |||
<!--[[Category:Electrochemistry]]--> | <!--[[Category:Electrochemistry]]--> | ||
Latest revision as of 14:12, 18 September 2026
Redox potentials of an electrochemical half-cell can be computed by thermodynamic integration, following the procedure of Jinnouchi, Karsai, and Kresse [1]. The Fe3+/Fe2+ half-cell serves as the example throughout (oxidized state/reduced state):
[math]\displaystyle{ \mathrm{Fe^{3+} + e^- \;\rightleftharpoons\; Fe^{2+}} }[/math]
The following two models are used:
- [math]\displaystyle{ [\mathrm{Fe}(\mathrm{H}_2\mathrm{O})_n]^{3+} }[/math] - Fe ion in 64 H2O bulk cell with 3 valence electrons put in the background charge
- [math]\displaystyle{ [\mathrm{Fe}(\mathrm{H}_2\mathrm{O})_n]^{2+} }[/math] - Fe ion in 64 H2O bulk cell with 2 valence electrons put in the background charge
Thermodynamic integration (TI) is used to calculate the free energy difference ΔA using the following basic equation:
[math]\displaystyle{ \Delta A = \int_0^1 \langle U_1 - U_0 \rangle_{\lambda} d\lambda - \mu n }[/math]
where U0 and U1 are the potential energies of the non-interacting and interacting states, respectively, λ is a coupling parameter, μ is the chemical potential (relative to vacuum), and n is the number of electrons in the half-cell reaction.
Three thermodynamic integrations with machine-learned force fields (MLFFs) are described from non-interacting (λ=0) to interacting (λ=1) systems:
| TI | λ=0 | λ=1 | Term |
| MLFF:MLFF | Oxidized (Fe3+) | Reduced (Fe2+) | [math]\displaystyle{ \Delta A_{\mathrm{ML}} }[/math] |
| MLFF:GGA | MLFF (Fe3+) | RPBE+D3 (Fe3+) | [math]\displaystyle{ \Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_0 }[/math] |
| MLFF:GGA | MLFF (Fe2+) | RPBE+D3 (Fe2+) | [math]\displaystyle{ \Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_1 }[/math] |
Finally, the absolute redox potential on a GGA-level [math]\displaystyle{ \Delta A^{\mathrm{FP}_\mathrm{GGA}} }[/math] is calculated:
- [math]\displaystyle{ \Delta A^{\mathrm{FP}_\mathrm{GGA}} = \Delta A ^{\mathrm{ML}} + \Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_1 - \Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_0 }[/math]
Input files
POSCARs
The chemical potential requires four different systems: bulk H2O with 64 molecules (64H2O_bulk), Fe3+ in 64 H2O (Fe3P_64H2O), Fe2+ in 64 H2O (Fe2P_64H2O), H2O slab with 128 molecules (128H2O_slab), that is, a bulk water reference, [math]\displaystyle{ [\mathrm{Fe}(\mathrm{H}_2\mathrm{O})_n]^{3+} }[/math], [math]\displaystyle{ [\mathrm{Fe}(\mathrm{H}_2\mathrm{O})_n]^{2+} }[/math], and the water slab.
Fe3P_64H2O), top center: Fe2+ in 64 H2O (Fe2P_64H2O), top right: bulk H2O with 64 molecules (64H2O_bulk), and bottom: H2O slab with 128 molecules (128H2O_slab).Click to reveal the
64H2O_bulk POSCARWater
1
12.42128700 0.00000000 0.00000000
0.00000000 12.42128700 0.00000000
0.00000000 0.00000000 12.42128700
H O
128 64
Direct
0.36257633 0.37953597 0.67785342
0.42947178 0.43286249 0.59010358
0.66870393 1.13989162 0.58842641
0.70823249 1.22041509 0.67807457
0.35637008 0.05839855 0.21517466
0.44580288 0.01834909 0.16004703
0.44380894 0.67759921 0.05971511
0.44857483 0.64047308 0.17235572
1.14522157 0.76633728 -0.07286828
1.26197238 0.78911446 -0.10391832
0.56790416 1.14251886 0.92143084
0.60754105 1.14691841 0.80551737
0.75633099 0.25670520 0.32073462
0.75786987 0.15087192 0.26471573
0.74642803 -0.03671116 1.07026855
0.68277568 -0.08917103 1.16765095
0.08323566 -0.02515907 -0.01412446
-0.00880432 1.02342521 0.04212898
0.68157988 0.24579303 0.03134658
0.60083186 0.23578159 0.12220439
1.14251010 0.79368279 0.23822678
1.03381447 0.81936982 0.30217442
1.45625953 0.82211861 -0.39364182
1.33249296 0.84676639 -0.38681569
0.73773505 0.37233458 0.47040761
0.62289024 0.34929795 0.43767209
0.63992674 0.70350215 0.78737637
0.67764626 0.66562038 0.90017103
0.38737624 0.10112290 0.97816413
0.44041019 0.20087678 1.04654040
-0.74229635 0.27326641 0.85574288
-0.86480570 0.30089547 0.83717409
0.98781733 0.58627453 1.16338342
0.98642719 0.45910015 1.19257026
0.32814541 0.15987161 0.50362237
0.25642761 0.26648016 0.48295405
0.69811093 -0.56389784 0.28108666
0.70770558 0.42946360 0.15518042
0.18565685 0.52373867 0.56827259
0.08659439 0.58375638 0.51565848
0.44991203 0.27342036 0.48425454
0.55101593 0.27228106 0.56380912
0.41624326 -0.17788427 0.80896103
0.50257275 -0.21360866 0.89707634
0.28998176 0.48372062 -0.01768589
0.19491021 0.39742192 -0.01269715
1.06654096 0.36410409 0.63682306
1.03975716 0.44960909 0.73041188
0.50170068 1.61380686 0.40736049
0.49942205 1.73793230 0.43266624
0.37390409 0.30295427 0.17857302
0.34033093 0.38510006 0.09520283
0.30807694 0.20122962 0.31494909
0.22529499 0.19805563 0.21739075
1.12308748 0.37150659 0.42477646
1.06276778 0.26307066 1.46962386
-0.15792914 0.29610530 0.85397054
-0.07308234 0.27252743 0.77021418
0.85228948 0.16062637 0.95526886
0.89963468 0.18253764 1.06178327
0.19181780 0.04626484 -0.19418583
0.15189262 -0.06556287 -0.18492234
-0.10238306 0.45933074 0.94465075
0.78170215 0.43034159 0.95991840
-0.04726476 0.53527628 0.38000776
-0.10379563 0.60547798 0.46538564
-0.00098972 1.02353900 0.73537767
-0.10239174 1.10640574 0.72833236
0.83303676 0.93517612 0.67263678
0.72416689 0.88061601 0.69336360
0.83396381 0.79727679 0.33276158
0.84935945 0.89986765 1.25490788
0.43767040 0.26474088 0.82690407
0.40610758 0.19731651 0.73463453
-0.22521691 0.99033821 0.52178482
-0.21732111 1.08901868 0.43461921
0.55979128 0.65441008 0.58905315
0.51739017 0.64402769 0.70843828
1.15263203 0.22553743 1.04375491
1.05410314 0.23287904 1.11317058
0.40787600 0.01178034 0.65057706
0.46239323 0.06947748 0.55789921
1.22298533 -0.28915950 0.52942510
0.19053464 -0.18268639 0.46366419
0.34645899 0.68043144 0.72737986
0.36594763 0.55141902 0.69243633
0.56017599 0.33134538 -0.03087852
0.61388449 0.39021229 -0.12522124
0.38664251 -0.13040994 -0.70002278
0.35405937 -0.19423015 -0.59774517
0.73338812 0.70036370 0.46781396
0.80345417 0.77305746 0.53896313
0.70280496 0.86052960 0.87581100
0.72653862 0.98582519 0.88541416
0.94986235 0.10516542 0.53628061
0.99735924 1.10316818 0.41487930
0.34586913 0.92208748 -0.04377027
0.27835954 0.97034641 0.04814922
0.88788268 1.29429166 0.20459804
-0.10649273 1.36074553 0.09165817
0.39371787 0.57519096 -0.11020542
0.46220136 0.48981436 -0.05946583
0.89286511 0.79370028 0.86701692
0.96099685 0.84339162 0.96334218
0.62514345 0.53133485 0.73825559
0.73640530 0.50006566 0.71024894
-0.12902917 0.72677131 1.16320660
-0.11814514 0.70432566 1.03349629
1.01129606 0.61898608 0.90675712
1.09968460 0.53447870 0.92387699
-0.14457334 0.35937352 0.63141164
-0.11499539 0.45568803 0.54735066
0.48813670 0.98589446 0.40760031
0.60110088 1.02927414 0.43954578
-0.37236788 0.57165925 0.04667849
-0.24946342 0.60016434 1.06579556
1.07450858 1.10382861 0.26264566
1.12554235 0.99644744 0.31325744
0.62490856 0.06519849 0.18845978
0.54885406 0.10569865 0.28792416
0.55694560 -0.17713139 0.27541275
0.56084272 -0.21506154 0.16885793
0.12891427 0.50587770 0.28000948
0.23833432 0.43434745 0.29542543
1.48346125 0.46743806 0.19449182
1.57481458 0.52369865 0.24680216
0.17873553 0.62275568 0.11942301
0.26852366 0.68872043 0.07818815
0.35693455 0.43925319 0.62530196
0.64792950 1.20400828 0.63184735
0.39412198 -0.01057009 0.20787038
0.43581616 0.71008914 0.13099142
1.19397105 0.77903373 0.85946991
0.63615555 1.12938446 0.87822715
0.80661353 0.20330608 0.29232873
0.73246578 -0.02838354 1.14624509
0.03244721 -0.04025903 0.03944908
0.63721217 0.28885881 0.08035039
1.11493495 0.82738114 0.30859190
1.39883698 0.85058088 -0.34459284
0.69698494 0.36107968 0.40464588
0.64192914 0.72485162 0.86089438
0.45396944 0.13747103 1.00222318
-0.81120011 0.28079526 0.89372357
1.01585028 0.52957121 1.21343372
0.31821591 0.22470258 0.45731626
0.70888882 -0.51685513 0.22019861
0.13832317 0.58769764 0.57254578
0.51028551 0.31315460 0.51155010
0.42617734 -0.19946999 0.88578770
0.21795347 0.46397897 0.02254897
1.03345078 0.37392845 0.70810584
0.53447426 1.66561969 0.45708230
0.40848000 0.35179327 0.12578669
0.30420794 0.19684477 0.23570844
1.10944135 0.32634895 1.48853144
-0.15068923 0.25159638 0.78858129
0.84350236 0.13666558 1.03077541
0.12621458 0.00799599 -0.18726449
-0.14969491 0.39803889 -0.04386545
-0.05257019 0.54926511 0.45604016
-0.07042394 1.04132825 0.69759096
0.77456016 0.89777696 0.63596725
0.87893918 0.82791839 1.27826204
0.37199961 0.24693838 0.78669458
-0.26629559 1.03776777 0.47472997
0.58383474 0.64578572 0.66290779
1.11758097 0.18829015 1.10755225
0.41441089 0.08360444 0.62404569
0.25236886 -0.21779138 0.50357805
0.38654211 0.61365045 0.74790977
0.54189953 0.36543129 -0.10123356
0.41572204 -0.15596095 -0.63075239
0.80055913 0.73724417 0.46650168
0.76508964 0.91308708 0.87755337
0.94189130 0.13529024 0.46400834
0.29460037 0.97883985 -0.02591476
-0.06126178 1.32786996 0.14953350
0.41172576 0.54824212 -0.03901993
0.94501130 0.77310663 0.92132598
0.69302524 0.50734691 0.77446848
-0.13207205 0.66750977 1.10613301
1.02911948 0.53980433 0.89023275
-0.17248923 0.40767816 0.57425135
0.52217208 1.04992599 0.43651343
-0.30113028 0.57147331 0.01268190
1.11295226 1.07452050 0.32297073
-0.41735391 0.12325009 0.21804170
0.61116931 -0.19691698 0.22556446
0.17536968 0.46512423 0.32953214
1.49704038 0.52822279 0.24728978
0.19113993 0.69295724 0.10392182
Click to reveal the [math]\displaystyle{ [\mathrm{Fe}(\mathrm{H}_2\mathrm{O})_n]^{3+} }[/math] POSCAR
Fe_64H2O
1
12.42282200 0.00000000 0.00000000
0.00000000 12.42282200 0.00000000
0.00000000 0.00000000 12.42282200
H O Fe
128 64 1
Direct
-0.38272023 0.47236734 0.69895078
-0.27203601 0.41866794 0.72914269
0.60916541 0.72796225 0.09306829
0.71070494 0.78725068 0.14392771
1.08865564 0.93820602 0.24600432
1.03755512 0.98120737 0.35350048
0.04276558 1.23801733 0.56678300
0.45616640 0.36926572 0.39344586
-0.16977976 0.40716591 0.53863440
-0.25237329 0.38130525 0.44755164
0.44169261 1.14516978 0.14529032
0.37457743 1.13530301 0.24819417
0.18407518 0.25222959 0.06812304
0.09998538 0.14744161 0.06444770
0.27831295 0.57479968 0.31380844
0.29633366 0.45926581 0.36724418
0.64763154 -0.53184911 0.87002912
0.63022496 0.40977773 0.99107908
0.47906160 0.58853267 0.63100265
0.59945353 0.62329058 0.64816927
0.87712851 0.41898140 1.07505748
0.26314503 -0.13253460 0.01243951
0.16153170 -0.19745195 0.05804701
0.34412600 -0.01810085 0.81692911
0.14554222 0.77130921 -0.25213594
0.21209078 0.79452634 -0.14442072
0.41500251 0.04461290 0.96837924
0.49008610 0.14359786 0.93019679
0.82125847 0.86254586 1.27147476
0.00873717 0.10267666 0.23211446
0.35980746 0.34237108 0.76169366
0.33485343 0.46345006 0.77046483
0.90320251 0.61221204 0.04264507
0.97690745 0.71802735 1.06017819
0.96030764 0.81465128 0.70019704
1.01749178 0.74674122 0.60432752
0.64007416 0.23981183 0.92947124
0.31788389 -0.12374136 0.74827741
0.99321033 0.93631524 0.55767630
0.15487454 0.51028490 -0.09542554
0.66815917 0.79686891 0.72718580
0.76202793 0.70002522 0.75366280
0.13591369 -0.02901335 0.00566347
0.08957800 1.06396642 -0.07143824
0.53392446 -0.11896063 0.23953882
0.48718378 -0.00799262 0.25860106
-0.07163033 1.26809135 0.50762252
0.67405502 0.12105566 0.95083205
0.74551359 0.62681866 0.34775553
0.66650260 0.64357368 0.25540936
0.38872827 0.12804968 0.55656518
0.33782256 0.21579745 0.49248801
0.92231233 -0.04307885 1.00089073
0.79411522 -0.04854941 1.01531498
-0.40705912 1.12847907 0.48347215
0.65420892 1.23981902 0.47957005
0.84599640 -0.09249550 0.84023452
0.80327236 0.97997213 0.74592535
0.88638541 0.47237368 0.72110537
0.88941805 0.35493190 0.74249811
0.55296221 0.93428417 0.55243420
0.51441228 0.92952738 0.43767406
0.73847683 0.51634987 1.03266322
0.40092134 0.63952229 -0.01737781
0.41176789 0.89134362 1.10615330
0.44444253 0.84405868 0.99393367
0.58073572 0.63347610 -0.04423652
0.60854638 0.72933752 -0.11607112
-0.16580359 0.22001380 0.89963070
-0.04391612 0.21085677 0.85654400
0.40533602 0.50000334 0.11555470
0.34998095 0.38355577 0.10448457
-0.24359234 0.77537730 0.45894937
-0.26614437 0.70560002 0.56106171
0.28303298 0.61126441 0.01597951
0.37578327 0.33380444 0.30290136
0.59618764 0.27699048 1.12423906
0.50319683 0.35108728 1.09799129
0.18989659 0.76741117 1.22566037
0.06153560 0.75549602 0.22936925
0.51785592 0.24354330 0.65223118
0.56708719 0.32005973 0.73245604
0.61206395 0.98022340 0.70759384
0.50678625 0.91620503 0.71723010
0.68602440 0.14707860 0.76412395
0.69294218 0.12172781 0.64036372
0.33488986 0.71525634 0.14956129
0.35814895 0.80854785 0.22920363
0.95945912 1.05245660 0.53679791
0.77994512 0.55222218 1.15138214
0.69626269 0.94689443 0.43758767
0.82244291 0.95270203 0.44158402
-0.08533410 0.30668689 0.99691887
0.88689642 0.15254779 0.21721075
0.96565778 0.66203563 0.74352603
0.91899641 0.64391575 0.86105003
0.60192845 0.32627053 0.34541527
0.54412777 0.22739074 0.29056180
0.33801876 0.24481339 -0.03610605
0.27498242 0.34868783 -0.07124872
0.63970642 0.49240211 0.35877270
0.61552019 0.46892953 0.48438447
0.30584874 0.64699933 -0.22951580
0.42787645 0.64881697 -0.19623458
0.05975532 0.12314654 0.75074831
-0.02832220 1.03112946 0.76086242
1.07382577 0.25251612 0.39457391
1.17074101 0.24004400 0.30909053
0.90882853 0.66787644 0.52270722
0.07399163 0.59412808 0.95304999
0.13906082 0.47905225 0.10518101
1.17121656 0.46180650 0.23562581
0.97971011 0.63429872 0.37278291
0.90631471 0.68858591 0.26904690
0.60516786 0.96463888 0.06395544
0.68233031 1.04347277 1.12012494
0.20897747 0.21380315 0.63498998
0.20507421 0.29490009 0.72649781
0.21069366 1.05104536 -0.68658423
0.27938942 1.09406910 0.40078460
0.70263174 0.14455648 1.25970691
0.75730930 0.04722448 1.29768294
0.91307779 0.40130091 0.34978752
0.94369056 0.50994390 0.28013594
0.08933991 0.43228293 0.60852163
1.17837154 0.48628664 0.53072146
1.00504493 0.58007882 0.54073092
0.86291931 -0.13828811 1.15173275
-0.34491075 0.43621986 0.75416516
0.63900568 0.75492759 0.15829659
1.07199795 1.00257563 0.28363155
0.38420131 0.34485841 0.37967695
-0.17342306 0.37643136 0.46438657
0.45349630 1.14005000 0.22505348
0.11869813 0.21939929 0.09604324
0.23622993 0.50759344 0.33909324
0.64875252 0.47771851 0.95286519
0.55590506 0.57693205 0.60050819
0.18736673 -0.15437724 -0.00677806
0.37523088 -0.06818877 0.76435686
0.21811495 0.76245384 -0.22055952
0.42220970 0.12606004 0.96558981
0.85780135 0.81574581 1.21531426
0.30433544 0.39124488 0.77849745
0.95758533 0.65697710 1.00747411
1.03085555 0.79681957 0.66923251
0.69067343 0.18643375 0.90833405
0.96523126 0.16955239 0.22035170
0.68612747 0.72230731 0.74299918
0.08987161 0.03581255 0.00348357
0.46357734 -0.08420990 0.24196675
-0.01222913 1.21468404 0.51335175
0.71162904 0.58691114 0.28804098
0.32268708 0.16090339 0.53975026
0.86234882 -0.08436224 0.98773807
-0.38355188 0.19092778 0.52248232
0.85291582 -0.08281162 0.75894361
0.86323773 0.40597333 0.69118967
0.54292990 0.98123241 0.48667605
0.80728528 0.53397589 1.07560692
0.39528125 0.89790682 1.02768167
0.56047647 0.70886600 -0.05694049
-0.09585039 0.25981162 0.89091800
0.40071565 0.42814971 0.14908255
-0.23238312 0.70278251 0.48864911
0.35169257 0.62880569 0.04279978
0.56212715 0.31299799 1.06345122
0.12066294 0.79773619 1.19966391
0.51154081 0.26173878 0.72839369
0.58450605 0.90652737 0.69203180
0.70957535 0.09399646 0.71195964
0.32002783 0.74390105 0.22241023
0.95247763 0.97831298 0.50380136
0.76168681 0.93342180 0.39565898
-0.07590969 0.35133103 1.07020111
0.90649162 0.63028257 0.78114555
0.61524001 0.25974801 0.30414551
0.30266553 0.31064613 -0.00899828
0.62120647 0.43621160 0.41233540
0.36878355 0.59976322 -0.20951626
0.00592307 0.09406140 0.79582533
1.10947213 0.28743919 0.33114811
0.13898678 0.54360995 0.97023057
1.15168897 0.42574337 0.16437693
0.96099089 0.62982868 0.29120309
0.66514203 1.00695060 0.04564722
0.15635030 0.24896214 0.68388955
0.26078411 1.11001890 0.32632346
0.75255117 0.09346553 1.23417762
0.93550375 0.42521425 0.27361930
1.12555194 0.49730673 0.58995753
0.97951694 0.64650296 0.50479031
1.02535902 0.32276265 0.18627127
Click to reveal the [math]\displaystyle{ [\mathrm{Fe}(\mathrm{H}_2\mathrm{O})_n]^{2+} }[/math] POSCAR
Fe_64H2O
1
12.42282200 0.00000000 0.00000000
0.00000000 12.42282200 0.00000000
0.00000000 0.00000000 12.42282200
H O Fe
128 64 1
Direct
0.12743711 0.48760461 0.68402834
0.07224030 0.59099573 0.69720079
0.39623895 0.90251919 -0.68562329
0.52396846 0.89222892 -0.69202626
0.42964083 0.63656491 0.92024856
0.40471536 0.51747508 0.92244566
0.50949529 0.22553243 0.82975356
0.79695561 -0.29210416 -0.14766098
0.00805766 0.01676094 0.93033834
-0.01969191 0.14077136 0.89733509
0.01773692 1.49221026 -0.45662142
-0.09754205 1.45083180 -0.48707747
0.28568895 0.07977712 0.14809775
0.40460027 0.11889958 0.12632577
-0.44823614 -0.23392047 0.91065326
-0.44655395 -0.27382665 1.03236410
0.55758309 0.63486700 0.51461921
0.54015662 1.52364880 0.45450703
0.67045916 0.85136844 0.58289035
0.71164261 0.89435637 0.47626762
1.16684705 0.13193713 0.89590371
0.33871812 0.20748811 0.69296855
0.38041592 0.31864036 0.75378916
-0.26840230 -0.60381695 0.57433173
0.25751301 1.09723808 -0.39131073
0.28561906 1.17135043 -0.49010853
-0.15559770 0.77405954 1.51947552
-0.03350167 0.77676575 1.52339975
0.75237735 1.02309221 1.12180723
0.04826055 0.15234086 0.05862687
-0.09230147 -0.06114088 1.02402855
-0.07779336 -0.14013315 0.92524744
0.68560087 0.63322667 1.36320530
0.79400335 0.69800156 2.34722324
1.01111679 0.16477327 0.53663862
0.98618018 0.28565957 0.49809698
0.20452861 0.81136804 0.83127978
-0.29824646 -0.48415938 0.54811270
0.37502040 1.26549068 0.41257460
1.30154976 0.64025841 0.11667748
1.75290818 0.59403116 0.98316689
1.65052923 0.63412803 0.91547295
0.75375838 0.00138371 -0.37473678
0.76130867 1.05807834 -0.25911110
-0.05696647 -0.29493528 0.76468746
0.02811791 -0.33267590 0.85098776
0.61065413 0.29830310 0.85951557
0.11763966 0.86822584 0.90729729
0.18919114 0.55109134 0.38993918
0.21099928 0.52226945 0.51012484
0.03269714 0.36292005 0.33650725
0.02620788 0.48115860 0.28181401
1.05541232 0.90157404 0.64524915
1.07849658 0.77961507 0.68820380
-0.62352298 1.63205141 0.30877689
0.43788457 1.73907737 0.26620496
0.64515164 -0.84752695 1.12986454
0.64730542 0.18633377 1.26064334
0.19923207 -0.08426170 0.07502044
0.22493320 -0.20749942 0.08281964
0.78950677 1.18291196 0.57758905
0.83466268 1.23338961 0.48038196
1.10490483 1.28683406 0.74494774
0.29928808 0.66864147 -0.17842015
0.69120317 0.84829619 1.11035412
0.67223805 0.75611954 1.18775106
-0.00106366 0.61277354 0.38525072
0.01326608 0.68609575 0.28998127
-0.07137268 0.05723237 0.64512696
0.03262162 0.06902330 0.71661499
0.77907316 0.40440365 0.75246685
0.73437154 0.47097142 0.84530449
0.22165260 0.78128841 0.57222522
0.24873347 0.68762843 0.48624453
0.21948415 0.65300265 -0.26571256
0.79089209 -0.20809083 -0.24704170
1.11107834 0.01381845 0.32957265
1.03239487 0.10538648 0.35890625
-0.22306466 1.43601128 1.11757892
-0.24593929 1.53904148 0.17695298
0.53330947 -0.13832394 0.74775872
0.48380538 -0.08752260 0.84569979
0.33958593 1.05428999 0.90703938
0.29513801 0.93047632 0.90464110
-0.10080201 -0.17071890 0.20726349
0.00812468 -0.20258160 0.13847835
0.08102303 0.64583254 0.03494449
-0.04079584 0.61193214 0.04209357
0.35121443 1.17006312 0.32377652
1.18017291 1.28895470 0.63748305
0.44815316 0.76165600 -0.42603898
0.46425067 0.68605942 -0.32598380
0.26710057 0.21583209 0.85066344
1.08895608 0.13400129 0.17359061
0.65188184 1.12136051 0.88104293
0.65005020 0.99400031 0.86041374
0.74276310 0.93357167 0.31617568
0.65193906 1.01050888 0.33337941
0.22155564 -0.06334224 -0.33974479
0.31584241 -0.05457389 -0.25072038
0.46701413 0.39534529 0.33741153
0.37353168 0.43618379 0.40490478
0.88534124 0.31114574 0.85608235
0.90593210 0.27739909 0.97341481
0.39857237 -0.04387750 0.56182123
0.47095275 1.02320312 0.49993905
1.25892990 0.45837652 0.17875477
1.15656687 0.39089427 0.22252886
0.51644837 0.17070272 0.40580457
1.27940834 0.60940694 0.99753791
0.42616197 0.40350440 0.06681017
1.44539498 0.30393771 -0.01376590
0.46562325 1.01534518 0.03284939
0.50319974 0.99064681 0.15470387
0.81037140 0.25082530 0.14657752
0.90914722 0.31376690 1.17885883
0.90887977 0.86774057 1.39942038
0.94187889 0.97203833 1.33441155
0.22890408 0.88765248 0.26506994
0.22943934 0.87502350 1.39740388
0.49252993 0.56873322 1.17531223
0.58322974 0.47973985 1.18073167
1.07353906 0.35801775 -0.07352972
1.10345201 0.46291908 -0.00569576
0.56233028 0.47301291 0.65701182
1.50111441 0.39245196 0.57504212
0.64733050 0.17609466 0.44047892
0.76155406 -0.00064089 0.99805907
0.11696365 0.55392898 0.64536596
0.45665944 0.88787932 -0.73431062
0.37333616 0.58173353 0.89358570
0.83287880 -0.23084103 -0.18836803
0.03157586 0.08326143 0.89500066
-0.01967168 1.44066290 -0.50120574
0.32649692 0.14675159 0.14354441
-0.47815402 -0.28476352 0.96183234
0.58825114 1.58190086 0.46693607
0.73530757 0.86486122 0.54696653
0.36850393 0.23869519 0.76098331
-0.24064441 -0.53169325 0.57977533
0.24867886 1.17176032 -0.42128557
-0.09114394 0.73086170 1.49415350
0.77082046 0.96821649 1.06941593
-0.03834292 -0.10477073 0.98185130
0.72904928 0.67702858 2.31155822
0.98643994 0.21068855 0.47648425
0.18404121 0.88355581 0.86864094
1.08703359 0.18844511 0.11444108
1.71734102 0.59616846 0.90840384
0.79123200 0.05674886 -0.33394462
0.02043659 -0.31156425 0.77570162
0.56069314 0.24338431 0.88537227
0.24569014 0.53057214 0.44182948
0.03448211 0.40131521 0.26357804
1.09819988 0.83512130 0.63612042
-0.60356192 0.67216851 0.24349363
0.68641610 -0.81762609 1.18817124
0.21611881 -0.14120573 0.12450673
0.76919341 1.22245249 0.51393968
1.14696413 1.32734585 0.69656507
0.65811040 0.77819298 1.11522722
0.04903142 0.62925399 0.32627031
0.01260756 0.05983530 0.64228866
0.74653920 0.39814097 0.82582336
0.26593357 0.76215637 0.50729285
0.25423481 0.71079817 -0.23075499
1.04656732 0.05170956 0.30685598
-0.27090311 1.50034802 1.11590884
0.55144479 -0.10997314 0.82138721
0.36595179 0.97842307 0.90355956
-0.04879022 -0.22282940 0.18797022
0.02303649 0.60374665 -0.00064438
0.36248716 1.18806896 0.40233898
0.50553342 0.72818132 -0.38174573
0.23339102 0.17685431 0.91622715
0.69519318 1.05764833 0.86986144
0.66244443 0.93542605 0.32313623
0.29157406 -0.03443908 -0.32197353
0.43576289 0.39304384 0.41288487
0.93290635 0.26807660 0.90192508
0.46947458 -0.04604768 0.53317357
1.23797107 0.38686102 0.19701795
1.24279488 0.62426866 1.06753929
1.39453671 0.36177382 0.00797544
0.49872289 1.04632858 0.09816681
0.86030136 0.31356352 0.11502915
0.87609952 0.92053833 1.35250286
0.23398309 0.92340184 1.33460122
0.50301930 0.49047944 1.17320678
1.13069333 0.39166539 -0.03166013
1.52690177 0.40275423 0.64871882
0.59258213 0.16054156 0.38717393
1.22868038 0.27296875 0.05742047
Click to reveal the
128H2O_slab POSCARSYSTEM
1
12.50000000 0.00000000 0.00000000
0.00000000 12.50000000 0.00000000
0.00000000 0.00000000 50.00000000
H O
256 128
Direct
-0.15047607 0.45334683 0.32730556
-0.07219808 -0.49010458 0.34403411
0.63805046 0.03123465 0.46999929
0.50838374 0.00754431 0.46372876
0.20742016 -0.23715390 0.31145651
0.21957350 0.63504600 0.31129793
0.49252565 0.97213243 0.52739381
0.46254946 0.86699495 0.51338367
0.46833659 0.30083696 0.47570094
0.56435815 0.28642228 0.46051829
0.73586863 -0.12693724 0.63650708
0.75181963 -0.25406949 0.63203814
0.41854353 0.84406884 0.44927074
0.33515626 0.90054023 0.46728072
0.34479509 0.57962028 0.37911763
0.29762780 -0.47903310 0.40660055
0.15115248 0.54108257 0.49957979
0.08260965 0.53074503 0.52635427
0.89266849 0.40696421 0.52478091
-0.06057907 0.40486920 0.55352216
0.55977900 -0.37296647 0.61748251
0.53858714 -0.26767121 0.59987163
0.71275400 0.42250550 0.56832405
0.78839559 0.34546846 0.58417890
0.67385329 -0.19610180 0.44987464
0.59318198 -0.16624091 0.42828513
0.12312508 0.37123899 0.28937208
0.05023666 0.27675971 0.29997236
-0.16681739 1.58406768 0.63195176
-0.21800579 1.57002107 0.66191007
0.20675683 -0.06108281 0.29221359
0.10312049 -0.09294129 0.30396903
0.20277637 0.06067826 0.60416144
0.19887396 -0.00744650 0.62890026
0.42062864 -0.24439336 0.35497596
0.36814211 -0.33879224 0.34193308
-0.46114333 0.32802380 0.39548517
0.45713828 0.41873416 0.39179412
1.05751346 0.29067408 0.55923161
1.04348507 0.18703047 0.57744764
0.21179010 0.87249263 0.56793233
0.09036024 0.85931851 0.57432053
0.60751605 1.01609789 0.65259722
0.72696147 1.05338974 0.65871807
0.44683747 0.21956247 0.29940199
0.48216169 0.19981749 0.26927194
0.10930269 0.40445023 0.47082959
0.17010928 0.50009382 0.45569496
0.78757796 -0.20083679 0.40459572
0.81784835 -0.20993279 0.37400609
0.63044955 -0.50216422 0.32507910
0.61823261 -0.56783464 0.29754461
0.38298254 1.64224442 0.61713447
0.27391580 1.65127025 0.62769264
-0.18349257 0.40823838 0.42743175
-0.11707568 0.32822140 0.44206697
0.37713689 0.15324147 0.36153587
0.47188832 0.15561964 0.38076476
0.84313290 0.66869553 0.54540293
0.78036268 0.55967659 0.53781156
0.62834565 0.07729786 0.39976222
0.58336683 0.13133194 0.42359675
1.57461032 1.43827952 0.71246905
1.47093630 1.44260924 0.69550961
0.56330564 1.09120276 0.50193580
0.61593811 1.12895503 0.52786251
0.73348421 0.31660428 0.46194097
0.75096480 0.20702765 0.47568711
0.18544389 -0.23514793 0.40166768
0.22931162 -0.14674533 0.42153910
0.36292308 0.50876490 0.64987458
0.33531338 0.38470088 0.66057777
1.49747974 -0.54677440 0.44932443
0.51766749 -0.42535318 0.46152212
1.85759374 0.75185320 0.76879203
0.88557936 0.66337125 0.74776712
0.66151301 0.68965700 0.41313588
0.60698800 0.57649179 0.41787791
-0.21783861 0.90704975 0.72268493
-0.18062377 1.02010056 0.73108294
0.79398077 0.91221618 0.53177815
0.89172292 0.88245064 0.51170593
1.21360535 1.19956259 0.58314768
1.23188220 0.10388812 0.56520461
0.69217913 -0.02288592 0.35644459
0.65013288 -0.02349422 0.32542268
0.49314736 0.09990072 0.55428120
0.40698133 0.13735440 0.57630192
0.06940342 0.10352406 0.53021711
0.12918295 0.01423295 0.54673459
1.23761207 0.03882752 0.37085803
0.23384130 0.05388881 0.33858930
-0.50295222 0.92300366 0.74813821
-0.38930037 0.97841032 0.74266531
0.13655853 0.16157403 0.39040653
0.15373516 0.09867462 0.41870070
0.35793276 0.15462364 0.64561652
0.44313351 0.08542197 0.66008140
0.12490452 1.60323574 0.67914603
0.10103912 0.67758518 0.66052594
-0.02910081 -0.04303719 0.33290797
-0.08397074 -0.06715582 0.30625602
0.47100963 -0.31135459 0.42347049
0.44678381 -0.21727665 0.40309666
0.17380547 1.23522872 0.64351497
0.28087835 1.29207879 0.63186950
0.72243965 1.04679669 0.60136565
0.77633776 1.15885707 0.60968516
0.20332595 1.15182338 0.29172154
0.31769260 1.09822536 0.28860721
0.90406410 0.53849761 0.41088426
0.78463589 0.58505950 0.41470974
1.12942702 1.53845314 0.71988868
1.20054880 1.43235266 0.71382259
-0.17158886 0.30160650 0.40269399
-0.07101085 0.26515881 0.38366272
0.74322462 1.32613449 0.64589263
0.77475704 1.41411198 0.62227898
0.14620565 -1.10616729 0.48008239
0.18160223 -0.03968360 0.50606301
-0.11186987 0.94149530 0.65599853
-0.13718975 1.02669042 0.63199678
0.56159298 0.19900069 0.59311202
0.55920469 0.18285975 0.62176711
0.90336505 0.15667383 0.52548512
0.98450787 0.25110265 0.52511525
0.97847305 -0.42779836 0.72340892
1.02887897 0.67051183 0.71108441
0.93053153 -0.03352705 0.43016463
0.94358325 0.08994529 0.42307038
-0.12458051 0.56502080 0.45311411
-0.04730932 0.59270705 0.47590220
-0.26829199 0.80830572 0.49363371
-0.21717255 0.70457633 0.48040899
0.43110180 0.68638080 0.68890088
0.32510025 0.62119216 0.70042335
0.12431384 0.78282478 0.61843684
0.20700484 0.81000752 0.64341301
0.70384786 0.11901625 0.29636532
0.58157041 0.12693352 0.30130944
0.75363478 0.23168055 0.54459190
0.75644483 0.12259613 0.56105372
-0.03231861 0.12955620 0.32499376
-0.09853097 0.23584383 0.32218754
1.35833840 -0.44225252 0.44591484
1.25776038 -0.35606860 0.43991425
0.02974229 -0.01448169 0.38634402
-0.04925278 -0.12297578 0.38637941
0.07972722 -0.25141011 0.36219086
0.18306699 -0.31067445 0.35551374
0.92747522 0.76466556 0.60868816
0.95937680 0.87571382 0.62123825
0.52932163 0.57770897 0.50403614
0.58437494 0.68163282 0.49342216
0.32115740 0.82996586 0.67667456
0.25755504 0.94154338 0.67730506
-0.22810590 0.37614398 0.49430267
-0.18576922 0.49671288 0.49149878
0.65189321 0.70221781 0.55887514
0.58469375 0.80139057 0.54950481
0.45751382 1.14320451 0.73814919
0.45757048 1.25211787 0.71984711
-0.16073265 0.53318711 0.58413419
-0.04460020 0.52255466 0.59636621
1.23763588 0.27151415 0.53650362
0.26869350 1.37267158 0.55194387
1.56531541 1.75888578 0.65789015
1.64170785 1.76839833 0.68268336
1.17071819 0.30450968 0.36193986
1.07205789 0.24647217 0.34988931
0.66268682 0.39143410 0.35576888
0.77470013 0.33753052 0.36083320
0.48701966 -0.62651414 0.32289219
0.45172427 -0.71444999 0.34206563
0.89274414 0.87693241 0.57141470
0.77760242 0.86753703 0.58449204
1.22348599 -0.38621798 0.54275880
1.26314524 -0.31955173 0.56641379
0.93813763 1.18794820 0.71316474
0.91809214 1.10192472 0.69334364
0.61064405 0.92402610 0.60719628
0.55192305 1.01039312 0.59021354
0.96908959 0.73238970 0.67284455
0.86559314 0.75071650 0.65769890
1.05951362 -0.45354408 0.38093143
1.05009801 -0.58200279 0.37590291
0.76172957 0.49270651 0.70861525
0.72438777 0.60310359 0.69947487
0.43493202 1.01642204 0.62370834
0.41157215 0.90717107 0.63919892
1.28709531 1.22977646 0.71568664
1.22240102 1.24760104 0.69023856
-0.36929002 0.30112790 0.68799301
0.73835598 1.24904006 0.69307962
1.08438837 1.16887216 0.75281234
1.10337181 1.25880709 0.73309939
0.00881530 0.23511529 0.66153556
0.01992074 0.10863332 0.65331960
-0.18762204 1.76544824 0.69589571
-0.22099427 0.73159748 0.72364889
0.47402075 0.21823941 0.52156939
0.39449064 0.32549396 0.52603852
0.49229702 0.00642229 0.70437200
0.49041094 -0.09636789 0.68559144
0.31519106 1.45346432 0.61090090
0.38603816 1.37161975 0.59824648
0.15729732 1.50568246 0.63546369
0.21559092 1.47629782 0.66066682
-0.19454501 1.01393222 0.39591120
-0.22080313 1.12719611 0.38222917
1.06810038 0.81254123 0.53206673
1.06352248 0.70747723 0.51307936
-0.12638475 -0.55537021 0.26026243
-0.05760421 -0.65321822 0.26577609
-0.09042409 1.07196610 0.58294900
0.02325129 1.03828409 0.59636367
0.32063486 0.92608540 0.72659258
0.35574911 0.80830229 0.72123409
0.10412373 0.41841079 0.59896831
0.14274925 0.51651990 0.58178437
-0.19591148 0.07044928 0.45554844
-0.19021572 0.03340561 0.48590483
1.07785299 0.53236279 0.30851784
1.17527108 0.48273981 0.32369409
-0.05242302 0.21986278 0.47826685
0.04602452 0.20858933 0.45856538
0.40878055 0.29190606 0.42949757
0.35385287 0.26786429 0.40109729
0.91361281 1.37763527 0.70485166
0.92601382 1.37974445 0.73573816
0.49180218 0.49176486 0.55630311
0.38888954 0.54455977 0.56839724
0.74382564 0.28498135 0.31337726
0.81303180 0.29929615 0.28888582
0.23876821 1.20570242 0.44319221
0.19992227 1.11701411 0.46333591
0.59407806 1.36278641 0.60873199
0.54524237 1.46654437 0.59522322
2.02346307 1.41217476 0.66399254
1.91147497 1.37915850 0.65989115
0.09845777 0.03688460 0.69616599
0.09766268 -0.08531166 0.69167377
1.25808011 0.44420386 0.36149889
1.31081579 0.38358222 0.33934553
0.07353811 0.81613159 0.44176424
-0.02019300 0.78679404 0.46293839
0.38984424 0.83369934 0.56338610
0.37381688 0.87922811 0.59234427
0.51453372 -0.09541985 0.35009748
0.39484305 -0.06645979 0.35017125
0.55593924 0.38753163 0.51869390
0.66004756 0.46315941 0.51627573
0.64848801 0.61180544 0.37008072
0.57832658 0.64084106 0.34599612
0.75386845 -0.15146447 0.33086236
0.73564205 -0.27427126 0.33331691
-0.09648741 0.50938984 0.32488838
0.56504329 0.06286545 0.46968741
0.24293994 0.69761578 0.32091286
0.48108748 0.89535134 0.53166188
0.49757936 0.32553516 0.45789152
0.76744020 -0.18340073 0.62487691
0.41144876 0.90246316 0.46247378
0.31794376 -0.48573743 0.38681989
0.12877684 0.58271134 0.51624414
0.95851509 0.39985886 0.53493477
0.52243754 -0.30282985 0.61791897
0.76895759 0.37346632 0.56560940
0.66848522 -0.18560993 0.43043112
0.09074693 0.30162769 0.28441615
-0.20154125 1.61789646 0.64720226
0.17181260 -0.12468412 0.29662573
0.15700557 0.00871911 0.61363921
0.42233976 -0.32994901 0.35514866
0.53286785 0.40322511 0.39184579
1.07727341 0.25676911 0.57497686
0.13716955 0.87720865 0.55959520
0.68216863 0.99627546 0.65137040
0.45119708 0.16614064 0.28412456
0.10973328 0.48171438 0.46781222
0.84388934 -0.20664059 0.39174052
0.58303852 -0.54048463 0.31286474
0.32808472 1.59763404 0.62452336
-0.19165309 0.33704380 0.43457786
0.41823919 0.20183034 0.37245093
0.78587144 0.61673307 0.55124293
0.56011355 0.11309492 0.40567704
1.51701295 1.39737040 0.70673631
0.54695764 1.10711476 0.52102138
0.72131688 0.27947218 0.47879616
0.18444949 -0.21210987 0.42113103
0.36668171 0.45488107 0.66482426
0.49064798 -0.46503727 0.44561439
0.82820218 0.69414688 0.75657380
0.65847198 0.61537373 0.40712496
0.76748450 0.96128069 0.73583267
0.81730697 0.91059549 0.51275491
1.27220645 1.14795719 0.57873009
0.64741928 -0.05620926 0.34353694
0.48391157 0.12047094 0.57316210
0.14376608 0.06922992 0.53329165
1.27850839 0.03499489 0.35317314
-0.46354066 0.99485951 0.74685190
0.12508520 0.09293059 0.39933758
0.43296574 0.12732681 0.64290755
0.09700574 1.67696831 0.67957377
-0.02662689 -0.02348385 0.31412699
0.44363826 -0.24009345 0.42143376
0.25025628 1.24855748 0.64562073
0.79823150 1.08209822 0.60445310
0.24357641 1.08550799 0.29686686
0.84487510 0.54791453 0.42274978
1.20665739 1.51259035 0.71454551
-0.14852600 0.27439958 0.38483199
0.79875463 1.36570466 0.63635127
0.20013661 -0.05770220 0.48695782
-0.08803369 1.01243312 0.64786941
0.60681036 0.21343304 0.60836199
0.96760512 0.18267756 0.51744002
1.01869709 -0.36164031 0.72861865
0.89279549 0.03025280 0.42811891
-0.12190623 0.58173239 0.47257350
-0.28364610 0.74658147 0.48027041
0.35461305 0.69479199 0.69387913
0.17636099 0.75539076 0.63132326
0.64567772 0.08234767 0.30555919
0.74693252 0.15345949 0.54261621
-0.02750156 0.20761939 0.32607759
1.27987775 -0.43232854 0.43958099
0.00284833 -0.07336365 0.37560644
0.13403374 -0.30077764 0.37038742
0.98462893 0.81761037 0.60965656
0.51924082 0.63814697 0.49321682
0.28461449 0.88126714 0.66475622
-0.20914894 0.44302019 0.50652580
0.58914466 0.74769132 0.56398020
0.41471483 1.19061656 0.72530883
-0.12231266 0.52176874 0.60092233
0.26021422 1.34539425 0.53435665
1.56662855 1.77066279 0.67736079
1.10461724 0.26953014 0.36735832
0.72369042 0.35764108 0.34663600
0.43147580 -0.68083095 0.32528778
0.81926513 0.85574301 0.56746752
1.25134221 -0.39325035 0.56126973
0.87947711 1.14636757 0.70571444
0.62353030 0.97359658 0.59258366
0.91315731 0.79159814 0.67046472
1.00887481 -0.51613717 0.38007575
0.71990656 0.52980508 0.69520722
0.42715157 0.93815001 0.61949356
1.22065693 1.25957506 0.70894987
0.68302181 0.25679220 0.67873356
1.04556333 1.22974551 0.74538666
0.05950376 0.17484591 0.65818294
-0.25231321 1.76973735 0.70816158
0.45404841 0.29350534 0.51621152
0.47580393 -0.01885257 0.68583286
0.31607039 1.37827042 0.60641490
0.13921945 1.47529964 0.65395905
-0.24949740 1.04980313 0.38476573
1.02198254 0.77527198 0.51809412
-0.13269955 -0.63172943 0.26362895
-0.01665243 1.05122555 0.57819040
0.34830926 0.86183449 0.73481007
0.08547124 0.49199314 0.59517301
-0.22931968 0.07879046 0.47296684
1.15374394 0.51400513 0.30712860
-0.02575549 0.24426930 0.46122795
0.33910288 0.28875076 0.42031444
0.89222545 1.42242644 0.72114703
0.45112644 0.49752204 0.57217511
0.78732813 0.24777689 0.30076703
0.17546728 1.17137686 0.45167807
0.60114018 1.43939441 0.60809185
1.96388141 1.37112395 0.67251971
0.13553192 -0.02618932 0.70064801
1.23640886 0.39518318 0.34779401
0.01704024 0.84626979 0.45444055
0.34009588 0.83637329 0.57759905
0.44929940 -0.11284502 0.35947285
0.58008850 0.46021347 0.52387715
0.65101061 0.61092524 0.34954978
0.79285789 -0.22059881 0.33571374
INCAR
The INCAR file differs from step to step. Each is given on the how-to page that the corresponding step links to.
KPOINTS
The KPOINTS file contains only the Γ point:
Gamma-point only 0 Monkhorst Pack 1 1 1 0 0 0
POTCAR
Standard POTCAR files are used throughout:
PAW_PBE H 15Jun2001PAW_PBE O 08Apr2002PAW_PBE Fe_sv 23Jul2007
This is different to in Ref. [1], where GW POTCARs are used.
Step-by-step instructions
Step 1: Preparing the MLFFs
First, train the MLFFs for the water slab (128H2O_slab) and the redox pairs ([math]\displaystyle{ [\mathrm{Fe}(\mathrm{H}_2\mathrm{O})_n]^{3+} }[/math] and [math]\displaystyle{ [\mathrm{Fe}(\mathrm{H}_2\mathrm{O})_n]^{2+} }[/math]).
Step 2: Calculating the vacuum potential
Run a molecular dynamics (MD) simulation, then calculate the local potential [math]\displaystyle{ \phi }[/math] of the water slab. Using this, the difference between the Oxygen 1s level between the water slab and the water molecules far away from the ions in the redox cells can be calculated. This provides the chemical potential relative to the vacuum [math]\displaystyle{ e \Delta \bar{\phi} }[/math].
| Important: This step requires using GGAs on many structures, so will be more time-consuming. |
Step 3: TI from oxidized to reduced state
Perform a thermodynamic integration from the oxidized state (λ=0; Fe3+) to the reduced state (λ=1; Fe2+).
| Mind: This is a thermodynamic integration between two MLFFs, so is very quick. |
Step 4: TI from MLFF to GGA
Perform a thermodynamic integration from MLFF (λ=0) to RPBE+D3 (λ=1). This should be done separately for Fe3+ and Fe2+.
| Important: This step requires using GGAs on many structures so will be more time consuming. |
| Mind: Thermodynamic perturbation theory (TPT) from RPBE+D3 to PBE0, for both Fe3+ and Fe2+, is not covered here because of its cost. It can be replicated using Δ-learning from RPBE to PBE0, using the same procedure as from MLFF to RPBE. |
Step 5: Analysis
Once each step is completed, you will have the following values:
- the chemical potential - [math]\displaystyle{ e \Delta \bar{\phi} }[/math]
- the free energy difference from MLFF to MLFF (with chemical potential reference)
- from Fe3+ to Fe2+ - [math]\displaystyle{ \Delta A_{\mathrm{ML}} }[/math]
- the free energy difference from MLFF to GGA
- Reduced: Fe2+ - [math]\displaystyle{ \Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_1 }[/math]
- Oxidized: Fe3+ - [math]\displaystyle{ \Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_0 }[/math]
- [math]\displaystyle{ \Delta A^{\mathrm{FP}_\mathrm{GGA}} = \Delta A ^{\mathrm{ML}} + \Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_1 - \Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_0 }[/math]
The vacuum alignment is already contained in [math]\displaystyle{ \Delta A^{\mathrm{ML}} }[/math]. Combining the remaining terms as above gives the free energy difference from Fe3+ to Fe2+ on the GGA level [math]\displaystyle{ \Delta A^{\mathrm{FP}_\mathrm{GGA}} }[/math]. Switch the sign to get the absolute redox potential [math]\displaystyle{ U^{\mathrm{FP}_\mathrm{GGA}}_{\mathrm{redox}} }[/math]. Finally, take the standard hydrogen electrode (SHE) (4.44 V [2]) away to compare to experiment [math]\displaystyle{ U^{\mathrm{FP}_\mathrm{GGA}}_{\mathrm{redox,SHE}} }[/math]. The results below match the literature to within 30 meV.
The results are summarized in the following table:
| Step | Term | Ours (eV) | Lit. [1] (eV) |
| 2 - Vacuum reference | [math]\displaystyle{ e \Delta \bar{\phi} }[/math] | -3.71 | -3.74 |
| 3 - TI (MLFF) | [math]\displaystyle{ \Delta A_{\mathrm{ML}} }[/math] | -4.98 | -4.95 |
| 4 - TI (GGA) | [math]\displaystyle{ \Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_1 - \Delta A^{\mathrm{FP}_\mathrm{GGA}-\mathrm{ML}}_0 }[/math] | -0.055 | -0.049 |
| Absolute redox potential | [math]\displaystyle{ U^{\mathrm{FP}_\mathrm{GGA}}_{\mathrm{redox}} }[/math] | 5.03 | 5.00 |
| Redox potential | [math]\displaystyle{ U^{\mathrm{FP}_\mathrm{GGA}}_{\mathrm{redox,SHE}} }[/math] | 0.59 | 0.56 |
Recommendations and advice
- The vacuum-reference step and the MLFF:GGA integrations are by far the most expensive parts of the workflow, since both evaluate a GGA on many structures. Training the MLFFs in the first place can also be time consuming. The MLFF:MLFF integration, by contrast, is very quick.
Related tags and articles
- How-tos
- Vacuum reference
- Thermodynamic integration between machine-learned force fields
- Thermodynamic integration between a machine-learned force field and a density functional
- Theory
- Tags