Pair your accounts.

Export articles to Mendeley

Get article recommendations from ACS based on references in your Mendeley library.

Pair your accounts.

Export articles to Mendeley

Get article recommendations from ACS based on references in your Mendeley library.

You’ve supercharged your research process with ACS and Mendeley!

STEP 1:
Click to create an ACS ID

Please note: If you switch to a different device, you may be asked to login again with only your ACS ID.

Please note: If you switch to a different device, you may be asked to login again with only your ACS ID.

Please note: If you switch to a different device, you may be asked to login again with only your ACS ID.

MENDELEY PAIRING EXPIRED
Your Mendeley pairing has expired. Please reconnect
ACS Publications. Most Trusted. Most Cited. Most Read
DFT + U Study of the Adsorption and Dissociation of Water on Clean, Defective, and Oxygen-Covered U3Si2{001}, {110}, and {111} Surfaces
My Activity

Figure 1Loading Img
  • Open Access
Article

DFT + U Study of the Adsorption and Dissociation of Water on Clean, Defective, and Oxygen-Covered U3Si2{001}, {110}, and {111} Surfaces
Click to copy article linkArticle link copied!

  • Ericmoore Jossou*
    Ericmoore Jossou
    Department of Mechanical Engineering, College of Engineering, University of Saskatchewan, 57 Campus Drive, Saskatoon S7N 5A9, Saskatchewan, Canada
    *E-mail: [email protected] (E.J.).
  • Linu Malakkal
    Linu Malakkal
    Department of Mechanical Engineering, College of Engineering, University of Saskatchewan, 57 Campus Drive, Saskatoon S7N 5A9, Saskatchewan, Canada
  • Nelson Y. Dzade
    Nelson Y. Dzade
    School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, U.K.
  • Antoine Claisse
    Antoine Claisse
    Westinghouse Electric Sweden AB, SE-721 63 Västerås, Sweden
  • Barbara Szpunar
    Barbara Szpunar
    Department of Physics and Engineering Physics, College of Art and Science, University of Saskatchewan, 116 Science Place, Saskatoon S7N 5E2, Saskatchewan, Canada
  • Jerzy Szpunar*
    Jerzy Szpunar
    Department of Mechanical Engineering, College of Engineering, University of Saskatchewan, 57 Campus Drive, Saskatoon S7N 5A9, Saskatchewan, Canada
    *E-mail: [email protected] (J.S.).
Open PDFSupporting Information (1)

The Journal of Physical Chemistry C

Cite this: J. Phys. Chem. C 2019, 123, 32, 19453–19467
Click to copy citationCitation copied!
https://doi.org/10.1021/acs.jpcc.9b03076
Published July 11, 2019

Copyright © 2019 American Chemical Society. This publication is licensed under CC-BY.

Abstract

Click to copy section linkSection link copied!

The interfacial interaction of U3Si2 with water leads to corrosion of nuclear fuels, which affects various processes in the nuclear fuel cycle. However, the mechanism and molecular-level insights into the early oxidation process of U3Si2 surfaces in the presence of water and oxygen are not fully understood. In this work, we present Hubbard-corrected density functional theory (DFT + U) calculations of the adsorption behavior of water on the low Miller indices of the pristine and defective surfaces as well as water dissociation and accompanied H2 formation mechanisms. The adsorption strength decreases in the order U3Si2{001} > U3Si2{110} > U3Si2{111} for both molecular and dissociative H2O adsorption. Consistent with the superior reactivity, dissociative water adsorption is most stable. We also explored the adsorption of H2O on the oxygen-covered U3Si2 surface and showed that the preadsorbed oxygen could activate the OH bond and speed up the dissociation of H2O. Generally, we found that during adsorption on the oxygen-covered, defective surface, multiple water molecules are thermodynamically more stable on the surface than the water monomer on the pristine surface. Mixed molecular and dissociative water adsorption modes are also noted to be stable on the {111} surface, whereas fully dissociative water adsorption is most stable on the {110} and {001} surfaces.

Copyright © 2019 American Chemical Society

1. Introduction

Click to copy section linkSection link copied!

There is an increasing interest in the development and use of metallic nuclear fuels such as U3Si2 given the enhanced thermophysical properties compared to traditional uranium dioxide fuel (UO2). (1) However, a key issue to deal with is the ease of oxidation of metallic nuclear fuels in the presence of water, oxygen, or combination of both, which requires the synthesis of U3Si2 in controlled environments. A significant number of earlier experimental works have assessed the behavior of U3Si2 in corrosion-susceptible environments in comparison to other metallic and UO2 fuels. Recently, Nelson et al. investigated the behavior of U3Si2 following exposure to pressurized H2O that is typical in light water reactors (LWRs) at temperatures ranging from 300 to 350 °C. Their results showed that both UN and U3Si5 rapidly pulverize in less than 50 h at 300 °C, while the behavior of U3Si2 was superior but still below the corrosion resistance of UO2 fuel. It is worth mentioning that the mechanism of pulverization of U3Si2 might be due to spallation of UO2 or due to internal hydriding. (2,3) The formation of U3Si2H1.8 with volumetric increase has been shown by experiment and density functional theory (DFT) calculations. (4) In this work, we unraveled the detailed mechanism of the early oxidation of U3Si2 from an atomistic point of view, which is difficult to realize experimentally. Earlier theoretical works have focused on the bulk properties of U3Si2 within the framework of DFT + U formalism and molecular dynamics simulations. (4−10) For instance, Middleburgh et al. (4) investigated the defect evolution in U3Si2 and also proposed a phase diagram capable of predicting fuel behavior during burnup. Furthermore, Wang et al. (7) studied the electronic, structural, and elastic properties showing detailed bonding characteristics of U3Si2 by electron density of states (DOS), charge density, and charge density difference analysis which corroborate the work of Remschnig et al. (11) as regards the metallic nature of U3Si2.
Surface science provides fundamental insight into the chemistry and physics of corrosion in materials, but such experiments are expensive and require dedicated facilities, especially in the study of an actinide containing compound. Hence, theoretical surface science provides an alternative tool for investigating oxidation mechanisms in nuclear fuel materials. Bo et al. (12) modeled the surface properties of the low-index NpO2(111), (110), and (100) surfaces as well as the adsorption and dissociation behaviors of water on these surfaces using DFT + U calculations in combination with ab initio atomistic thermodynamic simulations. Their results showed that water dissociation is enhanced by oxygen vacancy, while coverage of the water molecules plays no significant role during molecular adsorption. (13) More recently, Wellington and co-worker investigated the adsorption and dissociation of water on pristine and reduced UO2 and PuO2 surfaces using the periodic electrostatic-embedded cluster method and Hubbard-corrected periodic conditions implemented in the DFT-based code. Oxygen vacancies were shown to be easily formed on PuO2 compared to UO2, which is due to ease of Pu reduction in comparison to U metal ions. Their results also showed that dissociation is favored over molecular adsorption of water both on defect-free and oxygen deficit surfaces of UO2 and PuO2.
Given the success of first-principles DFT + U calculations in the modeling surfaces, we recently studied the adsorption of molecular and dissociated O2 on stoichiometric U3Si2 lower index surfaces. (14) To further our understanding of the underlying mechanism of oxidation, we have considered the synergetic effect of water and oxygen interaction with U3Si2 perfect surfaces. However, as real surfaces are never perfect but contain defects due to the method of synthesis or service conditions, in the present work, we have explored the effects of U and Si vacancies on the adsorption mechanisms of water and oxygen on the {001}, {110}, and {111} surfaces of U3Si2. Surface vacancies in UN, UO2, PuO2, and CeO2 have been shown in previous studies to affect the oxidation behavior of such surfaces in the presence of oxidizing and hydriding agents. (12,13,15−19) The fundamental aspects of oxygen and water adsorption, including the initial adsorption geometries, adsorption energies, structural parameters, and electronic properties, are presented. Our results reveal a chemical picture of the initial steps involved in the oxidation process of the U3Si2 surfaces in the presence of oxygen and water (considering both molecular and dissociative adsorption).

2. Computational Methodology

Click to copy section linkSection link copied!

The first-principles calculations were performed within the plane-wave pseudo-potential DFT technique, (20,21) as implemented in the Quantum ESPRESSO code. (22) The exchange–correlation functional potential was described by the generalized gradient approximation (GGA) in the Wu–Cohen (GGA-WC) formulation. (23) All our calculations deployed the norm-conserving Wu–Cohen (WC) pseudopotentials (where these potentials treat 7p06d15f37s1 as valence electrons for U and 3s23p6 as valence electrons for Si). Due to the onsite Coulomb repulsion among the localized U 5f electrons, we used the Hubbard (DFT + U) correction to account for the strong correlation effect. (24) We have used an effective U value of 1.5 eV, which has been shown to give an accurate description of the structural parameters and the electronic properties of U3Si2. (14)
The Fermi surface effects were treated by the smearing technique of Methfessel–Paxton, (25) using a smearing parameter of 0.02 Ry (0.27 eV). An energy threshold defining self-consistency of the electron density was set to 10–8 eV and a β mixing factor of 0.3. The Brillouin zone integration was performed using 7 × 7 × 10 and 5 × 5 × 1 Monkhorst–Pack (22)k-point grids (centered at the Γ point) for the bulk U3Si2 and the surface models, respectively. Structural relaxation was carried out to minimize the energy using the conjugate gradient method within the Broyden–Fletcher–Goldfarb–Shanno algorithm, (26) until the magnitude of the residual Hell–Feynman force on each relaxed atom reached 0.01 eV Å–1. Visualization and analysis of the structures were performed using the VESTA program. (27,28)
The {001}, {110}, and {111} surfaces were considered for the oxygen and water adsorption calculations as they are the dominant growth facets of U3Si2. The different surface structures were created from the fully optimized bulk tetragonal U3Si2 structure (Figure 1) to eliminate the presence of fictitious forces during surface relaxation using the METADISE code, (29) which ensures the creation of surfaces with a zero dipole moment perpendicular to the surface plane. A vacuum of 15 Å was introduced to the surface models in the z direction, which is large enough to avoid any spurious interactions between periodic slabs. Bader population analysis was carried out on all adsorbate–substrate systems to quantify any charge transfer between the U3Si2 surface and O2/H2O species. (30)

Figure 1

Figure 1. Unit cell of ordered uranium silicide containing two formula units of U3Si2 (10 atoms) (color scheme: U = gray and Si = blue).

3. Results and Discussion

Click to copy section linkSection link copied!

3.1. Defective Surface Models

3.1.1. Surface Defect Energies and Stability

Defective surfaces (vacancies) were created by removing one Si/U atom at a time from either the first or second layer (Figure 2) to access their energetics or stability.

Figure 2

Figure 2. Optimized surface geometry of (a) U3Si2{001}, (b) U3Si2{110}, and (c) U3Si2{111} with Si and U1 surface vacancy represented by a rectangular red box in the second and third rows respectively (color scheme: U = gray and Si = blue).

The silicon vacancies are denoted as VSi, whereas the two considered uranium vacancy sites are denoted as VU1 and VU2, respectively, on the surface. The method of Wellington et al. (19) was employed in the creation of vacancy such that the uranium and silicon atoms removed from the supercell, both from the surface and subsurface layers, are fully coordinated by other quantum mechanically treated atoms and lie close to the center of the supercell to avoid problems that may occur in forming a defect at the edge. (19) The vacancy formation energy was calculated using eq 1
(1)
where Esurface+Xvac is the energy of the relaxed slab with X = Si or U vacancy, Esurface is the energy of the relaxed stoichiometric slab of the same type, and EX is the reference energy for X = Si and U as the energy per atom in its fundamental face-centered cubic and orthorhombic structure, respectively. Presented in Table 1 are the calculated vacancy formation energies on the {001}, {110}, and {111} U3Si2 surfaces. We found higher Si vacancy formation for the {110} surface (2.62 eV) than for the {001} surface (0.16 eV) and {111} surface (1.54 eV), and a similar trend is observed in the subsurface Si vacancies formation. Clearly, Si vacancy on the {110} subsurface (2.62 eV) is close to 2.48 eV, calculated for bulk U3Si2 by Andersson et al. (31) This is not surprising, given that the {110} surface has the lowest surface energy and similar neighboring atom coordination. There are two U sites denoted here as U1 and U2 with vacancy defect formation energies of 1.64 and 2.65 eV, respectively, for the bulk. The first and second layer U vacancy energies were calculated as presented in Table 1. Generally, it is easier to form a vacancy in the surface layer compared to the subsurface, which is due to the difference in the nearest neighbor atoms and the coordination numbers.
Table 1. Uranium and Silicon Vacancy Formation Energies (eV) of the {001}, {110}, and {111} Surfaces of U3Si2
 first layersecond layer
surfaceVSiVU1VU2VSiVU1VU2
{001}0.161.721.871.742.132.88
{110}2.623.182.411.622.032.61
{111}1.543.813.702.441.463.10
The formation of Si and U vacancies resulted in a nonstoichiometric surface, which allows variation in the surface energies as the chemical potential, μ, of Si and U changes. The thermodynamic stability of a given surface, in general, depends on the specific chemical environment. To determine the stability of the surfaces due to vacancies, we calculate the surface energy, γ, as a function of the Si and U chemical potential, respectively. At zero temperature the surface energy of a crystal may be derived from a N-layer slab using eq 2
(2)
In this equation, Esurf is the slab energy with all atomic coordinates relaxed unconstrainedly; NEbulk is the energy of an equal number, N of bulk U3Si2 atoms; A is the area of the slab surface; and the factor of 2 reflects the fact that there are two surfaces for each slab. Due to the surface vacancy defects, the surface energy depends on the specific thermodynamic conditions, i.e., the reservoir with which the atoms of the compound are exchanged in a structural transition. Therefore, the chemical potential of the constituents enters the surface energy. The most stable surface structure is determined by the minimum of the free energy which at zero temperature is given by eq 3
(3)
In the case of U3Si2, eq 3 can be expressed by eq 4 to account for the chemical potential μ of U and Si atoms
(4)
The μU and μSi are bounded by a set of conditions. Assuming U and Si are in thermal equilibrium with the U3Si2 crystal, it would imply that
(5)
Furthermore, since there is no precipitation of U and Si on the U3Si2 surface, the following conditions must also be fulfilled
(6a,b)
Also, the heat of formation of bulk U3Si2, ΔHf (eV f.u.–1), is defined as
(7)
Combining eqs 5, 6a,b, and 7, we obtain a range for possible values of the Si chemical potential
(8)
Combining eqs 4 and 5, gives the surface energy as a function of μSi
(9)
Similar to eq 9, the surface energy can also be expressed as a function μU
(10)
At ambient temperature and pressure, Esurf can be taken as the total energy from DFT calculation neglecting contributions from configurational or vibrational entropies. In Figure 2, we show the results of the calculated surface energies of the {001}, {110}, and {111} surfaces of U3Si2 as functions of the change in silicon chemical potential μSi, ranging from Si4 (−7.46 eV) to that of an isolated Si atom (0 eV), which denotes lower and upper limit of the Si environment, data from NIST database. (32) The lower limit, which is the Si-poor environment, is defined by the decomposition of the silicide into U and Si. Meanwhile, the upper limit is the Si-rich state that corresponds to a situation where the gas phase is so rich in silicon that they condensed on the substrate. In the case of stoichiometric surfaces considered in this work, the surface free energy is independent of μSi as expected, since the coefficient of eq 6a,ba,b containing μSi vanishes. For the nonstoichiometric cases, the surface free energy is a linear function of μSi according to eq 7. It can be seen from Figure 3 that the equilibrium surface stoichiometry of the U3Si2 surfaces can be tuned by changing μSi through the use of different U3Si2 growth (Si vacancy is considered here) and annealing environments. It is clear from the increase in surface energies that the stoichiometric surfaces are found to be more stable than the silicon deficient (silicon vacancy) surfaces. This implies that surface adsorption of an oxidizing adsorbate is more stable on the defective surface in comparison to the stoichiometric surface.

Figure 3

Figure 3. Calculated surface energies of the {001}, {110}, and {111} surfaces of U3Si2 as functions of the change in silicon chemical potential μSi.

3.2. Adsorption of Water Molecule

The adsorption of the water molecule on the U3Si2 surfaces is an important starting step toward understanding its early oxidation mechanisms. The first interest of this study is therefore to determine the lowest-energy adsorption structures and modes of water on the low-index U3Si2{001}, {110}, and {111} surfaces, and to characterize the strength of their interaction and the extent of O–H bond activation during the dissociation process. Different H2O initial adsorption possibilities, including Hw and Ow head-on configurations, have been subjected to geometry optimization until the residual forces on all atoms were ≤0.03 eV Å–1. Prior to adsorption, we have calculated the reference energies, bond length, bond angle, and vibrational frequencies of one free H2O molecule and compared them with earlier theoretical results and available experimental data. The values of the O–H bond and the H–O–H angle of water are determined to be 0.970 Å and 104.5°, respectively, which are in good agreement with the previous experimental and theoretical values. (33,34) The calculated scissoring bend, asymmetric and symmetric stretching vibrational frequencies are calculated to be 1555, 3523, and 3635 cm–1, which agree well with the experimental results (1596, 3652, and 3756 cm–1), (35) thus making our calculation method reliable.
The adsorption energy (Eads), which characterizes the strength of water adsorption has been defined as the mean adsorption energy per molecule of H2O–U3Si2 interaction using eq 11
(11)
where Esurface+n(H2O) is the total energy of the substrate–adsorbate system in the equilibrium state, Esurface and nEH2O are the total energies for the surface and free molecules, respectively, and n is the number of water molecules in the cell. By this definition, a negative value of Eads indicates an exothermic and stable adsorption, whereas a positive value indicates an endothermic and unstable adsorption.

3.2.1. Water Adsorption and Dissociation on Clean U3Si2{001}

For the adsorption of molecular water on the U3Si2{001} surface, different adsorption sites and configurations were explored to determine the lowest-energy adsorption structures, as presented in Figure 4. The calculated adsorption energies and the optimized interatomic bond distances are summarized in Table 2. The calculated lowest-energy water adsorption structure is presented in Figure 4a, wherein the water molecule adsorbs through the O atom at the U site (U–O = 2.56 Å), releasing an adsorption energy of −3.70 eV. The adsorption of water at the Si site (O–Si = 2.376 Å) released an adsorption energy of −3.11 eV. When the water molecule is adsorbed with the hydrogen atoms pointing toward the surface Si site (Figure 4b), it moved away perpendicularly from the surface during energy minimization until the closest H–Si distance was 3.044 Å. The adsorption energy of this configuration was calculated to be −1.97 eV, while the α(H–O–H) bond angle and O–H bonds are obtained at 102.32° and 0.965/0.967 Å, respectively. Compared to the lowest-energy water adsorption on U3Si2{001}, Bo et al. calculated the adsorption energies of −2.07 and −1.27 eV for the lowest-energy adsorption structures of the water monomer on UN(001) and UO2(110), respectively, which suggest that U3Si2{001} is more reactive toward water adsorption than both UN(001) and UO2(110). (36,37) In all three adsorption modes, the O–H bond lengths were slightly elongated and the α(H–O–H) bonds bond angle was larger compared to the gas-phase free H2O molecule in a vacuum, suggesting that the O–H bonds are activated to some extent when water is adsorbed on U3Si2{001}. In the lowest-energy O–U structure (Figure 4a), the two O–H bonds are calculated at 0.975 and 0.996 Å compared to the gas-phase molecule at 0.970 Å and this is consistent with O–H bond stretching vibrational frequencies presented in Table 1, whereas the α(H–O–H) bonds bond angle is obtained at 107.6° compared to the gas-phase value of 104.5°, which is supported by the red and blue shifts in the bending vibrational frequencies as shown in Table 2. The stretched O–H bond lengths are indicative of weaker O–H bonds, resulting from the π-antibonding occupation. The activated O–H bonds suggest that these molecular adsorption states are likely precursors for H2O dissociation.

Figure 4

Figure 4. Top and side views of the relaxed adsorption structures of molecular H2O adsorbed at: (a) O–U, (b) H–Si, and (c) O–Si (color scheme: U = gray, Si = blue, H = white, and O = red).

Table 2. Calculated Adsorption Energy (Eads), Relevant Bond Distance (d), Vibrational Frequencies, and Variation of the Total Bader Charge of Molecular (H2O) and Dissociated (OH + H) Adsorbed on the {001}, {110}, and {111} Surfaces of U3Si2
surfaceadsorbateconfig.Eads. (eV)d(O–H1) (Å)d(O–H2) (Å)d(Si–O) (Å)d(U–O) (Å)d(U–H) (Å)d(Si–H) (Å)α(HOH) (deg)νb (cm–1)νs (cm–1)νas (cm–1)Δq (e)
  H2O 0.9700.970    104.490155535233635 
{001}H2OO–U–3.700.9750.996 2.622  107.6421534339036610.08
O–Si–3.110.9940.9762.3762.562  105.3991543355336750.09
H–Si–1.970.9750.977    102.3231540376738620.05
OH + HOH–Si: H on adj. U–6.810.983  2.598 1.551    2.13
OH–Si: H on Si–5.430.971 1.7412.618 1.501    1.21
OH–Si: H on U–3.610.965 1.644 2.119     1.17
{110}H2OO–Si–2.460.9960.9782.120   107.3431532341836730.15
O–U–0.600.9851.005    104.0411600332635900.16
H–Si0.550.9710.970    103.2321601342535710.05
OH + HOH–Si: H on U–3.240.979 1.676 2.3961.583    0.44
OH–U: H on Si–2.270.974  2.200 1.510    0.17
{111}H2OO–U–1.400.9770.995 2.548  105.4541543343236900.19
H–Si–0.200.9570.977    99.1811575367739480.05
O–Si1.020.9770.982    102.2051622358436910.01
OH + HOH–U: H on Si–3.300.971  2.194 1.697    2.18
OH–Si: H on U–1.210.969 1.702 2.181     0.34
Dissociative adsorption of H2O on the clean U3Si2{001} surface is found to be highly exothermic with adsorption energies of −6.81, −5.43, and −3.61 eV, respectively, for the different configurations presented in Figure 5 and Table 2. Consistent with their stronger adsorption, the collective amount of charge transfer gained by the dissociated H2O species is larger than those of the molecular adsorbed water systems. Due to its metallic nature, the U ions donate more electrons than Si ions during bond formation. The preference for dissociative over molecular adsorption of water can be attributed to the fact that the energy required to break an O–H bond requires less energy compared to the energy released in the formation of the Si–H and U–OH bonds on the U3Si2{001} surface. In like manner, on several metals and oxide surfaces, the dissociative state of H2O is thermodynamically more stable than the molecularly adsorbed state, for example, on Cu and Al2O3, because the breaking of an O–H bond is effectively balanced by the formation of a metal–O and another O–H bond with a surface oxygen. (38,39)

Figure 5

Figure 5. Top and side views of the relaxed adsorption structures of dissociated water configuration with (a) OH–Si: H on adj. U, (b) OH–Si: H on Si, and (c) OH–Si: H on U, on the U3Si2{001}-Si terminated surface (color scheme: U = gray, Si = blue, H = white, and O = red).

3.2.2. Water Adsorption and Dissociation on Clean U3Si2{110}

As on the U3Si2{001} surface, we have also explored different adsorption configurations of H2O on the {110} surface, including O adsorbed laterally at U or Si sites or H placed head-on Si sites. The optimized adsorption structures are shown in Figure 6, whereas the energies of adsorption, geometric parameters, and Bader charge transfer analysis are presented in Table 2. The lowest-energy water adsorption configuration at U3Si2{110} was calculated to be the O–Si structure (Figure 6a), wherein the water molecule is adsorbed at the Si site through the oxygen atom (O–Si = 2.120 Å), releasing an energy of −2.46 eV. In this structure, we observed elongation of the O–H bond calculated at 0.996 and 0.978 Å, and broadening of the α(H–O–H) bond angle (107.3°). The next stable configuration is the O–U structure (Figure 6b), which released an adsorption energy of −0.60 eV. The interacting O–U distance is calculated at 2.805 Å, and O–H bonds and α(H–O–H) bond angle are calculated at 0.985/1.005 Å and 104.0°, respectively. When the H atoms placed head-on Si (Figure 6c), the adsorption process is found to be endothermic by 0.55 eV and an average O–H bond length is 0.971 Å and there is a preferential Si–O bond formation with an interatomic distance of 1.932 Å after energy minimization. From Bader populations analysis, we found that the adsorption process is characterized by charge transfer from the surface to the water molecule, with the water molecule gaining a charge of 0.15e, 0.16e, and 0.05e in the O–Si, O–U, and H–Si configurations, respectively. We also observed that the change in the vibrational modes and bond strengths is related to the Bader charge transfer between the surface and the adsorbed water molecule (Table 1).

Figure 6

Figure 6. Top and side views of the relaxed adsorption structures of molecular oxygen adsorbed at (a) O–Si, (b) O–U, and (c) H–Si on the U3Si2{110} surface (color scheme: U = gray, Si = blue, H = white, and O = red).

Compared to the molecular adsorption, dissociative water adsorption is found to have average adsorption energies of −3.24 and −2.27 eV for OOH placed head-on Si and head-on U, respectively. For dissociative adsorption through the HOH, we initially placed an HOH above the Si surface. After optimization, a Si–O bond was formed through the OOH atom, which means that dissociative adsorption through the HOH is unstable. The dissociative adsorption structures are characterized by more electrons being transferred from the surface to the dissociating species as evident by the calculated Bader charges of 0.44e and 0.17e. The charge transfer from the U3Si2 surface to the OH group is a necessary condition for water to bind to the surface, which is consistent with water dissociation on metal surfaces. (40) Furthermore, the above energetics results indicate that the dissociative adsorption of water is notably stronger than the molecular adsorption on the U3Si2{110} surface. We also observed the formation of U–H with bond lengths of 2.396 and 2.275 Å in the dissociative configurations, which is required for the H+ to achieve stability. In comparison to previous work, the U–H distances are 2.30–2.32 Å in α and β-UH3 configurations. (41) The schematic representations of adsorption structures of dissociated water on the clean U3Si2{110} surface are shown in Figure 7.

Figure 7

Figure 7. Top and side views of the relaxed adsorption structures of dissociated H2O at (a) OH–Si: H on U and (b) OH–U: H on Si on the U3Si2{110} surface (color scheme: U = gray, Si = blue, H = white, and O = red).

3.2.3. Water Adsorption and Dissociation on Clean U3Si2{111}

The geometries of the lowest-energy adsorption structures obtained on the U3Si2{111} surface are shown in Figure 8a–c. Compared to the {001} and {110} surfaces, the water molecules are adsorbed least strongly on the {111} surface. The lowest-energy adsorption structure on the U3Si2{111} surface is calculated to be the O–U configuration (Figure 8a), which released an adsorption energy of −1.40 eV. In this structure, the water molecule binds via the O atom lying laterally above the U atom at an O–U distance of 2.548 Å, with the α(H–O–H) bond angle that slightly increased to 105.5° compared to the gas-phase value of 104.5°. When the H2O molecule is adsorbed with the H atoms pointing toward the surface Si atoms (Figure 8b), an adsorption energy of −0.20 eV was released. The closest distance between H and Si/U sites is calculated to be 2.632/2.942 Å ruling out the formation of silane and uranium hydrides. The O–H bond lengths are elongated within the range of 0.977–0.995 Å while the α(HOH) bond angle increased slightly by 0.964° as presented in Table 2. The adsorption of H2O through Ow is unstable with an endothermic energy of 1.20 eV and tends to move further away from Si after optimization as shown in Figure 8c.

Figure 8

Figure 8. Top and side views of the relaxed adsorption structures of molecular H2O adsorbed at (a) the O–U site, (b) the H–Si site, and (c) the O–Si site on the U3Si2{111} surface (color scheme: U = gray, Si = blue, H = white, and O = red).

For dissociative adsorption of water on U3Si2{111}, the OH radical and one H atom are placed on the surface with two different adsorption modes as shown in Figure 9a,b. In Figure 9a, the OH forms a bond with a surface U atom while the H atom bonded to a neighboring Si atom. The formed U–O, Si–O, and Si–H bond lengths are 2.655, 1.722, and 1.610 Å, respectively. The calculated adsorption energy for this configuration is −0.76 eV, whereas the adsorption of OOH on U leads to the formation of U–O (2.194 Å) and Si–H (1.697 Å) bonds after optimization requiring an energy of −3.30 eV. The larger adsorption energies suggest a preference for dissociative water adsorption over molecular adsorption on the U3Si2{111} surface, which is consistent with previous studies of uranium-bearing systems. (37,42,43)

Figure 9

Figure 9. Top and side views of the relaxed adsorption structures dissociated water configuration with (a) O–U: H and (b) O–Si: H, on the U3Si2{111}-Si terminated surface (color scheme: U = gray, Si = blue, H = white, and O = red).

3.3. Effects of Surface Coverage, O-Covered, and Surface Vacancy on Adsorption of Water

3.3.1. Water Adsorption at Higher Coverage

The effect of coverage on the adsorption properties of water on the U3Si2 surfaces was investigated by adsorbing up to four water molecules in a molecular, mixed, and partially dissociative modes. The optimized structures of the most stable adsorption modes are shown in Figure 10, whereas the adsorption energies are reported in Table 3. We found no clear trends in the adsorption energies with an increasing number of water molecules on the U3Si2 surfaces. At the {001} surface, the adsorption energies of one, two, three, and four water molecules representing a coverage of 0.25, 0.50, 0.75, and 1.00 monolayers (ML) are −3.70, −2.83, −3.02, and −3.10 eV, respectively. At the {110} surface, the adsorption energy increased from −2.46 eV for one water molecule to −3.23, −3.06, and −2.58 eV for two, three, and four water molecules, respectively. A decrease in the adsorption energy is observed at the {111} surface with an increasing number of water molecules; −3.30 for one water molecule compared to −1.42, −2.01, and −2.53 eV for the two, three, and four water molecules, respectively.
Table 3. Adsorption Energies (eV) per Water Molecule, Mixed Molecular and Dissociated Water on U3Si2 Surfacesa
adsorbate{001}{110}{111}
1M–3.70–2.46–3.30
2M–2.83–3.23–1.42
3M–3.02–3.06–2.01
4M–3.10–2.58–2.53
1D + 3M–3.25–2.68–2.69
2D + 2M–3.26–2.69–2.65
3D + 1M–4.12–2.72–2.48
4D–3.54–3.68–3.45
a

Where M = H2O and D = OH + H.

Apart from the molecularly adsorbed water molecules, we have also explored mixed adsorption modes, wherein some of the water molecules are dissociated and some remained molecularly adsorbed, and investigated the case of fully dissociated modes in which all of the water molecules are dissociated at full coverage. The optimized structures of the most stable adsorption states are shown in Figure 10, and the adsorption energies are reported in Table 3. The mixed adsorption modes of molecular and dissociative water results in the formation different oxides, hydroxides, and oxygen molecules on the U3Si2 surface depending on the ratio of dissociative to molecularly adsorbed water as shown in Figure 10 for the most stable adsorption energies. The remaining structures are given in the Supporting Information (Figures S1–S3). For instance, on the {001} surface, in the case of 75/25 mix adsorption, the oxygen molecule was formed with a bond length of 0.99 Å directly on the surface Si with a bond distance of 1.70 Å typical of SiO2 formation. The strongest and most stable adsorption modes at the {001} surface are predicted for 3D + 1M, which released an adsorption energy of −4.12 eV, compared to the 4D structures on the {110} and {111} surfaces, which released adsorption energies of −3.68 and −3.45 eV, respectively. Even though U3Si2 is a metallic fuel, it is important to draw a comparison between the silicide and urania, which is the standard fuel for LWRs. Hence, it is interesting to note that studies by Bo et al. (37) on UO2 and recent work by Tegner et al. (44) predicted the mixed 50/50 molecular/dissociative adsorption modes as the most stable adsorption configurations. The differences in the results can be attributed to differences in the crystallographic arrangement of the two systems and their electronic structures, UO2 is a semiconductor with an experimental band gap of 2.1 eV, (45,46) whereas U3Si2 is metallic. Note, however, that the energy difference between the fully dissociative case and the fully molecular case is −0.44, −1.10, and −0.92 eV on the {001}, {110}, and {111} surfaces, respectively. The increase in adsorption energies in the dissociative configuration is due to the formation of stronger intramolecular hydrogen bonds on the crowded surfaces compared to the purely molecular mode.

Figure 10

Figure 10. Top and side view of the relaxed adsorption structures of mixed (dissociative to molecular) 3:1 adsorption of H2O on (a) {001} and fully dissociative adsorption of H2O on (b) the {110} and (c) {111} U3Si2 surface (color scheme: U = gray, Si = blue, H= white, and O = red).

We also carried out a detailed analysis of the bond lengths for molecular and dissociative adsorption of multiple water molecules on the surfaces as shown in Table 4. We observed a shorter hydrogen bond length for the Si–H compared to U–H on the U3Si2{001}, {110}, and {111} surfaces, which suggest stronger hydrogen bonding in Si–H when it is formed on the surface. There is also strong adsorption of molecular H2O on the surface forming U–H2O and Si–H2O complexes with the exception of the {111} surface where the formation of Si–H2O is not observed.
Table 4. Selected Interatomic Distances (Å) for Molecular and Dissociative Water on the U3Si2{001}, {110}, and {111} Surfaces at Coverages from 0.5 to 1.0 Monolayers (ML)
bond type{001} 0.5–1.0 ML (Å){110} 0.5–1.0 ML (Å){111} 0.5–1.0 ML (Å)
Si–O1.70N/AN/A
Si–H1.50–1.711.57–1.60N/A
Si–OH1.66–1.681.67–1.681.62–1.70
Si–H2O1.89–2.112.02–2.06N/A
U–O2.13N/AN/A
U–H2.17–2.342.34–2.452.16–2.35
U–OHN/A2.31–2.782.48–2.71
U–H2O2.51–2.672.35–2.652.64–2.69

3.3.2. Water Adsorption and Dissociation on O-Covered U3Si2{100}, {110} and {111} Surfaces

Preadsorbed oxygen atoms play an important role in the activation of the O–H bond and further increase the rate of surface oxidation due to rapid water dissociation. (47,48) Hence, it is important to investigate the adsorption behavior of H2O on oxygen-covered U3Si2{001}, {110}, and {111} surfaces. Prior to investigating the adsorption of H2O on oxygen-covered U3Si2{001}, {110}, and {111} surfaces, we have systematically characterized the adsorption of atomic oxygen on different U3Si2 surfaces to elucidate their surface oxide formation. In our previous study, (14) we have provided comprehensive information regarding the mechanism of oxide layer formation on the U3Si2 surface by adsorbed O2 molecule. The adsorption of atomic O is found to be energetically more favorable at U sites than Si sites on the {001}, {110}, and {111} U3Si2 surfaces (Table S4 and Figure S5). The adsorption energies at the U and Si sites were calculated to be −2.44 and −0.67 eV on the {001} surface, −2.71 and −0.77 eV on the {110} surface, and −2.85 and −0.81 eV on the {111} surface. Similar results were predicted for molecular O2 at the different U3Si2 surfaces, which implies that the formation of UO2 on pristine U3Si2 is favored over SiO2. (14)
For the adsorption of water on the O-covered U3Si2 surfaces, various coadsorption structures were explored with the water adsorbed at the neighbor or distant sites from the preadsorbed O atoms. The coadsorption energies between the H2O and O on the U3Si2 surface is calculated as follows
where EH2O, EO2, Esurface, and Esurface+(H2O+O) are the total energy for the free molecule of water, molecular oxygen, the clean U3Si2 surface, and the coadsorbed (H2O + O) + U3Si2 surface systems, respectively. In most cases at the {001} surface, we found that the water molecule coadsorbed with oxygen dissociate to form OH ions due to the attractive force between the hydrogen atoms and preadsorbed O atoms. The most favorable coadsorption mode on the {001} surface is predicted for the configuration in which the water molecule is adsorbed at the U site near the O atom that is preadsorbed at the Si site. Due to the strong attractive force between the two adsorbates, the H2O molecule spontaneously dissociate to form two hydroxyl species, releasing an adsorption energy of −5.47 eV (Table 5), which is consistent with water dissociation on oxygen-covered metal surfaces. (34,49) The U–OH and Si–OH interaction bond lengths are calculated to be 2.339 and 1.641 Å, respectively. Consistent with the strong adsorption, the two OH species draws a combined charge of 1.95e from the interacting surface species. When the water molecule remained molecularly adsorbed at a Si site near a preadsorbed O atom (Figure 11b), an adsorption energy of −3.54 eV is released.

Figure 11

Figure 11. Top and side views of the relaxed adsorption structures of molecular H2O (a) H–U: O on U, (b) Ow–Si: O on Si, (c) Ow–Si: O on U, and (d) Ow–U: O on Si U3Si2{001} surface (color scheme: U = gray, Si = blue, H = white, and O = red).

Table 5. Calculated Adsorption Energy (Eads) and Relevant Bond Distances for H2O Coadsorbed with Atomic Oxygen on U3Si2{001}, U3Si2{110}, and U3Si2{111} Surfaces
surfaceadsorbateconfig.Eads (eV)d(Ow–H1) (Å)d(Ow–H2) (Å)d(O–H2) (Å)d(Si–O) (Å)d(Si–OH) (Å)d(U–O) (Å)d(U–OH) (Å)α(HOH) (deg)Δq (e)
{001}H2O + OHw–U: O on U–5.470.977 1.026 1.641 2.339 1.95
Ow–Si: O on Si–3.540.9700.988 1.729 2.525 106.7601.83
Ow–Si: O on U–3.490.975 0.953 1.680 2.525 1.79
Ow–U: O on Si–2.870.973 1.0991.675    1.75
{110}H2O + OOw–U: O on U bridge–7.410.9750.975 2.763 2.268 105.5231.21
Ow–Si: O on U–4.130.982    2.236  0.92
Hw–Si: O on U–3.670.978 1.154   1.80  
{111}H2O + OOw–Si: O on U–1.390.979 1.006 1.710 2.145 1.67
Ow–U: O on Si–0.810.9770.974 1.695 2.567 109.9751.87
Hw–Si: O on U–0.270.974    2.021   
We have also investigated the adsorption and dissociation of water on the U3Si2{110} surface in the presence of preadsorbed atomic oxygen. Three coadsorption modes have been explored with the H2O and O coadsorbed at the neighboring surface U and U bridge sites. The optimized adsorption geometries are shown in Figure 12a–c, while the coadsorption energies and the relevant optimized geometric parameters are summarized in Table 5. When water is adsorbed with the Ow head-on Si and the preadsorbed O atom at the nearest neighbor U–U bridge, the coadsorption energy is −7.41 eV, which is more negative than the sum of the separate adsorption energies (−3.23 eV), suggesting a strong interaction evident by the stretching of the OH bonds and significant transfer of electrons (1.21e) from Si to the water molecule. It is clear that the oxygen atom has a stronger affinity for the uranium forming a U–O–U complex with a bond angle of 97.65°. Hence, the water molecule preferred to bond with the surface Si atom rather than deprotonate to form hydroxyl ions as was observed on iron sulfide surfaces. (34) When water is coadsorbed with the Ow head-on a neighboring U site, the coadsorption energy is calculated to be −4.13 eV, which is also more negative than the sum of the separate adsorption energies (−3.31 eV), suggesting that the oxygen atom promotes the dissociation of water due to deprotonation resulting in the formation of silane (Si–H), uranium oxide (U–O), and hydroxyl (O–H) species (Figure 12b). In this case, where Hw is placed on the Si atom, we observed a repulsion between H and Si while the nearest O atom forms a bond with one of the Hw leading to the formation of hydroxyl ions as shown in Figure 12c releasing an energy of −2.67 eV.

Figure 12

Figure 12. Top and side views of the relaxed adsorption structures of molecular H2O (a) Ow–Si: O on the U bridge, (b) Ow–Si: O on U, and (c) Hw–Si: O on U on the U3Si2{110} surface (color scheme: U = gray, Si = blue, H = white, and O = red).

The most favorable coadsorbed configuration of the (H2O + O)/U3Si2 system on the {111} surface is found to be the structure with the water molecule adsorbed on Si while one of the hydrogen atoms pointing toward the preadsorbed O atom on the top-U site (Figure 13a) releasing an energy of −1.39 eV, which is more negative than the sum of the separate adsorption energies of 0.39 eV. This is indicative of attraction between the two species on the U3Si2{111} surface for this configuration. This accounts for the dissociation of the water molecule leading to the formation of two Si–OH and U–OH. Although in the adsorption of the O atom on the Si site (Figure 13b) leads to the O atom, preferentially forming a trihedral network with U and Si atoms, this is facilitated by drawing electrons from a metal surface that is easy compared to a Si atom while there is stretching of the OH bonds (0.977 and 0.974 Å) with the Ow forming bond with the surface U atom.

Figure 13

Figure 13. Top and side views of the relaxed adsorption structures of molecular H2O (a) Ow–U: O and (b) Ow–Si: O on the U3Si2{111} surface (color scheme: U = gray, Si = blue, H = white, and O = red).

3.3.3. Water Adsorption on Defective (Nonstoichiometric) Surfaces

Considering that real surfaces are never perfect under reaction conditions and chemical processes often occur at defect sites, we have also investigated water adsorption at the defective surfaces containing one Si or U vacancies at the {001}, {110}, and {111} surfaces. The optimized structure of defective U3Si2 with H2O is shown in Figure 14 and the calculated adsorption energies and optimized geometry parameters are listed in Table 6. At the defective {001}, {110}, and {111} surfaces containing one Si vacancy site (Si-1), the adsorption energies of the water monomer in the molecular state is calculated to be −7.52, −3.18, and −7.0 eV, respectively. In {001} and {110} adsorption structures, there is a complete dissociation of the water molecule which is consistent with the chemisorbed nature of water (16) while the OH bonds on the {111} surface stretch to 1.245 and 0.980 Å, respectively, suggesting the instability of the water molecule on the nonstoichiometric surface. The computed bond distance of U–O (2.163–2.321 Å) and Si–O (1.536–1.892 Å) suggest the formation of UO2 and SiO2.

Figure 14

Figure 14. Top and side views of the relaxed adsorption structures of Si and U vacancy-assisted molecular H2O adsorption (a) H2O–U3Si2–x {001}, (b) H2O–U3Si2–x {110}, (c) H2O–U3Si2–x {111}, (d) H2O–U3–xSi2{001}, (e) H2O–U3–xSi2{110}, and (f) H2O–U3–xSi2{111} surface (color scheme: U = gray, Si = blue, H = white, and O = red).

Table 6. Adsorption Energies (eV) of the Water Molecule on Nonstoichiometric U3Si2{001}, {110}, and {111} Surfaces
adsorption site{001}{110}{111}
Si vacancy–3.43–3.18–2.92
U1 vacancy–4.85–4.35–2.17
U2 vacancy–5.84–4.22–3.54
At the surfaces containing one U vacancy site (U-1), the dissociative adsorption of H2O on the defective U3–xSi2{110} is found to be highly exothermic (Eads = −5.84 eV), which suggests that the defective U3–xSi2{001} favors dissociative H2O adsorption rather than molecular adsorption, similar to the findings on the defective U3–xSi2{001} and {111} surfaces. Bader population analysis reveals that a significant amount of charge (∼1.32e) was transferred to dissociated water species from the interacting surface species. This is not surprising since vacancies result in the formation of dangling bonds (reactive sites) due to the availability of electrons for transfer to water molecules on the surface.

3.4. Electronic Structure and Bonding Mechanism

Fundamental understanding of the nature of interactions between the H2O molecule and the U3Si2 surfaces and any adsorption-induced changes in the electronic structures of U3Si2 were gained through projected density of states (PDOS) analysis (Figure 15). The empty (filled) electronic states near the Fermi level can accept (donate) more electrons to enhance the reduction (oxidation) reactions. We observed various degrees of shifts in the PDOS, indicative of a change in the electronic structure due to electron transfer from U and Si ions to the water molecules leading to OH bond elongation; a step before deprotonation of the water molecule. The charge gained by the H2O molecule in different adsorption complexes is calculated to be in the range of 0.05–0.11e at the water–U3Si2{001} surface, 0.04–0.05e at the water–U3Si2{111} surface, and 0.15–0.25e at the water–U3Si2{110} surface (Table 3).

Figure 15

Figure 15. (Right) PDOS for the interacting surface U f-states before and after the adsorption of H2O at the (a) water–U3Si2{001} and (b) for interacting surface Si p-states at the water–U3Si2{110} interface, and for the interacting surface U f-states before and after the adsorption of H2O at the (c) water–U3Si2{111} interface. (Left) the corresponding optimized water–U3Si2 interfaces with U–O and Si–O bond lengths.

There is a possibility of surface oxide formation as shown by the U–O bond lengths of 2.620 and 2.555 Å and Si–O bond length of 2.120 Å on the {001}, {111}, and {110} surfaces, respectively. This suggests that water–U3Si2 oxidation would proceed by any of the following chemical reactions (50,51)
The reaction energy of the different possible reaction pathways is found to be highly exothermic, which suggest that they are thermodynamically favorable and feasible reactions.
We further analyzed the PDOS of the isolated water molecule in the adsorbed state at the various surfaces and compared to the gas-phase molecule (Figure 16). The DOS for the gas-phase H2O molecule is shown in Figure 16a where the molecular orbitals (MOs) are labeled as 2a1, 1b2, 3a1, and 1b1, respectively, while those for the lowest-energy adsorption configurations at the water–U3Si2{001}, {110}, and {111} interfaces are shown in Figure 16b–d, respectively. Due to the strong water–U3Si2 hybridization, electron transfer from the interacting surface U f-states, we observed a shift or disappearance of the 3a1 and 1b1 MOs around the Fermi level of water adsorbed at the {001} and {111} surfaces. At the water–U3Si2{110} interface, we observe a splitting of the 1b1 MO, coupled with a shift toward lower energy levels (Figure 16d), which signifies stabilization of the water molecule via physisorption. The reactivity of the surface is influenced largely by the ease of electron transfer between the surface and the adsorbate. Previous, work by Li et al. showed that the electronic states of the valence band of the surface play a key role during the surface–adsorbate interaction due to the hybridization energy between the bonding and antibonding adsorbate states and the metal valence-bands. (52)

Figure 16

Figure 16. DOS for H2O in the (a) free state and adsorbed in the lowest-energy geometry at the water–U3Si2 interfaces (b–d).

As the preadsorbed oxygen species enhance the O–H bond activation to dissociate, a necessary step for the formation of surface oxides, we have also characterized the electronic structures of the O-covered U3Si2 surfaces and showed that the U atoms to which atomic oxygen is bound become more positive (1.12e, 1.20e, and 1.61e on the {001}, {110}, and {111} surfaces, respectively) compared to the clean surface U charge of +0.98e {001}, 0.96e {110}, and 1.01e {111}, which, from the qU4+/qU6+ ratio, is enough to suggest that they have been oxidized from U4+ to U6+ (see Table S4 and Figure S5, for full details of O–U3Si2 in terms of bond lengths and charge transfer). We further investigated the nature of the interaction of the bonding between surface U and atomic oxygen by analyzing the projected density of states (PDOS) of the interacting surface U f-states and O p-state, as shown in Figure 17. We observed strong hybridization between the interacting U f-states and O p-state, which is due to the charge transfer from surface U ions into the adsorbed oxygen π orbital, in good agreement with previous studies of oxygen adsorption on metal surfaces (53)

Figure 17

Figure 17. Partial DOS projected on the interacting surface U f-states and O p-states for adsorbed atomic oxygen on the top-U site on (a) {001}, (b) {110}, and (c) {111} surfaces.

The work function gives a picture of the nature of such electronic interactions; hence, we have calculated the work function of the clean and adsorbate containing U3Si2 surfaces to characterize the level of difficulty for an electron transferring from the surface to the vacuum. The work function is the minimum energy needed to remove an electron from the bulk of a material through a surface to a point outside the material and can be written as
(12)
The potential in the vacuum region (Vvacuum) and the Fermi energy (EFermi) were derived from the same calculation. In practice, this is the energy required at 0 K to remove an electron from the Fermi level of the material to the vacuum potential. We determined the work function of the clean surface and further probed the effect of oxygen on the electronic states of the surface and oxygen-assisted dissociation of H2O.
The work function depends on the crystallographic direction as shown by the results in Table 7. The anisotropic nature of the work function comes from the crystallographic arrangement of the surface planes, which determines the spreading of the electronic charge into the vacuum. (54) Furthermore, the adsorption acts to smoothen the surface electric charge distribution that lowers the work function. The work function is such that Φ(H2O + O) < Φ(H2O) < Φdry due to the partial transfer of electron charge from the substrate to the adsorbate and the resulting adsorption-induced surface dipoles (Table 7). (55)
Table 7. Calculated Work Functions of the Dry (Φdry), Hydrated (ΦH2O), and Coadsorbed H2O + O (ΦH2O+O) U3Si2 Surfaces
surfaceΦdry (eV)ΦH2O (eV)ΦH2O+O (eV)
{001}3.093.052.86
{110}1.891.481.21
{111}2.862.301.96

4. Summary and Conclusions

Click to copy section linkSection link copied!

We have investigated the adsorption and dissociation reactions of H2O on clean, oxygen-covered and defective U3Si2{001}, {110}, and {100} surfaces using density functional theory calculations. The adsorption energetics and characteristics of molecular water adsorption were compared with full dissociative and mixed (molecular and dissociative) adsorption modes on the U3Si2 surfaces. We showed from our calculations that the major interactions between the adsorbing water molecules and the U3Si2 surfaces occur through oxygen and the surface U or Si site. Compared to the clean surfaces, preadsorbed O atoms are shown to enhance the activation of the O–H bonds of water and their subsequent dissociation reactions to form surface hydroxyl species, which are driven by significant charge transfer from the surface to the adsorbing species. We demonstrated that surface vacancy defects enhance the adsorption and dissociation of H2O compared to a stoichiometric defect-free surfaces. The molecular-level insights derived from this work provide a fundamental understanding of the adsorption processes and mechanisms of the early stage of oxidation of U3Si2 in the presence of oxygen and water and might open new avenues for the rational design of oxidation resistance of metallic fuels for nuclear reactors.

Supporting Information

Click to copy section linkSection link copied!

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcc.9b03076.

  • Mixed mode water adsorption configurations and adsorption energies on stoichiometric U3Si2{001}, {110}, and {111} surfaces; calculated adsorption energy (Eads), relevant bond lengths (d) of atomic (O) oxygen on {001}, {110}, and {111} U3Si2 surfaces (PDF)

Terms & Conditions

Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.

Author Information

Click to copy section linkSection link copied!

  • Corresponding Authors
    • Ericmoore Jossou - Department of Mechanical Engineering, College of Engineering, University of Saskatchewan, 57 Campus Drive, Saskatoon S7N 5A9, Saskatchewan, CanadaOrcidhttp://orcid.org/0000-0002-4792-760X Email: [email protected]
    • Jerzy Szpunar - Department of Mechanical Engineering, College of Engineering, University of Saskatchewan, 57 Campus Drive, Saskatoon S7N 5A9, Saskatchewan, Canada Email: [email protected]
  • Authors
    • Linu Malakkal - Department of Mechanical Engineering, College of Engineering, University of Saskatchewan, 57 Campus Drive, Saskatoon S7N 5A9, Saskatchewan, Canada
    • Nelson Y. Dzade - School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, U.K.Orcidhttp://orcid.org/0000-0001-7733-9473
    • Antoine Claisse - Westinghouse Electric Sweden AB, SE-721 63 Västerås, Sweden
    • Barbara Szpunar - Department of Physics and Engineering Physics, College of Art and Science, University of Saskatchewan, 116 Science Place, Saskatoon S7N 5E2, Saskatchewan, Canada
  • Notes
    The authors declare no competing financial interest.

Acknowledgments

Click to copy section linkSection link copied!

The authors gratefully acknowledge Compute Canada (WestGrid) and University of Saskatchewan’s Research Cluster (Plato) for a generous amount of CPU time. We also acknowledge the financial support of the Canadian National Science and Engineering Research Council (NSERC) and the University of Saskatchewan’s International Dean’s Scholarship. N.Y.D. acknowledges the U.K. Engineering and Physical Sciences Research Council (EPSRC) for funding (Grant No. EP/S001395/1).

References

Click to copy section linkSection link copied!

This article references 55 other publications.

  1. 1
    White, J. T.; Nelson, A. T.; Dunwoody, J. T.; Byler, D. D.; Safarik, D. J.; McClellan, K. J. Thermophysical Properties of U3Si2 to 1773 K. J. Nucl. Mater. 2015, 464, 275280,  DOI: 10.1016/j.jnucmat.2015.04.031
  2. 2
    Wood, E. S.; White, J. T.; Grote, C. J.; Nelson, A. T. U3Si2 Behavior in H2O: Part I, Fl Owing Steam and the Effect of Hydrogen. J. Nucl. Mater. 2018, 501, 404412,  DOI: 10.1016/j.jnucmat.2018.01.002
  3. 3
    Nelson, A. T.; Migdisov, A.; Wood, E. S.; Grote, C. J. U3Si2 Behavior in H2O Environments: Part II, Pressurized Water with Controlled Redox Chemistry. J. Nucl. Mater. 2018, 500, 8191,  DOI: 10.1016/j.jnucmat.2017.12.026
  4. 4
    Middleburgh, S. C.; Claisse, A.; Andersson, D. A.; Grimes, R. W.; Olsson, P.; Mašková, S. Solution of Hydrogen in Accident Tolerant Fuel Candidate Material: U3Si2. J. Nucl. Mater. 2018, 501, 234237,  DOI: 10.1016/j.jnucmat.2018.01.018
  5. 5
    Beeler, B.; Baskes, M.; Andersson, D.; Cooper, M. W. D.; Zhang, Y. A modified Embedded-Atom Method interatomic potential for uranium-silicide. J. Nucl. Mater. 2017, 495, 267276,  DOI: 10.1016/j.jnucmat.2017.08.025
  6. 6
    Noordhoek, M. J.; Besmann, T. M.; Andersson, D.; Middleburgh, S. C.; Chernatynskiy, A. Phase Equilibria in the U-Si System from First-Principles Calculations. J. Nucl. Mater. 2016, 479, 216223,  DOI: 10.1016/j.jnucmat.2016.07.006
  7. 7
    Wang, T.; Qiu, N.; Wen, X.; Tian, Y.; He, J.; Luo, K.; Zha, X.; Zhou, Y.; Huang, Q.; Lang, J. First-Principles Investigations on the Electronic Structures of U3Si2. J. Nucl. Mater. 2016, 469, 194199,  DOI: 10.1016/j.jnucmat.2015.11.060
  8. 8
    Andersson, D. A.; Liu, X.-Y.; Beeler, B.; Middleburgh, S. C.; Claisse, A.; Stanek, C. R. Density Functional Theory Calculations of Self- and Xe Diffusion in U3Si2. J. Nucl. Mater. 2019, 515, 312325,  DOI: 10.1016/j.jnucmat.2018.12.021
  9. 9
    Beeler, B.; Baskes, M.; Andersson, D.; Cooper, M. W. D.; Zhang, Y. Molecular Dynamics Investigation of Grain Boundaries and Surfaces in U3Si2. J. Nucl. Mater. 2019, 514, 290298,  DOI: 10.1016/j.jnucmat.2018.12.008
  10. 10
    Middleburgh, S. C.; Grimes, R. W.; Lahoda, E. J.; Stanek, C. R.; Andersson, D. A. Non-Stoichiometry in U3Si2. J. Nucl. Mater. 2016, 482, 300305,  DOI: 10.1016/j.jnucmat.2016.10.016
  11. 11
    Remschnig, K.; Le Bihan, T.; Noël, H.; Rogl, P. Structural Chemistry and Magnetic Behavior of Binary Uranium Silicides. J. Solid State Chem. 1992, 97, 391399,  DOI: 10.1016/0022-4596(92)90048-Z
  12. 12
    Bo, T.; Lan, J.; Zhao, Y.; Zhang, Y.; He, C.; Chai, Z.; Shi, W. Surface Science Surface Properties of NpO2 and Water Reacting with Stoichiometric and Atomistic Thermodynamics. Surf. Sci. 2016, 644, 153164,  DOI: 10.1016/j.susc.2015.09.016
  13. 13
    Bo, T.; Lan, J.-H.; Zhao, Y.-L.; Zhang, Y.-J.; He, C.-H.; Chai, Z.-F.; Shi, W.-Q. Surface Properties of NpO2 and Water Reacting with Stoichiometric and Reduced NpO2 (111), (110), and (100) Surfaces from Ab Initio Atomistic Thermodynamics. Surf. Sci. 2016, 644, 153164,  DOI: 10.1016/j.susc.2015.09.016
  14. 14
    Jossou, E.; Eduok, U.; Dzade, N. Y.; Szpunar, B.; Szpunar, J. A. Oxidation Behaviour of U3Si2: An Experimental and First Principles Investigation. Phys. Chem. Chem. Phys. 2018, 20, 47084720,  DOI: 10.1039/C7CP07154J
  15. 15
    Molinari, M.; Parker, S. C.; Sayle, D. C.; Islam, M. S. Water Adsorption and Its Effect on the Stability of Low Index Stoichiometric and Reduced Surfaces of Ceria. J. Phys. Chem. C 2012, 116, 70737083,  DOI: 10.1021/jp300576b
  16. 16
    Bo, T.; Lan, J.-H.; Wang, C.-Z.; Zhao, Y.-L.; He, C.-H.; Zhang, Y.-J.; Chai, Z.-F.; Shi, W.-Q. First-Principles Study of Water Reaction and H2 Formation on UO2 (111) and (110) Single Crystal Surfaces. J. Phys. Chem. C 2014, 118, 2193521944,  DOI: 10.1021/jp503614f
  17. 17
    Fan, J.; Li, C.; Zhao, J.; Shan, Y.; Xu, H. The Enhancement of Surface Reactivity on CeO2 (111) Mediated by Subsurface Oxygen Vacancies. J. Phys. Chem. C 2016, 120, 2791727924,  DOI: 10.1021/acs.jpcc.6b07650
  18. 18
    Zhukovskii, Y.; Bocharov, D.; Gryaznov, D.; Kotomin, E. First Principles Simulations on Surface Properties and Oxidation of Nitride Nuclear Fuels; Advances in Nuclear Fuel; IntechOpen, 2012.
  19. 19
    Wellington, J. P. W.; Tegner, B. E.; Collard, J.; Kerridge, A.; Kaltsoyannis, N. Oxygen Vacancy Formation and Water Adsorption on Reduced AnO2{111}, {110}, and {100} Surfaces (An = U, Pu): A Computational Study. J. Phys. Chem. C 2018, 122, 71497165,  DOI: 10.1021/acs.jpcc.7b11512
  20. 20
    Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. 1964, 136, B864B871,  DOI: 10.1103/PhysRev.136.B864
  21. 21
    Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133A1138,  DOI: 10.1103/PhysRev.140.A1133
  22. 22
    Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G. L.; Cococcioni, M.; Dabo, I. QUANTUM ESPRESSO: A Modular and Open-Source Software Project for Quantum Simulations of Materials. J. Phys.: Condens. Matter 2009, 21, 395502  DOI: 10.1088/0953-8984/21/39/395502
  23. 23
    Wu, Z.; Cohen, R. E. More Accurate Generalized Gradient Approximation for Solids. Phys. Rev. B 2006, 73, 235116  DOI: 10.1103/PhysRevB.73.235116
  24. 24
    Liechtenstein, A. I.; Anisimov, V. I.; Zaanen, J. Density-Functional Theory and Strong Interactions: Orbital Ordering in Mott-Hubbard Insulators. Phys. Rev. B 1995, 52, R5467R5470,  DOI: 10.1103/PhysRevB.52.R5467
  25. 25
    Methfessel, M.; Paxton, A. T. High-Precision Sampling for Brillouin-Zone Integration in Metals. Phys. Rev. B 1989, 40, 36163621,  DOI: 10.1103/PhysRevB.40.3616
  26. 26
    Head, J. D.; Zerner, M. C. A Broyden—Fletcher—Goldfarb—Shanno Optimization Procedure for Molecular Geometries. Chem. Phys. Lett. 1985, 122, 264270,  DOI: 10.1016/0009-2614(85)80574-1
  27. 27
    Momma, K.; Izumi, F. VESTA: A Three-Dimensional Visualization System for Electronic and Structural Analysis. J. Appl. Crystallogr. 2008, 41, 653658,  DOI: 10.1107/S0021889808012016
  28. 28
    Momma, K.; Izumi, F. VESTA3 for Three-Dimensional Visualization of Crystal, Volumetric and Morphology Data. J. Appl. Crystallogr. 2011, 44, 12721276,  DOI: 10.1107/S0021889811038970
  29. 29
    Watson, G. W.; Kelsey, E. T.; de Leeuw, N. H.; Harris, D. J.; Parker, S. C. Atomistic Simulation of Dislocations, Surfaces and Interfaces in MgO. J. Chem. Soc., Faraday Trans. 1996, 92, 433438,  DOI: 10.1039/ft9969200433
  30. 30
    Bader, R. F. W. Atoms in Molecules. Acc. Chem. Res. 1985, 18, 915,  DOI: 10.1021/ar00109a003
  31. 31
    Andersson, A. D. Density Functional Theory Calculations of Defect and Fission Gas Properties in U-Si Fuels; Office of Scientific and Technical Information, 2016.
  32. 32
    Li, R.; Zhang, P.; Zhang, C.; Huang, X.; Zhao, J. Vacancy Trapping Mechanism for Multiple Helium in Monovacancy and Small Void of Vanadium Solid. J. Nucl. Mater. 2013, 440, 557561,  DOI: 10.1016/j.jnucmat.2013.03.068
  33. 33
    Lide, D. R. Handbook of Chemistry and Physics, 82nd ed.; CRC Press: Boca Raton, FL, 2001.
  34. 34
    Dzade, N. Y.; Roldan, A.; de Leeuw, N. H. DFT-D2 Study of the Adsorption and Dissociation of Water on Clean and Oxygen-Covered {001} and {011} Surfaces of Mackinawite (FeS). J. Phys. Chem. C 2016, 120, 2144121450,  DOI: 10.1021/acs.jpcc.6b06122
  35. 35
    Herzberg, G. Molecular Spectra and Molecular Structure. II. Infrared and Raman Spectra of Polyatomic Molecules; Lancaster Press: New York, 1946; p 365.
  36. 36
    Bo, T.; Lan, J.; Zhang, Y.; Zhao, Y.; He, C.; Chai, Z.; Shi, W. Surface of Uranium Mononitride: Energetics And. Phys. Chem. Chem. Phys. 2016, 18, 1325513266,  DOI: 10.1039/C6CP01175F
  37. 37
    Bo, T.; Lan, J.; Zhao, Y.; Zhang, Y.; He, C.; Chai, Z.; Shi, W. First-Principles Study of Water Adsorption and Dissociation on the UO2. J. Nucl. Mater. 2014, 454, 446454,  DOI: 10.1016/j.jnucmat.2014.09.001
  38. 38
    Shapovalov, V.; Truong, T. N. Ab Initio Study of Water Adsorption on α-Al2O3 (0001) Crystal Surface. J. Phys. Chem. B 2000, 104, 98599863,  DOI: 10.1021/jp001399g
  39. 39
    Pang, Z.; Duerrbeck, S.; Kha, C.; Bertel, E.; Somorjai, G. A.; Salmeron, M. Adsorption and Reactions of Water on Oxygen-Precovered Cu(110). J. Phys. Chem. C 2016, 120, 92189222,  DOI: 10.1021/acs.jpcc.6b00769
  40. 40
    Chen, L.; Lu, J.; Liu, P.; Gao, L.; Liu, Y.; Xiong, F.; Qiu, S.; Qiu, X.; Guo, Y.; Chen, X. Dissociation and Charge Transfer of H2O on Cu(110) Probed in Real Time Using Ion Scattering Spectroscopy. Langmuir 2016, 32, 1204712055,  DOI: 10.1021/acs.langmuir.6b03516
  41. 41
    Taylor, C. D.; Lookman, T.; Lillard, R. S. Ab Initio Calculations of the Uranium–Hydrogen System: Thermodynamics, Hydrogen Saturation of a -U and Phase-Transformation to UH3. Acta Mater. 2010, 58, 10451055,  DOI: 10.1016/j.actamat.2009.10.021
  42. 42
    Bo, T.; Lan, J.-H.; Zhang, Y.-J.; Zhao, Y.-L.; He, C.-H.; Chai, Z.-F.; Shi, W.-Q. Adsorption and Dissociation of H2O on the (001) Surface of Uranium Mononitride: Energetics and Mechanism from First-Principles Investigation. Phys. Chem. Chem. Phys. 2016, 18, 1325513266,  DOI: 10.1039/C6CP01175F
  43. 43
    Tian, X.; Wang, H.; Xiao, H.; Gao, T. Adsorption of Water on UO2 (111) Surface: Density Functional Theory Calculations. Comput. Mater. Sci. 2014, 91, 364371,  DOI: 10.1016/j.commatsci.2014.05.009
  44. 44
    Tegner, B. E.; Molinari, M.; Kerridge, A.; Parker, 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, 16751682,  DOI: 10.1021/acs.jpcc.6b10986
  45. 45
    Schoenes, J. Optical Properties and Electronic Structure of UO2. J. Appl. Phys. 1978, 49, 14631465,  DOI: 10.1063/1.324978
  46. 46
    Wang, J.; Ewing, R. C.; Becker, U. Electronic Structure and Stability of Hyperstoichiometric UO2+x under Pressure. Phys. Rev. B 2013, 88, 024109  DOI: 10.1103/PhysRevB.88.024109
  47. 47
    Henderson, M. A. The Interaction of Water with Solid Surfaces: Fundamental Aspects Revisited. Surf. Sci. Rep. 2002, 46, 1308,  DOI: 10.1016/S0167-5729(01)00020-6
  48. 48
    Ammon, C.; Bayer, A.; Steinrück, H.-P.; Held, G. Low-Temperature Partial Dissociation of Water on Cu(110). Chem. Phys. Lett. 2003, 377, 163169,  DOI: 10.1016/S0009-2614(03)01127-8
  49. 49
    Jiang, Z.; Fang, T. Dissociation Mechanism of H2O on Clean and Oxygen-Covered Cu(111) Surfaces: A Theoretical Study. Vaccum 2016, 128, 252258,  DOI: 10.1016/j.vacuum.2016.03.030
  50. 50
    Johnson, K.; Ström, V.; Wallenius, J.; Lopes, D. A. Oxidation of Accident Tolerant Fuel Candidates. J. Nucl. Sci. Technol. 2017, 54, 280286,  DOI: 10.1080/00223131.2016.1262297
  51. 51
    Johnson, K. D.; Wallenius, J.; Jolkkonen, M.; Claisse, A. Spark Plasma Sintering and Porosity Studies of Uranium Nitride. J. Nucl. Mater. 2016, 473, 1317,  DOI: 10.1016/j.jnucmat.2016.01.037
  52. 52
    Li, J.; Zhu, S.; Li, H.; Oguzie, E. E.; Li, Y.; Wang, F. Bonding Nature of Monomeric H2O on Pd: Orbital Cooperation and Competition. J. Phys. Chem. C 2009, 19311938,  DOI: 10.1021/jp809595y
  53. 53
    Akbay, T.; Staykov, A.; Druce, J.; Téllez, H.; Ishihara, T.; Kilner, J. A. The interaction of molecular oxygen on LaO terminated surfaces of La2NiO4. J. Mater. Chem. A 2016, 4, 13113,  DOI: 10.1039/c6ta02715f
  54. 54
    Fall, C. Ab Initio Study of the Work Functions of Elemental Metal Crystals; EPFL PP: Lausanne.
  55. 55
    Zhang, W.; Liu, L.; Wan, L.; Liu, L.; Cao, L.; Xu, F.; Zhao, J.; Wu, Z. Electronic Structures of Bare and Terephthalic Acid Adsorbed TiO2(110)-(1 × 2) Reconstructed Surfaces: Origin and Reactivity of the Band Gap States. Phys. Chem. Chem. Phys. 2015, 17, 2014420153,  DOI: 10.1039/C5CP01298H

Cited By

Click to copy section linkSection link copied!

This article is cited by 12 publications.

  1. Mohammad Asif, Faisal Khan, Kelly Hawboldt. Atomistic Analysis of the Microbial Influence on the Adsorption Characteristic of Sulfur, Hydrogen, and SO4on Iron Surfaces. Industrial & Engineering Chemistry Research 2023, 62 (48) , 20777-20788. https://doi.org/10.1021/acs.iecr.3c02592
  2. Mathew John Haskew, Benjamin Deacon, Chin Weng Yong, John George Hardy, Samuel Thomas Murphy. Atomistic Simulation of Water Incorporation and Mobility in Bombyx mori Silk Fibroin. ACS Omega 2021, 6 (51) , 35494-35504. https://doi.org/10.1021/acsomega.1c05019
  3. Linu Malakkal, Ericmoore Jossou, Jayangani Inoka Ranasinghe, Barbara Szpunar, Jerzy Szpunar. Density Functional Theory Study of Oxygen Adsorption and Dissociation on Lower Miller Index Surfaces of ThN. The Journal of Physical Chemistry C 2020, 124 (45) , 24849-24860. https://doi.org/10.1021/acs.jpcc.0c07770
  4. Wenji Wang. Physisorbed State Regulates the Dissociation Mechanism of H2O on Ni(100). The Journal of Physical Chemistry A 2020, 124 (42) , 8724-8732. https://doi.org/10.1021/acs.jpca.0c06130
  5. Chenglong Qin, Yushu Yu, Zihan Xu, Jiguang Du, Liang Zhao, Gang Jiang. Investigations on mechanical, thermodynamic and surface properties of U–Si alloys. Vacuum 2023, 216 , 112417. https://doi.org/10.1016/j.vacuum.2023.112417
  6. Anum Shahid Malik, Lisa A. Fredin. Untangling product selectivity on clean low index rutile TiO 2 surfaces using first-principles calculations. Physical Chemistry Chemical Physics 2023, 25 (3) , 2203-2211. https://doi.org/10.1039/D2CP04939B
  7. Shannon E. Cooney, Alex A. Fertig, Madeleine R. Buisch, William W. Brennessel, Ellen M. Matson. Coordination-induced bond weakening of water at the surface of an oxygen-deficient polyoxovanadate cluster. Chemical Science 2022, 13 (43) , 12726-12737. https://doi.org/10.1039/D2SC04843D
  8. Jiatu Liu, Patrick A. Burr, Joshua T. White, Vanessa K. Peterson, Pranesh Dayal, Christopher Baldwin, Deborah Wakeham, Daniel J. Gregg, Elizabeth S. Sooby, Edward G. Obbard. Structural and phase evolution in U3Si2 during steam corrosion. Corrosion Science 2022, 204 , 110373. https://doi.org/10.1016/j.corsci.2022.110373
  9. Xiaoyong Yang, Pavel A. Korzhavyi, Yang Liu, Qianglin Wei, Temirlan R. Arslanov, John P.A. Wärnå, Yu Yang, Ping Zhang. Structural, electronic, mechanical and thermodynamic properties of U–Si intermetallic compounds: A comprehensive first principles calculations. Progress in Nuclear Energy 2022, 148 , 104229. https://doi.org/10.1016/j.pnucene.2022.104229
  10. Faisal Khan, Mohmmad Asif, Kelly Hawboldt, Shams Anwar. Sulfide (H2s) Corrosion Modelling of Iron (Fe) Using Molecular Modeling Approach. SSRN Electronic Journal 2022, 31 https://doi.org/10.2139/ssrn.4003102
  11. Jiatu Liu, Patrick A. Burr, Joshua T. White, Vanessa K. Peterson, Pranesh Dayal, Christopher Baldwin, Deborah Wakeham, Daniel J. Gregg, Elizabeth S. Sooby, Edward G. Obbard. Structural and Phase Evolution in U3si2 During Steam Corrosion. SSRN Electronic Journal 2022, 484 https://doi.org/10.2139/ssrn.4050274
  12. Adrian Gonzales, Jennifer K. Watkins, Adrian R. Wagner, Brian J. Jaques, Elizabeth S. Sooby. Challenges and opportunities to alloyed and composite fuel architectures to mitigate high uranium density fuel oxidation: uranium silicide. Journal of Nuclear Materials 2021, 553 , 153026. https://doi.org/10.1016/j.jnucmat.2021.153026
Open PDF

The Journal of Physical Chemistry C

Cite this: J. Phys. Chem. C 2019, 123, 32, 19453–19467
Click to copy citationCitation copied!
https://doi.org/10.1021/acs.jpcc.9b03076
Published July 11, 2019

Copyright © 2019 American Chemical Society. This publication is licensed under CC-BY.

Article Views

2652

Altmetric

-

Citations

Learn about these metrics

Article Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.

Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.

The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated.

  • Abstract

    Figure 1

    Figure 1. Unit cell of ordered uranium silicide containing two formula units of U3Si2 (10 atoms) (color scheme: U = gray and Si = blue).

    Figure 2

    Figure 2. Optimized surface geometry of (a) U3Si2{001}, (b) U3Si2{110}, and (c) U3Si2{111} with Si and U1 surface vacancy represented by a rectangular red box in the second and third rows respectively (color scheme: U = gray and Si = blue).

    Figure 3

    Figure 3. Calculated surface energies of the {001}, {110}, and {111} surfaces of U3Si2 as functions of the change in silicon chemical potential μSi.

    Figure 4

    Figure 4. Top and side views of the relaxed adsorption structures of molecular H2O adsorbed at: (a) O–U, (b) H–Si, and (c) O–Si (color scheme: U = gray, Si = blue, H = white, and O = red).

    Figure 5

    Figure 5. Top and side views of the relaxed adsorption structures of dissociated water configuration with (a) OH–Si: H on adj. U, (b) OH–Si: H on Si, and (c) OH–Si: H on U, on the U3Si2{001}-Si terminated surface (color scheme: U = gray, Si = blue, H = white, and O = red).

    Figure 6

    Figure 6. Top and side views of the relaxed adsorption structures of molecular oxygen adsorbed at (a) O–Si, (b) O–U, and (c) H–Si on the U3Si2{110} surface (color scheme: U = gray, Si = blue, H = white, and O = red).

    Figure 7

    Figure 7. Top and side views of the relaxed adsorption structures of dissociated H2O at (a) OH–Si: H on U and (b) OH–U: H on Si on the U3Si2{110} surface (color scheme: U = gray, Si = blue, H = white, and O = red).

    Figure 8

    Figure 8. Top and side views of the relaxed adsorption structures of molecular H2O adsorbed at (a) the O–U site, (b) the H–Si site, and (c) the O–Si site on the U3Si2{111} surface (color scheme: U = gray, Si = blue, H = white, and O = red).

    Figure 9

    Figure 9. Top and side views of the relaxed adsorption structures dissociated water configuration with (a) O–U: H and (b) O–Si: H, on the U3Si2{111}-Si terminated surface (color scheme: U = gray, Si = blue, H = white, and O = red).

    Figure 10

    Figure 10. Top and side view of the relaxed adsorption structures of mixed (dissociative to molecular) 3:1 adsorption of H2O on (a) {001} and fully dissociative adsorption of H2O on (b) the {110} and (c) {111} U3Si2 surface (color scheme: U = gray, Si = blue, H= white, and O = red).

    Figure 11

    Figure 11. Top and side views of the relaxed adsorption structures of molecular H2O (a) H–U: O on U, (b) Ow–Si: O on Si, (c) Ow–Si: O on U, and (d) Ow–U: O on Si U3Si2{001} surface (color scheme: U = gray, Si = blue, H = white, and O = red).

    Figure 12

    Figure 12. Top and side views of the relaxed adsorption structures of molecular H2O (a) Ow–Si: O on the U bridge, (b) Ow–Si: O on U, and (c) Hw–Si: O on U on the U3Si2{110} surface (color scheme: U = gray, Si = blue, H = white, and O = red).

    Figure 13

    Figure 13. Top and side views of the relaxed adsorption structures of molecular H2O (a) Ow–U: O and (b) Ow–Si: O on the U3Si2{111} surface (color scheme: U = gray, Si = blue, H = white, and O = red).

    Figure 14

    Figure 14. Top and side views of the relaxed adsorption structures of Si and U vacancy-assisted molecular H2O adsorption (a) H2O–U3Si2–x {001}, (b) H2O–U3Si2–x {110}, (c) H2O–U3Si2–x {111}, (d) H2O–U3–xSi2{001}, (e) H2O–U3–xSi2{110}, and (f) H2O–U3–xSi2{111} surface (color scheme: U = gray, Si = blue, H = white, and O = red).

    Figure 15

    Figure 15. (Right) PDOS for the interacting surface U f-states before and after the adsorption of H2O at the (a) water–U3Si2{001} and (b) for interacting surface Si p-states at the water–U3Si2{110} interface, and for the interacting surface U f-states before and after the adsorption of H2O at the (c) water–U3Si2{111} interface. (Left) the corresponding optimized water–U3Si2 interfaces with U–O and Si–O bond lengths.

    Figure 16

    Figure 16. DOS for H2O in the (a) free state and adsorbed in the lowest-energy geometry at the water–U3Si2 interfaces (b–d).

    Figure 17

    Figure 17. Partial DOS projected on the interacting surface U f-states and O p-states for adsorbed atomic oxygen on the top-U site on (a) {001}, (b) {110}, and (c) {111} surfaces.

  • References


    This article references 55 other publications.

    1. 1
      White, J. T.; Nelson, A. T.; Dunwoody, J. T.; Byler, D. D.; Safarik, D. J.; McClellan, K. J. Thermophysical Properties of U3Si2 to 1773 K. J. Nucl. Mater. 2015, 464, 275280,  DOI: 10.1016/j.jnucmat.2015.04.031
    2. 2
      Wood, E. S.; White, J. T.; Grote, C. J.; Nelson, A. T. U3Si2 Behavior in H2O: Part I, Fl Owing Steam and the Effect of Hydrogen. J. Nucl. Mater. 2018, 501, 404412,  DOI: 10.1016/j.jnucmat.2018.01.002
    3. 3
      Nelson, A. T.; Migdisov, A.; Wood, E. S.; Grote, C. J. U3Si2 Behavior in H2O Environments: Part II, Pressurized Water with Controlled Redox Chemistry. J. Nucl. Mater. 2018, 500, 8191,  DOI: 10.1016/j.jnucmat.2017.12.026
    4. 4
      Middleburgh, S. C.; Claisse, A.; Andersson, D. A.; Grimes, R. W.; Olsson, P.; Mašková, S. Solution of Hydrogen in Accident Tolerant Fuel Candidate Material: U3Si2. J. Nucl. Mater. 2018, 501, 234237,  DOI: 10.1016/j.jnucmat.2018.01.018
    5. 5
      Beeler, B.; Baskes, M.; Andersson, D.; Cooper, M. W. D.; Zhang, Y. A modified Embedded-Atom Method interatomic potential for uranium-silicide. J. Nucl. Mater. 2017, 495, 267276,  DOI: 10.1016/j.jnucmat.2017.08.025
    6. 6
      Noordhoek, M. J.; Besmann, T. M.; Andersson, D.; Middleburgh, S. C.; Chernatynskiy, A. Phase Equilibria in the U-Si System from First-Principles Calculations. J. Nucl. Mater. 2016, 479, 216223,  DOI: 10.1016/j.jnucmat.2016.07.006
    7. 7
      Wang, T.; Qiu, N.; Wen, X.; Tian, Y.; He, J.; Luo, K.; Zha, X.; Zhou, Y.; Huang, Q.; Lang, J. First-Principles Investigations on the Electronic Structures of U3Si2. J. Nucl. Mater. 2016, 469, 194199,  DOI: 10.1016/j.jnucmat.2015.11.060
    8. 8
      Andersson, D. A.; Liu, X.-Y.; Beeler, B.; Middleburgh, S. C.; Claisse, A.; Stanek, C. R. Density Functional Theory Calculations of Self- and Xe Diffusion in U3Si2. J. Nucl. Mater. 2019, 515, 312325,  DOI: 10.1016/j.jnucmat.2018.12.021
    9. 9
      Beeler, B.; Baskes, M.; Andersson, D.; Cooper, M. W. D.; Zhang, Y. Molecular Dynamics Investigation of Grain Boundaries and Surfaces in U3Si2. J. Nucl. Mater. 2019, 514, 290298,  DOI: 10.1016/j.jnucmat.2018.12.008
    10. 10
      Middleburgh, S. C.; Grimes, R. W.; Lahoda, E. J.; Stanek, C. R.; Andersson, D. A. Non-Stoichiometry in U3Si2. J. Nucl. Mater. 2016, 482, 300305,  DOI: 10.1016/j.jnucmat.2016.10.016
    11. 11
      Remschnig, K.; Le Bihan, T.; Noël, H.; Rogl, P. Structural Chemistry and Magnetic Behavior of Binary Uranium Silicides. J. Solid State Chem. 1992, 97, 391399,  DOI: 10.1016/0022-4596(92)90048-Z
    12. 12
      Bo, T.; Lan, J.; Zhao, Y.; Zhang, Y.; He, C.; Chai, Z.; Shi, W. Surface Science Surface Properties of NpO2 and Water Reacting with Stoichiometric and Atomistic Thermodynamics. Surf. Sci. 2016, 644, 153164,  DOI: 10.1016/j.susc.2015.09.016
    13. 13
      Bo, T.; Lan, J.-H.; Zhao, Y.-L.; Zhang, Y.-J.; He, C.-H.; Chai, Z.-F.; Shi, W.-Q. Surface Properties of NpO2 and Water Reacting with Stoichiometric and Reduced NpO2 (111), (110), and (100) Surfaces from Ab Initio Atomistic Thermodynamics. Surf. Sci. 2016, 644, 153164,  DOI: 10.1016/j.susc.2015.09.016
    14. 14
      Jossou, E.; Eduok, U.; Dzade, N. Y.; Szpunar, B.; Szpunar, J. A. Oxidation Behaviour of U3Si2: An Experimental and First Principles Investigation. Phys. Chem. Chem. Phys. 2018, 20, 47084720,  DOI: 10.1039/C7CP07154J
    15. 15
      Molinari, M.; Parker, S. C.; Sayle, D. C.; Islam, M. S. Water Adsorption and Its Effect on the Stability of Low Index Stoichiometric and Reduced Surfaces of Ceria. J. Phys. Chem. C 2012, 116, 70737083,  DOI: 10.1021/jp300576b
    16. 16
      Bo, T.; Lan, J.-H.; Wang, C.-Z.; Zhao, Y.-L.; He, C.-H.; Zhang, Y.-J.; Chai, Z.-F.; Shi, W.-Q. First-Principles Study of Water Reaction and H2 Formation on UO2 (111) and (110) Single Crystal Surfaces. J. Phys. Chem. C 2014, 118, 2193521944,  DOI: 10.1021/jp503614f
    17. 17
      Fan, J.; Li, C.; Zhao, J.; Shan, Y.; Xu, H. The Enhancement of Surface Reactivity on CeO2 (111) Mediated by Subsurface Oxygen Vacancies. J. Phys. Chem. C 2016, 120, 2791727924,  DOI: 10.1021/acs.jpcc.6b07650
    18. 18
      Zhukovskii, Y.; Bocharov, D.; Gryaznov, D.; Kotomin, E. First Principles Simulations on Surface Properties and Oxidation of Nitride Nuclear Fuels; Advances in Nuclear Fuel; IntechOpen, 2012.
    19. 19
      Wellington, J. P. W.; Tegner, B. E.; Collard, J.; Kerridge, A.; Kaltsoyannis, N. Oxygen Vacancy Formation and Water Adsorption on Reduced AnO2{111}, {110}, and {100} Surfaces (An = U, Pu): A Computational Study. J. Phys. Chem. C 2018, 122, 71497165,  DOI: 10.1021/acs.jpcc.7b11512
    20. 20
      Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. 1964, 136, B864B871,  DOI: 10.1103/PhysRev.136.B864
    21. 21
      Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133A1138,  DOI: 10.1103/PhysRev.140.A1133
    22. 22
      Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G. L.; Cococcioni, M.; Dabo, I. QUANTUM ESPRESSO: A Modular and Open-Source Software Project for Quantum Simulations of Materials. J. Phys.: Condens. Matter 2009, 21, 395502  DOI: 10.1088/0953-8984/21/39/395502
    23. 23
      Wu, Z.; Cohen, R. E. More Accurate Generalized Gradient Approximation for Solids. Phys. Rev. B 2006, 73, 235116  DOI: 10.1103/PhysRevB.73.235116
    24. 24
      Liechtenstein, A. I.; Anisimov, V. I.; Zaanen, J. Density-Functional Theory and Strong Interactions: Orbital Ordering in Mott-Hubbard Insulators. Phys. Rev. B 1995, 52, R5467R5470,  DOI: 10.1103/PhysRevB.52.R5467
    25. 25
      Methfessel, M.; Paxton, A. T. High-Precision Sampling for Brillouin-Zone Integration in Metals. Phys. Rev. B 1989, 40, 36163621,  DOI: 10.1103/PhysRevB.40.3616
    26. 26
      Head, J. D.; Zerner, M. C. A Broyden—Fletcher—Goldfarb—Shanno Optimization Procedure for Molecular Geometries. Chem. Phys. Lett. 1985, 122, 264270,  DOI: 10.1016/0009-2614(85)80574-1
    27. 27
      Momma, K.; Izumi, F. VESTA: A Three-Dimensional Visualization System for Electronic and Structural Analysis. J. Appl. Crystallogr. 2008, 41, 653658,  DOI: 10.1107/S0021889808012016
    28. 28
      Momma, K.; Izumi, F. VESTA3 for Three-Dimensional Visualization of Crystal, Volumetric and Morphology Data. J. Appl. Crystallogr. 2011, 44, 12721276,  DOI: 10.1107/S0021889811038970
    29. 29
      Watson, G. W.; Kelsey, E. T.; de Leeuw, N. H.; Harris, D. J.; Parker, S. C. Atomistic Simulation of Dislocations, Surfaces and Interfaces in MgO. J. Chem. Soc., Faraday Trans. 1996, 92, 433438,  DOI: 10.1039/ft9969200433
    30. 30
      Bader, R. F. W. Atoms in Molecules. Acc. Chem. Res. 1985, 18, 915,  DOI: 10.1021/ar00109a003
    31. 31
      Andersson, A. D. Density Functional Theory Calculations of Defect and Fission Gas Properties in U-Si Fuels; Office of Scientific and Technical Information, 2016.
    32. 32
      Li, R.; Zhang, P.; Zhang, C.; Huang, X.; Zhao, J. Vacancy Trapping Mechanism for Multiple Helium in Monovacancy and Small Void of Vanadium Solid. J. Nucl. Mater. 2013, 440, 557561,  DOI: 10.1016/j.jnucmat.2013.03.068
    33. 33
      Lide, D. R. Handbook of Chemistry and Physics, 82nd ed.; CRC Press: Boca Raton, FL, 2001.
    34. 34
      Dzade, N. Y.; Roldan, A.; de Leeuw, N. H. DFT-D2 Study of the Adsorption and Dissociation of Water on Clean and Oxygen-Covered {001} and {011} Surfaces of Mackinawite (FeS). J. Phys. Chem. C 2016, 120, 2144121450,  DOI: 10.1021/acs.jpcc.6b06122
    35. 35
      Herzberg, G. Molecular Spectra and Molecular Structure. II. Infrared and Raman Spectra of Polyatomic Molecules; Lancaster Press: New York, 1946; p 365.
    36. 36
      Bo, T.; Lan, J.; Zhang, Y.; Zhao, Y.; He, C.; Chai, Z.; Shi, W. Surface of Uranium Mononitride: Energetics And. Phys. Chem. Chem. Phys. 2016, 18, 1325513266,  DOI: 10.1039/C6CP01175F
    37. 37
      Bo, T.; Lan, J.; Zhao, Y.; Zhang, Y.; He, C.; Chai, Z.; Shi, W. First-Principles Study of Water Adsorption and Dissociation on the UO2. J. Nucl. Mater. 2014, 454, 446454,  DOI: 10.1016/j.jnucmat.2014.09.001
    38. 38
      Shapovalov, V.; Truong, T. N. Ab Initio Study of Water Adsorption on α-Al2O3 (0001) Crystal Surface. J. Phys. Chem. B 2000, 104, 98599863,  DOI: 10.1021/jp001399g
    39. 39
      Pang, Z.; Duerrbeck, S.; Kha, C.; Bertel, E.; Somorjai, G. A.; Salmeron, M. Adsorption and Reactions of Water on Oxygen-Precovered Cu(110). J. Phys. Chem. C 2016, 120, 92189222,  DOI: 10.1021/acs.jpcc.6b00769
    40. 40
      Chen, L.; Lu, J.; Liu, P.; Gao, L.; Liu, Y.; Xiong, F.; Qiu, S.; Qiu, X.; Guo, Y.; Chen, X. Dissociation and Charge Transfer of H2O on Cu(110) Probed in Real Time Using Ion Scattering Spectroscopy. Langmuir 2016, 32, 1204712055,  DOI: 10.1021/acs.langmuir.6b03516
    41. 41
      Taylor, C. D.; Lookman, T.; Lillard, R. S. Ab Initio Calculations of the Uranium–Hydrogen System: Thermodynamics, Hydrogen Saturation of a -U and Phase-Transformation to UH3. Acta Mater. 2010, 58, 10451055,  DOI: 10.1016/j.actamat.2009.10.021
    42. 42
      Bo, T.; Lan, J.-H.; Zhang, Y.-J.; Zhao, Y.-L.; He, C.-H.; Chai, Z.-F.; Shi, W.-Q. Adsorption and Dissociation of H2O on the (001) Surface of Uranium Mononitride: Energetics and Mechanism from First-Principles Investigation. Phys. Chem. Chem. Phys. 2016, 18, 1325513266,  DOI: 10.1039/C6CP01175F
    43. 43
      Tian, X.; Wang, H.; Xiao, H.; Gao, T. Adsorption of Water on UO2 (111) Surface: Density Functional Theory Calculations. Comput. Mater. Sci. 2014, 91, 364371,  DOI: 10.1016/j.commatsci.2014.05.009
    44. 44
      Tegner, B. E.; Molinari, M.; Kerridge, A.; Parker, 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, 16751682,  DOI: 10.1021/acs.jpcc.6b10986
    45. 45
      Schoenes, J. Optical Properties and Electronic Structure of UO2. J. Appl. Phys. 1978, 49, 14631465,  DOI: 10.1063/1.324978
    46. 46
      Wang, J.; Ewing, R. C.; Becker, U. Electronic Structure and Stability of Hyperstoichiometric UO2+x under Pressure. Phys. Rev. B 2013, 88, 024109  DOI: 10.1103/PhysRevB.88.024109
    47. 47
      Henderson, M. A. The Interaction of Water with Solid Surfaces: Fundamental Aspects Revisited. Surf. Sci. Rep. 2002, 46, 1308,  DOI: 10.1016/S0167-5729(01)00020-6
    48. 48
      Ammon, C.; Bayer, A.; Steinrück, H.-P.; Held, G. Low-Temperature Partial Dissociation of Water on Cu(110). Chem. Phys. Lett. 2003, 377, 163169,  DOI: 10.1016/S0009-2614(03)01127-8
    49. 49
      Jiang, Z.; Fang, T. Dissociation Mechanism of H2O on Clean and Oxygen-Covered Cu(111) Surfaces: A Theoretical Study. Vaccum 2016, 128, 252258,  DOI: 10.1016/j.vacuum.2016.03.030
    50. 50
      Johnson, K.; Ström, V.; Wallenius, J.; Lopes, D. A. Oxidation of Accident Tolerant Fuel Candidates. J. Nucl. Sci. Technol. 2017, 54, 280286,  DOI: 10.1080/00223131.2016.1262297
    51. 51
      Johnson, K. D.; Wallenius, J.; Jolkkonen, M.; Claisse, A. Spark Plasma Sintering and Porosity Studies of Uranium Nitride. J. Nucl. Mater. 2016, 473, 1317,  DOI: 10.1016/j.jnucmat.2016.01.037
    52. 52
      Li, J.; Zhu, S.; Li, H.; Oguzie, E. E.; Li, Y.; Wang, F. Bonding Nature of Monomeric H2O on Pd: Orbital Cooperation and Competition. J. Phys. Chem. C 2009, 19311938,  DOI: 10.1021/jp809595y
    53. 53
      Akbay, T.; Staykov, A.; Druce, J.; Téllez, H.; Ishihara, T.; Kilner, J. A. The interaction of molecular oxygen on LaO terminated surfaces of La2NiO4. J. Mater. Chem. A 2016, 4, 13113,  DOI: 10.1039/c6ta02715f
    54. 54
      Fall, C. Ab Initio Study of the Work Functions of Elemental Metal Crystals; EPFL PP: Lausanne.
    55. 55
      Zhang, W.; Liu, L.; Wan, L.; Liu, L.; Cao, L.; Xu, F.; Zhao, J.; Wu, Z. Electronic Structures of Bare and Terephthalic Acid Adsorbed TiO2(110)-(1 × 2) Reconstructed Surfaces: Origin and Reactivity of the Band Gap States. Phys. Chem. Chem. Phys. 2015, 17, 2014420153,  DOI: 10.1039/C5CP01298H
  • Supporting Information

    Supporting Information


    The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcc.9b03076.

    • Mixed mode water adsorption configurations and adsorption energies on stoichiometric U3Si2{001}, {110}, and {111} surfaces; calculated adsorption energy (Eads), relevant bond lengths (d) of atomic (O) oxygen on {001}, {110}, and {111} U3Si2 surfaces (PDF)


    Terms & Conditions

    Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.