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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The approximate second order coupled-cluster method based on a size-consistent Brillouin–Wigner partitioning

We present a variant of the approximate second order coupled-cluster method (CC2) with a two-parameter size-consistent Brillouin–Wigner (BW-s) partitioning instead of a Møller–Plesset (MP) partitioning for the unperturbed Hamiltonian, which we refer to as BWs-CC2. The computational complexity of this model scales identically to CC2 with molecular size. Conventional CC2 and its regularized BWs-CC2 variants, as well as conventional MP2 and two of its regularized BW-s2 variants, were assessed on a 535 element database spanning thermochemistry, non-covalent interactions, barrier heights, and isomerization energies. To ensure a well-defined model chemistry, the assessment was performed using internally stable spin-polarized Hartree–Fock (HF) orbitals in the finite aug-cc-pVQZ basis without counterpoise corrections. As a result of using stable orbitals, contrary to conventional wisdom, we find that CC2 substantially outperforms MP2 on molecules with significantly spin contaminated reference orbitals without a significant increase in error on systems with a spin-pure reference, showing the value of its single substitutions. While no single choice of regularization parameters can be optimal for all datasets, we find that BWs-CC2 generally outperforms both CC2 and BW-s2 with a single judicious parameter choice. Additional tests on dipole moments and bond lengths of diatomics provide further support for the utility of this choice. Furthermore, the main outliers and poorest performing cases are associated with large amounts of spin-contamination in the HF reference, which is indicative of systems with either strong correlation or extensive artificial symmetry breaking. Overall, these findings argue that the perception of the quality of the CC2 ground state should be reevaluated and that it can be further improved upon by the soundly based BWs-CC2 variant with the recommended parameter choice.

Correlation energy↗
DOE OSTI2024

Fixed Source Sensitivity Calculations for Inertial Confinement Fusion Applications

A numerical code library was developed for the radiation transport code MCNP6.3 to calculate generalized response sensitivity coefficients for fixed source neutron transport problems with applications to inertial confinement fusion (ICF) experiments. The new MCNP6.3 dependency is used to generate a novel time convolution response that represents a neutron time-of-flight (nToF) signal. The traditional suite of macroscopic cross-section sensitivities and constrained fixed source probability distribution sensitivities are available for both the standard and the new response tallies in this library. However, novel sensitivity coefficients for the constrained hyperparameters of analytic fixed source probability distributions are emphasized in this work for their connection to ICF neutron transport models. Particularly, advanced Monte Carlo methods are developed for calculating the sensitivity of a nToF signal to perturbations in an ICF plasma’s ion temperature and burn history as well as perturbations in the target liner mass density and the shape parameters of the nToF detector’s impulse response function. Together, these capabilities form an advanced suite of computational tools that can be used to analyze and extract information from any ICF experimental platform.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Modeling Unconventional Superconductivity at the Crossover between Strong and Weak Electronic Interactions

High-temperature superconductivity emerges in many different quantum materials, often in regions of the phase diagram where the electronic kinetic energy is comparable to the electron-electron repulsion. Describing such intermediate-coupling regimes has proven challenging as standard perturbative approaches are inapplicable. In this study, we employ quantum Monte Carlo methods to solve a multiband Hubbard model that does not suffer from the sign problem and in which only repulsive interband interactions are present. In contrast to previous sign-problem-free studies, we treat magnetic, superconducting, and charge degrees of freedom on an equal footing. We find an antiferromagnetic dome accompanied by a metal-to-insulator crossover line in the intermediate-coupling regime, with a smaller superconducting dome appearing in the metallic region. Across the antiferromagnetic dome, the magnetic fluctuations change from overdamped in the metallic region to propagating in the insulating region. Our findings shed new light on the intertwining between superconductivity, magnetism, and charge correlations in quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

The inviscid incompressible limit of Kelvin–Helmholtz instability for plasmas

The Kelvin–Helmholtz Instability (KHI) is an interface instability that develops between two fluids or plasmas flowing with a common shear layer. KHI occurs in astrophysical jets, solar atmosphere, solar flows, cometary tails, planetary magnetospheres. Two applications of interest, encompassing both space and fusion applications, drive this study: KHI formation at the outer flanks of the Earth’s magnetosphere and KHI growth from non-uniform laser heating in magnetized direct-drive implosion experiments. Here, we study 2D KHI with or without a magnetic field parallel to the flow. We use both the GAMERA code, which solves the compressible Euler equations, and the STRATOSPEC code, which solves the Navier-Stokes equations under the Boussinesq approximation, coupled with the magnetic field dynamics. GAMERA is a global three-dimensional MHD code with high-order reconstruction in arbitrary nonorthogonal curvilinear coordinates, which is developed for a large range of astrophysical applications. STRATOSPEC is a three-dimensional pseudo-spectral code with an accuracy of infinite order (no numerical diffusion). Magnetized KHI is a canonical case for benchmarking hydrocode simulations with extended MHD options. An objective is to assess whether or not, and under which conditions, the incompressibility hypothesis allows to describe a dynamic compressible system. For comparing both codes, we reach the inviscid incompressible regime, by decreasing the Mach number in GAMERA, and viscosity and diffusion in STRATOSPEC. Here, we specifically investigate both single-mode and multi-mode initial perturbations, either with or without magnetic field parallel to the flow. The method relies on comparisons of the density fields, 1D profiles of physical quantities averaged along the flow direction, and scale-by-scale spectral densities. We also address the triggering, formation and damping of filamentary structures under varying Mach number or Atwood number, with or without a parallel magnetic field. Comparisons show very satisfactory results between the two codes. The vortices dynamics is well reproduced, along with the breaking or damping of small-scale structures. We end with the extraction of growth rates of magnetized KHI from the compressible regime to the incompressible limit in the linear regime assessing the effects of compressibility under increasing magnetic field. The observed differences between the two codes are explained either from diffusion or non-Boussinesq effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Polynomial-time preparation of low-temperature Gibbs states for two-dimensional toric code

In this work, we propose a polynomial-time algorithm for preparing the Gibbs state of the two-dimensional toric code Hamiltonian at any temperature, starting from any initial state, significantly improving upon prior estimates that suggested exponential scaling with inverse temperature. We prove that fast mixing at low temperature for the two-dimensional toric code can be achieved by augmenting local jump operators with simple global jump operators, which enable efficient transitions between logical sectors. To establish tight lower bounds on the spectral gap, we introduce a new reduction method that eventually maps the problem to estimating the spectral gap of a perturbed graph Laplacian on a stair graph. Our proof also shows that the Lindblad dynamics with a digitally implemented low-temperature local Davies generator is able to efficiently drive the quantum state toward the ground state manifold.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Hand Calculations for Nuclear Criticality Safety – Primer Revision

The document “ Hand Calculation Methods for Criticality Safety – A Primer, ” was published by the authors (Los Alamos National Laboratory report LA-14244-M) in December 2006 to provide an overview of the most common hand calculation methods used for criticality safety analysis. Hand calculation methods can be used as a starting point for more advanced calculations, and in many circumstances, they can provide sensitivity and perturbation information quicker than using a criticality code, such as KENO or MCNP. This report has been revised to augment existing basis information and example problems for each of the methods discussed: 1-group and modified 1- group diffusion theory, core-density conversions, buckling conversions, limiting surface density method, density analog method, surface density method, and solid angle method. This primer discusses the applicability of the various methods, illustrates how they are used, and provides an interpretation of the example problem results. The example problems provided are simple, practical problems, and the reference data to solve each example problem is provided. This paper will provide an overview of the Primer revision and how it can be useful for new Nuclear Criticality Safety practitioners.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Hand Calculations for Nuclear Criticality Safety — Primer Revision [Abstract]

The document “ Hand Calculation Methods for Criticality Safety –A Primer, ” was published by the authors (Los Alamos National Laboratory report LA-14244-M) in December 2006 to provide an overview of the most common hand calculation methods used for criticality safety analysis. Hand calculation methods can be used as a starting point for more advanced calculations, and in many circumstances, they can provide sensitivity and perturbation information quicker than using a criticality code, such as KENO or MCNP. This report has been revised to augment existing basis information and example problems for each of the methods discussed: 1-group and modified 1-group diffusion theory, core-density conversions, buckling conversions, limiting surface density method, density analog method, surface density method, and solid angle method. This primer discusses the applicability of the various methods, illustrates how they are used, and provides an interpretation of the example problem results. The example problems provided are simple, practical problems, and the reference data to solve each example problem is provided. This paper will provide an overview of the Primer revision and how it can be useful for new Nuclear Criticality Safety practitioners.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2020

Transverse momentum dependent PDFs at N3LO

We compute the quark and gluon transverse momentum dependent parton distribution functions at next-to-next-to-next-to-leading order (N 3 LO) in perturbative QCD. Our calculation is based on an expansion of the differential Drell-Yan and gluon fusion Higgs production cross sections about their collinear limit. This method allows us to employ cutting edge multiloop techniques for the computation of cross sections to extract these universal building blocks of the collinear limit of QCD. The corresponding perturbative matching kernels for all channels are expressed in terms of simple harmonic polylogarithms up to weight five. As a byproduct, we confirm a previous computation of the soft function for transverse momentum factorization at N 3 LO. Our results are the last missing ingredient to extend the q T subtraction methods to N 3 LO and to obtain resummed q T spectra at N 3 LL' accuracy both for gluon as well as for quark initiated processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Ultrafast Correlation Energy Estimator

A virtually no-cost method is proposed that can compute the correlation energies of general, covalently bonded, organic, and inorganic molecules (including conjugated π-electron systems) with a well-defined dominant Lewis structure at the accuracy of 99.5% of the near-exact values determined by the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] in the complete-basis-set (CBS) limit. This Correlation Energy Per Bond (CEPB) method assigns a partial correlation energy to each bond type (characterized by the identities of the two atoms forming the bond and its integer bond order) and to a lone pair, regardless of the bond length, bond angle, sp-hybridization, π-electron conjugation, ionicity, noncovalent interactions, etc. At its current stage, the method is mainly suitable for near-equilibrium geometries. The correlation energies per bond are determined by a fit to the CCSD(T)/CBS benchmarks. It can neither improve the equilibrium structures nor discern conformers or positional isomers, yet its accuracy for reaction energies rivals that of the second-order Møller–Plesset perturbation theory, which is far more expensive. Its promising performance underscores the possibility that surprisingly compact, chemically intuitive molecular fragments exist into which correlation energies can be partitioned, leading to various ultrafast correlation-energy estimators tailored to different purposes.

Binding energy↗
DOE OSTI2023

Vibrational heat-bath configuration interaction with semistochastic perturbation theory using harmonic oscillator or VSCF modals

Vibrational heat-bath configuration interaction (VHCI)—a selected configuration interaction technique for vibrational structure theory—has recently been developed in two independent works [J. H. Fetherolf and T. C. Berkelbach, J. Chem. Phys. 154, 074104 (2021); A. U. Bhatty and K. R. Brorsen, Mol. Phys. 119, e1936250 (2021)], where it was shown to provide accuracy on par with the most accurate vibrational structure methods with a low computational cost. Here, we eliminate the memory bottleneck of the second-order perturbation theory correction using the same (semi)stochastic approach developed previously for electronic structure theory. This allows us to treat, in an unbiased manner, much larger perturbative spaces, which are necessary for high accuracy in large systems. Stochastic errors are easily controlled to be less than 1 cm−1. We also report two other developments: (i) we propose a new heat-bath criterion and an associated exact implicit sorting algorithm for potential energy surfaces expressible as a sum of products of one-dimensional potentials; (ii) we formulate VHCI to use a vibrational self-consistent field (VSCF) reference, as opposed to the harmonic oscillator reference configuration used in previous reports. Our tests are done with quartic and sextic force fields, for which we find that with VSCF, the minor improvements to accuracy are outweighed by the higher computational cost associated the matrix element evaluations. We expect VSCF-based VHCI to be important for more general potential representations, for which the harmonic oscillator basis function integrals are no longer analytic.

Chemistry↗
DOE OSTI2023

A promising intersection of excited‐state‐specific methods from quantum chemistry and quantum Monte Carlo

Abstract We present a discussion of recent progress in excited‐state‐specific quantum chemistry and quantum Monte Carlo alongside a demonstration of how a combination of methods from these two fields can offer reliably accurate excited state predictions across singly excited, doubly excited, and charge transfer states. Both of these fields have seen important advances supporting excited state simulation in recent years, including the introduction of more effective excited‐state‐specific optimization methods, improved handling of complicated wave function forms, and ways of explicitly balancing the quality of wave functions for ground and excited states. To emphasize the promise that exists at this intersection, we provide demonstrations using a combination of excited‐state‐specific complete active space self‐consistent field theory, selected configuration interaction, and state‐specific variance minimization. These demonstrations show that combining excited‐state‐specific quantum chemistry and variational Monte Carlo can be more reliably accurate than either equation of motion coupled cluster theory or multi‐reference perturbation theory, and that it can offer new clarity in cases where existing high‐level methods do not agree. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Software > Quantum Chemistry

Otis, Leon↗
DOE OSTI2020

Assessing the sensitivity of land-atmosphere coupling strength to boundary and surface layer parameters in the WRF model over Amazon

Modeling tools can be used to diagnose regional land-atmosphere (L-A) coupling strength in the absence of sufficient observations, but subject to uncertainties associated with parameters in model physical parameterizations. Different sensitivity analysis (SA) approaches may lead to different conclusions about the underlying sensitivities. In this study, we quantify simulation uncertainties related to parameter perturbations, and use different approaches to conduct parameter SA on the WRF model pertaining to L-A coupling strength for simulations over the Amazon region. A total of twenty parameters from the Yonsei University (YSU) planetary boundary layer (PBL) and the revised MM5 surface layer (SL) schemes were selected in this analysis. Three different SA methods, the Morris One-at-A-Time (MOAT) method, the Multivariate Adaptive Regression Splines (MARS) method, and the Sobol’ method, were employed to analyze seven WRF-simulated variables and five L-A coupling metrics. Results show that 1) parameter perturbations cause large simulation uncertainties which are comparable to those in the observations; 2) three different SA methods give consistent L-A coupling strength outcomes; 3) six out of the twenty parameters contribute 80%–95% of the total variance in the metrics analyzed, and first-order effects dominate over interaction effects; 4) the twelve variables/metrics of interest show similar sensitivity patterns to the selected parameters, which is consistent across all the methods used. Physical mechanisms for how the sensitive parameters act in determining the L-A coupling strength and associated variables also are illustrated. Our results will help quantifying L-A coupling strength and establishing a basis for parameter calibration over the Amazon region.

Wang, Chen↗
DOE OSTI2026

A Space-Time Fluid (Unabridged)

Purpose: This essay is a retelling of general relativity in a language in which space-time geometry is expressed as a fluid. This trivial and useful reformulation gives 1) a non-perturbative covariant description of cosmological inhomogeneities and 2) a simple formula describing how cosmic inhomogeneities are generated on super-horizon scales. Methods: Equating the Ricci curvature with the associated matter stress-energy gives a description of space-time geometry in terms of fluid properties. These locally measurable (covariant) non-perturbative quantities are in some ways superior to commonly used "gauge invariant" quantities. The dynamics of a quantity (kurvature) which describes cosmological inhomogeneities is described in detail. A detailed comparison is made of space-time fluid dynamics with that of a classical (Newtonian physics) fluid. Results: The fluid lexicon permits an unambiguous definition of the velocity of space-time. The evolution of the space-time fluid is in many ways identical with that of the classical fluid when expressed in Lagrangian coordinates. Kurvature is a measure of the specific binding energy of the fluid and is a most useful covariant measure of cosmological inhomogeneities. For plausible matter models kurvature will increase, even on super-horizon scales, due to non-linear hydrodynamic effects rather than gravity. This phenomena is also exhibited by classical fluids. Conclusion: The space-time fluid representation of geometrodynamics gives a simple and useful description of the evolution of cosmological inhomogeneities.

Stebbins, Albert [Fermilab] (ORCID:000000023807725↗
DOE OSTI2021

An analysis of Bayesian estimates for missing higher orders in perturbative calculations

With current high precision collider data, the reliable estimation of theoretical uncertainties due to missing higher orders (MHOs) in perturbation theory has become a pressing issue for collider phenomenology. Traditionally, the size of the MHOs is estimated through scale variation, a simple but ad hoc method without probabilistic interpretation. Bayesian approaches provide a compelling alternative to estimate the size of the MHOs, but it is not clear how to interpret the perturbative scales, like the factorisation and renormalisation scales, in a Bayesian framework. Recently, it was proposed that the scales can be incorporated as hidden parameters into a Bayesian model. In this paper, we thoroughly scrutinise Bayesian approaches to MHO estimation and systematically study the performance of different models on an extensive set of high-order calculations. We extend the framework in two significant ways. First, we define a new model that allows for asymmetric probability distributions. Second, we introduce a prescription to incorporate information on perturbative scales without interpreting them as hidden model parameters. We clarify how the two scale prescriptions bias the result towards specific scale choice, and we discuss and compare different Bayesian MHO estimates among themselves and to the traditional scale variation approach. Finally, we provide a practical prescription of how existing perturbative results at the standard scale variation points can be converted to 68%/95% credibility intervals in the Bayesian approach using the new public code MiHO.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Identifying Differential Equations in Fourier Domain (FourierIdent)

We investigate identifying differential equations in the frequency domain. Fourier analysis is an important tool in theoretical analysis and numerical solvers of differential equations, yet there is limited work in exploring this connection in the identification of differential equations. This paper aims to identify the underlying differential equation in the frequency domain, from a given single realization of the differential equation perturbed by noise. Such setting imposes difficulties which are different from other identification methods where computation is carried out in the physical domain. We propose several ways to mitigate the challenges arising from noise in data and large differences in the magnitudes of frequency responses. The main takeaways are that identifying differential equations solely in the frequency domain is challenging, the method we propose is based on a form of domain partitions in the frequency domain, and this method shows benefits for complex data even with high level of noise. We introduce a Fourier feature denoising, and define the meaningful data region and the core regions of features to reduce the effect of noise in the frequency domain and to enhance the accuracy in coefficient identification. The proposed method is tested on various differential equations with linear, nonlinear, and high-order derivative feature terms, and shows advantages on complex data with many frequency modes, even under high level of noise.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Quantum adversarial learning for kernel methods

We show that hybrid quantum classifiers based on quantum kernel methods and support vector machines are vulnerable against adversarial attacks, namely small engineered perturbations of the input data can deceive the classifier into predicting the wrong result. Nonetheless, we also show that simple defense strategies based on data augmentation with a few crafted perturbations can make the classifier robust against new attacks. Our results find applications in security-critical learning problems and in mitigating the effect of some forms of quantum noise, since the attacker can also be understood as part of the surrounding environment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Ice Phase Classification Made Easy with Score-Based Denoising

Accurate identification of ice phases is essential for understanding various physicochemical phenomena. However, such classification for structures simulated with molecular dynamics is complicated by the complex symmetries of ice polymorphs and thermal fluctuations. For this purpose, both traditional order parameters and data-driven machine learning approaches have been employed, but they often rely on expert intuition, specific geometric information, or large training data sets. In this work, we present an unsupervised phase classification framework that combines a score-based denoiser model with a subsequent model-free classification method to accurately identify ice phases. Further, the denoiser model is trained on perturbed synthetic data of ideal reference structures, eliminating the need for large data sets and labeling efforts. The classification step utilizes the smooth overlap of atomic position (SOAP) descriptors as the atomic fingerprint, ensuring Euclidean symmetries and transferability to various structural systems. Our approach achieves a remarkable 100% accuracy in distinguishing ice phases of test trajectories using only seven ideal reference structures of ice phases as model inputs. This demonstrates the generalizability of the score-based denoiser model in facilitating phase identification for complex molecular systems. The proposed classification strategy can be broadly applied to investigate structural evolution and phase identification for a wide range of materials, offering new insights into the fundamental understanding of water and other complex systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗