Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time dependant density functional theory”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Using Neural Network Force Fields to Ascertain the Quality of Ab Initio Simulations of Liquid Water

Accurately simulating the properties of bulk water, despite the apparent simplicity of the molecule, is still a challenge. In order to fully understand and reproduce its complex phase diagram, it is necessary to perform simulations at the ab initio level, including quantum mechanical effects both for electrons and nuclei. This comes at a high computational cost, given that the structural and dynamical properties tend to require long timescales and large simulation cells. In this work, we evaluate the errors that density functional theory (DFT)-based simulations routinely incur into due time- and size-scale limitations. These errors are evaluated using neural-network-trained force fields that are accurate at the level of DFT methods. We compare different exchange and correlation potentials for properties of bulk water that require large timescales. Here we show that structural properties are less dependent on the system size and that dynamical properties such as the diffusion coefficient have a strong dependence on the simulation size and timescale. Our results facilitate comparisons of DFT-based simulation results with experiments and offer a path to discriminate between model and convergence errors in these simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Defect-Limited Mobility and Defect Thermochemistry in Mixed A-Cation Tin Perovskites: (CH3NH3)1-xCsxSnBr3

Hybrid organic-inorganic semiconductors crystallizing in the perovskite structure present a significant opportunity for realizing defect-tolerant semiconductors. In this work, we examine the solid solution, (CH3NH3)1-xCsxSnBr3, and identify the thermochemistry dictating the intrinsic carrier concentrations and how local structural distortions influences this electronic behavior. This family of compounds exhibits the expected systematic trend in decreasing optical gap with the cesium to methylammonium A-site mixing ratio, x, in the visible region. However, the carrier mobility, as determined from time-resolved microwave conductivity measurements, trends opposite to that expected from first-principles calculations combined with Boltzmann scattering theory calculations of the carrier mobility. We propose that this is a result of increasing carrier scattering with x in (CH3NH3)1-xCsxSnBr3 due to a significant increase in the carrier density with x. By examining the dependence of the carrier density as a function of x, we infer the compositional-dependence of the average enthalpy and nonconfigurational entropy per defect. While diffraction reveals a cubic aristotypic perovskite structure as a function of x at room temperature, the pair distribution functions obtained from synchrotron X-ray total scattering from these materials are better described by symmetry-adapted displacement modes of the Pm3m crystal structure, which we attribute to a large degree of anharmonic dynamics of the atom positions. This analysis shows that CH3NH3+-rich compositions retain mostly linear Sn-Br-Sn bonding environments through the displacements, while Cs+-rich compositions lead to more significantly bent Sn-Br-Sn environments. We propose that these bent bonding arrangements in Cs-rich compositions yield a higher propensity for defect formation. This also provides a rationale for carrier trapping that gives rise to anomalous microwave transients. Together, these results provide insight into the structure-dynamics-properties relationships in this highly anharmonic system with high amplitude atomic motions and low defect formation energies.

defect thermochemistry↗
DOE OSTI2024

Mixed Valence {Ni 2+ Ni 1+ } Clusters as Models of Acetyl Coenzyme A Synthase Intermediates

Acetyl coenzyme A synthase (ACS) catalyzes the formation and deconstruction of the key biological metabolite, acetyl coenzyme A (acetyl-CoA). The active site of ACS features a {NiNi} cluster bridged to a [Fe4S4] n+ cubane known as the A-cluster. The mechanism by which the A-cluster functions is debated, with few model complexes able to replicate the oxidation states, coordination features, or reactivity proposed in the catalytic cycle. In this work, we isolate the first bimetallic models of two hypothesized intermediates on the paramagnetic pathway of the ACS function. The heteroligated {Ni 2+ Ni 1+ } cluster, [K(12-crown-4) 2 ][1], effectively replicates the coordination number and oxidation state of the proposed “A red ” state of the A-cluster. Addition of carbon monoxide to [1] - allows for isolation of a dinuclear {Ni 2+ Ni 1+ (CO)} complex, [K(12-crown-2) n ][2] (n = 1–2), which bears similarity to the “A NiFeC ” enzyme intermediate. Structural and electronic properties of each cluster are elucidated by X-ray diffraction, nuclear magnetic resonance, cyclic voltammetry, and UV/vis and electron paramagnetic resonance spectroscopies, which are supplemented by density functional theory (DFT) calculations. Calculations indicate that the pseudo-T-shaped geometry of the three-coordinate nickel in [1] – is more stable than the Y-conformation by 22 kcal mol –1 , and that binding of CO to Ni 1+ is barrierless and exergonic by 6 kcal mol –1 . UV/vis absorption spectroscopy on [2] - in conjunction with time-dependent DFT calculations indicates that the square-planar nickel site is involved in electron transfer to the CO π*-orbital. Further, we demonstrate that [2] - promotes thioester synthesis in a reaction analogous to the production of acetyl coenzyme A by ACS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Local order of orthorhombic weberite-type Y 3 TaO 7 as determined by neutron total scattering and density functional theory calculations

Weberite-type oxides are a family of oxides with A 3 BO 7 stoichiometry that can adopt different space groups depending on chemical composition. There is a discrepancy in previous studies as to whether Y 3 TaO 7 is the orthorhombic C222 1 or the Ccmm space group. Here, we describe the short- and long-range structural properties of weberite-type Y 3 TaO 7 using neutron total scattering data, which has a high sensitivity to the oxygen sublattice. Simultaneous analysis of both short- and long-range structural data via conventional Rietveld and small box modelling demonstrates that Y 3 TaO 7 is best modeled as C222 1 . While the Ccmm describes equally well the long-range structure, pair distribution function analysis revealed that the local atomic configuration can be best modelled as C222 1 . This is corroborated by first-principles calculations that confirm the energetic preference of the C222 1 over Ccmm. Further, neutron total scattering data are reported for the first time for the Y 3 TaO 7 weberite-type ceramics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Vacancy-Dependent Diffusion Mechanism in Oxygen-Defective SrFeO 3 Perovskite Materials: First-Principles Density Functional Theory and Experimental Approach

Understanding oxygen diffusion at the atomic scale in SrFeO 3−δ perovskites is crucial for developing oxygen storage materials with optimal performance. Such materials are required to have high stability, corrosion resistance, and acceptable oxygen storage capacity at moderate operating temperatures and pressures. Here, in this study, we used first-principles density functional theory and thermogravimetric analysis to study the vacancy-dependent oxygen diffusion in oxygen-deficient SrFeO 3−δ (δ = 0, 0.065, 0.125, 0.25, 0.5) perovskites. The electronic structures, including the partial- and spin-resolved density of states, for different SrFeO 3−δ phases were calculated and compared with available experimental and theoretical results. By mapping the migration pathways, we investigated diffusion mechanisms and calculated the energy barriers for oxygen diffusion in cubic, orthorhombic, and brownmillerite phases of SrFeO 3−δ perovskites. Using the calculated energy barriers, we deduced the diffusion time scales and diffusion coefficients within SrFeO 3−δ . A diffusion coefficient on the order of 10 –8 m 2 /s was obtained for SrFeO 2.875 . We experimentally investigated the roles of temperature and oxygen partial pressures on the redox kinetics and deduced the kinetics rate and diffusion density, which agreed well with the calculated values for the density of diffusing oxygen vacancy in the lattice. Our results showed that the energy barrier tends to reduce at higher oxygen concentrations. Our results serve as an important guideline for designing oxygen storage materials with optimal redox kinetics.

chemical looping with oxygen uncoupling (CLOU)↗
DOE OSTI2021

Examination of How Well Long-Range-Corrected Density Functionals Satisfy the Ionization Energy Theorem

For this work, we calculated the vertical ionization energies (VIE) of 99 species in two ways to examine the accuracy of several long-range-corrected (LC) hybrid meta functionals in comparison with a gradient approximation (GA), global hybrids, and doubly hybrids. In the category of LC functionals, we examined both those with meta ingredients (i.e., that depend on the kinetic energy density) and those without them. The LC-hybrid meta functionals examined are M11, revM11, M11plus, and ωB97M-V. The reference data used to assess accuracy consist of 95 molecules and 4 atoms in the GW100 set. The two methods studied are the ΔSCF method (involving the difference of neutral and cation self-consistent field (SCF) energies) and the ionization energy theorem (involving the orbital energy of the highest occupied molecular orbital, HOMO). We calculated linear correlation coefficients (r 2 ) and mean absolute deviations (MADs) between each approach and the reference VIE value from the CCSD(T)/def2-TZVPP level of theory. We compared the new LC-hybrid meta calculations to calculations with the 10 functionals in a previous VIE study by Brémond et al. and to the calculations with LC-BLYP (LC-Becke, Lee–Yang–Parr), CAM-B3LYP (Coulomb-attenuating-method Becke-3-parameter Lee–Yang–Parr), LC-ωHPBE, and ωB97X-D. The results show that Minnesota LC-hybrid meta functionals have the smallest mean absolute deviation of ionization energy theorem VIEs with the reference data; the LC-ωHPBE functional also does quite well in this test. This is very encouraging and indicates that LC-hybrid meta functionals would be the best starting points for the tuning strategy that has been shown to be a very good procedure for improving time-dependent density functional calculations, and it also helps explain the good success of LC-hybrid meta functionals for molecular excitation energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Some Problems in Density Functional Theory

Though calculations based on density functional theory (DFT) are used remarkably widely in chemistry, physics, materials science, and biomolecular research and though the modern form of DFT has been studied for almost 60 years, some mathematical problems remain. From a physical science perspective, it is far from clear whether those problems are of major import. For context, we provide an outline of the basic structure of DFT as it is presented and used conventionally in physical sciences, note some unresolved mathematical difficulties with those conventional demonstrations, then pose several questions regarding both the time-independent and time-dependent forms of DFT that could benefit from attention in applied mathematics. Furthermore, progress on any of these would aid in development of better approximate functionals and in interpretation of DFT.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Density functional thermodynamic description of spin, phonon and displacement degrees of freedom in antiferromagnetic-to-paramagnetic phase transition in YNiO 3

This work herein demonstrates a direct density functional description of the finite-temperature thermodynamic properties of solids exhibiting phase transitions through positional and spin symmetry breaking degrees of freedom. A classic example addressed here is the rare-earth (R) nickelates RNiO 3 where the ground state is characterized by crystallographic and magnetic (e.g., antiferromagnetic) long-range order (LRO), whereas the higher temperature paramagnetic phase manifests a range of local spin and positional symmetry breaking motifs with short-range order (SRO). Unlike time-dependent simulations of spin and positional degrees of freedom, in the present work, phases are described via a superposition of static configurations constructed by populating a periodic base lattice supercell allowing for the formation of energy lowing distribution of positional and spin local motifs. The thermal populations of the configurations in such a superposition phase are obtained from the energy-minimized Density Functional Theory (DFT)-calculated partition functions at different temperatures. This approach offers flexible inclusion of different physical contributions to the free energy, such as elastic, electronic and phonon free energies, all obtained from the same underlying DFT total energy calculations of periodic structures. The thermodynamic and magnetic properties of both LRO and SRO crystallographic and spin phases, including antiferromagnetic (AFM) to paramagnetic (PM) Néel phase transition in YNiO3 are studied. Including spin and phonon contributions, we find a DFT-calculated Néel temperature to be 144 K in satisfactory agreement with the experimental value of 145 K; whereas omitting the phonon contribution, one obtains a Néel temperature of 81 K. We present phonon contributions to the DFT-calculated temperature-dependent SRO, heat capacities, and the polymorphous distribution of nonzero local magnetic moments in the PM phase. This approach thus extends to finite temperatures the symmetry-broken DFT description of both the AFM and PM phases, demonstrating that a thermodynamic superposition approach based on symmetry broken configurations evaluated by a mean-field like DFT is sufficient to obtain a consistent description of the thermal physics of the AFM, PM phases and their interconversion in 3d oxides illustrated by YNiO 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Diffusion-to-Imbibition Transition in Water Sorption in Nanoporous Media: Theoretical Studies

The ability to predict multiphase fluid transport in nanoporous rocks such as shales is critical for many geoscience applications, for example unconventional hydrocarbon production, geologic carbon sequestration, and nuclear waste disposal. When the pore sizes approach nanoscales, the impact of the molecular interaction forces between fluids and solids becomes increasingly important. These forces can alter macroscopic fluid phase behavior and control transport. Recent experimental studies have shown that capillary condensation and subsequent imbibition of liquid water can occur in hydrophilic nanoporous media even if the vapor phase is at a critical relative humidity (rhcrit) well below vapor saturation. This study presents a theoretical investigation of the processes controlling adsorption, capillary condensation and imbibition in nanoporous media, using the square-gradient classical density functional theory. The proposed theoretical model explicitly includes the relevant interaction forces among fluids and solids in macroscopic porous media. Application of the model to a relative-humidity-controlled water adsorption experiment is presented to demonstrate the impact of water-pore wall attractive forces on multiphase water behavior in a hydrophilic silicon nanoporous medium. Here, the model represents well the measured time-dependent evolution of the water imbibition front inside the nanoporous medium and also explains the diffusion-like water transport regimes observed at rh < rh crit and the imbibition-like flow regimes observed at rh > rh crit . The study furthermore gives an insight on hysteresis phenomenon in adsorption and desorption isotherms.

58 GEOSCIENCES↗
DOE OSTI2020

Solvent manipulation of the pre-reduction metal–ligand complex and particle-ligand binding for controlled synthesis of Pd nanoparticles

Understanding how to control the nucleation and growth rates is crucial for designing nanoparticles with specific sizes and shapes. Herein, we show that the nucleation and growth rates are correlated with the thermodynamics of metal–ligand/solvent binding for the pre-reduction complex and the surface of the nanoparticle, respectively. To obtain these correlations, we measured the nucleation and growth rates by in situ small angle X-ray scattering during the synthesis of colloidal Pd nanoparticles in the presence of trioctylphosphine in solvents of varying coordinating ability. The results show that the nucleation rate decreased, while the growth rate increased in the following order, toluene, piperidine, 3,4-lutidine and pyridine, leading to a large increase in the final nanoparticle size (from 1.4 nm in toluene to 5.0 nm in pyridine). Using density functional theory (DFT), complemented by 31 P nuclear magnetic resonance and X-ray absorption spectroscopy, we calculated the reduction Gibbs free energies of the solvent-dependent dominant pre-reduction complex and the solvent-nanoparticle binding energy. The results indicate that lower nucleation rates originate from solvent coordination which stabilizes the pre-reduction complex and increases its reduction free energy. At the same time, DFT calculations suggest that the solvent coordination affects the effective capping of the surface where stronger binding solvents slow the nanoparticle growth by lowering the number of active sites (not already bound by trioctylphosphine). The findings represent a promising advancement towards understanding the microscopic connection between the metal–ligand thermodynamic interactions and the kinetics of nucleation and growth to control the size of colloidal metal nanoparticles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Modeling dark photon oscillations in our inhomogeneous Universe

Not provided.A dark photon may kinetically mix with the Standard Model photon, leading to observable cosmological signatures. The mixing is resonantly enhanced when the dark photon mass matches the primordial plasma frequency, which depends sensitively on the underlying spatial distribution of electrons. Crucially, inhomogeneities in this distribution can have a significant impact on the nature of resonant conversions. We develop and describe, for the first time, a general analytic formalism to treat resonant oscillations in the presence of inhomogeneities. Our formalism follows from the theory of level crossings of random fields and only requires knowledge of the one-point probability density function (PDF) of the underlying electron number density fluctuations. We validate our formalism using simulations and illustrate the photon-to-dark photon conversion probability for several different choices of PDFs that are used to characterize the low-redshift Universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Statistical learning for predicting density–matrix-based electron dynamics

In this study, we consider the problem of learning density-dependent molecular Hamiltonian matrices from time series of electron density matrices, all in the context of Hartree–Fock theory. Prior work developed a solution to this problem for small molecular systems with density and Hamiltonian matrices of size at most 6 × 6. Here, using a battery of techniques, we scale prior methods to larger molecular systems with, for example, 29 × 29 matrices. This includes systems that either have more electrons or are expressed in large basis sets such as 6-311++G**. Scaling the method to larger systems enhances its relevance for realistic applications in chemistry and physics. To achieve this scaling, we apply dimensionality reduction, ridge regression and analytic computation of Hessians. Through the combination of these techniques, we are able to learn Hamiltonians by minimizing an objective function that encodes local propagation error. Importantly, these learned Hamiltonians can then be used to predict electron dynamics for thousands of steps: When we use our learned Hamiltonians to numerically solve the time-dependent Hartree–Fock equation, we obtain predicted dynamics that are in close quantitative agreement with ground truth dynamics. This includes field-off trajectories similar to the training data and field-on trajectories outside of the training data.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Dioxygen Activation Kinetics over Distinct Cu Site Types in Cu-Chabazite Zeolites

Cu-exchanged zeolites activate dioxygen to form active sites for partial methane oxidation (PMO), nitrogen oxide decomposition, and carbon monoxide oxidation. Apparent rates of O 2 activation depend both on the intrinsic kinetics of distinct Cu site types and the distributions of such sites within a given zeolite, which depend on the density and arrangement of the framework Al atoms. In this work, we use hydrothermal synthesis methods to control the arrangement of framework Al sites in chabazite (CHA) zeolites and, in turn, the distinct Cu site types formed. Time-resolved in situ Raman spectroscopy reveals the kinetics of O 2 adsorption and activation within these well-defined Cu-CHA materials and the concomitant structural evolution of copper–oxygen (Cu x O y ) complexes, which are interpreted alongside Cu(I) oxidation kinetics extracted from in situ X-ray absorption spectroscopy (XAS). Raman spectra of several plausible Cu x O y species simulated using density functional theory suggest that experimental spectra (λ ex = 532 nm) capture the formation of mono(μ-oxo)dicopper species (ZCuOCuZ). Transient experiments show that the timescales required to form Cu x O y structures that no longer change in Raman spectra correspond to the durations of oxidative treatments that maximize CH 3 OH yields in stoichiometric PMO cycles (approximately 2 h). Yet, these periods extend well beyond the timescales for the complete conversion of the initial Cu(I) intermediates to their Cu(II) states (<0.3 h, reflected in X-ray absorption near edge spectroscopy spectra), which demonstrates that Cu x O y complexes continue to evolve structurally following rapid Cu(I) oxidation. The dependence of ZCuOCuZ formation rates on O 2 pressure, H 2 O pressure, and temperature is consistent with a mechanism in which ZCuOH reduces to form ZCu + sites that bind molecular oxygen and form ZCu-O 2 intermediates. Subsequent reaction with proximate ZCu + forms bridging peroxo dicopper complexes that cleave O–O bonds to form ZCuOCuZ in steps facilitated by water. These data and interpretations provide evidence for the chemical processes that link rapid and kinetically irrelevant Cu oxidation steps (frequently probed by XAS and UV–vis spectroscopy) to the relatively slow genesis of reactive Cu complexes that form CH 3 OH during PMO. In doing so, we reveal previously unrecognized complexities in the processes by which Cu ions in zeolites activate O 2 to form active Cu x O y complexes, which underscore the insight afforded by judicious combinations of experimental and theoretical techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

First-principles calculation of the configurational energy density of states for a solid-state ion conductor with a variant of the Wang and Landau algorithm

In this work, a variant of the Wang and Landau algorithm for calculation of the configurational energy density of states is proposed. The algorithm was developed for the purpose of using first-principles simulations, such as density functional theory, to calculate the partition function of disordered sublattices in crystal materials. The expensive calculations of first-principles methods make a parallel algorithm necessary for a practical computation of the configurational energy density of states within a supercell approximation of a solid-state material. The algorithm developed in this work is tested with the two-dimensional (2d) Ising model to bench mark the algorithm and to help provide insight for implementation to a materials science application. Tests with the 2d Ising model revealed that the algorithm has good performance compared to the original Wang and Landau algorithm and the 1/$\textit{t}$ algorithm, in particular the short iteration performance. Further, a proof of convergence is presented within an adiabatic assumption, and the analysis is able to correctly predict the time dependence of the modification factor to the density of states. The algorithm was then applied to the lithium and lanthanum sublattice of the solid-state lithium ion conductor Li 0.5 La 0.5 TiO 3 . This was done to help understand the disordered nature of the lithium and lanthanum. The results find, overall, that the algorithm performs very well for the 2d Ising model and that the results for Li 0.5 La 0.5 TiO 3 are consistent with experiment while providing additional insight into the lithium and lanthanum ordering in the material. The primary result is that the lithium and lanthanum become more mixed between layers along the c axis for increasing temperature. In part, the simulation of the disordered Li 0.5 La 0.5 TiO 3 system serves as a benchmark for what size systems are currently and in the near future practical to calculate with density functional theory methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Biased tracers in redshift space in the EFTofLSS with exact time dependence

We study the effect of the Einstein-de Sitter (EdS) approximation on the one-loop power spectrum of galaxies in redshift space in the Effective Field Theory of Large-Scale Structure. The dark matter density perturbations and velocity divergence are treated with exact time dependence. Splitting the density perturbation into its different temporal evolutions naturally gives rise to an irreducible basis of biases. While, as in the EdS approximation, at each time this basis spans a seven-dimensional space, this space is a slightly different one, and the difference is captured by a single calculable time- and vec k-dependent function. We then compute the redshift-space galaxy one-loop power spectrum with the EdS approximation (P EdS-approx ) and without (P Exact ). For the monopole we find P 0 Exact /P 0 Ed S-approx ~ 1.003 and for the quadrupole P 2 Exact /P 2 EdS-approx ~ 1.007 at z=0.57, and sharply increasing at lower redshifts. Finally, we show that a substantial fraction of the effect remains even after allowing the bias coefficients to shift within a physically allowed range. This suggests that the EdS approximation can only fit the data to a level of precision that is roughly comparable to the precision of the next generation of cosmological surveys. Furthermore, we find that implementing the exact time dependence formalism is not demanding and is easily applicable to data. Both of these points motivate a direct study of this effect on the cosmological parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Correlating Crystallography, Magnetism, and Electronic Structure Across Anhysteretic First-Order Phase Transition in Pr 2 In

Here, temperature-dependent powder X-ray diffraction and magnetization measurements of Pr 2 In conclusively prove that the unusual anhysteretic first-order paramagnetic-ferromagnetic phase transition in the compound is related to concurrent changes in both the magnetic and crystallographic lattices. At the same time, the hexagonal Ni 2 In-type structure is stable at least between 6 and 298 K, including at T C = ~57 K. From the density functional theory calculations, the electronic structure of the compound is extraordinarily sensitive to minor changes in lattice parameters that occur across the phase transition, revealing the origin of strong magnetoelastic coupling. In the vicinity of T C , the maximum entropy change, ΔS Max = -16 J Kg -1 K -1 induced by a moderate magnetic field change of 20 kOe (ΔS Max = -20 J Kg -1 K -1 for 50 kOe magnetic field change) is comparable to other known potentially functional materials demonstrating large cryogenic magnetocaloric effect.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Realization of Organocerium-Based Fullerene Molecular Materials Showing Mott Insulator-Type Behavior

Abstract Electron-rich organocerium complexes (C 5 Me 4 H) 3 Ce and [(C 5 Me 5 ) 2 Ce(ortho-oxa)], with redox potentials E 1/2 = –0.82 V and E 1/2 = –0.86 V versus Fc/Fc + respectively, were reacted with fullerene (C 60 ) in different stoichiometries to obtain molecular materials. Structurally characterized co-crystals: [(C 5 Me 4 H) 3 Ce] 2 •C 60 (1) and [(C 5 Me 5 ) 2 Ce(ortho-oxa)] 3 •C 60 (2) of C 60 with cerium-based rare earth molecular precursors are reported for the first time. The extent of charge transfer in 1 and 2 was evaluated using a series of physical measurements: FT-IR, Raman, solidstate UV-vis-NIR spectroscopy, X-ray absorption near edge structure (XANES) spectroscopy, and magnetic susceptibility measurements. The physical measurements indicate that 1 and 2 comprise the cerium(III) oxidation state with formally neutral C 60 as a co-crystal in both cases. Pressure-dependent periodic density functional theory calculations were performed to study the electronic structure of 1. Inclusion of a Hubbard-U parameter removes Ce f states from the Fermi level, opens up a band gap, and stabilizes FM/AFM magnetic solutions that are isoenergetic because of the large distances between the Ce(III) cations. Furthermore, the electronic structure of this strongly correlated Mott insulator-type system is reminiscent of the well-studied Ce 2 O 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Fission Fragment Intrinsic Spins and Their Correlations

The intrinsic spins and their correlations are the least understood characteristics of fission fragments from both theoretical and experimental points of view. In many nuclear reactions, the emerging fragments are typically excited and acquire an intrinsic excitation energy and an intrinsic spin depending on the type of the reactions and interaction mechanism. Both the intrinsic excitation energies and the fragments’ intrinsic spins and parities are controlled by the interaction mechanism and conservations laws, which lead to their correlations and determines the character of their deexcitation mechanism. We outline here a framework for the theoretical extraction of the intrinsic spin distributions of the fragments and their correlations within the fully microscopic real-time density-functional theory formalism and illustrate it on the example of induced fission of 236 U and 240 Pu, using two nuclear energy density functionals. Here, these fission fragment intrinsic spin distributions display new qualitative features previously not discussed in literature. Within this fully microscopic framework, we extract for the first time the intrinsic spin distributions of fission fragments of 236 U and 240 Pu as well as the correlations of their intrinsic spins, which have been debated in literature for more than six decades with no definite conclusions so far.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗