Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “massively”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

DESI Massive Poststarburst Galaxies at z ~ 1.2 Have Compact Structures and Dense Cores

Poststarburst galaxies (PSBs) are young quiescent galaxies that have recently experienced a rapid decrease in star formation, allowing us to probe the fast-quenching period of galaxy evolution. In this work, we obtained Hubble Space Telescope (HST)/WFC3 F110W imaging to measure the sizes of 171 massive (log(M $\ast$ /M ⊙ ) ~ 11) spectroscopically identified PSBs at 1 < z 1.3 selected from the DESI Survey Validation luminous red galaxy sample. This statistical sample constitutes an order of magnitude increase from the ~20 PSBs with space-based imaging and deep spectroscopy. We perform structural fitting of the target galaxies with pysersic and compare them to quiescent and star-forming galaxies in the 3D-HST survey. We find that these PSBs are more compact than the general population of quiescent galaxies, lying systematically ~0.1 dex below the established size–mass relation. However, their central surface mass densities are similar to those of their quiescent counterparts (log(Σ 1kpc /(M ⊙ kpc -2 ))~10.1). These findings are easily reconciled by later ex situ growth via minor mergers or a slight progenitor bias. These PSBs are round in projection (b/a median ~ 0.8), suggesting that they are primarily spheroids, not disks, in 3D. We find no correlation between the time since quenching and light-weighted PSB sizes or central densities. This disfavors apparent structural growth due to the fading of centralized starbursts in this galaxy population. Instead, we posit that the fast quenching of massive galaxies at this epoch occurs preferentially in galaxies with preexisting compact structures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Strong Coupling of Hydrodynamics and Reactions in Nuclear Statistical Equilibrium for Modeling Convection in Massive Stars

We build on the simplified spectral deferred corrections (SDC) coupling of hydrodynamics and reactions to handle the case of nuclear statistical equilibrium (NSE) and electron/positron captures/decays in the cores of massive stars. Our approach blends a traditional reaction network on the grid with a tabulated NSE state from a very large, ${\mathcal O }(100)$ nuclei network. We demonstrate how to achieve second-order accuracy in the simplified-SDC framework when coupling NSE to hydrodynamics, with the ability to evolve the star on the hydrodynamics time step. We discuss the application of this method to convection in massive stars leading up to core collapse. We also show how to initialize the initial convective state from a 1D model in a self-consistent fashion. All of these developments are done in the publicly available Castro simulation code and the entire simulation methodology is fully GPU-accelerated.

Explosive nucleosynthesis↗
DOE OSTI2025

Impact of the Latest 22 Ne+ α Reaction Rates on Nucleosynthesis in Massive Stars and Galactic Chemical Evolution

In massive stars (initial mass of ≳9 M ⊙ ), the weak s- (slow neutron capture) process produces elements between Fe and Zr, enriching the Galaxy with these elements through core-collapse supernova explosions. The weak s-process nucleosynthesis is driven by neutrons produced in the 22 Ne(α, n) 25 Mg reaction during convective He core and C shell burning. The yields of heavy elements thus depend on the 22 Ne(α, n) 25 Mg and the competitive 22 Ne(α, γ) 26 Mg reaction rates, which are dominated by several narrow-resonance reactions. While the accuracy of these rates has been under debate for decades, recent experimental efforts, including ours, drastically reduced these uncertainties. In this work, we use a set of 280 massive star nucleosynthesis models calculated using different 22 Ne(α, n) 25 Mg and 22 Ne(α, γ) 26 Mg rates and a galactic chemical evolution (GCE) study to probe their impact on the weak s-process elemental abundances in the Galaxy. The GCE was computed with the OMEGA+ code, using the new sets of stellar yields with different 22 Ne+α rates. From GCE, we find that these rates are causing up to 0.45 dex of variations in the [Cu/Fe], [Ga/Fe], and [Ge/Fe] ratios predicted at solar metallicity. The greatest impact on the stellar nucleosynthesis and GCE results derives from uncertainties in the (α, n) strength (ωγ (α,n) ) of the E x = 11.32 MeV resonance. We show that variations observed in the GCE calculations for weak s-process elements become negligibly smaller than dispersions found in observations once the ωγ (α,n) is accurately determined within the uncertainty of 10%–20% (typically reported experimental errors for the resonance) in future nuclear physics experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

QRCODE: Massively parallelized real-time time-dependent density functional theory for periodic systems

We present a new software module, QRCODE (Quantum Research for Calculating Optically Driven Excitations), for massively parallelized real-time time-dependent density functional theory (RT-TDDFT) calculations of periodic systems in the open-source Qbox software package. Our approach utilizes a custom implementation of a fast Fourier transformation scheme that significantly reduces inter-node message passing interface (MPI) communication of the major computational kernel and shows impressive scaling up to 16,344 CPU cores. In addition to improving computational performance, QRCODE contains a suite of various time propagators for accurate RT-TDDFT calculations. As benchmark applications of QRCODE, we calculate the current density and optical absorption spectra of hexagonal boron nitride (h-BN) and photo-driven reaction dynamics of the ozone-oxygen reaction. We also calculate the second and higher harmonic generation of monolayer and multi-layer boron nitride structures as examples of large material systems. Our optimized implementation of RT-TDDFT in QRCODE enables large-scale calculations of real-time electron dynamics of chemical and material systems with enhanced computational performance and impressive scaling across several thousand CPU cores.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

OpenSn: A massively parallel, open-source simulation environment for discrete ordinates radiation transport

OpenSn is an open-source, massively parallel deterministic radiation transport code for solving the discrete-ordinates ( S N ) form of the Boltzmann transport equation on unstructured, arbitrary polyhedral meshes. It supports high-fidelity simulations involving steady-state, eigenvalue, and adjoint problems for neutral particles (e.g., neutrons, photons, multi-particles), using the multigroup approximation in energy. OpenSn combines angular discretization via discrete ordinates with a discontinuous Galerkin finite element method (DGFEM) in space, enabling accurate resolution of transport physics on arbitrary polyhedral cells, included locally refined spatial grids. It includes multiple angular quadrature types, including locally refined angular quadratures. Written in modern C++ with a Python API, OpenSn runs efficiently on platforms ranging from laptops to supercomputers. The transport sweep algorithm is implemented using a task-based, directed-acyclic-graph (DAG) approach for each angle and supports asynchronous parallelism across thousands of MPI ranks. Group-set aggregation improves compute intensity, and synthetic acceleration techniques (e.g., diffusion synthetic acceleration, second-moment method) enhance solver convergence. OpenSn has been verified on reactor physics problems and demonstrated excellent weak and strong scaling performance on more than 32,768 processes, making it a versatile and robust platform for large-scale transport simulations in complex geometries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Radiation asymmetry during shattered pellet and massive gas injection in DIII-D

Infrared thermography of the first wall in DIII-D is used to show the toroidal peaking of plasma radiation during mitigated disruptions with shattered pellet injection (SPI) and with massive gas injection (MGI). During MGI, the radiation peak location is shown to be due to the n = 1 magnetohydrodynamic (MHD) activity and continuously variable in toroidal phase based on externally applied error fields, consistent with previous experimental results. Furthermore, the measured toroidal peaking factor (TPF) is in agreement with the previous estimates based on radiometry, with a total TPF of 1.07 ± 0.05 when integrated over the entire duration of the disruption, and a value of 1.3 ± 0.1 during the thermal quench (TQ). For SPI, the location of the radiation peak is found to be determined predominantly by the particle source location, and thus can be varied on DIII-D by injecting from either of two toroidally separated injectors. 3D non-linear MHD simulations with the NIMROD code support this picture. The localized SPI particle source leads to higher peaking than for MGI, with a total TPF of 1.3 ± 0.1 over the full disruption and a TQ TPF of 1.9 +0.5/−0.3, a value also consistent with the NIMROD simulations. This TQ value is close to or potentially exceeding the allowable limit for the previously planned Be first wall in ITER, pointing to the importance of radiation asymmetries depending on the choice of first wall material.

Physics - Plasma physics↗
DOE OSTI2024

Extended-MHD simulations of disruption mitigation via massive gas injection in SPARC

Recent developments to the M3D-C1 code enable higher fidelity modeling of disruptions, and can be applied in the design verification of reactor-scale tokamaks. Among these new capabilities is a method to mesh conducting vessel structures such as coils and passive plates, packing of the toroidal mesh around gas injectors, as well as anisotropic resistivity inside the vessel structures. We present extended-magnetohydrodynamic (MHD) simulations of disruption mitigation via massive gas injection (MGI) in SPARC. The goal of this study is to inform the disruption mitigation layout of SPARC and aid in the design of an effective gas injector configuration. Fully three-dimensional simulations with M3D-C1 are carried out for various injector configurations with the primary goal of determining the effect of different MGI parameters on heat loads and vessel forces. The simulations include a model for impurity ionization, recombination, advection and radiation, as well as spatially resolved conducting structures around the plasma. A localized mixture of deuterium and neon with a small toroidal and poloidal width is injected in up to six locations. We demonstrate that M3D-C1 can model a rapid shutdown via MGI using narrow and more realistic gas plumes than in previous simulations. As a result of the q = 1 surface in the SPARC baseline case a sawtooth is observed early in the simulations. Despite the sawtooth and the onset of edge MHD instabilities, the impurity distribution remains localized around the injector locations, but enables a radiative shutdown of the plasma. We find that using the maximum of six gas injectors results in a lower peaking factor and leads to a more even distribution of radiation toroidally than using two injectors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Predicting synthetic mRNA stability using massively parallel kinetic measurements, biophysical modeling, and machine learning

Abstract mRNA degradation is a central process that affects all gene expression levels, though it remains challenging to predict the stability of a mRNA from its sequence, due to the many coupled interactions that control degradation rate. Here, we carried out massively parallel kinetic decay measurements on over 50,000 bacterial mRNAs, using a learn-by-design approach to develop and validate a predictive sequence-to-function model of mRNA stability. mRNAs were designed to systematically vary translation rates, secondary structures, sequence compositions, G-quadruplexes, i-motifs, and RppH activity, resulting in mRNA half-lives from about 20 seconds to 20 minutes. We combined biophysical models and machine learning to develop steady-state and kinetic decay models of mRNA stability with high accuracy and generalizability, utilizing transcription rate models to identify mRNA isoforms and translation rate models to calculate ribosome protection. Overall, the developed model quantifies the key interactions that collectively control mRNA stability in bacterial operons and predicts how changing mRNA sequence alters mRNA stability, which is important when studying and engineering bacterial genetic systems.

Cetnar, Daniel P.↗
DOE OSTI2025

Massively parallel reporter assays and mouse transgenic assays provide correlated and complementary information about neuronal enhancer activity

High-throughput massively parallel reporter assays (MPRAs) and phenotype-rich in vivo transgenic mouse assays are two potentially complementary ways to study the impact of noncoding variants associated with psychiatric diseases. Here, we investigate the utility of combining these assays. Specifically, we carry out an MPRA in induced human neurons on over 50,000 sequences derived from fetal neuronal ATAC-seq datasets and enhancers validated in mouse assays. We also test the impact of over 20,000 variants, including synthetic mutations and 167 common variants associated with psychiatric disorders. We find a strong and specific correlation between MPRA and mouse neuronal enhancer activity. Four out of five tested variants with significant MPRA effects affected neuronal enhancer activity in mouse embryos. Mouse assays also reveal pleiotropic variant effects that could not be observed in MPRA. Our work provides a catalog of functional neuronal enhancers and variant effects and highlights the effectiveness of combining MPRAs and mouse transgenic assays.

Kosicki, Michael↗
DOE OSTI2025

Exploring quantum statistics for massive Dirac and Majorana neutrinos using spinor-helicity techniques

Recently, there has been interest in the applicability of quantum statistics to distinguish Dirac from Majorana neutrinos in multineutrino final states. In particular, debate has arisen over the validity of the Dirac-Majorana confusion theorem in these processes, i.e., that any distinction between the Dirac and Majorana processes goes to zero as the neutrino mass goes to zero. Here we approach this problem equipped with spinor-helicity methods generalized for massive Dirac and Majorana fermions. We explicitly calculate all helicity amplitudes, and their squares, for the decay of a light scalar particle to two neutrinos and two oppositely charged leptons. This allows us to pinpoint the crucial steps which could lead to claims of a violation of the confusion theorem. We show that, if the correct antisymmetrization of Dirac to Majorana amplitudes is used, identification of which is clear in this framework, and all relevant contributions are appropriately summed, a scalar decay into two charged leptons and two neutrinos satisfies the Dirac-Majorana confusion theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

First displaced vertex search for electroproduced dark-sector strongly interacting massive particles by the HPS experiment

The Heavy Photon Search experiment (HPS) is a fixed-target, electron beam experiment designed to search for 𝑒 + ⁢𝑒 − mass resonances and displaced decays using a forward acceptance spectrometer. This paper details the search for naturally long-lived “dark” vector mesons (𝑉 𝐷 ) arising from a dark sector of beyond-Standard-Model strongly interacting massive particles, characterized by a QCD-like 𝑆⁢𝑈⁢(3) 𝐷 symmetry and coupled to the Standard Model photon via a new 𝑈⁢(1) 𝐷 gauge interaction mediated by the “heavy photon,” or 𝐴′. The results are based on an integrated luminosity of 10 608 nb −1 collected during the 2016 HPS engineering run. The displaced vertex search for 𝑉 𝐷 → 𝑒 +⁢ 𝑒 − in the 𝑒 +⁢ 𝑒 − invariant mass range 39–179 MeV showed no statistically significant evidence for signal above the QED background.

extensions of gauge sector↗
DOE OSTI2026

Massive tree-level splitting functions beyond kinematical limits

We present a compact form of the massive 1 → 3 tree-level QCD splitting functions and discuss a decomposition of the results in terms of lower-order expressions, scalar dipole antenna functions and pure higher-order remainders. The two-gluon radiator functions introduced in this context are novel and generalize expressions obtained from the double-soft approximation. Our results are obtained without reference to soft or quasicollinear limits.

Höche, Stefan [Fermi National Accelerator Laborato↗
DOE OSTI2025

Emergent Nonthermal Fluid from Jets in the Massive Schwinger Model Using Tensor Networks

We analyze the correlation between the energy, momentum, and spatial entanglement produced by two luminal jets in the massive Schwinger model. Using tensor network methods, we show that for m/g > 1/π, in the vicinity of the strong- to weak-coupling transition, a nearly perfect and chargeless effective fluid behavior appears around the midrapidity region with a universal energy-pressure relationship. The evolution of energy and pressure is strongly correlated with the rise of the spatial entanglement entropy, indicating a key role of quantum dynamics. Some of these observations may be used to analyze high multiplicity jet fragmentation events, energy-energy and energy-charge correlators at current collider energies.

Ab initio calculations↗
DOE OSTI2024

Astrophysical constraints from synchrotron emission on very massive decaying dark matter

If the cosmological dark matter (DM) couples to Standard Model (SM) fields, it can decay promptly to SM states in a highly energetic hard process, which subsequently showers and hadronizes to give stable particles including e ± , γ , p ± , and ν ν ¯ at lower energy. If the DM particle is very heavy, the high-energy e ± , due to the Klein-Nishina cross section suppression, preferentially lose energy via synchrotron emission which, in turn, can be of unusually high energies. Here, we present previously unexplored bounds on heavy decaying DM up to the Planck scale, by studying the synchrotron emission from the e ± produced in the ambient Galactic magnetic field. In particular, we explore the sensitivity of the resulting constraints on the DM decay width to (i) different SM decay channels, to (ii) the Galactic magnetic field configurations, and to (iii) various different DM density profiles proposed in the literature. We find that constraints from the synchrotron component complement and improve on constraints from very high-energy cosmic-ray and gamma-ray observatories targeting the prompt emission when the DM is sufficiently massive, most significantly for masses in excess of 10 12 GeV . Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Quasiparton distributions in massive QED2: Toward quantum computation

We analyze the quasiparton distributions of the lightest 𝜂′ meson in massive two-dimensional quantum electrodynamics (QED2) by exact diagonalization. The Hamiltonian and boost operators are mapped onto spin qubits in a spatial lattice with open boundary conditions. The lowest excited state in the exact diagonalization is shown to interpolate continuously between an anomalous 𝜂′ state at strong coupling, and a nonanomalous heavy meson at weak coupling, with a cusp at the critical point. The boosted 𝜂′ state follows relativistic kinematics but with large deviations in the luminal limit. The spatial quasiparton distribution function and amplitude for the 𝜂′ state are computed numerically for increasing rapidity both at strong and weak coupling, and compared to the exact light front results. The numerical results from the boosted form of the spatial parton distributions, compare fairly with the inverse Fourier transformation of the luminal parton distributions, derived in the lowest Fock space approximation. Our analysis points out some of the limitations facing the current lattice program for the parton distributions.

Lattice field theory↗
DOE OSTI2024

SPT clusters with DES and HST weak lensing. II. Cosmological constraints from the abundance of massive halos

We present cosmological constraints from the abundance of galaxy clusters selected via the thermal Sunyaev-Zel’dovich (SZ) effect in South Pole Telescope (SPT) data with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). The cluster sample is constructed from the combined SPT-SZ, SPTpol ECS, and SPTpol 500d surveys, and comprises 1,005 confirmed clusters in the redshift range 0.25–1.78 over a total sky area of 5200 deg 2 . We use DES Year 3 weak-lensing data for 688 clusters with redshifts 𝑧 < 0.95 and HST weak-lensing data for 39 clusters with 0.6 < 𝑧 < 1.7. The weak-lensing measurements enable robust mass measurements of sample clusters and allow us to empirically constrain the SZ observable-mass relation without having to make strong assumptions about, e.g., the hydrodynamical state of the clusters. For a flat Λ⁢ CDM cosmology, and marginalizing over the sum of massive neutrinos, we measure Ω m = 0.286 ± 0.032, 𝜎 8 = 0.817 ± 0.026, and the parameter combination 𝜎 8 ⁢(Ω m /0.3) 0.25 = 0.805 ± 0.016. Our measurement of 𝑆 8 ≡ 𝜎 8 ⁢$\sqrt{Ω_{m}/0.3}$ = 0.795 ± 0.029 and the constraint from Planck CMB anisotropies (2018 TT, TE, EE+lowE) differ by 1.1⁢𝜎. In combination with that Planck dataset, we place a 95% upper limit on the sum of neutrino masses ∑𝑚 𝜈 < 0.18 eV. When additionally allowing the dark energy equation of state parameter 𝑤 to vary, we obtain 𝑤 = −1.45 ± 0.31 from our cluster-based analysis. In combination with Planck data, we measure 𝑤 =−1.3⁢4$^{+0.22}_{−0.15}$, or a 2.2⁢𝜎 difference with a cosmological constant. We use the cluster abundance to measure 𝜎8 in five redshift bins between 0.25 and 1.8, and we find the results to be consistent with structure growth as predicted by the Λ⁢ CDM model fit to Planck primary CMB data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Quantum real-time evolution of entanglement and hadronization in jet production: Lessons from the massive Schwinger model

The possible link between entanglement and thermalization, and the dynamics of hadronization are addressed by studying the real-time response of the massive Schwinger model coupled to external sources. This setup mimics the production and fragmentation of quark jets, as the Schwinger model and quantum chromodynamics (QCD) share the properties of confinement and chiral symmetry breaking. By using simulations of quantum dynamics on classical hardware, we study the entanglement between the produced jets, and observe the growth of the corresponding entanglement entropy in time. This growth arises from the increased number of contributing eigenstates of the reduced density matrix with sufficiently large and close eigenvalues. We also investigate the physical nature of these eigenstates, and find that at early times they correspond to fermionic Fock states. We then observe the transition from these fermionic Fock states to mesonlike bound states as a function of time. In other words, we observe how hadronization develops in real time. At late times, the local observables at midrapidity (such as the fermion density and the electric field) approach approximately constant values, suggesting the onset of equilibrium and approach to thermalization. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗