Search NASA⌕ Search

SEARCH · Search NASA

Results for “THREE-BODY PROBLEM”

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
DOE OSTI2021

A Canonical Transformation to Eliminate Resonant Perturbations. I

We study dynamical systems that admit action-angle variables at leading order, which are subject to nearly resonant perturbations. If the frequencies characterizing the unperturbed system are not in resonance, the long-term dynamical evolution may be integrated by orbit-averaging over the high-frequency angles, thereby evolving the orbit-averaged effect of the perturbations. It is well known that such integrators may be constructed via a canonical transformation, which eliminates the high-frequency variables from the orbit-averaged quantities. An example of this algorithm in celestial mechanics is the von Zeipel transformation. However, if the perturbations are inside or close to a resonance, i.e., the frequencies of the unperturbed system are commensurate; these canonical transformations are subject to divergences. We introduce a canonical transformation that eliminates the high-frequency phase variables in the Hamiltonian without encountering divergences. This leads to a well-behaved symplectic integrator. We demonstrate the algorithm through two examples: a resonantly perturbed harmonic oscillator and the gravitational three-body problem in mean motion resonance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Computational approaches for three-nucleon systems

Highlights: • Pedagogical review on the computational aspects of the three-nucleon system. • Detailed numerical steps for solving the three-body problem. • Three-nucleon binding and scattering for two different separable potentials. • Triton binding with partial-wave projected chiral potentials. We revisit the three-nucleon system with two kinds of nucleon–nucleon interactions: separable potentials and chiral forces, showing the computational aspects in detail. We start with S-wave separable potentials for which there are simplified forms for the Faddeev equations describing the scattering and binding of three nucleons. We then discuss the partial-wave projected case with a chiral potential considering only one channel so that the computational details can be clearly shown.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Evolution of Efimov states

The Efimov phenomenon manifests itself as an emergent discrete scaling symmetry in the quantum three-body problem. In the unitarity limit, it leads to an infinite tower of three-body bound states with energies forming a geometric sequence. Here in this work, we study the evolution of these so-called Efimov states using relativistic scattering theory. We identify them as poles of the three-particle S matrix and trace their trajectories in the complex energy plane as they evolve from virtual states through bound states to resonances. We dial the scattering parameters toward the unitarity limit and observe the emergence of the universal scaling of energies and couplings - a behavior known from the non-relativistic case. Interestingly, we find that Efimov resonances follow unusual, cyclic trajectories accumulating at the three-body threshold and then disappear at some values of the two-body scattering length. We propose a partial resolution to this ?missing states? problem.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

On the stability of the open-string QED neutron and dark matter

We study the stability of a hypothetical QED neutron, which consists of a color-singlet system of two d quarks and a u quark interacting with the QED interaction. As a quark cannot be isolated, the intrinsic motion of the three quarks in the lowest-energy state may lie predominantly in 1 + 1 dimensions, as in a d-u-d open string. The attractive d-u and u-d QED interactions may overcome the weaker repulsive d-d QED interaction to bind the three quarks together. We examine the QED neutron in a phenomenological three-body problem in 1 + 1 dimensions with an effective interaction extracted from Schwinger’s exact QED solution in 1 + 1 dimensions. The phenomenological model in a variational calculation yields a stable QED neutron at 44.5 MeV. The analogous QED proton with two u quarks and a d quark has been found to be too repulsive to be stable and does not have a bound or continuum state, onto which the QED neutron can decay via the weak interaction. Consequently, the QED neutron is stable against the weak decay, has a long lifetime, and is in fact a QED dark neutron. It may be produced following the deconfinement-to-confinement phase transition of the quark gluon plasma in high-energy heavy-ion collisions. Because of the long lifetime of the QED dark neutron, self-gravitating assemblies of QED dark neutrons or dark antineutrons may be good candidates for a part of the primordial dark matter produced during the phase transition of the quark gluon plasma in the evolution of the early Universe.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Fast emulation of quantum three-body scattering

Here, we develop a class of emulators for solving quantum three-body scattering problems. They are based on combining the variational method for scattering observables and the recently proposed eigenvector continuation concept. The emulators are first trained by the exact scattering solutions of the governing Hamiltonian at a small number of points in its parameter space, and then employed to make interpolations and extrapolations in that space. Through a schematic nuclear-physics model with finite-range two and three-body interactions, we demonstrate the emulators to be extremely accurate and efficient. The computing time for emulation is on the scale of milliseconds (on a laptop), with relative errors ranging from 10 –13 to 10 –4 depending on the case. The emulators also require little memory. We argue that these emulators can be generalized to even more challenging scattering problems. Furthermore, this general strategy may be applicable for building the same type of emulators in other fields, wherever variational methods can be developed for evaluating physical models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Inclusion of electronic polarizability effect in high pressure structural properties of alloy of rare-earth antimonides

In the present paper, we have investigated the high-pressure structural phase transition of rare-earth antimonides (NdSb and DySb). We studied theoretically the structural properties of alloy of these compounds (NdSb and DySb) by using the three-body potential model with the effect of electronic polarizability (TBIPE{sub P}). These compounds exhibit first order crystallographic phase transition from NaCl (B{sub 1}) to CsCl (B{sub 2}) phase at 17.8 GPa and 22.6 GPa respectively. The study has been extended to mixed crystals and the effect of composition on transition pressure and volume change is investigated. The phase transition pressures and associated volume collapse obtained from present potential model (TBIPE{sub P}) show a good agreement with available experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

DYNAMICAL EVOLUTION OF MULTI-RESONANT SYSTEMS: THE CASE OF GJ 876

The GJ 876 system was among the earliest multi-planetary detections outside of the Solar System, and has long been known to harbor a resonant pair of giant planets. Subsequent characterization of the system revealed the presence of an additional Neptune mass object on an external orbit, locked in a three body Laplace mean motion resonance with the previously known planets. While this system is currently the only known extrasolar example of a Laplace resonance, it differs from the Galilean satellites in that the orbital motion of the planets is known to be chaotic. In this work, we present a simple perturbative model that illuminates the origins of stochasticity inherent to this system and derive analytic estimates of the Lyapunov time as well as the chaotic diffusion coefficient. We then address the formation of the multi-resonant structure within a protoplanetary disk and show that modest turbulent forcing in addition to dissipative effects is required to reproduce the observed chaotic configuration. Accordingly, this work places important constraints on the typical formation environments of planetary systems and informs the attributes of representative orbital architectures that arise from extended disk-driven evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Exomoons in Systems with a Strong Perturber: Applications to α Cen AB

The presence of a stellar companion can place constraints on occurrence and orbital evolution of satellites orbiting exoplanets, i.e., exomoons. In this work we revise earlier orbital stability limits for retrograde orbits in the case of a three-body system consisting of a star, planet, and satellite. The revised limit reads a{sub sat}{sup crit}≈0.668(1-1.236e{sub p}) for e {sub p} ≤ 0.8 in units of the Hill Radius and represents the lower critical orbit as a function of the planetary eccentricity e {sub p}. A similar formula is determined for exomoons hosted by planets in binary star systems, where e {sub p} is replaced with the components of free and forced eccentricity from secular orbit evolution theory. By exploring the dynamics of putative exomoons in α Centauri AB we find that the outer stability limit can be much less than half the Hill Radius due to oscillations in the planetary orbital eccentricity caused by the gravitational interaction with the binary star. We show, furthermore, how the resulting truncation of the outer stability limit can affect the outward tidal migration and potential observability of exomoons through transit-timing variations (TTVs). Typical TTV (rms) amplitudes induced by exomoons in binary systems are ≲10 minutes and appear more likely for planets orbiting the less massive stellar component.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Mitigating Green's function Monte Carlo signal-to-noise problems using contour deformations

The Green's function Monte Carlo (GFMC) method provides accurate solutions to the nuclear many-body problem and predicts properties of light nuclei starting from realistic two- and three-body interactions. Controlling the GFMC fermion sign problem is crucial, as the signal-to-noise ratio decreases exponentially with imaginary time, requiring significant computing resources. Inspired by similar scenarios in lattice quantum field theory and spin systems, in this work, we employ integration contour deformations to improve the GFMC signal-to-noise ratio. Machine learning techniques are used to select optimal contours with minimal variance from parametrized families of deformations. As a proof of principle, we consider the deuteron binding energies and Euclidean density response functions. Here, we only observe mild signal-to-noise improvement for the binding energy case. On the other hand, we achieve an order of magnitude reduction of the variance for Euclidean density response functions, paving the way for computing electron- and neutrino-nucleus cross sections of larger nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Mitigating Green's function Monte Carlo signal-to-noise problems using contour deformations

The Green's function Monte Carlo (GFMC) method provides accurate solutions to the nuclear many-body problem and predicts properties of light nuclei starting from realistic two- and three-body interactions. Controlling the GFMC fermion-sign problem is crucial, as the signal-to-noise ratio decreases exponentially with Euclidean time, requiring significant computing resources. Inspired by similar scenarios in lattice quantum field theory and spin systems, in this work, we employ integration contour deformations to improve the GFMC signal-to-noise ratio. Machine learning techniques are used to select optimal contours with minimal variance from parameterized families of deformations. As a proof of principle, we consider the deuteron binding energies and Euclidean density response functions. We only observe mild signal-to-noise improvement for the binding energy case. On the other hand, we achieve an order of magnitude reduction of the variance for Euclidean density response functions, paving the way for computing electron- and neutrino-nucleus cross-sections of larger nuclei.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Coupled-channel approach to isotensor π π π scattering from lattice QCD

The quest to understand three-body dynamics from first-principle QCD includes the study of nonresonant and resonant systems. The isospin I = 2 system is of particular interest having no three-body resonance but featuring a resonance in a subchannel, while also being a coupled-channel problem. In this study, we calculate the finite-volume spectrum from lattice QCD at two different pion masses, map the amplitude to the infinite volume through a generalized Finite-Volume Unitarity three-body quantization condition, investigate the limit of a narrow ρ , and compare with an effective Lagrangian prediction at leading order. Chiral extrapolations between different pion masses are performed.

Feng, Yuchuan [The George Washington University] (↗
DOE OSTI2025

Toward scalable bound-to-resonance extrapolations for few- and many-body systems

In open quantum many-body systems, the theoretical description of resonant states of many particles strongly coupled to the continuum can be challenging. Such states are commonplace in, for example, exotic nuclei and hadrons, and can reveal important information about the underlying forces at play in these systems. In this work, we demonstrate that the complex-augmented eigenvector continuation (CA-EC) method, originally formulated for the two-body problem with uniform complex scaling, can reliably perform bound-to-resonance extrapolations for genuine three-body resonances having no bound subsystems. Here, we first establish that three-body bound-to-resonance extrapolations are possible by benchmarking different few-body approaches, and we provide arguments to explain how the extrapolation works in the many-body case. We furthermore pave the way towards scalable resonance extrapolations in many-body systems by showing that the CA-EC method also works in the Berggren basis, studying a realistic application using the Gamow shell model.

Ab initio calculations↗
DOE OSTI2024

In-medium similarity renormalization group with flowing 3-body operators, and approximations thereof

Here, we explore the impact of retaining three-body operators within the in-medium similarity renormalization group (IMSRG), as well as various approximations schemes. After studying two toy problems, identical fermions with a contact interaction and the Lipkin-Meshkov-Glick model, we employ the valence-space formulation of the IMSRG to investigate the even- A carbon isotopes with a chiral two-body potential. We find that retaining only those commutators expressions that scale as N 7 provides an excellent approximation of the full three-body treatment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

How to renormalize coupled cluster theory

Coupled cluster theory is an attractive tool to solve the quantum many-body problem because its singles and doubles (CCSD) approximation is computationally affordable and yields about 90% of the correlation energy. Capturing the remaining 10%, e.g., via including triples, is numerically expensive. In this work, we assume that short-range three-body correlations dominate and—following Lepage (arXiv:nucl-th/9706029)—that their effects can be included within CCSD by renormalizing the three-body contact interaction. We renormalize this contact in 16 O and obtain systematically improved CCSD results for 24 O, 20–34 Ne, 40,48 Ca, 78 Ni, 90 Zr, and 100 Sn.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Benchmark calculations of infinite neutron matter with realistic two- and three-nucleon potentials

In this study, we present the equation of state of infinite neutron matter as obtained from highly realistic Hamiltonians that include nucleon-nucleon and three-nucleon coordinate-space potentials. We benchmark three independent many-body methods: Brueckner-Bethe-Goldstone (BBG), Fermi hypernetted chain/single-operator chain (FHNC/SOC), and auxiliary-field diffusion Monte Carlo (AFDMC). We find them to provide similar equations of state when the Argonne v 18 and the Argonne v' 6 nucleon-nucleon potentials are used in combination with the Urbana IX three-body force. Only at densities larger than about 1.5 the nuclear saturation density ( ρ 0 =0.16 fm ⁻3 ) the FHNC/SOC energies are appreciably lower than the other two approaches. The AFDMC calculations carried out with all of the Norfolk potentials fitted to reproduce the experimental trinucleon ground-state energies and nd doublet scattering length yield unphysically bound neutron matter, associated with the formation of neutron droplets. Including tritium β decay in the fitting procedure, as in the second family of Norfolk potentials, mitigates but does not completely resolve this problem. An excellent agreement between the BBG and AFDMC results is found for the subset of Norfolk interactions that do not make neutron-matter collapse, while the FHNC/SOC equations of state are moderately softer.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Variational Monte Carlo Calculations of A ≤ 4 Nuclei with an Artificial Neural-Network Correlator Ansatz

Here, the complexity of many-body quantum wave functions is a central aspect of several fields of physics and chemistry where nonperturbative interactions are prominent. Artificial neural networks (ANNs) have proven to be a flexible tool to approximate quantum many-body states in condensed matter and chemistry problems. In this work we introduce a neural-network quantum state ansatz to model the ground-state wave function of light nuclei, and approximately solve the nuclear many-body Schrodinger equation. Using efficient stochastic sampling and optimization schemes, our approach extends pioneering applications of ANNs in the field, which present exponentially scaling algorithmic complexity. We compute the binding energies and point-nucleon densities of A ≤ 4 nuclei as emerging from a leading-order pionless effective field theory Hamiltonian. We successfully benchmark the ANN wave function against more conventional parametrizations based on two- and three-body Jastrow functions, and virtually exact Green's function Monte Carlo results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Variational Monte Carlo calculations of n + H 3 scattering

A paramount goal in nuclear physics is to unify ab initio treatments of bound and unbound states. The position-space quantum Monte Carlo (QMC) methods have a long history of successful bound-state calculations in light systems but have seen minimal implementation in unbound systems. Here, we introduce a numerical method to improve the efficiency and accuracy of unbound-state calculations in QMC. As an initial application, we compute scattering observables for the smallest system available to probe three-body forces, the neutron-triton system, using variational Monte Carlo (VMC) wave functions. The method involves inferring long-range amplitudes in the wave function from integrals over the short-range region where all the particles interact. This approach using integral relations is well established in the literature; here, we develop it for the QMC framework. We validate our approach with a consistency check between short-range spectroscopic overlap functions computed from direct evaluation and from the integral relations; scattering amplitudes are long-range asymptotics of those overlaps. Comparison against published benchmark calculations using the same potential demonstrates that when applied to the current VMC wave functions, the integral method produces more accurate scattering observables than direct evaluation from the same variational wave function. However, it still differs noticeably from the exact results. Using additional interactions, we then present phase shifts and mixing parameters for the n + 3 H system. In particular, we present one of the first applications of the Norfolk family of local coordinate-space chiral potentials in unbound systems of A > 2. The Norfolk results accurately describe s-wave scattering but predict p-wave cross sections too large. Compared with previous QMC scattering calculations, the integral method avoids difficulties associated with the precise computation of energy differences and with convergence outside the interaction region, which is particularly severe in the variational calculation. Application of the integral method here paves the way for its use in Green's function Monte Carlo (GFMC) calculations. In GFMC, the wave functions are more accurate, but the high-precision convergence of their tails is slow, and there are additional difficulties in reading out amplitudes. Here, the integral methods will address both of those remaining problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Improved structure of calcium isotopes from ab initio calculations

The in-medium similarity renormalization group (IMSRG) is a powerful and flexible many-body method to compute the structure of nuclei starting from nuclear forces. Recent developments have extended the IMSRG from its standard truncation at the normal-ordered two-body level, the IMSRG(2), to a precision approximation including normal-ordered three-body operators, the IMSRG(3)-N 7 . This improvement provides a more precise solution to the many-body problem and makes it possible to quantify many-body uncertainties in IMSRG calculations. We explore the structure of 44,48,52 Ca using the IMSRG(3)-N 7 , focusing on understanding existing discrepancies of the IMSRG(2) to experimental results. We find a significantly better description of the first 2 + excitation energy of 48 Ca, improving the description of the shell closure at N=28. At the same time, we find that the IMSRG(3)-N 7 corrections to charge radii do not resolve the systematic underprediction of the puzzling large charge radius difference between 52 Ca and 48 Ca. We present estimates of many-body uncertainties of IMSRG(2) calculations applicable also to other systems based on the size extensivity of the method.

39 ≤ A ≤ 58↗