Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “invariant”

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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Hamiltonian reduction of Vlasov–Maxwell to a dark slow manifold

In this work, we show that non-relativistic scaling of the collisionless Vlasov–Maxwell system implies the existence of a formal invariant slow manifold in the infinite-dimensional Vlasov–Maxwell phase space. Vlasov–Maxwell dynamics restricted to the slow manifold recovers the Vlasov–Poisson and Vlasov–Darwin models as low-order approximations, and provides higher-order corrections to the Vlasov–Darwin model more generally. The slow manifold may be interpreted to all orders in perturbation theory as a collection of formal Vlasov–Maxwell solutions that do not excite light waves, and are therefore ‘dark’. We provide a heuristic lower bound for the time interval over which Vlasov–Maxwell solutions initialized optimally near the slow manifold remain dark. We also show how the dynamics on the slow manifold naturally inherits a Hamiltonian structure from the underlying system. After expressing this structure in a simple form, we use it to identify a manifestly Hamiltonian correction to the Vlasov–Darwin model. The derivation of higher-order terms is reduced to computing the corrections of the system Hamiltonian restricted to the slow manifold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A gauge-compatible Hamiltonian splitting algorithm for particle-in-cell simulations using finite element exterior calculus

A particle-in-cell algorithm is derived with a canonical Poisson structure in the formalism of finite element exterior calculus. The resulting method belongs to the class of gauge-compatible splitting algorithms, which exactly preserve gauge symmetries and their associated conservation laws via the momentum map. We numerically demonstrate this time invariance of the momentum map and its usefulness in establishing precise initial conditions with a desired initial electric field and fixed background charge. The restriction of this canonical, finite element Poisson structure to the 1X2P $1\frac {1}{2}$ -dimensional phase space is also considered and simulated numerically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Finite orbit width effects in large aspect ratio stellarators

New orbit-averaged equations for low collisionality neoclassical fluxes in large aspect ratio stellarators with mirror ratios close to unity are derived. The equations retain finite orbit width effects by employing the second adiabatic invariant $J$ as a velocity-space coordinate and they have been implemented in the orbit-averaged neoclassical code KNOSOS (Velasco et al. , J. Comput. Phys. , vol. 418, 2020, 109512; Velasco et al. , Nucl. Fusion , vol. 61, 2021, 116013). The equations are used to study the $1/\nu$ regime and the lower collisionality regimes. For generic large aspect ratio stellarators with mirror ratios close to unity, as the collision frequency decreases, the $1/\nu$ regime transitions directly into the $\nu$ regime, without passing through a $\sqrt {\nu }$ regime. An explicit formula for the neoclassical fluxes in the $\nu$ regime is obtained. The formula includes the effect of particles that transition between different types of wells. While these transitions produce stochastic scattering independent of the value of the collision frequency in velocity space, the diffusion in real space remains proportional to the collision frequency. The $\sqrt {\nu }$ regime is only recovered in large aspect ratio stellarators close to omnigeneity: large aspect ratio stellarators with large mirror ratios and optimized large aspect ratio stellarators with mirror ratios close to unity. Neoclassical transport in large aspect ratio stellarators with large mirror ratios can be calculated with the orbit-averaged equations derived by Calvo et al. ( Plasma Phys. Control. Fusion , vol. 59, 2017, 055014). In these stellarators, the $\sqrt {\nu }$ regime exists in the collisionality interval $(a/R) \ln (R/a) \ll \nu _{ii} R a/\rho _i v_{ti} \ll R/a$ . In optimized large aspect ratio stellarators with mirror ratios close to unity, the $\sqrt {\nu }$ regime occurs in an interval of collisionality that depends on the deviation from omnigeneity $\delta$ : $\delta ^{2} |\ln \delta | \ll \nu _{ii} R a/\rho _i v_{ti} \ll 1$ . Here, $\nu _{ii}$ is the ion–ion collision frequency, $\rho _i$ and $v_{ti}$ are the ion gyroradius and thermal speed, and $a$ and $R$ are the minor and major radii.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Detecting shearless phase-space transport barriers in global gyrokinetic turbulence simulations with test particle map models

In magnetically confined fusion plasmas, the role played by zonal 𝐸 x 𝐵 flow shear layers in the suppression of turbulent transport is relatively well understood. However, less is understood about the role played by the weak shear regions that arise in the non-monotonic radial electric field profiles often associated with these shear layers. In electrostatic simulations from the global total-𝑓 gyrokinetic particle-in-cell code XGC, we demonstrate how shearless regions with non-zero flow curvature form zonal ‘jets’ that, in conjunction with neighbouring regions of shear, can act as robust barriers to particle transport and turbulence spreading. By isolating quasi-coherent fluctuations radially localised to the zonal jets, we construct a map model for the Lagrangian dynamics of gyrokinetic test particles in the presence of drift waves. We identify the presence of shearless invariant tori in this model and verify that these tori act as partial phase-space transport barriers in the simulations. We also demonstrate how avalanches impinging on these shearless tori cause eddy detachment events that form ‘cold/warm core ring’ structures analogous to those found in oceanic jets, facilitating transport across the barriers without destroying them completely. We discuss how shearless tori may generically arise from tertiary instabilities or other types of discrete eigenmodes, suggesting their potential relevance to broader classes of turbulent fluctuations.

fusion plasma↗
DOE OSTI2021

Antiferromagnetic Order and Linear Magnetoresistance in Fe-Substituted Shandite Co 3 In 2 S 2

Co 3 In 2 S 2 is a shandite phase closely related to the ferromagnetic Weyl semimetal Co 3 Sn 2 S 2 and contains Co atoms on an ideal kagome net. Elemental substitutions in Co 3 Sn 2 S 2 have been shown to invariably suppress its magnetism, and it remains a rare example of a shandite phase with magnetic order. Here, we report a new magnetic shandite compound discovered by substituting Fe into Co 3 In 2 S 2 . Single crystals of (Co 1–x Fe x ) 3 In 2 S 2 were grown, and neutron diffraction, ac and dc magnetization, and heat capacity measurements reveal long-range antiferromagnetic order with an interesting noncollinear structure and transition temperatures from 11 to 17 K. Partially substituting Fe for Co was inspired by calculated electronic structures, which suggest that Co 3 In 2 S 2 is near a magnetic instability and shows a peak in the density of states that can be approached by hole doping. While the focus of the present work is primarily on magnetism, there is a nodal ring below the Fermi energy in Co 3 In 2 S 2 that may also provide interest from a topological perspective. Potentially related to this, we find large, linear, nonsaturating magnetoresistance in this compound. Furthermore, this work expands the scope of magnetic topological materials among shandites and motivates similar chemical manipulation studies to reveal hidden magnetic instabilities in this important family and other topological systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electrolyte-Induced Restructuring of Acid-Stable Oxygen Evolution Catalysts

Crystalline metal oxide catalysts operating under oxygen evolution reaction (OER) conditions invariably restructure, resulting in active sites with hydroxo/oxo species in an amorphous environment. An increase in the population of terminal hydroxo/oxo species (i.e., edge sites) facilitates proton-coupled electron-transfer (PCET) kinetics for oxygen generation and thus improves catalyst competency. While amorphous films benefit from a greater density of active sites, they suffer from diminished charge transport as compared to that of extended crystalline lattices. Managing this amorphous–crystalline dichotomy is essential when designing OER catalysts, which we highlight with the examination of electrodeposited PbO x materials, which historically are very poor OER catalysts. Along these lines, the presence of phosphate during PbO x electrodeposition truncates the growth of an extended lattice owing to its strong bonding to oxide surfaces to afford an amorphous catalyst film (A-PbO x ) with significant charge-transfer resistance (138 ± 42 Ω) and poor OER kinetics (420 ± 105 mV dec –1 Tafel slope). Conversely, electrodeposition of Pb 2+ in the presence of less coordinating electrolytes such as nitrate affords crystalline β-PbO2 with improved charge-transfer resistance (42.6 ± 1.1 Ω), though still poor OER kinetics (134 ± 36 mV dec –1 Tafel slope). By operating amorphous A-PbOx in less coordinating electrolytes, however, a new partially crystalline material can be generated (μc-PbO x ) with further reduced charge-transfer resistance (33.0 ± 1.4 Ω) and improved OER kinetics (70 ± 15 mV dec –1 Tafel slope). The enhanced OER activity of μc-PbO x is the result of coupling the high edge-site population of an amorphous PbOx phase with crystalline-like charge transport properties. Finally, the ability to use an electrolyte to induce OER activity in an inactive amorphous form of PbO x highlights the benefits of optimizing the amorphous–crystalline phase compositions in the design of active OER catalysts.

catalysts↗
DOE OSTI2022

Synergistic Enhancement of Lead and Selenate Uptake at the Barite (001)–Water Interface

The interactions of heavy metals with minerals influence the mobility and bioavailability of toxic elements in natural aqueous environments. The sorption of heavy metals on covalently bonded minerals is generally well described by surface complexation models (SCMs). However, understanding sorption on sparingly soluble minerals is challenging because of the dynamically evolving chemistry of sorbent surfaces. The interpretation can be even more complicated when multiple metal ions compete for sorption. In the present study, we observed synergistically enhanced uptake of lead and selenate on the barite (001) surface through two sorption mechanisms: lattice incorporation that dominates at lower coverages and two-dimensional monolayer growth that dominates at higher coverages. Further, we also observed a systematic increase in the sorption affinity with increasing co-sorbed ion coverages, different from the assumption of invariant binding constants for individual adsorption processes in classical SCMs. Computational simulations showed thermodynamically favorable co-incorporation of lead and selenate by simultaneously substituting for barium and sulfate in neighboring sites, resulting in the formation of molecular clusters that locally match the net dimension of the substrate lattice. These results emphasize the importance of ion–ion interactions at mineral–water interfaces that control the fate and transport of contaminants in the environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Investigation of Mass Transfer and Sorption in CO 2 /Brine/Rock Systems via In Situ FT-IR

CO 2 geological storage in deep saline formations is considered a promising method to mitigate anthropogenic CO 2 emissions and, thereby, minimize changes to the Earth’s atmosphere. A fundamental understanding of CO 2 mass transfer and sorption phenomena in brine-saturated reservoir formations is necessary to understand the long-term fate of injected CO 2 as it is subjected to different (physical, dissolution, and mineral) trapping mechanisms. In this work, we investigate CO 2 sorption in brine-saturated and dry Mt. Simon sandstone samples via in-situ Fourier Transform infrared spectroscopy (FT-IR) at elevated pressures, ranging from 0.3 MPa to 8.3 MPa, at a temperature of 50 ⁰C. The FT-IR spectra of bulk-phase CO 2 were simultaneously recorded under the same conditions. For bulk-phase CO 2 , we observed, in agreement with past studies, a doublet peak at 2361 cm -1 and 2336 cm -1 and another peak (ν 2 bending mode) at 667 cm -1 . With increasing pressure, the position of the peak at 667 cm -1 remains invariant, however, when crossing into the supercritical region the doublet peak degenerates onto a single peak at 2336 cm -1 with a barely visible shoulder at 2361 cm -1 . The bulk CO 2 data provide a perfect fit for Beer’s law for the whole range of pressure conditions. For the dry sample, the IR spectrum is experimentally indistinguishable from the bulk CO 2 spectrum, signifying that if physical adsorption occurs to any significant extent, the adsorbed CO 2 molecules are not substantially more rotationally constrained than the dense bulk CO 2 molecules. For the brine-saturated sample, we observe a strong band centered at 2342 cm -1 and a small companion peak at 2360-2361 cm -1 that degenerates into a barely visible shoulder peak at the higher pressures. The 2342 cm -1 band has been previously observed by other investigators for CO 2 dissolved in bulk water/brine as well during its adsorption on a variety of other wet natural porous media. We observe no peaks corresponding to bicarbonate or carbonate bulk species, which correlates well with the prior literature on similar low-pH aqueous solutions. The integrated peak area for the CO 2 sorbed in the brine-saturated sample correlates linearly with its solubility in the same bulk brine, as measured separately via a PVT-cell approach. Here, this validates the accuracy of both techniques, and the potential of the FT-IR method to be used in the study of mass transfer and adsorption in such systems. To that effect, a simple mathematical model is presented to analyze the FT-IR data to determine the CO 2 effective diffusivity in the brine-saturated sandstone sample.

58 GEOSCIENCES↗
DOE OSTI2022

Four-Coordinate Fe N 2 and Imido Complexes Supported by a Hemilabile NNC Heteroscorpionate Ligand

Inspired by mechanistic proposals for N 2 reduction at the nitrogenase FeMo cofactor, we report herein a new, strongly σ-donating heteroscorpionate ligand featuring two weak-field pyrazoles and an alkyl donor. This ligand supports four-coordinate Fe(I)-N 2 , Fe(II)-Cl, and Fe(III)-imido complexes, which we have characterized using a variety of spectroscopic and computational methods. Structural and quantum mechanical analysis reveal the nature of the Fe–C bonds to be essentially invariant between the complexes, with conversion between the (formally) low-valent Fe-N 2 and high-valent Fe-imido complexes mediated by pyrazole hemilability. This presents a useful strategy for substrate reduction at such low-coordinate centers and suggests a mechanism by which FeMoco might accommodate the binding of both π-acidic and π-basic nitrogenous substrates.

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

Periodic Trends within Actinyl(VI) Nitrates and Their Structures, Vibrational Spectra, and Electronic Properties

In this work, a series of actinyl(VI) nitrate salts of the form MAnO 2 (NO 3 ) 3 , where M = NH 4 + K + , Rb + , Cs + , and Me 4 N + and AnO 2 2+ = U, Np, Pu, and AnO 2 (NO 3 ) 2 (H 2 O) 2 ·H 2 O, and the uranyl tetranitrates M 2 UO 2 (NO 3 ) 4 have been synthesized from aqueous solution and their structures determined using single-crystal X-ray diffraction. Together, these complexes represent an isostructural series of actinide complexes among the salts crystallized with the same charge-compensating cation and have been studied using vibrational spectroscopy including Raman and Fourier-transform infrared. Periodic trends in both the structural properties of these complexes and their vibrational spectra are presented and discussed, in particular the invariant nature of the O≡An≡O asymmetric stretching frequencies observed across the actinyl series. Electronic structure calculations were performed at a variety of levels of theory to aid in the interpretation of the vibrational data and to correlate trends in the data with the underlying electronic properties of these molecules.

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

Commutative Algebra Modeling in Materials Science – A Case Study on Metal–Organic Frameworks (MOFs)

Metal-organic frameworks (MOFs) are a class of important crystalline and highly porous materials whose hierarchical geometry and chemistry hinder interpretable predictions in materials properties. Commutative algebra is a branch of abstract algebra that has been rarely applied in data and material sciences. We introduce the first ever commutative algebra modeling and prediction in materials science. Specifically, category-specific commutative algebra (CSCA) is proposed as a new framework for MOF representation and learning. It integrates element-based categorization with multiscale algebraic invariants to encode both local coordination motifs and global network organization of MOFs. These algebraically consistent, chemically aware representations enable compact, interpretable, and data efficient modeling of MOF properties such as Henry’s constants and uptake capacities for common gases. Compared to traditional geometric and graph-based approaches, CSCA achieves comparable or superior predictive accuracy while substantially improving interpretability and stability across data sets. By aligning commutative algebra with the chemical hierarchy, the CSCA establishes a rigorous and generalizable paradigm for understanding structure and property relationships in porous materials and provides a nonlinear algebra-based framework for data-driven material discovery.

Khaemba, Caleb S.↗
DOE OSTI2023

Accurate calculation of many-body energies in water clusters using a classical geometry-dependent induction model

We incorporate geometry-dependent distributed multipole and polarizability surfaces into an induction model that is used to describe the 3- and 4-body terms of the interaction between water molecules. The expansion is carried out up to hexadecapole with the multipoles distributed on the atom sites. Dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole distributed polarizabilities are used to represent the response of the multipoles to an electric field. We compare the model against two large databases consisting of 43,844 3-body terms and 3,603 4-body terms obtained from high level ab initio calculations previously used to fit the MB-pol and q-AQUA interaction potentials. The classical induction model with no adjustable parameters reproduces the ab-initio 3- and 4-body terms contained in these teo Databases with a Root-Mean-Square-Error (RMSE) of 0.104/0.058 and a Mean-Absolute-Error (MAE) of 0.054/0.026 kcal/mol, respectively, results that are on a par with those obtained 1 by fitting the same data using tens of thousands of Permutationally Invariant Polynomials (PIPs). This demonstrates the accuracy of this physically motivated model in describing the 3- and 4-body terms in the interactions between water molecules with no adjustable parameters. The triple-dipole-dispersion energy was included in the 3-body energy and was found to be small but not quite negligible. The model represents a practical, efficient and transferable approach for obtaining accurate non-additive interactions for multi-component systems without the need of performing tens of thousands of high level electronic structure calculations and fitting them with tens of thousands of PIPs.

molecular analysis water, aqueous↗
DOE OSTI2025

Identifying Band Inversions in Topological Materials Using Diffusion Monte Carlo

Topological insulators are characterized by insulating bulk states and robust metallic surface states. Band inversion is a hallmark of topological insulators. At time-reversal invariant points in the Brillouin zone, spin–orbit coupling (SOC) induces a swapping of orbital character at the bulk band edges. Reliably detecting band inversion in solid-state systems with many-body methods would aid in identifying possible candidates for spintronics and quantum computing applications and improve our understanding of the physics behind topologically nontrivial systems. Density functional theory (DFT) methods are a well-established means of investigating these interesting materials due to their favorable balance of computational cost and accuracy but often struggle to accurately model the electron–electron correlations present in the many materials containing heavier elements. In this work, we develop a novel method to detect band inversion within continuum quantum Monte Carlo (QMC) methods that can accurately treat the electron correlation and spin–orbit coupling that are crucial to the physics of topological insulators. Our approach applies a momentum-space-resolved atomic population analysis throughout the first Brillouin zone utilizing the Löwdin method and the one-body reduced density matrix produced with diffusion Monte Carlo (DMC). We integrate this method into QMCPACK, an open source ab initio QMC package, so that these ground-state methods can be used to complement experimental studies and validate prior DFT work on predicting the band structures of correlated topological insulators. Here, we demonstrate this new technique on the topological insulator bismuth telluride, which displays band inversion between its Bi-p and Te-p states at the Γ-point. We show an increase in charge on the bismuth-p orbital and a decrease in charge on the tellurium-p orbital when comparing band structures with and without SOC. Additionally, we use our method to compare the degree of band inversion present in monolayer Bi 2 Te 3 , which has no interlayer van der Waals interactions, to that seen in the bilayer and bulk. The method presented here will enable future many-body studies of band inversion that can shed light on the delicate interplay between correlation and topology in correlated topological materials.

Band structure↗
DOE OSTI2025

Modeling Strong Light-Matter Coupling in Correlated Systems: State-Averaged Cavity Quantum Electrodynamics Complete Active Space Self-Consistent Field Theory

The description of strongly correlated systems interacting with quantized cavity modes poses significant theoretical challenges due to the combinatorial scaling of electronic and photonic degrees of freedom. Recent advances addressing this complexity include cavity quantum electrodynamics (QED) generalizations of complete active space configuration interaction and density matrix renormalization group methods. In this work, we introduce a QED extension of state-averaged complete active space self-consistent field theory, which incorporates cavity-induced correlations through a second-order orbital optimization framework with robust convergence properties. The method is implemented using both photon number state and coherent state representations, with the latter showing robust origin invariance in the energies regardless of the completeness of the photonic Fock space. The implementation enables symmetry-free orbital relaxations to account for photon-mediated symmetry breaking in polaritonic systems. Numerical validation on lithium hydride, hydroxide anion, and magnesium hydride cation demonstrates that this method achieves significantly improved accuracy in modeling ground-state and polariton potential energy surfaces compared to QED-CASCI in a fixed orbital basis. In these studies, we reach sub-kcal/mol accuracy in potential energy surface in much smaller active spaces than are required for QED-CASCI. This advancement provides a more robust approach for studying cavity-altered chemical landscapes for ground and exited strongly coupled systems.

CASSCF↗
DOE OSTI2025

Field-Driven Simulations to Probe the Impact of Ionic Correlations on Solution Transport Coefficients in Binary, Ternary, and Reciprocal Quaternary Aqueous Electrolytes

Ionic correlations play a critical role in governing transport properties of mixed-salt aqueous electrolyte solutions, yet their quantitative characterization remains challenging, particularly for multicomponent aqueous systems. Here, we develop a general nonequilibrium molecular dynamics framework to efficiently compute Onsager transport coefficients and ionic correlations in mixed-salt aqueous solutions. Using field-driven simulations, we obtain accurate Onsager matrices for LiCl/ KCl, KCl/KBr, and LiBr/KCl electrolyte solutions with significantly reduced computational cost relative to equilibrium Green−Kubo methods. The framework enables direct assessment of how attractive cation−anion and repulsive like-ion interactions contribute to conductivity and salt diffusivity across compositions. While static ion pairing exhibits strong composition dependence, dynamic ion correlations remain nearly invariant, leading to constant deviations from Nernst−Einstein predictions. These results highlight the disconnect between static ion association and dynamic transport correlations, and they establish a transferable approach for analyzing ion transport in complex electrolyte environments relevant to separation processes and electrochemical systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Data-Driven Many-Body Models for Molecular Fluids: CO 2 /H 2 O Mixtures as a Case Study

In this study, we extend the scope of the many-body TTM-nrg and MB-nrg potential energy functions (PEFs), originally introduced for halide ion–water and alkali-metal ion–water interactions, to the modeling of carbon dioxide (CO2) and water (H2O) mixtures as prototypical examples of molecular fluids. Both TTM-nrg and MB-nrg PEFs are derived entirely from electronic structure data obtained at the coupled cluster level of theory and are, by construction, compatible with MB-pol, a many-body PEF that has been shown to accurately reproduce the properties of water. Although both TTM-nrg and MB-nrg PEFs adopt the same functional forms for describing permanent electrostatics, polarization, and dispersion, they differ in the representation of short-range contributions, with the TTM-nrg PEFs relying on conventional Born–Mayer expressions and the MB-nrg PEFs employing multidimensional permutationally invariant polynomials. By providing a physically correct description of many-body effects at both short and long ranges, the MB-nrg PEFs are shown to quantitatively represent the global potential energy surfaces of the CO 2 –CO 2 and CO 2 –H 2 O dimers and the energetics of small clusters, as well as to correctly reproduce various properties in both gas and liquid phases. Furthermore, building upon previous studies of aqueous systems, our analysis provides further evidence for the accuracy and efficiency of the MB-nrg framework in representing molecular interactions in fluid mixtures at different temperature and pressure conditions.

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

CASSCF with Extremely Large Active Spaces Using the Adaptive Sampling Configuration Interaction Method

The complete active space self-consistent field (CASSCF) method is the principal approach employed for studying strongly correlated systems. However, exact CASSCF can only be performed on small active spaces of ~20 electrons in ~20 orbitals due to exponential growth in the computational cost. Here, we show that employing the Adaptive Sampling Configuration Interaction (ASCI) method as an approximate Full CI solver in the active space allows CASSCF-like calculations within chemical accuracy (<1 kcal/mol for relative energies) in active spaces with more than ~50 active electrons in ~50 active orbitals, significantly increasing the sizes of systems amenable to accurate multiconfigurational treatment. The main challenge with using any selected CI-based approximate CASSCF is the orbital optimization problem; they tend to exhibit large numbers of local minima in orbital space due to their lack of invariance to active–active rotations (in addition to the local minima that exist in exact CASSCF). We highlight methods that can avoid spurious local extrema as a practical solution to the orbital optimization problem. We employ ASCI-SCF to demonstrate a lack of polyradical character in moderately sized periacenes with up to 52 correlated electrons and compare against heat-bath CI on an iron porphyrin system with more than 40 correlated electrons.

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

Interplay between London Dispersion, Hubbard U , and Metastable States for Uranium Compounds

High-throughput computational studies of lanthanide and actinide chemistry with density-functional theory are complicated by the need for Hubbard U corrections, which ensure localization of the f-electrons, but can lead to metastable states. This work presents a systematic investigation of the effects of both Hubbard U value and metastable states on the predicted structural and thermodynamic properties of four uranium compounds central to the field of nuclear fuels: UC, UN, UO 2 , and UCl 3 . We also assess the impact of the exchange-hole dipole moment (XDM) dispersion correction on the computed properties. Overall, the choice of Hubbard U value and inclusion of a dispersion correction cause larger variations in the computed geometric properties than result from metastable states. Further, the weak dependence of structure optimization on metastable states should simplify future high-throughput calculations on actinides. Conversely, addition of the dispersion correction is found to offset the repulsion introduced by the Hubbard U term and provides greatly improved agreement with experiment for both cell volumes and heats of formation. The XDM dispersion correction is largely invariant to the chosen U value, making it a robust dispersion correction for actinide systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗