Dissociation Mechanism of Water Molecules on the PuO2(110

Dec 20, 2017 - ... American Chemical Society. *E-mail: [email protected], *E-mail: [email protected]. .... Full Text HTML ... Email a Colleague...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIV OF NEW ENGLAND ARMIDALE

Article 2

Dissociation Mechanism of Water Molecules on the PuO (110) Surface: An Ab Intio Molecular Dynamics Study Cui Zhang, Yu Yang, and Ping Zhang

J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.7b08864 • Publication Date (Web): 20 Dec 2017 Downloaded from http://pubs.acs.org on December 21, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Dissociation Mechanism of Water Molecules on the PuO2 (110) Surface: An Ab Initio Molecular Dynamics Study Cui Zhang,† Yu Yang,∗,† and Ping Zhang∗,† Institute of Applied Physics and Computational Mathematics, PO Box 8009, Beijing 100088, China, and Center for Applied Physics and Technology, Peking University, Beijing 100871,China E-mail: [email protected]; [email protected]

∗ To

whom correspondence should be addressed of Applied Physics and Computational Mathematics, PO Box 8009, Beijing 100088, China ‡ Center for Applied Physics and Technology, Peking University, Beijing 100871,China † Institute

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract Reactions between water and plutonium oxide play a key role in determining the oxidation and corrosion rates of plutonium materials. We perform consecutive ab initio molecular dynamics simulations with the DFT + U approach, systematically studying the dissociation dynamics and mechanism of water molecule and small clusters on the PuO2 (110) surface. The dissociation of water on the surface is found to be a two-step hydroxylation process for both monomer and clusters, but different dissociation mechanisms are revealed. Water monomer dissociation, as the consequence of the hybridizations between the molecular orbitals of water and the electronic state of the surface, needs to overcome a reaction energy barrier of 0.18 eV. In contrast, the dissociations of hydrogen-bonded small water clusters are exothermic by -0.42 eV with no transition energy barrier. Hydrogen bonding interactions between water molecules facilitate the dissociation process and result in a more stable dissociation state with respect to that of single molecule. Such interaction also induces a hydrogen transfer reaction between formed surface hydroxyl group and remaining water molecule.

2

ACS Paragon Plus Environment

Page 2 of 19

Page 3 of 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

1 Introduction Plutonium chemistry is attracting much interest for its important role in military weapon security and storage. When exposed to air and moisture, the plutonium metal surface rapidly oxidizes and corrodes, forming a protective layer of dioxide. At room temperature, the corrosion rate of plutonium by water at its equilibrium vapor pressure is more than 100 times faster than that in dry air. Such moisture enhancement is increased to be about 105 times at boiling point. 1 To understand the dynamics and mechanism of such moisture-enhanced oxidation process is an essential scientific issue both in fundamental actinide science and nuclear energy application. Experimental studies of water interactions with PuO2 are very challenging. 2–9 In early work of Stakebak, 2 two distinct water desorption temperature range were found, one between 373−423 K and the other between 573−623 K, which suggested a strong chemisorbed first water layer and a weak physisorbed second layer. The higher temperature desorption was considered to be due to dissociative first layer water adsorption with an estimated adsorption energy of -2.94 eV, while the lower temperature was attributed to molecular second layer adsorption with an estimated adsorption energy of -0.88 eV. Paffet et al. 7 confirmed the process and revised the estimations, reporting adsorption energy values of first layer dissociative adsorption and second layer molecular adsorption at 371 K to be -1.82 and -1.11 eV, respectively. Reaction of water with plutonium oxide was also proposed to lead to the formation of higher oxide, a PuO2+x phase. 3,5,6,10 However, recent photoemission work excluded the stability of higher plutonium oxide as a bulk system by showing that PuO2 is only chemisorbed by an oxygen layer, which desorbs at elevated temperature. 8 Many problems involving PuO2 /H2 O interactions remain to be further explored. Plutonium oxides are systems involving strong electronic correlations driven by the 5 f electrons. 11–13 Standard exchange-correlation approximations in density functional theory (DFT), such as the Local Density Approximation (LDA) or the Generalized Gradient Approximation (GGA), fail to predict correct electronic structure and magnetic order of many actinide oxides, due to underestimating the strong on-site Coulomb repulsion of 5 f electrons. The DFT + U approach is one of the most popular methods allowing to correctly treat systems containing strongly correlated 3

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

electrons. By adding an energy term of the effective Hubbard parameter, the antiferromagnetic Mott insulator feature and structural parameters of PuO2 are well described. 13–22 Several previous studies have been focused on the properties of plutonium-oxide surfaces. 23–26 Jomard and Bottin investigated the influence of electronic correlations on thermodynamic stability of PuO2 surfaces in ref. 24 and also pointed out that the electronic structure of the insulating PuO2 surfaces can be well described within the DFT + U method, which is a prerequisite to explore chemical reactions on the PuO2 surfaces. Sun et al. further reported the effects of thickness and oxygen vacancy on the stability and chemical activity of the PuO2 surfaces, where the spin-orbit coupling (SOC) interaction was considered. Recently, the stabilities of clean and hydroxylated PuO2 surfaces were compared and discussed in ref. 26 . The first theoretical study on the interaction between water molecule and the PuO2 surfaces was carried out by Wu et al. within the DFT framework, where both the strong correlation effect of 5 f electrons in plutonium and the van der Waals interactions between the molecule and surfaces were not considered. 27 More recently, Jomard et al. applied the DFT + U approach to study the molecular and dissociative adsorptions of water molecules on the PuO2 (110) surface. 28 Adsorption energy were found to be -0.87 and -1.12 eV for the molecular and dissociative configurations, respectively, with a small energy barrier of 0.05 eV between two states. In ref., 26 a dramatic effect of the dissociated water on the relative stability of the PuO2 surfaces was reported. Tegner et al. presented results of water molecular and dissociative adsorptions at monolayer coverage on the PuO2 for all three surfaces. 29 They found a mixed molecular and dissociative water adsorption to be most stable on the (111) surface, whereas the fully dissociative configuration is most stable on the (110) and (100) surfaces. Wellington et al. reported a study of water adsorption on the PuO2 (110) and (111) surfaces, using hybrid DFT within the periodic electrostatic embedded cluster method. 30 Similar energies were found for water molecular and dissociative adsorption on the (111) surface, while a clear preference for dissociative adsorption was shown on the (110) surface. Most previous theoretical studies on the water interactions with the PuO2 surfaces have been focused on the adsorption geometries and energies. The dynamics and mechanism of water

4

ACS Paragon Plus Environment

Page 4 of 19

Page 5 of 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

dissociation on the PuO2 surfaces, which is a key factor in the corrosion and oxidation reactions on the plutonium surfaces, still remain unresolved. In this work, we perform ab initio molecular dynamic (AIMD) simulations, together with DFTbased electronic structure calculations, to investigate the dynamics and reactivity of water molecule and small clusters on the PuO2 (110) surface. The rest of this paper is organized as follows. Our methodological approach is presented in section 2, and our results are discussed in section 3. Section 4 contains our main conclusions.

2 Computational Methodology All DFT calculations were performed using Vienna ab initio simulation package (VASP), 31 a plane-wave DFT code using projector augmented wave (PAW) 32,33 potential to describe the electronion interactions. The PuO2 (110) surface model used in this work consists of six atomic monolayers, and a periodic (2×1) supercell with a vacuum spacing of 15 Å along surface normal direction (z) was adopted. The gradient corrected functional PBE, 34,35 with a Hubbard U correction for the 5 f electrons of plutonium, was used for the exchange correlation interactions. The Coulomb energy U and exchange energy J were selected with the following values: U = 4.75 eV and J = 0.75 eV, which were tested and used in previous work. 25 We adopted optPBE exchange functional for the van der Waals correlation to describe the dispersion forces within the adsorption system. 36–39 optPBE is an optimized PBE-style functional, yielding a low mean absolute deviation from the CCSD(T) results on the S22 set (a set of 22 weakly interacting dimers) and precise descriptions for the water hexamers. 38,39 Spin-polarization was included in all calculations. The calculated lattice parameter of bulk PuO2 crystal is 5.35 Å, in good agreement with the experimental value of 5.398 Å. 6 The atomic positions were allowed to relax until the total energy change was less than 10−4 eV per unit cell and all components of forces were smaller than 0.01 eV/Å. A plane-wave cutoff of 680 eV was used in all structural optimizations. The Brillouin zone integration was calculated

5

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 19

using 5×7×1 Monkhorst-Pack k-point meshes. The adsorption energy of water molecule is given by Ead = E(slab+mol) − Eslab − Emol ,

(1)

where Eslab+mol and Eslab are the energies of the slab with the adsorbed molecule and the clean slab, respectively; Emol is the total energy of the water molecule in the gas phase, calculated in a cubic cell with a dimension of 20 Å. Spin-orbit coupling was considered in the water adsorption energy calculations. All AIMD simulations were carried out in the NVT ensemble within a Born-Oppenheimer framework at 300 K. A plane-wave cutoff of 400 eV and a time step of 1 fs were used. Only Γ point was adopted to sample k space.

3 Results and Discussions We first investigate two types of water adsorption, molecular and dissociative, on the PuO2 (110) surface. Optimized structures of two types of adsorption were shown in Figure 1. Molecular adsorption of single water molecule occurs with the hydrogen atoms tilting towards the surface and the oxygen atom lying in the position above surface oxygen empty sites, as demonstrated in Figure 1(a). The Pu· · ·Ow (Ow denotes the oxygen atom in the adsorbed water molecule) distance for water molecularly adsorbing on the (110) surface is 2.61 Å. For dissociative adsorption, the adsorbed hydroxide (Ow H) locates right above a surface Pu atom with its hydrogen pointing against surface, while the other hydrogen from water is bonded with a surface oxygen atom and angled towards another surface oxygen, shown in Figure 1(c). A Pu−Ow bond distance of 2.14 Å is yielded in water dissociative adsorption on the PuO2 (110) surface. The adsorption energy calculated based on Eq.(1) are -0.89 and -1.17 eV for the molecular and dissociative geometries, respectively. It indicates that, on the PuO2 (110) surface, water dissociative adsorption is more energetically favorable than the molecular one. Such preference found in our work is consistent with previous studies. 26–30 In ref., 29 longer plutonium-oxygen distances, 2.72 and 2.44 Å, and lower adsorption 6

ACS Paragon Plus Environment

Page 7 of 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 1: Energy landscape of water molecule dissociation on the PuO2 (110) surface: (a) optimized molecular adsorption structure; (b) transition structure of the dissociation extracted from the MD simulation trajectory; and (c) optimized dissociative adsorption structure. energy, -0.88 and -1.14 eV, were reported for molecular and dissociative adsorption states, respectively. The small geometry and energy difference between our results and those reported by Tegner et al. can be attributed to the inclusion of dispersion corrections in our work. 30 To gain insight into the dynamics of water dissociation on the PuO2 (110) surface, we performed a series of AIMD simulations of 1−3 water molecules on the PuO2 (110) surface at ambient conditions. In the initial configuration of MD simulations, water molecules are randomly settled 4 Å above the surface. For the case of water monomer, as the molecule diffuses towards the surface, the water dissociation occurs by ejecting a hydrogen atom and forming surface hydroxyl group with an oxygen atom in the topmost surface layer. Following that, the remaining hydroxyl radical (Ow H) forms bond with surface plutonium within tens of femto seconds. A two-step dissociation of water molecule on the PuO2 (110) surface is observed. We extract a transition structure of water monomer dissociation on the PuO2 (110) surface from the MD simulation trajectory, where one of water hydrogen atoms breaks away and bonds to a surface oxygen, producing a hydroxyl group on the surface, as shown in Figure 1(b). We compute the normal modes of the dissociative water molecule in Figure 1(b), using density functional perturbation theory, and only one imaginary fre7

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

quency is found. At the transition state, an energy barrier of 0.18 eV is obtained with respect to the molecular adsorption state. The second step of the dissociation, where two surface hydroxyl groups are formed, is an exothermic process with an energy release of 0.46 eV. The projected density of states (PDOS) are then calculated to study the electronic interactions between water molecule and the PuO2 (110) surface in the dissociation path. The PDOS of two hydrogen and one oxygen atoms within water molecule are summed together to represent the molecular orbitals. In order to compare the molecular orbital energy levels of water in the molecular adsorption state and the isolated state, we place another water molecule in the molecular adsorption configuration of water on the PuO2 (110) surface (see Figure 1(a)), 15 Å away from the PuO2 surface, and increase the vacuum spacing to 30 Å along surface normal direction. As shown in Figure 2(a), three highest occupied molecular orbitals of the isolated water molecule, 1b1 , 3a1 , and 1b2 , declare themselves as three localized peaks at −10.09, −6.27 and −4.22 eV, respectively, where the Fermi level is shifted to zero. In the molecular adsorption state, the highest 1b2 orbital of the molecule is greatly broadened, as presented in Figure 2(b). The surface O-2p, Pu-5 f as well as Pu-6d states take part in the hybridizations with the 1b2 molecular orbital. Through Bader topological analysis, 40,41 we find 0.03 electron is donated from the surface to the water molecule during the electronic hybridization process. As the dissociation of water molecule proceeds, the electronic states of one surface oxygen atom increasingly interact with those of water molecule, as shown in Figure 2(c) and 2(d), which indicates the formation of new OH group between the surface oxygen atom and the dissociating hydrogen from water. In short, the dissociation of water monomer on the PuO2 (110) surface is the result of the hybridizations between three frontier orbitals of water molecule and the electronic states of the PuO2 surface. Such dissociation needs to conquer an energy barrier of 0.18 eV. In the case of water dimer, two water molecules are clustered by hydrogen bond (HB), as shown in Figure 3(a). When water dimer are adsorbed adjacent to the surface, the dissociation is initiated from the hydrogen atom in the HB acceptor of the dimer. The dissociating hydrogen forms hydroxyl group with surface oxygen atom and the process is followed by the bonding between the

8

ACS Paragon Plus Environment

Page 8 of 19

Page 9 of 19

4

(a)

2

1b1 3a1

1b2

0 -2 water molecule

-4 6

(b)

3

PDOS (arb. unit)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

0 -3 molecular state

-6 6

(c)

3 0 -3 transition state

-6 6

(d)

3 0 Osurf-p

-3

Pusurf-d

dissociative state

Pusurf-f

-6

water molecule

-12 -10 -8 -6 -4 -2

0

2

4

Energy (eV) Figure 2: The projected density of states (PDOS) for the (a) isolated water molecule, (b) molecular adsorption state, (c) transition state and (d) dissociative adsorption state of water molecule on the PuO2 (110) surface, respectively. Fermi level is shifted to zero. Black line represents the total PDOS of water molecule. Grey, blue and red lines represent the PDOS of surface O-2p, Pu-6d and Pu-5 f orbitals, respectively.

9

ACS Paragon Plus Environment

(a)

2.1 Å

(b)

1.9 Å

0 fs

70 fs

20 fs

0 fs

(d)

(c)

(e)

(f)

10 fs

20 fs

Figure 3: Snapshots of the dissociation process of water dimer (a-c) and trimer (d-f) on the PuO2 (110) surface. The dashed lines indicate hydrogen bonds (HBs) between water molecules. remaining hydroxide from water and a surface plutonium atom. It presents a similar two-step dissociation to that of water monomer on the PuO2 (110) surface. Furthermore, we compare the energy of the initial, transition and final configurations of the dimer dissociation obtained from the MD trajectory, as shown in the upper panel of Figure 3. The reaction is found to be exothermic by -0.42 eV with no energy barrier. It suggests a different type of dissociation mechanism of water dimer from that of water monomer, where an activation energy of 0.18 eV is required. Hydrogen bonding interaction between molecules facilitates the dimer dissociation on the PuO2 (110) surface, yielding a more energetically favorable process. Our simulation of three water molecules on the PuO2 (110) surface displays similar dissociation phenomena to that of water dimer. Hydrogen bonding interaction of water glues molecules in a cluster manner. Typical steps in the dissociation process is exhibited in the lower panel of Figure 3. To further demonstrate the effect of HB in water dissociation on the PuO2 (110) surface, we calculate the electron density difference of the initiation configure of water trimer dissociation on the surface, as displayed in Figure 4. We define the electron density difference as

10

Page 11 of 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 4: (a) Configuration of the initiation step of water trimer dissociation on the PuO2 (110) surface. (b) Electron density difference contour of the configuration shown in (a), as defined in the text. ∆ρ = ρsur f +3H2 O − ρsur f +H2 O − ρHHBA − ρHHBD , where ρsur f +3H2 O is the total electron density of wa2O 2O ter trimer adsorbed on the PuO2 (110) surface, ρsur f +H2 O is the electron density of the PuO2 (110) surface with the water molecule that is not directly hydrogen bonded to the dissociating molecule, ρHHBA and ρHHBD are the electron densities of the HB acceptor and donor molecules, respectively, 2O 2O highlighted in the blue square in Figure 4(a). Electron density difference contour in Figure 4(b) shows a clear electron redistribution around oxygen atom of the HB acceptor molecule. We also perform Bader atomic charge analysis and find no charge transfer between water molecules or between molecules and the PuO2 surface. It indicates that the electron redistribution occurs within the HB acceptor molecule and it may lead to a weakening of intramolecular OH bond and the initiation of the barrierless dissociation of water cluster. Interestingly, after the partial dissociation of the HB acceptor in water trimer, the HB donor molecule remains hydrogen bonded with the dangling hydroxyl group on the plutonium atom. As shown in Figure 5, a hydrogen transfer reaction takes place between HB donor molecule and surface dangling hydroxyl group. Because of the hydrogen bonding interaction, the hydroxyl group first separates from the surface plutonium to produce a free radical, and sequentially a hydrogen transfer arises from water molecule to the hydroxyl radical. The dehydrogenated water molecule turns to be a new hydroxide and bonds with adjacent surface plutonium atom, producing a new surface hydroxyl group. Such hydrogen transfer process consequently induces a hydroxyl group transfer on the PuO2 (110) surface. 11

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5: Hydrogen transfer process between water molecules on the PuO2 (110) surface.

4 Conclusions We have performed a series of AIMD simulations of water molecule and small clusters on the PuO2 (110) surface to investigate the dynamics and mechanism of water dissociation. Our results reveal that water dissociation on the PuO2 (110) surface is a two-step hydroxylation process both for single molecule and small clusters. The reaction is initiated by the dehydrogenation of water molecule to form hydroxyl group with surface oxygen atom, followed by successive surface hydroxylation on plutonium with the remaining hydroxide ion of the dissociating molecule. Despite the same dissociation steps, the mechanism of the reaction are different for water molecule and clusters. Water monomer dissociation, involving the hybridizations between the molecular orbitals of water and the electronic state of the top surface, requires an activation energy of 0.18 eV. In contrast, the dissociation of water dimer on the surface is found to be an exothermic reaction with no energy barrier. The process starts within the HB acceptor molecule in the cluster and results in a more stable dissociation state with respect to that of single molecule on the surface. Such facilitation in the dissociation of water clusters can be attributed to the hydrogen bonding interac-

12

ACS Paragon Plus Environment

Page 12 of 19

Page 13 of 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

tions between water molecules. Hydrogen bonding interaction can also induce a hydrogen transfer reaction between surface hydroxyl group and remaining water molecule, resulting in a hydroxyl group transfer along the surface.

5 Acknowledgement This work was supported by the Natural Science Foundation of China under Grant No. 21503019, the President Fund of China Academy of Engineering Physics under Grant No. 201402034 and the Foundations for Development of Science and Technology of China Academy of Engineering Physics under Grant No. 2015B0102020.

13

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

References (1) Haschke, J. M.; Allen, T. H.; Stakebake, J. L. Oxidation Kinetics of Plutonium in Air: Consequences for Environmental Dispersal. J. Alloys Compd. 1998, 271-273, 211–215. (2) Stakebake, J. L. Thermal Desorption Study of the Surface Interactions between Water and Plutonium Dioxide. J. Chem. Phys. 1973, 77, 581–586. (3) Stakebake, J. L.; Larson, D. T.; Haschke, J. M. Characterization of the Plutonium-water Reaction II: Formation of a Binary Oxide Containing Pu(VI). J. Alloys Compd. 1993, 202, 251–263. (4) Haschke, J. M.; Allen, T. H.; Stakebake, J. L. Reaction Kinetics of Plutonium with Oxygen, Water and Humid Air: Moisture Enhancement of the Corrosion Rate. J. Alloys Compd. 1996, 243, 23–35. (5) Haschke, J. M.; Ricketts, T. E. Adsorption of Water on Plutonium Dioxide. J. Alloys Compd. 1997, 252, 148–156. (6) Haschke, J. M.; Allen, T. H.; Morales, L. A. Reaction of Plutonium Dioxide with Water: Formation and Properties of PuO2+x . Science 2000, 247, 285–286. (7) Paffett, M. T.; Kelly, D.; Joyce, S.; Morris, J.; Veirs, K. A Critical Examination of the Thermodynamics of Water Adsorption on Actinide Oxide Surfaces. J. Nucl. Mater. 2003, 322, 45–56. (8) T. Gouder and A. Seibert and L. Havela and J. Rebizant, Search for Higher Oxides of Pu: A Photoemission Study. Surf. Sci. 2007, 601, L77. (9) McCleskey, T. M.; Bauer, E.; Jia, Q.; Burrell, A. K.; andSteven D. Conradson, B. L. S.; Mueller, A.; Roy, L.; Wen, X.; Scuseria, G. E.; Martin, R. L. Optical Band Gap of NpO2 and PuO2 from Optical Absorbance of Epitaxial Films. J. Appl. Phys. 2013, 113, 013515.

14

ACS Paragon Plus Environment

Page 14 of 19

Page 15 of 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(10) Petit, L.; Svane, A.; Szotek, Z.; Temmerman, W. M. First-Principles Calculations of PuO2 . Science 2003, 301, 498–501. (11) Sun, B.; Zhang, P.; Zhao, X.-G. First-principles Local Density Approximation+U and Generalized Gradient Approximation+U Study of Plutonium Oxides. J. Chem. Phys. 2008, 084705, 128. (12) Zhang, P.; Wang, B.-T.; Zhao, X.-G. Ground-state Properties and High-pressure Behavior of Plutonium Dioxide: Density Functional Theory Calculations. Phys. Rev. B 2010, 144110, 82. (13) Wen, X.-D.; Martin, R. L.; Henderson, T. M.; Scuseria, G. E. Density Functional Theory Studies of the Electronic Structure of Solid State Actinide Oxides. Chem. Rev. 2013, 113, 1063–1096. (14) Jomard, G.; Amadon, B.; Bottin, F.; Torrent, M. Structural, Thermodynamic, and Electronic Properties of Plutonium Oxides from First Principles. Phys. Rev. B 2008, 78, 075125. (15) Andersson, D. A.; Lezama, J.; Uberuaga, B. P.; Deo, C.; Conradson, S. D. Cooperativity Among Defect Sites in AO2+x and A4 O9 (A=U, Np, Pu): Density Functional Calculations. Phys. Rev. B 2009, 79, 024110. (16) Yang, Y.; Zhang, P. Chemical Bonds and Vibrational Properties of Ordered (U, Np, Pu) Mixed Oxides. J. Appl. Phys. 2013, 013501, 113. (17) Yang, Y.; Wang, B.-T.; Zhang, P. Electronic and Mechanical Properties of Ordered (Pu, U)O2 Compounds: A Density Functional Theory + U Study. J. Nucl. Mater. 2013, 433, 345–350. (18) Lu, Y.; Yang, Y.; Zhang, P. Charge States of Point Defects in Plutonium Oxide: A Firstprinciples Study. J. Alloys Compd. 2015, 649, 544–552. (19) Ao, B.; Qiu, R.; Lu, H.; Ye, X.; Shi, P.; Chen, P.; Wang, X. New Insights into the Formation of Hyperstoichiometric Plutonium Oxides. J. Phys. Chem. C 2015, 101-108, 119.

15

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(20) Ao, B.; Lu, H.; Qiu, R.; Ye, X.; Shi, P.; Chen, P.; Wang, X. First-Principles Energetics of Some Nonmetallic Impurity Atoms in Plutonium Dioxide. J. Phys. Chem. C 2015, 1487914889, 119. (21) Ao, B.; Qiu, R.; Lu, H.; Chen, P. Differences in the Existence States of Hydrogen in UO2 and PuO2 from DFT + U Calculations. J. Phys. Chem. C 2016, 18445-18451, 120. (22) Hernandez, S. C.; Holby, E. F. DFT+U Study of Chemical Impurities in PuO2 . J. Phys. Chem. C 2016, 13095-13102, 120. (23) Wu, X.; Ray, A. K. A Hybrid-density Functional Cluster Study of the Bulk and Surface Electronic Structures of PuO2 . Physica B 2001, 301, 359. (24) Jomard, G.; Bottin, F. Thermodynamic Stability of PuO2 Surface: Influence of Electronic Correlations. Phys. Rev. B 2011, 84, 195469. (25) Sun, B.; Liu, H.; Song, H.; Zhang, G.; Zheng, H.; Zhao, X.; Zhang, P. First-principles Study of Surface Properties of PuO2 : Effects of Thickness and O-vacancy on Surface Stability and Chemical Activity. J. Nucl. Mater. 2012, 426, 139–147. (26) Zs. Rak, R. E.; Becker, U. Hydroxylation-induced Surface Stability of AnO2 (An=U, Np, Pu) from First-principles. Surf. Sci. 2013, 608, 180. (27) Wu, X.; Ray, A. K. Density-functional Study of Water Adsorption on the PuO2 (110) Surface. Phys. Rev. B 2002, 65, 085403. (28) Jomard, G.; Bottin, F.; Geneste, G. Water Adsorption and Dissociation on the PuO2 (110) Surface. J. Nucl. Mater. 2014, 451, 28–34. (29) Tegner, B. E.; Molinari, M.; Kerridg, A.; Parke, S. C.; Kaltsoyannis, N. Water Adsorption on AnO2 (111), (110), and (100) Surfaces (An = U and Pu): A Density Functional Theory + U Study. J. Phys. Chem. C 2017, 121, 1675–1682.

16

ACS Paragon Plus Environment

Page 16 of 19

Page 17 of 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(30) Wellington, J. P.; Kerridge, A.; Austin, J.; Kaltsoyannis, N. Electronic Structure of Bulk AnO2 (An = U, Np, Pu) and Water Adsorption on the (111) and (110) Surfaces of UO2 and PuO2 from Hybrid Density Functional Theory within the Periodic Electrostatic Embedded Cluster Method. J. Nucl. Mater. 2016, 482, 124–134. (31) Kresse, G.; Furthmuller, ¨ J. Efficient Iterative Schemes for Ab Initio Total-energy Calculations using a Plane-wave Basis Set. Phys. Rev. B 1996, 54, 11169–11186. (32) Blochl, ¨ P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953–17979. (33) Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B 1999, 59, 1758–1775. (34) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. (35) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple [Phys. Rev. Lett. 77, 3865 (1996)]. Phys. Rev. Lett. 1997, 78, 1396. (36) Dion, M.; Rydberg, H.; Schröder, E.; Langreth, D. C.; Lundqvist, B. I. Van der Waals Density Functional for General Geometries. Phys. Rev. Lett. 2004, 92, 246401. (37) Román-Pérez, G.; Soler, J. M. Efficient Implementation of a van der Waals Density Functional: Application to Double-Wall Carbon Nanotubes. Phys. Rev. Lett. 2009, 103, 096102. (38) Klime˘s, J.; Bowler, D. R.; Michaelides, A. Chemical Accuracy for the van der Waals Density Functional. J. Phys.: Condens. Matter 2010, 22, 022201. (39) Klime˘s, J.; Bowler, D. R.; Michaelides, A. Van der Waals Density Functionals Applied to Solids. Phys. Rev. B 2011, 83, 195131. (40) Bader, R. F. W. Atoms in Molecules: A Quantum Theory; Oxford University Press: New York, 1990. 17

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(41) Tang, W.; Sanville, E.; Henkelman, G. A Grid-based Bader Analysis Algorithm Without Lattice Bias. J. Phys.: Condens. Matter 2009, 21, 084204.

18

ACS Paragon Plus Environment

Page 18 of 19

Page 19 of 19 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

19

ACS Paragon Plus Environment