Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonperturbative 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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Nonperturbative methods in HZE ion transport

A nonperturbative analytic solution of the high charge and energy (HZE) Green's function is used to implement a computer code for laboratory ion beam transport. The code is established to operate on the Langley Research Center nuclear fragmentation model used in engineering applications. Computational procedures are established to generate linear energy transfer (LET) distributions for a specified ion beam and target for comparison with experimental measurements. The code is highly efficient and compares well with the perturbation approximations.

Wilson, John W.↗
NASA NTRS1994

Nonperturbative methods in HZE propagation

An analytical solution to the perturbative multiple collision series of a fragmenting HZE ion beam has limited usefulness since the first collision term has several hundred contributions, the second collision term has tens of thousands of contributions, and each successive collision term progresses to unwieldy computational proportions. Our previous work has revealed the multiple collision terms in the straight-ahead approximation to be simple products of a spatially dependent factor times a linear energy-dependent factor of limited domain and unit normalization. The properties of these forms allow the development of the nonperturbative summation of the series to all orders assuming energy-independent nuclear cross sections as matrix products of a scaled Green's function described herein. This nonperturbative Green's function with multiple scattering correction factors compares well with experiments using 670 MeV/u neon-20 ion beams in thick water targets.

Non-NASA Center↗
NASA NTRS1996

A Green's function method for high charge and energy ion transport

A heavy-ion transport code using Green's function methods is developed. The low-order perturbation terms exhibiting the greatest energy variation are used as dominant energy-dependent terms, and the higher order collision terms are evaluated using nonperturbative methods. The recently revised NUCFRG database is used to evaluate the solution for comparison with experimental data for 625A MeV 20Ne and 517A MeV 40Ar ion beams. Improved agreements with the attenuation characteristics for neon ions are found, and reasonable agreement is obtained for the transport of argon ions in water.

Energy Transfer↗
NASA NTRS1995

A Green's function method for heavy ion beam transport

The use of Green's function has played a fundamental role in transport calculations for high-charge high-energy (HZE) ions. Two recent developments have greatly advanced the practical aspects of implementation of these methods. The first was the formulation of a closed-form solution as a multiple fragmentation perturbation series. The second was the effective summation of the closed-form solution through nonperturbative techniques. The nonperturbative methods have been recently extended to an inhomogeneous, two-layer transport media to simulate the lead scattering foil present in the Lawrence Berkeley Laboratories (LBL) biomedical beam line used for cancer therapy. Such inhomogeneous codes are necessary for astronaut shielding in space. The transport codes utilize the Langley Research Center atomic and nuclear database. Transport code and database evaluation are performed by comparison with experiments performed at the LBL Bevalac facility using 670 A MeV 20Ne and 600 A MeV 56Fe ion beams. The comparison with a time-of-flight and delta E detector measurement for the 20Ne beam and the plastic nuclear track detectors for 56Fe show agreement up to 35%-40% in water and aluminium targets, respectively.

NASA Discipline Radiation Health↗
DOE OSTI2020

Facile All-Optical Method for In Situ Detection of Low Amounts of Ammonia

As a key precursor for nitrogenous compounds and fertilizer, ammonia affects our lives in numerous ways. Rapid and sensitive detection of ammonia is essential, both in environmental monitoring and in process control for industrial production. Here we report a novel and nonperturbative method that allows rapid detection of ammonia at low concentrations, based on the all-optical detection of surface-enhanced Raman signals. We show that this simple and affordable approach enables ammonia probing at selected regions of interest with high spatial resolution, making in situ and operando observations possible.

09 BIOMASS FUELS↗
DOE OSTI2021

Convergence of Eigenvector Continuation

Eigenvector continuation is a computational method that finds the extremal eigenvalues and eigenvectors of a Hamiltonian matrix with one or more control parameters. It does this by projection onto a subspace of eigenvectors corresponding to selected training values of the control parameters. The method has proven to be very efficient and accurate for interpolating and extrapolating eigenvectors. However, almost nothing is known about how the method converges, and its rapid convergence properties have remained mysterious. In this Letter, we present the first study of the convergence of eigenvector continuation. In order to perform the mathematical analysis, we introduce a new variant of eigenvector continuation that we call vector continuation. We first prove that eigenvector continuation and vector continuation have identical convergence properties and then analyze the convergence of vector continuation. Our analysis shows that, in general, eigenvector continuation converges more rapidly than perturbation theory. The faster convergence is achieved by eliminating a phenomenon that we call differential folding, the interference between nonorthogonal vectors appearing at different orders in perturbation theory. From our analysis we can predict how eigenvector continuation converges both inside and outside the radius of convergence of perturbation theory. Further, while eigenvector continuation is a nonperturbative method, we show that its rate of convergence can be deduced from power series expansions of the eigenvectors. Our results also yield new insights into the nature of divergences in perturbation theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The anomalous magnetic moment of the muon in the Standard Model

We review the present status of the Standard Model calculation of the anomalous magnetic moment of the muon. This is performed in a perturbative expansion in the fine-structure constant α and is broken down into pure QED, electroweak, and hadronic contributions. The pure QED contribution is by far the largest and has been evaluated up to and including O(α 5 ) with negligible numerical uncertainty. The electroweak contribution is suppressed by (m μ / M W ) 2 and only shows up at the level of the seventh significant digit. It has been evaluated up to two loops and is known to better than one percent. Hadronic contributions are the most difficult to calculate and are responsible for almost all of the theoretical uncertainty. The leading hadronic contribution appears at O(α 2 ) and is due to hadronic vacuum polarization, whereas at O(α 3 ) the hadronic light-by-light scattering contribution appears. Given the low characteristic scale of this observable, these contributions have to be calculated with nonperturbative methods, in particular, dispersion relations and the lattice approach to QCD. The largest part of this review is dedicated to a detailed account of recent efforts to improve the calculation of these two contributions with either a data-driven, dispersive approach, or a first- principle, lattice-QCD approach. The final result reads $a$$^{SM}_{μ}$ = 116 591 810(43) × 10 −11 and is smaller than the Brookhaven measurement by 3.7σ. The experimental uncertainty will soon be reduced by up to a factor four by the new experiment currently running at Fermilab, and also by the future J-PARC experiment. Furthermore, this and the prospects to further reduce the theoretical uncertainty in the near future – which are also discussed here – make this quantity one of the most promising places to look for evidence of new physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Multireference Equation-of-Motion Driven Similarity Renormalization Group: Theoretical Foundations and Applications to Ionized States

We present a formulation and implementation of an equation-of-motion (EOM) extension of the multireference driven similarity renormalization group (MR-DSRG) formalism for ionization potentials (IP-EOM-DSRG). The IP-EOM-DSRG formalism results in a Hermitian generalized eigenvalue problem, delivering accurate ionization potentials for strongly correlated systems. The EOM step scales as O(N 5 ) with the basis set size N, allowing for efficient calculation of spectroscopic properties, such as transition energies and intensities. The IP-EOM-DSRG formalism is combined with three truncation schemes of the parent MR-DSRG theory: an iterative nonperturbative method with up to two-body excitations [MR-LDSRG(2)] and second- and third-order perturbative approximations [DSRG-MRPT2/3]. We benchmark these variants by computing (1) the vertical valence ionization potentials of a series of small molecules at both equilibrium and stretched geometries; (2) the spectroscopic constants of several low-lying electronic states of the OH, CN, N 2 + , and CO + radicals; and (3) the binding curves of low-lying electronic states of the CN radical. A comparison with experimental data and theoretical results shows that all three IP-EOM-DSRG methods accurately reproduce the vertical ionization potentials and spectroscopic constants of these systems. Notably, the DSRG-MRPT3 and MR-LDSRG(2) versions outperform several state-of-the-art multireference methods of comparable or higher cost.

Hamiltonians↗
DOE OSTI2023

Stochastic and mixed density functional theory within the projector augmented wave formalism for simulation of warm dense matter

Stochastic density functional theory (DFT) and mixed stochastic-deterministic DFT are burgeoning approaches for the calculation of the equation of state and transport properties in materials under extreme conditions. In the intermediate warm dense matter regime, a state between correlated condensed matter and kinetic plasma, electrons can range from being highly localized around nuclei to delocalized over the whole simulation cell. The plane-wave basis pseudopotential approach is thus the typical tool of choice for modeling such systems at the DFT level. Unfortunately, stochastic DFT methods scale as the square of the maximum plane-wave energy in this basis. To reduce the effect of this scaling and improve the overall description of the electrons within the pseudopotential approximation, we present stochastic and mixed DFT approaches developed and implemented within the projector augmented wave formalism. In conclusion, we compare results between the different DFT approaches for both single-point and molecular dynamics trajectories and present calculations of self-diffusion coefficients of solid density carbon from 1 to 50 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Large-momentum effective theory

Since the parton model was introduced by Feynman more than 50 years ago, much has been learned about the partonic structure of the proton through a large body of high-energy experimental data and dedicated global fits. However, limited progress has been made in calculating partonic observables such as the parton distribution function (PDFs) from the fundamental theory of strong interactions, quantum chromodynamics (QCD). Recently some advocated for a formalism, large-momentum effective theory (LaMET), through which one can extract parton physics from the properties of the proton traveling at a moderate boost factor such as γ ~ 2 – 5 . The key observation behind this approach is that Lorentz symmetry allows the standard formalism of partons in terms of light-front operators to be replaced by an equivalent one with large-momentum states and time-independent operators of a universality class. With LaMET, the PDFs, generalized PDFs or generalized parton distributions, transverse-momentum-dependent PDFs, and light-front wave functions can all be extracted in principle from lattice simulations of QCD (or other nonperturbative methods) through standard effective field theory matching and running. Future lattice QCD calculations with exascale computational facilities could help one to understand the experimental data related to the hadronic structure, including those from the upcoming electron-ion colliders dedicated to exploring the partonic landscape of the proton. In this work, the progress made in the past few years in the development of the LaMET formalism and its applications is reviewed, with an emphasis on a demonstration of its effectiveness from initial lattice QCD simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1980

A multiphoton treatment of near resonant light scattering in intense fields

The scattering by a two-state atom of monochromatic light of arbitrary intensity and detuning, is treated by solving Schroedinger's equation directly for the state of the scattered field. The calculation is performed in a basis of dressed atom states, which include to all orders the coupling of the incident field to the atom. Amplitudes describing sequences of scattering events are found by a nonperturbative method. A scattered intensity spectrum, directly related to the rate of change of energy in any given scattered field mode, is constructed in terms of bilinear products of the amplitudes representing all possible sequences of photon emission, and in the steady state limit the spectrum is found to be identical to the widely quoted results of Mollow. A discussion is devoted to understanding in physical terms the various quantities appearing in the expression for the intensity spectrum.

Ballagh, R. J.↗
DOE OSTI2024

Effective many-body interactions in reduced-dimensionality spaces through neural network models

Accurately describing properties of challenging problems in physical sciences often requires complex mathematical models that are unmanageable to tackle head on. Therefore, developing reduced-dimensionality representations that encapsulate complex correlation effects in many-body systems is crucial to advance the understanding of these complicated problems. However, a numerical evaluation of these predictive models can still be associated with a significant computational overhead. To address this challenge, in this paper we discuss a combined framework that integrates recent advances in the development of active-space representations of coupled cluster (CC) downfolded Hamiltonians with neural network approaches. The primary objective of this effort is to train neural networks to eliminate the computationally expensive steps required for evaluating hundreds or thousands of Hugenholtz diagrams, which correspond to multidimensional tensor contractions necessary for evaluating a many-body form of downfolded effective Hamiltonians. Using small molecular systems (the H 2 O and HF molecules) as examples, we demonstrate that training neural networks employing effective Hamiltonians for a few nuclear geometries of molecules can accurately interpolate or extrapolate their forms to other geometrical configurations characterized by different intensities of correlation effects. We also discuss differences between effective interactions that define CC downfolded Hamiltonians with those of bare Hamiltonians defined by Coulomb interactions in the active spaces. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Effect of the orientation of Rydberg atoms on their collisional ionization cross section

Collisional ionization between two Rydberg atoms in relative motion is examined. A classical trajectory Monte Carlo method is used to determine the cross sections associated with Penning ionization. The dependence of the ionization cross section on the magnitude and the direction of orbital angular momentum of the electrons and the direction of the Laplace-Runge-Lenz vector of the electrons is studied. For a given magnitude of angular momentum, there can exist a difference of a factor of up to ~2.5 in the ionization cross section between the orientations with the highest and the lowest ionization cross section. As a result, the case of exchange ionization is examined and its dependence on the magnitude of angular momentum is studied.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Strong-field nonsequential double photoionization using virtual-detector theory with path summation

We present an ab initio study of the nonsequential strong-field ionization dynamics of a model two-electron atom with helium character. Single- and double-ionization events are characterized and displayed using detector signals extracted at different points in the two-electron two-dimensional space. The double photoelectron momentum distribution is calculated via coherent path-summation over virtual-particle trajectories. A comparison is made between the momentum distributions obtained with the virtual-detector method, the Schrödinger wave function, and the time-dependent surface flux method developed by Tao and Scrinzi [New J. Phys. 14, 013021 (2012)]. Insights into different ionization and electron recollision pathways are gained from detailed virtual-particle tracking and energy-time readouts. Furthermore, this study demonstrates the extension of virtual-detector theory to strong-field multielectron quantum dynamics and highlights the importance of the evolving quantum phase in quasiclassical electron propagation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Interference in nonlinear Compton scattering using a Schrödinger-equation approach

The interference between Compton scattering and nonlinear Compton scattering from a two-color field in the X-ray regime is theoretically examined for bound electrons. Here, the underlying phase shifts are analyzed using a perturbative approach in the incoming classical field. The perturbative approach is bench marked with a non-perturbative approach in the classical field. The interference for different combinations of linear polarization of the two fields is examined when the Compton and the nonlinear Compton scattered waves have the same wave vector and polarization. Only two cases exhibit interference. When there is interference, the calculations reveal an intrinsic phase difference between the Compton scattered wave function and the nonlinear Compton scattered wave function of either 0 or pi depending on the scattering angle.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Exact time-dependent density-functional theory for nonperturbative dynamics of the helium atom

By inverting the time-dependent Kohn-Sham equation for a numerically exact dynamics of the helium atom, we show that the dynamical step and peak features of the exact correlation potential found previously in one-dimensional models persist for real three-dimensional systems. We demonstrate that the Kohn-Sham and true current-densities differ by a rotational component. The results have direct implications for approximate TDDFT calculations of atoms and molecules in strong fields, emphasizing the need to go beyond the adiabatic approximation, and highlighting caution in quantitative use of the Kohn-Sham current.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Simulations for x-ray imaging of wave-packet dynamics

Previous work on imaging wave-packet dynamics with x-ray scattering revealed that the scattering patterns deviate substantially from the notion of instantaneous momentum density of the wave packet. Here we show that scattering patterns can provide clear insights into the electron wave-packet dynamics if the final state of the scattered electron and the scattered photon momentum are determined simultaneously. The scattering probability is shown to be proportional to the modulus square of the Fourier transform of the instantaneous electronic spatial wave function weighted by the final state of the electron. Several cases for the choice of final state of the electron are explored. First, the case where the final state can be measured up to a given principal quantum number n and orbital angular momentum l are presented. Next, the case where the final states can only be determined up to a given energy is discussed. Lastly, the case of an initial wave packet consisting of a large amount of a known stationary state and a small amount of an unknown stationary state is examined. The scattering profile is used to determine the properties of the unknown state in the wave packet.

74 ATOMIC AND MOLECULAR PHYSICS↗