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

A novel fiber-fretting test for tribological characterization of the fiber/matrix interface

Ceramic matrix composites exhibit excellent high temperature properties and are candidate materials for structural applications from jet engines to nuclear fuel cladding. The strength and deformation behavior of these composites are intrinsically tied to the properties of the fiber/matrix interface. This research introduces a novel technique to evaluate the fundamental friction and wear characteristics of these interfaces. A case study on SiC f /PyC/SiC m composites applies in situ fiber fretting to investigate the kinetic friction coefficient and cyclic wear characteristics of the interphase. The experimental methodology is presented in explicit detail to encourage other researchers to explore the capabilities of this technique. Additionally, testing was carried out as a function of fiber roughness, PyC thickness (≈10, 500, 1200 nm), cycle count (10, 100, 1000), and frequency (1, 5, and 10 Hz). The friction coefficient and wear behavior changed with cycle length due to a transition from adhesive to abrasive mechanisms. SEM fractography provided insights to tribological evolution from adhesive sliding to three-body abrasive wear. The structural disorder at the tribo-surface relating to graphitic content and dangling bond density was investigated using Raman spectroscopy. Properties are contextualized with comparison to uniaxial tensile behavior of the parent composite.

36 MATERIALS SCIENCE↗
DOE OSTI2020

SympNets: Intrinsic structure-preserving symplectic networks for identifying Hamiltonian systems

Here, we propose new symplectic networks (SympNets) for identifying Hamiltonian systems from data based on a composition of linear, activation and gradient modules. In particular, we define two classes of SympNets: the LA-SympNets composed of linear and activation modules, and the G-SympNets composed of gradient modules. Correspondingly, we prove two new universal approximation theorems that demonstrate that SympNets can approximate arbitrary symplectic maps based on appropriate activation functions. We then perform several experiments including the pendulum, double pendulum and three-body problems to investigate the expressivity and the generalization ability of SympNets. The simulation results show that even very small size SympNets can generalize well, and are able to handle both separable and non-separable Hamiltonian systems with data points resulting from short or long time steps. In all the test cases, SympNets outperform the baseline models, and are much faster in training and prediction. We also develop an extended version of SympNets to learn the dynamics from irregularly sampled data. This extended version of SympNets can be thought of as a universal model representing the solution to an arbitrary Hamiltonian system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Exploring the NΛ–NΣ coupled system with high precision correlation techniques at the LHC

The interaction of Λ and Σ hyperons (Y) with nucleons (N) is strongly influenced by the coupled-channel dynamics. Due to the small mass difference of the NΛ and NΣ systems, the sizable coupling strength of the NΣ↔NΛ processes constitutes a crucial element in the determination of the NΛ interaction. In this letter we present the most precise measurements on the interaction of pΛ pairs, from zero relative momentum up to the opening of the NΣ channel. The correlation function in the relative momentum space for pΛ ⊕$\bar{p}$$\bar{Λ}$ pairs measured in high-multiplicity triggered pp collisions at $\sqrt{s}$=13 TeV at the LHC is reported. The opening of the inelastic NΣ channels is visible in the extracted correlation function as a cusp-like structure occurring at relative momentum k*=289 MeV/c. This represents the first direct experimental observation of the NΣ↔NΛ coupled channel in the pΛ system. The correlation function is compared with recent chiral effective field theory calculations, based on different strengths of the NΣ↔NΛ transition potential. A weaker coupling, as possibly supported by the present measurement, would require a more repulsive three-body NNΛ interaction for a proper description of the Λ in-medium properties, which has implications on the nuclear equation of state and for the presence of hyperons inside neutron stars.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Quenching of the quantum p -strength in light exotic nuclei

One-nucleon spectroscopic overlaps, their strength or spectroscopic factors (SFs), and nucleon removal cross sections for light nuclei in the mass range A≤12 were evaluated using the fully correlated Quantum Monte Carlo (QMC) wave functions (WFs). Harder (AV18+UX) and softer (NV2+3) bare interactions were used providing consistent results. The SFs were also taken from simple Shell Model (SM) calculations. This structure information was incorporated in the standard three-body Faddeev/Alt-Grassberger-Sandhas reaction formalism to evaluate the (p,pN) cross sections. The results further our understanding of the quenching of the quantum p-shell strength obtained from structure and reactions. We have found the ratios of the total QMC sums of SFs to the SM ones to be uniform and ∼ 3/4. The corresponding ratios of the sums Below Particle Threshold (BPT) of SFs, as well as of total cross sections, deviate strikingly from the uniform trend in some special cases. We find these ratios to be close to unit for 9 Li → 8 Li + n and 9 C → 8 B + p. In contrast, they are strongly reduced for 11 C → 10 C + n and 11 B → 10 Be + p mirror transitions, resulting from the quenching of the strength for low-lying 2 + final state transitions. The QMC theoretical cross section BPT for 11 C(p,pn) is about two times smaller than the experimental data.

(p,pN) reactions↗
DOE OSTI2025

Femtoscopic study of the proton-proton and proton-deuteron systems in heavy-ion collisions at the LHC

This work reports femtoscopic correlations of p – – p ($\overline{p}$ – – $\overline{p}$) and p – – d ($\overline{p}$ – – $\overline{d}$) pairs measured in Pb–Pb collisions at center-of-mass energy per nucleon $\sqrt{s_{NN}}$ = 5.02 TeV in the ALICE Collaboration. A fit to the measured proton-proton correlation functions allows one to extract the dependence of the nucleon femtoscopic radius of the particle-emitting source on the pair transverse mass (m T ) and on the average charge particle multiplicity $\langle$dN ch /dη$\rangle$ 1/3 for three centrality intervals (0–10%, 10 – – 30 %, 30 – – 50 %). In both cases, the expected power-law and linear scalings are observed, respectively. The measured p–d correlations can be described by both two- and three-body calculations, indicating that the femtoscopy observable is not sensitive to the short-distance features of the dynamics of the p-(p-n) system, due to the large inter-particle distances in Pb–Pb collisions at the LHC. Indeed, in this study, the minimum measured femtoscopic source sizes for protons and deuterons have a minimum value at 2.73$^{+0.05}_{–0.05}$ and 3.10$^{+1.04}_{–0.86}$ fm, respectively, for the 30–50% centrality collisions. Moreover, the m T -scaling obtained for the p–p and p–d systems is compatible within 1σ of the uncertainties. These findings provide new input for fundamental studies on the production of light (anti)nuclei under extreme conditions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Quasielastic $\overrightarrow{^{3}He}$ $(\overrightarrow{e},e′)$ asymmetry in the threshold region with a longitudinally polarized target

We report on the measurements of the double-spin asymmetry from electron- 3 He scattering in the threshold region of two- and three-body breakup of 3 He for 4-momentum transfer values of 0.1 and 0.2 (GeV/c) 2 . The measurement was done at Jefferson Lab using a polarized 3 He target with its spin aligned parallel to the beam direction. This asymmetry is sensitive to a combination of spin-dependent response functions not previously measured and may highlight differences in two body forces, meson-exchange currents and final-state interactions. The results of this measurement serve as a test of our understanding of few-body systems. When compared with calculations from plane wave impulse approximation and Faddeev calculations, we found that the latter, which use modern nuclear potentials and prescriptions for meson-exchange currents, demonstrate an overall good agreement with data.

Few-body nuclei↗
DOE OSTI2025

Dynamical coupled-channel models for hadron dynamics

Dynamical coupled-channel (DCC) approaches parametrize the interactions and dynamics of two and more hadrons and their response to different electroweak probes. The inclusion of unitarity, three-body channels, and other properties from scattering theory allows for a reliable extraction of resonance spectra and their properties from data. Here, we review the formalism and application of the ANL-Osaka, the Juelich-Bonn-Washington, and other DCC approaches in the context of light baryon resonances from meson, (virtual) photon, and neutrino-induced reactions, as well as production reactions, strange baryons, light mesons, heavy meson systems, exotics, and baryon-baryon interactions. Finally, we also provide a connection of the formalism to study finite-volume spectra obtained in Lattice QCD, and review applications involving modern statistical and machine learning tools.

Amplitude analysis↗
DOE OSTI2022

Laser-driven ionization mechanisms of aluminum for single particle aerosol mass spectrometry

Single particle aerosol mass spectrometry (SPAMS), an analytical technique for measuring the size and composition of individual micron-scale particles, is capable of analyzing atmospheric pollutants and bioaerosols much more efficiently and with more detail than conventional methods which require the collection of particles onto filters for analysis in the laboratory. Despite SPAMS’ demonstrated capabilities, the primary mechanisms of ionization are not fully understood, which creates challenges in optimizing and interpreting SPAMS signals. In this paper, we present a well-stirred reactor model for the reactions involved with the laser-induced vaporization and ionization of an individual particle. The SPAMS conditions modeled in this paper include a 248 nm laser which is pulsed for 8 ns to vaporize and ionize each particle in vacuum. The ionization of 1 μm, spherical Al particles was studied by approximating them with a 0-dimensional plasma chemistry model. The primary mechanism of absorption of the 248 nm photons was pressure-broadened direct photoexcitation to Al(y 2 D). Atoms in this highly excited state then undergo superelastic collisions with electrons, heating the electrons and populating the lower energy excited states. We found that the primary ionization mechanism is electron impact ionization of various excited state Al atoms, especially Al(y 2 D). Because the gas expands rapidly into vacuum, its temperature decreases rapidly. The rate of three-body recombination (e - + e - + Al + → Al + e - ) increases at low temperature, and most of the electrons and ions produced recombine within several μs of the laser pulse. The importance of the direct photoexcitation indicates that the relative peak heights of different elements in SPAMS mass spectra may be sensitive to the available photoexcitation transitions. We also discuss the effects of laser intensity, particle diameter, and expansion dynamics.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Size quantization of an exciton: A toy model of the “dead layer”

In this paper, size-quantization levels of an exciton in large nanocrystals is studied theoretically. For the nanocrystal size, L , much bigger than the Bohr radius, a B , the level positions do not depend on a B . The correction to the levels in a small parameter a B / L depends on the reflection phase of the exciton from the boundary. Calculation of this phase constitutes a three-body problem: electron, hole, and the boundary. This calculation can be performed analytically in the limit when the hole is much heavier than the electron. Physically, a slow motion of the hole towards the boundary takes place in the effective potential created by the fast motion of the electron orbiting the hole and touching the boundary. As a result, the hole is reflected before reaching the boundary. The distance of the closest approach of the hole to the boundary (the dead layer) exceeds a B parametrically .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

The characterization of wear-causing particles and silica sand in particular

Erodants and abradants are two types of wear-causing particles (WCP). Despite their effects on the severity and rate of wear, WCP are inadequately characterized in a surprising number of publications, especially those involving impingement erosion and three-body abrasion. That shortfall makes it difficult to correlate features of WCP with the details of worn surfaces, or to develop wear models that account for those features. It is argued that the documentation of WCP should go beyond simply reporting their composition and mean particle size. In late 1985, ASTM Committee G2 on Wear and Erosion established a task group on the characterization of WCP. At that time, image analysis was slower and less sophisticated than it is today. While that ASTM task group failed to produce a consensus standard, computerized particle characterization methods were developing in fields other than tribology, fields like geoscience, heavy sand mining, materials processing, and pharmaceuticals. A notable exception to this is ferrography, which is a widely-used diagnostic for lubricated tribosystems. In the context of dry wear, it is useful to identify which features of WCP would be beneficial to document, and to identify some techniques and scales of detail appropriate to particular tribosystems. In this paper, examples are presented for the morphological and mechanical characterization of silica sand grains, as prompted by a triboanalysis of the abrasive and erosive wear of biomass pre-processing equipment. The authors propose a minimum level of documentation for WCP, one that can enrich tribosystem analysis both in laboratory tests and field studies.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Permutationally Invariant Polynomial Expansions with Unrestricted Complexity

A general strategy is presented for constructing and validating permutationally invariant polynomial (PIP) expansions for chemical systems of any stoichiometry. Demonstrations are made for three categories of gas-phase dynamics and kinetics: collisional energy-transfer trajectories for predicting pressure-dependent kinetics, three-body collisions for describing transient van der Waals adducts relevant to atmospheric chemistry, and nonthermal reactivity via quasiclassical trajectories. In total, 30 systems are considered with up to 15 atoms and 39 degrees of freedom. Permutational invariance is enforced in PIP expansions with as many as 13 million terms and 13 permutationally distinct atom types by taking advantage of petascale computational resources. The quality of the PIP expansions is demonstrated through the systematic convergence of in-sample and out-of-sample errors with respect to both the number of training data and the order of the expansion, and these errors are shown to predict errors in the dynamics for both reactive and nonreactive applications. Here, the parallelized code distributed as part of this work enables the automation of PIP generation for complex systems with multiple channels and flexible user-defined symmetry constraints and for automatically removing unphysical unconnected terms from the basis set expansions, all of which are required for simulating complex reactive systems.

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

Parallel Implementation of Nonadditive Gaussian Process Potentials for Monte Carlo Simulations

A strategy is presented to implement Gaussian process potentials in molecular simulations through parallel programming. Attention is focused on the three-body nonadditive energy, though all algorithms extend straightforwardly to the additive energy. The method to distribute pairs and triplets between processes is general to all potentials. Results are presented for a simulation box of argon, including full box and atom displacement calculations, which are relevant to Monte Carlo simulation. Data on speed-up are presented for up to 120 processes across four nodes. A 4-fold speed-up is observed over five processes, extending to 20-fold over 40 processes and 30-fold over 120 processes.

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

Graph-EAM: An Interpretable and Efficient Graph Neural Network Potential Framework

The development of deep learning interatomic potentials has enabled efficient and accurate computations in quantum chemistry and materials science, circumventing computationally expensive ab initio calculations. However, the huge number of learnable parameters in deep learning models and their complex architectures hinder physical interpretability and affect the robustness of the derived potential. In this work, we propose graph-EAM, a lightweight graph neural network (GNN) inspired by the empirical embedded atom method to model the interatomic potential of single-element structures. Four material systems: platinum, niobium, silicon, and amorphous-carbon, for which quantum simulation data sets are publicly available, are examined to demonstrate that graph-EAM can achieve high energy and force prediction accuracy-comparable or better than existing state-of-the-art machine learning models-with much fewer parameters. It is also shown that the explicit inclusion of the angular information via three-body atomic density increases the prediction accuracy. In conclusion, the accuracy and efficiency of potentials obtained from graph-EAM can help accelerate the molecular dynamics simulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Machine Learning a Simple Interpretable Short-Range Potential for Silica

A wide array of models, spanning from computationally expensive ab initio methods to a spectrum of force-field approaches, have been developed and employed to probe silica polymorphs and understand growth processes and atomic-level dynamical transitions in silica. However, the quest for a model capable of making accurate predictions with high computational efficiency for various silica polymorphs is still ongoing. Recent developments in short-range machine-learned models, such as GAP and NNPScan, have shown promise in providing reasonable descriptions of silica, but their computational cost remains high compared to force fields such as BKS which are based on simple interpretable functional forms. Here, in this study, we build on the recent success of our reinforcement learning (RL) workflow to derive a new set of optimal parameters for a promising short-range BKS-based model proposed by Soules. We use RL to navigate the eight-dimensional parameter space of the Soules potential using an experimental training data set that includes both local and global structural features from approximately 21 experimentally realized silica polymorphs, including high density phases and porous zeolites. We compare the performance of our machine-learned ML-Soules model with other high quality models including our recent machine-learned parametrization of BKS (ML-BKS), a machine-learned potential (GAP), as well as predictions of ab initio calculations with the highly fidelity SCAN functional. The ML-Soules accurately captures the relative energetic ordering of various polymorphs as well as their structural features at a significantly reduced computational expense. The ML-Soules model also reasonably captures the structure, density, and elastic constants of quartz, as well as metastable silica polymorphs. We further discuss the limitations of the Soules functional form and propose potential enhancements, including the incorporation of additional three-body terms and/or the utilization of different short-ranged functional forms to achieve greater accuracy for both global and local features in the modeling of silica while retaining low computational cost.

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

Polynomial Scaling Localized Active Space Unitary Selective Coupled Cluster Singles and Doubles

We present a polynomial-scaling algorithm for the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. In this approach, cluster excitations are selected based on a threshold ϵ determined by the absolute gradients of the LAS-UCCSD energy with respect to cluster amplitudes. Using the generalized Wick’s theorem for multireference wave functions, we derive the gradient expression as a polynomial function of one-, two-, and three-body reduced density matrices and 1- and 2-electron integrals, valid for any multireference wave function. The resulting gradient implementation exhibits a memory scaling of 𝒪(N 6 ), with N spin orbitals in the combined active space of all fragments. The variational quantum eigensolver is used to optimize the selected cluster excitations on a quantum simulator. Furthermore, by plotting the energy error, defined as the difference between the LAS-USCCSD and corresponding CASCI energies, against the inverse cluster amplitude selection threshold (ϵ –1 ) for polyene chains containing 2 to 5 π-bond units, we establish a relationship between the energy error and the threshold. To further validate the accuracy of LAS-USCCSD, we computed the cis–trans isomerization energy of stilbene (a 20-qubit system) and the magnetic coupling constant of the tris-hydroxo-bridged chromium dimer [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ (evaluated as both 12- and 20-qubit systems) using the Qiskit-Qulacs simulator. Assessing such examples is important to determine the practical feasibility of quantum simulations for chemically realistic systems. Toward this goal, with the LAS-USCCSD algorithm we estimated the quantum resources required for simulating an active space of (30e,22o) in [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ , a size that remains beyond the reach of current quantum simulators for accurate treatment.

Algorithms↗
DOE OSTI2026

Hierarchical Reinforcement Learning of a Short-Range Bond-Order Potential for Silica: Analytic Embedding of Coordination with Classical Efficiency

Reinforcement learning (RL) has recently emerged as a data-efficient strategy to parametrize short-range interatomic potentials. Building on our past RL optimization of pairwise silica models, we extend the framework to a bond-order (Tersoff-type) potential that provides an analytic embedding of local coordination through a three-body term. A hierarchical RL workflow combining continuous-action Monte Carlo Tree Search and property-based rewards efficiently explores the 26-dimensional parameter space, sequentially optimizing lattice parameters, densities, angles, and cohesive energies of 21 silica polymorphs. The resulting models, Q-Tersoff and ML-Tersoff, reproduce the energetic ordering of low-energy phases and capture the angular correlations and amorphous structure factors of silica with improved fidelity over pairwise force fields, while remaining orders of magnitude faster than high-dimensional machine-learned potentials. Both models underperform for elastic constants and high-energy frameworks, delineating the limits of the current analytic form. The approach establishes a general and interpretable route to angle-aware, short-range potentials that bridge physics-based and machine-learned descriptions of silicate materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Accelerating Ab Initio Simulation via Nested Monte Carlo and Machine Learned Reference Potentials

As a corollary of the rapid advances in computing, ab initio simulation is playing an increasingly important role in modeling materials at the atomic scale. Two strategies are possible, ab initio Monte Carlo (AIMC) and molecular dynamics (AIMD) simulation. The former benefits from exact sampling from the correct thermodynamic distribution, while the latter is typically more efficient with its collective all-atom coordinate updates. In this study, using a relatively simple test model comprised of helium and argon, we show that AIMC can be brought up to, and even above, the performance levels of AIMD via a hybrid nested sampling/machine learning (ML) strategy. Here, ML provides an accurate classical reference potential (up to three-body explicit interactions) that can pilot long collective Monte Carlo moves that are accepted or rejected in toto à la nested Monte Carlo (NMC); this is in contrast to the single move nature of a naive implementation. Our proposed method only requires a small up front expense from evaluating the ab initio energies and forces of (100) random configurations for training. Importantly, our method does not totally rely on the trained, assuredly imperfect, interaction. We show that high performance and exact sampling at the desired level of theory can be realized even when the trained interaction has appreciable differences from the ab initio potential. Remarkably, at the highest levels of performance realized via our approach, a pair of statistically uncorrelated atomic configurations can be generated with (1) ab initio calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Differentiating and Quantifying Gas-Phase Conformational Isomers Using Coulomb Explosion Imaging

Conformational isomerism plays a crucial role in defining the physical and chemical properties and biological activity of molecules ranging from simple organic compounds to complex biopolymers. However, it is often a significant challenge to differentiate and separate these isomers experimentally as they can easily interconvert due to their low rotational energy barrier. Here, we use the momentum correlation of fragment ions produced after inner-shell photoionization to distinguish conformational isomers of 1,2-dibromoethane (C 2 H 4 Br 2 ). We demonstrate that the three-body breakup channel, C 2 H 4 + + Br + + Br + , contains signatures of both sequential and concerted breakup, which are decoupled to distinguish the geometries of two conformational isomers and to quantify their relative abundance. The sensitivity of our method to quantify these yields is established by measuring the relative abundance change with sample temperature, which agrees well with calculations. Finally, our study paves the way for using Coulomb explosion imaging to track subtle molecular structural changes.

Photoionization↗