Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite size”

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 235 records · Page 13

Requirements, Status, and Plans for Track Reconstruction at the sPHENIX Experiment

sPHENIX is a new experiment that is being constructed at the Relativistic Heavy Ion Collider at Brookhaven National Laboratory. The primary physics goals of sPHENIX are to measure jets, their substructure, and the upsilon resonances in p+p, p+Au and Au+Au collisions. To realize these goals, a tracking system composed of a time projection chamber and several silicon detectors will be used to identify tracks that correspond to jets and upsilon decays. However, the sPHENIX experiment will collect approximately 200 PB of data utilizing a finite-sized computing center; thus, performing track reconstruction in a timely manner is a challenge due to the large occupancy of heavy-ion collisions. The sPHENIX experiment, its track reconstruction, and the need for implementing faster track-fitting algorithms, such as that provided by the A Common Tracking Software package, into the sPHENIX software stack are discussed.

Osborn, Joe↗

A Flexible Power Method for Solving Infinite Dimensional Tensor Eigenvalue Problems

We propose a flexible power method for computing the leftmost, i.e., algebraically smallest, eigenvalue of an infinite dimensional tensor eigenvalue problem, $H x = \lambda x$, where the infinite dimensional symmetric matrix $H$ exhibits a translational invariant structure. We assume the smallest eigenvalue of $H$ is simple and apply a power iteration of $e^{-H}$ with the eigenvector represented in a compact way as a translational invariant infinite Tensor Ring (iTR). Hence, the infinite dimensional eigenvector can be represented by a finite number of iTR cores of finite rank. In order to implement this power iteration, we use a small parameter $t$ so that the infinite matrix-vector operation $e^{-Ht}x$ can efficiently be approximated by the Lie product formula, also known as Suzuki--Trotter splitting, and we employ a low rank approximation through a truncated singular value decomposition on the iTR cores in order to keep the cost of subsequent power iterations bounded. We also use an efficient way for computing the iTR Rayleigh quotient and introduce a finite size iTR residual which is used to monitor the convergence of the Rayleigh quotient and to modify the timestep $t$. In this paper, we discuss 2 different implementations of the flexible power algorithm and illustrate the automatic timestep adaption approach for several numerical examples.

Beeumen, Roel Van↗

Constrained curve fitting for semi-parametric models with radial basis function networks

Common to many analysis pipelines in lattice gauge theory and the broader scientific discipline is the need to fit a semi-parametric model to data. We propose a fit method that utilizes a radial basis function network to approximate the non-parametric component of such models. The approximate parametric model is fit to data using the basin hopping global optimization algorithm. Parameter constraints are enforced through Gaussian priors. The viability of our method is tested by examining its use in a finite-size scaling analysis of the $q$-state Potts model and $p$-state clock model with $q=2,3$ and $p=4,\infty$.

Peterson, Curtis T.↗

Bridging the time scale in exascale computing of chemical systems (Final Technical Report)

This report summarizes the work carried out with support of the United States Department of Energy under Award DE-SC0019441. The theme of this project was to develop and apply methods that allowed for the acceleration of atomistic calculations, particularly in challenging areas such as multiphase systems, electrified interfaces, uncertainty estimation, and applications requiring chemical accuracy, which tend to be applications where simulation time is severely bottlenecked by the computational time requirements. Much of the focus was on the application of emerging machine-learning methodologies, although a wide range of methodologies were employed. This report has two major sections. The first focuses on the methodological advances themselves. Within this part, we report a number of major advances, a few examples of which are described here. We report the first machine-learning scheme for the acceleration of electronically grand-canonical calculations (that is, those applicable to electrochemistry). We report new methods of performing transfer learning, in which physics-based priors can be used to provide predictions, often with uncertainty estimates, of images well outside of training sets; we also offer ways to fine-tune these transfer-learning models. We provide a new systematic means to generate and apply minimal training data sets to very large (10,000’s of atoms) systems, with only small training sets appropriate for electronic structure. We developed new methodologies to integrate surface vibrations into surface adsorption calculations. We made advances to the applicability of diffusion Monte Carlo methods to allow (learned) force prediction, finite-size error correction, and force-free means of searching for transition states. We integrated machine-learned atomistic predictions into mechanism generation codes. Additionally, we released new software including AmpTorch, a modernized version of our original atomistic machine-learning code Amp. The second part of this report focuses on the scientific applications that accompanied, and were often enabled by, the methodological advances described earlier. A few examples follow, but full details are in the individual chapters of the report. For example, we developed a general theory of phonon-induced friction on molecular adsorbates. We showed fundamentally how solvent influences the adsorption and desorption process and how it differs from the processes typically involved at the solid–gas interface, making aqueous-phase and electrocatalysis different from traditional thermocatalysis. We examined how metal–insulator and magnetic transitions can be probed, and accelerated exciton dynamics via Frenkel Hamiltonian parameters. We showed that the nearsighted force-training approach, developed within this project, can predict both the stability and reactivity of large nanoparticles, and can also lead to insights on catalyst coverage on binding energies and entropies. These applied studies, which generally integrated with our method development, allowed us to push forward the theoretical understanding of several reaction classes.

08 HYDROGEN↗

Investigating SU(3) with Nf=8 fundamental fermions at strong renormalized coupling

Lattice simulations have observed a novel strong coupling symmetric mass generation (SMG) phase for the SU(3) gauge system with $N_f=8$ fundamental fermions (represented by two sets of staggered fields) at very large renormalized coupling ($g^2_{GF} \gtrsim 25$). The results of Phys.Rev.D 106 (2022) 014513 suggest that the SMG phase is separated from the weak coupling, conformal phase by a continuous phase transition, implying that the SMG phase exists in the continuum limit. To scrutinize these findings, we are generating a set of large volume zero temperature ensembles using nHYP improved staggered fermions with additional Pauli-Villars fields to tame gauge field fluctuations. We consider the low-lying meson spectrum and verify the existence of the SMG phase. Based on a finite size scaling analysis we predict that the phase transition between the strong and weak coupling phases is likely governed by a merged fixed point that is ultraviolet in the strong coupling but infrared in the weak coupling side. This finding suggests that the SU(3) 8-flavor system sits at the opening of the conformal window

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Jefferson Lab High Precision Proton Radius Measurement Experiment: PRad

Proton is not a point-like particle; it has a finite size and an internal structure. The proton charge radius is usually measured through two different methods: using spectroscopy of hydrogen atoms, or through electron-proton (e - p) elastic scattering at low momentum transfer. In 2010, Phol et al [1] published a new measurement for proton charge radius obtained using the spectroscopy method from muonic hydrogen atoms. An ordinary hydrogen (H) atom is a bound state of an electron orbiting a proton. Likewise, a muonic hydrogen atom has a muon orbiting a proton. Since a muon has 200 times more mass than an electron, resulting in a much smaller orbiting radius, its orbit is much more sensitive to the proton charge distribution in space than the electron’s of an ordinary H-atom. Thus, the measurement from muonic hydrogen provides a 10 times more precise result than the previous methods. The radius from the muonic hydrogen measurement reported in [1] was significantly smaller than all the p

Bai, Xinzhan↗

Architectural scaling tradeoffs in modular 3D bosonic quantum processors

We propose a modular three-dimensional bosonic quantum processor built from repeatable coupled-cavity modules linked by configurable interconnect networks. Using hardware-motivated graph-theoretic measures, we compare nearest-neighbor, hub-based, and hybrid architectures in terms of interconnect count, communication distance, resource concentration, and implementation complexity. Rather than identifying a universally optimal topology, our analysis shows how these architectures redistribute the costs of scaling, including wiring and port requirements, nonlocal communication distance, exposure to shared resources, routing bottlenecks, and scheduling overhead. Case studies of a \(3\times3\) processor and a larger hierarchical architecture further distinguish finite-size performance from asymptotic scaling. The resulting framework provides a systematic basis for evaluating modular three-dimensional bosonic processors and for identifying the device-level parameters required for quantitative hardware design.

Zhu, Shaojiang [Fermilab] (ORCID:0000000293180092)↗

Characterizing the impact of finite matrix block size on conservative particle transport through three-dimensional fracture networks

Mass transfer of solutes between fractures and the surrounding rock matrix exerts a noticeable signature on the tail of travel time distributions. When the width of the matrix is assumed to be infinite and advective transport through the fracture is sufficiently fast, the tails of the travel time distributions exhibit a classically expected slope of ψ(t) ∝ t -3/2 . However, studies have yet to characterize how solute transfer between fractures via diffusion through finite matrix blocks influences the tail’s slope in three-dimensional fractured media. Here, in this study, we assess the impact of finite matrix block size on breakthrough curve shape at different spatio-temporal scales by con ducting particle tracking simulations in three-dimensional discrete fracture networks. We consider a variety of hydrodynamic and geostructural proper ties to determine their relative impact on the resulting travel time distributions. We observe that the impact of matrix diffusion through a finite block on travel time distributions is similar to that of an infinite matrix block when the fracture spacing is sufficiently large, matrix diffusion is relatively weak, or transport is considered at an early control plane distance. We observe that the converse of these conditions, results in deviations from the classical ψ(t) ∝ t -3/2 scaling. These results provide a first step toward developing a metric to assess when finite block size effects are expected to significantly influence transport.

58 GEOSCIENCES↗

Role of boundary conditions in quantum computations of scattering observables

Quantum computing may offer the opportunity to simulate strongly interacting field theories, such as quantum chromodynamics, with physical time evolution. This would give access to Minkowski-signature correlators, in contrast to the Euclidean calculations routinely performed at present. However, as with present-day calculations, quantum computation strategies still require the restriction to a finite system size, including a finite, usually periodic, spatial volume. In this work, we investigate the consequences of this in the extraction of hadronic and Compton-like scattering amplitudes. Using the framework presented in Briceño et al., we estimate the volume effects for various 1+1D Minkowski-signature quantities and show that these can be a significant source of systematic uncertainty, even for volumes that are very large by the standards of present-day Euclidean calculations. We then present an improvement strategy, based in the fact that the finite volume has a reduced symmetry. This implies that kinematic points, which yield the same Lorentz invariants, may still be physically distinct in the periodic system. As we demonstrate, both numerically and analytically, averaging over such sets can significantly suppress the unwanted volume distortions and improve the extraction of the physical scattering amplitudes. As the improvement strategy is based only in kinematics, it can be applied without detailed knowledge of the system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

How eigenmode self-interaction affects zonal flows and convergence of tokamak core turbulence with toroidal system size

Self-interaction is the process by which a microinstability eigenmode that is extended along the direction parallel to the magnetic field interacts non-linearly with itself. This effect is particularly significant in gyrokinetic simulations accounting for kinetic passing electron dynamics and is known to generate stationary $E\times B$ zonal flow shear layers at radial locations near low-order mode rational surfaces (Weikl et al. Phys. Plasmas, vol. 25, 2018, 072305). Here, we find that self-interaction, in fact, plays a very significant role in also generating fluctuating zonal flows, which is critical to regulating turbulent transport throughout the radial extent. Unlike the usual picture of zonal flow drive in which microinstability eigenmodes coherently amplify the flow via modulational instabilities, the self-interaction drive of zonal flows from these eigenmodes are uncorrelated with each other. It is shown that the associated shearing rate of the fluctuating zonal flows therefore reduces as more toroidal modes are resolved in the simulation. In simulations accounting for the full toroidal domain, such an increase in the density of toroidal modes corresponds to an increase in the toroidal system size, leading to a finite system size effect that is distinct from the well-known profile shearing effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transducer Resolution Effect on Pressure Fluctuations Beneath Hypersonic Turbulent Boundary Layers

The size of a pressure transducer is known to affect the accuracy of measurements of wall-pressure fluctuations beneath a turbulent boundary layer because of spatial averaging over the sensing area of the transducer. In this paper, the effect of finite transducer size is investigated by applying spatial averaging or wavenumber filters to a database of hypersonic wall pressure generated from a direct numerical simulation (DNS) that simulates the turbulent portion of the boundary layer over a sharp 7° half-angle cone at nominally Mach 8. Here, a good comparison between the DNS and the experiment in the Sandia Hypersonic Wind Tunnel at Mach 8 is achieved after spatial averaging is applied to the DNS data over an area similar to the sensing area of the transducer. The study shows that a finite sensor size similar to that of the PCB132 transducer can cause significant attenuation in the root-mean-square and power spectral density (PSD) of wall-pressure fluctuations, and the attenuation effect is identical between cone and flat plate configurations at the same friction Reynolds number. The Corcos theory is found to successfully compensate for the attenuated high-frequency components of the wall-pressure PSD.

42 ENGINEERING↗

Simulation and Modeling of Optical Properties of U, Th, Pb, and Co Nanoparticles of Interest to Nuclear Security Using Finite Element Analysis

In this work, the optical characteristics of uranium (U), lead (Pb), cobalt (Co), and thorium (Th) nanoparticles are fashioned and simulated employing the finite element analysis (FEA) approach concerning multiple particle sizes. Applying finite element analysis, it was found that the simulated absorption peaks of electronic excitations of nuclear nanoparticles are red-shifted from 365 nm to 555 nm for U; from 355 nm to 550 nm for Pb; from 415 nm to 610 nm for Co; and from 350 nm to 540 nm for Th, comparing expanding particle sizes from 60 nm to 100 nm (except for Co, which varied from 70 nm to 100 nm). The FEA-simulated optical band gap energies and far-field radiation patterns were also obtained for nuclear materials. The simulation approach in this research enables the prediction of optical properties and design of nuclear materials before manufacture for nuclear security applications.

36 MATERIALS SCIENCE↗

Electronic structure of intrinsic defects in c -gallium nitride: Density functional theory study without the jellium approximation

Here, we report the first nonjellium, systematic, density functional theory (DFT) study of intrinsic and extrinsic defects and defect levels in zinc-blende (cubic) gallium nitride. We use the local moment counter charge (LMCC) method, the standard Perdew-Becke-Ernzerhoff (PBE) exchange-correlation potential, and two pseudopotentials, where the Ga 3$\textit{d}$ orbitals are either in the core ($d^0$) or explicitly in the valence set ($d^{10}$). We studied 64, 216, 512, and 1000 atom supercells, and demonstrated convergence to the infinite limit, crucial for delineating deep from shallow states near band edges, and for demonstrating the elimination of finite cell-size errors. Contrary to common claims, we find that exact exchange is not required to obtain defect levels across the experimental band gap. As was true in silicon, silicon carbide, and gallium arsenide, the extremal LMCC defect levels of the aggregate of defects yield an effective LMCC defect band gap that is within 10% of the experimental gap (3.3 eV) for both pseudopotentials. We demonstrate that the gallium vacancy is more complicated than previously reported. There is dramatic metastability–a nearest-neighbor nitrogen atom shifts into the gallium site, forming an antisite, nitrogen vacancy pair, which is more stable than the simple vacancy for positive charge states. Our assessment of the $d^0$ and $d^{10}$ pseudopotentials yields minimal differences in defect structures and defect levels. The better agreement of the $d^0$ lattice constant with experiment suggests that the more computationally economical $d^0$ pseudopotentials are sufficient to achieve the fidelity possible within the physical accuracy of DFT, and thereby enable calculations in larger supercells necessary to demonstrate convergence with respect to finite size supercell errors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Analysis of Dynamical Peculiarities in Nanoalloys at Subsystems Level: Dynamical Degrees of Freedom, Temperature Differences, and the Chameleon Effect

Abstract A novel analysis of the dynamical behavior of nanoalloy systems, as represented by model Ni/Al 13‐atom clusters, over a broad range of energies that cover the stage‐wise transition of the systems from their solid‐like to liquid‐like state is presented. Conceptually, the analysis is rooted in partitioning the systems into judiciously chosen subsystems and characterizing the latter in terms of subsystem‐specific dynamical descriptors that include dynamical degrees of freedom , root‐mean‐square bond‐length fluctuation , and element‐specific subsystem temperature . The analysis reveals a host of intriguing new peculiarities in the dynamical behavior of the Ni/Al 13‐mers, among which are what we call the chameleon effect and the difference in the temperatures of the Ni and Al subsystems at high energies, a difference that strongly depends on the cluster composition and also changes with energy. These do not have an analog in pure Ni 13 and Al 13 and are explained in terms of the coupled effects of the difference between the masses of the Ni and Al atoms (the mass effect) and of the difference in the anharmonicity of the overall interaction potential as experienced by the Ni and Al subsystems of the clusters (the potential effect).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of sample size on the maximum value distribution of fatigue driving forces in metals and alloys

An analytical framework is presented to predict the effects of sample size on the maximum value distribution (MVD) of the driving forces for fatigue crack formation in metals and alloys. The distribution of the maximum driving force for fatigue crack formation over the domain follows the generalized extreme value theory in the limit as the domain size increases to infinity. Here, a simulation-based analysis of microstructure influences on fatigue resistance for polycrystalline metals and alloys is very costly, and reaching those limits is intractable. This work models the MVD of Fatigue Indicator Parameters (FIPs), which serve as surrogate measures for the driving force for fatigue crack formation, at finite sample sizes prior to their convergence to a limiting extreme value distribution. Large-scale crystal plasticity finite element (CPFE) simulations of FCC Al 7075-T6 with microstructure realizations of various sizes are incorporated to calibrate and evaluate the developed framework, and a total of ∼6.5 million grains of Al 7075-T6 are examined. The calibrated analytical solution agrees well with the brute force Monte Carlo simulation results extracted from the CPFE simulations. Furthermore, the developed formulation can predict the MVD of FIPs for different sample sizes using a size-dependent parameter, and it is capable of accurately extrapolating the MVD of FIPs for much larger microstructure sample sizes than the size used for its calibration.

Crystal plasticity↗

Many-body quantum dynamics slows down at low density

We study quantum many-body systems with a global U(1) conservation law, focusing on a theory of N interacting fermions with charge conservation, or N interacting spins with one conserved component of total spin. We define an effective operator size at finite chemical potential through suitably regularized out-of-time-ordered correlation functions. The growth rate of this density-dependent operator size vanishes algebraically with charge density; hence we obtain new bounds on Lyapunov exponents and butterfly velocities in charged systems at a given density, which are parametrically stronger than any Lieb-Robinson bound. We argue that the density dependence of our bound on the Lyapunov exponent is saturated in the charged Sachdev-Ye-Kitaev model. We also study random automaton quantum circuits and Brownian Sachdev-Ye-Kitaev models, each of which exhibit a different density dependence for the Lyapunov exponent, and explain the discrepancy. We propose that our results are a cartoon for understanding Planckian-limited energy-conserving dynamics at finite temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Generalized kinetic equation for tokamak plasma equilibrium distribution function

A generalised kinetic equation for the equilibrium distribution function in a finite beta, arbitrary tokamak plasma is derived. The equation is correct to second order in ρ/L (ρ is the particle Larmor radius and L is the system size). Resolving finite Larmor radius length scales with no restriction on the ratio of poloidal to total equilibrium magnetic field, B $\vartheta$ /B, it generalises the drift kinetic theory of [Hazeltine Phys. Plasmas 15 (1973) 77] to the limit of B $\vartheta$ /B ~ 1 (e.g. to ensure validity for spherical tokamaks). Two cases are considered. The first provides the equilibrium distribution function, consistent with the generalised gyrokinetic formalism of [Dudkovskaia et al. Plasma Phys. Control. Fusion 65 (2023) 045010], derived specifically to capture neoclassical equilibrium currents in gyrokinetic stability analyses in strong gradient regions. The second assumes short length scales in the direction perpendicular to the magnetic field, which can occur as a result of small coherent magnetic structures in the plasma, such as neoclassical tearing mode magnetic islands close to threshold. This then extends the drift island equations of [Dudkovskaia et al. Nucl. Fusion 63 (2023) 016020] for the plasma response to magnetic islands to a spherical tokamak plasma configuration. In conclusion, resolving ρ ~ ρ $\vartheta$ (or B $\vartheta$ ~ B), where ρ $\vartheta$ is the particle poloidal Larmor radius, is also expected to influence calculations of the magnetic island propagation frequency and the associated contributions to the island onset conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗