Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalization”

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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The role of peer review and science for advancing public understanding

We live in a time of rapid change. The Covid-19 pandemic has united the world to understand and develop solutions to protect the health and wellbeing of society. At the same time, access to information on the internet and social media has exploded. People are exposed to information with many claims. Some of the information is relevant and accurate. A significant fraction may be inaccurate, biased, and potentially dangerous. With Covid-19, the information might point to cures that are toxic or risks that are not relevant. The scientific community has developed a process to develop and convey knowledge. The validation of results and their interpretation is critical to advancing knowledge. The core of the scientific process is based on developing a hypothesis-driven concept, designing and executing an experimental approach to prove or disprove the hypothesis, conducting a statistical analysis of the results to confirm their validity, and interpreting their meaning and impact. When their work is completed the scientist prepares to publish the work. Publication requires peer review. Hypotheses, experimental approaches, results, and interpretation of results are independently and anonymously reviewed by a scientific peer. The reviewers are selected based on their expertise in the field. The peer review process is essential for the integrity of the scientific process. Without a robust scientific process, society cannot discriminate between validated results and interpretations and biased, inaccurate, and potentially dangerous conjecture. With inaccurate science society and policy makers may make poor decisions that could hurt people and the world around them. While scientists should not censor information available to the public, scientists have a responsibility to inform the public. We must articulate what the scientific process is and why it is essential when considering the validity of information. Scientist are responsible for pointing out when public information has not been vetted through a scientific process and therefore cannot be considered conclusive. Scientists are responsible for highlighting information that has been generated and vetted with the scientific process. Over time, our goal is for the public to value the scientific process based on peer review. The public could become enlightened to the necessity for peer review to support decisions and action. Scientists need to contribute to advancing public knowledge. When invited to serve as a peer reviewer, scientists have a responsibility to contribute to the process. We all are extremely busy making peer review a distraction. While peer review does not directly benefit our own work and career, it benefits our chosen field of science. This is very important for the research area that our journal, Water-Energy Nexus, is focusing on. We cover the convergence of science and technology integrating across large spatial ranges (micro to global) in environmental pollution and the balance of remediation issues with energy use. Water-Energy Nexus wants to provide scientists with a platform for scientific debate and discussion on this rich and complicated topic. With the manuscripts submitted to WEN cover integrated research topics across a broad range of subjects from environment sciences through energy engineering, the role of peer review is essential. Peer review on a new submission to WEN should not be biased towards a specific science field. It should be done considering the complex nature of issues in the water-energy nexus. Scientists need to contribute to the public dialog. To accomplish these goals, we need to be honest, accurate, articulate, and use understandable language. The impact of WEN is very diverse affecting society in many ways. The scientific audience for WEN is diverse and readers may not fully understand the content of the papers outside their research expertise. Therefore, scientists, especially those who want to contribute to WEN, should work with communications specialists to convey message to benefit society. The audience includes the general public and researchers outside the author’s core specialty. While the peer-reviewed journal article is the cornerstone of our careers and critical for advancing scientific knowledge, a general interest story could benefit society even more. The author cannot lose sight that our work has two audiences, the science community, and more importantly, our family, friends, neighbors, and policy makers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A learned closure method applied to phase mixing in a turbulent gradient-driven gyrokinetic system in simple geometry

We present a new method for formulating closures that learn from kinetic simulation data. We apply this method to phase mixing in a simple gyrokinetic turbulent system – temperature-gradient-driven turbulence in an unsheared slab. The closure, called the learned multi-mode (LMM) closure, is constructed by, first, extracting an optimal basis from a nonlinear kinetic simulation using singular value decomposition. Subsequent nonlinear fluid simulations are projected onto this basis and the results are used to formulate the closure. We compare the closure with other closures schemes over a broad range of the relevant two-dimensional parameter space (collisionality and gradient drive). We find that the turbulent kinetic system produces phase-mixing rates much lower than the linear expectations, which the LMM closure is capable of capturing. We also compare radial heat fluxes. A Hammett–Perkins closure, generalized to include collisional effects, is quite successful throughout the parameter space, producing ${\sim }14\,\%$ root-mean-square (r.m.s.) error. The LMM closure is also very effective: when trained at three (two) points (in a 35 point parameter grid), the LMM closure produces $8\,\%$ ( $12\,\%$ ) r.m.s. errors. The LMM procedure can be readily generalized to other closure problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

On the Thermodynamic Condition for Adsorption Azeotropes

Adsorption azeotropy is a common phenomenon in mixed-gas adsorption equilibria. Since the first literature report of binary adsorption azeotropes in 1933, numerous studies investigated the thermodynamic conditions for adsorption azeotropes but failed to reach definitive conclusions. Based on the generalized Langmuir isotherm model for multicomponent adsorption equilibria which takes into account the vacant site as part of the adsorbed phase, this study presents the thermodynamic condition for adsorption azeotropes as derived from the generalized Langmuir isotherm to be the equality of the ratios of adsorbed phase activity coefficient γi and adsorption equilibrium constant $K$$^{o}_{i}$ for the two adsorbates in the binary 1–2 adsorption system, i.e., γ 1 /$K$$^{o}_{1}$ = γ 2 /$K$$^{o}_{2}$. This adsorption azeotropic condition is analogous to the vapor–liquid equilibrium azeotropic condition, i.e., the equality of the products of liquid phase activity coefficient and saturation vapor pressure Pisat for the two components, i.e., γ 1 $P$$^{sat}_{1}$ = γ 2 $P$$^{sat}_{2}$. We validated the thermodynamic condition for adsorption azeotropes with 14 azeotrope-forming adsorption systems in this study. As a result, we investigated the effects of the pressure, temperature, and adsorbed phase nonideality on azeotrope formation.

42 ENGINEERING↗
DOE OSTI2026

Signatures of Thermoreversible Associations in X-ray and Neutron Scattering from Dilute Polyzwitterion Solutions

In aqueous solutions of polyzwitterions (PZs), an interplay between dipole–dipole interactions and hydration of zwitterionic groups can lead to thermoreversible associations, which have been difficult to detect in experiments. Here, in this study, we investigated dilute aqueous solutions of poly(1-(3-sulphonatopropyl)-2-vinylpyridinium) (P2VPPS) using small-angle X-ray and neutron scattering (SAXS and SANS) to probe the structure and neutron spin-echo (NSE) spectroscopy to probe dynamics. The SAXS and SANS data show that the correlation length increases with an increase in the concentration of P2VPPS for three different molecular weights. The addition of 0.1 M NaCl to one of the solutions led to almost no dependence of the correlation length on the concentration. Such a concentration dependence of the correlation length suggests the formation of clusters driven by thermoreversible associations in the solutions. The NSE measurements show that the solutions with a larger correlation length display slower relaxation, reflecting the reduced mobility of larger clusters. These results should be considered as signatures of thermoreversible associations in the solutions of P2VPPS. To establish a quantitative link between local structure (clusters) and dynamics in dilute solutions of P2VPPS, we combined a thermoreversible gelation theory for the structure of PZ solutions (Li, S.-F.; Muthukumar, M. Theory of Thermoreversible Gelation and Anomalous Concentration Fluctuations in Polyzwitterion Solutions. J. Chem. Phys. 2024, 161, 024903) with a model for the dynamics of the clusters (generalized Zimm model), developed in this work. Using such a theoretical framework, we have predicted the distribution of clusters in the solutions probed with SANS and SAXS. With the distributions, the generalized Zimm model has been used to extract the diffusion constant of the clusters and their characteristic size from the NSE data, where the latter agrees with the values estimated from the SAXS/SANS data. These findings confirm the presence of thermoreversible associations and establish a quantitative link between local structure (clusters) and dynamics in dilute solutions of P2VPPS. Furthermore, with growing interest in technological applications, this work can provide useful insights into the structural and dynamical properties of other PZ solutions.

field theory↗
DOE OSTI2024

Emulator of PR‐DNS: Accelerating Dynamical Fields With Neural Operators in Particle‐Resolved Direct Numerical Simulation

Abstract Particle‐resolved direct numerical simulations (PR‐DNS) play an increasing role in investigating aerosol‐cloud‐turbulence interactions at the most fundamental level of processes. However, the high computational cost associated with high resolution simulations poses considerable challenges for large domain or long duration simulation using PR‐DNS. To address these issues, here we present an emulator of the complex physics‐based PR‐DNS developed by use of the data‐driven Fourier Neural Operator (FNO) method. The effectiveness of the method is showcased by presenting turbulence and temperature fields in a two‐dimensional space. The results demonstrate high accuracy at various resolutions and the emulator is two orders of magnitude cheaper in terms of computational demand compared to the physics‐based PR‐DNS model. Furthermore, the FNO emulator exhibits strong generalization capabilities for different initial conditions and ultra‐high‐resolution without the need to retrain models. These findings highlight the potential of the FNO method as a promising tool to simulate complex fluid dynamics problems with high accuracy, computational efficiency, and generalization capabilities, enhancing our understanding of the aerosol‐cloud‐precipitation system.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

New quantum physics, solving puzzles of Wheeler’s delayed choice and a particle’s passing N slits simultaneously and quantum oscillator in experiments

Abstract This paper discovers new quantum physics, and gives solutions to puzzles of Wheeler’s delayed choice and a particle’s passing many slits simultaneously by exact quantum physics expressions. We further show new quantum control, new quantum oscillation, new quantum control experiments and new quantum oscillator being able to be installed in quantum communication network etc. We discover that the ability of a photon to hit electrons out in photoelectric effect is complementarily equivalent to the ability of wave of a photon to simultaneously pass through many slits in wave-particle duality. Objective criterion for distinguishing classical and quantum particles is found, and this paper gives applicable realm of quantum theories and new quantum physics expressions of wave-particle duality. All these studies above should be classified as classical and quantum particles, then classical particle and quantum particle wave cannot and can pass many slits, respectively. This paper discovers wave-particle duality’s origin of displaying both wave property from plane wave part of the general Fourier expansion and particle property from the general Fourier expansion coefficients with the particle’s global property and spins etc. We give the superposition state representation of wave-particle duality, further find the collapse of the duality superposition state to wave or particle state. The collapsed wave or particle state is related to the measure of wave or particle property. Then, we explain why sometimes it's a wave or a particle. Our achieved results are truly tested, and we discover new measured attractive state and quantum wave collapse velocity expression.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

How well do self-interaction corrections repair the overestimation of static polarizabilities in density functional calculations?

Herein, we examine the effect of removing self-interaction error (SIE) on the calculation of molecular polarizabilities in the local spin density (LSDA) and generalized gradient approximations (GGA). To this end, we utilize a database of 132 molecules taken from a recent benchmark study [Hait and Head-Gordon, Phys. Chem. Chem. Phys. , 2018, 20 , 19800] to assess the influence of SIE on polarizabilities by comparing results with accurate reference data. Our results confirm that the general overestimation of molecular polarizabilities by these density functional approximations can be attributed to SIE. However, removing SIE using the Perdew–Zunger self-interaction-correction (PZ-SIC) method, implemented using the Fermi–Löwdin Orbital SIC approach, leads to an underestimation of molecular polarizabilities, showing that PZ-SIC overcorrects when combined with LSDA or GGA. Application of a recently proposed locally scaled SIC [Zope, et al., J. Chem. Phys. , 2019, 151 , 214108] is found to provide more accurate polarizabilities. We attribute this to the ability of the local scaling scheme to selectively correct for SIE in the regions of space where the correction is needed most.

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

Adsorption of ethylenediamine on Cu surfaces: attributes of a successful capping molecule using first-principles calculations

The shape-controlled synthesis of Cu nanocrystals can benefit a wide range of applications, though challenges exist in achieving high and selective yields to a particular shape. Capping agents play a pivotal role in controlling shape, but their exact role remains ambiguous. In this study, the adsorption of ethylenediamine (EDA) on Cu(100) and Cu(111) was investigated with quantum density functional theory (DFT) to reveal the complex roles of EDA in promoting penta-twinned Cu nanowire growth. We find EDA has stronger binding on Cu(100) than on Cu(111), which agrees the general expectation that penta-twinned Cu nanowires express facets with stronger capping-molecule binding. Despite this stronger binding, ab initio thermodynamics reveals the surface energy of EDA-covered Cu(111) is lower than that EDA-covered Cu(100) at all solution-phase EDA chemical potentials, so there is no thermodynamic driving force for penta-twinned nanowires. We also investigated the capability of EDA to protect Cu surfaces from oxidation in water by quantifying energy barriers for a water molecule to diffuse through EDA layers on Cu(100) and Cu(111). The energy barrier on Cu(100) is significantly lower, which supports observations of faster oxidation of Cu(100) in electrochemical experiments. Thus, we elucidate another possible function of a capping agent – to enable selective oxidation of crystal facets. This finding adds to the general understanding of successful attributes of capping agents for shape-selective nanocrystal growth.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Comparative analysis of atmospheric parameters from high-resolution spectroscopic sky surveys: APOGEE, GALAH, Gaia -ESO

SDSS-IV APOGEE-2, GALAH, and Gaia-ESO are high-resolution, ground-based, multi-object spectroscopic surveys providing fundamental stellar atmospheric parameters and multiple elemental abundance ratios for hundreds of thousands of stars of the Milky Way. Data from these and other surveys contribute to investigations of the history and evolution of the Galaxy. We undertake a comparison between the most recent data releases of these surveys to investigate the accuracy and precision of derived parameters by placing the abundances on an absolute scale. We also discuss the correlations in parameter and abundance differences as a function of main parameters. Uncovering the variants provides a basis to continue the efforts of future sky surveys. Quality samples from the APOGEE-GALAH (15 537 stars), APOGEE-GES (804 stars), and GALAH-GES (441 stars) overlapping catalogs were collected. We investigated the mean variants between the surveys, and linear trends were also investigated. We compared the slope of correlations and mean differences with the reported uncertainties. The average and scatter of v rad , T eff , log g, [M/H], and v micro , along with numerous species of elemental abundances in the combined catalogs, show that in general there is a good agreement between the surveys. We find large radial velocity scatters ranging from 1.3 km s -1 to 4.4 km s -1 when comparing the three surveys. We observe some weak trends (e.g., in ΔT eff vs. Δlog g for the APOGEE-GES stars) and a clear correlation in the v micro - Δv micro planes in the APOGEE-GALAH common sample. For [α/H], [Ti/H] (APOGEE-GALAH giants), and [Al/H] (APOGEE-GALAH dwarfs) potential strong correlations are discovered as a function of the differences in the main atmospheric parameters, and we also find weak trends for other elements. In general we find good agreement between the three surveys within their respective uncertainties. However, there are certain regimes in which strong variants exist, which we discuss. There are still offsets larger than 0.1 dex in the absolute abundance scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Galaxy cluster matter profiles - I. Self-similarity, mass calibration, and observable-mass relation validation employing cluster mass posteriors

We present a study of the weak lensing inferred matter profiles ΔΣ(R) of 698 South Pole Telescope (SPT) thermal Sunyaev-Zel’dovich effect (tSZE) selected and MCMF optically confirmed galaxy clusters in the redshift range 0.25 < z < 0.94 that have associated weak gravitational lensing shear profiles from the Dark Energy Survey (DES). Rescaling these profiles to account for the mass dependent size and the redshift dependent density produces average rescaled matter profiles ΔΣ(R/R200c)/(ρcritR200c) with a lower dispersion than the unscaled ΔΣ(R) versions, indicating a significant degree of self-similarity. Galaxy clusters from hydrodynamical simulations also exhibit matter profiles that suggest a high degree of self-similarity, with RMS variation among the average rescaled matter profiles with redshift and mass falling by a factor of approximately six and 23, respectively, compared to the unscaled average matter profiles. We employed this regularity in a new Bayesian method for weak lensing mass calibration that employs the so-called cluster mass posterior P(M200|ζ̂, λ̂, z), which describes the individual cluster masses given their tSZE (ζ̂) and optical (λ̂, z) observables. This method enables simultaneous constraints on richness λ-mass and tSZE detection significance ζ-mass relations using average rescaled cluster matter profiles. We validated the method using realistic mock datasets and present observable-mass relation constraints for the SPT×DES sample, where we constrained the amplitude, mass trend, redshift trend, and intrinsic scatter. Our observable-mass relation results are in agreement with the mass calibration derived from the recent cosmological analysis of the SPT×DES data based on a cluster-by-cluster lensing calibration. Our new mass calibration technique offers a higher efficiency when compared to the single cluster calibration technique. We present new validation tests of the observable-mass relation that indicate the underlying power-law form and scatter are adequate to describe the real cluster sample but that also suggest a redshift variation in the intrinsic scatter of the λ-mass relation may offer a better description. In addition, the average rescaled matter profiles offer high signal-to-noise ratio (S/N) constraints on the shape of real cluster matter profiles, which are in good agreement with available hydrodynamical ΛCDM simulations. This high S/N profile contains information about baryon feedback, the collisional nature of dark matter, and potential deviations from general relativity.Key words: gravitational lensing: weak / galaxies: clusters: general / large-scale structure of Universe

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A multifrequency study of sub-parsec jets with the Event Horizon Telescope

Context. The 2017 observing campaign of the Event Horizon Telescope (EHT) delivered the first very long baseline interferometry (VLBI) images at the observing frequency of 230 GHz, leading to a number of unique studies on black holes and relativistic jets from active galactic nuclei (AGN). In total, eighteen sources were observed, including the main science targets, Sgr A* and M 87, and various calibrators. Sixteen sources were AGN.Aims. We investigated the morphology of the sixteen AGN in the EHT 2017 data set, focusing on the properties of the VLBI cores: size, flux density, and brightness temperature. We studied their dependence on the observing frequency in order to compare it with the Blandford-Königl (BK) jet model. In particular, we aimed to study the signatures of jet acceleration and magnetic energy conversion.Methods. We modeled the source structure of seven AGN in the EHT 2017 data set using linearly polarized circular Gaussian components (1749+096, 1055+018, BL Lac, J0132–1654, J0006–0623, CTA 102, and 3C 454.3) and collected results for the other nine AGN from dedicated EHT publications, complemented by lower frequency data in the 2–86 GHz range. Combining these data into a multifrequency EHT+ data set, we studied the dependences of the VLBI core component flux density, size, and brightness temperature on the frequency measured in the AGN host frame (and hence on the distance from the central black hole), characterizing them with power law fits. We compared the observations with the BK jet model and estimated the magnetic field strength dependence on the distance from the central black hole.Results. Our observations spanning event horizon to parsec scales indicate a deviation from the standard BK model, particularly in the decrease of the brightness temperature with the observing frequency. Only some of the discrepancies may be alleviated by tweaking the model parameters or the jet collimation profile. Either bulk acceleration of the jet material, energy transfer from the magnetic field to the particles, or both are required to explain the observations. For our sample, we estimate a general radial dependence of the Doppler factor δ ∝ r≤0.5. This interpretation is consistent with a magnetically accelerated sub-parsec jet. We also estimate a steep decrease of the magnetic field strength with radius B ∝ r−3, hinting at jet acceleration or efficient magnetic energy dissipation.Key words: techniques: interferometric / galaxies: active / galaxies: jets / galaxies: nuclei / quasars: general / quasars: supermassive black holes⋆⋆ Deceased.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

A real-time search for Type Ia Supernovae with late-time interactions with circumstellar material in ZTF data

While it is generally accepted that Type Ia supernovae (SNe Ia) are the terminal explosions of white dwarfs (WDs), the nature of their progenitor systems and the mechanisms that lead up to these explosions remain widely debated. In rare cases, the SN ejecta interact with circumstellar material (CSM) that had previously been ejected from the progenitor system. The longer the delay between the creation of the CSM and the SN explosion, the greater the distance between the SN explosion site and the CSM and the later the onset of the interaction. The unknown distance between the CSM and SN explosion site makes it impossible to predict when the interaction will start. If the time between the SN explosion and the onset of the CSM interaction is of the order of several months to years, the SN has generally faded and it is no longer actively followed up on. This makes it even more difficult to detect the interaction while it is happening. In this work, we report on a real-time monitoring programme running between 13 November 2023 and 9 July 2024. It monitored 6914 SNe Ia for signs of late-time rebrightening using the Zwicky Transient Facility (ZTF). Flagged candidates were rapidly followed up on with photometry and spectroscopy to confirm the late-time excess and its position. We report the discovery of a ∼50 day rebrightening event in SN 2020qxz around 1200 rest-frame days after the peak of its light curve. SN 2020qxz exhibited signs of an early CSM interaction, but had faded from view over two years before its reappearance. Initial follow-up spectroscopy revealed the presence of four emission lines, while later follow-up spectroscopy showed that these had faded shortly after the end of the ZTF-detected rebrightening event. Our best match for these emission lines are H β (blueshifted by ∼5900 km s −1 ) and Ca II λ8542 , N I λ8567 , and K I λλ8763, 8767 (all blueshifted by 5100 km s −1 ; although we note that the line identifications are uncertain). This shows that catching and following up on late-time interactions as they occur can offer new clues on the nature of the progenitor systems that produce these SNe by putting constraints on the possible type of donor star. The only way to do this systematically is to use large sky surveys such as ZTF and the upcoming Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST) to monitor a large sample of objects for the rare events that reappear long after the object has faded from view.

circumstellar matter↗
DOE OSTI2021

Linear convergence of accelerated conditional gradient algorithms in spaces of measures

A class of generalized conditional gradient algorithms for the solution of optimization problem in spaces of Radon measures is presented. The method iteratively inserts additional Dirac-delta functions and optimizes the corresponding coefficients. Under general assumptions, a sub-linear [see formula in PDF] rate in the objective functional is obtained, which is sharp in most cases. To improve efficiency, one can fully resolve the finite-dimensional subproblems occurring in each iteration of the method. We provide an analysis for the resulting procedure: under a structural assumption on the optimal solution, a linear [see formula in PDF] convergence rate is obtained locally.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Re-evaluation of the TSL for Yttrium Hydride

Yttrium hydride (YH x ) is of interest as a high-temperature moderator material in advanced nuclear reactor systems because of its superior ability to retain hydrogen at elevated temperatures. Thermal neutron scattering laws (TSL) for hydrogen bound in yttrium hydride (H-YH 2 ) and yttrium bound in yttrium hydride (Y-YH 2 ) were previously evaluated by Naval Nuclear Laboratory using the ab initio approach and released in ENDF/B-VIII.0. In that work, density functional theory, incorporating the generalized gradient approximation (GGA) for the exchange-correlation energy, was used to simulate the face-centered cubic structure of YH 2 and calculate the interatomic Hellmann-Feynman forces for a 2×2×2 supercell containing 96 atoms. Lattice dynamics calculations using PHONON were used to determine the phonon density of states. The calculated phonon density of states for H and Y in YH 2 were then used to prepare H-YH 2 and Y-YH 2 TSL evaluations, in the incoherent approximation, using the LEAPR module of NJOY. In addition, elastic scattering was assumed to be incoherent for both H and Y. While the incoherent elastic scattering approximation is appropriate for H-YH 2 , it introduces an undesirable approximation for Y-YH 2 . In this work, we re-evaluate the TSL for Y-YH 2 using FLASSH (Full Law Analysis Scattering System Hub). Y-YH 2 is evaluated using the FLASSH generalized coherent elastic scattering capability in order to capture the Bragg peaks associated with the YH 2 crystal structure which were neglected in the prior NJOY-based evaluation due to limitations in LEAPR. An experimental approach to validate the Y-YH 2 TSL using neutron transmission measurements is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Analytical gradients for molecular-orbital-based machine learning

We report molecular-orbital-based machine learning (MOB-ML) enables the prediction of accurate correlation energies at the cost of obtaining molecular orbitals. Here, we present the derivation, implementation, and numerical demonstration of MOB-ML analytical nuclear gradients, which are formulated in a general Lagrangian framework to enforce orthogonality, localization, and Brillouin constraints on the molecular orbitals. The MOB-ML gradient framework is general with respect to the regression technique (e.g., Gaussian process regression or neural networks) and the MOB feature design. We show that MOB-ML gradients are highly accurate compared to other ML methods on the ISO17 dataset while only being trained on energies for hundreds of molecules compared to energies and gradients for hundreds of thousands of molecules for the other ML methods. The MOB-ML gradients are also shown to yield accurate optimized structures at a computational cost for the gradient evaluation that is comparable to a density-corrected density functional theory calculation.

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

A multi-dimensional Child–Langmuir law for any diode geometry

While prior theoretical studies of multi-dimensional space-charge limited current (SCLC) assumed emission from a small patch on infinite electrodes, none have considered emission from an entire finite electrode. In this paper, we apply variational calculus (VC) and conformal mapping, which have previously been used to derive analytic solutions for SCLC density (SCLCD) for nonplanar one-dimensional geometries, to obtain mathematical relationships for any multi-dimensional macroscopic diode with finite cathode and anode. We first derive a universal mathematical relationship between space-charge limited potential and vacuum potential for any diode and apply this technique to determine SCLCD for an eccentric spherical diode. We then apply VC and the Schwartz–Christoffel transformation to derive an exact equation for SCLCD in a general two-dimensional planar geometry with emission from a finite emitter. Particle-in-cell simulations using VSim agreed within 4%–13% for a range of ratios of emitter width to gap distance using the thinnest electrodes practical for the memory constraints of our hardware, with the difference partially attributed to the theory's assumption of infinitesimally thin electrodes. After generalizing this approach to determine SCLCD for any orthogonal diode as a function of only the vacuum capacitance and vacuum potential, we derive an analytical formulation of the three-dimensional Child–Langmuir law for finite parallel rectangular and disk geometries. These results demonstrate the utility for calculating SCLCD for any diode geometry using vacuum capacitance and vacuum potential, which are readily obtainable for many diode geometries, to guide experiment and simulation development.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Semi-empirical quantum optics for mid-infrared molecular nanophotonics

Nanoscale infrared (IR) resonators with sub-diffraction limited mode volumes and open geometries have emerged as new platforms for implementing cavity quantum electrodynamics at room temperature. The use of IR nanoantennas and tip nanoprobes to study strong light–matter coupling of molecular vibrations with the vacuum field can be exploited for IR quantum control with nanometer spatial and femtosecond temporal resolution. In order to advance the development of molecule-based quantum nanophotonics in the mid-IR, we propose a generally applicable semi-empirical methodology based on quantum optics to describe light–matter interaction in systems driven by mid-IR femtosecond laser pulses. The theory is shown to reproduce recent experiments on the acceleration of the vibrational relaxation rate in infrared nanostructures. It also provides physical insights on the implementation of coherent phase rotations of the near-field using broadband nanotips. We then apply the quantum framework to develop general tip-design rules for the experimental manipulation of vibrational strong coupling and Fano interference effects in open infrared resonators. We finally propose the possibility of transferring the natural anharmonicity of molecular vibrational levels to the resonator near-field in the weak coupling regime to implement intensity-dependent phase shifts of the coupled system response with strong pulses and develop a vibrational chirping model to understand the effect. The semi-empirical quantum theory is equivalent to first-principles techniques based on Maxwell’s equations, but its lower computational cost suggests its use as a rapid design tool for the development of strongly coupled infrared nanophotonic hardware for applications ranging from quantum control of materials to quantum information processing.

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

Ab initio calculations in atoms, molecules, and solids, treating spin–orbit coupling and electron interaction on an equal footing

We incorporate explicit, non-perturbative treatment of spin–orbit coupling into ab initio auxiliary-field quantum Monte Carlo (AFQMC) calculations. The approach allows a general computational framework for molecular and bulk systems in which material specificity, electron correlation, and spin–orbit coupling effects can be captured accurately and on an equal footing, with favorable computational scaling vs system size. We adopt relativistic effective-core potentials that have been obtained by fitting to fully relativistic data and that have demonstrated a high degree of reliability and transferability in molecular systems. This results in a two-component spin-coupled Hamiltonian, which is then treated by generalizing the ab initio AFQMC approach. We demonstrate the method by computing the electron affinity in Pb, the bond dissociation energy in Br 2 and I 2 , and solid Bi.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗