Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Adiabatic process”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Adiabatic and radiative cooling are both important causes of aerosol activation in simulated fog events in Europe

Aerosol–fog interactions affect the visibility in, and life cycle of, fog and are difficult to represent in weather and climate models. Here we explore processes that impact the simulation of fog droplet number concentrations (N d ) at sub-kilometer scale horizontal grid resolutions in the UK Met Office Unified Model. We modify the parameterization of aerosol activation to include droplet activation by radiative cooling in addition to adiabatic cooling and determine the relative importance of the two cooling mechanisms. We further test the sensitivity of simulated N d to: (a) interception of droplets by trees and buildings, (b) overestimation of updrafts in temperature inversions (which leads to artificially high N d values), and (c) potential mechanisms for droplet deactivation due to downward fluctuations in supersaturation. We evaluate our model against observations from the ParisFog and LANFEX field campaigns, building on evaluation described in the companion paper. Including radiative cooling in the activation mechanism improves how accurately we represent the liquid water path and the vertical structure of the fog in our LANFEX case study. However, with radiative cooling, the N d are overestimated for most of the ParisFog cases and for the LANFEX case. The time-averaged overestimate exceeds a factor of three (the normalized mean bias factor exceeds 2.0) in 4 out of 11 ParisFog cases. Our sensitivity studies demonstrate how these overestimates can be mitigated. Assuming the overestimate affects both radiative and adiabatic cooling, we find that although radiative cooling is more often the dominant source, both cooling sources can sometimes dominate activation.

Ghosh, Pratapaditya [Carnegie Mellon University, P↗
DOE OSTI2025

Outer Radiation Belt Dynamics During the October 2012 Storm Revisited: Rapid Inward Radial Transport From a Dynamic Outer Boundary

Earth's outer radiation belt electron flux is highly variable and can be enhanced by over an order of magnitude over timescales less than one day, as observed during the October 2012 storm. Previous studies of this storm (e.g., Reeves et al., 2013, https://doi.org/10.1126/science.1237743) have invoked local acceleration to explain this. However, here, we argue that the observations can instead be explained by fast inward radial transport. One method often invoked to distinguish between these two acceleration processes is the existence of local peaks in electron phase space density (PSD) as a function of L* at fixed first, M, and second, K, adiabatic invariants. However, this method relies on the assumption that the evolution of the PSD as a function of L* occurs over timescales slower than the satellite orbital period. Here, high spatiotemporal resolution data from the Global Positioning System (GPS) spacecraft constellation is used to show that enhancements in the PSD occur during the October 2012 storm over short timescales not resolvable by the Van Allen Probes. In addition, Geostationary Operational Environmental Satellite spacecraft data also indicate that these enhancements are consistent with relativistic electron injections. A radial diffusion model is shown to reproduce the PSD dynamics observed by the Van Allen Probes, once rapid variations at the simulation outer boundary are included, consistent with GPS data. This verifies that apparently “locally growing” peaks in PSD along high apogee satellite orbits can be produced by fast inward radial transport without requiring the action of any local acceleration processes.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Scattering wave packets of hadrons in gauge theories: Preparation on a quantum computer

Quantum simulation holds promise of enabling a complete description of high-energy scattering processes rooted in gauge theories of the Standard Model. A first step in such simulations is preparation of interacting hadronic wave packets. To create the wave packets, one typically resorts to adiabatic evolution to bridge between wave packets in the free theory and those in the interacting theory, rendering the simulation resource intensive. In this work, we construct a wave-packet creation operator directly in the interacting theory to circumvent adiabatic evolution, taking advantage of resource-efficient schemes for ground-state preparation, such as variational quantum eigensolvers. By means of an ansatz for bound mesonic excitations in confining gauge theories, which is subsequently optimized using classical or quantum methods, we show that interacting mesonic wave packets can be created efficiently and accurately using digital quantum algorithms that we develop. Specifically, we obtain high-fidelity mesonic wave packets in the Z 2 and U(1) lattice gauge theories coupled to fermionic matter in 1+1 dimensions. Our method is applicable to both perturbative and non-perturbative regimes of couplings. The wave-packet creation circuit for the case of the Z 2 lattice gauge theory is built and implemented on the Quantinuum H1-1 trapped-ion quantum computer using 13 qubits and up to 308 entangling gates. The fidelities agree well with classical benchmark calculations after employing a simple symmetry-based noise-mitigation technique. This work serves as a step toward quantum computing scattering processes in quantum chromodynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Post-inflationary enhancement of adiabatic perturbations in modular cosmology

We show that multi-field inflationary models with negligible turning in field space during inflation can lead to an effective sourcing of adiabatic from entropic perturbations afterthe end of inflation. We illustrate this general phenomenon with a detailed analysis of an inflationary model whose scalar potential is determined by modular invariance. Its entropic perturbations are frozen during inflation, but instead, they are converted into adiabatic perturbations in the first post-inflationary e-folds. The curvature power spectrum, giving rise to CMB fluctuations, reaches a novel and enhanced plateau in this process; we address the implications for the inflationary observables As , n s and r.

cosmological perturbation theory↗
DOE OSTI2024

Development of Steady-State and Dynamic Mass and Energy Constrained Neural Networks for Distributed Chemical Systems Using Noisy Transient Data

The paper presents the development of algorithms for mass and energy constrained neural network models that can exactly conserve the overall mass and energy of distributed chemical process systems, even though the noisy transient data used for optimal model training violate the same. In contrast to approximately satisfying mass and energy balance constraints of a system by soft penalization of objective function, algorithms have been developed for solving equality-constrained nonlinear optimization problems, thus providing the guarantee of exactly satisfying the system mass and energy conservation laws. For developing dynamic mass-energy constrained network models for distributed systems, hybrid series and parallel dynamic-static neural networks have been leveraged. The developed algorithms for solving both the training and forward problems are validated using both steady-state and dynamic data in the presence of various noise characteristics. The developed data-driven algorithms are flexible to exactly satisfy mass and energy balance constraints for dynamic chemical processes if the system holdup information is available. The proposed network structures and algorithms are applied to the development of data-driven lumped and distributed models of an adiabatic superheater/reheater system, a nonisothermal continuous stirred tank reactor, as well as an electrically heated plug-flow reactor system where one form of energy gets transformed to another. It has been observed that the mass-energy constrained neural networks yield a root mean squared error of <1% with respect to the system truth for the case studies evaluated in this work.

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

Inverse kinematics study of the energy levels of 21 Ne populated with the 20 Ne +𝑑 reaction

In recent years there has been significant experimental effort aimed at studying the impact of 16 O neutron poisoning on the weak 𝑠 process in rotating massive stars. Improving the understanding of energy levels in 21 Ne is crucial to reducing uncertainties in the rates of the 𝛼-induced reactions on 17 O that determine the overall efficiency of the weak 𝑠 process. This paper reports on one such experiment: a study of the 20 Ne (𝑑,𝑝)⁢ 21 Ne reaction in inverse kinematics. Deduced spin-parity assignments were made based on adiabatic distorted-wave approximation analysis and neutron partial widths were estimated by comparison with a previous experiment. Of particular significance for nuclear astrophysics is the resulting neutron partial width for the 7820-keV energy level, estimated to be 12 200 ⁢(2900)⁢ eV, and an upper limit of ≤ 9300 eV for that of the 7749-keV energy level. These results are in tension with previous studies and therefore this paper also discusses the current state of the research into the astrophysically relevant energy levels and highlights both areas of agreement and areas of disagreement between this and various other studies that have investigated this nucleus.

direct reactions↗
DOE OSTI2024

Study on the magnetothermal properties of the D y 1 - x T b x A l 2 series of compounds

Here, in this work, we developed a theoretical model Hamiltonian, in the mean field approximation, to describe the magnetic and magnetocaloric behavior of the series of compounds Dy 1-x Tb x Al 2 (x = 0.00, 0.15, 0.25, 0.30, 0.35, 0.40 and 0.75). We adjusted the exchange parameters λ DyDy , λ TbTb , and λ DyTb to obtain the spin reorientation temperatures (T SR ) and the critical temperature (T C ) for each compound in the series. The results obtained by the Hamiltonian model agree satisfactorily with the experimental results. The heat capacity curves with and without an applied magnetic field, adiabatic temperature variation and isothermal entropy variation were modeled and compared with the experimental data. As the experimental results show, our model was also able to reproduce the change in the spin reorientation process: a first order spin reorientation transition appears for concentrations x = 0.15, 0.25, 0.30, 0.35, and no spin reorientation transitions after x = 0.40.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Time-resolved vacuum-ultraviolet photoelectron spectroscopy of the à 1 A u state of acetylene

Ultrafast time-resolved photoelectron spectra are reported for the vacuum-ultraviolet (VUV) photoionization of acetylene following excitation to the à 1 A u state via UV absorption at 200 nm. The excitation energy lies above the lowest dissociation threshold to C 2 H X̃ 2 Σ + + H, as well as above the threshold for adiabatic dissociation of the à 1 A u state to form C 2 H (à 2 Π) + H. The time-dependent mass spectra and photoelectron spectra provide insight into the intramolecular decay processes of the à 1 A u state. In addition, photoelectron spectra of the à 1 A u state with VUV light access both the X̃ 2 Π u and à 2 Σ g + states of the ion, as well as the predicted, but previously unobserved, 1 2 Π g state, which corresponds to a two-hole, one-particle configuration that lies in close proximity to the à 2 Σ g + state. The 1 2 Π g state is split into 2 A 2 + 2 B 2 and 2 A g + 2 B g states in the cis and trans configurations, respectively. In conclusion, electronic structure calculations, along with trajectory calculations, reproduce the principal features of the experimental data and confirm the assignment of the 1 2 Π g state.

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

Initial in-orbit operation of the soft X-ray spectrometer Resolve onboard the X-ray imaging and spectroscopy mission satellite

The X-Ray Imaging and Spectroscopy Mission satellite was launched on September 6, 2023 (UT). Its Resolve instrument is a high-resolution X-ray spectrometer enabled by a microcalorimeter array thermally anchored to a 50-mK heat sink. Many sensitive, critical sub-systems comprise Resolve, including a multistage cryogenic cooling system, thin-film aperture filters, low-noise electronics, on-board signal processing, and several sources of X-rays for calibration. We summarize the initial on-orbit power-on and checkout of Resolve that commenced immediately after launch. Soon after launch, the cryocoolers were activated, and their operation was successfully established. On October 9, 2023, the first cycle of the adiabatic demagnetization refrigerator was carried out, bringing the sensors to their steady-state operational temperatures. Following this, the energy resolution at 5.9 keV was successfully measured. The energy scale of the system is highly sensitive to the thermal environment surrounding both the sensors and their analog electronics. Gain correction was performed using reference X-ray lines from onboard calibration sources. To optimize cooler frequency settings, noise spectra were collected across a range of frequencies, and the most suitable frequency pair was selected based on the in-orbit environment. During the final phase of the checkout, an attempt was made to open the gate valve, which is designed to protect the Dewar’s interior from external pressure during ground operations and launch. Unfortunately, this attempt was unsuccessful. As a result, the checkout process was temporarily paused, and a stable operational strategy was subsequently developed to enable Resolve to function effectively with the gate valve remaining closed.

X-ray micro-calorimeter↗
DOE OSTI2025

Investigation of Isobaric Mixing as a Mechanism for Boundary‐Layer Cloud Formation

This study investigates the potential role of isobaric mixing in the formation of marine boundary layer clouds. Cloud formation theory emphasizes uplift and adiabatic cooling, but recent observations support the existence of small clouds forming at various altitudes, even below the lifting condensation level. Isobaric mixing of air with different thermodynamic properties can generate localized supersaturation. A Gaussian mixing model is employed to simulate this process, considering the correlation between temperature and water vapor. Cloud droplet size distributions from aircraft measurements show a persistent and prominent mode of small droplets at 9 m, and the size of this mode compares favorably with predictions from the model. The results suggests that isobaric mixing plausibly contributes to the formation of clouds, particularly those observed at multiple altitudes with narrow droplet size distributions. This finding highlights the importance of considering isobaric mixing processes in understanding and modeling cloud formation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Cryogenic Photoelectron Spectroscopic and Theoretical Study of the Electronic and Geometric Structures of Undercoordinated Osmium Chloride Anions OsCl n – ( n = 3–5)

A series of anionic transition metal halides OsCl n - (n=3–5) have been investigated using a newly developed, home-constructed, cryogenic anion cluster photoelectron spectroscopy. Here, the target anionic species are generated through collision-induced dissociation in a two-stage ion funnel. The measured vertical detachment energies (VDEs) are 3.48 eV, 4.54 eV, and 4.81 eV for n = 3, 4, and 5, respectively. Density functional theory calculations at the B3LYP-D3(BJ)//aug-cc-pVTZ (-pp) level predict the lowest energy structures of OsCl n - (n = 3–5) to be a quintet triangle, quartet square, and quintet square-based pyramid, respectively. The CCSD(T)-calculated VDEs and corresponding adiabatic detachment energies (ADEs) agree well with our experimental measurements. Analysis of the corresponding frontier molecular orbitals (FMOs) and charge density differences suggests that the d-orbitals of the transition metal Os play a primary role in the single-photon detachment processes, and the detached electrons originating from different molecular orbitals are distinguishable.

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

Efficient state preparation for the Schwinger model with a theta term

We present a comparison of different quantum state preparation algorithms and their overall efficiency for the Schwinger model with a theta term. While adiabatic state preparation is proved to be effective, in practice it leads to large gate counts to prepare the ground state. The quantum approximate optimization algorithm (QAOA) provides excellent results while keeping the counts small by design, at the cost of an expensive classical minimization process. We introduce a “blocked” modification of the Schwinger Hamiltonian to be used in the QAOA that further decreases the length of the algorithms as the size of the problem is increased. The rodeo algorithm (RA) provides a powerful tool to efficiently prepare any eigenstate of the Hamiltonian, as long as its overlap with the initial guess is large enough. We obtain the best results when combining the blocked QAOA ansatz and the RA, as this provides an excellent initial state with a relatively short algorithm without the need to perform any classical steps for large problem sizes. Published by the American Physical Society 2025

Bazavov, Alexei (ORCID:0000000321411901)↗
DOE OSTI2025

Assessment of metadynamic recrystallization in single copper particle impacts by focused ion beam tomography

We study single Cu-on-Cu impacts relevant to cold spray deposition and quantitatively analyze the metadynamic recrystallization (mDRX) that takes place after the impact by virtue of lingering impact adiabatic heat. Unlike prior studies, the current full 3D tomographic analysis of the mDRX volume shows that mDRX is extremely common in such impacts, although it is often missed when examining 2D sections. We also report an unexpected trend: there is a “sour spot” for mDRX at velocities about 20–40 % above the velocity for particle adhesion. This non-monotonic trend is contrary to the expectations based on increasing adiabatic heating with velocity. With a schematic model, we show that the trend can be explained on the basis of heat transfer: cooling of the heat-affected region is limited by transport through the bonded regions at the particle-substrate interface. Thus, bonding has a prominent role in the heat dissipation process and the best bonded particles most rapidly bulk quench, avoiding mDRX. Here, the developed semi-empirical model aligns with the experimental findings and may help inform microstructural evolution during cold spray and post-spray processing.

FIB-SEM tomography↗
DOE OSTI2025

Energetic and Electronic Properties of UX +/0/– for X = Li and Be and Comparison of the Properties of the Uranium Atom Binding to 2nd Row Elements Li–F

The bonding and spectroscopic properties of ULi +/0/– and UBe +/0/– to complete the series for UX +/0/– for X = Li to F were investigated by high-level ab initio SO-CASPT2 and CCSD(T) electronic structure calculations. The low-lying spin–orbit states were obtained at the SA-CASPT2/aQ-PP level; bond dissociation energies (BDEs), ionization energies (IEs), adiabatic electronic affinities (AEAs), and vertical detachment energies (VDEs) were calculated at the Feller-Peterson-Dixon (FPD) level. A dense manifold of low-lying states was predicted for ULi +/0/– and UBe +/0/– . Here, the calculated BDEs for ULi (37.7 kJ/mol) and UBe (8.0 kJ/mol) show that UBe is weakly bound. For redox processes, the BDEs increased for ULi + (109.3 kJ/mol), ULi – (47.4 kJ/mol), UBe + (35.6 kJ/mol), and UBe – (72.3 kJ/mol). The IE(ULi) = 4.650 eV is lower than IE(Li); the IE(UBe) = 5.901 eV is close to the IE(U) and to the IEs of UB, UC, UN, UO, and UF. The AEAs of ULi (0.708 eV) and UBe (0.989 eV) are lower than those for UB, UC, UN, and UO but higher than that for EA(UF). Natural bond orbital (NBO) calculations show that ULi has the 5f 3 6d 1 7s 2 configuration for U and 2s 1 for Li, with a small partial negative charge slightly delocalized on U. UBe arises from the U(5f 3 6d 1 7s 2 ) and Be(2s 2 ) electron configurations with no charge separation. The same calculations were made for WX (X = Li, Be, C–F) to enable detailed comparisons of the properties for UX with WX (X = Li–F). For WX, BDE(WX) is higher than that for UX for X = Li to N and lower than BDE(UX) for X = O and F, mostly due to the higher IE of W than U as ionic character becomes more important going from Li to F.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Asynchronicity in opposed-piston RCMs: Does it matter?

Rapid Compression Machines (RCMs) are widely utilized to study combustion phenomena at engine-relevant conditions, and significant efforts are typically made to create a quiescent environment, particularly for investigations of autoignition chemistry. Opposed-piston configurations can be advantageous due to shorter compression times and reduced surface area to volume ratios. Each side must be actuated simultaneously, but this can be challenging in practice. These devices, like most RCMs, utilize hydraulics for actuation, speed control and arrestation of the piston at the end of the stroke; there is no mechanical control or linkage of the two piston trajectories. To quantify the magnitudes and effects of piston asynchronous behavior, this work employs both detailed experimental measurements and, for the first time, high-fidelity, Direct Numerical Simulation (DNS). The boundary conditions are carefully considered applying insight from high-resolution linear variable differential transformer (LVDT) measurements of the piston trajectory and a zero-dimensional kinematics model of the piston-shaft assembly. Sufficient resolution in the piston crevice region is used. The complicated fluid dynamical behavior that can evolve during piston compression and the ensuing delay processes due to offset timings from t offset = 0-10 ms is elucidated. It is found that near t offset = 6 ms and beyond, the boundary layer on the face of the first-seating piston can be sufficiently perturbed, due initially to reemergence of gas from the crevice of the firstseating piston, so that the adiabatic core can become degraded at long ignition delay times. Substantial mixing of colder gas into the interior of the reaction chamber can alter the measurements, similar to effects previously observed for improper piston crevice configuration. In conclusion, experimental techniques to mitigate asynchronous behavior are discussed and demonstrated.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Spin-Controllable Dynamics in Defect-Engineered Carbon Nanotubes as Single Photon Emitters: Data-Driven Modeling and Computations

Quantum technologies, such as quantum computing and sensing, require efficient single-photon emission (SPE) sources that operate at room temperature in telecom wavelengths. While several materials can serve as SPE sources, no single platform meets all the criteria for efficiency, ambient operation, and scalability. Single-walled carbon nanotubes (SWCNTs) with covalently attached molecules offer a promising solution. Their SPE can be easily tuned via modifications of the SWCNT's diameter, chirality, and bonded molecules, enabling emission across near-IR to telecom wavelengths at ambient conditions. However, to fully realize the potential of SWCNTs and unlock their quantum capabilities, a deeper understanding of how structural defects from molecular adducts affect their emission and competing photoexcited processes is essential. To address this gap in our knowledge, this project combined quantum chemistry calculations with data-driven methods of cheminformatics (QSAR) and machine learning (ML). The developed computational approaches have provided several design strategies for covalent functionalization of SWCNTs to improve their optical response. The collaboration with Los Alamos National Lab (LANL) enabled direct comparison of computational and experimental data, facilitating method validation. This partnership was enhanced through access to LANL's Center for Integrated Nanotechnologies (CINT) utilizing User Facility Program and summer internships, which provided three NDSU graduate students with hands-on experience at LANL. The outcomes of this project included (1) Advancing the current stage of computational methods in accurate modeling of non-adiabatic spin-dependent photoexcited dynamics and its applicability to nanosystems consisting of thousands of atoms, realized as open-access codes linked to existing DFT-based software; (2) Establishing the relationship between the structure of adducts and SWCNTs and intrinsic excitonic and spin properties of defect states for guiding novel synthetic strategies and experimental probes of chemically functionalized SWCNTs as near-IR emitting materials; (3) Generating virtual libraries of hypothetical functionalized SWCNTs for virtual screening of their chemical structures and optical properties, leveraging new functionalities of SWCNTs; (4) Offering a unique experience for NDSU graduate students that prepared them for future scientific careers related to materials modeling and big data processing. These results were summarized in 12 published journal papers and 3 recently submitted papers. One of a key finding is that the position of defect sites on the SWCNT surface primarily drives the emission redshift (up to 100 meV), while the polarity of the defect-inducing molecules has a much smaller effect (~10 meV). However, the electron-donating or withdrawing properties of a molecule influence selecting reactivity of defect sites. These insights important for optimizing synthetic protocols for desired emissions in SWCNTs. We also revealed that the interaction between two defects at various positions on the SWCNT enhances the redshift and optical activity of states, favoring strong near-IR emission. This suggests that manipulations in defect concentrations is a promising strategy for controlling efficient emission. Mostly important, the defect position was found controllable by the spin states of photoexcited intermediates: Excited aromatic molecules form ortho defects with SWCNTs at their singlet states in the presence of oxygen, while oxygen-free conditions favor para defects via the triplet-state mechanism. Additionally, a heat-activated [2+2] cycloaddition reaction facilitates divalent defect formation with fewer bonding positions that narrows emission bands. These groundbreaking findings have been experimentally validated and significantly advance our understanding of defect chemistry in SWCNTs. Using a novel encoding technique and 3D-MoRSE descriptors, we developed highly accurate ML/QSAR models to predict both the 3D structure and optical properties of SWCNTs with chemical defects. This model enabled the creation of a virtual library of 125,556 structures, providing new insights into the relationship between SWCNT-defect structure and emission.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

High sensitivity of simulated fog properties to parameterized aerosol activation in case studies from ParisFog

Aerosols influence fog properties such as visibility and lifetime by affecting fog droplet number concentrations (N d ). Numerical weather prediction (NWP) models often represent aerosol–fog interactions using highly simplified approaches. Incorporating prognostic size-resolved aerosol microphysics from climate models could allow them to simulate N d and aerosol–fog interactions without incurring excessive computational expense. However, microphysics code designed for coarse spatial resolution may struggle with sub-kilometer-scale grid spacings. Here, we test the ability of the UK Met Office Unified Model to simulate aerosol and fog properties during case studies from the ParisFog field campaign in 2011. We examine the sensitivity of fog properties to variations in N d caused by modifications to simulated aerosol activation. Our model, with a 500 m horizontal resolution and interactive aerosol and cloud microphysics, significantly underpredicts N d , although it only slightly underestimates the cloud condensation nuclei concentration. With an updated version of the Abdul-Razzak and Ghan (2000) activation scheme, we produce N d that are more consistent with those predicted by a cloud parcel model under fog-like conditions. We activate droplets only by adiabatic cooling. We incorporate more realistic hygroscopicities for sulfate and organic aerosols and explore the sensitivity of simulated N d to unresolved updrafts. We find that both N d and simulated fog liquid water content are very sensitive to the updated activation scheme but remain less affected by the update to hygroscopicities. Our improvements offer insights into the physical processes regulating N d in stable conditions, potentially laying foundations for improved operational fog forecasts that incorporate interactive aerosol simulations or aerosol climatologies.

Ghosh, Pratapaditya [Carnegie Mellon University, P↗
DOE OSTI2024

ELECTRONIC STRUCTURE METHODS AND PROTOCOLS WITH APPLICATION TO DYNAMICS, KINETICS AND THERMOCHEMISTRY

Hydrocarbon combustion involves the reaction dynamics of a tremendous number of species beginning with many-component fuel mixtures and proceeding via a complex system of intermediates to form primary and secondary products. Combustion conditions corresponding to new advanced engines and/or alternative fuels rely increasingly on autoignition and low-temperature-combustion chemistry. In these regimes various transient radical species such as HO2, ROO·, ·QOOH, HCO, NO2, HOCO, and Criegee intermediates play important roles in determining the detailed as well as more general dynamics. A clear understanding and accurate representation of these processes is needed for effective modeling. Given the difficulties associated with making reliable experimental measurements of these systems, computation can play an important role in developing these energy technologies. Accurate calculations have their own challenges since even within the simplest dynamical approximations such as transition state theory, the rates depend exponentially on critical barrier heights and these may be sensitive to the level of quantum chemistry. Moreover, it is well-known that in many cases it is necessary to go beyond statistical theories and consider the dynamics. Quantum tunneling, resonances, radiative transitions, and non-adiabatic effects governed by spin-orbit or derivative coupling can be determining factors in those dynamics. Building upon progress made during a period of prior support through the DOE Early Career Program, this project combines developments in the areas of potential energy surface (PES) fitting and multistate multireference quantum chemistry to allow spectroscopically and dynamically/kinetically accurate investigations of key molecular systems (such as those mentioned above), many of which are radicals with strong multireference character and have the possibility of multiple electronic states contributing to the observed dynamics. An ongoing area of investigation is to develop general strategies for robustly convergent electronic structure theory for global multichannel reactive surfaces including diabatization of energy and other relevant surfaces such as dipole transition. Combining advances in ab initio methods with automated interpolative PES fitting allows the construction of high-quality PESs (incorporating thousands of high-level data) to be done rapidly through parallel processing on high-performance computing (HPC) clusters. In addition, new methods and approaches to electronic structure theory will be developed and tested through applications. This project will explore limitations in traditional multireference calculations (e.g., MRCI) such as those imposed by internal contraction, lack of high-order correlation treatment and poor scaling. Methods such as DMRG-based extended active-space CASSCF and various Quantum Monte Carlo (QMC) methods will be applied (including VMC/DMC and FCIQMC). Insight into the relative significance of different orbital spaces and the robustness of application of these approaches on leadership class computing architectures will be gained. Synergy with other components of this research program such as automated PES fitting and multireference quantum chemistry will be used to address challenges encountered by the standard approaches to computational thermochemistry (those being single-reference quantum chemistry and perturbative treatments of the anharmonic vibrational energy, which break down for some cases of electronic structure or floppy strongly coupled vibrational modes).

74 ATOMIC AND MOLECULAR PHYSICS↗