Stability and optical properties of plutonium monoxide from first-principle calculation.

The resolution of questions about the existence of condensed plutonium monoxide (PuO) has long been hindered by lack of thermochemical data. Here we perform first-principles calculation to investigate the reaction Pu2O3 + Pu → 3 PuO and find that PuO is thermodynamically unstable under ambient pressure. We also find that pressure could stabilize PuO by strengthening the hybridization between Pu-5f/6d and O-2p states. Moreover, the dynamical stability of NaCl-type PuO is verified by the phonon calculation. Optical properties such as reflectivity are also predicted for the detection of metallic PuO.

Introduction

Despite the technological importance and scientific interest of metallic plutonium (Pu), its phase diagram with oxygen is still controversial. Now it is recognized that the metal is covered by a layer of trivalent sesquioxide (Pu2O3), followed by a layer of dioxide (PuO2)[1,2]. The existence of Pu2O3 is inevitable because the reduction of PuO2 by Pu (i.e., 3 PuO2 + Pu → 2 Pu2O3) is thermodynamically favourable[1]. In addition, PuO2− is shown to exist between PuO2 and Pu2O3 [3]. However, the existence of plutonium monoxide (PuO) between Pu2O3 and Pu or in the bulk form has not been proven conclusively. Although several investigators reported its preparation including the surface[4-6] and bulk phase[3,7-10], there are also strong contradictory evidences[11-13].

Sooner after the preparation of metallic plutonium in 1944, PuO was first observed by Zachariasen[4] as a silver layer on the surface of plutonium metal, and it was identified using X-ray diffraction (XRD) as the NaCl-type structure with the lattice parameter being 4.948 ± 0.002 Å[4,14]. Another value of lattice parameter reported by Coffinberry and Ellinger was given as 4.96 ± 0.01 Å[3,5]. In 1960, Akimoto[7] reported the preparation of PuO in the bulk form by the reaction of molten plutonium with the oxygen generated from Ag2O. The measured lattice parameter was given by 4.960 ± 0.003 Å. Later Chikalla et al.[3,8]. prepared PuO by reacting PuO2 and carbon at 1800°C in helium and at 1550°C in vacuum. However, Mulford et al.[11]. stated that all attempts to prepare PuO in bulk at Los Alamos have been unsuccessful. In addition, the authors found that the solid solution phase, Pu(C, O), could not exist with a high oxygen content and thus argued that PuO does not exist as a pure binary phase. In fact, Chikalla et al.[3,8]. already suspected that their compound may be a carbon-stabilized form or an oxycarbide. Furthermore, the high-temperature XRD study of Terada et al.[6], which was performed on the PuO2-coated metal, showed that the formation reaction of PuO,begins only at high temperature of about 250°C. The authors also noted that the free energy change for reaction (1) is positive as calculated from Oetting’s table[12], i.e., PuO is metastable. The landmark experiment was performed by Larson and Haschke[13], who used X-ray photoelectron spectroscopy (XPS) and Auger electron spectroscopy (AES) to identify the surface phase of vacuum-heated-treated PuO2-coated metal. They showed that the surface phase previously identified as PuO is actually PuOC.

The story of PuO was supposed to end here. However in 1990, Haschke[9] reopen this question by claiming that the residue from the decomposition of plutonium monoxide monohydride (PuOH)[15,16] may be the metastable PuO. Because of the extreme pyrophoric nature of the residue, even XRD measurements were not attempted. Afterwards Russian scientist Reshetnikov[10] reported the preparation of pure PuO in gramme amounts by reacting plutonium oxychloride with calcium. But there is no any characterization of this compound.

Since this question has existed for such a long time, it is necessary to find the solution from a new point of view. In recent years, first-principles calculations have emerged as a powerful and indispensable tool for quantitatively predicting the properties of materials such as actinide oxides[17]. In addition, plutonium and its compounds are sensitive to variables such as impurity content, temperature, and pressure[1]. In particular, pressure is often used as an efficient tool to stabilize the materials or explore the new phase. Larson and Haschke[13] already argued that PuO may be synthesized under high pressure, by noting that several lanthanide elements (La-Nd[18], Sm[18,19], Yb[20]) react with their higher oxides to form NaCl-type monoxides under pressure. For actinide elements, ThO was proved to be stable under pressure through first-principles calculation[21]. Therefore, in this work we employ the density-functional theory (DFT) and DFT + U to study the reaction (1) under pressure as well as the physical properties of PuO.

From available literature, α-Pu is stable up to around 40 GPa[22], whereas no report is available on the high-pressure study of α-Pu2O3. In addition, for our estimate of formation enthalpy of Pu, we used the δ-phase, which is only 0.047 eV/atom[12] higher than the (most stable) α-phase. The first-principles calculation of α-Pu is not a trivial task since there is a significant site-selective electronic correlation in this complicated metal[23]. Similar treatment was used in the previous literature such as the calculation of formation energies of plutonium oxides[24].

Methods

The electronic structure problems related to first-principles calculation were solved using the code Vienna Ab-initio Simulation Package (VASP)[25]. Within the framework of density-functional theory (DFT)[26,27], the exchange correlation potential in the Kohn-Sham equation[27] was described by the local-density approximation (LDA)[28] and the generalized gradient approximation (GGA) of Perdew, Becke and Ernzerhof (PBE)[29]. Moreover, the rotationally invariant DFT + U approach[30] was also utilized to capture the possible localization effect of 5f electrons coming from the strong electron-electron interaction. Within DFT + U, the “full localized limit” double counting expression and the typical parameters: U = 4.0 eV and J = 0.7 eV, are chosen here. These choices were already proved to give a reasonable description of Pu[24] and Pu2O3 [31,32]. In addition, scalar relativistic calculations were performed using spin-orbit coupling (SOC) but SOC was not included in the enthalpy calculation because it is already found that the inclusion of SOC has insignificant effect on the bulk properties[24,31-36]. For example, within GGA, inclusion of SOC increases the equilibrium lattice parameter of PuO by only 0.14%. However, GGA + SOC is used for electronic structure calculation to illustrate the spin-orbit splitting in the 5f states.

Furthermore, the projector augmented-wave (PAW) method[37,38] was chosen to solve the Kohn-Sham equation. The PAW pseudopotentials for Pu and O were generated using the 6s 27s 26p 66d 25f  4 and 2s 22p 4 valence electronic configurations, respectively. These pseudopotentials are implemented in the VASP code[25]. All computational parameters were carefully selected to yield a total energy convergence better than 1 meV/atom. Here the cut-off energy for plane wave basis set was set as 700 eV and the Brillouin zone of δ-Pu, PuO, and α-Pu2O3 are sampled by 8 × 8 × 8, 8 × 8 × 8, and 4 × 4 × 4 Monkhorst-Pack k point sampling mesh[39], respectively. Here we also considered three possibilities for the magnetic states: nonmagnetic (NM), antiferromagnetic (AFM), and ferromagnetic (FM). Within the above computational schemes, the lowest-energy magnetic states of δ-Pu and PuO are AFM, while that of α-Pu2O3 is FM for pure DFT and AFM for DFT + U. These findings are in accord with the previous literature[24,31-36,40].

Results and Discussion

Thermodynamical stability under ambient pressure

The reaction enthalpy of (1) was defined bywhere H PuO, H Pu, and are the formation enthalpy of NaCl-type PuO, δ-Pu, and α-Pu2O3, respectively. At ambient pressure, the lattice parameter of PuO, a PuO, the formation enthalpies of PuO and Pu2O3, and the reaction enthalpy calculated using LDA/GGA and LDA/GGA + U are listed in Table 1. For comparison, the available experimental data and the previous theoretical results using self-interaction correction (SIC)[41] and full-potential linearized augmented plane wave (FLAPW)[34] are also shown.

Table 1

The lattice parameter a PuO and formation enthalpy H PuO of PuO, formation enthalpy of Pu2O3, , and the reaction enthalpy defined by equation (2) under ambient pressure. For comparison, we also show the FLAPW results in the square brackets, the SIC result and the experimental values.

Method	a PuO (Å)	HPuOf (eV/f.u.)	HPu_2O_3f (eV/f.u.)	ΔH(eV)
LDA	4.79 [4.80[34]]	−5.37	−15.61	−0.15
GGA	4.92 [4.96[34]]	−4.82	−14.54	0.50
LDA + U	4.89 [4.99[34]]	−5.76	−17.13	0.41
GGA + U	5.03 [5.12[34]]	−5.01	−15.77	1.49
LDA + SIC	5.13[41]
Exp.	4.96 ± 0.01[5]	−5.85 ± 0.22[12]	−18.65 ± 0.17[12]	1.10 ± 0.83
			−16.39[1]

For a PuO, the calculated values are consistent with the previous theoretical results[34], and GGA and LDA + U reproduce well the experimental data. The deviation between our work and the FLAPW results[34] may due to the different magnetic states and different schemes of DFT + U. The lattice parameters predicted by DFT + U are larger than those from pure DFT. This is because DFT + U scheme favours the localization of Pu 5f electrons, consequently Pu 5f electrons participate less to the bonding which leads to an increase of lattice parameter. In addition, since the overbinding is usually found in the LDA, the LDA underestimate the lattice parameter. Overall, GGA and LDA + U gives a better description while pure LDA underestimates the lattice parameter and GGA + U overestimates it. Similar deviations of LDA/GGA and LDA/GGA + U have been found in their description of Pu[24,40] and PuO2 [24,42,43].

The calculated formation enthalpies of PuO and Pu2O3 are in reasonable agreement with the experimental values. This indicates that our calculations are consistent and reliable. The larger formation enthalpies found in LDA(+U) with respect to GGA are consistent with the overbinding in LDA(+U). For same reason, the reaction (1) is more exothermic in LDA(+U) than in GGA(+U). Under ambient pressure, the calculated reaction enthalpy is positive only except the LDA result. However, the lattice parameter of PuO predicted by LDA deviates the experiment by 0.17 Å, which corresponds to a larger negative pressure ~ −13 GPa. Similar underestimate by LDA is also found in the other plutonium compounds[24,40,42,43]. We will see that is positive at a small negative pressure of ~ −1 GPa (see Fig. 1(a)). Thus it is reasonable to reckon the modified value of in LDA by enlarging the cell as positive. Overall, our first-principles calculation predicted that PuO is thermodynamically unstable under ambient pressure. As discussed in the introduction, this finding agrees with the high temperature XRD[6] and XPS-AES[13] experiment, whereas disagrees with the recent preparation experiment[10].

Figure 1

Thermodynamical and dynamical stability of PuO. (a) The reaction enthalpy and (b) the phonon density of states of PuO at 15 GPa obtained from LDA/GGA and LDA/GGA + U.

Thermodynamical and dynamical stability under high pressure

In Fig. 1(a) we plot the reaction enthalpy as a function of pressure which ranges from −5 GPa to 15 GPa. It is clear that the reaction enthalpy becomes negative when pressure is increased above a critical value P c. That is to say, PuO is thermodynamically stable above this critical pressure P c. From Fig. 1(a), the critical pressures predicted by GGA and GGA + U are 6 and 12 GPa, respectively. Thus it is likely that P c lies within this range 6–12 GPa.

In the following, we will take a value of 15 GPa as the pressure after the transition. The dynamical stability of PuO under pressure is examined by the phonon calculations. Figure 1(b) plot the phonon density of states (DOS) of PuO using LDA/GGA and LDA/GGA + U. No imaginary frequency was found for the entire Brillouin zone. The large frequency gap observed in Fig. 1(b) was attributed to the large mass ratio between Pu and O, i.e., Pu atoms contribute the lowest acoustic modes (0–5 THz) and O atoms contribute the higher modes (10–13 THz). In a word, PuO is dynamically stable in rocksalt structure.

Electronic structure

To obtain insights into the stabilization of PuO under pressure, we investigated and compared the electronic structures at 0 GPa and 15 GPa. The total density of states (DOS) and projected DOS (pDOS) for Pu-5f, Pu-6d and O-2p states are plotted in Fig. 2. Note that SOC is included in the following. Apparently, PuO is metallic, which is consistent with the previous theoretical result[34].

Figure 2

Total and projected density of states within the GGA + SOC formalism under pressure of (a) 0 GPa and (b) 15 GPa. The Fermi level is denoted by the dashed line.

It was well-known that Pu-6d orbital is quasi-degenerate with Pu-5f orbital and thus the hybridization between Pu-5f and Pu-6d states are very strong in the plutonium compounds. This feature is apparent in Fig. 2 and those bands could be denoted by Pu-5f/6d states. We can also see that the conduction bands and the semicore bands (energy bands between −8 eV and −4 eV) are mainly contributed by Pu-5f/6d and O-2p states, respectively. The sharp-peak feature of Pu-5f pDOS indicates the localization of the Pu-5f states. In the conduction and semicore bands, the shape of Pu-5f and Pu-6d pDOS follow that of O-2p, which implies that there is a considerable hybridization between O-2p and Pu-5f/6d states. This also indicates the presence of Pu-O bond. Through Bader charge analysis[44], the charge transfer from Pu to O is evaluated to be about 1.3 e, which approximates half of that in PuO2 [43]. This suggests the formal oxidation state of Pu in PuO is indeed II.

Next let us focus the effect of pressure. Under pressure, the bandwidth of conduction bands and semicore bands is increased and the band gap among them is decreased. This implies that the Pu-5f states become relatively itinerant and participate more and more to the hybridization and bonding. To clarify this, let us make a close look at the band structures. The band structures of PuO along -X-L--W directions are shown in Fig. 3. The band shift under pressure is highlighted and it is clear that pressure forces Pu-5f/6d bands and O-2p to mix together at point. So we conclude that hybridization among Pu-5f, Pu-5d, and O-2p serves as the dominant stabilizing factor.

Figure 3

Band structures around the Fermi level within the GGA + SOC formalism under pressure of (a) 0 GPa and (b) 15 GPa. The Fermi level is denoted by the dashed line. The red lines show the crossing band under pressure.

Optical properties

For the synthesis of LnO (Ln = Lanthanide) under high pressure[18], the mixture of Ln and Ln2O3 were introduced into a compressible gasket apparatus of the belt type calibrated by using the polymorphic transition of bismuth and barium. After being unloaded, the products were identified from XRD and chemical analysis. For the synthesis of PuO, modern in situ diamond-anvil cell x-ray/laser-heating experiment may provide more information than the older methods. For example, XPS could also be used to identify the phase of PuO. In particular, optical detection is often used in the characterization of new metallic phase in the high pressure experiment. For example, the important evidence of the metallic hydrogen at 495 GPa is the high reflectivity[45].

Here we also calculate the optical properties of PuO including the reflectivity and absorption coefficient, as shown in Fig. 4. The optical conductivity is computed from the Kubo-Greenwood formula[46,47] within the independent-particle approximation. And then the optical properties are calculated using Kramers-Kronig relation. The technical details are described in ref.[48]. The peak observed at 2 eV is mainly contributed by inter-band transition, i.e., the transition from Pu-5f states at the bottom of the conduction band to the Pu-5f states at the top of the conduction band. And the peaks at 6 and 10 eV are contributed from the O-2p states to the conduction band. The peaks at higher energy are due to transition from semicore states to the conduction band. This information is useful for the optical detection of the plutonium monoxide in the high pressure experiment.

Figure 4

The calculated reflectivity (left) and absorption coefficient (right) of PuO at 15 GPa using the GGA + SOC scheme.

Conclusion

In this work, we have calculated the reaction enthalpy of reaction Pu2O3 + Pu → 3 PuO in the framework of DFT and DFT + U. The results suggest that PuO is thermodynamically unstable at ambient pressure and probably decompose into α-Pu2O3 and Pu. However, PuO may be stabilized by the other environment variables such as impurity and temperature. Here we find that PuO is likely to exist at the pressure above 12 GPa by calculating the reaction enthalpies and phonon density of states under high pressure. The origin of the stabilization under pressure lies in the delocalization of Pu-5f states and the increasing hybridization between Pu-5f/6d and O-2p states. Furthermore, the optical properties of PuO are also predicted and may provide a useful reference for future experiment.