Search NASA⌕ Search

SEARCH · Search NASA

Results for “three-body system”

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 OSTI2023

A Formulation of the Many-Body Expansion (MBE) for Periodic Systems: Application to Several Ice Phases

A new formulation of the many-body expansion (MBE) for periodic systems, based on their inherent translational symmetry, is introduced and applied to 7 ice polymorphs. This new formulation is built via a hierarchical procedure connecting gas phase unit cells over finite supercells to infinite solids. For finite supercells, it is demonstrated that this method successfully recovers the energetics while reducing the scaling of the calculation of the many-body terms by a factor of ~N, where N is the system size. Furthermore, the proposed framework delineates a straightforward way to obtain properties of macroscopic systems in the limit of an infinite cell. For periodic systems, the success of this approach is demonstrated by showing that the lattice energies computed (up to the 4-th order in the MBE) reproduce the lattice energies obtained using periodic boundary conditions with an Ewald summation for 7 polymorphs of ice (Ih, II, VIII, IX, XIII, XIV, XV). This development makes it possible to quantify, for the first time, the many-body contributions to the lattice energy of various ice polymorphs. The many-body (three-body and higher) interactions were found to vary significantly among the 7 ice phases, amounting between 7-24% of the total lattice energies. So, this development opens the door for obtaining insights into solid-state properties, while leveraging the computational benefits of the MBE.

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

Many-Body Potential for Simulating Self-Assembly of Polymer-Grafted Nanoparticles in a Polymer Matrix

Many-body interactions between polymer-grafted nanoparticles (NPs) play a key role in promoting their assembly into low-dimensional structures within polymer melts, even when the particles are spherical and isotropically grafted. However, capturing such interactions in simulations of NP assembly is very challenging because explicit modeling of the polymer grafts and melt chains is highly computationally expensive, even using coarse-grained models. Here, we develop a many-body potential for describing the effective interactions between spherical polymer-grafted NPs in a polymer matrix through a machine-learning approach. The approach involves using permutationally invariant polynomials to fit two- and three-body interactions derived from the potential of mean force calculations. The potential developed here reduces the computational cost by several orders of magnitude, thereby, allowing us to explore assembly behavior over large length and time scales. We show that the potential not only reproduces previously known assembled phases such as 1D strings and 2D hexagonal sheets, which generally cannot be achieved using isotropic two-body potentials, but can also help discover interesting phases such as networks, clusters, and gels. We demonstrate how each of these assembly morphologies intrinsically arises from a competition between two- and three-body interactions. Our approach for deriving many-body effective potentials can be readily extended to other colloidal systems, enabling researchers to make accurate predictions of their behavior and dissect the role of individual interaction energy terms of the overall potential in the observed behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Descriptors of water aggregation

For this work, we rely on a total of 23 (cluster size, 8 structural, and 14 connectivity) descriptors to investigate structural patterns and connectivity motifs associated with water cluster aggregation. In addition to the cluster size n (number of molecules), the 8 structural descriptors can be further categorized into (i) one-body (intramolecular): covalent OH bond length (r OH ) and HOH bond angle (θ HOH ), (ii) two-body: OO distance (r OO ), OHO angle (θ OHO ), and HOOX dihedral angle ($\phi$ HOOX ), where X lies on the bisector of the HOH angle, (iii) three-body: OOO angle (θ OOO ), and (iv) many-body: modified tetrahedral order parameter (q) to account for two-, three-, four-, five-coordinated molecules (q m , m = 2, 3, 4, 5) and radius of gyration (R g ). The 14 connectivity descriptors are all many-body in nature and consist of the AD, AAD, ADD, AADD, AAAD, AAADD adjacencies [number of hydrogen bonds accepted (A) and donated (D) by each water molecule], Wiener index, Average Shortest Path Length, hydrogen bond saturation (% HB), and number of non-short-circuited three-membered cycles, four-membered cycles, five-membered cycles, six-membered cycles, and seven-membered cycles. We mined a previously reported database of 4 948 959 water cluster minima for (H 2 O) n , n = 3–25 to analyze the evolution and correlation of these descriptors for the clusters within 5 kcal/mol of the putative minima. It was found that r OH and % HB correlated strongly with cluster size n, which was identified as the strongest predictor of energetic stability. Marked changes in the adjacencies and cycle count were observed, lending insight into changes in the hydrogen bond network upon aggregation. A Principal Component Analysis (PCA) was employed to identify descriptor dependencies and group clusters into specific structural patterns across different cluster sizes. The results of this study inform our understanding of how water clusters evolve in size and what appropriate descriptors of their structural and connectivity patterns are with respect to system size, stability, and similarity. The approach described in this study is general and can be easily extended to other hydrogen-bonded systems.

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

Light nuclei with semilocal momentum-space regularized chiral interactions up to third order

In this work, we present a systematic investigation of few-nucleon systems and light nuclei using the current Low Energy Nuclear Physics International Collaboration interactions comprising semilocal momentum-space regularized two- and three-nucleon forces up to third chiral order (N 2 LO). Following our earlier study utilizing the coordinate-space regularized interactions, the two low-energy constants entering the three-body force are determined from the triton binding energy and the differential cross-section minimum in elastic nucleon-deuteron scattering. Predictions are made for selected observables in elastic nucleon-deuteron scattering and in the deuteron breakup reactions, for properties of the A = 3 and A = 4 nuclei, and for spectra of p-shell nuclei up to A = 16. A comprehensive error analysis is performed including an estimation of correlated truncation uncertainties for nuclear spectra. The obtained predictions are generally found to agree with experimental data within errors. Similarly to the coordinate-space regularized chiral interactions at the same order, a systematic overbinding of heavier nuclei is observed, which sets in for A~10 and increases with A.

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 OSTI2023

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↗
DOE OSTI2024

Magnetic structure of A ≤ 10 nuclei using the Norfolk nuclear models with quantum Monte Carlo methods

Here we present quantum Monte Carlo calculations of magnetic moments, form factors, and densities of A ≤ 10 nuclei within a chiral effective field theory approach. We use the Norfolk two- and three-body chiral potentials and their consistent electromagnetic one- and two-nucleon current operators. We find that two-body contributions to the magnetic moment can be large (up to ≈ 33% in A = 9 systems). We study the model dependence of these observables and place particular emphasis on investigating their sensitivity to using different cutoffs to regulate the many-nucleon operators. Calculations of elastic magnetic form factors for A ≤ 10 nuclei show excellent agreement with the data out to momentum transfers q ≈ 3 fm -1 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Quantum Monte Carlo calculations of electron scattering from 12 C in the short-time approximation

The short-time approximation is a method introduced to evaluate electroweak nuclear response for systems with A ≥ 12, extending the reach of first-principle many-body quantum Monte Carlo calculations. Using realistic two- and three-body nuclear interactions and consistent one- and two-body electromagnetic currents, we calculate longitudinal and transverse response densities and response functions of 12 C. Here, we compare the resulting cross sections with experimental data for electron-nucleus scattering, finding good agreement.

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

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↗
DOE OSTI2024

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↗
DOE OSTI2024

Unitary coupled-channel three-body amplitude with pions and kaons

Three-body dynamics above threshold is required for the reliable extraction of many amplitudes and resonances from experiment and lattice QCD. The S-matrix principle of unitarity can be used to construct dynamical coupled-channel approaches in which three particles scatter off each other, rearranging two-body subsystems by particle exchange. This paper reports the development of a three-body coupled-channel, amplitude including pions and kaons. The unequal-mass amplitude contains two-body S- and P-wave subsystems (“isobars”) of all isospins, 𝐼 = 0,1/2,1,3/2,2 , and it also allows for transitions within a given isobar. The 𝑓 0 ⁡(500)⁢(𝜎),𝑓 0 ⁡(980),𝜌⁡(700),𝐾$^{*}_{0}$⁡(700)⁢(𝜅), and 𝐾*⁡(892) resonances are included, apart from repulsive isobars. Different methods to evaluate the amplitude for physical momenta are discussed. Production amplitudes for 𝑎 1 quantum numbers are shown as a proof of principle for the numerical implementation.

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 OSTI2025

Investigating the p--π ± and p--p--π ± dynamics with femtoscopy in pp collisions at $\sqrt{s}$ = 13 TeV

The interaction between pions and nucleons plays a crucial role in hadron physics. It represents a fundamental building block of the low-energy QCD dynamics and is subject to several resonance excitations. This work studies the p--π ± dynamics using femtoscopic correlations in high-multiplicity pp collisions at $\sqrt{s}$ = 13 TeV measured by ALICE at the LHC. As the final-state interaction between protons and pions is well constrained by scattering experiments and the study of pionic hydrogen, the results give access to information on the particle-emitting source in pp collisions using the femtoscopy methods. The scaling of the source size of primordial protons and pions against their pair transverse mass is extracted. The results are compared with the source sizes studied with p–p, p--K + , and π ± –π ± pairs by ALICE in the same collision system and are found to be in agreement for the different particle pairs. This reinforces recent findings by ALICE of a common emission source for all hadron-pairs in pp collisions at LHC energies. Furthermore, the p--p--π ± systems are studied using three-particle femtoscopy in pp collisions at TeV. The presence of three-body effects is analyzed utilizing the cumulant expansion method. In this formalism, the known two-body interactions are subtracted in order to isolate the three-body effects. For both, p--p--π + and p--p--π – , a non-zero cumulant is found, indicating effects beyond pairwise interactions. These results give information on the coupling of the pion to multiple nucleons.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗