SEARCH · Engineering Papers

Results for “statistical analytics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Scaling Theory of Neutral Sequence-Specific Polyampholytes

The primary sequence of charged monomers in a polyampholyte has a profound effect on its conformational properties. In this work, theory and simulations are used to predict how the sequence influences the single-chain behavior of globally neutral polyampholytes (PAs) under salt-free conditions. We consider PAs with Markov statistics of charges, where charge “blockiness” is defined by the correlation λ along the chain. These PAs cover a wide spectrum of primary sequences ranging from alternating (λ = –1) to ideally random (λ = 0) to diblock (λ → 1). In a theta solvent, sufficiently long PAs of any primary sequence form globules, but their internal structure, as well as the position and sharpness of the coil–globule transition, crucially depends on λ. Further, using scaling arguments, we demonstrate that the excess charge of the globule correlation volume (blob) increases with increasing λ, thereby inducing stronger Coulomb attractions between neighboring blobs of opposite charge and leading to higher density PA globules. Scaling predictions are supported by the results of the random phase approximation (RPA) and coarse-grained molecular dynamics simulations. A generalization to good and poor solvents is presented that enables the construction of the resulting scaling diagram of PA conformations in terms of solvent quality and charge blockiness, thereby providing a comprehensive molecular understanding of the role of the PA primary sequence in the chain structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Near-canopy horizontal concentration heterogeneity of semivolatile oxygenated organic compounds and implications for 2-methyltetrols primary emissions

Semivolatile oxygenated organic compounds (SV-OVOCs) are important atmospheric species, in particular for the production and chemistry of atmospheric particulate matter and related impacts on air quality and climate. In this study, SV-OVOCs were collected in the horizontal plane of the roughness layer over the tropical forest in the central Amazon during the wet season of 2018. A sampler mounted to a copter-type, hovering unmanned aerial vehicle was used. Underlying the collection region, a plateau forest transitioned into a slope forest across several hundred meters. The concentrations of pinonic and pinic acids, which are monoterpene oxidation products, had no statistical difference over the two forests. By comparison, across the study period, differences in the concentration of 2-methyltetrols, which are products of isoprene oxidation, ranged from -70% to +480% over the two forests. The chemical lifetime of 2-methyltetrols in the atmosphere is sufficiently long that heterogeneity in the isoprene emission rate from the two forests followed by atmospheric oxidation does not explain the concentration heterogeneity of 2-methyltetrols. Standing waves and local meteorology also cannot account for the heterogeneity. Of the possibilities considered, the most plausible explanation is the direct emission from the forest of 2-methyltetrols produced through biological processes within the plants. Under this explanation, the rate of direct SV-OVOC emissions should be modulated by forest type and related environmental stressors. Direct emissions of SV-OVOCs should be more broadly considered for constraining and improving models of atmospheric composition, transport, and chemistry over tropical forests.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

A crossed molecular beams and computational study on the unusual reactivity of banana bonds of cyclopropane (c-C 3 H 6 ; X 1 A$^{'}_{1}$ ) through insertion by ground state carbon atoms (C( 3 P j ))

The mechanism and chemical dynamics of the reaction of ground electronic state atomic carbon C( 3 P j ) with cyclopropane c-C 3 H 6 (X 1 A$^{'}_{1}$) have been explored by combining crossed molecular beams experiments with electronic structure calculations of the pertinent triplet C 4 H 6 potential energy surface and statistical computations of product branching ratios under single-collision conditions. The experimental findings suggest that the reaction proceeds via indirect scattering dynamics through triplet C 4 H 6 reaction intermediate(s) leading to C 4 H 5 product(s) plus atomic hydrogen via a tight exit transition state, with the overall reaction exoergicity evaluated as 231 ± 52 kJ mol –1 . The calculations indicate that C( 3 P j ) can easily insert into one of the three equivalent C–C ‘banana’ bonds of cyclopropane overcoming a low barrier of only 2 kJ mol –1 following the formation of a van der Waals reactant complex stabilized by 15 kJ mol –1 . The carbon atom insertion into one of the six C–H bonds is also feasible via a slightly higher barrier of 5 kJ mol –1 . Here, these results highlight an unusual reactivity of cyclopropane's banana C–C bonds, which behave more like unsaturated C–C bonds with a π-character than saturated σ C–C bonds, which are known to be generally unreactive toward the ground electronic state atomic carbon such as in ethane (C 2 H 6 ). The statistical theory predicts the overall product branching ratios at the experimental collision energy as 50% for 1-butyn-4-yl, 33% for 1,3-butadien-2-yl, i-C 4 H 5 , and 11% for 1,3-butadien-1-yl, n-C 4 H 5 , with i-C 4 H 5 (230 kJ mol –1 below the reactants) favored by the C–C insertion providing the best match with the experimentally observed reaction exoergicity. The C( 3 P j ) + c-C 3 H 6 reaction is predicted to be a source of C 4 H 5 radicals under the conditions where its low entrance barriers can be overcome, such as in planetary atmospheres or in circumstellar envelopes but not in cold molecular clouds. Both i- and n-C 4 H 5 can further react with acetylene eventually producing the first aromatic ring and hence, the reaction of the atomic carbon with c-C 3 H 6 can be considered as an initial step toward the formation of benzene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Entanglement phases in large- N hybrid Brownian circuits with long-range couplings

Here, we develop solvable models of large-N hybrid quantum circuits on qubits and fermions with tunable long-range power-law interactions and continuous local monitoring. These models provide analytical access to the entanglement phase diagram and error-correcting properties of many-body entangled nonequilibrium states generated by such dynamics. In one dimension, the long-range couplings are irrelevant for α > 3/2, where α is the power-law exponent, and the models exhibit a conventional measurement-induced phase transition between volume- and area-law entangled phases. For 1/2 < α < 3/2, the long-range couplings become relevant, leading to a nontrivial dynamical exponent at the measurement-induced phase transition. More interestingly, for α < 1, the entanglement pattern receives a subvolume correction for both area-law and volume-law phases, indicating that the phase realizes a quantum error correcting code whose code distance scales as L 2–2α . While the entanglement phase diagram is the same for both the interacting qubit and fermionic hybrid Brownian circuits, we find that long-range free-fermionic circuits exhibit a distinct phase diagram with two different fractal entangled phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Compressive Response and Energy Absorption of Additively Manufactured Elastomers with Varied Simple Cubic Architectures

Additive manufacturing, and particularly the vat photopolymerization process, enables the fabrication of complex geometries at high resolution and small length scales, making it well-suited for fabricating cellular structures (e.g., foams and lattices). Among these, elastomeric cellular structures are of growing interest due to their tunable compliance and energy dissipation. However, comprehensive data on the compressive behavior of these structures remains limited, especially for investigating the structure-property effects from changing the density and distribution of material within the cellular structure. This study explores how the mechanical response of polyurethane-based simple cubic structures changes when varying volume fraction, unit cell length, and unit cell patterning, which have not been systematically investigated previously in additively manufactured elastomers. Increasing volume fraction from 10% to 50% yielded significant changes in compressive stress–strain performance (decreasing strain at 0.5 MPa by 41.6% and increasing energy absorption density by 3962.5%). Although changing the unit cell length between 2.5 and 7 mm in ~30 mm parts did not result in statistically different stress–strain responses, modifying the configuration of struts of different thicknesses across designs with 30% volume fraction altered the stress–strain behavior (differences of 12.5% in strain at 0.5 MPa and 109.4% for energy absorption density). Power law relationships were developed to understand the interactions between volume fraction, unit cell length, and elastic modulus, and experimental data showed strong fits (R 2 > 0.91). These findings enhance the understanding of how multiple structural design aspects influence the performance of elastomeric cellular materials, providing a foundation for informing strategic design of tailorable materials for diverse mechanical applications.

36 MATERIALS SCIENCE↗

DOE OSTI2023

All-Polyhydroxyalkanoate Triblock Copolymers via a Stereoselective-Chemocatalytic Route

Biodegradable polyhydroxyalkanoate (PHA) homopolymers and statistical copolymers are ubiquitous in microbially produced PHAs, but the step-growth polycondensation mechanism the biosynthesis operates on presents a challenge to access well-defined block copolymers (BCPs), especially higher-order tri-BCP PHAs. We report a stereoselective-chemocatalytic route to produce discrete hard–soft–hard ABA all-PHA tri-BCPs based on the living chain-growth ring-opening polymerization of racemic (rac) 8-membered diolides (rac-8DL R ; R denotes the two substituents on the ring). Depending on the composition of the soft B block, originated from rac-8DL R (R = Et, n Bu), and its ratio to the semicrystalline, high-melting hard A block, derived from rac-8DL Me , the resulting all-PHA tri-BCPs with high molar mass (M n up to 238 kg mol –1 ) and low dispersity (&#208: = 1.07) exhibit tunable mechanical properties characteristic of a strong and tough thermoplastic, elastomer, or a semicrystalline thermoplastic elastomer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Emulation and detection of physical faults and cyber-attacks on building energy systems through real-time hardware-in-the-loop experiments

The increasing use of remote or mobile access, integrated wearable technologies, data exchange, and cloud-based data analytics in modern smart buildings is steering the building industry towards open communication technologies. The increased connectivity and accessibility could lead to more cyber-attacks in smart buildings. On the other hand, physical faults (e.g., HVAC -heating, ventilation, and air-conditioning faults) may have similar adverse impacts as those from the cyber-attacks on building energy systems, such as occupant discomfort, energy wastage, and equipment downtime. However, current physical behavior-based anomaly detection methods fail to differentiate between cyber-attacks and physical faults in building energy systems. Moreover, the challenge in collecting real-world threat data with ground truth has led researchers to rely on numerical models with user-defined assumptions, which may not accurately reflect real-world conditions due to the lack of in-situ experimental datasets. To address these challenges and gaps, this paper presents a flexible hardware-in-the-loop (HIL) testbed for generating cyber-attack and physical fault datasets and demonstrating threat detection algorithms in a real building automation system (BAS) environment. This testbed combines hardware (i.e., real BAS with local HVAC controllers and a physical network) with software (i.e., high-fidelity models to represent behaviors of building envelope and HVAC energy systems), enabling emulations of realistic threats. Five HIL experiments, including one baseline without any threats, two with physical faults, and two with cyber-attacks, were conducted to generate datasets containing detailed network traffic and system states. A joint classification framework, incorporating a network analyzer and a physical HVAC fault detector, was proposed to automatically detect cyber-physical abnormalities on BAS at both the network and the physical HVAC levels. The network analyzer comprises a conditional random fields (CRF) based command validator and a statistics-based detection strategy. The fault detector employs a weather and schedule-based pattern matching and feature-based principal component analysis (WPM-FPCA) method. Evaluation of the classification using four metrics from the multi-class confusion matrix revealed an average accuracy of 90.2%, recall of 89.7%, precision of 88.5% and F1-score of 89.2%. Finally, these results demonstrate that the proposed joint classification framework can effectively differentiate between specific types of cyber-attacks (e.g., device reinitialization attack, network Denial-of-Service attack) and physical faults (e.g., air handling unit operational fault, cooling coil valve stuck) in real time for improved building energy management.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2020

Effects of cation stoichiometry on surface morphology and crystallinity of ZnGeN 2 films grown on GaN by metalorganic chemical vapor deposition

ZnGeN 2 films were grown on GaN-on-sapphire templates via metalorganic chemical vapor deposition. Energy dispersive X-ray spectroscopy was used to estimate the Zn/(Zn+Ge) composition ratio in the films. This ratio decreased with increase in growth temperature but increased with increase in total reactor pressure or Zn/Ge precursor flow rate ratio. Systematic mapping of these key growth parameters has allowed us to identify the growth window to achieve ZnGeN 2 with stoichiometric cation composition. Compositional and statistical analyses performed on data acquired from atom probe tomography provided insight into the local compositional homogeneity. The cations Zn and Ge did not demonstrate segregation or clustering at the sub-nanometer level. Based on X-ray diffraction 2θ-ω scan profiles and transmission electron microscope nano-diffraction patterns, the films with near-stoichiometric cation ratios were single crystalline with planar surfaces, whereas zinc-rich or zinc-poor films were polycrystalline with nonplanar surfaces. The growth direction of the single crystalline ZnGeN 2 films on GaN templates was along the c-axis. Room temperature Raman spectra showed features associated with the phonon density of states, indicating the presence of cation disorder in the lattice. A cathodoluminescence peak associated with transitions involving deep level defects was observed around 640 nm. The intensity of this peak increased by almost 2.5 times as the temperature was reduced to 77 K from room temperature. A similar peak was observed in photoluminescence spectra collected at 80 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Nanostructures for Electrical Energy Storage (NEES) (2020 Final Technical Report)

Nanostructures for Electrical Energy Storage (NEES, www.efrc.umd.edu) was an Energy Frontier Research Center supported by the DOE Office of Science, Basic Energy Sciences, from 8/1/2009 to 7/31/2020. Led by the University of Maryland, NEES enjoyed extensive collaborations with its funded partners, including two DOE Laboratories and six universities. The NEES vision has been to reveal a set of scientific insights and design principles that can underpin a next-generation electrical energy storage approach, building on advances in nanoscale science and technology to achieve simultaneous high power and high energy over extended charge/discharge cycling. The vision is motivated by the recognition that scaling into the nano regime opens the door to new physical phenomena and that the tools enlisted in nanoscale research provide major new opportunities for the synthesis not only of materials at molecular scale but for structures at nano scale and above. NEES has translated this vision into its research program based on two observations. First, while the behavior of ions and electrons in electrolytes and in electrode materials is crucial to electrical energy storage (or more appropriately electrochemical energy storage), it is the transport of ion and electron charge between different structural components of a storage device that ultimately determine its performance. With it well recognized that the choice of electrode materials typically constrain ion transport kinetics as well as maximum ion concentration, the search for better electrode materials has been a primary driver of battery research. At the same time the synthesis of electrodes is typically based on aggregation of particles with varying size, shape, and orientation in the electrode. Together with the presence of additional materials to impart electrical conductivity and cohesion to the composite electrode, change in electrode materials is necessarily accompanied by structural changes at the nano/micro scale that are difficult to categorize and manage. From the beginning, NEES’ vision has been to create and study simpler, highly controlled spatial arrangements of known materials as battery components (electrodes, current collectors, and electrolyte) and to understand how design and structure above the molecular scale determines the energy storage performance available from known materials. Second, advances in nanoscience dramatically expanded the portfolio of synthesis methods, structural motifs, and new phenomena available for research. Some of these gave rapid access to new building blocks at the deep nanoscale (e.g., carbon nanotubes grown by self-assembly, nanoscale arrays formed by electrochemical self-alignment, monolayer films controlled by self-limiting reaction). Such advances served as the enabler for the NEES vision to be pursued experimentally through study of 3D structures created and controlled at the nano, micro, and meso scales. Here, we use meso as in the BES MESO Report, implying not only intermediate or varying length scales, but very much the way behavior is influenced by other factors including aggregation of nanocomponents at different densities and spatial configurations, statistical variations in the aggregates, hierarchical architectures in which they can be assembled, or local 3D configurations that result from the architectures. Over its life cycle, NEES has pursued two overarching goals: (1) to understand the scientific fundamentals of electrochemical storage from the nanoscale to the mesoscale; and (2) to create and learn from innovative, controlled, heterogeneous nanostructures, where such nanostructures can enable the first goal and serve as models for future paradigms in energy storage. Specific goals have included: Synthesize heterogeneous nanostructures comprised of multiple materials arranged in controlled fashion and characterize their behavior; Demonstrate and elucidate design principles for achieving simultaneous high power and high energy; Develop materials processes which enable precision control of thin layers and 3D structures; Investigate the impact of artificial interphases on electrode stability during ion insertion/deinsertion; Create dense arrays of nanostructures to understand how the architecture of these assemblies, along with nanostructure design, influences energy storage behavior at the mesoscale; Identify and understand the consequences of nanoconfinement and local inhomogeneities in 3D mesoscale arrays; Develop and apply computational models to stimulate, guide and interpret experiments.

25 ENERGY STORAGE↗

DOE OSTI2023

Inference of the optical depth to reionization τ from Planck CMB maps with convolutional neural networks

The optical depth to reionization, τ, is the least constrained parameter of the cosmological Λ cold dark matter (ΛCDM) model. To date, its most precise value is inferred from large-scale polarized cosmic microwave background (CMB) power spectra from the High Frequency Instrument (HFI) aboard the Planck satellite. These maps are known to contain significant contamination by residual non-Gaussian systematic effects, which are hard to model analytically. Therefore, robust constraints on τ are currently obtained through an empirical cross-spectrum likelihood built from simulations. In this paper, we present a likelihood-free inference of τ from polarized Planck HFI maps which, for the first time, is fully based on neural networks (NNs). NNs have the advantage of not requiring an analytical description of the data and can be trained on state-of-the-art simulations, combining the information from multiple channels. By using Gaussian sky simulations and Planck SRoll2 simulations, including CMB, noise, and residual instrumental systematic effects, we trained, tested, and validated NN models considering different setups. We inferred the value of τ directly from Stokes Q and U maps at ~4° pixel resolution, without computing angular power spectra. On Planck data, we obtained τ NN = 0.0579 ± 0.0082, which is compatible with current EE cross-spectrum results but with a ~30% larger uncertainty, which can be assigned to the inherent nonoptimality of our estimator and to the retraining procedure applied to avoid biases. While this paper does not improve on current cosmological constraints on τ, our analysis represents a first robust application of NN-based inference on real data, and highlights its potential as a promising tool for complementary analysis of near-future CMB experiments, also in view of the ongoing challenge to achieve the first detection of primordial gravitational waves.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Toward quantum Monte Carlo forces on heavier ions: Scaling properties

Quantum Monte Carlo (QMC) forces have been studied extensively in recent decades because of their importance with spectroscopic observables and geometry optimization. Here, we benchmark the accuracy and computational cost of QMC forces. The zero-variance zero-bias (ZVZB) force estimator is used in standard variational and diffusion Monte Carlo simulations with mean-field based trial wavefunctions and atomic pseudopotentials. Statistical force uncertainties are obtained with a recently developed regression technique for heavy tailed QMC data [P. Lopez Rios and G. J. Conduit, Phys. Rev. E 99, 063312 (2019)]. By considering selected atoms and dimers with elements ranging from H to Zn (1 ≤ Z eff ≤ 20), we assess the accuracy and the computational cost of ZVZB forces as the effective pseudopotential valence charge, Z eff , increases. We find that the costs of QMC energies and forces approximately follow simple power laws in Zeff. The force uncertainty grows more rapidly, leading to a best case cost scaling relationship of approximately Z e ff 6.5 ( 3 ) for diffusion Monte Carlo. We find that the accessible system size at fixed computational cost scales as Z e ff − 2 , insensitive to model assumptions or the use of the “space warp” variance-reduction technique. Our results predict the practical cost of obtaining forces for a range of materials, such as transition metal oxides where QMC forces have yet to be applied, and underscore the importance of further developing force variance-reduction techniques, particularly for atoms with high Z eff .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Charge Transfer and Spin Dynamics in a Zinc Porphyrin Donor Covalently Linked to One or Two Naphthalenediimide Acceptors

We report quantum coherence effects on charge transfer and spin dynamics in a system having two degenerate electron acceptors are studied using a zinc 5,10,15-tri(n-pentyI)-20-phenylporphyrin (ZnP) electron donor covalently linked to either one or two naphthalene-1,8:4,5-bis(dicarboximide) (NDI) electron acceptors using an anthracene (An) spacer, ZnP-An-NDI (1) and ZnP-An-NDI 2 (2), respectively. Following photoexcitation of 1 and 2 in toluene at 295 K, femtosecond transient absorption spectroscopy shows that the electron transfer (ET) rate constant for 2 is about three times larger than that of 1, which can be accounted for by the statistical nature of incoherent ET as well as the electron couplings for the charge separation reactions. In contrast, the rate constant for charge recombination (CR) of 1 is about 25% faster than that of 2. Using femtosecond transient infrared spectroscopy and theoretical analysis, we find that the electron on NDI 2 ∙- in 2 localizes onto one of the two NDIs prior to CR, thus precluding electronically coherent CR from NDI 2 ∙- . Conversely, CR in both 1 and 2 is spin coherent as indicated by the observation of a resonance in the 3* ZnP yield following CR as a function of applied magnetic field, giving spin-spin exchange interaction energies of 2J = 210 and 236 mT, respectively, where the line width of the resonance for 2 is greater than 1. These data show that while CR is a spin-coherent process, incoherent hopping of the electron between the two NDIs in 2, consistent with the lack of delocalization noted above, results in greater spin decoherence in 2 relative to 1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Structure and Reactivity of Binuclear Cu Active Sites in Cu-CHA Zeolites for Stoichiometric Partial Methane Oxidation to Methanol

Aluminosilicate zeolites exchanged with copper ions facilitate partial methane oxidation (PMO) to methanol in stoichiometric oxidation and reduction cycles, yet the identities of active Cu sites and details of the reaction mechanism remain debated. Here, we use the high-symmetry chabazite (CHA) zeolite framework as a model support to probe the relationship between bulk composition, Cu speciation, and response to various oxidizing and reducing treatments. Density functional theory and first-principles thermodynamics combined with statistical models reveal that Cu speciation and composition depend strongly on Al configuration and external gas conditions. Cu-CHA samples were synthesized to survey broad regions of Si/Al and Cu/Al composition space and framework Al proximity. Characterization by in situ X-ray absorption and UV–visible spectroscopy during exposure to different oxidation conditions reveal that the extent of Cu oxidation is sensitive to activation conditions and thus that both kinetic and thermodynamic factors influence Cu oxidizability in a given material. Similar characterizations during CO reduction reveal that CO titrates Cu 2+ in amounts suggesting the presence of both O- and O 2 -bridged species. In contrast, CH 4 and autoreduction (He) treatments reduce similar but smaller numbers of Cu sites than CO, implicating O 2 -bridged Cu dimers as a potential common intermediate in the former reduction pathways. Furthermore, a systematic increase in methanol yields (per Cu) in stoichiometric PMO cycles increase with the fraction of binuclear O x -bridged Cu sites suggests these species as active sites, as depicted in an updated PMO reaction mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Multiple Reaction Pathways for the Oxygen Evolution Reaction May Contribute to IrO 2 (110)’s High Activity

Density functional theory calculations in conjunction with statistical mechanical arguments are performed on the rutile IrO 2 (110) facet in order to characterize multiple reaction pathways on the surface at the highest active limit (the stoichiometric surface with all metal sites available) and at the lowest active limit (the oxygen-terminated surface). Alternative pathways to the oxygen evolution reaction (OER) are found, with multiple pathways determined at each step of the four proton-coupled electron transfer reaction. Of particular interest is the detailed characterization of a co-adsorption pathway utilizing neighboring, adsorbed O, OH species in order to evolve oxygen; activation energies of this pathway are <0.5 eV and therefore easily surmountable at the high operating potentials of OER. We also determined that surface Ir atoms can potentially participate in deprotonating an OOH* intermediate; the activation energy to this is 0.67 eV on the oxygen-terminated surface. These theoretical findings explain in part the high activity present in iridium oxide catalysts and also provide insight into the mechanistic pathways available on metal oxide catalysts, which may require the concerted interaction of nearest neighbor co-adsorbates to produce chemicals of interest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

One Collision—Two Substituents: Gas–Phase Preparation of Xylenes under Single–Collision Conditions

The fundamental reaction pathways to the simplest dialkylsubstituted aromatics—xylenes (C 6 H 4 (CH 3 ) 2 )—in high-temperature combustion flames and in low-temperature extraterrestrial environments are still unknown, but critical to understand the chemistry and molecular mass growth processes in these extreme environments. Exploiting crossed molecular beam experiments augmented by state-of-the-art electronic structure and statistical calculations, this study uncovers a previously elusive, facile gas-phase synthesis of xylenes through an isomer-selective reaction of 1-propynyl (methylethynyl, CH 3 CC) with 2-methyl-1,3-butadiene (isoprene, C 5 H 8 ). The reaction dynamics are driven by a barrierless addition of the radical to the diene moiety of 2-methyl-1,3-butadiene followed by extensive isomerization (hydrogen shifts, cyclization) prior to unimolecular decomposition accompanied by aromatization via atomic hydrogen loss. This overall exoergic reaction affords a preparation of xylenes not only in high-temperature environments such as in combustion flames and around circumstellar envelopes of carbon-rich Asymptotic Giant Branch (AGB) stars, but also in low-temperature cold molecular clouds (10 K) and in hydrocarbon-rich atmospheres of planets and their moons such as Triton and Titan. Furthermore, our study established a hitherto unknown gas-phase route to xylenes and potentially more complex, disubstituted benzenes via a single collision event highlighting the significance of an alkyl-substituted ethynyl-mediated preparation of aromatic molecules in our Universe.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

One Collision—Two Substituents: Gas–Phase Preparation of Xylenes under Single–Collision Conditions

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Heterogeneous and Framework-Bound Copper Species Contribute to Catalytic Partial Methane Oxidation in Cu–Chabazite Zeolites

The relationship between continuous partial methane oxidation (PMO) rates and Cu site speciation in Cu-CHA zeolite catalysts is explored through differential rate measurements across a series of samples of varying compositions combined with density functional theory, first-principles thermodynamics, and statistical models that characterize Cu speciation. Under continuous PMO conditions (573 K, 0.07 kPa O 2 , 3 kPa H 2 O), Cu ions are shown to anchor to the CHA framework in both monomeric (Z 2 Cu and Z 2 CuH 2 O) and dimeric (O- and OH-bridged Cu) forms that are sensitive to the identity of the local framework anchoring site. Consequently, across the studied compositional range, Cu-CHA catalysts are predicted to contain a mixture of monomeric and dimeric Cu sites. Cu-normalized CH 3 OH formation rates extrapolated to zero conversion reflect contributions from multiple site types. Predicted CH 3 OH formation rates indicate that the Cu site reactivity toward CH 4 is influenced by zeolite composition and is likely limited by the reduction half-cycle.

03 NATURAL GAS↗