Hua Y. Geng

Structural behavior of uranium dioxide under pressure by LSDA + U calculations

The structural behavior of UO2 under high pressure up to 300GPa has been studied by first-principles calculations with LSDA+U approximation. The results show that a pressure-induced structural transition to the cotunnite-type (orthorhombic Pnma) phase occurs at 38GPa. It agrees well with the experimentally observed ~42 GPa. An isostructural transition following that is also predicted to take place from 80 to 130GPa, which has not yet been observed in experiments. Further high compression beyond 226GPa will result in a metallic and paramagnetic transition. It corresponds to a volume of 90A^3 per cell, in good agreement with a previous theoretical analysis in the reduction of volume required to delocalize 5f states.

Interplay of defect cluster and the stability of xenon in uranium dioxide from density functional calculations

Self-defect clusters in bulk matrix might affect the thermodynamic behavior of fission gases in nuclear fuel such as uranium dioxide. With first-principles LSDA+U calculations and taking xenon as a prototype, we find that the influence of oxygen defect clusters on the thermodynamics of gas atoms is prominent, which increases the solution energy of xenon by a magnitude of 0.5 eV, about 43% of the energy difference between the two lowest lying states at 700 K. Calculation also reveals a thermodynamic competition between the uranium vacancy and tri-vacancy sites to incorporate xenon in hyper-stoichiometric regime at high temperatures. The results show that in hypo-stoichiometric regime neutral tri-vacancy sites are the most favored position for diluted xenon gas, whereas in hyper-stoichiometric condition they prefer to uranium vacancies even after taking oxygen self-defect clusters into account at low temperatures, which not only confirms previous studies but also extends the conclusion to more realistic fuel operating conditions. The observation that gas atoms are ionized to a charge state of Xe+ when at a uranium vacancy site due to strong Madelung potential implies that one can control temperature to tune the preferred site of gas atoms and then the bubble growth rate. A solution to the notorious meta-stable states difficulty that frequently encountered in DFT+U applications, namely, the quasi-annealing procedure, is also discussed.

Stability mechanism of cuboctahedral clusters inUO2+x: First-principles calculations

The stability mechanism of cuboctahedral clusters in nonstoichiometric uranium dioxide is investigated by first-principles LSDA+U method. Calculations reveal that the structural stability is inherited from U6O12 molecular cluster whereas the energy gain through occupying its center with an additional oxygen makes the cluster win out by competition with point oxygen interstitials. Local displacement of the center oxygen along direction also leads the cluster 8-folded degeneracy and increases relatively the concentration at finite temperatures. But totally, elevation of temperature, i.e., the effect of entropy, favors point interstitial over cuboctahedral clusters.

Ab initio investigation on oxygen defect clusters in UO2+x

Oxygen defect clustering in uranium dioxide had been indicated in powder neutron diffraction measurements, and an empirical clustering mechanism had been proposed to explain the data. However, using first-principles LSDA+U calculations, we find that this empirical model, in fact, cannot work. A more physically reasonable model is proposed based on a thermodynamical competition between point defects and cuboctahedral clusters. This mechanism interprets the puzzled origin of the observed asymmetric interstitial O′ and O″ naturally. It also gives a good and consistent agreement with all available experimental data, except the high occupation of the O″ site.

First-principles study on oxidation effects in uranium oxides and high-pressure high-temperature behavior of point defects in uranium dioxide

Formation Gibbs free energy of point defects and oxygen clusters in uranium dioxide at high-pressure high-temperature conditions are calculated from first principles, using the LSDA+U approach for the electronic structure and the Debye model for the lattice vibrations. The phonon contribution on Frenkel pairs is found to be notable, whereas it is negligible for the Schottky defect. Hydrostatic compression changes the formation energies drastically, making defect concentrations depend more sensitively on pressure. Calculations show that, if no oxygen clusters are considered, uranium vacancy becomes predominant in overstoichiometric UO2 with the aid of the contribution from lattice vibrations, while compression favors oxygen defects and suppresses uranium vacancy greatly. At ambient pressure, however, the experimental observation of predominant oxygen defects in this regime can be reproduced only in a form of cuboctahedral clusters, underlining the importance of defect clustering in UO2+x. Making use of the point defect model, an equation of state for non-stoichiometric oxides is established, which is then applied to describe the shock Hugoniot of UO2+x. Furthermore, the oxidization and compression behavior of uranium monoxide, triuranium octoxide, uranium trioxide, and a series of defective UO2 at zero Kelvin are investigated. The evolution of mechanical properties and electronic structures with an increase of the oxidation degree are analyzed, revealing the transition of the groundstate of uranium oxides from metallic to Mott insulator and then to charge-transfer insulator due to the interplay of strongly correlated effects of 5f orbitals and the shift of electrons from uranium to oxygen atoms.

High-pressure behavior of dense hydrogen up to 3.5 TPa from density functional theory calculations

Structural behavior and equation of state (EOS) of atomic and molecular crystal phases of dense hydrogen at pressures up to 3.5 TPa are systematically investigated with density functional theory. The results indicate that the Vinet EOS model that fitted to low-pressure experimental data overestimates the compressibility of dense hydrogen drastically when beyond 500 GPa. Metastable multi-atomic molecular phases with weak covalent bonds are observed. When compressed beyond about 2.8 TPa, these exotic low-coordinated phases become competitive with the ground state and other high-symmetry atomic phases. Using nudged elastic band method, the transition path and the associated energy barrier between these high-pressure phases are evaluated. In particular for the case of dissociation of diatomic molecular phase into the atomic metallic Cs-IV phase, the existent barrier might raise the transition pressure about 200 GPa at low temperatures. Plenty of flat and broad basins on the energy surface of dense hydrogen ha...

Hybrid cluster expansions for local structural relaxations

A model is constructed in which pair potentials are combined with the cluster expansion method in order to better describe the energetics of structurally relaxed substitutional alloys. The effect of structural relaxations away from the ideal crystal positions and the effect of ordering is described by interatomic-distance-dependent pair potentials, while more subtle configurational aspects associated with correlations of three or more sites are described purely within the cluster expansion formalism. Implementation of such a hybrid expansion in the context of the cluster variation method or Monte Carlo method gives improved ability to model phase stability in alloys from first-principles.

Shock-induced order-disorder transformation in Ni3Al

The Hugoniot of Ni3Al with L12 structure is calculated with an equation of state (EOS) based on a cluster expansion and variation method from first principles. It is found that an order-disorder transition occurs at a shock pressure of 205GPa, corresponding to 3750K in temperature. On the other hand, an unexpected high melting temperature about 6955K is obtained at the same pressure, which is completely different from the case at ambient pressure where the melting point is slightly lower than the order-disorder transition temperature, implying the high pressure phase diagram has its own characteristics. The present work also demonstrates the configurational contribution is more important than electronic excitations in alloys and mineral crystals within a large range of temperature, and an EOS model based on CVM is necessary for high pressure metallurgy and theoretical Earth model.

Cluster expansion of electronic excitations: Application to fcc Ni–Al alloys

The cluster expansion method is applied to electronic excitations and a set of effective cluster densities of states (ECDOS) is defined, analogous to effective cluster interactions (ECIs). The ECDOSs are used to generate alloy thermodynamic properties as well as the equation of state (EOS) of electronic excitations for the fcc Ni-Al systems. When parent clusters have a small size, the convergence of the expansion is not so good but the electronic density of state (DOS) is well reproduced. However, the integrals of the DOS such as the cluster expanded free energy, entropy, and internal energy associated with electronic excitations are well described at the level of the tetrahedron-octahedron cluster approximation, indicating that the ECDOS is applicable to produce electronic ECIs for cluster variation method (CVM) or Monte Carlo calculations. On the other hand, the Gruneisen parameter, calculated with first-principles methods, is no longer a constant and implies that the whole DOS profile should be considered for EOS of electronic excitations, where ECDOS adapts very well for disordered alloys and solid solutions.

Order-disorder effects on the equation of state for fcc Ni-Al alloys

Order-disorder effects on equation of state EOS properties of substitutional binary alloys are investigated with the cluster variation method CVM based on ab initio effective cluster interactions ECI. Calculations are applied to the fcc based system. Various related quantities are shown to vary with concentration around stoichiometry with a surprising “W shape,” such as the thermal expansion coefficient, the heat capacity, and the Gruneisen parameter, due to configurational ordering effects. Analysis shows that this feature originates from the dominated behavior of some elements of the inverse of Hessian matrix, and relates to antisite defects occurring around stoichiometric compositions. This kind of strong compositional effects on EOS properties highlights the importance of subtle thermodynamic behavior of order-disorder systems.