Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Greens function”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Coordinate-space solver for superfluid many-fermion systems with the shifted conjugate-orthogonal conjugate-gradient method

Self-consistent approaches to superfluid many-fermion systems in three dimensions (and their subsequent use in time-dependent studies) require a large number of diagonalizations of very large dimension Hermitian matrices, which results in enormous computational costs. We present an approach based on the shifted conjugate-orthogonal conjugate-gradient (COCG) Krylov method for the evaluation of the Green's function, from which we subsequently extract various densities (particle number, spin, current, kinetic energy, anomalous, etc.) of a nuclear system. The approach eschews the determination of the quasiparticle wave functions and their corresponding quasiparticle energies, which never explicitly appear in the construction of a single-particle Hamiltonian or needed for the calculation of various static nuclear properties, which depend only on densities. As benchmarks this study presents calculations for nuclei with axial symmetry, including the ground state of spherical (magic or semimagic) and axially deformed nuclei, the saddle point in the $^{240}Pu$ constrained fission path, and a vortex in the neutron star crust, and demonstrate the superior efficiency of the shifted COCG Krylov method over traditional approaches.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Quantum Monte Carlo calculations of light nuclei with local chiral two- and three-nucleon interactions

Local chiral effective field theory interactions have recently been developed and used in the context of quantum Monte Carlo few- and many-body methods for nuclear physics. In this paper, we go over detailed features of local chiral nucleon-nucleon interactions and examine their effect on properties of the deuteron, paying special attention to the perturbativeness of the expansion. We then turn to three-nucleon interactions, focusing on operator ambiguities and their interplay with regulator effects. We then discuss the nuclear Green's function Monte Carlo method, going over both wave-function correlations and approximations for the two- and three-body propagators. Finally, following this, we present a range of results on light nuclei: Binding energies and distribution functions are contrasted and compared, starting from several different microscopic interactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Quantum Monte Carlo calculations of weak transitions in A = 6 – 10 nuclei

{\it Ab initio} calculations of the Gamow-Teller (GT) matrix elements in the $\beta$ decays of $^6$He and $^{10}$C and electron captures in $^7$Be are carried out using both variational and Green's function Monte Carlo wave functions obtained from the Argonne $v_{18}$ two-nucleon and Illinois-7 three-nucleon interactions, and axial many-body currents derived from either meson-exchange phenomenology or chiral effective field theory. The agreement with experimental data is excellent for the electron captures in $^7$Be, while theory overestimates the $^6$He and $^{10}$C data by $\sim 2\%$ and $\sim 10\%$, respectively. We show that for these systems correlations in the nuclear wave functions are crucial to explain the data, while many-body currents increase by $\sim 2$--$3\%$ the one-body GT contributions. These findings suggest that the longstanding $g_A$-problem, {\it i.e.}, the systematic overprediction ($\sim 20 \%$ in $A\le 18$ nuclei) of GT matrix elements in shell-model calculations, may be resolved, at least partially, by correlation effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Relativistic effects in ab initio electron-nucleus scattering

The electromagnetic responses obtained from Green's function Monte Carlo (GFMC) calculations are based on realistic treatments of nuclear interactions and currents. The main limitations of this method comes from its nonrelativistic nature and its computational cost, the latter hampering the direct evaluation of the inclusive cross sections as measured by experiments. We extend the applicability of GFMC in the quasielastic region to intermediate momentum transfers by performing the calculations in a reference frame that minimizes nucleon momenta. Additional relativistic effects in the kinematics are accounted for employing the two-fragment model. In addition, we developed a novel algorithm, based on the concept of first-kind scaling, to compute the inclusive electromagnetic cross section of 4 He through an accurate and reliable interpolation of the response functions. A very good agreement is obtained between theoretical and experimental cross sections for a variety of kinematical setups. Finally, this offers a promising prospect for the data analysis of neutrino-oscillation experiments that requires an accurate description of nuclear dynamics in which relativistic effects are fully accounted for.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Computing the dipole polarizability of Ca 48 with increased precision

Here, we compute the electric dipole polarizability of 48 Ca with an increased precision by including more correlations than in previous studies. Employing the coupled-cluster method we go beyond single and double excitations and include leading-order three-particle-three-hole (3p-3h) excitations for the ground state, excited states, and the similarity-transformed operator. We study electromagnetic sum rules, such as the bremsstrahlung sum rule m 0 and the polarizability sum rule αD using interactions from chiral effective field theory. To gauge the quality of our coupled-cluster approximations we perform several benchmarks with the effective interaction hyperspherical harmonics approach in 4 He and with self consistent Green's function in 16 O. We compute the dipole polarizability of 48 Ca employing the chiral interaction N 2 LO sat and the 1.8/2.0 (EM). We find that the effect of 3p-3h excitations in the ground state is small for 1.8/2.0 (EM) but non-negligible for N 2 LO sat . The addition of these new correlations allows us to improve the precision of our 48 Ca calculations and reconcile the recently reported discrepancy between coupled-cluster results based on these interactions and the experimentally determined αD from proton inelastic scattering in 48 Ca. For the computation of electromagnetic and polarizability sum rules, the inclusion of leading-order 3p-3h excitations in the ground state is important, while it is less so for the excited states and the similarity-transformed dipole operator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Perturbative quantum evolution of the gravitational state and dressing in general backgrounds

This paper sets up a perturbative treatment of the evolving quantum state of a gravitational system, in a Schrödinger-like picture, working about a general background. This connects gauge symmetry, the constraints, gravitational dressing, and evolution. Starting with a general time slicing, we give a simple derivation of the relation between the constraints, the Hamiltonian, and its well-known boundary term. Among different approaches to quantization with constraints, we focus on a “gauge-invariant canonical quantization,” which is developed perturbatively in the gravitational coupling. The leading-order solution of the constraints (including the Wheeler-DeWitt equation) for perturbations about the background is given in terms of an explicit construction of gravitational dressings built using certain generalized Green’s functions; different such dressings corresponding to adding propagating gravitational waves to a particular solution of the constraints. Dressed operators commute with the constraints, expressing their gauge invariance, and have an algebraic structure differing significantly from the undressed operators of the underlying field theory. These operators can act on the vacuum to create dressed states, and evolution of general such states is then generated by the boundary Hamiltonian, and alternately may be characterized using other relational observables. This provides a concrete approach to studying perturbative time evolution, including the leading gravitational backreaction, of quantum states of black holes with flat or anti–de-Sitter asymptotics, for example on horizon-crossing slices. This description of evolution in turn provides a starting point for investigating possibly important corrections to quantum evolution, that go beyond quantized general relativity. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Tetraquarks made of sufficiently unequal-mass heavy quarks are bound in QCD

Tetraquarks, bound states composed of two quarks and two antiquarks, have been the subject of intense study but are challenging to understand from first principles. We apply variational and Green’s function Monte Carlo methods to compute tetraquark ground-state energies in potential nonrelativistic QCD using a wide range of color and spatial wave functions. We find no evidence for bound tetraquarks composed of equal-mass quarks and antiquarks. Conversely, we find clear evidence for the existence of bound tetraquarks for sufficiently unequal quark/antiquark mass ratios at all overall mass scales where our effective theory results are applicable. We predict the critical mass ratios for bound state formation and study tetraquark bound states’ spatial and color structure at leading order and next-to-leading order in potential nonrelativistic QCD. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Composite vector particles in external electromagnetic fields

We report that lattice quantum chromodynamics (QCD) studies of electromagnetic properties of hadrons and light nuclei, such as magnetic moments and polarizabilities, have proven successful with the use of background field methods. With an implementation of nonuniform background electromagnetic fields, properties such as charge radii and higher electromagnetic multipole moments (for states of higher spin) can additionally be obtained. This can be achieved by matching lattice QCD calculations to a corresponding low-energy effective theory that describes the static and quasistatic responses of hadrons and nuclei to weak external fields. With particular interest in the case of vector mesons and spin-1 nuclei such as the deuteron, we present an effective field theory of spin-1 particles coupled to external electromagnetic fields. To constrain the charge radius and the electric quadrupole moment of the composite spin-1 field, the single-particle Green’s functions in a linearly varying electric field in space are obtained within the effective theory, providing explicit expressions that can be used to match directly onto lattice QCD correlation functions. Lastly, the viability of an extraction of the charge radius and the electric quadrupole moment of the deuteron from the upcoming lattice QCD calculations of this nucleus is discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Multiple scattering theory for dense plasmas

Dense plasmas occur in stars, giant planets, and in inertial fusion experiments. Accurate modeling of the electronic structure of these plasmas allows for prediction of material properties that can in turn be used to simulate these astrophysical objects and terrestrial experiments. But modeling them remains a challenge. Here we explore the Korringa-Kohn-Rostoker Green's function (KKR-GF) method for this purpose. We find that it is able to predict equation of state in good agreement with other state-of-the-art methods, where they are accurate and viable. In addition, it is shown that the computational cost does not significantly change with temperature, in contrast with other approaches. Moreover, the method does not use pseudopotentials—core states are calculated self consistently. We conclude that KKR-GF is a very promising method for dense plasma simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Short-Range Correlations and the EMC Effect in Effective Field Theory

We propose that the empirical linear relation between the magnitude of the EMC effect in deep inelastic scattering on nuclei and the short-range correlation scaling factor a 2 extracted from high-energy quasielastic scattering at x ≥ 1 is a natural consequence of scale separation and derive the relationship using effective field theory. While the scaling factor a 2 is a ratio of nuclear matrix elements that individually depend on the calculational scheme, we show that the ratio is independent of this choice. We perform Green’s function Monte Carlo calculations with both chiral and Argonne-Urbana potentials to verify this and determine the scaling factors for light nuclei. The resulting values for 3 He and 4 He are in good agreement with experimental values. We additionally present results for 9 Be and 12 C extracted from variational Monte Carlo calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Light-Nuclei Spectra from Chiral Dynamics

In recent years local chiral interactions have been derived and implemented in quantum Monte Carlo methods in order to test to what extent the chiral effective field theory framework impacts our knowledge of few- and many-body systems. Here in this Letter, we present Green’s function Monte Carlo calculations of light nuclei based on the family of local two-body interactions presented by our group in a previous paper in conjunction with chiral three-body interactions fitted to bound- and scattering-state observables in the three-nucleon sector. These interactions include Δ intermediate states in their two-pion-exchange components. We obtain predictions for the energy levels and level ordering of nuclei in the mass range A=4–12, accurate to ≤2% of the binding energy, in very satisfactory agreement with experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Atom Pairing in Optical Superlattices

We study the pairing of fermions in a one-dimensional lattice of tunable double-well potentials using radio-frequency spectroscopy. The spectra reveal the coexistence of two types of atom pairs with different symmetries. Our measurements are in excellent quantitative agreement with a theoretical model, obtained by extending the Green’s function method of Orso et al., [Phys. Rev. Lett. 95, 060402 (2005)], to a bichromatic 1D lattice with non-zero harmonic radial confinement. The predicted spectra comprise hundreds of discrete transitions, with symmetry-dependent initial state populations and transition strengths. Our work provides an understanding of the elementary pairing states in a superlattice, paving the way for new studies of strongly interacting many-body systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Coherent Superconductivity with a Large Gap Ratio from Incoherent Metals

Here, a mysterious incoherent metallic (IM) normal state with T-linear resistivity is ubiquitous among strongly correlated superconductors. Recent progress with microscopic models exhibiting IM transport has presented the opportunity for us to study new models that exhibit direct transitions into a superconducting state out of IM states within the framework of connected Sachdev-Ye-Kitaev “quantum dots.” Here, local Sachdev-Ye-Kitaev interactions within a dot produce IM transport in the normal state, while local attractive interactions drive superconductivity. Through explicit calculations, we find two features of superconductivity arising from an IM normal state. First, despite the absence of quasiparticles in the normal state, the superconducting state still exhibits coherent superfluid transport. Second, the nonquasiparticle nature of the IM Green’s functions produces a large enhancement in the ratio of the zero-temperature superconducting gap Δ and transition temperature T SC , 2Δ/T SC , with respect to its BCS value of 3.53.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Nuclear Charge Radii of B 10 , 11

The first laser spectroscopic determination of the change in the nuclear charge radius for a five-electron system is reported. This is achieved by combining high-accuracy ab initio mass-shift calculations and a high-accuracy measurement of the isotope shift in the 2s 2 2p 2 P 1/2 → 2s 2 3s 2 S 1/2 ground state transition in boron atoms. Accuracy is increased by orders of magnitude for the stable isotopes 10,11 B and the results are used to extract their difference in the mean-square charge radius < r$^{2}_{c}$ > 11– < r$^{2}_{c}$ >10= –0.49(12) fm 2 . The result is qualitatively explained by a possible cluster structure of the boron nuclei and quantitatively used to benchmark new ab initio nuclear structure calculations using the no-core shell model and Green’s function Monte Carlo approaches. These results are the foundation for a laser spectroscopic determination of the charge radius of the proton-halo candidate 8 B.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Fermion Pair Dynamics in Open Quantum Systems

Three-body decay is a rare decay mode observed in a handful of unbound rare isotopes. The angular and energy correlations between emitted nucleons are of particular interest, as they provide invaluable information on the interplay between structure and reaction aspects of the nuclear open quantum system. To study the mechanism of two-nucleon emission, we developed a time-dependent approach that allows us to probe emitted nucleons at long times and large distances. We successfully benchmarked the new method against the Green’s function approach and applied it to low-energy two-proton and two-neutron decays. In particular, we studied the interplay between initial-state nucleon-nucleon correlations and final-state interaction. We demonstrated that the time evolution of the two-nucleon wave function is strongly impacted by the diproton/dineutron dynamics and that the correlations between emitted nucleons provide invaluable information on the dinucleon structure in the initial state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Criticality and Phase Classification for Quadratic Open Quantum Many-Body Systems

Here we study the steady states of translation-invariant open quantum many-body systems governed by Lindblad master equations, where the Hamiltonian is quadratic in the ladder operators, and the Lindblad operators are either linear or quadratic and Hermitian. These systems are called quasifree and quadratic, respectively. We find that steady states of one-dimensional systems with finite-range interactions necessarily have exponentially decaying Green’s functions. For the quasifree case without quadratic Lindblad operators, we show that fermionic systems with finite-range interactions are noncritical for any number of spatial dimensions and provide bounds on the correlation lengths. Quasifree bosonic systems can be critical in D > 1 dimensions. Last, we address the question of phase transitions in quadratic systems and find that, without symmetry constraints beyond invariance under single-particle basis and particle-hole transformations, all gapped Liouvillians belong to the same phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Quantum Monte Carlo Method in the Steady State

We present a numerically exact steady-state inchworm Monte Carlo method for nonequilibrium quantum impurity models. Rather than propagating an initial state to long times, the method is directly formulated in the steady state. This eliminates any need to traverse the transient dynamics and grants access to a much larger range of parameter regimes at vastly reduced computational costs. We benchmark the method on equilibrium Green’s functions of quantum dots in the noninteracting limit and in the unitary limit of the Kondo regime. We then consider correlated materials described with dynamical mean field theory and driven away from equilibrium by a bias voltage. We show that the response of a correlated material to a bias voltage differs qualitatively from the splitting of the Kondo resonance observed in bias-driven quantum dots.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗