Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sim model”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Multibody for Everybody (M4E) - A Linearization Approach to Enable Frequency Domain Analysis, Time Integration and Control Co-Design

1.1 Background/Objectives: Marine energy represents a promising yet underexploited source of power. To increase the harvested power, significant efforts have been made to improve wave energy converter (WEC) modeling capabilities and optimize power take-off (PTO) performance; however, these efforts have often treated WEC dynamics, PTO design, and controller development sequentially. In contrast, control co-design (CCD) is emerging as a promising strategy to address these issues directly, creating a growing need for fast analysis tools suitable for repeated simulation and parametric studies [1]. To support this need, this work presents the Multibody for Everybody (M4E) [2] linearization module, which employs a symbolic toolbox to provide deeper insight of WEC design parameters. The objective is to demonstrate that a minimal-coordinate linearization of articulated WEC dynamics can provide accurate wave response predictions and substantial computational savings relative to nonlinear time-domain simulation, while preserving compatibility with broader wave-energy analysis workflows, enabling CCD. 1.2 Approach/Activities: The proposed approach linearizes the equations of motion, generated by M4E, in minimal coordinates about a selected operating point and combines the resulting system with frequencydomain hydrodynamic terms to incorporate the reduced mass, damping, stiffness, and forcing operators. The linearized model is used for both impedance-based response amplitude operator (RAO) prediction and rapid regular-wave time integration. The methodology is demonstrated on a single-flap device and a FOSWEC configuration, with linearized M4E responses compared against the corresponding nonlinear M4E simulations and WEC-Sim results. Regular-wave time histories, RAO trends, and runtime differences are assessed. The framework is also compatible with broader wave-energy workflows, including coupling to WecOptTool, although that capability is not the focus of this work [3]. 1.3 Results/Lessons: The linearized M4E model reproduces key regularwave response characteristics such as integration and Response Amplitude over multiple frequencies. This module matches nonlinear M4E and WEC-Sim results while substantially reducing integration cost. Thus, the proposed framework can serve as a rapid analysis layer for articulated WEC design, parameter studies, and controls-oriented workflows. The analysis is most appropriate in the near-equilibrium regime, about the linearization point.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

A new constraint on galaxy–halo connections of [O ii ] emitters via HOD modelling with angular clustering and luminosity functions from the Subaru HSC survey

ABSTRACT Establishing a robust connection model between emission-line galaxies (ELGs) and their host dark haloes is of paramount importance in anticipation of upcoming redshift surveys. We propose a novel halo occupation distribution (HOD) framework that incorporates galaxy luminosity, a key observable reflecting ELG star-formation activity, into the galaxy occupation model. This innovation enables prediction of galaxy luminosity functions (LFs) and facilitates joint analyses using both angular correlation functions (ACFs) and LFs. Using physical information from luminosity, our model provides more robust constraints on the ELG–halo connection compared to methods relying solely on ACF and number density constraints. Our model was applied to $\rm [O\, {\small II}]$-emitting galaxies observed at two redshift slices at $z=1.193$ and 1.471 from the Subaru Hyper Suprime-Cam PDR2. Our model effectively reproduces observed ACFs and LFs observed in both redshift slices. Compared to the established Geach et al. HOD model, our approach offers a more nuanced depiction of ELG occupation across halo mass ranges, suggesting a more realistic representation of ELG environments. Our findings suggest that ELGs at $z\sim 1.4$ may evolve into Milky-Way-like galaxies, as their inferred halo masses evolve accordingly based on the extended Press–Schechter formalism, highlighting their role as potential building blocks in galaxy formation scenarios. By incorporating the LF as a constraint linking galaxy luminosity to halo properties, our HOD model provides a more precise understanding of ELG-host halo relationships. Furthermore, this approach facilitates the generation of high-quality ELG mock catalogues for future surveys. As the LF is a fundamental observable, our framework is potentially applicable to diverse galaxy populations, offering a versatile tool for analysing data from next-generation galaxy surveys.

Ishikawa, Shogo (ORCID:0000000221184211)↗
DOE OSTI2024

Atomistic simulations for investigation of substrate and salt effects on lipid in-source fragmentation in secondary ion mass spectrometry: A follow-up study

In-source fragmentation (ISF) poses a significant challenge in secondary ion mass spectrometry (SIMS). These fragment ions increase the spectral complexity and can lead to incorrect annotation of fragments as intact species. The presence of salt that is ubiquitous in biological samples can influence the fragmentation and ionization of analytes in a significant manner, but their influences on SIMS have not been well characterized. To elucidate the effect of substrates and salt on ISF in SIMS, we have employed experimental SIMS in combination with atomistic simulations of a sphingolipid on a gold surface with various NaCl concentrations as a model system. Our results revealed that a combination of bond dissociation energy and binding energy between N-palmitoyl-sphingomyelin and a gold surface is a good predictor of fragment ion intensities in the absence of salt. However, ion-fragment interactions play a significant role in determining fragment yields in the presence of salt. Additionally, the charge distribution on fragment species may be a major contributor to the varying effects of salt on fragmentation. This study demonstrates that atomistic modeling can help predict ionization potential when salts are present, providing insights for more accurate interpretations of complex biological spectra.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Reducing the matrix effect in mass spectral imaging of biofilms using flow-cell culture

The interactions between soil microorganisms and soil minerals play a crucial role in the formation and evolution of minerals and the stability of soil aggregates. Due to the heterogeneity and diversity of the soil environment, the under-standing of the functions of bacterial biofilms in soil minerals at the microscale is limited. A soil mineral-bacterial biofilm system was used as a model in this study, and it was analyzed by time-of-flight secondary ion mass spectrometry (ToF-SIMS) to acquire molecular level information. Static culture in multi-wells and dynamic flow-cell culture in microfluidics of biofilms were investigated. Our results show that more characteristic molecules of biofilms can be observed in SIMS spectra of the flow-cell culture. In contrast, biofilm signature peaks are buried under the mineral components in SIMS spectra in the static culture case. Spectral overlay was used in peak selection prior to performing Principal component analysis (PCA). Comparisons of the PCA results between the static and flow-cell culture show more pronounced molecular features and higher loadings of organic peaks of the dynamic cultured specimens. For example, fatty acids secreted from bacterial biofilm extracellular polymeric substance are likely to be responsible for biofilm dispersal due to mineral treatment up to 48 h. Such findings suggest that the use of microfluidic cells to dynamically culture biofilms be a more suitable method for reducing the matrix effect arisen from the growth medium and minerals as a perturbation fac-tor for improved spectral and multivariate analysis of complex mass spectral data in ToF-SIMS. These results show that the interaction mechanism between biofilms and soil minerals at the molecular level can be better studied using the flow-cell culture and advanced mass spectral imaging techniques like ToF-SIMS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Shadows in the Dark: Low-surface-brightness Galaxies Discovered in the Dark Energy Survey

We present a catalog of 23,790 extended low-surface-brightness galaxies (LSBGs) identified in $\sim 5000\,{\deg }^{2}$ from the first three years of imaging data from the Dark Energy Survey (DES). Based on a single-component Sérsic model fit, we define extended LSBGs as galaxies with g-band effective radii ${R}_{\mathrm{eff}}(g)\gt 2\buildrel{\prime\prime}\over{.} 5$ and mean surface brightness ${\bar{\mu }}_{\mathrm{eff}}(g)\gt 24.2\,\mathrm{mag}\,{\mathrm{arcsec}}^{-2}$. We find that the distribution of LSBGs is strongly bimodal in (g-r) versus (g-i) color space. We divide our sample into red (g- i ≥ 0.60) and blue (g- i < 0.60) galaxies and study the properties of the two populations. Redder LSBGs are more clustered than their blue counterparts and are correlated with the distribution of nearby (z < 0.10) bright galaxies. Red LSBGs constitute ~33% of our LSBG sample, and $\sim 30 \% $ of these are located within 1° of low-redshift galaxy groups and clusters (compared to ~8% of the blue LSBGs). For nine of the most prominent galaxy groups and clusters, we calculate the physical properties of associated LSBGs assuming a redshift derived from the host system. In these systems, we identify 41 objects that can be classified as ultradiffuse galaxies, defined as LSBGs with projected physical effective radii ${R}_{\mathrm{eff}}\gt 1.5\,\mathrm{kpc}$ and central surface brightness ${\mu }_{0}(g)\gt 24.0\,\mathrm{mag}\,{\mathrm{arcsec}}^{-2}$. The wide-area sample of LSBGs in DES can be used to test the role of environment on models of LSBG formation and evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Multiwavelength Observations of 2HWC J1928+177: Dark Accelerator or New TeV Gamma-Ray Binary?

2HWC J1928+177 is a Galactic TeV gamma-ray source detected by the High Altitude Water Cherenkov (HAWC) Observatory up to ~56 TeV. The HAWC source, later confirmed by the High Energy Stereoscopic System, still remains unidentified as a dark accelerator since there is no apparent supernova remnant or pulsar wind nebula detected in the lower-energy bands. The radio pulsar PSR J1928+1746, coinciding with the HAWC source position, has no X-ray counterpart. In this work, our SED modeling shows that inverse Compton scattering in the putative pulsar wind nebula can account for the TeV emission only if the unseen nebula is extended beyond $r\sim 4^{\prime} $. Alternatively, TeV gamma-rays may be produced by hadronic interactions between relativistic protons from an undetected supernova remnant associated with the radio pulsar and a nearby molecular cloud G52.9+0.1. NuSTAR and Chandra observations detected a variable X-ray point source within the HAWC error circle, potentially associated with a bright infrared (IR) source. The X-ray spectra can be fitted with an absorbed power-law model with N H = (9 ± 3) × 10 22 cm –2 and Γ X = 1.6 ± 0.3 and exhibit long-term X-ray flux variability over the last decade. If the X-ray source, possibly associated with the IR source (likely an O star), is the counterpart of the HAWC source, it may be a new TeV gamma-ray binary powered by collisions between the pulsar wind and stellar wind. Follow-up X-ray observations are warranted to search for diffuse X-ray emission and determine the nature of the HAWC source.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Constraints on Inflationary Gravitational Waves with Two Years of SPT-3G Data

We present a measurement of the $B$-mode polarization power spectrum of the cosmic microwave background anisotropies at 32 $\le$$\ell$$<$ 502 for three bands centered at 95, 150, and 220 GHz using data from the SPT-3G receiver on the South Pole Telescope. This work uses SPT-3G observations from the 2019 and 2020 winter observing seasons of a $\sim$1500 deg$^2$ patch of sky that directly overlaps with fields observed with the BICEP/Keck family of telescopes, and covers part of the proposed Simons Observatory and CMB-S4 deep fields. Employing new techniques for mitigating polarized atmospheric noise, the SPT-3G data demonstrates a white noise level of 9.3 (6.7) $\mu$K-arcmin at $\ell \sim 500$ for the 95 GHz (150 GHz) data, with a $1/\ell$ noise knee at $\ell$=128 (182). We fit the observed six auto- and cross-frequency $B$-mode power spectra to a model including lensed $\Lambda$CDM $B$-modes and a combination of Galactic and extragalactic foregrounds. This work characterizes foregrounds in the vicinity of the BICEP/Keck survey area, finding foreground power consistent with that reported by the BICEP/Keck collaboration within the same region, and a factor of $\sim$ 3 higher power over the full SPT-3G survey area. Using SPT-3G data over the BICEP/Keck survey area, we place a 95% upper limit on the tensor-to-scalar ratio of $r < 0.25$ and find the statistical uncertainty on $r$ to be $\sigma(r) = 0.067$.

Zebrowski, J. A. [Chicago U., KICP; Chicago U., As↗
DOE OSTI2020

Insights into Na + Diffusion in Silicon Modules under Operating Conditions: Measuring Low Concentrations by D-SIMS

Sodium induced shunting under an electric field is a challenging reliability issue in crystalline Si solar modules. The source of this Potential-Induced Degradation of the Shunting type (PID-s) is well understood and its influence on power loss has been intensively studied based on phenomenological models on cell or module level relating the experimental power-loss to stressing parameters (time, temperature, voltage). However, little is known about the Na ion migration kinetics, responsible for PID on a microscopic level, and its quantitative relation to the efficiency degradation. In this paper we present our investigations of sodium ion migration in Ethylene-Vinyl Acetate (EVA) and silicon through Dynamic Secondary Ion Mass Spectroscopy (D-SIMS). Each sample was annealed at field relevant temperatures from 60-90 °C to address typical migration mechanisms of common PV installations. Analysis of the SIMS migration profiles revealed a diffusivity constant D 0,EVA = 0.09 ± 0.14 cm 2 /s and an activation energy E A,EVA = 0.85 ± .04 eV for Na in EVA and diffusivities higher than extrapolated literature values in silicon (D 0,Si = (3.03 ± 2.42)x10 -5 cm 2 /s, and E A,Si = 0.98 ± 0.02 eV). The new insight will be included in a drift-diffusion based degradation model accounting for the partition coefficient across all relevant interfaces. This model can assist in predicting PID-failure in the field based on the given model stack and the diffusion of Na + through each material. Furthermore, this tool can be used for process optimization as well as material selection significantly reducing the cost and time to validate a technology.

14 SOLAR ENERGY↗
DOE OSTI2025

On the optimal sizing of power take-off systems for wave energy converters

The power take-off (PTO) system in a wave energy converter (WEC) is the means by which the energy in an ocean wave is converted into useful energy. There has been a general lack of technology convergence in the wave energy industry and PTO design has been largely device-specific with the design process determined by the WEC developer. An optimized PTO design can significantly affect the efficiency, reliability, performance, and overall cost and viability of a device. There remains a need in the marine energy industry for best practices and efficient design processes to be identified as it could greatly benefit developers and researchers in designing an optimal WEC PTO. This research utilizes the open-source Matlab-Simulink based software, WEC-Sim, which was developed by the National Renewable Energy Laboratory, and the Eagle, a high-performance computing system, to model and optimize two WECs of different archetypes and then investigates if any trends exist that can be exploited for greater efficiency in optimal PTO design. The results from this research indicate that there may be an optimal sizing for WEC PTOs for both power and PTO force rating and that substantial reductions can be made in the cost of a device without significant loss in the energy produced. Furthermore this optimal sizing may be independent of WEC type and deployment location. With an optimal, or near optimal, approach to WEC PTO sizing demonstrated, a methodology is proposed to address the challenge of nameplate ratings in the wave energy industry.

13 HYDRO ENERGY↗
DOE OSTI2023

Vertical control of DIII-D discharges with strong negative triangularity

Through predictive modeling validated by a series of experiments on DIII-D, the vertical stability of low β diverted plasmas with strong negative triangularity (NT) ($\delta\sim-0.6$) is assessed. As a result of their unique magnetic geometry, NT plasmas feature larger Shafranov shifts and more elongated inner flux surfaces than positive triangularity counterparts, typically leading to enhanced vertical instability growth rates. However, coupling with the non-conformal vessel DIII-D wall reduces these growth rates to controllable values, providing a path forward for stabilizing strongly NT plasma with a diverted geometry. These discharges are used to validate GSdesign (part of the TokSys code suite) stability calculations in strong NT plasmas on DIII-D, with errors of no more than ~20% observed between modeled and experimentally measured growth rates. Additions of diagnostic noise and power supply tuning to the TokSys model are needed to accurately capture the time dependence of DIII-D NT discharges, assisting with the design of control schemes specific to the DIII-D poloidal field coils. Lastly, implications of these results on a future NT reactor are briefly described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

An Efficient Time-Domain Model to Simulate Parametric Resonances in a Floating Body Free to Move in Six Degrees of Freedom: Preprint

We present a computationally efficient time-domain model capable of simulating parametric resonances in a floating body in waves. The model assumes all wave forces to be linear, but the inertia and restoring forces acting on the body are expanded to second order in body motions. The simulation speed on a standard computer is approximately 40 times faster than real time. The model is applied to a soft-moored floating axisymmetric body which absorbs energy through heave, but is otherwise free to move in six degrees of freedom. Under certain conditions, we show that the body responds parametrically with large amplitudes not only in surge and pitch, but also in sway, roll, and yaw, provided it is given some small initial displacement in one of these out-of-plane modes. The predictions are confirmed by simulations using state-of-the-art nonlinear Froude-Krylov and computational fluid dynamics models.

parametric resonance↗
DOE OSTI2026

Predicting core transport in ITER baseline discharges with neon injections

Achieving self-consistent performance predictions for ITER requires integrated modeling of core transport and divertor power exhaust under realistic impurity conditions. We present results from a systematic power-flow and impurity-content study for the ITER 15 MA baseline scenario constrained directly by existing SOLPS-ITER neon-seeded divertor solutions. Using the OMFIT STEP workflow, stationary temperature and density profiles are predicted with TGYRO for $1.5 \unicode{x2A7D} Z_\textrm{eff} \unicode{x2A7D} 2.5$, and the corresponding power crossing the separatrix $P_\textrm{sep}$ is evaluated. We find that $P_\textrm{sep}$ varies by more than a factor of 1.7 across this scan and matches the ${\sim}100$ MW SOLPS-ITER prediction when $Z_\textrm{eff} \simeq 1.6$ or when auxiliary heating is reduced to ${\sim}75\%$ of nominal. Rotation-sensitivity studies show that plausible variations in toroidal flow magnitude modify $P_\textrm{sep}$ by $\lesssim 20\%$, while AURORA modeling confirms that charge-exchange radiation inside the separatrix is dynamically negligible under predicted ITER neutral densities. These results identify a restricted compatibility window, $Z_\textrm{eff} \approx 1.6$ –1.75 and $0.75 \lesssim f_{P_\textrm{aux}} \unicode{x2A7D} 1.0$, in which core transport predictions remain aligned with neon-seeded divertor protection targets. This self-consistent, model-constrained framework provides actionable guidance for impurity control and auxiliary-heating scheduling in early ITER operation and supports future whole-device scenario optimization.

ITER↗
DOE OSTI2026

Beyond Point Estimates: Benchmarking Uncertainty Quantification Methods on the AION-1 Astronomical Foundation Model

Foundation models for astronomical surveys offer powerful learned representations that can be transferred to downstream regression tasks such as galaxy property estimation. However, point predictions alone are insufficient for scientific inference; reliable uncertainty quantification (UQ) is essential. We compare seven UQ methods on galaxy property regression using frozen AION-1 foundation-model embeddings, predicting redshift, stellar mass, stellar-population age, gas-phase metallicity, and specific star-formation rate, from Legacy Survey photometry/imaging and DESI spectra, with PROVABGS-derived labels. Distribution-free conformal methods achieve marginal coverage within $\sim$1 pp of the nominal 90% across all properties, while non-conformal baselines (Deep Ensembles, MC~Dropout) fail to calibrate reliably. Among conformal approaches, Conformalized Quantile Regression (CQR) delivers the best coverage in the bin with the poorest model predictions. More importantly, only the Locally Valid and Discriminative (LVD) framework -- particularly when operating on AION-1 embeddings -- also provides finite-sample \emph{local validity}, producing intervals that adapt to each galaxy's local prediction difficulty rather than relying on marginal guarantees alone. These results establish conformal prediction, and LVD in particular, as the preferred UQ framework for uncertainty-aware inference on foundation-model embeddings in astrophysics.

Tame-Narvaez, Karla [Fermilab] (ORCID:000000022249↗
DOE OSTI2025

Search for supersymmetry in final states with missing transverse momentum and charm-tagged jets using 139 fb−1 of proton-proton collisions at s = 13 TeV with the ATLAS detector

The paper presents a search for supersymmetric particles produced in proton-proton collisions at s$$ \sqrt{s} $$ = 13 TeV and decaying into final states with missing transverse momentum and jets originating from charm quarks. The data were taken with the ATLAS detector at the Large Hadron Collider at CERN from 2015 to 2018 and correspond to an integrated luminosity of 139 fb−1. No significant excess of events over the expected Standard Model background expectation is observed in optimized signal regions, and limits are set on the production cross-sections of the supersymmetric particles. Pair production of charm squarks or top squarks, each decaying into a charm quark and the lightest supersymmetric particle χ~10$$ {\overset{\sim }{\chi}}_1^0 $$, is excluded at 95% confidence level for squarks with masses up to 900 GeV for scenarios where the mass of χ~10$$ {\overset{\sim }{\chi}}_1^0 $$ is below 50 GeV. Additionally, the production of leptoquarks with masses up to 900 GeV is excluded for the scenario where up-type leptoquarks decay into a charm quark and a neutrino. Model-independent limits on cross-sections and event yields for processes beyond the Standard Model are also reported.

Experimental Nuclear Physics↗
DOE OSTI2022

Constraining the baryonic feedback with cosmic shear using the DES Year-3 small-scale measurements

ABSTRACT We use the small scales of the Dark Energy Survey (DES) Year-3 cosmic shear measurements, which are excluded from the DES Year-3 cosmological analysis, to constrain the baryonic feedback. To model the baryonic feedback, we adopt a baryonic correction model and use the numerical package baccoemu to accelerate the evaluation of the baryonic non-linear matter power spectrum. We design our analysis pipeline to focus on the constraints of the baryonic suppression effects, utilizing the implication given by a principal component analysis on the Fisher forecasts. Our constraint on the baryonic effects can then be used to better model and ameliorate the effects of baryons in producing cosmological constraints from the next-generation large-scale structure surveys. We detect the baryonic suppression on the cosmic shear measurements with a ∼2σ significance. The characteristic halo mass for which half of the gas is ejected by baryonic feedback is constrained to be $M_c \gt 10^{13.2} \, h^{-1} \, \mathrm{M}_{\odot }$ (95 per cent C.L.). The best-fitting baryonic suppression is $\sim 5{{\ \rm per\ cent}}$ at $k=1.0 \, {\rm Mpc}\ h^{-1}$ and $\sim 15{{\ \rm per\ cent}}$ at $k=5.0 \, {\rm Mpc} \ h^{-1}$. Our findings are robust with respect to the assumptions about the cosmological parameters, specifics of the baryonic model, and intrinsic alignments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

FLAMINGO: combining kinetic SZ effect and galaxy–galaxy lensing measurements to gauge the impact of feedback on large-scale structure

Energetic feedback processes associated with accreting supermassive black holes can expel gas from massive haloes and significantly alter various measures of clustering on $\sim$Mpc scales, potentially biasing the values of cosmological parameters inferred from analyses of large-scale structure (LSS) if not modelled accurately. Here, we use the state-of-the-art FLAMINGO suite of cosmological hydrodynamical simulations to gauge the impact of feedback on large-scale structure by comparing to Planck + ACT stacking measurements of the kinetic Sunyaev–Zel’dovich (kSZ) effect of SDSS BOSS galaxies. We make careful like-with-like comparisons to the observations, aided by high precision KiDS and DES galaxy–galaxy lensing measurements of the BOSS galaxies to inform the selection of the simulated galaxies. In qualitative agreement with several recent studies using dark matter only simulations corrected for baryonic effects, we find that the kSZ effect measurements prefer stronger feedback than predicted by simulations which have been calibrated to reproduce the gas fractions of low redshift X-ray-selected groups and clusters. We find that the increased feedback can help to reduce the so-called $S_8$ tension between the observed and CMB-predicted clustering on small scales as probed by cosmic shear (although at the expense of agreement with the X-ray group measurements). However, the increased feedback is only marginally effective at reducing the reported offsets between the predicted and observed clustering as probed by the thermal SZ (tSZ) effect power spectrum and tSZ effect – weak lensing cross-spectrum, both of which are sensitive to higher halo masses than cosmic shear.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Production and Persistence of Extreme Two-temperature Plasmas in Radiative Relativistic Turbulence

Turbulence is a predominant process for energizing electrons and ions in collisionless astrophysical plasmas, and thus is responsible for shaping their radiative signatures (luminosity, spectra, and variability). To better understand the kinetic properties of a collisionless radiative plasma subject to externally driven turbulence, we investigate particle-in-cell simulations of relativistic plasma turbulence with external inverse Compton cooling acting on the electrons. We find that ions continuously heat up while electrons gradually cool down (due to the net effect of radiation), and hence the ion-to-electron temperature ratio T i /T e grows in time. Furthermore, we show that Ti/Te is limited only by the size and duration of the simulations (reaching ${T}_{i}/{T}_{e}\sim {10}^{3}$), indicating that there are no efficient collisionless mechanisms of electron–ion thermal coupling. This result has implications for models of radiatively inefficient accretion flows, such as observed in the Galactic center and in M87, for which so-called two-temperature plasmas with ${T}_{i}/{T}_{e}\gg 1$ have been invoked to explain their low luminosity. Additionally, we find that electrons acquire a quasi-thermal distribution (dictated by the competition of turbulent particle energization and radiative cooling), while ions undergo efficient nonthermal acceleration (acquiring a harder distribution than in equivalent nonradiative simulations). There is a modest nonthermal population of high-energy electrons that are beamed intermittently in space, time, and direction; these beamed electrons may explain rapid flares in certain high-energy astrophysical systems (e.g., in the Galactic center). These numerical results demonstrate that extreme two-temperature plasmas can be produced and maintained by relativistic radiative turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Fate of supernova progenitors in massive binary systems

ABSTRACT How massive stars end their lives depends on the core mass, core angular momentum, and hydrogen envelopes at death. However, these key physical facets of stellar evolution can be severely affected by binary interactions. In turn, the effectiveness of binary interactions itself varies greatly depending on the initial conditions of the binaries, making the situation much more complex. We investigate systematically how binary interactions influence core–collapse progenitors and their fates. Binary evolution simulations are performed to survey the parameter space of supernova progenitors in solar metallicity binary systems and to delineate major evolutionary paths. We first study fixed binary mass ratios ($q=M_2/M_1$ = 0.5, 0.7, and 0.9) to elucidate the impacts of initial mass and initial separation on the outcomes, treating separately Type Ibc supernova, Type II supernova, accretion-induced collapse (AIC), rapidly rotating supernova (Ibc-R), black hole formation, and long gamma ray burst (long GRB). We then conduct 12 binary population synthesis model calculations, varying the initial condition distributions and binary evolution parameters, to estimate various supernova fractions. We obtain a Milky Way supernova rate $R_{\rm SN} = (1.78$–$2.47) \times 10^{-2} \, {\rm yr}^{-1}$ which is consistent with observations. We find the rates of AIC, Ibc-R, and long GRB to be $\sim 1/100$ the rate of regular supernovae. Our estimated long GRB rates are higher than the observed long GRB rate and close to the low luminosity GRB rate, although care must be taken considering our models are computed with solar metallicity. Furthering binary modelling and improving the inputs one by one will enable more detailed studies of these and other transients associated with massive stars.

Astronomy & Astrophysics↗