Search NASA⌕ Search

SEARCH · Search NASA

Results for “low order model”

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 OSTI2022

Hydrostatic pressure effect on the Co-based honeycomb magnet BaCo 2 (AsO 4 ) 2

Here the honeycomb antiferromagnet BaCo 2 (AsO 4 ) 2 , in which small in-plane magnetic fields (H1 ≈ 0.26T and H2 ≈ 0.52T at T = 1.8 K < T N ≈ 5.4K) induce two magnetic phase transitions, has attracted attention as a possible candidate material for the realization of Kitaev physics based on the 3d element Co 2+ . Here, we report on the change in the transition temperature T N and the critical fields H1 and H2 of BaCo 2 (AsO 4 ) 2 with hydrostatic pressure up to ~20 kbar, as determined from magnetization and specific heat measurements. Within this pressure range, a marginal increase in the magnetic ordering temperature is observed. At the same time, the critical fields are changed significantly (up to ~25 –35%). Specifically, we find that H1 is increased with hydrostatic pressure, i.e., the antiferromagnetic state is stabilized, whereas H2, which was previously associated with a transition into a proposed Kitaev spin-liquid state, decreases with increasing pressure. We discuss to what extent these results are compatible with suggested models with sizable third-nearest-neighbor exchange. Overall, the results put constraints on the magnetic models that are used to describe the low-temperature magnetic properties of BaCo 2 (AsO 4 ) 2 .

2-dimensional systems↗
DOE OSTI2024

Anthropogenic extremely low volatility organics (ELVOCs) Govern the Growth of Molecular Clusters over the Southern Great Plains during the Springtime

New particle formation (NPF) and growth govern cloud condensation nuclei (CCN) concentrations in many regions. The mechanisms governing the nucleation of molecular clusters vary substantially in different regions of the atmosphere. Additionally, the growth of these clusters from ~2 to 20 nm sizes is often governed by the availability of extremely low volatility organic vapours (ELVOCs). While the pathways to ELVOC formation from the oxidation of biogenic monoterpenes with ozone is better understood, the chemical and mechanistic pathways for ELVOC formation from oxidation of anthropogenic organics are not well understood. We integrate measurements and three-dimensional regional model simulations with the Weather Research and Forecasting Model coupled to chemistry (WRF-Chem) to understand the processes governing new particle formation and growth and secondary organic aerosol (SOA) formation during the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystems (HI-SCALE) field campaign at the Southern Great Plains (SGP) observatory in Oklahoma, and contrast it with a site within the Bankhead National Forest (BNF), Alabama in Southeast USA, where 5-year long measurements will begin in 2024. Simulations show that nucleation rates are at least an order of magnitude higher at SGP compared to BNF during the springtime days (April 28 and May 14, 2016), largely due to lower H2SO4 concentrations at BNF, which are needed for nucleation. In addition, the larger CS at BNF (compared to SGP) increase the loss of molecular clusters by coagulation to pre-existing particles. Among the 8 different nucleation mechanisms in WRF-Chem, we find that the amine+H2SO4 nucleation mechanism dominates at the SGP site, while the pure organic ion induced nucleation mechanism dominates over BNF. Through various WRF-Chem sensitivity simulations, we find that anthropogenic ELVOCs are critical for explaining the growth of newly formed particles and the resulting number size distribution observed near the surface at the SGP site during the daytime. In addition, we show that treating organic particles as semisolid, with strong diffusion-limited uptake of organic vapours, brings model predictions into closer agreement with the observed evolution of particle size distribution. Simulations also predict that anthropogenic SOA, formed by the oxidation of aromatic volatile organic compounds (VOCs), is the dominant organic aerosol component at SGP, while biogenic SOA dominates particle composition at the BNF site in Southeast USA on these days.

Shrivastava, ManishKumar B.↗
DOE OSTI2025

How does ion temperature gradient turbulence depend on magnetic geometry? Insights from data and machine learning

Magnetic geometry has a significant effect on the level of turbulent transport in fusion plasmas. Here, we model and analyse this dependence using multiple machine learning methods and a dataset of >200 000 nonlinear gyrokinetic simulations of ion-temperature-gradient turbulence in diverse non-axisymmetric geometries. The dataset is generated using a large collection of both optimised and randomly generated stellarator equilibria. At fixed gradients and other input parameters, the turbulent heat flux varies between geometries by several orders of magnitude. Trends are apparent among the configurations with particularly high or particularly low heat flux. Regression and classification techniques from machine learning are then applied to extract patterns in the dataset. Due to a symmetry of the gyrokinetic equation, the heat flux and regressions thereof should be invariant to translations of the raw features in the parallel coordinate, similar to translation invariance in computer vision applications. Multiple regression models including convolutional neural networks (CNNs) and decision trees can achieve reasonable predictive power for the heat flux in held-out test configurations, with highest accuracy for the CNNs. Using Spearman correlation, sequential feature selection and Shapley values to measure feature importance, it is consistently found that the most important geometric lever on the heat flux is the flux surface compression in regions of bad curvature. The second most important geometric feature relates to the magnitude of geodesic curvature. These two features align remarkably with surrogates that have been proposed based on theory, while the methods here allow a natural extension to more features for increased accuracy. The dataset, released with this publication, may also be used to test other proposed surrogates, and we find that many previously published proxies do correlate well with both the heat flux and stability boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The MicroBooNE Single-Photon Low-Energy Excess Search (Public Note 1087)

MicroBooNE is a short baseline neutrino experiment at Fermilab designed to address the low energy excess observed by the MiniBooNE experiment. This note describes and presents preliminary results for the MicroBooNE analysis developed to address this excess as a single photon plus one or zero protons in the final state. The analysis assumes neutrino neutral current Δ resonance production followed by Δ radiative decay on argon (NC Δ → Nγ ) as the "signal model"; event reconstruction and selection have been developed and optimized in order to maximize efficiency and reduce cosmogenic and other beam-related backgrounds to the NC Δ → Nγ signal. We present the analysis methodology and validation checks performed on limited-statistics open data sets, corresponding to 5 x10 19 protons on target (POT), following a blind analysis, as well as the projected sensitivities for testing the Standard Model (SM) predicted rate for the NC Δ → Nγ process and for testing the interpretation of the previously observed MiniBooNE low energy excess as NC Δ → Nγ events, using the full anticipated MicroBooNE data set of 12.25 x 10 20 POT.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A higher-order finite-element implementation of the nonlinear Fokker–Planck collision operator for charged particle collisions in a low density plasma

Collisions between particles in a low density plasma are described by the Fokker–Planck collision operator. In applications, this nonlinear integro-differential operator is often approximated by linearised or ad-hoc model operators due to computational cost and complexity. In this work, we present an implementation of the nonlinear Fokker–Planck collision operator written in terms of Rosenbluth potentials in the Rosenbluth–MacDonald–Judd (RMJ) form. The Rosenbluth potentials may be obtained either by direct integration or by solving partial differential equations (PDEs) similar to Poisson's equation: we optimise for performance and scalability by using sparse matrices to solve the relevant PDEs. We represent the distribution function using a tensor-product continuous-Galerkin finite-element representation and we derive and describe the implementation of the weak form of the collision operator. We present tests demonstrating a successful implementation using an explicit time integrator and we comment on the speed and accuracy of the operator. Finally, we speculate on the potential for applications in the current and next generation of kinetic plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Identification of SPR Caverns with Multiple Oil Layers

This analysis shows that when lower density crude oil is injected into the top of an underground salt storage cavern containing more dense crude, separate oil phases can form and coexist indefinitely. This has been observed at the U.S. Strategic Petroleum Reserve in spite of geothermal heating and natural convection, which tend to mix the contents of containers with significant vertical extent subjected to wall and bottom heating. Such persistent layering can create operational challenges for meeting delivery specifications if high-value, low-vapor pressure oil becomes trapped below incoming low-density, high-vapor pressure oil, effectively blocking access to the lower layers until the top layer is removed. Previous conceptual models assumed that the oil injection process mixed incoming oil with resident oil in a storage cavern, forming a single oil phase with relatively homogeneous properties. Here, a review of historical data from the Strategic Petroleum Reserve revealed that several caverns contain multiple oil layers. As a result, oil layering needs to be another variable considered when planning oil movements at SPR in order to optimize low-vapor pressure oil availability to assist in oil delivery blending.

02 PETROLEUM↗
DOE OSTI2022

Anomalies in the topology of the temperature fluctuations in the cosmic microwave background: An analysis of the NPIPE and FFP10 data releases

We present a topological analysis of the temperature fluctuation maps from the Planck 2020 Data Release 4 NPIPE dataset and the Planck 2018 Data Release 3 FFP10 dataset. We performed a multiscale analysis in terms of the homology characteristics of the maps, invoking relative homology to account for the analysis in the presence of masks. We performed our analysis for a range of smoothing scales spanning sub- and super-horizon scales corresponding to a full width at half maximum (FWHM) of 5',10',20',40',80',160',320', and 640', and employed simulations based on the standard model for comparison, which assumes the initial fluctuation field to be an isotropic and homogeneous Gaussian random field. Examining the behavior of topological components, represented by the 0D homology group, we find the observations to be approximately 2σ or less deviant from the simulations for all resolutions and scales for the NPIPE dataset. For the FFP10 dataset, we detect a 2.96σ deviation between the observations and simulations at N = 128, FWHM = 80'. For the topological loops, represented by the first homology group, the simulations and observations are consistent within 2σ for most resolutions and scales for both the datasets. However, for the NPIPE dataset, we observe a high deviation between the observation and simulations in the number of loops at FWHM = 320', but at a low dimensionless threshold ν = –2.5. Under a Gaussian assumption, this would amount to a deviation of ~4σ. However, the distribution in this bin is manifestly non-Gaussian and does not obey Poisson statistics either. In the absence of a true theoretical understanding, we simply note that the significance is higher than what may be resolved by 600 simulations, yielding an empirical p-value of at most 0.0016. Specifically in this case, our tests indicate that the numbers arise from a statistically stable regime, despite being based on small numbers. For the FFP10 dataset, the differences are not as strong as for the NPIPE dataset, indicating a 2.77σ deviation at this resolution and threshold. The Euler characteristic, which is the alternating sum of the ranks of relative homology groups, reflects the deviations in the components and loops. To assess the significance of combined levels for a given scale, we employed the empirical and theoretical versions of the χ 2 test as well as the nonparametric Tukey depth test. Although all statistics exhibit a stable distribution, we favor the empirical version of the χ 2 test in the final interpretation, as it indicates the most conservative differences. For the NPIPE dataset, we find that the components and loops differ at more than 95%, but agree within the 99% confidence level with respect to the base model at N = 32, FWHM = 320'. The Euler characteristic at this resolution displays a per mil deviation. In contrast, the FFP10 dataset shows that the observations are consistent with the base model within the 95% confidence level, at this and smaller scales. This is consistent with the observations of the Planck analysis pipeline via Minkowski functionals. For the largest smoothing scale, N = 16, FWHM = 640', both datasets exhibit an anomalous behavior of the loops, where FFP10 data exhibit a deviation that is larger by an order of magnitude than that of the NPIPE dataset. In contrast, the values for the topological components and the Euler characteristic agree between observations and model to within a confidence level of 99%. However, for the largest scales, the statistics are based on low numbers and may have to be regarded with caution. Even though both datasets exhibit mild to significant discrepancies, they also exhibit contrasting behaviors at various instances. Therefore, we do not find it feasible to convincingly accept or reject the null hypothesis. Disregarding the large-scale anomalies that persist at similar scales in WMAP and Planck, observations of the cosmic microwave background are largely consistent with the standard cosmological model within 2σ.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Contributions of Solar + Storage to Future Power Needs

The costs of solar photovoltaics (PV) have been dropping in recent years, leading to increasing installations of solar PV systems and growing interest in how the technology will impact the electric grid. Additionally, as battery costs decline, there is increasing interest in understanding the benefits and limitations of battery storage for the grid as well as potential co-benefits of coupling battery storage and photovoltaics. This paper presents recent methodological developments to more accurately represent the value and limitations of coupled PV and battery systems in capacity expansion models. These are demonstrated using the Resource Planning Model (RPM), which co-optimizes capacity investments, transmission investments, and reduced-order dispatch through 2045. We simulate the evolution of the generation and transmission system in the Western Interconnection both with and without a coupled PV and battery technology option under two core scenarios - a baseline scenario and a high renewable penetration scenario - coupled with sensitivities assuming low and mid PV and battery cost projections. We analyze the degree to which a variety of configurations of coupled PV + Battery systems are able to meet future electricity needs, under which system conditions this coupled technology becomes beneficial to the electric grid and identify the displaced technologies when coupled solar and storage is introduced into the model. We find that capturing the benefits of co-locating these systems becomes increasingly important as the penetration of solar PV rises, and that in particular it is important to accurately capture the ability of these technologies to shift energy, provide firm capacity, and reduce expected curtailment.

14 SOLAR ENERGY↗
DOE OSTI2020

Calculation of mode lifetimes in weakly anharmonic solids using self-consistent ensemble eigenstates of the Liouvillian

Recent schemes for the efficient calculation of vibrational mode frequencies and lifetimes in weakly anharmonic systems such as solid crystal structures have shown promise, with simple relationships being observed between low order moments of the Liouvillian and the observed lifetime. While only shown empirically, two parameters seemed sufficient to determine the lifetime to high accuracy over a large range of temperatures for a few simple low degree of freedom models and a Lennard–Jones solid. In this paper we analytically reproduce the numerically observed relationship using a few generally applicable assumptions about the density of states. This advance is made by studying the autocorrelation based on an “ensemble eigenstate” of the Liouvillian. We propose the use of this ensemble eigenstate as the basis for autocorrelations to simplify calculations in richer systems. Finally, an analysis is also performed of the relationship between canonical and microcanonical autocorrelation functions and a long standing error in the long time behavior of such functions originally presented by Sen, Sinkovits, and Chakravarti is corrected.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Microscopic investigation of the 8 Li(n,γ)Li 9 reaction

Background: The 8 Li (n, γ) 9 Li reaction plays an important role in several astrophysics scenarios. It cannot be measured directly and indirect experiments have so far provided only cross section limits. Theoretical predictions differ by an order of magnitude. Purpose: In this work we study the properties of 9 Li bound states and low-lying resonances and calculate the 8 Li(n, γ) 9 Li cross section within the no-core shell model with continuum (NCSMC) with chiral nucleon-nucleon and three-nucleon interactions as input. Methods: Here, the NCSMC is an ab initio method applicable to light nuclei that provides a unified description of bound and scattering states well suited to calculate low-energy nuclear scattering and reactions. For the capture cross section calculation, we adjust calculated thresholds to experimental values. Results: Our calculations reproduce the experimentally known bound states as well as the lowest 5/2 – resonance of 9 Li. We predict a 3/2 – spin-parity assignment for the resonance observed at 5.38 MeV. In addition to the very narrow 7/2 – resonance corresponding presumably to the experimental 6.43 MeV state, we find several other broad low-lying resonances. Conclusions: Our calculated 8 Li(n, γ) 9 Li cross section is within the limits derived from the 1998 National Superconducting Cyclotron Laboratory Coulomb-dissociation experiment. However, it is higher than cross sections obtained in recent phenomenological studies. It is dominated by a direct E1 capture to the ground state with a resonant contribution at ≈ 0.2 MeV due to E2/M1 radiation enhanced by the 5/2 – resonance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Time-dependent-bases with local CUR decomposition method for accelerating turbulent combustion simulations

Here, this study presents a novel reduced-order modeling framework, Time-Dependent Bases with Local CUR decomposition (TDB-L-CUR), designed to efficiently and accurately approximate the species transport equations in reacting flow simulations. The method extends the existing TDB-CUR approach for chemically reacting flows (Jung et al. Comput. Methods Appl. Mech. Engrg. 437 (2025) 117758), which leverages matrix decomposition techniques to form a global-in-space, time-dependent low-dimensional manifold. While TDB-CUR performs well in homogeneous systems, it may be less well-suited to spatially heterogeneous systems such as turbulent flames, where higher-rank approximations are typically required. The proposed TDB-L-CUR framework introduces two methodological extensions to the baseline approach. First, it applies unsupervised clustering to partition the physical domain into distinct regions, enabling spatially localized manifold construction, thereby reducing the rank required for the reduced-order representation. Second, it incorporates a computational singular perturbation (CSP)-based scheme for identifying and penalizing fast species, allowing for spatio-temporally adaptive mitigation of chemical stiffness. The proposed framework is validated on a hierarchy of test cases, including a one-dimensional premixed flame, a two-dimensional nonpremixed ignition case with vortex interaction, and a three-dimensional turbulent premixed flame. TDB-L-CUR significantly improves accuracy over TDB-CUR while further reducing computational cost. The fully on-the-fly formulation of TDB-L-CUR (i.e., requiring no offline training or prior knowledge) makes it a robust and scalable tool for reduced-order modeling of reactive flows.

Local manifold↗
DOE OSTI2020

The MicroBooNE Single-Photon Low-Energy Excess Search

MicroBooNE is a short baseline neutrino experiment at Fermilab designed to address the low energy excess observed by the MiniBooNE experiment. This note describes and presents preliminary results for the MicroBooNE analysis developed to address this excess as a single photon plus one or zero protons in the final state. The analysis assumes neutrino neutral current Δ resonance production followed by Δ radiative decay on argon (NC Δ → $N_γ$) as the “signal model”; event reconstruction and selection have been developed and optimized in order to maximize efficiency and reduce cosmogenic and other beam-related backgrounds to the NC Δ → $N_γ$ signal. We present the analysis methodology and validation checks performed on limited-statistics open data sets, corresponding to 5×10 19 protons on target (POT), following a blind analysis, as well as the projected sensitivities for testing the Standard Model (SM) predicted rate for the NC Δ → $N_γ$ process and for testing the interpretation of the previously observed MiniBooNE low energy excess as NC Δ → $N_γ$ events, using the full anticipated MicroBooNE data set of 12.25×10 20 POT.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Effect of electron–electron interactions on the emission characteristics of ultracold, nanoscale photoemission-based electron sources

For nearly 60 years, cold field emitters have been the source of choice for electron microscopy due to their high brightness and relatively low energy spread. In this paper, we have examined an alternative: nanoscale electron sources based on near-threshold photoemission. While these sources have not yet been realized, they hold the potential to produce significantly brighter beams than cold field emitters. We model electron–electron Coulomb interactions in beams emitted from such sources to calculate the impact of these interactions on the brightness and the energy spread. Our results show that these sources can theoretically deliver more than an order of magnitude brighter electron beams compared to cold field emitters along with more than an order of magnitude smaller energy spread, before being limited by the Coulomb interactions. Electron sources with such high brightness and low energy spread would be transformational for electron microscopy, enabling electron energy loss-based vibrational spectroscopy at the sub-nanometer scale.

Gevorkyan, G S↗
DOE OSTI2024

Recalibration of missing low-frequency variability and trends in the North Atlantic Oscillation

Abstract Multi-decadal trends in the wintertime North Atlantic Oscillation (NAO) are under-represented by coupled general circulation models (CGCMs), consistent with a lack of autocorrelation in their NAO index series. This study proposes and tests two simple “reddening” approaches for correcting this problem in simulated indices based on simple one parameter short-term (AR; Auto-Regressive order 1) and long-term (FD; Fractional-Difference) time series filters. Using CGCMs from the Coupled Model Intercomparison Project Phase 6 (CMIP6), the FD filter successfully improves the autocorrelation structure of the NAO, and in turn the simulation of extreme trends, while the AR filter is less successful. The 1963–1993 NAO trend is the maximum 31-year trend in the historical period. Raw CGCMs underestimate the likelihood of this trend by a factor of ten but this discrepancy is corrected after reddening. CMIP6 future projections show that long-term (2024–2094) NAO ensemble mean trends systematically increase with the magnitude of radiative forcing: -2.4 to 3.5 hPa/century for low-to-high forcing after reddening (more than double the range using raw output). The related likelihood of future maximum 31year trends comparable to 1963–1993 ranges from 3 to 7% whereas none of these CMIP6 projections simulate this without reddening. Near-term projections of the next 31 years (2024–2054) are less sensitive than long term trends to the future scenario, showing weak-to-no forced trend. However, reddening increases the ensemble range by 74% (to +/-1 standard deviation/decade), which could increase/decrease regional climate change signals in the Northern Hemisphere by magnitudes that are underestimated when using raw CGCM output.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Ultra-low concentration terbium (Tb) adsorption on garlic peels biosorbent and its application for $\mathrm{Nd}$-$\mathrm{Fe}$-$\mathrm{B}$ scraps recovery

Motivated by the strategic value of middle-heavy rare earth elements (MHREEs), we proposed a bio-adsorption method to recover ultra-low concentration terbium (Tb(Ⅲ)) from industrial wastewater that originated from the Nd-Fe-B scraps recovery process. It could be found that Tb(Ⅲ) ions as low as 5.34 ng·mL -1 (ppb) could be adsorbed onto Ca(Ⅱ)-modified garlic peels (Ca-GP) within 10 min with the adsorption efficiency of 99.2% at pH 3.5. The adsorption behavior of Tb(Ⅲ) onto Ca-GP conformed the Langmuir isotherm model and pseudo-second-order kinetic model. Attenuated total reflection Fourier Transform Infrared Spectroscopy (ATR-FTIR), X-ray photoelectron spectroscopy (XPS) and density-functional theory (DFT) calculations results showed that Tb(Ⅲ) ions could be ion-exchanged with cations in Ca-GP such as Ca(Ⅱ) and -COOH. In application, column experiment results showed that the maximum bed adsorption capacity for Tb, praseodymium (Pr), neodymium (Nd), dysprosium (Dy), and Total rare earth elements (REEs) ions calculated by the Thomas model was 0.06, 3.50, 9.65, 7.34, and 27.37 µg·g -1 , respectively, with the initial concentration of 0.37 ng·mL -1 Tb, 20 ng·mL -1 Pr, 44 ng·mL -1 Nd, 67 ng·mL -1 Dy and 170 ng·mL -1 total REEs ions of actual solutions (below the solubility of rare-earth hydroxide). Here this manuscript thus provided a novel cost-effective and efficient approach to the enrichment and recovery of MHREEs from ultra-low concentration REEs ions-containing solutions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Confinement and String Breaking in the Compact Abelian Higgs Model

While real-time simulation of Quantum Chromodynamics remains technologically out of reach, simplified models for studying elements of QCD phenomenology abound. This work presents a simple model, a spin-1 truncation of the Compact Abelian Higgs Model simulated on qutrit sites, in which confinement and string breaking is accessible to current simulation methods. In the low-energy regime of 1+1D scalar electrodynamics, the heavy modes are integrated out, producing a spin chain effective Hamiltonian in which Gauss' law is implicitly satisfied. We study the spectrum of string-like excitations using DMRG methods on the order of 100 sites. We demonstrate that an added, local chemical potential, playing a role analogous to external charges, permits parameter-dependent measurements of physical features of interest like the string tension and effective meson mass. Varying the chemical potential also permits a characterization of string stability not assessed in prior studies of confining lattice models.

Senseman, Blake [Iowa U.]↗
DOE OSTI2020

The convergence of quadrupole rotational invariants from the nuclear shell model

Nuclei exhibit both single-particle and collective degrees of freedom, with the latter often subdivided into vibrational and rotational motions. Experimentally identifying the relative roles of these collective modes is extremely challenging, particularly in the face of possible shape coexistence. Model-independent, invariant quantities describing the deformation of a nucleus in the intrinsic frame have long been known but their determination potentially requires a large quantity of experimental data to achieve convergence. Through comparison with the nuclear shell model, the question of convergence is addressed. Shell-model calculations performed in the sd- and pf-shell model spaces are used to determine electric-quadrupole matrix elements for a multitude of low-lying states using the first 40 states of the relevant spins. Relative contributions to the rotationally invariant quantities from multiple states can therefore be determined. It is found that, on average, the inclusion of four intermediate states results in the leading-order invariant, $\langle\hat{Q^2}\rangle$, converging to within 10% of its true value and the triaxiality term, cos (3δ), converging to its true value, though some variance remains. Higher-order quantities relating to the softness of the nuclear shape are found to converge more slowly. The convergence of quadrupole rotationally invariant sum rules was quantified in the sd- and pf-shell model spaces and indicates the challenge inherent in a full determination of nuclear shape. The present study is limited to relatively small valence spaces. Finally, larger spaces, such as the rare-earth region, potentially offer faster convergence.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Ab initio generalized Langevin equation

We introduce a machine learning–based approach called ab initio generalized Langevin equation (AIGLE) to model the dynamics of slow collective variables (CVs) in materials and molecules. In this scheme, the parameters are learned from atomistic simulations based on ab initio quantum mechanical models. Force field, memory kernel, and noise generator are constructed in the context of the Mori–Zwanzig formalism, under the constraint of the fluctuation–dissipation theorem. Combined with deep potential molecular dynamics and electronic density functional theory, this approach opens the way to multiscale modeling in a variety of situations. Here, we demonstrate this capability with a study of two mesoscale processes in crystalline lead titanate, namely the field-driven dynamics of a planar ferroelectric domain wall, and the dynamics of an extensive lattice of coarse-grained electric dipoles. In the first case, AIGLE extends the reach of ab initio simulations to a regime of noise-driven motions not accessible to molecular dynamics. In the second case, AIGLE deals with an extensive set of CVs by adopting a local approximation for the memory kernel and retaining only short-range noise correlations. The scheme is computationally more efficient than molecular dynamics by several orders of magnitude and mimics the microscopic dynamics at low frequencies where it reproduces accurately the dominant far-infrared absorption frequency.

97 MATHEMATICS AND COMPUTING↗