Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “particle data”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Self-consistent analysis for the η c → γ γ process

The next-to-next-to-leading-order (NNLO) perturbative QCD (pQCD) predictions for both the decay width and the transition form factor in the η c → γ γ process, based on nonrelativistic QCD (NRQCD), deviate from precise experimental measurements. These significant discrepancies have cast doubt on the applicability of NRQCD to charmonium processes. In this paper, we analyze the η c → γ γ process by applying the principle of maximum conformality (PMC), a systematic method for eliminating renormalization scheme-and-scale ambiguities. The PMC renormalization scales are determined by absorbing the nonconformal β -terms that govern the behavior of the QCD running coupling via the renormalization group equation. We obtain the PMC scale Q ⋆ = 4.49 m c for the η c → γ γ decay width. Even after using the PMC method, the convergence of the pQCD series is still poor, which indicates the importance of uncalculated next-to-next-to-next-to-leading-order and higher-order terms. The resulting value for Γ η c → γ γ is in agreement with the Particle Data Group’s reported value of Γ η c → γ γ = 5.1 ± 0.4 keV within the bounds of uncertainties. Moreover, the transition form factor obtained using the PMC is also in good agreement with precise experimental measurements. The application of the PMC suggests a potential resolution to η c → γ γ puzzle and supports the applicability of NRQCD to charmonium processes. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Performance of the NA62 trigger system

The NA62 experiment at CERN targets the measurement of the ultra-rare $K^{+}→π^{+}ν\bar{ν}$ decay, and carries out a broad physics programme that includes probes for symmetry violations and searches for exotic particles. Data were collected in 2016–2018 using a multi-level trigger system, which is described highlighting performance studies based on 2018 data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

In-pixel integration of signal processing and AI/ML based data filtering for particle tracking detectors

We present the first physical realization of in-pixel signal processing with integrated AI-based data filtering for particle tracking detectors. Building on prior work that demonstrated a physics-motivated edge-AI algorithm suitable for ASIC implementation, this work marks a significant milestone toward intelligent silicon trackers. Our prototype readout chip performs real-time data reduction at the sensor level while meeting stringent requirements on power, area, and latency. The chip is taped-out in 28nm TSMC CMOS bulk process, which has been shown to have sufficient radiation hardness for particle experiments. This development represents a key step toward enabling fully on-detector edge AI, with broad implications for data throughput and discovery potential in high-rate, high-radiation environments such as the High-Luminosity LHC.

Parpillon, Benjamin [Fermilab; Illinois U., Chicag↗
DOE OSTI2024

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Measurement of the production cross section of prompt $ {\Xi}_{\mathrm{c}}^0 $ baryons at midrapidity in pp collisions at $ \sqrt{s} $ = 5.02 TeV

The transverse momentum (p T ) differential cross section of the charm-strange baryon $ {\Xi}_{\mathrm{c}}^0 $ is measured at midrapidity (|y| < 0.5) via its semileptonic decay into e + Ξ - ν e in pp collisions at $\sqrt{s}$ = 5.02 TeV with the ALICE detector at the LHC. The ratio of the p T -differential $ {\Xi}_{\mathrm{c}}^0 $-baryon and D 0 -meson production cross sections is also reported. The measurements are compared with simulations with different tunes of the PYTHIA 8 event generator, with predictions from a statistical hadronisation model (SHM) with a largely augmented set of charm-baryon states beyond the current lists of the Particle Data Group, and with models including hadronisation via quark coalescence. The p T -integrated cross section of prompt $ {\Xi}_{\mathrm{c}}^0 $-baryon production at midrapidity is also reported, which is used to calculate the baryon-to-meson ratio $ {\Xi}_{\mathrm{c}}^0 $/D 0 = 0.20 ± 0.04 (stat.)$^{+0.08}_{-0.07}$ (syst.). These results provide an additional indication of a modification of the charm fragmentation from e + e - and e - p collisions to pp collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Domain decomposition in the GPU-accelerated Shift Monte Carlo code

The GPU solver within the Shift continuous-energy Monte Carlo neutron transport code has been extended to provide domain decomposition in addition to domain replication to enable the solution of problems with memory requirements exceeding the capacity of a single GPU. The strategy follows the Multiple Set, Overlapping Domain (MSOD) approach that is used in Shift’s CPU solver and integrates into the event-based algorithm used for Shift’s GPU solver. Furthermore, the ability to assign processors to spatial domains non-uniformly has been maintained. In this work, two different approaches for communicating particle data between domains are considered, and multiple criteria for load balancing problems have been investigated. Numerical results are presented for both fresh and depleted small modular nuclear reactor (SMR) cores. A parallel efficiency of approximately 80% was achieved with up to 16 spatial domains measured relative to full domain replication. A scaling study on the Summit supercomputer demonstrates a weak scaling parallel efficiency of over 90% on over 24000 GPUs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Code modernization strategies for short-range non-bonded molecular dynamics simulations

Modern HPC systems are increasingly relying on greater core counts and wider vector registers. Thus, applications need to be adapted to fully utilize these hardware capabilities. One class of applications that can benefit from this increase in parallelism are molecular dynamics simulations. In this paper, we describe our efforts at modernizing the ESPResSo++ simulation package for molecular dynamics by restructuring its particle data layout for efficient memory accesses and applying vectorization techniques to benefit the calculation of short-range non-bonded forces, which results in an overall three times speedup and serves as a baseline for further optimizations. We also implement fine-grained parallelism for multi-core CPUs through HPX, a C++ runtime system which uses lightweight threads and an asynchronous many-task approach to maximize concurrency. Our goal is to evaluate the performance of an HPX-based approach compared to the bulk-synchronous MPI-based implementation. This requires the introduction of an additional layer to the domain decomposition scheme that defines the task granularity. On spatially inhomogeneous systems, which impose a corresponding load-imbalance in traditional MPI-based approaches, we demonstrate that by choosing an optimal task size, the efficient work-stealing mechanisms of HPX can overcome the overhead of communication resulting in an overall 1.4 times speedup compared to the baseline MPI version.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Indra: a public computationally accessible suite of cosmological N -body simulations

ABSTRACT Indra is a suite of large-volume cosmological N-body simulations with the goal of providing excellent statistics of the large-scale features of the distribution of dark matter. Each of the 384 simulations is computed with the same cosmological parameters and different initial phases, with 10243 dark matter particles in a box of length 1 h−1 Gpc, 64 snapshots of particle data and halo catalogues, and 505 time-steps of the Fourier modes of the density field, amounting to almost a petabyte of data. All of the Indra data are immediately available for analysis via the SciServer science platform, which provides interactive and batch computing modes, personal data storage, and other hosted data sets such as the Millennium simulations and many astronomical surveys. We present the Indra simulations, describe the data products and how to access them, and measure ensemble averages, variances, and covariances of the matter power spectrum, the matter correlation function, and the halo mass function to demonstrate the types of computations that Indra enables. We hope that Indra will be both a resource for large-scale structure research and a demonstration of how to make very large data sets public and computationally accessible.

Falck, Bridget↗
DOE OSTI2023

The evolution of HCO + in molecular clouds using a novel chemical post-processing algorithm

Modelling the chemistry of molecular clouds is critical to accurately simulating their evolution. To reduce computational cost, 3D simulations generally restrict their chemistry to species with strong heating and cooling effects. Time-dependent information about the evolution of other species is therefore often neglected. We address this gap by post-processing tracer particles in the SILCC-Zoom molecular cloud simulations. Using a chemical network of 39 species and 301 reactions (including freeze-out of CO and H 2 O) and a novel algorithm to reconstruct a density grid from sparse tracer particle data, we produce time-dependent density distributions for various species. We focus upon the evolution of HCO + , which is a critical formation reactant of CO but is not typically modelled on the fly. We find that ∼ 90 per cent of the HCO + content of the cold molecular gas forms in situ around n HCO + ∼ 10 3 –10 4 cm −3 , over a time-scale of approximately 1 Myr. The remaining ∼ 10 per cent forms at high extinction sites, with minimal turbulent mixing out into the less dense gas. We further show that the dominant HCO + formation pathway is dependent on the visual extinction, with the reaction H 3 + + CO contributing 90 per cent of the total HCO + production above A V, 3D = 3. We produce the very first maps of the HCO + column density, N(HCO + ), and show that it reaches values as high as 10 15 cm −2 . We find that 50 per cent of the HCO + mass is located within AV ∼ 10–30 in a density range of 10 3.5 –10 4.5 cm −3 . Our maps of N(HCO + ) are shown to be in good agreement with recent observations of the W49A star-forming region.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Precise 𝐾* ⁡(892) → 𝐾⁡𝜋 branching fractions

Although discovered more than sixty years ago, direct measurement of the 𝐾*⁡(892) → 𝐾⁡𝜋 branching fractions is a formidable challenge that has not been attempted. Typically they are assumed to obey the isospin limit in hundreds of particle data measurements. We show that an abundance of recent amplitude analyses and other data, however, enables recovery of the ratios ℬ⁡(𝐾 *+ → 𝐾 + ⁢𝜋 0 )/ℬ⁡(𝐾 *+ → 𝐾$^0_S$ ⁢𝜋 + ) and 4⁢ℬ⁡(𝐾 *0 → 𝐾$^0_S$ ⁢𝜋 0 )/ℬ⁡(𝐾 *0 → 𝐾 +⁢ 𝜋 − ) at ∼5% precision.

Gligorov, Vladimir V. [Sorbonne Université, Paris ↗
DOE OSTI2025

Generalized susceptibilities and the properties of charm degrees of freedom across the QCD crossover temperature

We study the generalized charm susceptibilities in 2+1 flavor QCD on the lattice at several lattice spacings. We show that, below the chiral crossover, these susceptibilities are well described by the hadron resonance gas (HRG) model if charmed hadrons not listed in tables of the Particle Data Group are included. However, the HRG description abruptly breaks down just above the chiral crossover. To understand this, we use a model for the charm pressure in which it is expressed as the sum of partial pressures from charmed baryons, charmed mesons, and charm quarks. We present continuum estimates of these partial pressures and find that, while the partial pressures of charmed mesons and baryons drop below their respective HRG predictions, the charm quark pressure becomes nonzero above the chiral crossover.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

symfind : Addressing the Fragility of Subhalo Finders and Revealing the Durability of Subhalos

Abstract A major question in ΛCDM is what this theory actually predicts for the properties of subhalo populations. Subhalos are difficult to accurately simulate and to find within simulations, and this propagates into uncertainty in theoretical predictions for satellite galaxies. We present Symfind , a new particle-tracking-based subhalo finder, and demonstrate that it can track subhalos to orders-of-magnitude lower masses than commonly used halo-finding tools, with a focus on Rockstar and consistent-trees . These longer survival times mean that at a fixed peak subhalo mass, we find ≈ 15%–40% more subhalos within the virial radius, R vir , and ≈35%–120% more subhalos within R vir /4 in the Symphony dark-matter-only simulation suite. More subhalos are found as the resolution is increased, in contrast to the Rockstar halo finder, which appears to be converged at smaller subhalo counts. We perform extensive numerical testing. In agreement with idealized simulations, we show that the v max , the maximum circular velocity, is systematically biased low until high resolutions ( n peak ≳ 3 × 10 4 ) are achieved, but that mass loss itself can be resolved at much more modest resolutions ( n peak ≳ 4 × 10 3 ). We show that Rockstar converges to false solutions for the mass function, radial distribution, and disruption masses of subhalos. We argue that our new method can trace resolved subhalos until the point of typical galaxy disruption without invoking post hoc orphan modeling. We outline a concrete set of steps for determining whether other subhalo finders meet the same criteria. We publicly release Symfind catalogs and particle data for the Symphony simulation suite at http://web.stanford.edu/group/gfc/symphony .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The Dynamics of Debris Disk Creation in Neutron Star Mergers

Abstract The detection of GW170817/AT2017gfo inaugurated an era of multimessenger astrophysics, in which gravitational-wave and multiwavelength photon observations complement one another to provide unique insight into astrophysical systems. A broad theoretical consensus exists, in which the photon phenomenology of neutron star mergers largely rests upon the evolution of the small amount of matter left on bound orbits around the black hole or massive neutron star remaining after the merger. Because this accretion disk is far from inflow equilibrium, its subsequent evolution depends very strongly on its initial state, yet very little is known about how this state is determined. Using both snapshot and tracer particle data from a numerical relativity/MHD simulation of an equal-mass neutron star merger that collapses to a black hole, we show how gravitational forces arising in a nonaxisymmetric, dynamical spacetime supplement hydrodynamical effects in shaping the initial structure of the bound debris disk. The work done by hydrodynamical forces is ∼10 times greater than that due to time-dependent gravity. Although gravitational torques prior to remnant relaxation are an order of magnitude larger than hydrodynamical torques, their intrinsic sign symmetry leads to strong cancellation; as a result, hydrodynamical and gravitational torques have a comparable effect. We also show that the debris disk’s initial specific angular momentum distribution is sharply peaked at roughly the specific angular momentum of the merged neutron star’s outer layers, a fewr g c, and identify the regulating mechanism.

Astronomy & Astrophysics↗
DOE OSTI2024

A Measurement of the Eta Meson Radiative Decay Width via the Primakoff Effect

The ? meson is an interesting tool to study fundamental symmetries in Quantum Chromodynamics (QCD). In particular, its radiative decay width, ? p? Ñ ??q, is an important quantity that can be predicted in the framework of Chiral Perturbation Theory. A precision measurement of this quantity would provide critical inputs to understanding the mixing of the ? and ?1 mesons and extracting constants with wide-ranging applications in low-energy QCD. This decay width has been measured in the past using two different experimental techniques. The more popular technique utilized e`e´ collisions to produce ? mesons through electromagnetic interactions. Today, the Particle Data Group (PDG) averages the results of five such experiments to obtain their currently-accepted value of the decay width as: 0.515?0.018 keV. However the first measurement of this quantity was obtained from a fixed-target experiment that measured the cross section for photoproduction of ? mesons on a nuclear target via the Primakoff effect. Their result of 0.324?0.046 keV shows strong tension with the average of the collider measurements, motivating a new, high precision measurement using the Primakoff method. For this purpose, the PrimEx-eta experiment was conducted in Hall D of the Thomas Jefferson National Accelerator Facility (Jefferson Lab or JLab). The data is currently being analyzed to measure the differential cross section for the photoproduction of ? mesons on a liquid, 4He target. Preliminary results obtained from the analysis of the first phase of the PrimEx-eta experiment show reasonable agreement with the currently-accepted PDG value of the radiative decay width. However, as will be discussed, there are many challenges to this precision measurement which must be studied before any results can be finalized and compared with previous measurements. In parallel to the ? decay width measurement, the PrimEx-eta experiment measured the total cross section for the fundamental, Quantum Electrodynamics (QED) process of Compton scattering from the atomic electrons inside the target. The results obtained from this measurement are in strong agreement with the next-to-leading order QED calculations, and the total combined uncertainties are below 3% for incident photon energies between 7-10 GeV. In addition to providing the first precision measurement of the total Compton scattering cross section within this beam energy range, this measurement verifies the capability of the PrimEx-eta experimental setup to perform absolute cross section measurements at forward angles, and serves as a reference process for the calibration of systematic uncertainties.

Smith, Andrew↗
DOE OSTI2021

ARM Data for Examining the Ice-Nucleating Particles from SGP Part II (ExINP-SGP II) Aerodynamic Particle Sizer, Condensation Particle Counter, and Meteorological Instrument Data Analysis

Knowledge of airborne particulate matter (PM), especially the particles that have supermicron diameters, is key for understanding ice-nucleating particles (INPs). Supported by the Atmospheric Radiation Measurement (ARM) user facility, we sampled PM and surface soil materials at the Southern Great Plains observatory (SGP; 36&deg; 36&prime; 18&Prime; N, 97&deg; 29&prime; 6&Prime; W) to systematically compare the INP abundance and ice nucleation efficiency of different SGP samples (i.e., airborne versus surface materials). The field campaign, named Examing INP from SGP II (ExINP-SGP II), was conducted from 20 January to 20 April, 2021. Our data analysis products include (1) physical surface sorption characterization data (i.e., BET, pore volume) of two sets of samples -- SGP Soil and USDA, (2) X-ray diffraction spectra for the same two samples, (3) immersion freezing assay-based ice nucleation active mass density data as a function of freezing temperature for both ambient and surface samples, and (4) time-series data of ambient meteorological conditions, concentration particle counter-derived aerosol particle concentration, and aerosol particle sizer measurement during ExINP-SGP II.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Systems and methods for screening particle source manufacturing and development test data

A computing system obtains test data for a particle source. The test data was generated by the particle source when the particle source was caused to emit particles. The test data comprises a first set of measurements of a first type and a second set of measurements of a second type. The computing system applies a data agnostic predictive model to the test data. The data agnostic predictive model is generated without a parametric analysis of variables of the first type and variables of the second type. The data agnostic predictive model outputs, based upon the test data, a value that is indicative of whether or not the test data is abnormal. Based upon the value, the computing system outputs an indication that the particle source was operating sub-optimally when emitting the particles.

Multari, Rosalie A.↗
DOE OSTI2023

Real-time mixed reality display of dual particle radiation detector data

Radiation source localization and characterization are challenging tasks that currently require complex analyses for interpretation. Mixed reality (MR) technologies are at the verge of wide scale adoption and can assist in the visualization of complex data. Herein, we demonstrate real-time visualization of gamma ray and neutron radiation detector data in MR using the Microsoft HoloLens 2 smart glasses, significantly reducing user interpretation burden. Radiation imaging systems typically use double-scatter events of gamma rays or fast neutrons to reconstruct the incidence directional information, thus enabling source localization. The calculated images and estimated ’hot spots’ are then often displayed in 2D angular space projections on screens. By combining a state-of-the-art dual particle imaging system with HoloLens 2, we propose to display the data directly to the user via the head-mounted MR smart glasses, presenting the directional information as an overlay to the user’s 3D visual experience. We describe an open source implementation using efficient data transfer, image calculation, and 3D engine. We thereby demonstrate for the first time a real-time user experience to display fast neutron or gamma ray images from various radioactive sources set around the detector. We also introduce an alternative source search mode for situations of low event rates using a neural network and simulation based training data to provide a fast estimation of the source’s angular direction. Using MR for radiation detection provides a more intuitive perception of radioactivity and can be applied in routine radiation monitoring, education & training, emergency scenarios, or inspections.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗