Search NASA⌕ Search

SEARCH · Search NASA

Results for “strong scaling”

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

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

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Multiscale ACI Satellite Database

The SATELLITE_EAGLES_PNNL NetCDF dataset contains a suite of satellite- and reanalysis-derived atmospheric and surface parameters on a regular latitude–longitude grid. The dataset includes core geophysical fields such as land fraction, aerosol optical depth at multiple wavelengths (465, 550, 667, and 865 nm), sea surface temperature, estimated inversion strength, and various thermodynamic and dynamic quantities (e.g., relative humidity, vertical velocity, boundary-layer height, and surface fluxes) from both MERRA and ERA reanalysis products, provided as daily-mean and instantaneous values. A major component of the dataset consists of MODIS-retrieved cloud microphysical properties, including cloud droplet number concentration, cloud effective radius, optical thickness, and liquid water path, provided for three compositing regimes (“All,” “Q06,” and “G18”). Corresponding cloud-top parameters—temperature, height, and pressure—along with total and domain-mean cloud fraction fields are also included. The file further integrates additional satellite data from AMSR-E (for cloud water, rain water, and surface precipitation retrievals) and CERES (for top-of-atmosphere radiative fluxes, cloud fractions, and albedo). This dataset is designed to evaluate aerosol–cloud interactions in warm clouds, emphasizing the use of MODIS for deriving cloud droplet number concentration and liquid water path statistics. The complementary satellite and reanalysis fields are co-located and time-matched to the same instantaneous MODIS observations, enabling consistent comparisons between cloud properties, aerosol loading, and large-scale meteorological conditions. The dataset is recently featured in Christensen et al. (2025), Machine Learning Reveals Strong Grid-Scale Dependence in the Satellite Nd–LWP Relationship, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2025-3850, 2025.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Disentangling the Impacts of Microtopography and Shrub Distribution on Snow Depth in a Subarctic Watershed: Toward a Predictive Understanding of Snow Spatial Variability

Snow plays a critical role in carbon cycling, vegetation dynamics, and permafrost hydrology at high latitudes by influencing surface energy exchange. Predicting snow distribution patterns is essential for understanding the evolution of Arctic ecosystems, yet scaling process-level knowledge to landscape predictions remains challenging. Here, we analyze snow depth (2019 and 2022), terrain elevation, and vegetation height from a watershed on the Seward Peninsula, Alaska, to examine how topography and shrubs shape snow redistribution across spatial scales. We find that snow depth is strongly coupled to terrain at scales below ∼60 m but becomes increasingly decoupled at larger scales. The topographic model of snow depth variation, which transforms terrain data to align with these scale-dependent snow patterns, is well correlated with local snow depth variations (linear fit R 2 > 0.5 for 85% of 100-m patches). A machine learning reconstruction of shrub canopy snow trapping reveals a simple exponential relationship between canopy structure and snow accumulation ( R 2 = 0.59), highlighting the combined influence of topography and vegetation on snow distribution. Together, these empirical relationships capture much of the observed snow variability in the watershed ( R 2 = 0.49, root mean square error (RMSE) = 30 cm), though systematic limitations persist in areas of strong scour and at coarser scales where wind-terrain interactions are more complex. These findings provide a framework for more efficient snow depth prediction and offer insights to improve snow-vegetation feedback representation in Earth System Models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

ALEGRA Parallel Scaling for Shock in a Heterogeneous Structure

We investigate the strong and weak parallel scaling performance of the ALEGRA multiphysics finite element program when solving a problem involving shock propagation through a heterogeneous material. We determine that ALEGRA scales well over a wide range of problem sizes, cores, and element sizes, and that scaling generally improves as the minimum element size in the mesh increases.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Symmetric Mass Generation with Four SU(2) Doublet Fermions

We study a single exactly massless staggered fermion in the fundamental representation of an SU(2) gauge group. We utilize an nHYP-smeared fermion action supplemented with additional heavy Pauli-Villars fields that serve to decrease lattice artifacts. The phase diagram exhibits a clear two-phase structure with a conformal phase at weak coupling and a novel new phase, the symmetric mass generation (SMG) phase, appearing at strong coupling. The SMG phase is confining with all states gapped and chiral symmetry unbroken. Our finite size scaling analysis provides strong evidence that the phase transition between these two phases is continuous, which would allow for the existence of a continuum SMG phase. Furthermore, the renormalization group flows are consistent with a 𝛽 function that vanishes quadratically at the new fixed point, suggesting that the 𝑁 𝑓 = 4 flavor SU(2) gauge theory lies at the opening of the conformal window.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Extending the Dynamic Range of Galaxy Outflow Scaling Relations: Massive Compact Galaxies with Extreme Outflows

We investigate galactic winds in the HizEA galaxies, a collection of 46 late-stage galaxy mergers at z = 0.4–0.8, with stellar masses of $\mathrm{log}({M}_{* }/{M}_{\odot })=10.4\mbox{--}11.5$ , star formation rates (SFRs) of 20–500 M ⊙ yr -1 , and ultra-compact (a few 100 pc) central star-forming regions. We measure their gas kinematics using the Mg II λλ 2796,2803 absorption lines in optical spectra from MMT, Magellan, and Keck. We find evidence of outflows in 90% of targets, with maximum outflow velocities of 550–3200 km s -1 . We combine these data with ten samples from the literature to construct scaling relations for outflow velocity versus SFR, star formation surface density (Σ SFR ), M * , and SFR/M * . The HizEA galaxies extend the dynamic range of the scaling relations by a factor of ~2–4 in outflow velocity and an order of magnitude in SFR and Σ SFR . The ensemble scaling relations exhibit strong correlations between outflow velocity, SFR, SFR/R, and Σ SFR , and weaker correlations with M * and SFR/M * . The HizEA galaxies are mild outliers on the SFR and M* scaling relations, but they connect smoothly with more typical star-forming galaxies on plots of outflow velocity versus SFR/R and Σ SFR . These results provide further evidence that the HizEA galaxies' exceptional outflow velocities are a consequence of their extreme star formation conditions rather than hidden black hole activity, and they strengthen previous claims that Σ SFR is one of the most important properties governing the velocities of galactic winds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Smoother Surfaces Enhance Diffusion of Nanorods in Entangled Polymer Melts

Coarse-grained molecular dynamics simulations are used to study the diffusion of thin nanorods in entangled polymer melts for varying nanorod length and roughness. While prior studies observed a nanorod parallel diffusion constant scaling inversely with rod length D ∥ ~ l –1 , here, we show that this scaling is not universal and depends sensitively on the nanorod surface roughness. We observe D ∥ ~ l –k , where k < 1 and decreases with decreasing surface roughness. The weaker scaling is driven by the non-Gaussian diffusion of nanorods due to the emergence of an intermittent hopping process that becomes more pronounced with decreasing roughness at the monomer scale. Analysis shows that the mean hop size grows for smoother rods but shows little to no variation with rod length. The mean hopping frequency shows no dependence on either rod length or roughness, suggesting it originates from the polymer melt environment. Further, our results show that the small-scale features of the nanorod surface strongly influence the large-scale and long-time transport of nanorods in polymer matrices, creating new material design opportunities for precisely engineered nanocomposites.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Montana Community-Scale Solar Strategy Project (Final Technical Report)

Grid-tied solar photovoltaic (PV) development began gradually in Montana, initially taking root in the early 2000s, largely in the form of small residential systems, rarely larger than two kilowatts. Utility support for solar technology through public benefit funding triggered the installation of the state’s first solar PV installations on community-owned facilities including schools, libraries, fire stations and local government buildings in the mid-2000s. By 2016, solar PV installations totaled approximately 28 megawatts of installed capacity, generating less than 0.05% of electricity sold in Montana. Montana’s solar fleet in 2016 included 17 megawatts of utility-scale solar farms, and three shared solar projects developed by rural electric cooperatives that ranged from 25 to 50 kilowatts in size. However, with only 3,000 households utilizing solar PV, either through an on-site installation or participation in a shared solar project, the potential benefits and reach of solar technology was limited. The Montana Solar Community Project (MSCP) was developed by the Montana Energy Office (MEO) at the Department of Environmental Quality in order to implement a “community-scale” solar energy strategy for the state that would cost-effectively expand Montanans’ access to solar PV. This project defines three types of community-scale solar projects: shared solar, community-sited solar, and group purchasing programs. The MSCP project conducted research and stakeholder engagement that helped to inform and develop model community-scale initiatives and an overarching community-scale solar strategy for Montana. Specifically, the project included a solar market assessment, stakeholder meetings across Montana, development of a menu of community-scale solar options, a study of solar potential on schools, community solar market research, and community-scale project implementation support. MSCP activities were successful in implementing solar projects and supporting solar development in communities across Montana. Project activities realized 280 kW of newly installed solar from two solarize programs, a shared solar array at a rural electric cooperative, and a solar array atop a Bozeman high school. The community-scale project implementation support and Solar on Schools components of the project resulted directly in the analysis and design of nearly 2 MW of distributed solar systems across the state. The Solar on Schools analysis alone produced solar designs for 25 schools, representing 920.8 kW of solar potential. Many of the deliverables from this project provide guidance on future development of community-scale solar in Montana, and resulted in a number of important conclusions about community-scale solar in Montana. There is strong interest in community-scale solar in Montana, as noted during stakeholder meetings, in response to solar market research surveys, and as seen by applications for implementation support grants provided by MEO. Project activities also suggest that Montanans demonstrated a desire and need for more information about solar energy technology, solar financing options, and the state’s policy and regulatory framework. Additionally, the project found that schools will be key partners for future community-scale solar programs. In addition to being excellent hosts for solar arrays, schools have the opportunity to use the installation of a solar array to develop educational materials and curricula to educate Montana’s future leaders about these technologies, which can be shared and replicated across the state. Project activities laid a clear path forward for solar energy in Montana. Certain activities have concrete next steps. For example, numerous solar arrays were designed for schools and other public buildings around the state. Many of those designs will be used as soon as funding becomes available. Other activities have less developed next steps, though lay the foundation for future work. The solar market assessment, solar market research, and community stakeholder meetings are good examples. These documents can be used as stepping stones to increase renewable energy education, increase stakeholder outreach, and support communities as they seek to develop community-scale solar projects.

14 SOLAR ENERGY↗
DOE OSTI2024

Harmonizing tau positron emission tomography in Alzheimer's disease: The CenTauR scale and the joint propagation model

Abstract INTRODUCTION Tau‐positron emission tomography (PET) outcome data of patients with Alzheimer's disease (AD) cannot currently be meaningfully compared or combined when different tracers are used due to differences in tracer properties, instrumentation, and methods of analysis. METHODS Using head‐to‐head data from five cohorts with tau PET radiotracers designed to target tau deposition in AD, we tested a joint propagation model (JPM) to harmonize quantification (units termed “CenTauR” [CTR]). JPM is a statistical model that simultaneously models the relationships between head‐to‐head and anchor point data. JPM was compared to a linear regression approach analogous to the one used in the amyloid PET Centiloid scale. RESULTS A strong linear relationship was observed between CTR values across brain regions. Using the JPM approach, CTR estimates were similar to, but more accurate than, those derived using the linear regression approach. DISCUSSION Preliminary findings using the JPM support the development and adoption of a universal scale for tau‐PET quantification. Highlights Tested a novel joint propagation model (JPM) to harmonize quantification of tau PET. Units of common scale are termed “CenTauRs”. Tested a Centiloid‐like linear regression approach. Using five cohorts with head‐to‐head tau PET, JPM outperformed linearregressionbased approach. Strong linear relationship was observed between CenTauRs values across brain regions.

Neurosciences & Neurology↗
DOE OSTI2025

Advancing Quantum Many-Body GW Calculations on Exascale Supercomputing Platforms

Advanced ab initio materials simulations face growing challenges as increasing systems and phenomena complexity requires higher accuracy, driving up computational demands. Quantum many-body GW methods are state-of-the-art for treating electronic excited states and couplings but often hindered due to the costly numerical complexity. Here, we present innovative implementations of advanced GW methods within the BerkeleyGW package, enabling large-scale simulations on Frontier and Aurora exascale platforms. Our approach demonstrates exceptional versatility for complex heterogeneous systems with up to 17,574 atoms, along with achieving true performance portability across GPU architectures. We demonstrate excellent strong and weak scaling to thousands of nodes, reaching double-precision core-kernel performance of 1.069 ExaFLOP/s on Frontier (9,408 nodes) and 707.52 PetaFLOP/s on Aurora (9,600 nodes), corresponding to 59.45% and 48.79% of peak, respectively. Our work demonstrates a breakthrough in utilizing exascale computing for quantum materials simulations, delivering unprecedented predictive capabilities for rational designs of future quantum technologies.

Zhang, Benran [University of Southern California, ↗
DOE OSTI2023

Exploring the parameter space of MagLIF implosions using similarity scaling. II. Current scaling

Magnetized liner inertial fusion (MagLIF) is a magneto-inertial-fusion (MIF) concept, which is presently being studied on the Z pulsed power facility. The MagLIF platform has achieved interesting plasma conditions at stagnation and produced significant fusion yields in the laboratory. Given the relative success of MagLIF, there is a strong interest to scale the platform to higher peak currents. However, scaling MagLIF is not entirely straightforward due to the large dimensionality of the experimental input parameter space and the numerous physical processes involved in MIF implosions. Here, in this work, we propose a novel method to scale MagLIF loads to higher currents. Our method is based on similarity (or similitude) scaling and attempts to preserve much of the physics regimes already known or being studied on today's Z pulsed-power driver. By avoiding significant deviations into unexplored and/or less well-understood regimes, the risk of unexpected outcomes on future scaled-up experiments is reduced. Using arguments based on similarity scaling, we derive the scaling rules for the experimental input parameters characterizing a MagLIF load (as functions of the characteristic current driving the implosion). We then test the estimated scaling laws for various metrics measuring performance against results of 2D radiation–magneto-hydrodynamic hydra simulations. Agreement is found between the scaling theory and the simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Scaling neural simulations in STACS

Abstract As modern neuroscience tools acquire more details about the brain, the need to move towards biological-scale neural simulations continues to grow. However, effective simulations at scale remain a challenge. Beyond just the tooling required to enable parallel execution, there is also the unique structure of the synaptic interconnectivity, which is globally sparse but has relatively high connection density and non-local interactions per neuron. There are also various practicalities to consider in high performance computing applications, such as the need for serializing neural networks to support potentially long-running simulations that require checkpoint-restart. Although acceleration on neuromorphic hardware is also a possibility, development in this space can be difficult as hardware support tends to vary between platforms and software support for larger scale models also tends to be limited. In this paper, we focus our attention on Simulation Tool for Asynchronous Cortical Streams (STACS), a spiking neural network simulator that leverages the Charm++ parallel programming framework, with the goal of supporting biological-scale simulations as well as interoperability between platforms. Central to these goals is the implementation of scalable data structures suitable for efficiently distributing a network across parallel partitions. Here, we discuss a straightforward extension of a parallel data format with a history of use in graph partitioners, which also serves as a portable intermediate representation for different neuromorphic backends. We perform scaling studies on the Summit supercomputer, examining the capabilities of STACS in terms of network build and storage, partitioning, and execution. We highlight how a suitably partitioned, spatially dependent synaptic structure introduces a communication workload well-suited to the multicast communication supported by Charm++. We evaluate the strong and weak scaling behavior for networks on the order of millions of neurons and billions of synapses, and show that STACS achieves competitive levels of parallel efficiency.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Lidar-Based Evaluation of HRRR Performance in California’s Diablo Range

The performance of the NOAA High-Resolution Rapid Refresh (HRRR) model for capturing low-level winds near a wind energy production site during summer 2019 is evaluated. This study catalogs the ability of HRRR to predict boundary layer dynamics relevant to wind energy interests over complex terrain, which has presented challenges for weather and energy forecasting. Performance is evaluated by comparing HRRR output to wind-profiling Doppler lidars at Lawrence Livermore National Laboratory Site 300. HRRR captured the diurnal profile of horizontal winds in the observed 150-m layer, despite strong underpredictions (∼4 m s −1 ) during evening and nighttime hours. These underpredictions may be a result of local speedup flows observed by the lidars, which were unresolved in HRRR due to their small spatial extent. HRRR bias magnitude relative to observations was found to be minimal during days with synoptic-scale troughs and strong 850-hPa geopotential gradients, while bias magnitude was maximal during days with synoptic ridging and weak 850-hPa geopotential gradients. To translate wind speed predictions to energy forecasting, generic turbine models were used to estimate power generation for turbines characteristic of the nearby Altamont Pass Wind Resource Area. Results show that HRRR-based energy estimates predicted daytime power generation adequately relative to lidar-based estimates with an 18-h lead time (bias magnitude < 0.4 MW from 0900 to 1400 LT) but overpredicted power during the rest of the diurnal cycle (bias > 1 MW). These results demonstrate conditions under which HRRR performs well for wind energy applications in complex terrain, while highlighting biases that require further investigation to support usage of a high-resolution model for wind energy forecasts.

Boundary layer↗
DOE OSTI2022

Turbulent Regimes in Collisions of 3D Alfvén-wave Packets

Using three-dimensional gyrofluid simulations, we revisit the problem of Alfvén-wave (AW) collisions as building blocks of the Alfvénic turbulent cascade and their interplay with magnetic reconnection at magnetohydrodynamic (MHD) scales. Depending on the large-scale value of the nonlinearity parameter χ 0 (the ratio between the AW linear propagation time and nonlinear turnover time), different regimes are observed. For strong nonlinearities (χ 0 ~ 1), turbulence is consistent with a dynamically aligned, critically balanced cascade—fluctuations exhibit a scale-dependent alignment $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1/4}$, resulting in a ${k}_{\perp }^{-3/2}$ spectrum and ${k}_{\parallel }\propto {k}_{\perp }^{1/2}$ spectral anisotropy. At weaker nonlinearities (small χ 0 ), a spectral break marking the transition between a large-scale weak regime and a small-scale ${k}_{\perp }^{-11/5}$ tearing-mediated range emerges, implying that dynamic alignment occurs also for weak nonlinearities. At χ 0 < 1 the alignment angle ${\theta }_{{k}_{\perp }}$ shows a stronger scale dependence than in the χ 0 ~ 1 regime, namely $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1/2}$ at χ 0 ~ 0.5, and $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1}$ at χ 0 ~ 0.1. Dynamic alignment in the weak regime also modifies the large-scale spectrum, scaling approximately as ${k}_{\perp }^{-3/2}$ for χ 0 ~ 0.5 and as ${k}_{\perp }^{-1}$ for χ 0 ~ 0.1. A phenomenological theory of dynamically aligned turbulence at weak nonlinearities that can explain these spectra and the transition to the tearing-mediated regime is provided; at small χ 0 , the strong scale dependence of the alignment angle combines with the increased lifetime of turbulent eddies to allow tearing to onset and mediate the cascade at scales that can be larger than those predicted for a critically balanced cascade by several orders of magnitude. Such a transition to tearing-mediated turbulence may even supplant the usual weak-to-strong transition.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Scaled Vecchia Approximation for Fast Computer-Model Emulation

Many scientific phenomena are studied using computer experiments consisting of multiple runs of a computer model while varying the input settings. Gaussian processes (GPs) are a popular tool for the analysis of computer experiments, enabling interpolation between input settings, but direct GP inference is computationally infeasible for large datasets. We adapt and extend a powerful class of GP methods from spatial statistics to enable the scalable analysis and emulation of large computer experiments. Specifically, we apply Vecchia’s ordered conditional approximation in a transformed input space, with each input scaled according to how strongly it relates to the computer-model response. The scaling is learned from the data by estimating parameters in the GP covariance function using Fisher scoring. Our methods are highly scalable, enabling estimation, joint prediction, and simulation in near-linear time in the number of model runs. In several numerical examples, our approach substantially outperformed existing methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Towards Scaling Law Analysis For Spatiotemporal Weather Data

Compute-optimal scaling laws are relatively well studied for NLP and CV, where objectives are typically single-step and targets are comparatively homogeneous. Weather forecasting is harder to characterize in the same framework: autoregressive rollouts compound errors over long horizons, outputs couple many physical channels with disparate scales and predictability, and globally pooled test metrics can disagree sharply with per-channel, late-lead behavior implied by short-horizon training. We extend neural scaling analysis for autoregressive weather forecasting from single-step training loss to long rollouts and per-channel metrics. We quantify (1) how prediction error is distributed across channels and how its growth rate evolves with forecast horizon, (2) if power law scaling holds for test error, relative to rollout length when error is pooled globally, and (3) how that fit varies jointly with horizon and channel for parameter, data, and compute-based scaling axes. We find strong cross-channel and cross-horizon heterogeneity: pooled scaling can look favorable while many channels degrade at late leads. We discuss implications for weighted objectives, horizon-aware curricula, and resource allocation across outputs.

Kiefer Jr, Alexander [ORNL] (ORCID:000000025398874↗
DOE OSTI2022

The Astronomy Commons Platform: A Deployable Cloud-based Analysis Platform for Astronomy

Abstract We present a scalable, cloud-based science platform solution designed to enable next-to-the-data analyses of terabyte-scale astronomical tabular data sets. The presented platform is built on Amazon Web Services (over Kubernetes and S3 abstraction layers), utilizes Apache Spark and the Astronomy eXtensions for Spark for parallel data analysis and manipulation, and provides the familiar JupyterHub web-accessible front end for user access. We outline the architecture of the analysis platform, provide implementation details and rationale for (and against) technology choices, verify scalability through strong and weak scaling tests, and demonstrate usability through an example science analysis of data from the Zwicky Transient Facility’s 1Bn+ light-curve catalog. Furthermore, we show how this system enables an end user to iteratively build analyses (in Python) that transparently scale processing with no need for end-user interaction. The system is designed to be deployable by astronomers with moderate cloud engineering knowledge, or (ideally) IT groups. Over the past 3 yr, it has been utilized to build science platforms for the DiRAC Institute, the ZTF partnership, the LSST Solar System Science Collaboration, and the LSST Interdisciplinary Network for Collaboration and Computing, as well as for numerous short-term events (with over 100 simultaneous users). In a live demo instance, the deployment scripts, source code, and cost calculators are accessible. 4 4 http://hub.astronomycommons.org/

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Energy balance and Alfvén Mach numbers in compressible magnetohydrodynamic turbulence with a large-scale magnetic field

Energy equipartition is a powerful theoretical tool for understanding astrophysical plasmas. It is invoked, for example, to measure magnetic fields in the interstellar medium (ISM), as evidence for small-scale turbulent dynamo action, and, in general, to estimate the energy budget of star-forming molecular clouds. Here, in this study, we motivate and explore the role of the volume-averaged root-mean-squared (rms) magnetic coupling term between the turbulent, δB , and large-scale, B 0 , fields, $\langle$(δB·B 0 ) 2 $\rangle$$_{V}^{1/2}$. By considering the second moments of the energy balance equations we show that the rms coupling term is in energy equipartition with the volume-averaged turbulent kinetic energy for turbulence with a sub-Alfvénic large-scale field. Under the assumption of exact energy equipartition between these terms, we derive relations for the magnetic and coupling term fluctuations, which provide excellent, parameter-free agreement with time-averaged data from 280 numerical simulations of compressible magnetohydrodynamic (MHD) turbulence. Furthermore, we explore the relation between the turbulent mean field and total Alfvén Mach numbers, and demonstrate that sub-Alfvénic turbulence can only be developed through a strong, large-scale magnetic field, which supports an extremely super-Alfvénic turbulent magnetic field. This means that the magnetic field fluctuations are significantly subdominant to the velocity fluctuations in the sub-Alfvénic large-scale field regime. Throughout our study, we broadly discuss the implications for observations of magnetic fields and understanding the dynamics in the magnetized ISM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

On the Jet–Ejecta Interaction in 3D GRMHD Simulations of a Binary Neutron Star Merger Aftermath

Short γ-ray burst (sGRB) jets form in the aftermath of a neutron star merger, drill through disk winds and dynamical ejecta, and extend over four to five orders of magnitude in distance before breaking out of the ejecta. We present the first 3D general-relativistic magnetohydrodynamic sGRB simulations to span this enormous scale separation. They feature three possible outcomes: jet+cocoon, cocoon, and neither. Typical sGRB jets break out of the dynamical ejecta if (i) the bound ejecta's isotropic equivalent mass along the pole at the time of the BH formation is ≲10 -4 M ⊙ , setting a limit on the delay time between the merger and BH formation, otherwise, the jets perish inside the ejecta and leave the jet-inflated cocoon to power a low-luminosity sGRB; (ii) the postmerger remnant disk contains a strong large-scale vertical magnetic field, ≳1015 G; and (iii) if the jets are weak (≲10 50 erg), the ejecta's isotropic equivalent mass along the pole must be small (≲10 -2 M ⊙ ). Generally, the jet structure is shaped by the early interaction with disk winds rather than the dynamical ejecta. As long as our jets break out of the ejecta, they retain a significant magnetization (≲1), suggesting that magnetic reconnection is a fundamental property of sGRB emission. The angular structure of the outflow isotropic equivalent energy after breakout consistently features a flat core followed by a steep power-law distribution (slope ≳3), similar to hydrodynamic jets. In the cocoon-only outcome, the dynamical ejecta broadens the outflow angular distribution and flattens it (slope ~1.5).

Jets↗