Search NASA⌕ Search

SEARCH · Search NASA

Results for “Particle identication methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Masked Particle Modeling on Sets: Towards Self-Supervised High Energy Physics Foundation Models

Abstract We propose masked particle modeling (MPM) as a self-supervised method for learning generic, transferable, and reusable representations on unordered sets of inputs for use in high energy physics (HEP) scientific data. This work provides a novel scheme to perform masked modeling based pre-training to learn permutation invariant functions on sets. More generally, this work provides a step towards building large foundation models for HEP that can be generically pre-trained with self-supervised learning and later fine-tuned for a variety of down-stream tasks. In MPM, particles in a set are masked and the training objective is to recover their identity, as defined by a discretized token representation of a pre-trained vector quantized variational autoencoder. We study the efficacy of the method in samples of high energy jets at collider physics experiments, including studies on the impact of discretization, permutation invariance, and ordering. We also study the fine-tuning capability of the model, showing that it can be adapted to tasks such as supervised and weakly supervised jet classification, and that the model can transfer efficiently with small fine-tuning data sets to new classes and new data domains.

Heinrich, Lukas (ORCID:0000000240487584)↗

Electromagnetic Scattering by Spheroidal Volumes of Discrete Random Medium

We use the superposition T-matrix method to compare the far-field scattering matrices generated by spheroidal and spherical volumes of discrete random medium having the same volume and populated by identical spherical particles. Our results fully confirm the robustness of the previously identified coherent and diffuse scattering regimes and associated optical phenomena exhibited by spherical particulate volumes and support their explanation in terms of the interference phenomenon coupled with the order-of-scattering expansion of the far-field Foldy equations. We also show that increasing non-sphericity of particulate volumes causes discernible (albeit less pronounced) optical effects in forward and backscattering directions and explain them in terms of the same interference/multiple-scattering phenomenon.

Electromagnetic scattering; Multi-particle groups;↗

QCD–Gravity Double Copy in Regge Asymptotics: From \(2\rightarrow n\) Amplitudes to Radiation in Shockwave Collisions

This paper discusses multi-particle production in QCD and in gravity at ultrarelativistic energies, their double-copy relations, and strong parallels in emergent shockwave dynamics. Dispersive techniques are applied to derive the BFKL equation for multi-gluon production in Regge asymptotics. Identical methods apply in gravity and are captured by a gravitational Lipatov equation. The building blocks in both cases are Lipatov vertices and reggeized propagators satisfying double-copy relations; in gravity, Weinberg’s soft theorem is recovered as a limit of the Lipatov framework. BFKL evolution in QCD generates wee parton states of maximal occupancy characterized by an emergent semi-hard saturation scale. Renormalization group equations in the Color Glass Condensate (CGC) EFT describe wee parton correlations and their rapidity evolution. A shockwave picture of deeply inelastic scattering and hadron–hadron collisions follows, with multi-particle production described by Cutkosky’s rules in strong time-dependent fields. Gluon radiation in the CGC EFT has a double copy in gravitational shockwave collisions, with a similar correspondence applicable between gluon and graviton shockwave propagators. Possible extensions of this semi-classical double copy are outlined for computing multi-particle production in gravitational shockwave collisions, self-force and tidal contributions, and classical and quantum noise in the focusing of geodesics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Assessing Metal Ion Assignment Accuracy in Protein Data Bank Models via Elemental Spectroscopy

Accurate representation of metal ions in macromolecular structures is critical for chemical interpretation, computational modeling, and machine-learning methods that rely on Protein Data Bank (PDB) entries. However, the elemental identity of metals modeled in crystallographic structures is often inferred indirectly and rarely validated experimentally. Here, we combine Particle Induced X-ray Emission (PIXE) and X-ray Fluorescence Spectroscopy (XRFS) to determine the elemental composition of protein samples used to generate 70 deposited metalloprotein crystal structures. By analyzing the original protein material employed for crystallization, but before the addition of crystallization buffer solutions, we assess whether the modeled metal ions in deposited structures are consistent with experimentally detectable elemental content. We find that in a majority of cases, the metals modeled in the corresponding PDB entries are inconsistent with the metals present in the protein samples before crystallization, or that additional metals are present but not represented in the structural models. Spectroscopic results were integrated with automated crystallographic validation metrics, including real-space Z-difference (RSZD) analysis and systematic rerefinement, to evaluate atomic-number mismatch at metal sites. PIXE and XRFS show strong agreement for dominant elemental signals and provide complementary, scalable approaches for identifying suspect metal assignments. This work does not address physiological or functional metalation but instead highlights a widespread data integrity issue in deposited macromolecular structures, PDB-wide. These results establish an experimentally corroborated link between elemental identity and crystallographic validation metrics, enabling the large-scale detection of chemically inconsistent annotations in structural databases used for computational modeling and machine learning.

Crystallization↗

Solving complex nanostructures with ptychographic atomic electron tomography

Abstract Transmission electron microscopy (TEM) is essential for determining atomic scale structures in structural biology and materials science. In structural biology, three-dimensional structures of proteins are routinely determined from thousands of identical particles using phase-contrast TEM. In materials science, three-dimensional atomic structures of complex nanomaterials have been determined using atomic electron tomography (AET). However, neither of these methods can determine the three-dimensional atomic structure of heterogeneous nanomaterials containing light elements. Here, we perform ptychographic electron tomography from 34.5 million diffraction patterns to reconstruct an atomic resolution tilt series of a double wall-carbon nanotube (DW-CNT) encapsulating a complex ZrTe sandwich structure. Class averaging the resulting tilt series images and subpixel localization of the atomic peaks reveals a Zr 11 Te 50 structure containing a previously unobserved ZrTe 2 phase in the core. The experimental realization of atomic resolution ptychographic electron tomography will allow for the structural determination of a wide range of beam-sensitive nanomaterials containing light elements.

36 MATERIALS SCIENCE↗

Simulations of classical three-body thermalization in one dimension

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

1-dimensional systems↗

Solar system applications of Mie theory and of radiative transfer of polarized light

A theory of the multiple scattering of polarized light is discussed using the doubling method of van de Hulst. The concept of the Stokes parameters is derived and used to develop the form of the scattering phase matrix of a single particle. The diffuse reflection and transmission matrices of a single scattering plane parallel atmosphere are expressed as a function of the phase matrix, and the symmetry properties of these matrices are examined. Four matrices are required to describe scattering and transmission. The scattering matrix that results from the addition of two identical layers is derived. Using the doubling method, the scattering and transmission matrices of layers of arbitrary optical thickness can be derived. The doubling equations are then rewritten in terms of their Fourier components. Computation time is reduced since each Fourier component doubles independently. Computation time is also reduced through the use of symmetry properties.

Whitehill, L. P.↗

A phase-space fluid simulation of a two-component narrow planetary ring - Particle size segregation, edge formation, and spreading rates

The Krook kinetic equation for identical planetary ring particles is presently generalized for the case of two-component systems, and the equations are numerically solved on the basis of Brophy and Esposito's (1989) phase-space CFD method. Attention is given to the simulation results obtained for a two-component narrow ring, in which the large particles are eight times as massive as the small particles. This ring's unconstrained edge dynamics are resolved by the simulation, and are found to exhibit a sharpening that would not have been expected in single-component rings.

Brophy, Thomas G.↗

Comparison of Surface Tension Generation Methods in Smoothed Particle Hydrodynamics for Dynamic Systems

Developing robust numerical models of dynamic surface tension dominated multiphase systems is an ongoing challenge, especially in scenarios with large density and viscosity ratios. This is critical to the design and understanding of various physical and engineering systems, such as fluidized beds, fuel injectors, and drug delivery schemes. Much of the computational work in surface tension dominated multiphase flows has employed the continuum surface force method (CSF) of Brackbill et al. [1], which recasts surface tension from a surface force to a volumetric force that can be imposed in the vicinity of an interface. The CSF method produces accurate results across a variety of systems, however it relies on the identication of surface normals, which can be unreliable under certain conditions. Alternative methods of simulating surface tension have been proposed. Here the advantages and disadvantages of the CSF method in comparison to a pairwise forces (PF) method proposed by Tartakovsky and Panchenko [2] are explored. The CSF and PF methods are used in a smoothed particle hydrodynamics (SPH) framework to model dynamic systems. Results are compared to existing test cases from the literature and to analytic solutions derived from fundamental normal mode behavior of bubbles and droplets. A physical system for which the PF method is more stable and physically appropriate than the CSF method is identied.

smoothed particle hydrodynamics, multiphase flow↗

Mesoscale simulation of granular materials under weak shock compaction–pore size distribution effects

This research established a systematic method to generate various pore-size distributions (PSDs) and studied the effect of PSDs on the shock compaction response of granular materials using two-dimensional mesoscale simulations under identical porosity. Simulations utilized various PSDs for three particle shapes (circle, ellipse, and square). Contacting particle configurations using three PSDs, characterized by spatially uniform distributed pores to heterogeneous distributed pores, and non-contacting particle configurations under a single case of PSD were tested. The PSD of generated particle sets was characterized using coordination number, mean diameter, and bimodality coefficient as statistical metrics. Mesoscale simulations showed that regardless of the conditions of pore distributions, shock compaction of granular materials consistently demonstrates a precursor, shock compaction front, and end. However, the shock compaction velocity of contacting particles was dependent on the PSDs despite the constant initial porosity. The compaction velocity was faster in particle configurations with relatively uniform pore distributions than in heterogeneous pore distributions, which our study demonstrated can be attributed to particle rearrangement during compaction. Circular-shaped particles had high sensitivity in shock compaction response to the various PSDs. Furthermore, a contacting particle configuration tended to propagate the shock compaction wave relatively faster than particles that were in a non-contact configuration. This study established the relative importance of considering PSD as a metric over the coordination number in studies of the shock compaction response of granular materials. Further, insights are provided on the evolving shock substructure to characterize the shock compaction response of granular materials.

36 MATERIALS SCIENCE↗

Devising novel methods for the controlled synthesis with morphology and size control of scintillator materials

Nanoscale scintillators represent an important class of materials underlying the functioning of high resolution, high ‘light yield’ radiation detection equipment, applicable for medicine in addition to security and environmental monitoring. This current study has focused on the controlled synthesis of Cs 2 ZnX 4 (X = Cl, Br) not only due to its very fast and promising radioluminescence decay time but also because no reliable synthesis protocols have been reported for producing nanoscale motifs of this material, to date. Solution-based methods can be used to synthesize products, possessing different, discrete morphologies with the ability to dictate particle size and shape. Specifically, we have demonstrated that a hot injection technique can be used to fabricate nanoscale rods and plates. Using this method, we systematically studied the effect of varying reaction variables, such as (i) reaction temperature, (ii) surfactants, and (iii) reaction time. As synthesized nanorods of Cs 2 ZnCl 4 possessed lengths of 101.3 ± 24.4 nm with widths of 17.3 ± 4.5 nm, whereas Cs 2 ZnBr 4 , created under identical conditions, were characterized by lengths of 206 ± 50 nm and widths of 26.7 ± 12.6 nm. Ligand-assisted re-precipitation (LARP) was used to produce the corresponding micron-sized particles. Under these conditions, we successfully synthesized spindles and spherical particles of Cs 2 ZnCl 4 with sizes of 1.58 ± 0.16 μm and 1.03 ± 0.34 μm, respectively, whereas plates of Cs 2 ZnBr 4 were prepared with average sizes of 0.76 ± 0.21 μm. High resolution synchrotron mapping, using the hard X-ray nanoprobe at NSLS II, was used to confirm the expected spatial distribution of Cs, Zn, and Cl within the as-produced crystals. Nanorods of Cs 2 ZnCl 4 and Cs 2 ZnBr 4 were subsequently tested for scintillation light yield and decay times. Preliminary results showed that Cs 2 ZnCl 4 maintained a light yield of 100–300 photons per MeV with a fast decay component of 2.6 ns. Finally, the Cs 2 ZnBr 4 analogue did not give rise to any perceptible scintillation emission.

36 MATERIALS SCIENCE↗

Representations of relativistic particles of arbitrary spin in Poincaré, Lorentz, and Euclidean covariant formulations of relativistic quantum mechanics

Background: Relativistic treatments of quantum mechanical systems are important for understanding hadronic structure and dynamics at subnucleon scales. Relativistic invariance of a quantum system means that there is an underlying unitary representation of the Poincaré group. This is equivalent to the requirement that the quantum observables (probabilities, expectation values, and ensemble averages) for equivalent measurements performed in different inertial reference frames are identical. Different representations are used in practice, including Poincaré covariant forms of dynamics, representations based on Lorentz covariant wave functions, Euclidean covariant representations, and representations generated by Lorentz covariant fields. Purpose: The purpose of this work is to illustrate the relation between the different equivalent representations of states in relativistic quantum mechanics. Method: The starting point is a description of a particle of mass m and spin j using irreducible representations of the Poincaré group. Since any unitary representation of the Poincaré group can be decomposed into a direct integral of irreducible representations, these are the basic building blocks of any relativistically invariant quantum theory. The equivalence is established by constructing equivalent Lorentz covariant irreducible representations from Poincaré covariant irreducible representations and constructing equivalent Euclidean covariant irreducible representations from Lorentz covariant irreducible representations. Results: Equivalent descriptions for positive mass representations of arbitrary spin are presented in each of these frameworks. Dynamical realizations of the different representations are briefly discussed. Conclusion: Poincaré covariant, Lorentz covariant, and Euclidean covariant realizations of relativistic dynamics are shown to be equivalent by explicitly relating the positive-mass positive-energy irreducible representations of the Poincaré group that appear in the direct integral.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Stochastic GW -GPU: Rapid Quasi-Particle Energies for Molecules beyond 10,000 Atoms

StochasticGW is a code for computing accurate quasi-particle (QP) energies of molecules and material systems in the GW approximation. StochasticGW utilizes the stochastic Resolution of the Identity (sROI) technique to enable a massively parallel implementation with computational costs that scale semilinearly with system size, allowing the method to access systems with tens of thousands of electrons. Here, we introduce a new implementation, StochasticGW-GPU, for which the main bottleneck steps have been ported to GPUs and give substantial performance improvements over previous versions of the code. We showcase the new code by computing band gaps of hydrogenated silicon clusters (Si x H y ) containing up to 10,001 atoms and 35,144 electrons, and we obtain individual QP energies with a statistical precision of better than ±0.03 eV with times-to-solution of less than 1 h.

Thomas, Phillip S. [Lawrence Berkeley National Lab↗

Breaking the curse of dimensionality: Solving configurational integrals for crystalline solids by tensor networks

Accurately evaluating configurational integrals for dense solids remains a central and difficult challenge in the statistical mechanics of condensed systems. Here, we present a tensor network approach that reformulates the high-dimensional configurational integral for identical-particle crystals into a sequence of computationally efficient summations. We represent the integrand as a high-dimensional tensor and apply tensor-train (TT) decomposition together with a custom TT-cross interpolation. This approach circumvents the need to explicitly construct the full tensor. We introduce tailored rank-1 and rank-2 schemes optimized for sharply peaked Boltzmann probability densities, typical for identical-particle crystals. When applied to the calculation of internal energy and pressure-temperature curves for crystalline Cu and Ar at high (GPa) pressures, as well as the alpha-to-beta phase transition diagram of Sn, our method accurately reproduces molecular dynamics simulation results using tight-binding, machine learning, hierarchical interacting particle–neural network, and modified embedded atom method potentials,all within seconds of computation time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Polarized Radiative Transfer of a Cirrus Cloud Consisting of Randomly Oriented Hexagonal Ice Crystals: The 3 x 3 Approximation for Non-Spherical Particles

The reflection and transmission of polarized light for a cirrus cloud consisting of randomly oriented hexagonal columns were calculated by two very different vector radiative transfer models. The forward peak of the phase function for the ensemble-averaged ice crystals has a value of order 6 x 10(exp 3) so a truncation procedure was used to help produce numerically efficient yet accurate results. One of these models, the Vectorized Line-by-Line Equivalent model (VLBLE), is based on the doubling- adding principle, while the other is based on a vector discrete ordinates method (VDISORT). A comparison shows that the two models provide very close although not entirely identical results, which can be explained by differences in treatment of single scattering and the representation of the scattering phase matrix. The relative differences in the reflected I and Q Stokes parameters are within 0.5 for I and within 1.5 for Q for all viewing angles. In 1971 Hansen showed that for scattering by spherical particles the 3 x 3 approximation is sufficient to produce accurate results for the reflected radiance I and the degree of polarization (DOP), and he conjectured that these results would hold also for non-spherical particles. Simulations were conducted to test Hansen's conjecture for the cirrus cloud particles considered in this study. It was found that the 3 x 3 approximation also gives accurate results for the transmitted light, and for Q and U in addition to I and DOP. For these non-spherical ice particles the 3 x 3 approximation leads to an absolute error 2 x 10(exp -6) for the reflected and transmitted I, Q and U Stokes parameters. Hence, it appears to be an excellent approximation, which significantly reduces the computational complexity and burden required for multiple scattering calculations.

Polarization↗

An ultra-fast method for generating synthetic down-scattered neutron data for inertial confinement fusion implosions

In inertial confinement fusion experiments at the National Ignition Facility, asymmetries are probed by a variety of neutron diagnostics, including neutron imaging systems, real-time neutron activation diagnostics (RTNADs), and neutron spectrometers. It is often useful to generate synthetic data based on these diagnostics to validate and tune models. However, current methods of doing so using Monte Carlo particle tracing are time-consuming. In this paper, an ultra-fast method is presented for generating synthetic neutron images, RTNAD data, and spectrometry data using line integrals and 3D convolutions. While it does not contain as much physics as particle tracing codes, it is thousands of times faster and produces nearly identical data. This enables analysis techniques that depend on generating large amounts of synthetic data, which will prove very useful for the study of asymmetries going forward.

Deuterium↗

Comparison of Eulerian Bin and Lagrangian Particle-Based Microphysics in Simulations of Nonprecipitating Cumulus

A single nonprecipitating cumulus congestus setup is applied to compare droplet spectra grown by the diffusion of water vapor in Eulerian bin and particle-based Lagrangian microphysics schemes. Bin microphysics represent droplet spectral evolution applying the spectral density function. In the Lagrangian microphysics, computational particles referred to as superdroplets are followed in time and space with each superdroplet representing a multiplicity of natural cloud droplets. The same cloud condensation nuclei (CCN) activation and identical representation of the droplet diffusional growth allow the comparison. The piggybacking method is used with the two schemes operating in a single simulation, one scheme driving the dynamics and the other one piggybacking the simulated flow. Piggybacking allows point-by-point comparison of droplet spectra predicted by the two schemes. Overall, the results show the impact of inherent limitations of the two microphysics simulation methods, numerical diffusion in the Eulerian scheme and a limited number of superdroplets in the Lagrangian scheme. Numerical diffusion in the Eulerian scheme results in a more dilution of the cloud upper half and thus smaller cloud droplet mean radius. The Lagrangian scheme typically has larger spatial fluctuations of droplet spectral properties. A significantly larger mean spectral width in the bin microphysics across the entire cloud depth is the largest difference between the two schemes. A fourfold increase of the number of superdroplets per grid volume and a twofold increase of the spectral resolution and thus the number of bins have small impact on the results and provide only minor changes to the comparison between simulated cloud properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Absorption and scattering of light by nonspherical particles

Using the example of the polarization of scattered light, it is shown that the scattering matrices for identical, randomly ordered particles and for spherical particles are unequal. The spherical assumptions of Mie theory are therefore inconsistent with the random shapes and sizes of atmospheric particulates. The implications for corrections made to extinction measurements of forward scattering light are discussed. Several analytical methods are examined as potential bases for developing more accurate models, including Rayleigh theory, Fraunhoffer Diffraction theory, anomalous diffraction theory, Rayleigh-Gans theory, the separation of variables technique, the Purcell-Pennypacker method, the T-matrix method, and finite difference calculations.

Bohren, C. F.↗