Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Perturbative methods”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Application of deep learning methods for beam size control during user operation at the Advanced Light Source

Past research at the Advanced Light Source (ALS) provided a proof-of-principle demonstration that deep learning methods could be effectively employed to compensate for the significant perturbations to the transverse electron beam size induced by user-controlled adjustments of the insertion devices. However, incorporating these methods into the ALS’ daily operations has faced notable challenges. The complexity of the system’s operational requirements and the significant upkeep demands has restricted their sustained application during user operation. Here, we introduce the development of a more robust neural network (NN)-based algorithm that utilizes a novel online fine-tuning approach and its systematic integration into the day-to-day machine operations. Our analysis emphasizes the process of NN model selection, demonstrates the superior performance of the NN-based method over traditional feedback methods, and examines the effectiveness and resilience of the new algorithm during user-operation scenarios. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Long-Range Interaction between Molecules with Electronic Degeneracy: Beyond the Born–Oppenheimer Approximation

Long-range intermolecular interaction plays an essential role in (ultra)cold collisions. For monomers with degenerate electronic states, an accurate description of the long-range interaction can be challenging due to the potential breakdown of the Born–Oppenheimer approximation. In this work, we developed a general method to construct long-range interaction potential energy surfaces (IPESs) between arbitrary molecules using the degenerate perturbation theory of Schrödinger. To further simplify the application of this method, a group-theoretical method is introduced to identify nonzero components of electrostatic properties of monomers, especially for the nonlinear molecules. To demonstrate this approach, we report the long-range IPESs of the O( 3 P)–OH( 2 Π) and OH( 2 Π)–OH( 2 Π) systems. The excellent agreement with ab initio calculations in large separations, even near electronic degeneracies, confirms the accuracy of our approach. The method is expected to lay the groundwork for the accurate investigation of (ultra)cold collisional dynamics involving open-shell species.

ab initio calculations↗
DOE OSTI2021

A diagrammatic approach for automatically deriving analytical gradients of tensor hyper-contracted electronic structure methods

In this work, we introduce a diagrammatic approach to facilitate the automatic derivation of analytical nuclear gradients for tensor hyper-contraction (THC) based electronic structure methods. The automatically derived gradients are guaranteed to have the same scaling in terms of both operation count and memory footprint as the underlying energy calculations, and the computation of a gradient is roughly three times as costly as the underlying energy. The new diagrammatic approach enables the first cubic scaling implementation of nuclear derivatives for THC tensors fitted in molecular orbital basis (MO-THC). Furthermore, application of this new approach to THC-MP2 analytical gradients leads to an implementation, which is at least four times faster than the previously reported, manually derived implementation. Finally, we apply the new approach to the 14 tensor contraction patterns appearing in the supporting subspace formulation of multireference perturbation theory, laying the foundation for developments of analytical nuclear gradients and nonadiabatic coupling vectors for multi-state CASPT2.

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

First principles thermal transport modeling in GaN and related materials

Gallium nitride is a wide bandgap material utilized in a variety of technologies, including high-power electronics and light-emitting diodes, partly due to its favorable thermal properties. This chapter describes modern first-principles-based modeling of phonons and lattice thermal conductivity (k) of GaN, III-nitrides and related materials. In particular, we describe the theoretical underpinnings of calculating phonon dispersions, intrinsic phonon interactions, and other lattice dynamical properties from quantum perturbation theory and density functional theory (DFT) methods. Description of how these methods are then coupled with the Peierls-Boltzmann transport (PBT) equation to determine phonon distributions and lifetimes relevant for thermal transport is given. These theoretical and numerical methods have demonstrated quantitative accuracy and predictive power for calculating pristine k and defect-limited k from first principles for a variety of materials. We present a review of the literature utilizing DFT-PBT methods to understand novel k behaviors in III-nitrides and related materials.

Lindsay, Lucas↗
DOE OSTI2021

Towards complete assignment of the infrared spectrum of the protonated water cluster H + (H 2 O) 21

The spectroscopic features of protonated water species in dilute acid solutions have been long sought after for understanding the microscopic behavior of the proton in water with gas-phase water clusters H + (H 2 O) n extensively studied as bottom-up model systems. We present a new protocol for the calculation of the infrared (IR) spectra of complex systems, which combines the fragment-based Coupled Cluster method and anharmonic vibrational quasi-degenerate perturbation theory, and demonstrate its accuracy towards the complete and accurate assignment of the IR spectrum of the H + (H 2 O) 21 cluster. The site-specific IR spectral signatures reveal two distinct structures for the internal and surface four-coordinated water molecules, which are ice-like and liquid-like, respectively. The effect of inter-molecular interaction between water molecules is addressed, and the vibrational resonance is found between the O-H stretching fundamental and the bending overtone of the nearest neighboring water molecule. The revelation of the spectral signature of the excess proton offers deeper insight into the nature of charge accommodation in the extended hydrogen-bonding network underpinning this aqueous cluster.

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

Collisional excitation of interstellar PN by H 2 : New interaction potential and scattering calculations

Here, rotational excitation of interstellar PN molecules induced by collisions with H 2 is investigated. We present the first ab initio four-dimensional potential energy surface (PES) for the PN–H 2 van der Waals system. The PES was obtained using an explicitly correlated coupled cluster approach with single, double, and perturbative triple excitations [CCSD(T)-F12b]. The method of interpolating moving least squares was used to construct an analytical PES from these data. The equilibrium structure of the complex was found to be linear, with H 2 aligned at the N end of the PN molecule, at an intermolecular separation of 4.2 Å. The corresponding well-depth is 224.3 cm -1 . The dissociation energies were found to be 40.19 cm -1 and 75.05 cm -1 for complexes of PN with ortho-H 2 and para-H 2 , respectively. Integral cross sections for rotational excitation in PN–H 2 collisions were calculated using the new PES and were found to be strongly dependent on the rotational level of the H 2 molecule. These new collisional data will be crucial to improve the estimation of PN abundance in the interstellar medium from observational spectra.

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

LATTE: open-source, high-performance traveltime computation, tomography and source location in acoustic and elastic media

Traveltime-based tomography and source location are fundamental approaches for imaging subsurface structures and understanding the spatiotemporal distribution of seismicity from local to global scales. We present an open-source, high-performance framework integrating eikonal equation solvers and adjoint-state theory for traveltime computation, velocity tomography, source location and joint tomography-location in 2-D/3-D acoustic and elastic media. We introduce novel regularization schemes based on total generalized p-variation, structural similarity and multitask machine learning to enhance the fidelity and interpretability of inverted models and source locations. Key features of our implementation also include the ability to leverage both absolute-difference and double-difference traveltime misfits for high-fidelity velocity tomography and source parameter estimation; support for traveltime computation and inversion in diverse 2-D/3-D scenarios with arbitrary source and receiver distributions; and a perturbation-based optimal step-size estimation method to reduce computational costs. In addition, our implementation employs shared-memory and distributed-memory parallelization to provide an efficient solution for traveltime computation, tomography, and source location. In conclusion, we validate the efficacy and accuracy of our approach through multiple synthetic data examples.

58 GEOSCIENCES↗
DOE OSTI2020

Two-particle correlation functions in cluster perturbation theory: Hubbard spin susceptibilities

Cluster perturbation theory (CPT) is a computationally economic method commonly used to estimate the momentum- and energy-resolved single-particle Green's function. It has been used extensively in direct comparisons with experiments that effectively measure the single-particle Green's function, e.g., angle-resolved photoemission spectroscopy. However, many experimental observables are given by two-particle correlation functions. CPT can be extended to compute two-particle correlation functions by approximately solving the Bethe-Salpeter equation. We implement this method and focus on the transverse spin susceptibility, measurable via inelastic neutron scattering or with optical probes of atomic gases in optical lattices. Finally, we benchmark the method with the one-dimensional Fermi-Hubbard model by comparing with known results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Tensor train continuous time solver for quantum impurity models

The simulation of strongly correlated quantum impurity models is a significant challenge in modern condensed matter physics that has multiple important applications. Thus far, the most successful methods for approaching this challenge involve Monte Carlo techniques that accurately and reliably sample perturbative expansions to any order. However, the cost of obtaining high precision through these methods is high. Recently, tensor train decomposition techniques have been developed as an alternative to Monte Carlo integration. In this study, we apply these techniques to the single-impurity Anderson model at equilibrium by calculating the systematic expansion in power of the hybridization of the impurity with the bath. Furthermore, we demonstrate the performance of the method in a paradigmatic application, examining the first-order phase transition on the infinite-dimensional Bethe lattice, which can be mapped to an impurity model through dynamical mean field theory. Our results indicate that using tensor train decomposition schemes allows the calculation of finite-temperature Green's functions and thermodynamic observables with unprecedented accuracy. The methodology holds promise for future applications to frustrated multiorbital systems, using a combination of partially summed series with other techniques pioneered in diagrammatic and continuous time quantum Monte Carlo.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Modest set of simulations of local-type primordial non-Gaussianity in the DESI era

A measurement of a primordial non-Gaussianity (PNG) signal through late- or early-Universe probes has the potential to transform our understanding of the physics of the primordial Universe. While large-scale structure observables in principle contain vital information, interpreting these measurements is challenging due to poorly understood astrophysical effects. Luckily, N-body simulations, such as the abacuspng set presented in this study, consisting of 9 boxes, each of size L box = 2 Gpc/h and particle mass of 1.01 × 10 10 M ⊙ /h, provide a viable path forward. As validation, we find good agreement between the simulations and our expectations from one-loop perturbation theory (PT) and the "separate universe" method for the matter bispectrum, matter power spectrum and the halo bias parameter associated with PNG, b φ . As a science application, here we investigate the link between halo assembly bias and b φ for halo properties known to play a vital role in accurately predicting galaxy clustering: concentration, shear (environment), and accretion rate. We find a strong response for all three parameters, suggesting that the connection between b φ and the assembly history of halos needs to be taken into account by future PNG analyses. We further perform the first study of the b φ parameter from fits to early DESI data of the luminous red galaxy (LRG) and quasistellar object (QSO) samples and comment on the effect on f NL constraints for the allowed galaxy-halo models (note that σ[f NL ]∝ $\frac{σ[b_φ]}{b_φ}$). We find that the error on f NL is 15, 8, 7 for the LRGs at z = 0.5 and z = 0.8 and QSOs at z = 1.4, respectively, suggesting that a thorough understanding of galaxy assembly bias is warranted so as to perform robust high-precision analysis of local-type PNG with future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Bending effects and optical properties of WSe 2 nanoribbons of topological phase

A WSe 2 monolayer of 1T' phase is a large band gap quantum spin Hall insulator, supporting dissipationless charge and spin transports through the topologically protected edge states. Here, in this work, we explore the nanoribbon forms of 1T' phase WSe 2 by first-principles density functional calculations and the many-body perturbation GW and Bethe-Salpeter equation method. We found that the 1T' WSe 2 nanoribbon can show topological edge states with a ribbon width of ~4-6 nm. Those edge bands show crossing through the Fermi level an odd number of times, with one kind of spin-polarization connecting the valence band continuum and conduction band continuum. The topological features of the edge bands hold even under small and medium bending in the nanoribbon, while large bending induces large band splitting, resulting in a topological switch-off in the edge bands. The semiconducting 1T' WSe2 nanoribbon shows a large tunability with bending in optical absorption spectra and exciton states. The lowest-energy exciton is changed from optically dark in the flat nanoribbon to bright in the bent nanoribbons. These properties in the 1T' WSe 2 nanoribbons suggest potential applications in controllable quantum electronics and exciton-based quantum information processes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Maximum Power Reference Tracking Algorithm for Power Curtailment of Photovoltaic Systems

This paper presents an algorithm for power curtailment of photovoltaic (PV) systems under fast solar irradiance intermittency. Based on the Perturb and Observe (P&O) technique, the method contains an adaptive gain that is compensated in real-time to account for moments of lower power availability. In addition, an accumulator is added to the calculation of the step size to reduce the overshoot caused by large irradiance swings. A testbed of a three-phase single-stage, 500 kVA PV system is developed on the OPAL-RT eMEGAsim real-time simulator. Field irradiance data and a regulation signal from PJM (RTO) are used to compare the performance of the proposed method with other techniques found in the literature. Results indicate an operation with smaller overshoot, less dc-link voltage oscillations, and improved power reference tracking capability.

Paduani, Victor↗
DOE OSTI2026

Stochastic tensor contraction for quantum chemistry

Many computational methods in ab initio quantum chemistry are formulated in terms of high-order tensor contractions, whose cost determines the size of system that can be studied. We introduce stochastic tensor contraction to perform such operations with greatly reduced cost, and present its application to the gold-standard quantum chemistry method, coupled cluster theory with up to perturbative triples. For total energy errors more stringent than chemical accuracy, we reduce the computational scaling to that of mean-field theory, while starting to approach the mean-field absolute cost, thereby challenging the existing cost-to-accuracy landscape. Benchmarks against state-of-the-art local correlation approximations further show that we achieve an order-of-magnitude improvement in both total computation time and error, with significantly reduced sensitivity to system dimensionality and electron delocalization. We conclude that stochastic tensor contraction is a powerful computational primitive to accelerate a wide range of quantum chemistry.

Chemical Physics (physics.chem-ph)↗
DOE OSTI2020

Monte Carlo Global QCD Analyses of the Pion Parton Distribution Functions

As the lightest hadron, the pion presents itself as a dichotomy. While being the pseudo Goldstone boson associated with chiral symmetry breaking, it is simultaneously regarded as the lightest pseudoscalar meson typically composed of constituent quark and antiquark. As the chiral symmetry shares the same global symmetry as QCD, studying the pion and its internal structure is key in understanding one of the four fundamental forces of nature, namely the quantum chromodynamic strong interaction. We study the pion's parton distribution functions (PDFs), which are universal quantities that describe the structure of the pion in terms of its constituent quarks, antiquarks, and gluons. Our goal is to discover what the available data reveal for these universal quantities, namely the PDFs. Through the use of factorization theorems and perturbative QCD, we use Monte Carlo (MC) methods to extract the PDFs from the available Drell-Yan (DY) and leading neutron (LN) data. While the DY process involves two hadrons colliding, the LN electroproduction involves an electron beam incident on a target nucleon. We can constrain well the pion PDFs at large momentum fraction using DY data, and at low momentum fraction using LN data. We also use threshold resummation in the DY process to gather predictable higher order terms associated with the soft gluon radiation because their contributions to the cross section are nontrivial. In this dissertation, we parametrize the pion PDFs and fit those parameters to the available data. We make use of Bayesian inference and use MC techniques to sample the parameter space. Three sets of results are presented. We first extract the pion PDFs by fitting the DY and LN data. Then, we include transverse momentum dependent DY data in conjunction with the DY and LN data to extract the pion PDFs. Finally, we apply various methods of threshold resummation to the DY cross section and extract pion PDFs. The pion PDFs presented here are the first global QCD analyses performed as well as the first MC extracted pion PDFs, which are at the forefront of the pion PDF community.

Barry, Patrick↗
DOE OSTI2021

Reduced scaling formulation of CASPT2 analytical gradients using the supporting subspace method

We present a reduced scaling and exact reformulation of state specific complete active space second-order perturbation (CASPT2) analytical gradients in terms of the MP2 and Fock derivatives using the supporting subspace method. This work follows naturally from the supporting subspace formulation of the CASPT2 energy in terms of the MP2 energy using dressed orbitals and Fock builds. For a given active space configuration, the terms corresponding to the MP2-gradient can be evaluated with O(N5) operations, while the rest of the calculations can be computed with O(N3) operations using Fock builds, Fock gradients, and linear algebra. When tensor-hyper-contraction is applied simultaneously, the computational cost can be further reduced to O(N4) for a fixed active space size. The new formulation enables efficient implementation of CASPT2 analytical gradients by leveraging the existing graphical processing unit (GPU)-based MP2 and Fock routines. We present benchmark results that demonstrate the accuracy and performance of the new method. Example applications of the new method in ab initio molecular dynamics simulation and constrained geometry optimization are given.

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

Sensitivity of the pseudo-global warming method under flood conditions: a case study from the northeastern US

Abstract. Intensified extreme precipitation and corresponding floods are the most relevant consequences of climate change over the northeastern US (NEUS). To evaluate the impacts of climate change or certain climate perturbations on future extreme weather events which are dynamically similar to historic analogs, the pseudo-global warming (PGW) method has been frequently employed; however, this method lacks precise definition and guidelines, thus limiting its application. More specifically, three key questions related to the application of the PGW method remain unanswered: at what spatial scale should climate perturbations be applied? Among the different meteorological variables available, which ones should be perturbed? And will PGW projections vary significantly when different perturbations are applied? To address these questions, we examine the sensitivity and robustness of conclusions drawn from the PGW method over the NEUS by conducting multiple PGW experiments with varied perturbation spatial scales and choice of perturbed meteorological variables. The results show that the projections of precipitation and other essential variables at the regional mean scale are consistent across the PGW simulations, with a relative difference of much less than 10 %; however, different perturbation modifications can cause significant displacements of the storm events being simulated. Several previously assumed advantages of modifying only the temperature at regional mean scale, such as the preservation of geostrophic balance, do not appear to hold. Also, for these experiments, we find the regional mean perturbation produces a positive precipitation bias because it ignores the land–ocean warming contrast, which is a robust regional response to global warming. Overall, PGW experiments with perturbations from temperature or the combination of temperature and wind at the grid point scale are both recommended, depending on the research questions. The first approach can isolate the spatially dependent thermodynamic impact, and the latter incorporates both the thermodynamic and dynamic impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Unified analysis of finite-size error for periodic Hartree-Fock and second order Møller-Plesset perturbation theory

Despite decades of practice, finite-size errors in many widely used electronic structure theories for periodic systems remain poorly understood. For periodic systems using a general Monkhorst-Pack grid, there has been no comprehensive and rigorous analysis of the finite-size error in the Hartree-Fock theory (HF) and the second order Møller-Plesset perturbation theory (MP2), which are the simplest wavefunction based method, and the simplest post-Hartree-Fock method, respectively. Such calculations can be viewed as a multi-dimensional integral discretized with certain trapezoidal rules. Due to the Coulomb singularity, the integrand has many points of discontinuity in general, and standard error analysis based on the Euler-Maclaurin formula gives overly pessimistic results. The lack of analytic understanding of finite-size errors also impedes the development of effective finite-size correction schemes. We propose a unified analysis to obtain sharp convergence rates of finite-size errors for the periodic HF and MP2 theories. Our main technical advancement is a generalization of the result of Lyness [Math. Comp. 30 (1976), pp. 1–23] for obtaining sharp convergence rates of the trapezoidal rule for a class of non-smooth integrands. Our result is applicable to three-dimensional bulk systems as well as low dimensional systems (such as nanowires and 2D materials). Our unified analysis also allows us to prove the effectiveness of the Madelung-constant correction to the Fock exchange energy, and the effectiveness of a recently proposed staggered mesh method for periodic MP2 calculations (see X. Xing, X. Li, and L. Lin [J. Chem. Theory Comput. 17 (2021), pp. 4733–4745]). In conclusion, our analysis connects the effectiveness of the staggered mesh method with integrands with removable singularities, and suggests a new staggered mesh method for reducing finite-size errors of periodic HF calculations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Assessing the Sensitivity of the Tropical Cyclone Boundary Layer to the Parameterization of Momentum Flux in the Community Earth System Model

Recent studies have demonstrated that high-resolution (~25 km) Earth System Models (ESMs) have the potential to skillfully predict tropical cyclone (TC) occurrence and intensity. However, biases in ESM TCs still exist, largely due to the need to parameterize processes such as boundary layer (PBL) turbulence. Building on past studies, we hypothesize that the depiction of the TC PBL in ESMs is sensitive to the configuration of the PBL parameterization scheme, and that the targeted perturbation of tunable parameters can reduce biases. The Morris one-at-a-time (MOAT) method is implemented to assess the sensitivity of the TC PBL to tunable parameters in the PBL scheme in an idealized configuration of the Community Atmosphere Model, version 6 (CAM6). The MOAT method objectively identifies several parameters in an experimental version of the Cloud Layers Unified by Binormals (CLUBB) scheme that appreciably influence the structure of the TC PBL. We then perturb the parameters identified by the MOAT method within a suite of CAM6 ensemble simulations and find a reduction in model biases compared to observations and a high-resolution, cloud-resolving model. Importantly, we demonstrate that the high-sensitivity parameters are tied to PBL processes that reduce turbulent mixing and effective eddy diffusivity, and that in CAM6 these parameters alter the TC PBL in a manner consistent with past modeling studies. In this way, we provide an initial identification of process-based input parameters that, when altered, have the potential to improve TC predictions by ESMs.

54 ENVIRONMENTAL SCIENCES↗