Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three body”

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 37 records · Page 2

Newton versus the machine: solving the chaotic three-body problem using deep neural networks

ABSTRACT Since its formulation by Sir Isaac Newton, the problem of solving the equations of motion for three bodies under their own gravitational force has remained practically unsolved. Currently, the solution for a given initialization can only be found by performing laborious iterative calculations that have unpredictable and potentially infinite computational cost, due to the system’s chaotic nature. We show that an ensemble of converged solutions for the planar chaotic three-body problem obtained using an arbitrarily precise numerical integrator can be used to train a deep artificial neural network (ANN) that, over a bounded time interval, provides accurate solutions at a fixed computational cost and up to 100 million times faster than the numerical integrator. In addition, we demonstrate the importance of training an ANN using converged solutions from an arbitrary precise integrator, relative to solutions computed by a conventional fixed precision integrator, which can introduce errors in the training data, due to numerical round-off and time discretization, that are learned by the ANN. Our results provide evidence that, for computationally challenging regions of phase space, a trained ANN can replace existing numerical solvers, enabling fast and scalable simulations of many-body systems to shed light on outstanding phenomena such as the formation of black hole binary systems or the origin of the core collapse in dense star clusters.

Breen, Philip G.↗

Unitary Bethe-Salpeter Methods in Two- and Three-Body Systems

The study of resonances, unstable particles formed in particle scattering, is motivated by the questions about the strong interaction and the composition of hadrons. In order to advance this study in the intermediate energy region, a variety of techniques are used. In this thesis, we present work done on two different analyses that use the principle of unitarity and the Bethe-Salpeter equation to study the resonances Lambda(1405), Sigma(1385), and a1(1260). Firstly, we perform a simultaneous analysis of Sand P-waves of the Strangeness = -1 two-body meson-baryon scattering amplitude using all low-energy data. For the first time, differential cross section data are included for chiral unitary coupled-channel models. From this model, S- and P-wave amplitudes are extracted and we observe both well-known I(JP)=0(1/2^-) S-wave states as well as a new I(JP)=1(1/2^+) state absent in quark models and lattice QCD results. Multiple statistical and phenomenological tests suggest that, while the data clearly require an I =1 P-wave resonance, the new state just accounts for the absence of the decuplet Sigma(1385)3/2+ in the model which is subsequently included in the parameterization. Secondly, we formulate the final state interaction of the a1(1260) resonance decay in a manifestly three-body-unitary parameterization and fit it to the a1(1260) lineshape measured by the ALEPH experiment. Dalitz plots calculated from this fit are presented. The work demonstrates the feasibility to numerically solve a previously derived amplitude and its generalization to isobars with spin and coupled channels. The model is a good test case because it can also be applied to other meson decays including exotic states and modified for the finite-volume problem as it arises in lattice QCD due to its manifest unitarity. The understanding of the resonances being analyzed in both the two-body work and the three-body work gives a better understanding of intermediate energy QCD.

Sadasivan, Daniel↗

6Li as a three-body system in the (p, 3He) reaction at astrophysical energies

Several astrophysical processes are governed by the occurrence of nuclear reactions involving light nuclei at energies below the Coulomb barrier. Among other effects, their understanding is challenged by the appearance of clustered structures in the ground-state configuration of some of these nuclei, which may have a significant impact on the reaction cross section. In this contribution, we focus on the Li 6 (p, He 3 ) He 4 reaction, to probe the role of clustered configurations of Li 6 . In particular, we consider a three-body ab-initio calculation, based on the hyperspherical harmonics (HH) method, of the Li 6 wave function (WF), together with a more phenomenological three-body model. We observe that the HH WF entails a degree of clustering much larger than obtained from the phenomenological WFs. However, the corresponding reaction cross section, evaluated as a direct two-nucleon transfer in distorted-wave Born approximation, still follows the scaling with the clustering strength already pointed out in a previous work [1] and exhibits an energy trend very similar to that obtained with realistic phenomenological WFs. This opens up interesting perspectives towards constraining

Perrotta, S.↗

Exploring the Strong Interaction of Three-Body Systems at the LHC

Deuterons are atomic nuclei composed of a neutron and a proton held together by the strong interaction. Unbound ensembles composed of a deuteron and a third nucleon have been investigated in the past using scattering experiments, and they constitute a fundamental reference in nuclear physics to constrain nuclear interactions and the properties of nuclei. In this work, K + -d and p-d femtoscopic correlations measured by the ALICE Collaboration in proton-proton (pp) collisions at $\sqrt{s}$ = 13 TeV at the Large Hadron Collider (LHC) are presented. It is demonstrated that correlations in momentum space between deuterons and kaons or protons allow us to study three-hadron systems at distances comparable with the proton radius. The analysis of the K + -d correlation shows that the relative distances at which deuterons and protons or kaons are produced are around 2 fm. The analysis of the p-d correlation shows that only a full three-body calculation that accounts for the internal structure of the deuteron can explain the data. In particular, the sensitivity of the observable to the short-range part of the interaction is demonstrated. These results indicate that correlations involving light nuclei in pp collisions at the LHC will also provide access to any three-body system in the strange and charm sectors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Bound states in the B -matrix formalism for the three-body scattering

In this work, we consider a model of relativistic three-body scattering with a bound state in the two-body subchannel. We show that the naiïve K-matrix type parametrization, here referred to as the B-matrix, has nonphysical singularities near the physical region. We show how to eliminate such singularities by using dispersion relations and also show how to reproduce unitarity relations by taking into account all relevant open channels.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Symmetrizing relativistic three-body partial wave amplitudes

S matrix principles and symmetries impose constraints on three-particle scattering amplitudes, which can be formulated as a class of integral equations for their partial wave projections. However, these amplitudes are typically expressed in an asymmetric basis, where one of the initial and final state particles is singled out, and all quantum numbers are defined relative to this spectator. In this work, we show how to construct symmetric partial wave amplitudes, which have been symmetrized over all possible spectator combinations, using their asymmetric counterparts and sets of recoupling coefficients. We derive these recoupling coefficients for arbitrary angular momentum and isospin for arbitrary systems of spinless particles with SU(2) flavor symmetry. We propose a simple intensity observable suitable for visualizing the structure of three-body dynamics in Dalitz distributions. Finally, we provide some numerical examples of Dalitz distributions relevant for future lattice QCD calculations of 3⁢𝜋 systems and provide numerical evidence that the symmetrization procedure is consistent with expected symmetries of Dalitz plots.

Lattice QCD↗

Lattice model of heavy-light three-body system

We present a study of a $1+1$ dimensional heavy-light three-body system in finite volume. The heavy-light system is simulated by a coupled-channel $\phi^4$ type lattice model, and both ground state and excited states of multiparticle energy spectra are measured on various lattices. The lattice simulation data analysis is performed based on variational approach.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Tidal Origin for a Three-body Resonance in Kepler-221

Over the course of the last two decades, traditional models of planet formation have been repeatedly challenged by the emerging census of extrasolar planets. Key among them is the orbital architecture problem: while standard models of orbital migration predict resonant orbits for short-period objects, most planets do not appear to lie in orbital resonances. Here, we show that the four-planet system Kepler-221, not previously recognized to have active orbital resonances, has a three-body commensurability relation unique within the Kepler sample. Using a suite of numerical experiments as well as a perturbative analysis, we demonstrate that this system likely began as a resonant chain and proceeded to undergo large-scale divergence away from resonance, under the action of tidal dissipation. Our results further indicate that obliquity tides, driven by a secular spin–orbit resonance and mutual inclination, are an excellent candidate for driving this orbital divergence, and that the high tidal luminosity may also explain the anomalous size of planet b, which lies within the Fulton radius gap.

79 ASTRONOMY AND ASTROPHYSICS↗

Dalitz plots and lineshape of a 1 ( 1260 ) from a relativistic three-body unitary approach

We formulate the final-state interaction of the a 1 ( 1260 ) resonance decay in a manifestly three-body unitary parametrization and fit it to the a 1 ( 1260 ) line shape measured by the ALEPH experiment. Dalitz plots calculated from this fit are presented. The work demonstrates the feasibility to numerically solve a previously derived amplitude and its generalization to isobars with spin and coupled channels. The model can also be applied to other meson decays and modified for the finite-volume problem as it arises in lattice QCD due to its manifest unitarity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Randomized low-rank decompositions of nuclear three-body interactions

First-principles simulations of many-fermion systems are commonly limited by the computational requirements of processing large data objects. As a remedy, we propose the use of low-rank approximations of three-body interactions, which are the dominant such limitation in nuclear physics. We introduce a randomized decomposition technique to handle the excessively large matrix dimensions and study the sensitivity of low-rank properties to interaction details. The developed low-rank three-nucleon interactions are benchmarked in ab initio simulations of few- and many-body systems. Exploiting low-rank properties provides a promising route to extend the microscopic description of atomic nuclei to large systems where storage requirements exceed the computational capacities of the most advanced high-performance computing facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantum Monte Carlo studies of a trimer scaling function with microscopic two- and three-body interactions

Here, we present an energy scaling function to predict, in a specific range, the energy of bosonic trimers with large scattering lengths and finite range interactions, which is validated by quantum Monte Carlo calculations using microscopic Hamiltonians with two- and three-body potentials. The proposed scaling function depends on the scattering length, effective range, and a reference energy, which we chose as the trimer energy at unitarity. We obtained the scaling function as a limit cycle from the solution of the renormalized zero-range model with effective range corrections. We proposed a simple parametrization of the energy scaling function. Besides the intrinsic interest in theoretical and experimental investigations, this scaling function allows one to probe Efimov physics with only the trimer ground states, which may open opportunities to identify Efimov trimers whenever access to excited states is limited.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Simulations of classical three-body thermalization in one dimension

One-dimensional systems, such as nanowires or electrons moving along strong magnetic field lines, have peculiar thermalization physics. The binary collision of pointlike particles, typically the dominant process for reaching thermal equilibrium in higher-dimensional systems, cannot thermalize a 1D system. We study how dilute classical 1D gases thermalize through three-body collisions. We consider a system of identical classical point particles with pairwise repulsive inverse power-law potential V ij ∝ 1/|x i –x j | n or the pairwise Lennard-Jones potential. Using Monte Carlo methods, we compute a collision kernel and use it in the Boltzmann equation to evolve a perturbed thermal state with temperature T toward equilibrium. We explain the shape of the kernel and its dependence on the system parameters. Additionally, we implement molecular dynamics simulations of a many-body gas and show agreement with the Boltzmann evolution in the low-density limit. For the inverse power-law potential, the rate of thermalization is proportional to ρ 2 ⁢T$\frac{1}{2}$ – $\frac{1}{n}$, where ρ is the number density. Furthermore, the corresponding proportionality constant decreases with increasing n.

1-dimensional systems↗

Measurement of the absolute branching fraction of the three-body decay Λ c + → Ξ 0 K + π 0 and search for Λ c + → n K + π 0 , Σ 0 K + π 0 , and Λ K + π 0

The Cabbibo-favored decay Λ c + → Ξ 0 K + π 0 is studied for the first time using 6.1 fb − 1 of e + e − collision data at center-of-mass energies between 4.600 and 4.840 GeV, collected with the BESIII detector at the BEPCII collider. With a double-tag method, the branching fraction of the three-body decay Λ c + → Ξ 0 K + π 0 is measured to be ( 7.79 ± 1.46 ± 0.95 ) × 10 − 3 , where the first and second uncertainties are statistical and systematic, respectively. The branching fraction of the two-body decay Λ c + → Ξ ( 1530 ) 0 K + is ( 5.99 ± 1.04 ± 0.32 ) × 10 − 3 , which is consistent with the previous result of ( 5.02 ± 0.99 ± 0.31 ) × 10 − 3 . In addition, the upper limit on the branching fraction of the doubly Cabbibo-suppressed decay Λ c + → n K + π 0 is 7.1 × 10 − 4 at the 90% confidence level. The upper limits on the branching fractions of Λ c + → Σ 0 K + π 0 and Λ K + π 0 are also determined to be 1.8 × 10 − 3 and 2.0 × 10 − 3 , respectively. Published by the American Physical Society 2024

Ablikim, M.↗

Partial-wave projection of the one-particle exchange in three-body scattering amplitudes

As the study of three-hadron physics from lattice QCD matures, it is necessary to develop proper analysis tools in order to reliably study a variety of phenomena, including resonance spectroscopy and nuclear structure. Reconstructing the three-particle scattering amplitude requires solving integral equations, which can be written in terms of data-constrained dynamical functions and physical on shell quantities. The driving term in these equations is the so-called one-particle exchange, which leads to a kinematic divergence for particles on mass shell. A vital component in defining three-particle amplitudes with definite parity and total angular momentum, which are used in spectroscopic studies, is to project the one-particle exchange into definite partial waves. We present a general procedure to construct exact analytic partial-wave projections of the one-particle exchange contribution for any system composed of three spinless hadrons. Our result allows one full control over the analytic structure of the projection, which we explore for some low-lying partial waves with applications to three pions. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Impact of three-body forces on elastic nucleon-nucleus scattering observables

In a previous series of papers we investigated the domain of applicability of chiral potentials to the construction of a microscopic optical potential (OP) for elastic nucleon-nucleus scattering. The OP was derived at the first order of the spectator expansion of the Watson multiple scattering theory and its final expression was a folding integral between the nucleon-nucleon ( N N ) t matrix and the nuclear density of the target. In the calculations N N and three-nucleon ( 3 N ) chiral interactions were used for the target density and only the N N interaction for the N N t matrix. The purpose of this study is to achieve another step towards the calculation of a more consistent OP introducing the 3 N force also in the dynamic part of the OP. The full treatment of the 3 N interaction is beyond our present capabilities. Thus, here it is approximated with a density dependent N N interaction obtained after the averaging over the Fermi sphere. In practice, in our model the 3 N force acts as a medium correction of the bare N N interaction used to calculate the t matrix. Even if the 3 N force is treated in an approximate way, this method naturally extends our previous model of the OP and allows a direct comparison of our present and previous results. We consider as case studies the elastic scattering of nucleons off C 12 and O 16 . We present results for the differential cross section and the spin observables for different values of the projectile energy. From the comparison with the experimental data and with the results of our previous model we assess the importance of the 3 N interaction in the dynamic part of the OP. Conclusions: Our analysis indicates that the contribution of the 3 N force in the t matrix is small for the differential cross section and it is sizable for the spin observables, in particular, for the analyzing power. We find that the two-pion exchange term is the major contributor to the 3 N force. A chiral expansion order-by-order analysis of the scattering observables confirms the convergence of our results at the next-to-next-to-next-to-leading-order, as already established in our previous work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantum dynamics of positron-hydrogen scattering and three-body bound state formation with an assisting laser field: predictions of a reduced-dimensionality model

Abstract We investigate the quantum dynamics of target excitation and positronium formation in the positron-hydrogen atom scattering without and with an external assisting laser field within a reduced-dimensional quantum model. Strong interference fringes between the incident and reflected positron wave packets are observed in the reaction region. We further investigate the critical behavior of transition probabilities near the channel-opening thresholds for hydrogen excitation and positronium formation and find a strong competition between channels with similar threshold energies, but different parities. The transmission ratios of the incident positron in different reaction channels are calculated, and it is shown that only positronium formation in the ground state prefers forward scattering. Our simulation of the positron-hydrogen scattering with an assisting laser field indicates that the three-particle bound states can be formed during the collisions due to the photon emission induced by the external laser field.

Optics↗

Role of three-body dynamics in nucleon-deuteron correlation functions

Correlation functions of hadrons can be accessed in high-energy collisions of atomic nuclei, revealing information about the underlying interaction. This work complements experimental efforts to study nucleon-deuteron Nd—with N=p (proton) or N=n (neutron)—correlations with theory evaluations using different techniques. The correlation functions C nd and C pd are calculated based on a scattering wave function, extending previous benchmarks for the Nd scattering matrix to this new observable. We use hyperspherical harmonics and Faddeev techniques with one of the widely used nucleon-nucleon (NN) interactions, the Argonne v18 potential. Moreover, in the low-energy region we perform additional calculations in the framework of pionless effective field theory. The pd correlation function is computed in the large-energy region to make contact with a recent measurement by the ALICE Collaboration. We show that the scattering wave function has the proper dynamical input to describe an initial rise and subsequent oscillations of C pd as a function of the energy. Effects on the observables using different NN and three-nucleon potentials are evaluated with the conclusion that variations of around 2% are observed. Although these effects are small, future measurements can go beyond this accuracy, allowing for new detailed studies of strong interaction in light nuclear systems. Furthermore, the present study supports the current efforts devoted to the measurement of correlation functions in systems dominated by the strong interactions, such as pd, ppp, Λd, and ppΛ.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗