Search NASA⌕ Search

SEARCH · Search NASA

Results for “FORMAL 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Classical eikonal from Magnus expansion

In a classical scattering problem, the classical eikonal is defined as the generator of the canonical transformation that maps in-states to out-states. It can be regarded as the classical limit of the log of the quantum S-matrix. In a classical analog of the Born approximation in quantum mechanics, the classical eikonal admits an expansion in oriented tree graphs, where oriented edges denote retarded/advanced worldline propagators. The Magnus expansion, which takes the log of a time-ordered exponential integral, offers an efficient method to compute the coefficients of the tree graphs to all orders. We exploit a Hopf algebra structure behind the Magnus expansion to develop a fast algorithm which can compute the tree coefficients up to the 12th order (over half a million trees) in less than an hour. In a relativistic setting, our methods can be applied to the post-Minkowskian (PM) expansion for gravitational binaries in the worldline formalism. We demonstrate the methods by computing the 3PM eikonal and find agreement with previous results based on amplitude methods. Importantly, the Magnus expansion yields a finite eikonal, while the naïve eikonal based on the time-symmetric propagator is infrared-divergent from 3PM on.

Black Holes↗
DOE OSTI2022

General, Rigorous Approach for the Treatment of Interfragment Covalent Bonds

Here, a generalized, projection-based transformation of the method-agnostic Fock operator in various ab initio fragment-based quantum chemistry methods has been developed for the treatment of interfragment covalent bonds. This transformation freezes the relevant localized molecular orbital associated with each interfragment bond, thereby restricting the variational subspace of the fragment wave functions, in order to maintain the proper physical characteristics of the involved covalent bonds. In addition, sets of orbitals that would lead to multiple occupancy of certain orbitals are explicitly removed from the variational space. The transformation is developed for the specific case of mutually orthonormal frozen and unfrozen orbitals within each fragment. The newly developed approach is then used to study model systems with two popular ab initio fragment-based methods, and the results of these calculations are compared to those obtained by existing methodologies. Analysis is focused on both quantitative and qualitative accuracy as well as computational scalability and stability. Other methods for which the developed formalisms are appropriate are outlined, and future extensions of the methods are discussed.

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

Training and projecting: A reduced basis method emulator for many-body physics

Here, we present the reduced basis method as a tool for developing emulators for equations with tun able parameters within the context of the nuclear many-body problem. The method uses a basis expansion informed by a set of solutions for a few values of the model parameters and then projects the equations over a well-chosen low-dimensional subspace. We connect some of the results in the eigenvector continuation literature to the formalism of reduced basis methods and show how these methods can be applied to a broad set of problems. As we illustrate, the possible success of the formalism on such problems can be diagnosed beforehand by a principal component analysis. We apply the reduced basis method to the one-dimensional Gross-Pitaevskii equation with a harmonic trap ping potential and to nuclear density functional theory for 48 Ca, achieving speed-ups of more than x150 and x250, respectively, when compared to traditional solvers. The outstanding performance of the approach, together with its straightforward implementation, show promise for its application to the emulation of computationally demanding calculations, including uncertainty quantification.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Comparison of second-order spectral line widths formulae

Frequently used second-order spectral line width formulae from the projector operator and kinetic theory methods have not been formally compared previously. Here it is shown that a systematic second-order expansion of the projection operator expression including initial correlations agrees with the second-order kinetic theory result. The agreement assumes a common approximation in the projector operator method that introduces a screened radiator–perturber interaction to account for neglected electron–electron correlations. On the other hand, it is shown that the usual width expression from the projection operator approach neglecting initial correlations differs from kinetic theory. The differences, however, are at least third order in the radiator–perturber interaction. The comparisons suggest using the more compact width expression from kinetic theory, which includes initial correlations and a systematic screening of the radiator–perturber interactions, as the starting point for second-order width calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Bringing discrete-time Langevin splitting methods into agreement with thermodynamics

In light of the recently published complete set of statistically correct Grønbech–Jensen (GJ) methods for discrete-time thermodynamics, we revise a differential operator splitting method for the Langevin equation in order to comply with the basic GJ thermodynamic sampling features, namely, the Boltzmann distribution and Einstein diffusion, in linear systems. This revision, which is based on the introduction of time scaling along with flexibility of a discrete-time velocity attenuation parameter, provides a direct link between the ABO splitting formalism and the GJ methods. This link brings about the conclusion that any GJ method has at least weak second order accuracy in the applied time step. It further helps identify a novel half-step velocity, which simultaneously produces both correct kinetic statistics and correct transport measures for any of the statistically sound GJ methods. Explicit algorithmic expressions are given for the integration of the new half-step velocity into the GJ set of methods. Finally, numerical simulations, including quantum-based molecular dynamics (QMD) using the QMD suite Los Alamos Transferable Tight-Binding for Energetics, highlight the discussed properties of the algorithms as well as exhibit the direct application of robust, time-step-independent stochastic integrators to QMD.

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

Lowering the Scaling of Self-Consistent Field Methods by Combining Tensor Hypercontraction and a Density Difference Ansatz

We present the tensor hypercontraction difference self-consistent field (SCF) method, an approach that reduces the formal computational scaling of traditional naive self-consistent field methods from 𝑂(𝑁 4 ) to 𝑂(𝑁 3 ) with system size 𝑁. The scaling reduction is achieved by developing a new technique for constructing the tensor hypercontraction decomposition based on the fundamental approximation made in density fitting. Combining this scheme with the difference self-consistent field methodology, we achieve a method that enables 𝑂(𝑁 3 ) scaling SCF calculations with only 𝑂(𝑁 2 ) storage requirements. In conclusion, our proof-of-concept numerical tests demonstrate robust performance with errors in total energies below 8 × 10 –4 E h and with sub 1 kcal/mol errors for relative energies.

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

Efficient general method for numerically modeling laser pulse propagation, overlap, and lifetime effects in amplifiers

An efficient numerical time-dependent general method is developed to address incoherent pulse overlap and lifetime effects in laser amplifiers. The alternating propagation-population laser energetics method (APPLE) has been validated against a semi-discrete coupled rate equation numerical method (SDRE) and analytic formalisms in bounding cases. APPLE is based on decoupled rates applied to a time-dependent framework where both space-time-dependent populations and pulse energetics are consistently updated in each time step. A significant advantage of APPLE lies in its conceptual simplicity, ease of implementation, and relatively small computational cost. SDRE tracks the populations through coupled rates and uses the method of lines to discretize the hyperbolic partial differential transport equations allowing for use of ordinary differential equation solvers. With reasonably sized mesh, we report both energetic and power pulse shape relative differences on the order of one percent between the models over a large range of initial conditions.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Giant Wilson loops and AdS 2 /dCFT 1

The 1/2-BPS Wilson loop in N = 4 supersymmetric Yang-Mills theory is an important and well-studied example of conformal defect. In particular, much work has been done for the correlation functions of operator insertions on the Wilson loop in the fundamental representation. In this paper, we extend such analyses to Wilson loops in the large-rank symmetric and antisymmetric representations, which correspond to probe D3 and D5 branes with AdS 2 × S 2 and AdS 2 × S 4 worldvolume geometries, ending at the AdS 5 boundary along a one-dimensional contour. We first compute the correlation functions of protected scalar insertions from supersymmetric localization, and obtain a representation in terms of multiple integrals that are similar to the eigenvalue integrals of the random matrix, but with important differences. Using ideas from the Fermi Gas formalism and the Clustering method, we evaluate their large N limit exactly as a function of the ’t Hooft coupling. The results are given by simple integrals of polynomials that resemble the Q-functions of the Quantum Spectral Curve, with integration measures depending on the number of insertions. Next, we study the correlation functions of fluctuations on the probe D3 and D5 branes in AdS. We compute a selection of three- and four-point functions from perturbation theory on the D-branes, and show that they agree with the results of localization when restricted to supersymmetric kinematics. We also explain how the difference of the internal geometries of the D3 and D5 branes manifests itself in the localization computation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Perspective on Many-Body Methods for Molecular Polaritonic Systems

Recent advances in strong light–matter interactions have revealed a wealth of new physical phenomena in molecules embedded in optical cavities, including modified chemical reactivity, altered excitation spectra, and novel quantum correlations. To describe these effects from first-principles, the field of ab initio quantum electrodynamics (QED) has emerged as a compelling extension of quantum chemistry that treats electronic and photonic degrees of freedom on equal footing. In this Perspective, we review the growing landscape of many-body QED methods, including Hartree–Fock, density functional theory (QEDFT), time-dependent DFT (QED-TDDFT), configuration interaction (QED-CI), complete active space (QED-CASSCF), coupled cluster (QED-CC), quantum Monte Carlo (QED-QMC), and density matrix renormalization group (QED-DMRG), highlighting recent developments and implementations. We further explore real-time methods, gradient and Hessian formalisms, and the integration of nonadiabatic nuclear dynamics. Applications range from benchmark simulations of polaritonic chemistry to quantum simulations on emerging quantum hardware. We conclude by outlining future directions for theory development and interdisciplinary efforts at the interface of quantum chemistry, condensed matter, and quantum optics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Numerically exact configuration interaction at quadrillion-determinant scale

The combinatorial growth of configuration interaction (CI) has long limited this formally exact quantum chemistry method to only the smallest molecules. Here, we report a numerically exact CI calculation exceeding one quadrillion (10 15 ) determinants, made possible by a lossless categorical compression strategy within the small-tensor-product distributed active space (STP-DAS) framework. This approach overcomes the traditional memory bottlenecks of CI by a numerically exact compression of the wavefunction representation and reformulating the most computationally demanding matrix–vector operations. Using this method, we performed a fully relativistic CI calculation of the ground state of HBrTe with over 10 15 complex-valued determinants in just 34.5 h on 1000 computing nodes—the largest CI calculation ever reported. We further achieved fast computation for systems with hundreds of billions of determinants on only a few compute nodes. Extensive benchmarks confirm that the method retains full numerical exactness while cutting memory and computational cost by orders of magnitude. Compared to previous state-of-the-art CI calculations, this work achieves a 1000 times increase in CI space, a 10 6 -fold increase in floating-point operations performed, and a 10 6 -fold improvement in computational speed.

Computational chemistry↗
DOE OSTI2021

A Bayesian analysis of classical shadows

Abstract The method of classical shadows proposed by Huang, Kueng, and Preskill heralds remarkable opportunities for quantum estimation with limited measurements. Yet its relationship to established quantum tomographic approaches, particularly those based on likelihood models, remains unclear. In this article, we investigate classical shadows through the lens of Bayesian mean estimation (BME). In direct tests on numerical data, BME is found to attain significantly lower error on average, but classical shadows prove remarkably more accurate in specific situations—such as high-fidelity ground truth states—which are improbable in a fully uniform Hilbert space. We then introduce an observable-oriented pseudo-likelihood that successfully emulates the dimension-independence and state-specific optimality of classical shadows, but within a Bayesian framework that ensures only physical states. Our research reveals how classical shadows effect important departures from conventional thinking in quantum state estimation, as well as the utility of Bayesian methods for uncovering and formalizing statistical assumptions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Exact restoration of Galilei invariance in density functional calculations with quantum Monte Carlo

We investigate some aspects of the connection between mean-field and ab-initio calculations for any-nucleon systems employing the Gogny effective interaction. In particular, we present variational Monte Carlo calculations of 4 He and 16 O, for which we extended the formalism of the linear method to the case of a density-dependent Hamiltonian. Monte Carlo calculations using an uncorrelated basis reduce to Hartree-Fock ones, but Galilean invariance can be preserved, contrary to the usual treatment. An analysis of the results obtained with the Gogny force when used in conjunction with correlated wave functions has also been performed. Here, the calculations show the collapse of the deuteron, leading to a set of unphysical predictions for the structure of heavier nuclei. This behavior remarks the different role of the Gogny effective interaction as compared to bare interactions typically used in fully microscopic calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Cosmology at high redshift — a probe of fundamental physics

An observational program focused on the high redshift (2<z<6) Universe has the opportunity to dramatically improve over upcoming LSS and CMB surveys on measurements of both the standard cosmological model and its extensions. Using a Fisher matrix formalism that builds upon recent advances in Lagrangian perturbation theory, we forecast constraints for future spectroscopic and 21-cm surveys on the standard cosmological model, curvature, neutrino mass, relativistic species, primordial features, primordial non-Gaussianity, dynamical dark energy, and gravitational slip. We compare these constraints with those achievable by current or near-future surveys such as DESI and Euclid, all under the same forecasting formalism, and compare our formalism with traditional linear methods. Here, our Python code FishLSS — used to calculate the Fisher information of the full shape power spectrum, CMB lensing, the cross-correlation of CMB lensing with galaxies, and combinations thereof — is publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Fast, differentiable, and extensible big bang nucleosynthesis package

Here, we introduce light isotope nucleosynthesis with JAX (LINX), a new differentiable public big bang nucleosynthesis code designed for fast parameter estimation. By leveraging JAX, LINX achieves both speed and differentiability, enabling the use of Bayesian inference, including gradient-based methods. We discuss the formalism used in LINX for rapid primordial elemental abundance predictions and give examples of how LINX can be used. When combined with differentiable cosmic microwave background power spectrum emulators, LINX can be used for joint cosmic microwave background and big bang nucleosynthesis analyses without requiring extensive computational resources, including on personal hardware.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Predicting magnetic properties of single-molecule magnets from self-interaction-free density-functional theory (Final Report)

In this project we investigated electronic structure of an intermediate-sized copper-based molecule and magnetic and hyperfine properties of several small non-magnetic and magnetic molecules including transition-metal elements by applying self-interaction corrections to density functional theory (DFT). In a sufficient number of cases, DFT-calculated electronic structure and magnetic properties of single-molecule magnets qualitatively differ from corresponding experimental data. This is partly due to self-interacting electrons within the DFT formalism. Recently, an efficient method to correct the self-interactions was proposed, i.e., Fermi-Lowdin prbital (FLO) based self-interaction corrected (SIC) methodology, within DFT. Henceforth, this method is referred to as FLO-SIC method which exists in FLOSIC code. We used this FLO-SIC method for our studies of electronic structure and magnetic and hyperfine properties of small magnetic molecules and non-magnetic molecules. Our study will provide insight into predictions of magnetic properties of single-molecule magnets where self-interaction corrections play a critical role. There are two components of this project. In the first work, we studied the electronic structure of a planar mononuclear Cu-based molecule in two oxidation states, using DFT with the FLO-SIC method. We chose this system because it is small enough and it includes a transition metal element. We found that the standard FLO-SIC method takes too much compute time even for the small transition-metal molecule and so we slightly modified the method in order to expedite the process. In the dianionic state, we found that the FLO-SIC spin density agrees quantitatively with accurate quantum chemistry methods, while DFT spin density without self-interaction corrections are severely deviated from the quantum chemistry methods. We also showed that the energy gap between the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of the dianionic state is larger than that of the monoanionic state. This result is consistent with experimental data. In the second work, we investigated how the interaction between the electron spin and nuclear spin is affected by electron self-interactions within small non-magnetic and magnetic molecules. Such an interaction is called hyperfine interaction. For molecules without significant orbital angular momentum, the hyperfine interaction consists of Fermi contact and dipolar interaction terms. Since the Fermi contact term depends on electron spin density at the nuclear site, it would be highly affected by self-interaction corrections. Therefore, we calculated the hyperfine interaction for the small molecules using DFT with the slightly modified expedited FLO-SIC method which was obtained in the first work, and compared the results to experimental data and DFT calculations without self-interaction corrections. We found significant improvement of the Fermi contact term computed using the FLO-SIC method for small magnetic molecules. Overall, the first and second work provided positive outlook of application of the FLO-SIC method to magnetic molecules and systems including transition-metal elements.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mixed quantum–classical approach to model non-adiabatic electron–nuclear dynamics: Detailed balance and improved surface hopping method

We develop a density matrix formalism to describe coupled electron–nuclear dynamics. To this end, we introduce an effective Hamiltonian formalism that describes electronic transitions and small (quantum) nuclear fluctuations along a classical trajectory of the nuclei. Using this Hamiltonian, we derive equations of motion for the electronic occupation numbers and for the nuclear coordinates and momenta. We show that, in the limit, when the number of nuclear degrees of freedom coupled to a given electronic transition is sufficiently high (i.e., the strong decoherence limit), the equations of motion for the electronic occupation numbers become Markovian. Furthermore, the transition rates in these (rate) equations are asymmetric with respect to the lower-to-higher energy transitions and vice versa. In thermal equilibrium, such asymmetry corresponds to the detailed balance condition. We also study the equations for the electronic occupations in the non-Markovian regime and develop a surface hopping algorithm based on our formalism. To treat the decoherence effects, we introduce additional “virtual” nuclear wave packets whose interference with the “real” (physical) wave packets leads to the reduction in coupling between the electronic states (i.e., decoherence) as well as to the phase shifts that improve the accuracy of the numerical approach. Remarkably, the same phase shifts lead to the detailed balance condition in the strong decoherence limit.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Conservation laws in coupled cluster dynamics at finite temperature

We extend the finite-temperature Keldysh non-equilibrium coupled cluster theory (Keldysh-CC) [A. F. White and G. K.-L. Chan, J. Chem. Theory Comput. 15, 6137–6253 (2019)] to include a time-dependent orbital basis. When chosen to minimize the action, such a basis restores local and global conservation laws (Ehrenfest’s theorem) for all one-particle properties while remaining energy conserving for time independent Hamiltonians. We present the time-dependent Keldysh orbital-optimized coupled cluster doubles method in analogy with the formalism for zero-temperature dynamics, extended to finite temperatures through the time-dependent action on the Keldysh contour. To demonstrate the conservation property and understand the numerical performance of the method, we apply it to several problems of nonequilibrium finite-temperature dynamics: a 1D Hubbard model with a time-dependent Peierls phase, laser driving of molecular H2, driven dynamics in warm-dense silicon, and transport in the single impurity Anderson model.

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

Scattering Observables from One- and Two-body Densities: Formalism and Application to $\pmb \gamma $ ${}^3\hbox {He}$ Scattering

We introduce the transition-density formalism, an efficient and general method for calculating the interaction of external probes with light nuclei. One- and two-body transition densities that encode the nuclear structure of the target are evaluated once and stored. They are then convoluted with an interaction kernel to produce amplitudes, and hence observables. Here, by choosing different kernels, the same densities can be used for any reaction in which a probe interacts perturbatively with the target. The method therefore exploits the factorisation between nuclear structure and interaction kernel that occurs in such processes. We study in detail the convergence in the number of partial waves for matrix elements relevant in elastic Compton scattering on 3 He. The results are fully consistent with our previous calculations in Chiral Effective Field Theory. But the new approach is markedly more computationally efficient, which facilitates the inclusion of more partial-wave channels in the calculation. We also discuss the usefulness of the transition-density method for other nuclei and reactions. Calculations of elastic Compton scattering on heavier targets like 4 He are straightforward extensions of this study, since the same interaction kernels are used. And the generality of the formalism means that our 3 He densities can be used to evaluate any 3 He elastic-scattering observable with contributions from one- and two-body operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗