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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Merging ab initio theory and few-body approach for (d, p) reactions

A framework for A(d, p)B reactions is introduced by merging the microscopic approach to computing the properties of the nucleon-target systems and the three-body n + p + A reaction formalism, thus providing a consistent link between the reaction cross sections and the underlying microscopic structure. In this first step toward a full microscopic description, we focus on the inclusion of the neutron-target microscopic properties. The properties of the neutron-target subsystem are encapsulated in the Green's function which is computed with the coupled Cluster theory using a chiral nucleon-nucleon and three-nucleon interactions. Subsequently, this many-body information is introduced in the few-body Green's Function Transfer approach to (d, p) reactions. Our benchmarks on stable targets 40,48 Ca show an excellent agreement with the data. We then proceed to make specific predictions for (d, p) on neutron rich 52,54 Ca isotopes. These predictions are directly relevant to testing the new magic numbers N = 32, 34 and are expected to be feasible in the first campaign of the projected FRIB facility.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 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 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 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 OSTI2023

High-pressure CO 2 dissociation with nanosecond pulsed discharges

Herein we investigate the conversion of CO 2 into CO and O 2 with nanosecond repetitively pulsed (NRP) discharges in a high-pressure batch reactor. Stable discharges are obtained at up to 12 bar. By-products are measured with gas chromatography. The energy efficiency is determined for a range of processing times, pulse energy, and fill pressures. It is only weakly sensitive to the plasma operating parameters, i.e the extent of CO 2 conversion is almost linearly-dependent on the specific energy invested. A conversion rate as high as 14% is achieved with an energy efficiency of 23%. For long processing times, saturation in the yield and a drop in efficiency are observed, due to the increasing role of three-body recombination reactions, as described by zero-dimensional detailed kinetic modeling. The modeling reveals the presence of three-stage kinetics between NRP pulses, controlled by electron-impact CO 2 dissociation, vibrational relaxation, and neutral elementary kinetics. Transport effects are shown to be important for CO 2 conversion at high pressures. For fill pressures beyond 10 bar, CO 2 may locally transit into supercritical states. The supercritical plasma kinetics may bypass atomic oxygen pathways and directly convert CO 2 into O 2 . This work provides a detailed analysis of plasma-based high-pressure CO 2 conversion, which is of great relevance to future large-scale sustainable carbon capture, utilization, and storage.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Probing Few-Body Nuclear Dynamics via H 3 and He 3 ( e , e ' p ) pn Cross-Section Measurements

Here, we report the first measurement of the (e, e' p) three-body breakup reaction cross sections in helium-3 ( 3 He) and tritium ( 3 H) at large momentum transfer ($\langle{Q^2}\rangle$ ≈ 1.9 (GeV/c) 2 ) and x B > 1 kinematics, where the cross section should be sensitive to quasielastic (QE) scattering from single nucleons. The data cover missing momenta 40 ≤ p miss ≤ 500 MeV/c that, in the QE limit with no rescattering, equals the initial momentum of the probed nucleon. The measured cross sections are compared with state-of-the-art ab-initio calculations. Overall good agreement, within ±20%, is observed between data and calculations for the full p miss range for 3 H and for 100 ≤ p miss ≤ 350 MeV/c for 3 He. Including the effects of rescattering of the outgoing nucleon improves agreement with the data at p miss > 250 MeV/c and suggests contributions from charge-exchange (SCX) rescattering. The isoscalar sum of 3 He plus 3 H, which is largely insensitive to SCX, is described by calculations to within the accuracy of the data over the entire pmiss range. This validates current models of the ground state of the three-nucleon system up to very high initial nucleon momenta of 500 MeV/c.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Finite- and infinite-volume study of DDπ scattering

We develop a comprehensive framework for extracting the pole position and properties of the doubly-charmed tetraquark $T$$^{+}_{cc}$(3875) from lattice QCD data using the relativistic three-particle formalism. This approach incorporates the effect of the one-pion exchange diagram in DDπ and DD ∗ scattering, making it applicable at energies coinciding with the left-hand cut in the partial-wave projected DD ∗ amplitude. We present an example application of this framework to existing lattice QCD data at mπ = 280 MeV. We solve the integral equations describing the DDπ reaction, use LSZ reduction to determine the corresponding DD ∗ amplitude, and find the values of the infinite-volume two- and three-body K matrices that lead to agreement with lattice DD ∗ phase shifts within their uncertainties. Using these K matrices in the three-particle quantization condition, we describe the finite- volume DD ∗ spectrum and find good agreement with the lattice QCD energies. Our results suggest that, at this pion mass, the tetraquark appears as a pair of subthreshold complex poles whose precise location strongly depends on the value of the DDπ three-particle K matrix.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Photo- and Electro-Induced Hadron Production from Nuclei at Jefferson Laboratory

Understanding many-body knockout processes is crucial for nuclear physics, particularly in photo- and electro-induced reactions. In turn, understanding two- and three-body forces, including higher-order forces, is vital for a complete understanding of atoms. We present photo-induced many-proton knockout processes, with multiplicities from 1 to 6, using 12C, CH2, and C4H9OH targets in the g9a FROST dataset. Our analysis covers photon energies from 600 to 4500 MeV, significantly expanding current world data. Comparing our experimental data to the state-of-the-art GiBUU model offers a new challenge in the model’s theoretical description of many-body processes. GiBUU reasonably describes the data at lower photon energies but struggles at higher energies and missing masses, likely due to missing processes, such as initial 3-pion photoproduction. Our results will inform future developments in describing proton knockout processes, indicating GiBUU’s overall reasonable description of many-proton knockout data up to around 2.2 GeV. We also assess various electro-induced reactions using 2D, 12C, and 40Ar targets in the RGM dataset. Our results, obtained at electron beam energies of 2, 4, and 6 GeV, are compared in detail to GENIE and GiBUU, two widely used theory models in neutrino oscillation experiments. Discrepancies between model predictions and experimental data underscore the need for refining the two theoretical models. Despite discrepancies, GiBUU provides a more accurate modelling of electro-induced reactions, especially for 40Ar - crucial for future neutrino oscillation facilities such as DUNE. Understanding the fundamental nuclear physics involved in neutrino-nuclei interactions is essential for reducing the systematic uncertainties in extracting neutrino oscillation parameters. Many-body processes significantly contribute to the background processes observed in neutrino-nuclei interactions, hence the results from both analyses are crucial for developing the theoretical framework for the underlying nuclear physics.

Williams, Rhidian↗
DOE OSTI2024

Measurement of jet substructure in boosted $t\overline{t}$ events with the ATLAS detector using 140 fb -1 of 13 TeV $pp$ collisions

Measurements of the substructure of top-quark jets are presented, using 140 fb -1 of 13 TeV pp collision data recorded with the ATLAS detector at the LHC. Top-quark jets reconstructed with the anti-k t algorithm with a radius parameter R = 1.0 are selected in top-quark pair ($t\overline{t}$) events where one top quark decays semileptonically and the other hadronically, or where both top quarks decay hadronically. The top-quark jets are required to have transverse momentum p T > 350 GeV, yielding large samples of data events with jet p T values between 350 and 600 GeV. One- and two-dimensional differential cross sections for eight substructure variables, defined using only the charged components of the jets, are measured in a particle-level phase space by correcting for the smearing and acceptance effects induced by the detector. The differential cross sections are compared with the predictions of several Monte Carlo simulations in which top-quark pair-production quantum chromodynamic matrix-element calculations at next-to-leading-order precision in the strong coupling constant α S are passed to leading-order parton shower and hadronization generators. The Monte Carlo predictions for measures of the broadness, and also the two-body structure, of the top-quark jets are found to be in good agreement with the measurements, while variables sensitive to the three-body structure of the top-quark jets exhibit some tension with the measured distributions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Ab initio calculation of the β-decay spectrum of 6 He

Here, we calculate the β spectrum in the decay of 6 He using quantum Monte Carlo methods with nuclear interactions derived from chiral effective field theory and consistent weak vector and axial currents. We work at second order in the multipole expansion, retaining terms suppressed by O(q 2 /m$^{2}_{π}$), where q denotes low-energy scales such as the reaction’s Q value or the electron energy, and m π is the pion mass. We go beyond the impulse approximation by including the effects of two-body vector and axial currents. We estimate the theoretical error on the spectrum by using four potential models in the Norfolk family of local two- and three-nucleon interactions, which have different cutoffs, fit two-nucleon data up to different energies, and use different observables to determine the couplings in the three-body force. We find the theoretical uncertainty on the β spectrum, normalized by the total rate, to be well below the permille level, and to receive contributions of comparable size from first- and second-order corrections in the multipole expansion. We consider corrections to the β decay spectrum induced by beyond-standard-model charged-current interactions in the standard model effective field theory, with and without sterile neutrinos, and discuss the sensitivity of the next generation of experiments to these interactions.

6 ≤ A ≤ 19↗
DOE OSTI2020

Shell-model study of calcium isotopes toward their drip line

Here, we report in this paper a study in terms of the nuclear shell model about the location of the calcium isotopes drip line. The starting point is considering the realistic two-body potential derived by Entem and Machleidt within chiral perturbation theory at next-to-next-to-next-to-leading order (N 3 LO), as well as a chiral three-body force at next-to-next-to-leading order (N 2 LO) whose structure and low-energy constants are consistent with the two-body potential. Then we construct the effective single-particle energies and residual interaction needed to diagonalize the shell-model Hamiltonian. The calculated two-neutron separation energies agree nicely with experiment until 56 Ca, which is the heaviest isotope whose mass has been measured, and do not show any sign of two-neutron emission until 70 Ca . We discuss the role of the choice of the model space in determining the neutron drip line, and also the dependence of the results on the parameters of the shell-model Hamiltonian.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Spinbox: tools for many-body quantum systems in a Monte Carlo context

Spinbox is a piece of software that facilitates quantum mechanical calculations relevant to Monte Carlo simulation of atomic nuclei. At the front lines of research on the nuclear many-body problem are a large number of supercomputer-scale simulation codes. These codes produce valuable results but can be hard to understand, especially for those without intimate knowledge of the relevant theoretical methods. Thus, tools that fill pedagogical roles are extremely valuable. Spinbox makes it easy for one to replicate and analyze the computational processes relevant to a Quantum Monte Carlo (QMC) simulation that may be difficult to understand/debug/analyze due to the scale of the corresponding simulation software. Spinbox is written in Python using other state-of-the-art Python modules for numerical calculations. While a number of Python libraries exist that are suited to general quantum many-body calculations, the motivation of Spinbox is quite particular. In Diffusion Monte Carlo methods (DMC, GFMC, AFDMC), the central calculation is the imaginary-time propagation of individual samples of the many-body wavefunction. Although quantum wavefunctions generally must be described by a probability distribution over a basis, DMC imbues particles (within one sample) with classical spatial coordinates. This method is unusual, so other Python packages are typically not set up to do this easily. Furthermore, the software has built-in options for nuclear systems assuming isospin symmetry, which can be set up with other libraries but is a nontrivial process to do so. Features: - numerical representation of samples of the many-body wavefunctions, including tensor-product states (used in AFDMC) - numerical representation of many-body operators, including tensor-product operators: general, spin, imaginary-time propagation, etc. - the correct associated arithmetic and algebra, implemented as class methods - classes for representing realistic nuclear two- and three-body Hamiltonians (e.g. Argonne V18, Illinois NNN) - large-scale parallel integration over random variables, crucial for the AFDMC method My goal is to make this package open source so that anyone may use it and contribute to it, particularly other researchers doing AFDMC calculations

Fox, Jordan↗
DOE OSTI2023

Study of Nucleon Resonance Electroexcitations at Jefferson Lab

Nucleons are the most fundamental three-body systems in Nature. During the last decade crucial progress has been achieved in the exploration of the spectrum and the structure of nucleon resonances. Studies of exclusive ?N, ?+??p, K?, and K? electroproduction off protons on the evolution of the electrocouplings with photon virtuality Q2 using the CLAS spectrometer has revealed a complex interplay between an inner core of three dressed quarks surrounded by an external meson-baryon cloud. CLAS12 is the only facility capable of extending these results on the electrocouplings into the unexplored Q2 range, spanning the domain of quark momenta where ? 50% of hadron mass is generated. The increase of the CEBAF energy and the upgrade of the CLAS12 spectrometer to measure exclusive electroproduction channels at higher luminosity and photon virtualities will offer a unique opportunity to explore how the dominant part of hadron mass and the bound three-quark structure of the resonant states emerge from QCD.

Achenbach, Patrick↗
DOE OSTI2024

Can Third-Body Stabilization of Bimolecular Collision Complexes in Cold Molecular Clouds Happen?

For more than half a century, networks of radiative association, dissociative recombination, and bimolecular reactions have been postulated to drive the low-temperature chemistry of cold molecular clouds. Third-body stabilizations of collision complexes have been assumed to be 'irrelevant' due to short lifetime of such complexes. Here, we conduct crossed molecular beam studies of ground state atomic silicon with diacetylene in combination with electronic structure calculations and microcanonical kinetics models operating under cold molecular cloud conditions. Our combined experimental, electronic structure, and microcanonical kinetics modelling investigations provide compelling evidence that three-body collisions of molecular hydrogen with long-lived reaction intermediates accessed through intersystem crossing are prevalent deep inside molecular clouds. This concept might be exportable to reactions involving polycyclic aromatic hydrocarbons thus affording a versatile machinery to complex organics via third-body stabilizations of bimolecular collision complexes deep inside cold molecular clouds.

astrochemistry↗
DOE OSTI2021

Search for CP violation in $D^{+}_{(s)}\to h^{+} \pi^{0}$ and $D^{+}_{(s)}\to h^{+}\eta$ decays

Searches for CP violation in the two-body decays $D^{+}_{(s)}\to h^{+} \pi^{0}$ and $D^{+}_{(s)}\to h^{+}\eta$ (where h + denotes a π + or K + meson) are performed using pp collision data collected by the LHCb experiment corresponding to either 9 fb –1 or 6 fb –1 of integrated luminosity. The π 0 and η mesons are reconstructed using the e + e – γ final state, which can proceed as three-body decays π 0 → e + e – γ and η → e + e – γ, or via the two-body decays π 0 → γγ and η → γγ followed by a photon conversion. The measurements are made relative to the control modes $D^{+}_{(s)}\to K^{0}_{S}h^{+}$ to cancel the production and detection asymmetries. The CP asymmetries are measured to be A CP (D + →π + π 0 )=(–1.3±0.9±0.6)%, A CP (D + →K + π 0 )=(–3.2±4.7±2.1)%, A CP (D + →π + η)=(–0.2±0.8±0.4)%, A CP (D + →K + η)=(–6±10±4)%, A CP (D$^{+}_{s}$→K + π 0 )=(–0.8±3.9±1.2)%, A CP (D$^{+}_{s}$→π + η)=(0.8±0.7±0.5)%, A CP (D$^{+}_{s}$→K + η)=(0.9±3.7±1.1)%, where the first uncertainties are statistical and the second systematic. These results are consistent with no CP violation and mostly constitute the most precise measurements of A CP in these decay modes to date.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Development of a Transferable Density-Functional Tight-Binding Model for Organic Molecules at the Water/Platinum Interface

A computationally efficient and transferable approach for modeling reactions at metal/water interfaces could significantly accelerate our understanding and ultimately the development of new catalytic transformations, particularly in the context of the emerging field of biomass conversion. Here, we present a parametrization of Pt–X (X = H, O, C) density-functional tight-binding (DFTB) for addressing this need. We first constructed Pt–H, Pt–O, and Pt–C repulsive potential splines. These pairwise parameters were then augmented to include many-body interactions using the Chebyshev Interaction Model for Efficient Simulation (ChIMES). We compare the geometrical and energetic performances of both DFTB and DFTB/ChIMES methods with DFT reference data across a variety of organic molecules at the platinum surface from nanoparticles to single-crystal surfaces. DFTB shows limited transferability between extended crystal surfaces and small nanoparticles. This transferability is significantly improved through the introduction of three-body interactions with Pt in DFTB/ChIMES, which provides consistent results across various systems, with reductions in the RMSD from around 30 kcal/mol in DFTB to around 10 kcal/mol. We demonstrate the stability and reliability of the obtained parameters by performing metadynamic simulations for the adsorption of phenol on Pt(111). We observe that DFTB itself is undersolvating the surface, leading to only one or two chemisorbed water molecules in a c(4 × 6) unit cell. In contrast, DFTB/ChIMES leads to a coverage of about 0.5 ML and successfully captures the chemisorbed mode of phenol at both the solid/liquid and the solid/gas interfaces. Furthermore, in agreement with experimental measurements, the adsorption at the solid/liquid interface is significantly weaker than that at the solid/gas interface. As a result, we highlight that even with DFTB, where we can accumulate dynamics for more than 1 ns for a given system, the simulations are not fully converged.

Adsorption↗
DOE OSTI2020

Nucleon-Nucleon Scattering Up to N5LO in Chiral Effective Field Theory

During the past few decades a large effort has been made toward describing the NN interaction in the framework of chiral Effective Field Theory (EFT). The main idea is to exploit the symmetries of QCD to obtain an effective theory for low energy nuclear systems. In 2003, the first accurate charge-dependent NN potential in this scheme was developed and it has been applied to many ab-initio calculations, opening the possibility to study nuclear systems in a systematic and accurate way. It was shown that the fourth order (N3LO) was necessary and sufficient to describe the NN scattering data with a χ 2 /d.o.f on the order of so-called high precision potentials. However the systematics of chiral EFT also allow to relate two- and many-body interactions in a well-defined way. Since many-body forces make their first appearance at higher order, they are substantially smaller than their two-body counterparts, but may never-the-less be crucial for some processes. Thus, there are observables where they can have a big impact and, for example, there are indications that they solve the long standing Ay puzzle of N-d scattering. The last few years, have also seen substantial progress toward higher orders of chiral EFT which was motivated by the fact that only three-body forces of rather high order may solve some outstanding issues in microscopic nuclear structure and reactions. In this chapter we will review the latest contributions of the authors to development of chiral EFT based potentials up to N4LO as well as first calculations conducted for NN scattering at N5LO.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗