Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupled cluster methods”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

NASA Tech Briefs, July 2011

Topics covered include: 1) Collaborative Clustering for Sensor Networks; 2) Teleoperated Marsupial Mobile Sensor Platform Pair for Telepresence Insertion Into Challenging Structures; 3) Automated Verification of Spatial Resolution in Remotely Sensed Imagery; 4) Electrical Connector Mechanical Seating Sensor; 5) In Situ Aerosol Detector; 6) Multi-Parameter Aerosol Scattering Sensor; 7) MOSFET Switching Circuit Protects Shape Memory Alloy Actuators; 8) Optimized FPGA Implementation of Multi-Rate FIR Filters Through Thread Decomposition; 9) Circuit for Communication Over Power Lines; 10) High-Efficiency Ka-Band Waveguide Two-Way Asymmetric Power Combiner; 11) 10-100 Gbps Offload NIC for WAN, NLR, and Grid Computing; 12) Pulsed Laser System to Simulate Effects of Cosmic Rays in Semiconductor Devices; 13) Flight Planning in the Cloud; 14) MPS Editor; 15) Object-Oriented Multi Disciplinary Design, Analysis, and Optimization Tool; 16) Cryogenic-Compatible Winchester Connector Mount and Retaining System for Composite Tubes; 17) Development of Position-Sensitive Magnetic Calorimeters for X-Ray Astronomy; 18) Planar Rotary Piezoelectric Motor Using Ultrasonic Horns; 19) Self-Rupturing Hermetic Valve; 20) Explosive Bolt Dual-Initiated from One Side; 21) Dampers for Stationary Labyrinth Seals; 22) Two-Arm Flexible Thermal Strap; 23) Carbon Dioxide Removal via Passive Thermal Approaches; 24) Polymer Electrolyte-Based Ambient Temperature Oxygen Microsensors for Environmental Monitoring; 25) Pressure Shell Approach to Integrated Environmental Protection; 26) Image Quality Indicator for Infrared Inspections; 27) Micro-Slit Collimators for X-Ray/Gamma-Ray Imaging; 28) Scatterometer-Calibrated Stability Verification Method; 29) Test Port for Fiber-Optic-Coupled Laser Altimeter; 30) Phase Retrieval System for Assessing Diamond Turning and Optical Surface Defects; 31) Laser Oscillator Incorporating a Wedged Polarization Rotator and a Porro Prism as Cavity Mirror; 32) Generic, Extensible, Configurable Push-Pull Framework for Large-Scale Science Missions; 33) Dynamic Loads Generation for Multi-Point Vibration Excitation Problems; 34) Optimal Control via Self-Generated Stochasticity; 35) Space-Time Localization of Plasma Turbulence Using Multiple Spacecraft Radio Links; 36) Surface Contact Model for Comets and Asteroids; 37) Dust Mitigation Vehicle; 38) Optical Coating Performance for Heat Reflectors of the JWST-ISIM Electronic Component; 39) SpaceCube Demonstration Platform; 40) Aperture Mask for Unambiguous Parity Determination in Long Wavelength Imagers; 41) Spaceflight Ka-Band High-Rate Radiation-Hard Modulator; 42) Enabling Disabled Persons to Gain Access to Digital Media; 43) Cytometer on a Chip; 44) Principles, Techniques, and Applications of Tissue Microfluidics; and 45) Two-Stage Winch for Kites and Tethered Balloons or Blimps.

Source record↗
DOE OSTI2021

Convergent close-coupling calculations of electrons scattering on electronically excited molecular hydrogen

We use the adiabatic-nuclei molecular convergent close-coupling method to perform calculations of 0.01–1000 eV electrons scattering on the c 3 Π u , a 3 $Σ$$^{+}_{g}$, B 1 $Σ$$^{+}_{u}$, C 1 Π u , and EF 1 $Σ$$^{+}_{g}$ states of H 2 in the v=0 vibrational level. Elastic, superelastic, ionization, and grand-total cross sections are presented, as well as cross sections for excitation of the n=2–3 singlet and triplet states of H 2 (where n is the atomic-limit principle quantum number). Comparison with available theoretical results is made. Good agreement is found with the recent R-matrix results [J. Phys. B 53, 245203 (2020)] for most of the exchange and dipole-forbidden transitions, but not for the dipole-allowed transitions. Here, the sources of disagreement were found to be an unconverged partial-wave expansion and the utilization of the fixed-nuclei approximation (as opposed to adiabatic-nuclei) in the R-matrix calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

A multi-dimensional parametric study of variability in multi-phase flow dynamics during geologic CO 2 sequestration accelerated with machine learning

Successful geologic CO 2 storage projects depend on numerical simulations to predict reservoir performance during site selection, injection verification, and post-injection monitoring phases of the project. These numerical simulations solve non-linear sets of coupled partial differential equations, while accounting for multi-phase fluid dynamics on the basis of constitutive equations that are embedded into the solution scheme. As a consequence, individual simulations often require tens to hundreds of hours to complete on high-performance computing clusters. Moreover, laboratory experiments reveal that parametric functions for capillary pressure and relative permeability exhibit substantial variability, even within the same rock type. This combination of computational expense and wide-ranging parametric variability means that there remains substantial uncertainty in the behavior of multi-phase CO 2 -water systems, particularly in the context of feedbacks between relative permeability and capillary pressure. To bridge this knowledge gap, here we develop a novel workflow that utilizes physics-based numerical simulation to train an artificial neural network (ANN) emulator for interrogating the multivariate parameter space that governs both capillary pressure and relative permeability. With this approach, the ANN is trained to emulate both fluid pressure distribution and CO 2 saturation, which are then interrogated quantitatively to generate parametric response surface mappings with high-fidelity resolution. Results from this study initially show that capillary entry pressure is the dominant control on both CO 2 plume geometry and fluid pressure propagation when considering the combined effects of capillary pressure and relative permeability, particularly when phase interference is low and residual CO 2 saturation is high. Moreover, the ANN emulator provides tremendous computational speed-up by computing 2691 individual simulations in several minutes; whereas, the same simulation ensemble would have required ~3 years of simulation time using only physics-based simulation methods (25,000 times speed up).

58 GEOSCIENCES↗
DOE OSTI2026

A dendritic strontium river isoscape for fisheries applications in the Sacramento River basin, California, USA

Objective Understanding the origins and movements of fish is fundamental to effective conservation and fisheries management. Strontium isotope ratios ( 87 Sr/ 86 Sr) in otoliths provide a powerful tracer of natal origin and migratory pathways. However, existing 87 Sr/ 86 Sr isoscapes for the Sacramento River basin, an ecosystem that supports ecologically and economically important salmon populations, rely on discrete classification approaches that overlook unsampled habitats and do not incorporate spatial uncertainty. Our objective was to develop a continuous, network-explicit 87 Sr/ 86 Sr isoscape with quantified uncertainty to fill in data gaps and enable probabilistic assignments of fish origin and movement. Methods We used river water 87 Sr/ 86 Sr data from 106 sites (1997–2021) to develop spatial stream network models that use dendritic connectivity and watershed characteristics (lithology, bedrock age, and land cover) to predict river water 87 Sr/ 86 Sr throughout the basin. Models were fitted using maximum and restricted likelihood and were evaluated via Akaike’s information criterion and leave-one-out cross validation. We produced both historical (pre-dam) and present-day (below-dam) isoscapes, delineated uncertainty-informed isotopic ranges using k -means clustering, and applied a proof-of-concept Bayesian assignment to estimate natal origins and early rearing habitats for two endangered winter-run Chinook Salmon Oncorhynchus tshawytscha. Results Cross validation indicated strong performance of the 87 Sr/ 86 Sr model (leave-one-out cross validation: R 2 = 0.91; root mean square error = 0.0005). Uncertainty-informed clustering identified 19 isotopic “suites” (reaches with indistinguishable 87 Sr/ 86 Sr values) in present-day anadromous habitats and 25 suites in the historical network. Example natal and early rearing assignments included predictions that challenged expectations for juvenile salmon migration based on predicted river 87 Sr/ 86 Sr compositions. Conclusions This study developed a continuous, network-explicit 87 Sr/ 86 Sr isoscape that integrates existing river data to predict 87 Sr/ 86 Sr in unsampled reaches and the likely achievable range and resolution of otolith-based origin and life history inference. The resulting river isoscape provides a valuable tool to predict salmon movements and identify habitats supporting their survival and growth that otherwise might remain undetected. Coupling these predictions with complementary approaches that ground-truth juvenile presence (e.g., targeted fish surveys) represents an important step toward science-informed restoration and management of critical habitats throughout the Sacramento River basin.

Environmental sciences↗
DOE OSTI2025

Thermal Weight Determination and Interstate Coupling in State-Averaged ADAPT-VQE

Characterizing electronic thermal states at low temperatures is an important but challenging task in quantum chemistry and condensed matter physics, making it a prime candidate for a useful application in quantum computing. One of the most successful methods for state preparation on quantum computers is the Adaptive, Problem-Tailored (ADAPT) Variational Quantum Eigensolver (VQE), which has recently been generalized to treat excited states within a state-averaged framework as well as Gibbs states. In this work, we introduce Helmholtz-Optimized Thermal (HOT) ADAPT-VQE, an ancilla-free strategy for preparing Gibbs states that directly minimizes the Helmholtz free energy by targeting the dominant eigenstates of the thermal ensemble. We demonstrate the usefulness of HOT-ADAPT-VQE by predicting the free energy of two model systems with strongly correlated ground states: (1) the Fe 2+ cation in a magnetic field and (2) a [Cu 2 O 7 ] 10– fragment of the Mott insulator La 2 CuO 4 . Our results demonstrate that HOT-ADAPT-VQE significantly improves upon Gibbs-state estimates from multistate variants of ADAPT-VQE, often with substantially shallower quantum circuits, making it a promising candidate for thermal-state calculations.

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

Coupled Cluster Study of the Interactions of AnO 2 , AnO 2 + , and AnO 2 2+ (An = U, Np) with N 2 and CO

We report thermochemical and spectroscopic properties for actinyl complexes involving UO 2 2+/1+/0 and NpO 2 2+/1+/0 with N 2 and CO, together with the UO 2 -O 2 , UO 2 + -O 2 , and UO 2 + -NO complexes, have been studied for the first time using an accurate composite coupled cluster approach. Two general bonding motifs were investigated, end-on (η 1 ) and side-on (η 2 ) relative to the metal center of the actinyls. For end-on CO complexes, both C-coordinated (An-C) and O-coordinated (An-O) structures were examined, with the former always being lower in energy. All of the η 1 complexes were calculated to be stable, with dissociation energies ranging from 2 to 36 kcal/mol, except for that of UO 2 + -O 2 (the η 1 orientation for UO 2 + -NO was not amenable to single reference coupled cluster). In agreement with a previous study, the η 2 structure for UO 2 + -O 2 was calculated to be relatively strongly bound, by 22.3 kcal/mol in this work. The closely related NO complex, however, had a calculated dissociation energy of just 4.0 kcal/mol. The binding energy of O 2 to neutral UO 2 in a η 2 orientation was calculated to be very strong, 75.4 kcal/mol, and strongly resembled a UO 2 + (O 2 – ) complex at equilibrium. The N–N and C–O bonds were found to be somewhat activated for all the side-on (η 2 ) neutral An(IV) complexes, with stretching frequencies of N 2 or CO being red-shifted by as much as 480 cm –1 with a 0.06 Å bond length elongation. Dissociation energies for the η 1 complexes are strongly correlated with the extent of electron transfer from ligand to actinyl. The nature of bonding in the actinyl complexes is examined using natural resonance theory (NRT). The correlation between bonding motif and small molecule activation is in agreement with experiments in condensed phases.

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

NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts

In this work, we present NENCI-2021, a benchmark database of ~8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation–π and anion–π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion–π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (E int ) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The E int values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) E int value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent—these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)- based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation–π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the E int values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive E int values, and maximum errors of ±0.2–0.3 kcal/mol (i.e., ~±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.

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

Magnetic structure and properties of the compositionally complex perovskite (Y 0.2 La 0.2 Pr 0.2 Nd 0.2 Tb 0.2 )MnO 3

Large configurational disorder in compositionally complex ceramics can lead to unique functional properties that deviate from traditional rules of alloy mixing. In recent years, compositionally complex oxides (CCOs) have shown intriguing magnetic behavior including long-range order, enhanced magnetic exchange couplings, and mixed phase magnetic structures. This work focuses on how large local spin disorder affects magnetic ordering in a CCO. Specifically, we investigated the A-site alloyed perovskite, (Y 0.2 La 0.2 Pr 0.2 Nd 0.2 Tb 0.2 )MnO 3 , or (5A)MnO 3 , using a combination of bulk magnetometry, synchrotron X-ray diffraction, and temperature-dependent neutron diffraction. The five A-site ions have an average spin and ionic radius nearly equal to that of Nd 3+ ions, which minimizes structural distortions and allows for an understanding of the local spin disorder effects through a direct comparison with NdMnO 3 . Our magnetometry data show that (5A)MnO 3 exhibits two distinct phase transitions associated with the A-site and B-site sublattices, as seen in NdMnO 3 , as well as the presence of domain pinning and exchange bias at low temperature, suggesting a mixed phase magnetic ground state, as seen in other magnetic CCOs. Neutron powder diffraction shows clear long-range antiferromagnetic ordering below 67 K and refines to a Pn'ma' magnetic structure at low temperature, in excellent agreement with the well-studied behavior of NdMnO 3 . The two most notable differences in (5A)MnO 3 magnetism apparent from our data are a slight suppression of the B-site ordering temperature, which is explained by a smaller Mn–O–Mn bond angle in (5A)MnO 3 than NdMnO 3 , and the presence of a magnetic susceptibility transition above the B-site ordering, which could indicate the formation of a cluster glass but requires further study. Finally, this work demonstrates a general method of isolated investigation of size and spin disorder in CCOs and motivates future work using local structure probes to better understand the effects of nanoscale clustering and local spin disorder in magnetic CCOs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Trace Elemental Abundances in Calcium-Aluminum-Rich Inclusions in CV Chondrites

Introduction: Calcium-aluminum-rich inclusions (CAIs), are the first formed solids that define the age of the Solar System [1,2]. CAIs are thought to have condensed from nebular gas [3,4] within the first <1 Ma of Solar System formation [5,6]. CAIs have experienced numerous early Solar System processes including condensation, evaporation, melting, recrystallization, and aqueous alteration [e.g., 7]. The chemical, mineralogical, and textural diversity among CAIs results from a range of chemical and physical processes recorded during nebular and parent body epoch. This study aims to explore the mineralogical, textural, and chemical compositions of CAIs including the trace elemental abundances in CAI phases to determine the early Solar System processes recorded in them. Samples and Analytical Methods: We analyzed one CAI each from CV3 chondrites Northwest Africa (NWA) 5508 designated as ‘Saguaro’, and Northwest Africa (NWA) 12772 designated as ‘Hoopoe’. Back-scatter electron (BSE) images were collected using a Phenom XL scanning electron microscope (SEM) at the Lunar and Planetary Institute (LPI) and the JEOL JXA-8530F electron probe microanalyzer (EPMA) at Johnson Space Center (JSC)-NASA. Additionally, energy dispersive X-ray spectrometry (EDS) elemental maps of select areas for these samples were collected using a 15.0kV beam energy and a 40µA emission current. Using the EPMA, wavelength-dispersive X-ray spectroscopy (WDS) quantitative data were collected. In-situ trace element measurements for both CAIs were determined at JSC-NASA using a Photon Machines 193nm laser ablation system and a Thermo-Scientific Element-XR inductively coupled plasma mass spectrometer (ICP-MS). Analyses consisted of 30s ablations at 10Hz, spot sizes of 20-25µm, and a fluence of 6.0 J/cm2 for anorthite and melilite, and a 3.5 J/cm2 fluence for all other phases. NIST612 was used to correct for instrument drift, while BHVO-2g was used as a primary calibration standard. BCR-2g and in-house mineral standards were regularly measured as unknowns to ensure accuracy. Results: Saguaro is a coarse-grained CAI, ~11 x 6 mm in dimensions. Saguaro contains spinel, Al-rich pyroxene, anorthite, Mg-rich melilite, and minor perovskite in its interior and is therefore classified as a Type B CAI. Individual melilite grains shows normal compositional zoning with an Ak content ranging from ~24 to 54 with no apparent trend from the core to the edge of the CAI. The spinel appears euhedral and occurs both as clusters and as spinel palisades [8]. Two rim sequences surround most of the sample: the inner rim being a Wark-Lovering (WL) rim (~10-35 µm) containing pyroxene, spinel, and melilite (or anorthite), and the outer rim is a finer-grained, thicker (~100 µm), accretionary rim (Fig. 1). The mineral phases in Saguaro record an overall flat REE pattern with an average negative Eu anomaly in pyroxene, and an average positive Eu anomaly in anorthite and melilite respectively. Anorthite, melilite, and pyroxene have a minor depletion in Tm (Fig. 2). The Hoopoe CAI is a compact, coarse-grained ~6 × 4 mm in size. The major mineralogy includes hibonite, spinel, melilite, anorthite, and perovskite. Therefore, it is classified as a compact transitional type A and B. (?)zoning was observed in some hibonites. Individual melilite grains show both reverse and normal zoning, where the Ak content ranges from ~6- to 28. Melilite shows two distinct textures. One texture consisted of smooth melilite that appeared homogenous, while the second appeared to consist of many fine fractures. The spinel also often appears clustered. The WL-rim sequence surrounding Hoopoe is ~25 µm thick and composed of spinel, perovskite, hibonite, and melilite/anorthite. It is then partially surrounded by an outer accretionary rim (~75µm). Like before, refractory metal nuggets appeared concentrated near the WL rims. Other metal assemblages rich in Fe and Ni were also observed. All major mineral phases in Hoopoe display relatively flat REE patterns, except for varying Eu and Tm between phases (Fig. 2). There is a prominent negative Eu anomaly in perovskite and an average positive Eu anomaly in spinel, anorthite, and hibonite respectively (Fig. 2, 3). The mixed phases along the rim of the CAI also display a negative Eu anomaly, and all phases the CAI were depleted in Pb. Discussion: The CV3 CAIs analyzed in this study were classified based on their mineralogy and textures into Type A versus Type B CAIs [10]. Hibonite appears to be pseudomorphically replacing the spinel, (i.e., is hibonite in composition, but appears in the shape of spinel). Spinel palisades. The presence of spinel palisades present in Saguaro are consistent with the melting and recrystallization experienced by this CAI. Trace elemental analyses. Saguaro and Hoopoe display similar trace element patterns to each other, with both appearing generally flat, with anomalies in Eu, and Tm. Melilite and anorthite display positive Eu anomalies in both CAIs, in addition to the hibonite in Hoopoe (Fig. 2). The phases that are depleted in Eu are pyroxene and perovskite in both Saguaro and Hoopoe, respectively (Fig. 3). Given that Eu is volatile in reducing environments [11], this could possibly indicate reducing conditions at the time anorthite and melilite crystallized, with the gas they formed from containing Eu. As these CAIs continued to form, this gas as a result would become depleted in Eu. This also could be supported by the propensity of anorthite and melilite to take up Eu from its surroundings and incorporate it into their structure [12]. In addition, analyzing the assemblage of the phases in the Saguaro, melilite and anorthite (Eu enriched) often surround the pyroxene (Eu depleted) as they are crystallized. This intergrowth of phases and the proximity of the phases would support that the Eu is being incorporated into some phases, preventing it from incorporating into other. Trace elemental analyses of the CAI rims will be evaluated in more detail, as they are complicated by the transient signal being composed of a mixture of mineral phases. Broadly, however, the patterns in the rims of both CAIs are comparable to each other, and for Hoopoe, to the mixed phase patterns in the core (Fig. 3). Other studies have found that CAI rims can be depleted in Ce and Yb [13], however we did not observe these anomalies in the two CAIs discussed here. Given their similarity, the trace elemental analyses of the mixed interior (i.e. core) and rim phases could be interpreted as forming from similar, if not the same, reservoirs. The REE abundance between the rim and core of Hoopoe are also similar, indicating they may have formed from a gas of the same or similar composition. Acknowledgments: We thank the ASU Center for Meteorite Studies for loaning the samples used in this work and Tabb Prissel for his assistance with the analysis. Mouti Al-Hashimi thanks Sam Crossley and Cyrena Goodrich for their help with the LPI SEM training. This work was supported by the LPI Summer Intern Program in Planetary Science and the LPI Cooperative Agreement. References: [1] Connelly J.N. (2012) Science, 338, 651-655. [2] MacPherson G. J. (2014) Treatise on Geochem., 2, 139-179. [3] Grossman L. (1972) GCA, 36, 597-619. [4] Ebel, D.S. (2006) Meteorites and the Early Solar System II (D. S. Lauretta & H. Y. McSween, Eds.) 253-277. [5] MacPherson G. J. (2012) Earth Planet. Sci. Lett., 331-332, 43-54. [6] MacPherson G.J. (2017) GCA, 201, 65-82. [7] Krot A.N. (1995) Meteoritics & Planet. Sci., 30, 748-775. [8] Wark and Lovering (1982) GCA, 46, 2595-2607. [9] Palme H. and Jones A. (2003) Treatise on Geochemistry (H. D. Holland and K. K. Turekian Eds.), 1, 41-61. [10] Grossman L. (1980) Ann. Rev. Earth Planet. Sci., 8, 559-608. [11] Floss C. et al. (1996) GCA, 60, 1975-1997. [12] Mason B. and Martin P. M. (1974) Earth Planet. Sci. Lett., 22, 141-144. [13] Wark B. and Boynton W. V. (2001) Meteoritics & Planet. Sci., 36, 1135-1166.

X Mouti↗
NASA NTRS1996

Self-Consistent Simulation of the Brownian Stage of Dust Growth

It is a widely accepted view that in proto-planetary accretion disks the collision and following sticking of dust particles embedded in the gas eventually leads to the formation of planetesimals (coagulation). For the smallest dust grains, Brownian motion is assumed to be the dominant source of their relative velocities leading to collisions between these dust grains. As the dust grains grow they eventually couple to the turbulent motion of the gas which then drives the coagulation much more efficiently. Many numerical coagulation simulations have been carried out to calculate the fractal dimension of the aggregates, which determines the duration of the ineffective Brownian stage of growth. Predominantly on-lattice and off-lattice methods were used. However, both methods require simplification of the astrophysical conditions. The aggregates found by those methods had a fractal dimension of approximately 2 which is equivalent to a constant, mass-independent friction time. If this value were valid for the conditions in an accretion disk, this would mean that the coagulation process would finally 'freeze out' and the growth of a planetesimal would be impossible within the lifetime of an accretion disk. In order to investigate whether this fractal dimension is model independent, we simulate self-consistently the Brownian stage of the coagulation by an N-particle code. This method has the advantage that no further assumptions about homogeneity of the dust have to be made. In our model, the dust grains are considered as aggregates built up of spheres. The equation of motion of the dust grains is based on the probability density for the diffusive transport within the gas atmosphere. Because of the very low number density of the dust grains, only 2-body-collisions have to be considered. As the Brownian stage of growth is very inefficient, the system is to be simulated over long periods of time. In order to find close particle pairs of the system which are most likely to undergo a collision, we use a particle-in-cell (PIC) method for the early stages of the simulation where the system is still very homogeneous and a tree method later when the particles are more clustered.

Kempf, S.↗
DOE OSTI2024

Complex orders and chirality in the classical Kitaev-Γ model

It is well recognized that the low-energy physics of many Kitaev materials is governed by two dominant energy scales, the Ising-type Kitaev coupling 𝐾 and the symmetric off-diagonal Γ coupling. An understanding of the interplay between these two scales is therefore the natural starting point toward a quantitative description that includes subdominant perturbations that are inevitably present in real materials. This study focuses on the classical 𝐾−Γ model on the honeycomb lattice, with a specific emphasis on the region 𝐾< 0 and Γ > 0 , which is the most relevant for the available materials and which remains enigmatic in both quantum and classical limits, despite much effort. We employ large-scale Monte Carlo simulations on specially designed finite-size clusters and unravel the presence of a complex multisublattice magnetic order in a wide region of the phase diagram, whose structure is characterized in detail. We show that this order can be quantified in terms of a coarse-grained scalar-chirality order, featuring a counterrotating modulation on the two spin sublattices. Here, we also provide a comparison to previous studies and discuss the impact of quantum fluctuations on the phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Quantum simulations of excited states with active-space downfolded Hamiltonians

Many-body techniques based on the double unitary coupled cluster ansatz (DUCC) can be used to downfold electronic Hamiltonians into low-dimensional active spaces. It can be shown that the resulting dimensionality reduced Hamiltonians are amenable for quantum computing. Recent studies performed for several benchmark systems using quantum phase estimation (QPE) algorithms demonstrated that these formulations can recover a significant portion of ground-state dynamical correlation effects that stem from the electron excitations outside of the active space. These results have also been confirmed in studies of ground-state potential energy surfaces using quantum simulators. In this letter, we study the effectiveness of the DUCC formalism in describing excited states. We also emphasize the role of the QPE formalism and its stochastic nature in discovering/identifying excited states or excited-state processes in situations when the knowledge about the true configurational structure of a sought after excited state is limited or postulated (due to the specific physics driving excited-state processes of interest). In this context, we can view the QPE algorithm as an engine for verifying various hypotheses for excited-state processes and providing statistically meaningful results that correspond to the electronic state(s) with the largest overlap with a postulated configurational structure. We illustrate these ideas on examples of strongly cor- related molecular systems, characterized by small energy gaps and high density of quasi-degenerate states around the Fermi level.

quantum chemistry, quantum computing, computationa↗
DOE OSTI2022

Newton trust-region methods with primary variable switching for simulating high temperature multiphase porous media flow

Coupling multiphase flow with energy transport due to high temperature heat sources introduces significant new challenges since boiling and condensation processes can lead to dry-out conditions with subsequent re-wetting. The transition between two-phase and single-phase behavior can require changes to the primary dependent variables adding discontinuities as well as extending constitutive nonlinear relations to extreme physical conditions. Practical simulations of large-scale engineered domains lead to Jacobian systems with a very large number of unknowns that must be solved efficiently using iterative methods in parallel on high-performance computers. Performance assessment of potential nuclear repositories, carbon sequestration sites and geothermal reservoirs can require numerous Monte-Carlo simulations to explore uncertainty in material properties, boundary conditions, and failure scenarios. Due to the numerical challenges, standard NR iteration may not converge over the range of required simulations and require more sophisticated optimization method like trust-region. In this study, we use the open-source simulator PFLOTRAN for the important practical problem of the safety assessment of future nuclear waste repositories in the U.S. DOE geologic disposal safety assessment Framework. The simulator applies the PETSc parallel framework and a backward Euler, finite volume discretization. We demonstrate failure of the conventional NR method and the success of trust-region modifications to Newton’s method for a series of test problems of increasing complexity. Trust-region methods essentially modify the Newton step size and direction under some circumstances where the standard NR iteration can cause the solution to diverge or oscillate. Furthermore, we show how the Newton Trust-Region method can be adapted for Primary Variable Switching (PVS) when the multiphase state changes due to boiling or condensation. The simulations with high-temperature heat sources which led to extreme nonlinear processes with many state changes in the domain did not converge with NR, but they do complete successfully with the trust-region methods modified for PVS. This implementation effectively decreased weeks of simulation time needing manual adjustments to complete a simulation down to a day. Finally, we show the strong scalability of the methods on a single node and multiple nodes in an HPC cluster.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

How to detect lensing rotation

Gravitational lensing rotation of images is predicted to be negligible at linear order in density perturbations, but can be produced by the post-Born lens-lens coupling at second order. This rotation is somewhat enhanced for Cosmic Microwave Background (CMB) lensing due to the large source path length, but remains small and very challenging to detect directly by CMB lensing reconstruction alone. We show the rotation may be detectable at high significance as a cross-correlation signal between the curl reconstructed with Simons Observatory (SO) or CMB-S4 data, and a template constructed from quadratic combinations of large-scale structure (LSS) tracers. Equivalently, the lensing rotation-tracer-tracer bispectrum can also be detected, where LSS tracers considered include the CMB lensing convergence, galaxy density, and the Cosmic Infrared Background (CIB), or optimal combinations thereof. We forecast that an optimal combination of these tracers can probe post-Born rotation at the level of 5.7σ–6.1σ with SO and 13.6σ–14.7σ for CMB-S4, depending on whether standard quadratic estimators or maximum a posteriori iterative methods are deployed. We also show possible improvement up to 21.3σ using a CMB-S4 deep patch observation with polarization-only iterative lensing reconstruction. However, these cross-correlation signals have non-zero bias because the rotation template is quadratic in the tracers, and exists even if the lensing is rotation free. We estimate this bias analytically, and test it using simple null-hypothesis simulations to confirm that the bias remains subdominant to the rotation signal of interest. Detection and then measurement of the lensing rotation cross-spectrum is therefore a realistic target for future observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Benchmarking an Embedded Adaptive Sampling Configuration Interaction Method for Surface Reactions: H 2 Desorption from and CH 4 Dissociation on Cu(111)

Embedded (emb-) correlated wavefunction (CW) theory enables accurate assessments of both ground- and excited-state reaction mechanisms involved in heterogeneous catalysis. Embedded multireference second-order perturbation theory (emb-MRPT2) based on reference wavefunctions generated via embedded complete active space self-consistent field (emb-CASSCF) theory is currently state-of-the-art. However, the factorial scaling of CASSCF limits the size of active space and the complexity of systems that can be studied. In this work, we assess the efficacy of an alternative CW method, adaptive sampling configuration interaction (ASCI)–which enables large active spaces to be used–for studying surface reactions. We couple ASCI with density functional embedding theory (DFET) and benchmark its performance for two reactions: H 2 desorption from and CH 4 dissociation on the Cu(111) surface. Unlike embedded complete active space second-order perturbation theory (emb-CASPT2) that accurately reproduces a measured H 2 desorption barrier, embedded ASCI, using a very large active space (though one that still comprises a small portion of the full set of orbitals) fails to do so. Adding an extra correlation term from embedded Møller–Plesset second-order perturbation theory (emb-MP2) improves the desorption barrier and endothermicity predictions. Thus, the inaccuracy of embedded ASCI comes from the missing dynamic correlation from the many other electrons and orbitals not included in the active space. For CH 4 dissociation, again embedded ASCI overestimates the dissociation barrier compared to emb-CASPT2 predictions. Adding dynamic correlation from emb-MP2 helps correct the barrier. However, this composite approach suffers from double counting of correlation within embedded ASCI followed by emb-MP2 calculations. We therefore conclude that the state-of-the-art emb-MRPT2 based on reference wavefunctions generated via emb-CASSCF remains the method of choice for studying surface reactions. emb-ASCI is useful when large active spaces beyond the limit of emb-CASSCF are essential, such as to study complex surface reactions with significant multiconfigurational character (static correlation) but weak dynamic correlation.

08 HYDROGEN↗
DOE OSTI2025

Do We Know How to Model Reionization?

I compare the power spectra of the radiation fields from two recent sets of fully-coupled simulations that model cosmic reionization: “Cosmic Reionization On Computers” (CROC) and “Thesan”. While both simulations have similar power spectra of the radiation sources, the power spectra of the photoionization rate are significantly different at the same values of cosmic time or the same values of the mean neutral hydrogen fraction. However, the power spectra of the photoionization rate can be matched at large scales for the two simulations when the matching snapshots are allowed to vary independently. I.e., on large scales, the clustering of the radiation field in two simulations evolves similarly, but the exact timing of this evolution is different in different simulations and is not parameterized by an easily interpretable physical quantity like the mean neutral fraction or the mean free path. On small scales, large differences are present and remain partially unexplained. Both CROC and Thesan use the Variable Eddington Tensor approximation for modeling radiative transfer, but adopt different closure relations (optically thin OTVET versus M1). The role of this key difference is tested by using smaller simulations with a new cosmological simulation code that implements both closure relations in a controlled environment (the same hydro, cooling, and gravity solvers and the star formation recipe). In these controlled tests, both the M1 closure and the OTVET ansatz follow the expected behavior from a simple analytical approximation, demonstrating that the differences in the 2-point function of the radiation field induced by the choice of the Eddington tensor are not dominant.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Topology-Informed Design Rules for Deconstructable Thermoset Copolymer Networks

Existing models of thermoset deconstruction facilitated by incorporating cleavable comonomers rely on a mean-field reverse gel point paradigm, which predicts network dissolution once cleavable bonds reach a critical stoichiometric threshold, but does not account for where those bonds reside within the network architecture. Using reactive coarse-grained molecular dynamics simulations coupled with graph-theoretic analysis, we extend this stoichiometric picture to show that deconstructability is governed by the curing-imprinted network topology rather than stoichiometry alone. This topological organization is hierarchical: at the local scale, the elastic effectiveness of cross-link junctions determines which cross-links constitute the load-bearing scaffold; at the mesoscale, the cross-linking rate kinetically templates that scaffold into topologically modular communities─densely cross-linked clusters connected by sparse bridging strands that sustain network connectivity. Using betweenness centrality to identify nodes that disproportionately lie on intercommunity shortest paths, we demonstrate that effective deconstruction of the network into macromolecular fragments requires cleavable comonomers to intercept these high-centrality bridging strands. We further find that under uniform, disassortative comonomer incorporation, this topological requirement provides a mechanistic basis for extending the reverse gel point to incorporate network topology. We also show that modularity imposes a fundamental limit on fragment uniformity that persists even when the centrality requirement is met. Finally, we demonstrate that chain stiffness provides a nearly independent lever to suppress mechanically redundant cross-links and raise the glass transition temperature without significantly altering the deconstruction outcome. Together, these findings reframe the thermoset design space around network topology and provide actionable guidelines for engineering thermoset copolymers with predictable deconstructability and targeted thermomechanical performance.

coarse-grained molecular dynamics↗
DOE OSTI2025

Hydrogen Bond Benchmark: Focal‐Point Analysis and Assessment of DFT Functionals

We performed a hierarchical, convergent ab initio benchmark study and systematically analyzed the performance of density functional approximations for describing hydrogen bonds in small neutral, cationic, and anionic complexes, as well as in larger systems involving amide, urea, deltamide, and squaramide moieties. Focal point analyses (FPA), extrapolating to the ab initio limit, were carried out using correlated wave function methods up to CCSDT(Q) for the small complexes and CCSD(T) for the larger systems, together with correlation-consistent Gaussian basis sets up to the complete basis set limit. Optimized geometries and vibrational frequencies were obtained at the CCSD(T) level. The resulting FPA hydrogen-bond energies converge within a few tenths of a kcal mol −1 . These reference data were used to evaluate 60 density functionals (including 12 dispersion-corrected), spanning the local-density approximation (LDA), generalized gradient approximations (GGAs), meta-GGAs, hybrids, meta-hybrids, double-hybrids, and range-separated hybrids. Overall, the meta-hybrid M06-2X provides the best performance for both hydrogen bond energies and geometries, while the dispersion-corrected GGAs BLYP-D3(BJ) and BLYP-D4 also yield accurate hydrogen-bond data and can serve as cost-effective options for studying large and complex systems.

coupled cluster theory↗