I. INTRODUCTION

Uranium nitrides are considered as promising fuel for the fast nuclear Generation IV reactors. Compared to many uranium and plutonium oxide nuclear fuels, uranium nitrides have several advantages [1, 2], such as higher thermal conductivity, melting point, metal density, and smaller lattice constant. Many research groups have identified an important role of UN in anti-corrosion application since the 1960s [3, 4, 5], and investigations on uranium nitride compounds renew great attention due to enviormental pollution and increasing interest for development of new nuclear fuels. The bulk properties of actinide nitrides have been investigated in these works, especially the elastic and magnetic properties.

Contrary to a number of available experimental data, some theoretical work have paid much attention on the pure and defective UN, physical properties of various defects (such as vacancy, O impurities, grain boundary) [6-13]. Basic bulk properties have been considered in these studies for uranium nitrides with emphasis on elastic and magnetic properties [14-18]. Petit et al. [9] have clarified the partial localization of 5f electrons in UN and reproduced the experimental total magnetic moment. In Ref. [7], the all-electron relativistic spin-polarized DFT calculations were performed to evaluate the total energies, optimized geometries, and electronic and thermodynamic properties of perfect stoichiometric UN and UN₂ single crystals. However, high reactivities of uranium nitrides with hydrogen, carbon, and oxygen at ambient atmosphere can affect the fuel fabrication process and fuel performance [1]. Experimental studies also clearly showed that oxygen in contact with the surface of uranium mononitride can result in growth of the oxide compound and, in initial stages, can lead to the formation of a surface layer structurally similar to oxynitrides UO_xN_y [19].

To predict UN fuel performance under different operating conditions and understand the material properties at the microscopic scale, it is crucial to investigate the surface properties and chemisorption process based on electronic structure calculations. However, a few papers have focused on the chemisorption behaviors of molecules or atoms on the surface of uranium nitrides.

Only recently [19-23], some groups have simulated the reactivity of molecules/atoms with the surface of uranium nitrides. These reports have indicated that O₂ molecule would spontaneously dissociate after chemisorption on the UN(001) surface, then the produced O atoms exhibit a strong chemisorption behavior.

According to Tasker’s analysis, the (001) surface must have the lowest surface energy for the rock-salt compounds [24]. At the same time, the (111) surface is usually the most dense plane for face-centered cubic (fcc) structure. Therefore, we consider two representative planes for the UN lattice–UN(001) and UN(111) surface. In this work, we perform electronic structure calculations within the framework of density functional theory (DFT) method to deeply understand chemisorption mechanisms of H, C, and O atoms on the UN(001) and (111) surfaces, shed light on the mechanisms of hydrogenization, oxidation, and cabonization of UN in the air at the electronic and atomic level, and analyse the chemical bonding of U 6d, 5f states with H 1s, C 2s, 2p, and O 2s, 2p states.

II. METHODOLOGY

Uranium mononitride (UN) has NaCl-type structure (face-centered cubic) with the experimental lattice constant a=0.4889 nm. In this paper, we use the generalized gradient approximation (GGA)+U (Hubbard U parameter to represent a correction to Coulomb repulsion interaction, separating f manifold into lower and upper Hubbard bands and removing f degrees of freedom from the Fermi level, U=4.0 eV) within the framework of density functional theory (DFT) and revised Perdew-Burke-Ernzerhof (RPBE) exchange-correlation functional with a periodic $\sqrt{2}\times \sqrt{2}$ three-layer slab model and single-sided chemisorption mode (an atom is placed on one side of the slab model, namely 0.5 ML adsorptivity) to simulate chemisorption behaviors of H, C, and O atoms on UN(001) and (111) surfaces in all calculations. A vacuum layer of 2.0 nm is added to a unit cell of the layers. The vacuum height test will be described in Sec. III in detail. We do not yet consider fully relativistic effects instead by using the scalar-relativistic approach. However, all other relativistic kinematic effects, such as mass-velocity, Darwin, and higher order terms are retained. Wang and Hay et al. [25, 26] have found that one can adequately describe the electronic and geometric properties of actinide complexes without treating spin-orbit effects explicitly, since we are interested in the chemisorption energies, defined as the difference in total energies. Therefore, we expect the inclusion of other relativistic effects, such as spin-orbit coupling, spin polarization, and orbital polarization effects, will not alter the main qualitative and quantitative conclusions of our work, as discussed in other works [27].

The outer fourteen electrons (6s²6p⁶5f⁶6d¹7s²) of the U atom are treated as the valence electrons and the remaining seventy-eight electrons are treated as the core electrons. DFT semi-core pseudopotentials (DSPP) and a double numerical basis set with polarization functions (DNP) have been used to treat the core electrons and the valence electrons, respectively. All electron basis sets are used for H, C, N, and O atoms. A 8×8×2 Monkhorst-Pack k-point mesh is applied for the Brillouin zone (BZ) integration. The convergence of a self-consistent field (SCF) is less than 1.0×10^-5 eV/atom. A plane-wave cutoff energy is fixed at E_cut=500 eV, which is enough for convergence of chemisorption calculations. Nonmagnetic configuration is appropriate for the light actinide U element from the point of total energy, so U 5f electrons are in the delocalized 5f³ electronic configuration in the present work.

Single atom, one per unit cell, is allowed to approach the UN(001) surface along four different symmetrical positions, namely (I) on the middle of the two nearest U atoms (Bridge position), (II) the adsorption atom sees a U atom located on the layer directly below the surface hollow site (Hollow position), (III) directly on top of a U atom (U-top position), and (IV) directly on top of a N atom (N-top position), as shown in Fig. 1, whereas only the former three symmetrical chemisorption positions exist for the UN(111) surface (not shown here). Chemisorption energy E_C is optimized with respect to the height R of the chemisorbed atom above the surface, and is given by [28]

Fig. 1.Full-screen View

UN(001) surface has four symmetrical chemisorption positions, namely Bridge, Hollow, U-top and N-top, whereas only the former three symmetrical chemisorption positions exist for UN(111) surface.

where E(M) is total energy of a bare UN(001) or UN(111) slab, E(X) is the total energy of the isolated atom in reference crystal structure, and E(M+X) is the total energy of the entire chemisorption system.

The relative change for bond length is utilized to describe the change for U-N bond length as a result of atom chemisorption

where Ri is the bond length for the central N atom on the UN(001) surface and the i^th U atom, R₀ is the original U-N bond length, and Δ denotes the relative change for Ri. The serial numbers i are shown in Fig. 2.

Fig. 2.Full-screen View

A configuration model for relaxation calculation. Ri (i=1, 2) and dij (i=1, 2, j=1, 2) represent U-N bond length in the intralayer and interlayer, respectively.

III. CONFIGURATIONS FOR UN(001) AND (111) SURFACES

To test the validity of the computational parameters, we first cleave the UN(001) surface, then check the convergence of total energy of the UN(001) surface with different vacuum thicknesses. We consider a vacuum thickness test to be convergent as long as change for total energy ΔE is less than 10 meV. Test results indicate that the total energy of a system is convergent when vacuum thickness is larger than 1.8 nm (Fig. 3). Therefore, we add a vacuum layer of 2.0 nm onto the unit cell of the periodic slab in order to reduce the influence of the boundary condition on the computation procedure. A model for the UN(001) surface is plotted in Fig. 4.

Fig. 3.Full-screen View

Convergence test of the vacuum thickness for the UN(001) surface.

Fig. 4.Full-screen View

A calculation model for the UN(001) surface.

The interactions of the surface atoms with UN matrix atoms will be unstable due to the absence of an adjacent atom, which is contrary to the matrix atoms. At the same time, these non-equilibrium interactions may cause the surface atoms to relax, reconstruct, find a new equilibrium site, and finally lower the total energy of the system. Moreover, this relaxation behavior may also change U-N bond length. A configuration model for relaxation calculation is shown in Fig. 2. In our work, relaxation is defined as relative change of U-N bond length, and the calculation result is listed in Table 1. From this table, we can see that the relative relaxation is relatively small (the maximum value at most 4.527%), so we fix the atoms in two low-lying layer in the following section, otherwise particularly declaration. For clarity, the configuration model and the relaxation calculation for the UN(111) surface are not shown here.

IV. RESULTS AND DISCUSSION

The energy minimum principle demonstrates that the higher the system symmetry, the lower the system energy, and the more stable the system. Therefore, we prefer an atom to be chemisorbed onto the high symmetrical position in the crystal surface. The UN crystal has fcc structure and several high symmetrical chemisorption positions exist on the crystal surface. As discussed above, we consider four representative positions in this work, namely Bridge, Hollow, U-top, and N-top. The chemisorbed atom is directly placed on the top of individual positions to study chemisorption behavior of the atom on these positions, as shown in Fig. 1. Chemisorption parameters not only include position, but also the orientation and height of the atom from the UN(001) and (111) surfaces, where chemisorption orientation is associated with the molecule structure and chemisorption height mainly depends on stability of the system.

A. H atom

We define chemisorption height h as the nearest distance for an H atom from UN(001) and UN(111) surfaces, and express it in terms of fractional coordinates. We optimize the H-UN(001)/UN(111) system with the minimum total energies (Table 2). The results for H chemisorption energies and Mulliken charges of individual configurations are listed in Table 3. As shown in Table 2, chemisorption heights of Bridge position are the same as Hollow position, and chemisorption energy differences are within several MeV. We think Bridge position to be equivalent with Hollow position. To verify this statement, we still consider Bridge, Hollow, and U-Top positions for C and O chemisorptions. The most stable position for H chemisorptions on the UN(001) and UN(111) surfaces are Hollow position (chemisorption energy 13.06 kJ/mol) and Bridge/Hollow position (chemisorption energy 11.532 kJ/mol), respectively.

TABLE 2.Full-screen View

Total energies of H-UN(001)/UN(111) systems with different chemisorption positions and heights. Numerical values for H-UN(111) system are listed in parentheses. Meanings for Bridge, Hollow, U-top and N-Top configurations are discussed in the text

Configuration	h (nm)	Esystem(eV)
Bridge	0.263 (0.282)	-11569.569 (-3094.786)
	0.264 (0.283)	-11569.569 (-3094.776)
	0.265 (0.284)	-11569.576 (-3094.786)
	0.266 (0.285)	-11569.575 (-3094.775)
	0.267 (0.286)	-11569.571 (-3094.773)
Hollow	0.251 (0.282)	-11570.033 (-3094.776)
	0.252 (0.283)	-11570.037 (-3094.776)
	0.253 (0.284)	-11570.038 (-3094.786)
	0.254 (0.285)	-11570.037 (-3094.775)
	0.255 (0.286)	-11570.033 (-3094.773)
U-top	0.274 (0.324)	-11569.733 (-3093.689)
	0.275 (0.325)	-11569.742 (-3093.691)
	0.276 (0.326)	-11569.746 (-3093.691)
	0.277 (0.327)	-11569.744 (-3093.688)
	0.278 (0.328)	-11569.737 (-3093.687)
N-top	0.237	-11569.194
	0.238	-11569.229
	0.239	-11569.245
	0.240	-11569.244
	0.241	-11569.231

Show more

TABLE 3.Full-screen View

Chemisorption configurations, chemisorption energies E_C (in kJ/mol) and Mulliken charges Q (in e) for H atom. Numerical values for H-UN(111) system are listed in parentheses. Meanings for Bridge, Hollow, U-top and N-Top configurations are discussed in the text

Configuration	EC (kJ/mol)	Q(e)
Bridge	12.023 (11.532)	-0.229 (-0.341)
Hollow	13.058 (11.531)	-0.368 (-0.341)
U-top	11.893 (11.242)	-0.144 (-0.314)
N-top	10.602	0.206

Show more

To further understand interaction of H atoms with the UN(001) surface, we analyze the projected density of states (PDOS) before and after H chemisorption in terms of electronic structure calculations, as shown in Fig. 5(a)–5(c). In fact, the ground state valence electronic configurations of U and N atoms are 5f³6d¹7s² and 2s²2p³, respectively. From the point of ionic bonding behavior, the U6d and 7s electrons fill the N2p states and the three U5f electrons form the highest occupied molecular orbital (HOMO) [29]. The characteristic peak of H 1s PDOS in the energy range of 17.157 eV and 20.271 eV shifts towards lower energy band induced by H chemisorption (Fig. 5(a)) and expands its peak area, indicating that H atoms gain electrons as a result of chemisorption, which is in agreement with the Mulliken charge analysis (the third row in Table 3). We neglect other s, p orbitals in U atoms because 6d, 5f orbitals dominate the structural and electronic properties of the U element.

Fig. 5.Full-screen View

Projected density of states (PDOS) of H 1s orbital (a), U 6d, 5f orbitals (b) and H 1s, U 6d, 5f orbitals (c) before and after H chemisorption on the Hollow position of UN(001) surface. The Fermi energy (dash line) stands at 0 eV.

We can see that H chemisorption causes the shape of U 5f PDOS with energy of -1.026∼0 eV to slightly change and peak position of U 5f PDOS shifts from 0.353 eV to 0.502 eV (Fig. 5(b)). From this plot we can see that peak value decreases from 247.337 states/eV to 211.083 states/eV and peak area of U 6d PDOS obviously diminishes with energies of -1.042 eV and 1.515 eV, which show that U 6d, 5f orbitals both lose electrons, also in agreement with the Mulliken charge analysis (Table 3). Fig. 5(c)) depicts H 1s and U 6d, 5f PDOS. H 1s orbitals remarkably overlap with U 6d orbitals in the energy of -5.043 eV and 1.161 eV. However, H 1s orbitals almost separate from U 6f orbitals with energy intervals of -4.152 eV and -1.241 eV. It is known that the larger the overlapping area of PDOS, the higher the hybridized bonding. Therefore, the H atom will hybridize with the U atom as a consequence of H chemisorption. Electronic charges of U 6d, 5f orbitals (mainly U 5f orbitals) transfer to the H 1s orbital, which is also consistent with the Mulliken charge analysis (Table 3).

Similarly, PDOSs before and after H chemisorption in Bridge or Hollow positions for the UN(111) surface are shown in Figs. 6(a)–6(c). H chemisorption induces an H 1s peak around the Fermi level to the energy band in the energy range of -7.374 eV to -2.502 eV (Fig. 6(a)). However, H chemisorption has little effect on U 6d, 5f PDOSs. U 6d states clearly hybridize with U 5f states between -2.215 eV and 2.00 eV, as shown in Fig. 6(b). H 1s states clearly overlap with U 6d, 5f states (Fig. 6(c)), and acquire electrons from U atom, which is consistent with Mulliken analysis (Table 3).

Fig. 6.Full-screen View

Projected density of states (PDOS) of H 1s orbital (a), U 6d, 5f orbitals (b) and H 1s, U 6d, 5f orbitals (c) before and after H chemisorption on the Bridge or Hollow positions of UN(111) surface. The Fermi energy (dash line) stands at 0 eV.

B. C atom

Previous report [2] has shown that U carbides can effectively retard corrosion of U metal. The crystal structure of UN is the same as that of UC (NaCl type, face-centered cubic structure) and the lattice constant is also very close to the latter (lattice constants of UN and UC are 0.4965 nm and 0.4889 nm, respectively). Therefore, we wonder whether chemisorption behavior of a C atom on UN(001) and UN(111) surfaces is also similar with that of a U monocarbide.

Next we perform electronic structure calculations of total energy for the C-UN(001)/UN(111) systems with different chemisorption positions and heights using the method discussed above and relax different configurations to find the optimal chemisorption positions and heights. The total energy results of the C-UN(001)/UN(111) systems are listed in Table 4. According to the chemisorption energy formula (Eq. (1)), a configuration with the maximum chemisorption energy is the most stable configuration for chemisorption on the UN(001) and UN(111) surface. Therefore, we optimize configurations with the minimum total energy for different chemisorption positions and the optimal configuration is the most stable configuration for certain chemisorption positions.

TABLE 4.Full-screen View

Total energies of C-UN(001)/UN(111) systems with different chemisorption positions and heights. Numerical values for C-UN(111) system are listed in parentheses. Meanings for Bridge, Hollow, U-top and N-Top configurations are discussed in the text

Configuration	h (nm)	Esystem(eV)
Bridge	0.253 (0.265)	-12587.704 (-3232.383)
	0.254 (0.266)	-12587.720 (-3232.386)
	0.255 (0.267)	-12587.727 (-3232.387)
	0.256 (0.268)	-12587.726 (-3232.384)
	0.257 (0.269)	-12587.719 (-3232.379)
Hollow	0.249 (0.265)	-12587.983 (-3232.383)
	0.250 (0.266)	-12587.995 (-3232.386)
	0.251 (0.267)	-12588.001 (-3232.387)
	0.252 (0.268)	-12588.000 (-3232.384)
	0.253 (0.269)	-12587.993 (-3232.379)
U-top	0.276 (0.314)	-12586.574 (-3229.448)
	0.277 (0.315)	-12586.578 (-3229.456)
	0.278 (0.317)	-12586.579 (-3229.459)
	0.279 (0.318)	-12586.574 (-3229.456)
	0.280 (0.319)	-12586.564 (-3229.447)
N-top	0.246	-12589.483
	0.247	-12589.528
	0.248	-12589.541
	0.249	-12589.527
	0.250	-12589.490

Show more

The chemisorption energy and Mulliken charge analysis results are presented in Table 5. From this table we can see that chemisorption energies reach the maximum value when the C atom chemisorbs or at the Bridge and N-top positions of the UN(001) surface (chemisorption energies are 25.513 kJ/mol and 25.502 kJ/mol, respectively). The chemisorption energy of the Bridge configuration is very close to that of the N-top configuration. The relaxtion calcultion shows that the former is unstable, and will transform to the latter. Therefore, the most stable configuration for C chemisorption on the UN(001) surface is N-top configuration. Similarly, relaxation result shows that the most stable configuration for C chemisorption on the UN(111) surface is the Bridge or Hollow configurations.

TABLE 5.Full-screen View

Chemisorption configurations, chemisorption energies E_C (in kJ/mol) and Mulliken charges Q (in e) for C atom. Numerical values for C-UN(111) system are listed in parentheses. Meanings for Bridge, Hollow, U-top and N-Top configurations are discussed in the text

Configuration	EC (kJ/mol)	Q(e)
Bridge	25.513 (24.132)	-0.400 (-0.543)
Hollow	18.81 (24.132)	-0.142 (-0.542)
U-top	15.69 (14.147)	-0.103 (-0.011)
N-top	25.502	-0.392

Show more

To further understand the interactions of C with the UN(001) and (111) surfaces, we analyze PDOS before and after C chemisorption in terms of electronic structure calculations, as shown in Figs. 7(a)–7(d) and Figs. 8(a)–8(d). C 2s and 2p PDOSs shift towards a lower energy band after chemisorption (Fig. 7(a)), and the peak area of C 2s and 2p PDOSs obviously increase, indicating that C atoms gain electrons from the adsorbate (UN(001) surface in the present work) as a result of chemisorption, which is in agreement with the Mulliken charge analysis (Table 5). U 6d, 5f PDOSs before and after chemisorption are plotted in Fig. 7(b). C chemisorption has a negligible effect on U 6d PDOS, however, the peak value of the U 6d PDOS decreases (Fig. 7(b)), which shows that U 6d orbitals lose electrons. At the same time, C chemisorption also has a negligible effect on the U 5f PDOS in the energy interval of 0∼1.520 eV. However, a new characteristic peak appears in the energy range of -1.0∼0.50 eV, showing that U 5f orbitals hybridize with C 2s or 2p orbitals. As shown in Fig. 7(b), peak position of U 5f PDOS shifts from 0.508 eV to 0.572 eV and peak value decreases from 241.746 states/eV to 227.908 states/eV, indicating that U 5f orbitals also lose electrons. Figs. 7(c)–7(d)) depict interactions of U 6d, 5f orbitals with C 2s, 2p orbitals. C 2s orbitals and U 6d orbitals produce a remarkable hybridization/mixing peak (Fig. 7(c)), while the C 2s PDOS almost separates from the U 5f PDOS. However, C 2p orbitals obviously overlap with U 6d, 5f orbitals (Fig. 7(d)), showing that C 2p orbitals hybrid with U 6d, 5f orbitals to form a covalent bond, which is in sharp contrast with Fig. 7(c). Therefore, electronic charges of U 6d orbitals transfer to C 2s, 2p orbitals as a result of C chemisorption on the UN(001) surface and U 5f orbitals transfer to C 2p orbitals, these results are consistent with the Mulliken charge analysis (Table 5).

Fig. 7.Full-screen View

Projected density of states (PDOS) of C 2s, 2p orbital (a), U 6d, 5f orbital (b), C 2s, U 6d, 5f orbital (c) and C 2p, U 6d, 5f orbital (d) before and after C chemisorption on N-top position of the UN(001) surface. The Fermi energy (dash line) stands at 0 eV.

Fig. 8.Full-screen View

Projected density of states (PDOS) of C 2s, 2p orbital (a), U 6d, 5f orbital (b), C 2s, U 6d, 5f orbital (c) and C 2p, U 6d, 5f orbital (d) before and after C chemisorption on Bridge or Hollow positions of the UN(111) surface. The Fermi energy (dash line) stands at 0 eV.

As shown in Fig. 8(a), C chemisorption widens the C 2p peak around the Fermi level, and new C 2s states appear in the energy interval of -5.241∼-2.632 eV to -21.352∼-18.471 eV (Fig. 8(a)). However, C chemisorption also has a negligible effect on U 6d, 5f PDOSs (Fig. 8(b)). C 2s states slightly overlap with U 6d, 5f states (Fig. 8(c)). However, C 2p states obviously hybridize with U 6d, 5f states (Fig. 8(d)) and acquire a majority of electrons (about 2/3) from the U atom, which is consistent with Mulliken analysis (Table 5).

V. CONCLUSION

We also performed electronic structure calculations for O chemisorption on the UN(001) and UN(111) surfaces, and depicted the PDOS before and after chemisorption in Figs. 9(a)–9(d) and Figs. 10(a)–10(d). C 2s and 2p PDOSs shift towards a lower energy band after chemisorption (Fig. 9(a)) and the peak area of the C 2p PDOS increases, indicating that C atoms gain electrons from the UN(001) as a result of chemisorption, which is in agreement with the Mulliken charge analysis (the third row in Table 7). Peak areas of the U 6d, 5f PDOS in the energy of -1.5∼1.3 eV diminish, which shows that U 6d, 5f electrons transfer to the O 2s or 2p orbitals. However, peak areas of the U 6d, 5f PDOS in the energy of -5.0∼-1.50 eV increase, indicating that U 6d, 5f orbitals may hybridize with O 2s or 2p orbitals (Fig. 9(b)). The peak position of U 5f PDOS shifts from 0.354 eV to 0.418 eV, and the peak value decreases from 247.361 states/eV to 237.268 states/eV, implying that U 5f orbitals lose electrons. U 6d orbitals and O 2s orbitals form a relatively small overlapping peak, while the U 5f PDOS is almost separate from the O 2s PDOS (Fig. 9(c)). However, O 2p PDOS obviously overlap with U 6d, 5f PDOSs (Fig. 9(d)), particularly the U 6d PDOS, which is in sharp contrast with Fig. 9(c). It is shown that the smaller the overlapping area, the weaker the chemical bonding and the lower the transferred charge. Therefore, electronic charges of U 6d orbitals transfer to O 2s and 2p orbitals (mainly O 2p orbitals) and U 5f orbitals transfer to O 2p orbitals, which is in agreement with the Mulliken charge analysis (Table 7).

Fig. 9.Full-screen View

Projected density of states (PDOS) of H 1s, O 2s, 2p orbital (a), U 6d, 5f orbital (b), O 2s, U 6d, 5f orbital (c) and O 2p, U 6d, 5f orbital (d) before and after O₂ chemisorption on the BU position of the UN (001) surface. The Fermi energy (dash line) stands at 0 eV.

Fig. 10.Full-screen View

Projected density of states (PDOS) of H 1s, O 2s, 2p orbital (a), U 6d, 5f orbital (b), O 2s, U 6d, 5f orbital (c) and O 2p, U 6d, 5f orbital (d) before and after H₂O chemisorption on Bridge or Hollow positions of the UN(111) surface. The Fermi energy (dash line) stands at 0 eV.

From Fig. 10(a), we can see that O chemisorption shifts the O 2s, 2p peak to a lower energy, and widens O 2p PDOS. However, O chemisorption has an ignorable effect on U 6d, 5f PDOSs (Fig. 10(b)). O 2s states slightly overlap with U 6d, 5f states (Fig. 10(c)). However, O 2p states remarkably hybridize with U 6d, 5f states (Fig. 10(d)) and acquire electrons from the U atom, which is consistent with Mulliken analysis (Table 7).

In the present work, we performed density functional theory calculations for H, C, and O chemisorptions on UN(001) and UN(111) surfaces using the generalized gradient approximation (GGA) and revised Perdew-Burke-Ernzerhof (RPBE) exchange-correlation functional at a non-spin polarized level with a periodic slab model. The results show that Hollow, N-top, and Hollow chemisorption sites are the most stable sites for H, C, and O atoms chemisorbed on the UN(001) suraface, respectively. The calculated electronic density of states (DOS) demonstrate electronic charge transfer between s, p orbitals in chemisorbed atoms and U 6d, 5f orbitals and transferred electronic charges agree with the Mulliken charge analysis. Electronic structure calculations of H, C and O atom chemisorption on Bridge, Hollow, and U-Top positions of the UN(111) surface indicate that Bridge position is equivalent with Hollow position, two symmetrical chemisorption positions exist on UN(111) surface, namely Bridge (or Hollow) and U-Top, and the most stable chemisorption position for H, C, and O atoms is Bridge (or Hollow) position.

In the future, we plan to investigate the chemisorption of other atoms or molecules on the surface of actinide compounds, especially promising nuclear fuels such as Pu and U compounds, providing corresponding anti-corrosion techniques, and improving operational life and efficiency for nuclear fuel.