Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized functions”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

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

Electrochemical Hydrogenation of CO on Cu(100): Insights from Accurate Multiconfigurational Wavefunction Methods

Copper (Cu) remains the most efficacious electrocatalyst for electrochemical CO 2 reduction (CO 2 R). Its activity and selectivity are highly facet-dependent. We recently examined the commonly proposed rate-limiting CO hydrogenation step on Cu(111) via embedded correlated wavefunction (ECW) theory and demonstrated that only this higher-level theory yields predictions consistent with potential-dependent experimental kinetics. Here, to understand the differing activities of Cu(111) and Cu(100) in catalyzing CO 2 R, we explore CO hydrogenation on Cu(100) using ECW theory. We predict that the preferred pathway involves the reduction of adsorbed CO (*CO) to *COH via proton-coupled electron transfer (PCET) at working potentials, although *CHO also may form with a kinetically accessible but higher barrier. In contrast, our earlier work on Cu(111) concluded that *COH and *CHO formation via PCET are equally feasible. Furthermore, this work illustrates one possible origin of the facet dependence of CO 2 R mechanisms and products on Cu electrodes and sheds light on how the selectivity of CO 2 R electrocatalysts can be controlled by the surface morphology.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Multidimensional four-wave-mixing spectroscopy with squeezed light

A set of nonlinear spectroscopic measurements with quantum light is presented here. It involves a strong pump with frequency ω pu and a weak probe atω pr interacting with a solid state target to generate a conjugated beam ω c =2ω pu -ω pr via nondegenerate four-wave mixing. The $\chi^{(3)}$ susceptibility can be measured by the noise spectra of the intensity difference of the squeezed beams. We discuss three spectroscopic setups based on squeezed light: first, single four-wave mixing in a single crystal; second, cascading scheme involving two crystals; third, an $\textit{SU}$(1, 1) interferometer based on two separate four-wave mixing processes. We further investigate the microscopic noise and optical losses in all three setups. Simulations are presented for silicon-vacancy color centers in diamond.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Density dependence of the excitation gaps in an undoped Si/SiGe double-quantum-well heterostructure

In this work, we report low-temperature magneto-transport measurements of an undoped Si/SiGe asymmetric double quantum well heterostructure. The density in both layers is tuned independently utilizing top and bottom gates, allowing the investigation of quantum wells at both imbalanced and matched densities. Integer quantum Hall states at total filling factor v T =1 and v T =2 are observed in both density regimes, and the evolution of their excitation gaps is reported as a function of the density. The v T =1 gap evolution departs from the behavior generally observed for valley splitting in the single layer regime. Furthermore, by comparing the v T =2 gap to the single particle tunneling energy, Δ SAS , obtained from Schrödinger–Poisson (SP) simulations, evidence for the onset of spontaneous interlayer coherence is observed for a relative filling fraction imbalance smaller than ~50%.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Toward the core-edge coupling of delta-f and total-f gyrokinetic models

The coupling of core reduced delta-f and edge total-f gyrokinetic models could enable a significant speed-up of the whole device modeling simulations. It is thus interesting to compare these models in the core where they will be coupled. The main difference between these models is the presence of the zeroth order term on the right-hand side of the total-f gyrokinetic equation. This term is responsible for driving a radial electric field that is also associated with Geodesic acoustic mode-like oscillations. To investigate the coupling between these models, the subsequent large-scale gyrokinetic 3D turbulence simulations have been performed starting from an already saturated gyrokinetic axisymmetric equilibria (collisions are not included). This work has also been extended to couple different axisymmetric and turbulent models available in XGC. It is shown that the radial electric field and its drive have to be taken care of consistently while coupling different models together.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Computational evolution of high-performing unfused non-fullerene acceptors for organic solar cells

Materials optimization for organic solar cells (OSCs) is a highly active field, with many approaches using empirical experimental synthesis, computational brute force to screen a subset of chemical space, or generative machine learning methods that often require significant training sets. While these methods may find high-performing materials, they can be inefficient and time-consuming. Genetic algorithms (GAs) are an alternative approach, allowing for the “virtual synthesis” of molecules and a prediction of their “fitness” for some property, with new candidates suggested based on good characteristics of previously generated molecules. In this work, a GA is used to discover high-performing unfused non-fullerene acceptors (NFAs) based on an empirical prediction of power conversion efficiency (PCE) and provides design rules for future work. The electron-withdrawing/donating strength, as well as the sequence and symmetry, of those units are examined. The utilization of a GA over a brute-force approach resulted in speedups up to 1.8 × 10 12 . New types of units, not frequently seen in OSCs, are suggested, and in total 5426 NFAs are discovered with the GA. Of these, 1087 NFAs are predicted to have a PCE greater than 18%, which is roughly the current record efficiency. While the symmetry of the sequence showed no correlation with PCE, analysis of the sequence arrangement revealed that higher performance can be achieved with a donor core and acceptor end groups. Future NFA designs should consider this strategy as an alternative to the current A-D-A'-D-A architecture.

14 SOLAR ENERGY↗
DOE OSTI2022

Quasilinear theory of general electromagnetic fluctuations including discrete particle effects for magnetized plasmas: General analysis

The general quasilinear Fokker–Planck kinetic equation for the gyrophase-averaged plasma particle distribution functions in magnetized plasmas is derived, making no restrictions on the energy of the particles and on the frequency of the electromagnetic fluctuations and avoiding the often made Coulomb approximation of the electromagnetic interactions. Here, the inclusion of discrete particle effects breaks the dichotomy of nonlinear kinetic plasma theory divided into the test particle and the test fluctuation approximation because it provides expression of both the non-collective and collective electromagnetic fluctuation spectra in terms of the plasma particle distribution functions. Within the validity of the quasilinear approach, the resulting full quasilinear transport equation can be regarded as a determining nonlinear equation for the time evolution of the plasma particle distribution functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond

Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within “chemical accuracy” (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob’s ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 − N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.

Chemistry↗
DOE OSTI2024

Erosion and impurity transport for the edge localized mode suppression window in KSTAR

A 3D analysis of plasma wall interactions and global impurity transport for the edge localized mode suppression window in KSTAR during H-Mode discharges has shown that carbon erosion at the divertor plates is a strong function of the resonant magnetic perturbation (n = 1) coil current and relative phasing. The Generalized Perturbation Equilibrium Code was used to determine a realistic initial perturbed plasma equilibrium, and EMC3-EIRENE was used to calculate the resulting scrape-off layer plasma used in this study as a fixed background for the ERO2.0 plasma–material interaction model. The resulting transport leads to deposition of impurities along the targets positioned at the high-field side of the device. An attempt at calculating the resulting effective charge state has demonstrated a similar dependence on the perturbation coil current and has been able to determine a window for the experimentally observed values of Zeff by including contributions of all ionized carbon charge states and deuterium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Al-alkyls as acceptor dopant precursors for atomic-scale devices

Atomically precise ultradoping of silicon is possible with atomic resists, area-selective surface chemistry, and a limited set of hydride and halide precursor molecules, in a process known as atomic precision advanced manufacturing (APAM). It is desirable to expand this set of precursors to include dopants with organic functional groups and in this work we consider aluminium alkyls, to expand the applicability of APAM. We explore the impurity content and selectivity that results from using trimethyl aluminium and triethyl aluminium precursors on Si(001) to ultradope with aluminium through a hydrogen mask. Comparison of the methylated and ethylated precursors helps us understand the impact of hydrocarbon ligand selection on incorporation surface chemistry. Combining scanning tunneling microscopy and density functional theory calculations, we assess the limitations of both classes of precursor and extract general principles relevant to each.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Robust wind farm layout optimization

Wake interactions in wind farms cause losses in annual energy production (AEP) on the order of 10%. Wind farm designers optimize the layout of the farm to mitigate wake losses, especially in the dominant site-specific wind directions. As wind turbines and wind farms grow in scale, optimization becomes more complex. Offshore wind farms regularly comprise more than 100 wind turbines and are characterized by complex boundaries due to shipping lanes, neighboring wind farms, and other constraints. Layout optimization methods are broadly split between gradient-based and gradient-free approaches. Gradient-based approaches can converge quickly and perform well for smaller, academic problems but are often sensitive to initial conditions and tuning parameters and require expert knowledge to use. On the other hand, gradient-free approaches can be more robust to problem complexities. We present a robust layout optimization approach based on a random search algorithm. The algorithm is intended for those who are not optimization experts and has few tuning parameters that need specification to achieve satisfactory results. Unlike off-the-shelf methods, which use generally available, non-domain-specific optimization routines that accept as inputs an optimization function and constraint definitions, this approach takes advantage of the relative computational costs of the different evaluations by evaluating cheaper computations first (boundary and minimum distance constraints) and running expensive AEP evaluations only if all other checks pass. Moreover, an outer genetic algorithm allows multiple solutions to evolve in parallel, enabling rapid solution development on high-performance computers. We discuss the relative ease of selecting necessary tuning parameters and demonstrate the efficacy of the genetic random search on a complex layout problem consisting of placing 70 turbines in a nonconvex and unconnected boundary region.

17 WIND ENERGY↗
DOE OSTI2023

Certifiably Robust Neural ODE With Learning-Based Barrier Function

Neural Ordinary Differential Equations (ODEs) have gained traction in many applications. While recent studies have focused on empirically increasing the robustness of neural ODEs against natural or adversarial attacks, certified robustness is still lacking. In this work, we propose a framework for training a neural ODE using barrier functions and demonstrate improved robustness for classification problems. Finally, we further provide the first generalization guarantee of robustness against adversarial attacks using a wait-and-judge scenario approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Bloom filters for molecules

Abstract Ultra-large chemical libraries are reaching 10s to 100s of billions of molecules. A challenge for these libraries is to efficiently check if a proposed molecule is present. Here we propose and study Bloom filters for testing if a molecule is present in a set using either string or fingerprint representations. Bloom filters are small enough to hold billions of molecules in just a few GB of memory and check membership in sub milliseconds. We found string representations can have a false positive rate below 1% and require significantly less storage than using fingerprints. Canonical SMILES with Bloom filters with the simple FNV (Fowler-Noll-Voll) hashing function provide fast and accurate membership tests with small memory requirements. We provide a general implementation and specific filters for detecting if a molecule is purchasable, patented, or a natural product according to existing databases at https://github.com/whitead/molbloom .

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

Safety Evaluation of Nitric Acid Reactions with Non-Polysaccharide Organic Materials

The hazards of non-polysaccharide materials in a transuranic (TRU) waste drum exposed to nitric acid (HNO 3 ) and metal nitrate salts have been evaluated, focusing on sorbents and on resins used in ion-exchange chromatography. The range of sorbent materials can be grouped into two general categories: 1) a variety of polyacrylate and polyacrylamide compounds that incorporate polar carbonyl functional groups designed to sorb protic (i.e., acidic) substrates and solutions; 2) hydrocarbon based polymers that include polystyrene, polybutadiene, and polyethylene derivatives, which are designed to sorb non-polar organic substrates and solutions. The ion-exchange resins are constructed of hydrocarbon based polymeric networks equipped with pendant ionizable functional groups designed to reversibly sorb/desorb ionic species that are targeted for separation from a mobile phase.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Realistic Noise Generation to Enhance Realism of Virtual Lidar Scans

Many real-world phenomena corrupt light detection and ranging (lidar) measurements, such as laser energy attenuation, variations in aerosol concentration and composition with height, and hard target returns. Accurate studies of lidar scans using virtual lidar methods should include some realistic model of these corrupting effects to generate more realistic simulations of lidar scans. We present a simple model that characterizes noise caused by energy attenuation and aerosol stratification. The model requires limited inputs and is developed for a Halo Photonics Streamline XR lidar but is readily generalizable for other lidar systems. A critical component of this model is a model of the standard deviation of measured wind speed as a function of the backscattered signal’s signal-to-noise ratio. We derive a general model for this behavior that can be adapted to different scan settings.

17 WIND ENERGY↗
DOE OSTI2021

Computational Fermi Level Engineering and Doping-Type Conversion of Ga2O3 Via Three-Step Processing

Ga2O3 is being actively explored for power electronics, deep-ultraviolet optoelectronics, and other applications due to its ultra-wide bandgap and low projected fabrication cost of large-size and high-quality crystals. N-type doping of Ga2O3 can be achieved and tuned, but p-type doping faces fundamental obstacles due to deep character of acceptor levels and polaron transport of resulting holes. However, successful engineering of Ga2O3 based devices requires critical control of doping density, Fermi level position, and free carrier concentration, providing opportunities for predictive process simulation. We use first-principles defect theory and defect equilibrium calculations to simulate a 3-step growth-annealing-quench protocol for hydrogen assisted Mg doping in Ga2O3, taking into account the hydrogen-oxygen-water gas phase equilibrium. We predict type conversion to a net p-type regime following O-rich annealing after growth under reducing conditions in the presence of H2. This process is similar to the Mg acceptor activation by H removal in GaN. We show that there is an optimal temperature that maximizes the net acceptor density during the equilibrium annealing step for a given Mg doping level. Quenching of non-equilibrium annealed samples then results in a Fermi level EF below mid-gap down to about EV +1.5 eV, creating a significant number of uncompensated neutral MgGa0 acceptors. The resulting free hole concentration in quenched samples is very low (109 – 1012 cm-3) due to deep energy level of these Mg acceptors, but this type converted Ga2O3 material can create a significant built-in field in a p-n junction with an adjoining n-type material. Additionally, the electron concentration is greatly suppressed in such acceptor-doped Ga2O3, which could enable the use as a current blocking layer to fabricate normally-off (enhancement-mode) vertical Ga2O3 based metal-oxide-semiconductor field effect transistors (MOSFETs), and as a guard ring for edge termination in MOSFETs and Schottky barrier diodes (SBDs) with increased breakdown voltage.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,M↗