Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy predictions”

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 559 records · Page 31

Subseasonal Tropical Cyclone Prediction and Modulations by MJO and ENSO in CESM2

Subseasonal tropical cyclone (TC) reforecasts from the Community Earth System Model version 2 (CAM6) subseasonal prediction system are examined in this study. Here, we evaluate the modeled TC climatology and the probabilistic forecast skill of basin-wide TC genesis at weekly temporal resolution. Prediction skill is calculated using the Brier skill score relative to a constant annual mean climatology and to a monthly varying seasonal climatology during TC season. The model captures the observed basin-wide climatological TC seasonality and spatial distributions at weeks 1–6, but TC genesis is largely underestimated from Week 2 onward. For some basins and lead times, the predicted TC genesis is primarily controlled by the number of TC “seeds” and the mean-state climate condition. The model has good prediction skill relative to the constant climatology across all the basins and lead times, but is only skillful in the eastern Pacific, North Indian Ocean, and Southern Hemisphere at Week 1 when compared to the seasonal climatology, indicating limited skill in predicting deviations from the seasonal cycle. We find strong modulations of the predicted TC genesis at up to 3 weeks of forecast lead time by the Madden-Julian Oscillation. The interannual variability of predicted TC genesis and accumulated cyclone energy are skillfully predicted in the North Atlantic and the Northwestern Pacific, with a strong modulation by the El Nino-Southern Oscillation.

54 ENVIRONMENTAL SCIENCES↗

Teaching a neural network to attach and detach electrons from molecules

Abstract Interatomic potentials derived with Machine Learning algorithms such as Deep-Neural Networks (DNNs), achieve the accuracy of high-fidelity quantum mechanical (QM) methods in areas traditionally dominated by empirical force fields and allow performing massive simulations. Most DNN potentials were parametrized for neutral molecules or closed-shell ions due to architectural limitations. In this work, we propose an improved machine learning framework for simulating open-shell anions and cations. We introduce the AIMNet-NSE (Neural Spin Equilibration) architecture, which can predict molecular energies for an arbitrary combination of molecular charge and spin multiplicity with errors of about 2–3 kcal/mol and spin-charges with error errors ~0.01e for small and medium-sized organic molecules, compared to the reference QM simulations. The AIMNet-NSE model allows to fully bypass QM calculations and derive the ionization potential, electron affinity, and conceptual Density Functional Theory quantities like electronegativity, hardness, and condensed Fukui functions. We show that these descriptors, along with learned atomic representations, could be used to model chemical reactivity through an example of regioselectivity in electrophilic aromatic substitution reactions.

36 MATERIALS SCIENCE↗

Nickase fidelity drives EvolvR-mediated diversification in mammalian cells

Abstract In vivo genetic diversifiers have previously enabled efficient searches of genetic variant fitness landscapes for continuous directed evolution. However, existing genomic diversification modalities for mammalian genomic loci exclusively rely on deaminases to generate transition mutations within target loci, forfeiting access to most missense mutations. Here, we engineer CRISPR-guided error-prone DNA polymerases (EvolvR) to diversify all four nucleotides within genomic loci in mammalian cells. We demonstrate that EvolvR generates both transition and transversion mutations throughout a mutation window of at least 40 bp and implement EvolvR to evolve previously unreported drug-resistantMAP2K1variants via substitutions not achievable with deaminases. Moreover, we discover that the nickase’s mismatch tolerance limits EvolvR’s mutation window and substitution biases in a gRNA-specific fashion. To compensate for gRNA-to-gRNA variability in mutagenesis, we maximize the number of gRNA target sequences by incorporating a PAM-flexible nickase into EvolvR. Finally, we find a strong correlation between predicted free energy changes underlying R-loop formation and EvolvR’s performance using a given gRNA. The EvolvR system diversifies all four nucleotides to enable the evolution of mammalian cells, while nuclease and gRNA-specific properties underlying nickase fidelity can be engineered to further enhance EvolvR’s mutation rates.

Science & Technology - Other Topics↗

Atomistic learning in the electronically grand-canonical ensemble

Abstract A strategy is presented for the machine-learning emulation of electronic structure calculations carried out in the electronically grand-canonical ensemble. The approach relies upon a dual-learning scheme, where both the system charge and the system energy are predicted for each image. The scheme is shown to be capable of emulating basic electrochemical reactions at a range of potentials, and coupling it with a bootstrap-ensemble approach gives reasonable estimates of the prediction uncertainty. The method is also demonstrated to accelerate saddle-point searches, and to extrapolate to systems with one to five water layers. We anticipate that this method will allow for larger length- and time-scale simulations necessary for electrochemical simulations.

36 MATERIALS SCIENCE↗

Toward fully automated UED operation using two-stage machine learning model

To demonstrate the feasibility of automating UED operation and diagnosing the machine performance in real time, a two-stage machine learning (ML) model based on self-consistent start-to-end simulations has been implemented. This model will not only provide the machine parameters with adequate precision, toward the full automation of the UED instrument, but also make real-time electron beam information available as single-shot nondestructive diagnostics. Furthermore, based on a deep understanding of the root connection between the electron beam properties and the features of Bragg-diffraction patterns, we have applied the hidden symmetry as model constraints, successfully improving the accuracy of energy spread prediction by a factor of five and making the beam divergence prediction two times faster. The capability enabled by the global optimization via ML provides us with better opportunities for discoveries using near-parallel, bright, and ultrafast electron beams for single-shot imaging. It also enables directly visualizing the dynamics of defects and nanostructured materials, which is impossible using present electron-beam technologies.

36 MATERIALS SCIENCE↗

New thermal decomposition pathway for TATB

Abstract Understanding the thermal decomposition behavior of TATB (1,3,5-triamino-2,4,6-trinitrobenzene) is a major focus in energetic materials research because of safety issues. Previous research and modelling efforts have suggested benzo-monofurazan condensation producing H 2 O is the initiating decomposition step. However, early evolving CO 2 (m/z 44) along with H 2 O (m/z 18) evolution have been observed by mass spectrometric monitoring of head-space gases in both constant heating rate and isothermal decomposition studies. The source of the CO 2 has not been explained, until now. With the recent successful synthesis of 13 C 6 -TATB ( 13 C incorporated into the benzene ring), the same experiments have been used to show the source of the CO 2 is the early breakdown of the TATB ring, not adventitious C from impurities and/or adsorbed CO 2 . A shift in mass m/z 44 (CO 2 ) to m/z 45 is observed throughout the decomposition process indicating the isotopically labeled 13 C ring breakdown occurs at the onset of thermal decomposition along with furazan formation. Partially labeled (N 18 O 2 ) 3 -TATB confirms at least some of the oxygen comes from the nitro-groups. This finding has a significant bearing on decomposition computational models for prediction of energy release and deflagration to detonation transitions, with respect to conditions which currently do not recognize this oxidation step.

36 MATERIALS SCIENCE↗

Sub 20 cm -1 computational prediction of the CH bond energy – a case of systematic error in computational thermochemistry

Here, the bond dissociation energy of methylidyne, D 0 (CH), is studied using an improved version of the High-Accuracy Extrapolated ab initio Thermochemistry (HEAT) approach as well as the Feller–Peterson–Dixon (FPD) model chemistry. These calculations, which include basis sets up to nonuple (aug-cc-pCV9Z) quality, are expected to be capable of providing results substantially more accurate than the ca. 1 kJ mol -1 level that is characteristic of standard high-accuracy protocols for computational thermochemistry. The calculated 0 K CH bond energy (27 954 ± 15 cm -1 for HEAT and 27 956 ± 15 cm -1 for FPD), along with equivalent treatments of the CH ionization energy and the CH + dissociation energy (85 829 ± 15 cm -1 and 32 946 ± 15 cm -1 , respectively), were compared to the existing benchmarks from Active Thermochemical Tables (ATcT), uncovering an unexpected difference for D 0 (CH). This has prompted a detailed reexamination of the provenance of the corresponding ATcT benchmark, allowing the discovery and subsequent correction of a systematic error present in several published high-level calculations, ultimately yielding an amended ATcT benchmark for D 0 (CH). Finally, the current theoretical results were added to the ATcT Thermochemical Network, producing refined ATcT estimates of 27 957.3 ± 6.0 cm -1 for D0(CH), 32 946.7 ± 0.6 cm -1 for D 0 (CH + ), and 85 831.0 ± 6.0 cm -1 for IE(CH).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solvation at metal/water interfaces: An ab initio molecular dynamics benchmark of common computational approaches

Determining the influence of the solvent on electrochemical reaction energetics is a central challenge in our understanding of electrochemical interfaces. To date, it is unclear how well existing methods predict solvation energies at solid/liquid interfaces, since they cannot be assessed experimentally. Ab initio molecular dynamics (AIMD) simulations present a physically highly accurate, but also a very costly approach. In this work, we employ extensive AIMD simulations to benchmark solvation at charge-neutral metal/water interfaces against commonly applied continuum solvent models. We consider a variety of adsorbates including *CO, *CHO, *COH, *OCCHO, *OH, and *OOH on Cu, Au, and Pt facets solvated by water. The surfaces and adsorbates considered are relevant, among other reactions, to electrochemical CO 2 reduction and the oxygen redox reactions. We determine directional hydrogen bonds and steric water competition to be critical for a correct description of solvation at the metal/water interfaces. As a consequence, we find that the most frequently applied continuum solvation methods, which do not yet capture these properties, do not presently provide more accurate energetics over simulations in vacuum. We find most of the computed benchmark solvation energies to linearly scale with hydrogen bonding or competitive water adsorption, which strongly differ across surfaces. Furthermore, we determine solvation energies of adsorbates to be non-transferable between metal surfaces, in contrast to standard practice.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Connected three-body terms in single-reference unitary many-body theories: Iterative and perturbative approximations

This work introduces various approaches to include connected three-body terms in unitary many-body theories, focusing on the driven similarity renormalization group (DSRG). Starting from the least approximate method—the linearized DSRG truncated to one-, two-, and three-body operators [LDSRG(3)]—we develop several approximate LDSRG(3) models with reduced computational cost. Through a perturbative analysis, we motivate a family of iterative LDSRG(3)-n and -n' (n = 1, 2, 3, 4) methods that contain a subset of the LDSRG(3) diagrams. Among these variants, the LDSRG(3)-2 scheme has the same computational complexity of coupled cluster theory with singles, doubles, and triples (CCSDT), but it outperforms CCSDT in the accuracy of the predicted correlation energies. We also propose and implement two perturbative triples corrections based on the linearized DSRG truncated to one- and two-body operators augmented with recursive semi-quadratic commutators [qDSRG(2)]. Overall, the resulting qDSRG(2)+(T) approach matches the accuracy of the “gold-standard” coupled cluster theory with singles, doubles, and perturbative triples model on the energetics of twenty-eight closed-shell atoms and small molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparing GGA, GGA+ U , and meta-GGA functionals for redox-dependent binding at open metal sites in metal–organic frameworks

Metal–organic frameworks (MOFs) with open metal sites have been widely investigated for the selective adsorption of small molecules via redox mechanisms where charge transfer can take place between the binding site and the adsorbate of interest. Quantum-chemical screening methods based on density functional theory have emerged as a promising route to accelerate the discovery of MOFs with enhanced binding affinities toward various adsorbates. However, the success of this approach is linked to the accuracy of the underlying density functional approximations (DFAs). In this work, we compare commonly used generalized gradient approximation (GGA), GGA+U, and meta-GGA exchange-correlation functionals in modeling redox-dependent binding at open metal sites in MOFs using O 2 and N 2 as representative small molecules. We find that the self-interaction error inherent to the widely used Perdew, Burke, and Ernzerhof (PBE) GGA predicts metal sites that are artificially redox-active, as evidenced by their strong binding affinities, short metal–adsorbate bond distances, and large degree of charge transfer. The incorporation of metal-specific, empirical Hubbard U corrections based on the transition metal oxide literature systematically reduces the redox activity of the open metal sites, often improving agreement with experiment. Additionally, the binding behavior shifts from strong chemisorption to weaker physisorption as a function of U. The M06-L meta-GGA typically predicts binding energies between those of PBE-D3(BJ) and PBE-D3(BJ)+U when using empirically derived U values from the transition metal oxide literature. Despite the strong sensitivity of the binding affinities toward a given DFA, the GGA, GGA+U, and meta-GGA approaches often yield the same qualitative trends and structure–property relationships.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An ab initio exciton model for singlet fission

Here we present an ab initio exciton model that extends the Frenkel exciton model and includes valence, charge-transfer, and multiexcitonic excited states. It serves as a general, parameter-free, yet computationally efficient and scalable approach for simulation of singlet fission processes in multichromophoric systems. A comparison with multiconfigurational methods confirms that our exciton model predicts consistent energies and couplings for the pentacene dimer and captures the correct physics. Calculations of larger pentacene clusters demonstrate the computational scalability of the exciton model and suggest that the mixing between local and charge-transfer excitations narrows the gap between singlet and multiexcitonic states. Local vibrations of pentacene molecules are found to facilitate singlet–multiexcitonic state-crossing and hence are important for understanding singlet fission. The exciton model developed in this work also sets the stage for further implementation of the nuclear gradients and nonadiabatic couplings needed for first principles nonadiabatic quantum molecular dynamics simulations of singlet fission.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First principles reactive simulation for equation of state prediction

The high cost of density functional theory (DFT) has hitherto limited the ab initio prediction of the equation of state (EOS). In this article, we employ a combination of large scale computing, advanced simulation techniques, and smart data science strategies to provide an unprecedented ab initio performance analysis of the high explosive pentaerythritol tetranitrate (PETN). Comparison to both experiment and thermochemical predictions reveals important quantitative limitations of DFT for EOS prediction and thus the assessment of high explosives. In particular, we find that DFT predicts the energy of PETN detonation products to be systematically too high relative to the unreacted neat crystalline material, resulting in an underprediction of the detonation velocity, pressure, and temperature at the Chapman–Jouguet state. The energetic bias can be partially accounted for by high-level electronic structure calculations of the product molecules. Furthermore we demonstrate a modeling strategy for mapping chemical composition across a wide parameter space with limited numerical data, the results of which suggest additional molecular species to consider in thermochemical modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A molecular dynamics study of laser-excited gold

The structural evolution of laser-excited systems of gold has previously been measured through ultrafast MeV electron diffraction. However, there has been a long-standing inability of atomistic simulations to provide a consistent picture of the melting process, leading to large discrepancies between the predicted threshold energy density for complete melting, as well as the transition between heterogeneous and homogeneous melting. We make use of two-temperature classical molecular dynamics simulations utilizing three highly successful interatomic potentials and reproduce electron diffraction data presented by Mo et al. [Science 360, 1451–1455 (2018)]. We recreate the experimental electron diffraction data, employing both a constant and temperature-dependent electron–ion equilibration rate. In all cases, we are able to match time-resolved electron diffraction data, and find consistency between atomistic simulations and experiments, only by allowing laser energy to be transported away from the interaction region. This additional energy-loss pathway, which scales strongly with laser fluence, we attribute to hot electrons leaving the target on a timescale commensurate with melting.

36 MATERIALS SCIENCE↗

Plateau–Rayleigh instability with a grain boundary twist

We demonstrate using theoretical, computational, and experimental studies a morphological instability, in which a polycrystalline nanorod breaks up at grain boundaries (GBs) into an array of isolated domains. Our theoretical model is used to establish a neutral stability surface demarcating stable and unstable perturbations. It is shown that GBs play a destabilizing role in which the critical wavelength for the instability decreases with the increase in the GB energy. We carry out phase field simulations, which reveal accelerated pinch-off kinetics with the increase in the GB energy and predict temporal evolution of interfacial profiles in quantitative agreement with experimental observations.

Hussein, Omar↗

Identification and characterization of deep nitrogen acceptors in β-Ga 2 O 3 using defect spectroscopies

The ability to achieve highly resistive beta-phase gallium oxide (β-Ga 2 O 3 ) layers and substrates is critical for β-Ga 2 O 3 high voltage and RF devices. To date, the most common approach involves doping with iron (Fe), which generates a moderately deep acceptor-like defect state located at E C -0.8 eV in the β-Ga 2 O 3 bandgap. Recently, there has been growing interest in alternative acceptors, such as magnesium (Mg) and nitrogen (N), due to their predicted deeper energy levels, which could avoid inadvertent charge modulation during device operation. In this work, a systematic study that makes direct correlations between the introduction of N using ion implantation and the observation of a newly observed deep level at E C -2.9 eV detected by deep-level optical spectroscopy (DLOS) is presented. The concentration of this state displayed a monotonic dependence with N concentration over a range of implant conditions, as confirmed by secondary ion mass spectrometry (SIMS). With a near 1:1 match in absolute N and E C -2.9 eV trap concentrations from SIMS and DLOS, respectively, which also matched the measured removal of free electrons from capacitance-voltage studies, this indicates that N contributes a very efficiently incorporated compensating defect. Density functional theory calculations confirm the assignment of this state to be an N (0/-1) acceptor with a configuration of N occupying the oxygen site III [N O(III) ]. The near ideal efficiency for this state to compensate free electrons and its location toward the midgap region of the β-Ga 2 O 3 bandgap demonstrates the potential of N doping as a promising approach for producing semi-insulating β-Ga 2 O 3 .

36 MATERIALS SCIENCE↗

Large-database cross-verification and validation of tokamak transport models using baselines for comparison

State-of-the-art 1D transport solvers ASTRA and TRANSP are verified, then validated across a large database of semi-randomly selected, time-dependent DIII-D discharges. Various empirical models are provided as baselines to contextualize the validation figures of merit using statistical hypothesis tests. For predicting plasma temperature profiles, no statistically significant advantage is found for the ASTRA and TRANSP simulators over a baseline empirical (two-parameter) model. For predicting stored energy, a significant advantage is found for the simulators over a baseline empirical model based on confinement time scaling. Uncertainty in the results due to diagnostic and profile fitting uncertainties is approximated and determined to be insignificant due in part to the large quantity of discharges employed in the study. Advantages are discussed for validation methodologies like this one that employ (1) large databases and (2) baselines for comparison that are specific to the intended use-case of the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Unnuclear physics: Conformal symmetry in nuclear reactions

Significance Symmetry plays a key role in modern physics as a guiding principle for fundamental theories of nature. We use the nonrelativistic conformal symmetry to predict universal energy spectra of final states in reactions of quantum particles at vastly different length scales. Our work extends Georgi’s “unparticle” proposal to nonrelativistic spin-1/2 fermions in the unitary regime of strong interactions. These systems form so-called unparticles, which behave characteristically different from normal particles. They can be engineered with ultracold atoms and occur naturally in nuclear reactions with multiple neutrons—a situation we dub “unnuclear physics.” We present the general scenario of unnuclear physics, apply it to identify unnuclei in nuclear reactions, and highlight the opportunity to measure unnuclei at radioactive beam facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Morphogenesis of spin cycloids in a noncollinear antiferromagnet

Pattern formation in spin systems with continuous-rotational symmetry (CRS) provides a powerful platform to study emergent complex magnetic phases and topological defects in condensed-matter physics. However, its understanding and correlation with unconventional magnetic order along with high-resolution nanoscale imaging are challenging. Here, we employ scanning nitrogen vacancy (NV) magnetometry to unveil the morphogenesis of spin cycloids at both the local and global scales within a single ferroelectric domain of (111)-oriented BiFeO 3 , which is a noncollinear antiferromagnet, resulting in formation of a glassy labyrinthine pattern. We find that the domains of locally oriented cycloids are interconnected by an array of topological defects and exhibit isotropic energy landscape predicted by first-principles calculations. We propose that the CRS of spin-cycloid propagation directions within the (111) drives the formation of the labyrinthine pattern and the associated topological defects such as antiferromagnetic skyrmions. Unexpectedly, reversing the as-grown ferroelectric polarization from [$\overline{1}$ $\overline{1}$ $\overline{1}$] to [111] produces a noncycloidal NV image contrast which ¯ could be attributed to either the emergence of a uniformly magnetized state or a reversal of the cycloid polarity. These findings highlight that (111)-oriented BiFeO 3 is not only important for studying the fascinating subject of pattern formation but could also be utilized as an ideal platform for integrating novel topological defects in the field of antiferromagnetic spintronics.

magnetoelectric↗