Oxidation of water on a Pt(111) surface¶

Bulk Pt¶

We begin by setting up a geometry optimization job to build and optimize bulk Pt(111).

To create the Pt lattice structure

Import geometry from xyz fileBuild a new structure

1. Download the xyz file of bulk Pt

2. In AMSinput select File → Import Coordinates and select the downloaded file

For the geometry optimization we will use the revPBE-D4(EEQ) XC-functional:

In the Main panel

1. Select XC functional → GGA-D → revPBE-D4(EEQ)

2. Select Task → Geometry Optimization

3. Set Relativity (ZORA) to Scalar

4. Click the button next to Task or in the panel bar open Details → Geometry Optimization

5. Check the Yes checkbox next to Optimize lattice

We can now run the calculation using File → Run. This will prompt you to save the job somewhere if you have not saved it already. The progress can be viewed in the log file viewer (SCM → Logfile)

Once the calculation has finished, a pop-up will appear asking you to read in the new coordinates. Click Yes. Under Model → Lattice you can find the new lattice constant. Check that it is close to the literature value of (3.92 Å). You can try to see for yourself that omission of relativistic effects will cause the lattice to expand quite significantly. For more accurate relativistic effects, one could employ spin-orbit coupling effects as well, however these are more computationally demanding. You can visualize the band structure and density of states (DOS) by clicking on SCM → Band Structure and SCM → DOS respectively while in the output viewer.

OH and H₂O on a Pt(111) surface¶

Next we cut a (111) surface from the optimized bulk structure.

Import geometry from xyz fileBuild the Pt slab yourself

1. Download the xyz file of the Pt(111) surface

2. In AMSinput select File → Import Coordinates and select the downloaded file

Save this as a new file and run a single point calculation (change the Task from Geometry Optimization to Single Point) using the same settings as before.

Now we will put one monolayer (ML) of OH on the surface in the top position and set up the run.

To set up the geometry:

Import geometry from xyz fileBuild the OH monolayer yourself

1. Download the xyz file of the OH monolayer on a Pt(111) surface

2. In AMSinput select File → Import Coordinates and select the downloaded file

And to set up the calculation:

1. Go to the Main panel and select Geometry Optimization as Task

2. Click the button next to Task or in the panel bar open Details → Geometry Optimization

3. Deselect the Yes checkbox next to Optimize lattice

4. Select the Pt atoms by clicking them while holding Shift

5. In the panel bar, select Model → Geometry Constraints and PES Scan and press the top (fixed position) button to freeze the atoms in place

Run the job. It will take a few cycles to finish. In the meantime we will prepare the optimization for a monolayer of water on Pt(111).

Import geometry from xyz fileBuild the H₂O monolayer yourself

1. Download the xyz file of the H₂O monolayer on a Pt(111) surface

2. In AMSinput select File → Import Coordinates and select the downloaded file

Note that neither the 1 ML H₂O nor the OH geometries are physically (chemically) meaningful. Water will typically adopt a hexagonal bilayer structure. Check out the references at the end of the tutorials for a starting point on more realistic geometries.

Save as a different file and run the job. The OH calculation should be done by now, but the H₂O calculation will take slightly longer than the OH one. In the meantime, we will prepare calculations on free H₂, H₂O and OH, so that we can calculate their thermodynamic properties using ADF.

H₂, H₂O and OH molecules¶

We now continue with calculations on free H2, H2O and OH molecules. The general workflow is as follows:

1. Open a new AMSinput window and ensure that you are in the panel

2. Click the button and click anywhere in the molecule drawing area to place a hydrogen atom

3. With the first hydrogen selected, press anywhere else to create H₂

4. Preoptimize the molecule by pressing the button at the bottom of the molecule drawing area. By right-clicking you may choose a different engine for preoptimization

5. In the Main Panel, select XC-functional → GGA-D → revPBE-D4(EEQ)

6. Select Relativity → Scalar

7. Set Task to Geometry Optimization and check the Frequencies checkbox below Task

Selecting Frequencies will ensure that vibrational and thermodynamic properties of the systems are also calculated.

Run the job and repeat for H₂O and OH by changing the structure of the molecule. For the OH calculation, please note that OH is a radical (doublet) species and you should ensure that you select the Unrestricted checkbox in the Main panel and enter 1 in the Spin polarization field.

Estimated water oxidation potential on Pt(111) surface from thermodynamical properties¶

We now have all the pieces available to estimate the minimum potential at which water can get oxidized over a Pt(111) surface:

where the starred species denote adsorbed molecules on the Pt(111) surface. Following Nørskov et al. (J. Phys. Chem. B, 108, 17886-17892 (2004)) we choose the standard hydrogen electrode as a reference potential, such that

is in equilibrium. We get the overall reaction

The free energy of the reaction is then:

The electrode potential \(U\) is incorporated as a linear energy shift \(qU = -U \text{eV}\) for every electron, giving

As a first approximation for the Gibbs free energies of OH* and H₂O* we will just use the formation energies (from the atomic fragments). By far the largest contribution to the (difference in the) thermodynamic properties of the adsorbed species will be the zero-point energies (ZPEs), which can be crudely estimated from the gas phase frequency calculations for OH and H₂O.

From the log file of the OH and H₂O on Pt(111), write down the energy of formation near the end. (If necessary, look in the list of jobs in the AMSJobs window and click on the triangle in front of the job name so that all related files are listed. Double click on the .logfile files to open the appropriate logfile.)

from the H₂ geometry optimization:

For the rotation and vibration contributions to \(\Delta G(\text{H}_2)\), we examine the more detailed output file (.out or Output from the SCM menu). The menu bars for the output viewer enable quick jumping through the calculation output. In this case we want the Gibbs free energy from Statistical Thermal Analysis which is located in the Other Properties menu.

and thus:

Consequently, the reaction over a Pt(111) layer will only take place when \(\Delta G < 0\). This is only true when we set the potential to at least \(U = 0.1494 \text{ V}\).

The calculated value is quite far from the experimental value of 0.8 V (Markovic and Ross, Surf. Sci. Rep. 45, 117 (2002)). Will including ZPE effects improve things? No: look in the output of the frequency calculations Other Properties → Zero-Point Energy:

so that the estimated \(\Delta\text{ZPE} = -0.3265 \text{ eV}\) for the gas phase species. Therefore, including \(\text{ZPE}\) would make it easier to oxidize water over Pt(111), reducing the oxidation potential. However, if you have some time left at the end of this tutorial, you can convince yourself that lower coverages will improve agreement with experiment:

Supercells for lower adsorbate coverage¶

Alternatively, you could have started from the optimized OH* and H₂O* 1 ML structures, making supercells from these and then removing two adsorbate molecules in the super cell. In this way we create a cell which has a third of the coverage of the 1ML layer.

These calculations should take about 2 and 6 hours to complete, so just run them overnight and read back the last formation energies from the logfile when you return tomorrow. These values (-77.9892 and -82.2107 eV for OH and H₂O respectively) predict a much better water oxidation potential of about 0.8685 eV (with the same estimated gas phase \(\Delta \text{ZPE}\)).

Note¶

that this tutorial was just a short introduction to illustrate the calculation of thermodynamical properties of electrochemical half-reactions. In general for these kind of calculations you would want to converge carefully with respect to basis set, k-point sampling, and slab thickness. One should also carefully consider the adsorption geometry, e.g. water will want to form hexagon overlayers (necessitating √3 x √3 supercells).

All this goes too far to cover in a short tutorial session, but there are plenty of good literature examples where they have carefully considered these and more subtleties involved in these kind of calculations. Useful starting point are the papers by Nørskov and coworkers, for instance addressing oxygen reduction (2004) and water redox reactions on Pt(111) (2006).

BAND can also perform frequency calculations (more computationally demanding than ADF), so that ZPE effects can also be included for the adsorbed species rather than estimated from the gas phase. Finally, note that in BAND, uniquely, both solvation effects (COSMO) and homogeneous electric fields can be included properly (true 2D periodic systems). Both these effects will influence the adsorbate energies, and thereby the calculated oxidation potential.