Atomistic modeling of the negative thermal expansion in δ-plutonium based on the two-state description

The δ phase of plutonium with the fcc structure exhibits an unusual negative thermal expansion (NTE) over its narrow temperature range of stability, 593-736 K. An accurate description of the anomalous high-temperature volume effect of plutonium goes beyond the current capability of electronic-structure calculations. We propose an atomistic scheme to model the thermodynamic properties ofδ-Pu based on the two-state model of Weiss for the Invar alloys, inspired by the simple free-energy analysis previously conducted by Lawson et al. The two-state mechanism is incorporated into the atomistic description of a many-body interacting system. Two modified embedded atom method potentials are employed to represent the binding energies of two competing electronic states inδ-Pu. We demonstrate how the NTE takes place inδ-Pu by means of Monte Carlo simulations implemented with the two-state mechanism.