Computational investigations of Dienes defect- and vacancy-induced changes in the electronic and vibrational properties of carbon fiber structural units.

Carbon fiber (CF) is a promising lightweight alternative to steel and is of significant interest for energy applications. As CF continues to find new uses and is exposed to new external conditions, a noninvasive method of monitoring its structural integrity is critical. Raman spectroscopy is a commonly used method for this monitoring; however, it is highly inferential, and the interpretation of the data is not always straightforward. In this work, we perform density functional theory (DFT) calculations to investigate changes in the vibrational properties of CF structural units (i.e., graphene and graphite) caused by monovacancy and Dienes defects as a foundation for modeling more complex defects that move our model toward that of realistic CF. Using large computational supercells, we can understand how these defects change the electronic structure and vibrational properties of graphene and graphite for interdefect distances near those of the lower experimental limit. The monovacancy opens an electronic bandgap at the K point. Although no such electronic gap is opened by the Dienes defect, both defects introduce flat defect bands near the Fermi energy. The Dienes defect creates long-range deviations of the phonons, leading to substantial broadening of the highest frequency optical modes in the band structure compared to that of the pristine material. In contrast, the phonon changes caused by the monovacancy are short range, and only minor changes in the band structure or phonon density of states were observed. These findings can assist in the interpretation of experimental results by providing atomic-scale insight into key electronic and vibrational features.

Computationally Guided Investigation of the Optical Spectra of Pure ß-UO3.

Single-phase ß-UO3 is synthesized by flash heating UO2(NO3)·6H2O in air to 450 °C and annealing for 60 h under the same conditions. For the first time, we report the Raman spectra of pure ß-UO3. To facilitate the assignment of Raman and infrared vibrational modes, we use density functional theory with density functional perturbation theory. By employing a novel analysis scheme that includes the mode frequencies as well as a quantitative analysis of the mode eigenvectors, we assign the observed spectral features to individual chemical modes. In particular, the density functional theory optimized structure, observed Raman spectrum, and eigenvector analysis suggest the presence of four crystallographically distinct uranyl ions, one more than has previously been suggested.

Interaction of hydrogen with actinide dioxide (011) surfaces.

The corrosion and oxidation of actinide metals, leading to the formation of metal-oxide surface layers with the catalytic evolution of hydrogen, impacts the management of nuclear materials. Here, the interaction of hydrogen with actinide dioxide (AnO2, An = U, Np, or Pu) (011) surfaces by Hubbard corrected density functional theory (PBEsol+U) has been studied, including spin-orbit interactions and non-collinear 3k anti-ferromagnetic behavior. The actinide dioxides crystalize in the fluorite-type structure, and although the (111) surface dominates the crystal morphology, the (011) surface energetics may lead to more significant interaction with hydrogen. The dissociative adsorption of hydrogen on the UO2 (0.44 eV), NpO2 (-0.47 eV), and PuO2 (-1.71 eV) (011) surfaces has been calculated. It is found that hydrogen dissociates on the PuO2 (011) surface; however, UO2 (011) and NpO2 (011) surfaces are relatively inert. Recombination of hydrogen ions is likely to occur on the UO2 (011) and NpO2 (011) surfaces, whereas hydroxide formation is shown to occur on the PuO2 (011) surface, which distorts the surface structure.

Elucidation of the Structure and Vibrational Spectroscopy of Synthetic Metaschoepite and Its Dehydration Product.

We confirm that synthetic uranyl hydroxide hydrate metaschoepite [(UO)24O(OH)6]·5H2O is unstable against dehydration under dry conditions, and we present a structural and vibrational spectroscopic study of synthetic metaschoepite and its ambient temperature dehydration product. Complementary structural (X-ray diffraction and neutron diffraction) and vibrational spectroscopic techniques (Raman spectroscopy, infrared spectroscopy, and inelastic neutron scattering) are used to probe different components of these species. Analysis of the dehydration product suggests that it contains both pentagonally coordinated and hexagonally coordinated uranyl ions, necessitating that some uranyl ions undergo a coordination change during the dehydration of pentagonally coordinated metaschoepite. Vibrational spectra of metaschoepite and its dehydration product are interpreted with power spectra generated from ab initio molecular dynamics trajectories, allowing assignment of all major features. We identify the uranyl symmetric stretching modes of the four distinct uranyl ions in synthetic metaschoepite and clarify the assignment of lower energy Raman modes in both structures. The coanalysis of experimental and computational data reveals a strong coupling between the uranyl stretching modes and hydroxide bending modes in the anhydrous structure, leading to the presence of several high-energy combination bands in the inelastic neutron scattering data.

Magnetic structure of UO2 and NpO2 by first-principle methods.

The magnetic structure of the actinide dioxides (AnO2) remains a field of intense research. A low-temperature experimental investigation of the magnetic ground-state is complicated by thermal energy released from the radioactive decay of the actinide nuclei. To establish the magnetic ground-state, we have employed high-accuracy computational methods to systematically probe different magnetic structures. A transverse 1k antiferromagnetic ground-state with Fmmm (No. 69) crystal symmetry has been established for UO2, whereas a ferromagnetic (111) ground-state with R3[combining macron]m (No. 166) has been established for NpO2. Band structure calculations have been performed to analyse these results.

Interaction of hydrogen with actinide dioxide (111) surfaces.

The interaction of atomic and molecular hydrogen with actinide dioxide (AnO2, An = U, Np, Pu) (111) surfaces has been investigated by DFT+U, where noncollinear 3k antiferromagnetic behaviour and spin-orbit interactions are considered. The adsorption of atomic hydrogen forms a hydroxide group, coupled to the reduction of an actinide ion. The energy of atomic hydrogen adsorption on the UO2 (0.82 eV), NpO2 (-0.10 eV), and PuO2 (-1.25 eV) surfaces has been calculated. The dissociation of molecular hydrogen is not observed, shown to be due to kinetic rather than thermodynamic factors. As a barrier to the formation of a second hydroxyl group, an unusual charge distribution has been shown. This could be a limitation of a (1·1) unit cell method or an artefact of the systems. The recombination of hydrogen ions on the AnO2 (111) surfaces is favoured over hydroxide formation.

Hidden magnetic order in plutonium dioxide nuclear fuel.

A thorough understanding of the chemistry of PuO2 is critical to the design of next-generation nuclear fuels and the long-term storage of nuclear materials. Despite over 75 years of study, the ground-state magnetic structure of PuO2 remains a matter of much debate. Experimental studies loosely indicate a diamagnetic (DM) ground-state, whereas theoretical methods have proposed either a collinear ferromagnetic (FM) or anti-ferromagnetic (AFM) ground-state, both of which would be expected to cause a distortion from the reported Fm3[combining macron]m symmetry. In this work, we have used accurate calculations based on the density functional theory (DFT) to systematically investigate the magnetic structure of PuO2 to resolve this controversy. We have explicitly considered electron-correlation, spin-orbit interaction and noncollinear magnetic contributions to identify a hereto unknown longitudinal 3k AFM ground-state that retains Fm3[combining macron]m crystal symmetry. Given the broad interest in plutonium materials and the inherent experimental difficulties of handling this compound, the results presented in this paper have considerable implications for future computational studies relating to PuO2 and related actinide structures. As the crystal structure is coupled by spin-orbit interactions to the magnetic state, it is imperative to consider relativity when creating computational models.

Investigating the role of phonons in the phase stability of uranium-based Laves phases.

Laves phase alloys possess unique thermal and electrical conduction properties, yet the factors governing phase stability in these systems remain an open question. The influence of phonons in particular has been broadly overlooked. Here, we investigate the UCo2x Ni2(1-x) chemical space using density functional theory, which offers a unique opportunity to explore the factors influencing Laves phase stability as all three primary Laves phases (C14, C15, C36) can be stabilized by changing the ratio of Co to Ni. Calculations of the thermodynamic and dynamical stability of pure UCo2 and UNi2 in each of three primary Laves phases confirm the stability of experimentally known Laves phases for UNi2 and UCo2. A decrease in bonding strength is identified in UNi2 compared to UCo2, aligned with redshifts observed in the UNi2 phonon density of states and a decoupling of the U and Ni vibrational modes. Phonon calculations of C14 UCo2 reveal dynamical instabilities. Efforts to remove the unstable mode at the Γ point in UCo2 via atomic displacements break the symmetry of the C14 phase, revealing a lower energy P2/c structure. Vibrational contributions to the free energy were calculated and did not change the thermodynamically stable Laves phase below 1000 K. The temperature-dependent free energies of single phase UCo2 and UNi2 were used to interpolate the relative stability of ternary UCo2x Ni2(1-x) in each of the three Laves phases at varying temperatures and stoichiometries. The ternary C36 phase is only predicted to be thermodynamically stable over a narrow stoichiometric range below 600 K.

Analysis of Water Coupling in Inelastic Neutron Spectra of Uranyl Fluoride.

Inelastic neutron scattering (INS) is uniquely sensitive to hydrogen due to its comparatively large thermal neutron scattering cross-section (82 b). Consequently, the inclusion of water in real samples presents significant challenges to INS data analysis due directly to the scattering strength of hydrogen. Here, we investigate uranyl fluoride (UO2F2) with inelastic neutron scattering. UO2F2 is the hydrolysis product of uranium hexafluoride (UF6), and is a hygroscopic, uranyl-ion containing particulate. Raman spectral signatures are commonly used for inferential understanding of the chemical environment for the uranyl ion in UO2F2, but no direct measurement of the influence of absorbed water molecules on the overall lattice dynamics has been performed until now. To deconvolute the influence of waters on the observed INS spectra, we use density functional theory with full spectral modeling to separate lattice motion from water coupling. In particular, we present a careful and novel analysis of the Q-dependent Debye-Waller factor, allowing us to separate spectral contributions by mass, which reveals preferential water coupling to the uranyl stretching vibrations. Coupled with the detailed partial phonon densities of states calculated via DFT, we infer the probable adsorption locations of interlayer waters. We explain that a common spectral feature in Raman spectra of uranyl fluoride originates from the interaction of water molecules with the uranyl ion based on this analysis. The Debye-Waller analysis is applicable to all INS spectra and could be used to identify light element contributions in other systems.

Formation of a uranyl hydroxide hydrate via hydration of [(UO2F2)(H2O)]7·4H2O.

Hydrated uranyl fluoride, [(UO2F2)(H2O)]7·4H2O, is not stable at elevated water vapor pressure, undergoing a complete loss of fluorine to form a uranyl hydroxide hydrate. Powder X-ray diffraction data of the resultant uranyl hydroxide species is presented for the first time, along with Raman and infrared (IR) spectra. The new uranyl hydroxide species is structurally similar to the layered uranyl hydroxide hydrate minerals schoepite and metaschoepite, but has a significantly expanded interlayer spacing (c = 15.12 vs. 14.73 Å), suggesting that additional H2O molecules may be present between the uranyl layers. Comparison of the Raman and IR spectra of this new uranyl hydroxide hydrate and synthetic metaschoepite ([(UO2)4O(OH)6]·5H2O) suggests that the equatorial environment of the uranyl ion may differ and that H2O molecules in the new species participate in stronger hydrogen bonds. In addition, the interlayer spacing of both this new uranyl hydroxide species and synthetic metaschoepite is shown to be sensitive to the environmental humidity, contracting and re-expanding with desiccation and rehydration. Structural distinction between the new uranyl hydroxide species and synthetic metaschoepite is confirmed by a comparison of the thermal behavior; unlike metaschoepite, the new hydrate does not form α-UO2(OH)2 upon dehydration.

Comparison of ab initio and DFT electronic structure methods for peptides containing an aromatic ring: effect of dispersion and BSSE.

We establish that routine B3LYP and MP2 methods give qualitatively wrong conformations for flexible organic systems containing pi systems and that recently developed methods can overcome the known inadequacies of these methods. This is illustrated for a molecule (a conformer of the Tyr-Gly dipeptide) for which B3LYP/6-31+G(d) and MP2/6-31+G(d) geometry optimizations yield strikingly different structures [Mol. Phys. 2006, 104, 559-570]: MP2 predicts a folded "closed-book" conformer with the glycine residue located above the tyrosine ring, whereas B3LYP predicts a more open conformation. By employing different levels of theory, including the local electron correlation methods LMP2 (local MP2) and LCCSD(T0) (local coupled cluster with single, double, and noniterative local triple excitations) and large basis sets (aug-cc-pVnZ, n=D, T, Q), it is shown that the folded MP2 minimum is an artifact caused by large intramolecular BSSE (basis set superposition error) effects in the MP2/6-31+G(d) calculations. The B3LYP functional gives the correct minimum, but the potential energy apparently rises too steeply when the glycine and tyrosine residues approach each other, presumably due to missing dispersion effects in the B3LYP calculations. The PWB6K and M05-2X functionals, designed to give good results for weak interactions, remedy this to some extent. The reduced BSSE in the LMP2 calculations leads to faster convergence with increasing basis set quality, and accurate results can be obtained with smaller basis sets as compared to canonical MP2. We propose LMP2 as a suitable method to study interactions with pi-electron clouds.