Roman GrögerChainofstates atomistic calculations of the pathways and energetics of structural phase transformations in plutoniumUnderstanding structural phase transformations in plutonium is an important step toward obtaining thermodynamic and kinetic descriptions of this material. At ambient pressures, plutonium exists in six crystal structures that range from the highly symmetric facecentered cubic δ phase with four atoms per unit cell to one of the least symmetric structures: the monoclinic lowtemperature α phase with 16 atoms per unit cell. Just above room temperature, plutonium transforms to its bodycentered monoclinic β phase whose unit cell is experimentally predicted to consist of 34 atoms, an unusual number. Some insight into the nature of the phase transformations between various allotropes of plutonium has been gained recently using the Modified Embedded Atom Method (MEAM). These calculations not only correctly predict the order of the crystal structures and their stabilities but also reproduce the large 25% volume difference between the α and δ phases. However, the time scale required for observing a phase transformation between two phases is typically beyond the reach of conventional molecular dynamics simulations. To probe the structural phase transformations in plutonium, we use a chainofstates calculation implemented within the socalled nudged elastic band (NEB) method. We consider a chain of states {0, 1, 2, …, M+1} of M+2 replicas (images) of the system, that is a discrete representation of the transformation pathway. The atoms in the images 0 and M+1 are fixed and correspond to the initial and the target states of the system, while the remaining M movable images form the socalled elastic band. The atomic structure of each image is determined by the 3N positions of atoms, each of which is connected to its counterpart in the neighboring image by a linear spring to obtain a discrete representation of the transformation pathway. The coordinates of atoms in the M intermediate images are then determined by minimizing the total potential energy, represented by the objective function with respect to the positions of atoms in images. The potential energy E of every image i can be calculated readily by using empirical or semiempirical potentials, in our case the MEAM potential for plutonium. The relaxed configuration corresponds to zero forces on all atoms due to the interatomic potential and is largely independent of the choice of the spring constant k. The obtained coordinates of atoms in images 1 to M can be used to search systematically for the symmetry, space group, and Wyckoff positions by using the Rietveld refinement. Importantly, the variation of the potential energy E along the elastic band provides a direct measure of the energy barrier that has to be surmounted to transform between the given crystal structures.
The main goal of this project is to couple the NEB method with the semiempirical MEAM potential for plutonium to elucidate the mechanism of structural phase transformations in elemental plutonium and its dilute alloys. The current calculations aim to determine the energy barriers for α ↔ β, β ↔ γ, γ ↔ δ, and δ ↔ ε phase transformations that are required for reliable formulation of the activation enthalpies and the freeenergy functional of plutonium. The calculations of the energy barriers for α ↔ β with the β structure containing both 32 and 34 atoms and the subsequent comparison of the calculated photoemission spectra with experimental measurements will shed new light on the still controversial structure of the β phase.
Moreover, Los Alamos researchers Turab Lookman, Avadh Saxena, and Robert Albers of the Statistical Physics and Condensed Matter Group have recently proposed a phonon mechanism for the α( α') ↔ δ phase transformation in plutonium that is based purely on symmetry relations (group–subgroup) and which predicts an existence of intermediate trigonal and simple hexagonal phases. The existence of these intermediate structures, and thus the validity of this model, can be checked by examining the calculated transformation pathway between the α and δ phases. The same method will be used to investigate the influence of alloying on the transformation pathway and shape of the energy barrier between the α' and δ phases. Next: Material Disposition Program Overview

