Search NASA⌕ Search

SEARCH · Search NASA

Results for “self-consistent field 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 19 records

Cartesian equivariant representations for learning and understanding molecular orbitals

Qualitative and quantitative orbital properties such as bonding/antibonding character, localization, and orbital energies are critical to how chemists understand reactivity, catalysis, and excited-state behavior. Despite this, representations of orbitals in deep learning models have been very underdeveloped relative to representations of molecular geometries and Hamiltonians. Here, we apply state-of-the-art equivariant deep learning architectures to the task of assigning global labels to orbitals, namely energies characterizations, given the molecular coefficients from Hartree–Fock or density functional theory. The architecture we have developed, the Cartesian Equivariant Orbital Network (CEONET), shows how molecular orbital coefficients are readily featurized as equivariant node features common to all graph-based machine-learned potentials. We find that CEONET performs well at predicting difficult quantitative labels such as the orbital energy and orbital entropy. Furthermore, we find that the CEONET representation provides an intuitive latent space for differentiating orbital character for the qualitative assignment of e.g. bonding or antibonding character. In addition to providing a useful representation for further integrating deep learning with electronic structure theory, we expect CEONET to be useful for automatizing and interpreting the results of advanced electronic structure methods such as complete active space self-consistent field theory. In particular, the ability of CEONET to infer multireference character via the orbital entropy paves the way toward the machine-learned selection of active spaces.

chemical reactions↗

Bridging the gap between molecules and materials in quantum chemistry with localized active spaces

The number of materials that “bridge the gap” between single molecules and extended solids, such as metal-organic frameworks and organic semiconductors, has been increasing. Consequently, there is a growing need for modeling approaches that effectively integrate the real-space molecular perspective employed by computational chemists and the reciprocal-space dispersive perspective employed by computational physicists. Here, we propose the localized active space (LAS) approach as a promising method to successfully bridge this gap. The LAS approach extends the active space concept from multiconfigurational methods such as complete active space self-consistent field theory to multiple molecular fragments via a product-form wave function ansatz. Here, we apply this method to solid state phenomena by treating each unit cell as a fragment with different sets of local quantum numbers (e.g., charge and excitation number). State interaction between these LAS states (LASSI) thus provides a comprehensive basis for the study of charge and energy transfer, meeting and surpassing the capabilities of single-reference fragmentation approaches such as constrained density functional theory (cDFT). Most centrally, we show how combining this LASSI approach with multiconfigurational pair-density functional theory (MC-PDFT) provides an elegant and efficient method to compute band structures that capture multiconfigurational character. We apply the LASSI band structure approach to the computation of band gaps in stretched hydrogen chain, polyacetylene, and bulk nickel oxide (NiO), finding good or excellent quantitative agreement with reference values in all cases. Additionally, we use the LAS basis in one-dimensional model systems to demonstrate its ability to treat difficult solid-state phenomena such as exciton transfer and excitation at p-n junctions.

method development↗

Radiation-induced segregation in dilute Fe–Cr: A rate-theory framework for the Cr enrichment–depletion transition at the grain boundary

Radiation-induced segregation (RIS) poses a significant challenge for ferritic Fe–Cr alloys under irradiation, as it can compromise mechanical integrity and increase susceptibility to intergranular corrosion. Yet, the mechanisms governing Cr segregation remain incompletely understood. Here, in this study, we present a physics-based rate-theory model parameterized using self-consistent mean field theory-based Onsager transport coefficients to investigate RIS at the grain boundary in dilute Fe-(0.1 at. %) Cr. Under equal production rates of vacancies and self-interstitial atoms (SIA), and their equal absorption rates by bulk dislocations, the model simulates the experimentally observed transition from Cr enrichment at low temperatures to depletion at higher temperatures. Under these unbiased conditions, systematic investigation reveals that while temperature-dependent transport properties dictate the segregation direction, dose rate, grain size, and dislocation density only influence the magnitude and spatial extent of Cr segregation. However, under more realistic conditions of preferential vacancy production within damage cascade and/or preferential SIA absorption by bulk dislocations, the enrichment-to-depletion transition shifts to lower temperatures. Our findings demonstrate that RIS predictions based solely on transport coefficients are valid only under symmetric point defect flux conditions, and that biases in defect production and absorption must be considered for accurate predictions. This work provides a mechanistic framework for understanding RIS in ferritic alloys and informs alloy design for advanced nuclear systems.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Substrate-Mediated Evaporation and Stochastic Evolution of Supported Au Nanoparticles

Here, we use in situ transmission electron microscopy with automated tracking to study supported gold nanoparticles (NPs) during high-temperature vacuum annealing. The average mass loss per NP is governed by a flat, nearly size-independent substrate-mediated evaporation profile. On top of this mean shrinkage, individual NPs show significant fluctuations in apparent growth or shrinkage, and NP volume follows a random-walk-like trajectory. To rationalize both the ensemble-mean behavior and the particle-resolved variability, we develop a self-consistent theory that couples substrate-mediated evaporation to collective 2D Ostwald-type mass exchange through a shared adatom field, described in terms of a renormalized screening length and background concentration. In the experimentally relevant regime, the theory predicts an approximately size-independent mean shrinkage rate and clarifies how net mass loss suppresses classical coarsening. Superimposed on this deterministic drift, we quantify stochastic volume trajectories and capture their fluctuation spectrum with a minimal Langevin description consistent with intermittent adatom attachment and detachment events. In addition, we characterize the lateral diffusive motion of NPs, which is responsible for their coalescence. Altogether, our results highlight that stochasticity is intrinsic at the nanoscale and that predicting the evolution of supported NPs at early and intermediate times requires a unified framework combining substrate-mediated evaporation, collective mass exchange, and stochastic fluctuations.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Structure of domain walls in chiral spin liquids

The chiral spin liquid is a canonical state of quantum spins combining topological and symmetry-breaking order, and possible experimental realizations have attracted growing interest. We examine the physics at interfaces between chiral spin liquid domains of opposite chirality. We show that a self-consistent mean-field description of spinons remains possible in the vicinity of a domain wall and use this to formulate a Ginzburg–Landau theory of the domain wall. The bulk of a chiral spin liquid contains gapped spinon excitations and gauge fluctuations, set by a finite spinon mass and a nonzero spinon Chern number. A third class of excitations consists of amplitude fluctuations of the spinon hoppings, which admit a geometric interpretation in terms of effective vielbein fields. These fluctuations are usually neglected because they are irrelevant for a homogeneous chiral spin liquid and are suppressed in standard large- N treatments. Going beyond the purely topological Chern–Simons limit, we incorporate these fluctuations into an effective field theoretic framework and show that they generate momentum-dependent corrections, including Chern–Simons-like linear terms and higher-order contributions, beyond the universal topological limit. We then analyze nontopological properties, including domain wall tension and edge velocity, and explain how they modify observables relative to the uniform case. These results connect measurable, nonuniversal quantities such as domain wall width, domain wall tension, and edge velocity to microscopic parameters and provide concrete targets for experiments.

domain walls↗

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established nonequilibrium Green's function (GF) approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal GFs, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz, and exact diagonalization for small systems using second-order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single-particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the GFs, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich nonequilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and band gap renormalization, which are largely washed out in self-consistent approaches yet are present in exact solutions. Furthermore, our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Electronic structure↗

Nonadiabatic dynamics of photoexcited thiopyridone isomers: An interplay between El-Sayed’s conditions and energy gap law

We investigated the non-adiabatic dynamics of photoexcited thiopyridone systems across their ortho-, meta-, and para-isomeric forms. The relaxation pathways of the three isomers in both gas phase and solvent environments are mapped using surface hopping dynamics based on time-dependent density functional theory. Our analysis highlights the influence of isomeric structures on photophysical behavior, offering insights into design principles to control photochemical phenomena. The simulations suggest a systematic reduction in the rate of intersystem crossing (ISC) from ortho- to meta- to para-isomer. Comparisons with multiconfigurational wave function methods in the gas phase further demonstrate how electronic structure influences the predicted dynamical pathways. The simulated dynamics demonstrates that the spin–orbit coupling strength alone does not determine the rate of ISC, as both state energetics and underlying electronic and structural features play decisive roles. These aspects explain the much slower ISC in the para-isomer, as well as the non-negligible role of the El-Sayed forbidden pathway in the computed ISC dynamics.

Complete-active space self-consistent field↗

Redox Potentials for Hydrogen Species in Acetonitrile and Tetrahydrofuran

The properties of the fundamental hydrogen species, proton (H+), hydrogen atom (H•), and hydride anion (H-) are critical to a vast range of chemical processes, yet their thermodynamic properties in nonaqueous solvents are not well established. A hybrid supermolecule-continuum approach is used to predict the Gibbs free energies of solvation (?G°solv) and standard redox potentials (E°) for the 2H+/H2 and H•/H- couples in acetonitrile (MeCN) and tetrahydrofuran (THF) following the approach previously used for water. Gas phase geometry optimizations were done using density functional theory (DFT) with the ?B97XD functional, Møller-Plesset second-order perturbation theory (MP2) with augmented correlation-consistent basis sets and composite-correlated molecular orbital theory (G3MP2 and G3MP2B3). Additional single point calculations in gas phase were performed for H- affinity for MeCN and THF at the coupled cluster CCSD(T) level using MP2 geometries. Solvation was included using the self-consistent reaction field SMD model. For H+, the ?G°solv values are predicted to be -252.2 kcal/mol in MeCN and -261.0 kcal/mol in THF. The redox potentials for the 2H+/H2 couple are +0.49 V in MeCN and +0.11 V in THF relative to the aqueous SHE, consistent with the available experimental data. The ?G°solv for H- is predicted to be -80 kcal/mol in MeCN and -68 in THF, and for H•, ?G°solv is predicted to be -1.3 kcal/mol in THF and -2.0 kcal/mol in MeCN. These solvation energies yield respective calculated redox potentials for the H•/H- couple of -0.14 V in MeCN and -0.68 V in THF relative to the aqueous SHE. H2 solvation is slightly exergonic in MeCN (-3.3 kcal/mol) and near thermoneutral in THF (+1.3 kcal/mol). The pKa(H2) is predicted to be 45.4 in MeCN and 47.8 in THF. The proton-coupled electron transfer (PCET) accounting term CG was evaluated giving values on the ferrocene scale within approximately 2 kcal/mol of literature in both solvents. These results provide benchmarks for hydrogen redox chemistry in MeCN and THF and improve consistency across E°, ?G°solv, pKa, and CG. These results provide benchmarks for hydrogen redox chemistry in nonaqueous media challenging previous estimates based on indirect extrapolations and assumptions.

Duda, Damian P.↗

A tale of two polarization paradoxes: The diamagnetic polarization paradox

An accurate calculation of the total polarization charge density in a plasma is essential for a self-consistent determination of the electric field. Yet, the “diamagnetic polarization paradox” refers to the fact that there is a paradoxical factor of 1/2 difference between the pressure-driven “diamagnetic polarization” density calculated using real space drift theory vs action-angle space guiding center and gyrokinetic theory that has not been explained before. Here, we show that both results can be made consistent with one another. Half of the diamagnetic polarization is due to the transformation from the guiding center density to the real space density. The other half is due to the fact that, within the drift kinetic ordering assumptions, the guiding center density should be expressed as the gyroaverage of the density in the limit of vanishing Larmor radius. A comprehensive review of polarization is presented to complete the derivation, and then we derive results that are required for the polarization calculations to agree. Expressions for the diamagnetic polarization density are given that are accurate to first order in amplitude and all orders in gyroradius within the gyrokinetic theory for a constant magnetic field. Applications to Maxwell–Boltzmann particle distribution functions (PDFs), including anisotropic temperature, are presented. Local invariants, like total energy and toroidal momentum, do not generate net polarization effects; the electric and thermodynamic polarizations must precisely cancel. In contrast, anisotropic dependence on the magnetic moment generates a net polarization proportional to the temperature anisotropy. Finally, when sources are present, the equilibrium PDF is approximately the ratio of two orbit averages.

adiabatic theory↗

A self-consistent electrostatic electrified shallow water model for charged liquid surface dynamics

The dynamics of an electrified liquid surface are investigated using a shallow water model that is self-consistently coupled with an electrostatic solver. To account for the spatiotemporal variation of a curved liquid surface, the electric field on curved surfaces is calculated by solving a two-dimensional electrostatic equation using the weighted least squares (WLSQ) interpolation method. First, the WLSQ implementation is verified with analytical theory obtained from a Laplace solution assuming a half-cylinder liquid surface profile. Then, the coupled electrified shallow water model is used to study the instability of liquid surface perturbation as a function of the potential drop between an electrode and conducting liquid, surface tension, and gravity. We present a linear dispersion theory of liquid surface instability for long-wavelength perturbations, including the effects of liquid viscosity. Furthermore, the effects of multiple sinusoidal surface perturbations on the electrified liquid surface instability are investigated.

Capillary waves↗

Instant folded strings, dark energy and a cyclic bouncing universe

We present a wholly self-consistent, complete cyclic bouncing cosmology based on components drawn from string theory and constructed in a way that is under perturbative control throughout (e.g., with temperature much less than the string scale and string coupling g s ≪ 1 at all times). The cyclic evolution is governed by standard dilaton-gravity in (3+1)-dimensions with a perturbatively generated potential and a coupling between the dilaton and a second field that becomes massless at ϕ = ϕ ESP , resulting in an enhanced symmetry point (ESP) that prevents the dilaton from running all the way to zero coupling. A central role is played by instant folded strings (IFSs) — fundamental strings with the unusual property of being much lighter than the string mass while extending far beyond the string length, and violating the Null Energy Condition (NEC). IFSs are produced classically when the string coupling grows with time, which occurs at two critical points in each cycle. In turn, they fulfill a dual function: Enabling cosmological bounces and initiating transient epochs of dark-energy domination that naturally transition into slow contraction. The resulting cosmology eliminates the cosmic singularity and multiverse problems of big bang inflationary models and robustly predicts time-varying IFS-induced dark energy and the absence of primordial B-mode polarization in the cosmic microwave background.

alternatives to inflation↗

Datasets for Custom-trained Machine-learning Interatomic Potentials: Nitric Acid Aqueous Solution

This dataset was generated using an iterative active learning strategy with the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials (MLIPs) for aqueous nitric acid. Each active-learning cycle consisted of three stages: (1) training, (2) exploration, and (3) labeling. The initial training set comprised approximately 800 randomly selected configurations from a previous study by Lewis et al. (https://doi.org/10.1021/jp205510q), which investigated nitric acid solutions at 2, 3, 4, and 5 mol/L. For all configurations, single-point calculations of atomic forces and total energies were performed at the quantum density functional theory BLYP-D2 and PBE-D3 levels of theory using the CP2K Quickstep module. Valence electrons were treated explicitly, while core electrons on all atoms were represented by norm-conserving Goedecker–Teter–Hutter (GTH) pseudopotentials. Long-range dispersion interactions were accounted for using Grimme dispersion corrections. Wave functions were expanded in a mixed Gaussian-and-plane-wave scheme using TZV2P-MOLOPT basis sets for all elements and an 800 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent field convergence was accelerated using orbital transformation and Direct Inversion in the Iterative Subspace, with a convergence threshold of 10^{-6}. All single-point calculations were carried out in periodic orthorhombic cells whose dimensions match those of the molecular configurations sampled from earlier trajectories. The CELL_REF keyword in CP2K was used to define a fixed reference cell, ensuring consistency in the reference data used for MLIP training, particularly when cell fluctuations are present in NpT simulations. The resulting high-fidelity energies and forces constitute the ground-truth labels used to train the MLIPs contained in this dataset.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗

Quantum critical collapse abhors a naked singularity

Classical critical collapse yields naked singularities from smooth initial data, challenging cosmic censorship, and shaping the spectrum of primordial black holes. We show that one-loop vacuum polarization near the threshold qualitatively changes this outcome by dressing the singularity with a horizon within a controlled semiclassical regime. In analytically tractable Einstein-scalar critical spacetimes, a one-loop $s$-wave treatment linearized around self-similar backgrounds shows that regularity uniquely selects an asymptotically Minkowskian, vacuum-polarization state. Its renormalized stress tensor carries a universal quantum growing mode that competes with the classical unstable mode, shifts the critical point, and generates a trapped surface along with a finite mass gap at the new threshold, thereby enforcing horizon formation even under arbitrary fine-tuning. In primordial collapse, the threshold shift enters exponentially into the formation fraction, while the mass gap truncates the low-mass tail, suggesting potentially important consequences for the predicted mass spectrum. Furthermore, these results provide a self-consistent semiclassical treatment of critical collapse and yield sharp predictions within the one-loop, near-critical, linearized regime.

Anomalies↗

Cyclotron breaking: a mechanism for parallel ion cyclotron waves to heat the fast solar wind

The Parker Solar Probe mission has observed near-continuous power in parallel ion cyclotron waves (PICWs) in the young, fast solar wind. These waves are unlikely to be directly produced by the turbulent cascade and are likely born of a local instability; yet, they are observed to both cool – and heat – the plasma. We propose that these observations can be self-consistently explained as the natural consequence of PICWs propagating in the inhomogeneous solar wind after they have been driven unstable. In this work, we argue that strong proton heating by a turbulent cascade of oblique ICWs will result in PICWs being driven unstable in a process known as quasi-linear focusing. Because the power in the turbulent cascade is concentrated at scales above the turbulent transition region, PICWs will be driven unstable within a range of wavenumbers parallel to the background magnetic field, 𝑘 ∥ , that is bounded from above by 𝑘$^{∗}_{∥P}$, corresponding to the start of the transition region. As unstable PICWs propagate away from the Sun to regions of lower proton density, their 𝑘 ∥ , multiplied by the proton inertial length 𝑑 p , increases. Eventually, 𝑘$^{∗}_{∥P}$ of the PICWs becomes larger than 𝑘$^{∗}_{∥P}$⁢𝑑 p and the waves damp, heating the solar wind. We call this effect ‘cyclotron breaking’, in analogy with ocean waves breaking on the shore. We then discuss the testable predictions of the theory, including a distinct heating signature in which PICWs cool fast protons and heat slow protons at any given heliocentric distance 𝑟. Finally, we conjecture that cyclotron breaking can lead to net heating by PICWs if the power emitted as PICWs decreases sufficiently rapidly with 𝑟 that local emission of PICWs is overwhelmed by the local damping of PICWs generated closer to the Sun.

plasma heating↗

Unveiling horizons in quantum critical collapse

Critical gravitational collapse offers a unique window into regimes of arbitrarily high curvature, culminating in a naked singularity arising from smooth initial data — thus providing a dynamical counterexample to weak cosmic censorship. Near the critical regime, quantum effects from the collapsing matter are expected to intervene before full quantum gravity resolves the singularity. Despite its fundamental significance, a self-consistent treatment has so far remained elusive. In this work, we perform a one-loop semiclassical analysis using the robust anomaly-based method in the canonical setup of Einstein gravity minimally coupled to a free, massless scalar field. Focusing on explicitly solvable near-critical solutions in both 2 + 1 and 3 + 1 dimensions, we analytically solve the semiclassical Einstein equations and obtain controlled, quantitative results for several long-standing questions within the dominant s-wave sector. We find that regularity uniquely selects a Boulware-like quantum state, encoding genuine vacuum polarization effects from the collapsing matter. Remarkably, the resulting quantum corrections manifest as a growing mode. Horizon-tracing analyses, incorporating both classical and quantum modes, reveal the emergence of a finite mass gap, signaling a phase transition from classical Type II to quantum-modified Type I behavior, thereby providing a quantum enforcement of the weak cosmic censorship. The most non-trivial aspect of our analysis involves dealing with non-conformal matter fields in explicitly time-dependent critical spacetimes. Along the way, we uncover intriguing and previously underexplored features of quantum field theory in curved spacetime.

2D Gravity↗

Atomistic modeling of lanthanide diffusion in refractory body-centered cubic molybdenum

Lanthanide fission products can strongly interact with candidate cladding alloys, but their transport properties in refractory metals remain poorly understood. Here, in this work, we investigate the atomic-scale diffusion behavior of La, Ce, Pr, and Nd in body-centered cubic (bcc) molybdenum, a potential candidate for advanced nuclear cladding. Self-consistent mean-field transport modeling is performed to evaluate the fission product transport and vacancy mobility, informed by first-principles and nudged elastic band calculations of vacancy formation energies, migration barriers, and solute–vacancy binding characteristics. Compared with bcc Fe, lanthanide solutes in bcc Mo exhibit slower tracer diffusion due to higher vacancy formation and migration energies. Furthermore, the calculations reveal that the influence of fission products on migration barriers in bcc Mo are not as extensive in range compared to bcc Fe. Among the studied lanthanides, La exhibits the strongest vacancy binding while also being the fastest diffuser in Mo. These findings highlight how refractory bcc alloys can reduce fission product infiltration, offering valuable insight into the development of durable cladding systems for advanced reactors.

36 - MATERIALS SCIENCE↗

Quantum solver for single-impurity Anderson models with particle-hole symmetry

Quantum embedding methods, such as dynamical mean-field theory (DMFT), provide a powerful framework for investigating strongly correlated materials. A central computational bottleneck in DMFT is in solving the Anderson impurity model (AIM), whose exact solution is classically intractable for large bath sizes. In this work, we benchmark a quantum-classical hybrid solver tailored for particle-hole symmetric AIMs, using the variational quantum eigensolver to prepare the ground state of the model with shallow quantum circuits. The solver uses shallow quantum ansätze and one set of variational parameters to prepare the ground state and its particle and hole excitations, enabling the construction of the impurity Green’s function through a continued-fraction expansion. We evaluate the performance of this approach across a few bath sizes and interaction strengths under noisy, shot-limited conditions. We compare three optimization routines (COBYLA, Adam, and L-BFGS-B) in terms of convergence and fidelity, assess the benefits of estimating a quantum-computed moment correction to the variational energies, and benchmark the approach by comparing the density of states computed from the impurity Green’s function against that obtained using a classical pipeline. Our results demonstrate the feasibility of Green’s function construction on near-term devices and establish practical benchmarks for quantum impurity solvers embedded within self-consistent DMFT loops.

Karabin, Mariia [ORNL]↗

Fractional charge and spin errors in self-consistent Green’s function theory

We examine fractional charge and spin errors in self-consistent Green’s function theory within a second-order approximation (GF2). For GF2, it is known that the summation of diagrams resulting from the self-consistent solution of the Dyson equation removes the divergences pathological to second-order Møller-Plesset (MP2) theory for strong correlations. In the language often used in density functional theory contexts, this means GF2 has a greatly reduced fractional spin error relative to MP2. The natural question then is what effect, if any, does the Dyson summation have on the fractional charge error in GF2? To this end, we generalize our previous implementation of GF2 to open-shell systems and analyze its fractional spin and charge errors. We find that like MP2, GF2 possesses only a very small fractional charge error, and consequently minimal many electron self-interaction error. This shows that GF2 improves on the critical failings of MP2, but without altering the positive features that make it desirable. Furthermore, we find that GF2 has both less fractional charge and fractional spin errors than typical hybrid density functionals as well as random phase approximation with exchange.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗