Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated emission”

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 361 records · Page 20

Informed Multi‐Scale Approach Applied to the British Columbia Fires of Late Summer 2017

Abstract Pyrocumulonimbus (PyroCb) clouds have a complex origin dependent on fire dynamics and meteorological conditions. When a pyrocumulonimbus cloud develops and is maintained over a period of time, it can inject significant aerosol into the troposphere and lower stratosphere, resulting in a longer‐term (months to years) occurrence of aerosol in the stratosphere. In this work, we investigate the British Columbia wildfires on 12–13 August 2017 using a multi‐scale simulation framework. We use the output of a physics‐based wildfire model (FIRETEC) with parameterized energy, particle, and gas emissions to drive the upper atmospheric aerosol mass injection within a regional cloud resolving model (HIGRAD). We demonstrate that vertical motions produced by latent heat release of the condensation of ice and cloud particles within the PyroCbs induce another 5 km of lifting of the simulated aerosol plume. Primary black carbon and organic aerosols (OAs) alone may not be enough to explain the observed aerosol burden, thus we show that secondary OA produced via condensation of gases by the fires, ash, and possibly dust can enhance lofted aerosol mass. A simulation with all emission mechanisms active, driven by the observed fuel load and environmental conditions, reasonably reproduces an aerosol profile inferred from observational data.

54 ENVIRONMENTAL SCIENCES↗

Monte Carlo Simulations of Polarized Radiative Transfer in Neutron Star Atmospheres

Soft X-ray emission from neutron stars affords powerful diagnostic tools for uncovering their surface and interior properties, as well as their geometric configurations. In the atmospheres of neutron stars, the presence of magnetic fields alters the photon-electron scattering cross sections, resulting in nontrivial angular dependence of intensity and polarization of the emergent signals. This paper presents recent developments of our Monte Carlo simulation, MAGTHOMSCATT, which tracks the complex electric field vector for each photon during its transport. Our analysis encompasses the anisotropy and polarization characteristics of X-ray emission for field strengths ranging from nonmagnetic to extremely magnetized regimes that are germane to magnetars. In the very low field domain, we reproduced the numerical solution to the radiative transfer equation for nonmagnetic Thomson scattering, and provided analytical fits for the angular dependence of the intensity and the polarization degree. These fits can be useful for studies of millisecond pulsars and magnetic white dwarfs. By implementing a refined injection protocol, we show that, in the magnetar regime, the simulated intensity and polarization pulse profiles of emission from extended surface regions becomes invariant with respect to the ratio of photon (ω) and electron cyclotron (ω B ) frequencies once ω/ω B ≲ 0.01. This circumvents the need for simulations pertinent to really high magnetic field strengths, which are inherently slower. Our approach will be employed elsewhere to model observational data to constrain neutron star geometric parameters and properties of emitting hot spots on their surfaces.

79 ASTRONOMY AND ASTROPHYSICS↗

A kinetic study of fusion burn waves in compressed deuterium–tritium and proton–boron plasmas

We present particle-in-cell simulations with Monte Carlo collisions of fusion burn waves in compressed deuterium–tritium and proton–boron plasmas. We study the energy balance in the one-dimensional expansion of a hot-spot by simulating Coulomb collisions, fusion reactions, and bremsstrahlung emission with a Monte Carlo model and inverse bremsstrahlung absorption using a new PIC model. This allows us to self-consistently capture the alpha particle heating and radiative losses in the expanding hot-spot and surrounding cold fuel. After verifying our model in a code-to-code comparison with both kinetic and fluid codes for the case of a deuterium–tritium hot-spot, we simulate the expansion of a proton–boron hot-spot initialized at 200 keV and 1,000 g/cm 3 . Our model predicts that energy radiated by the hot-spot is recaptured by the surrounding high-density opaque fuel reducing the expansion work done by the propagating burn wave. As a result, we find the net fusion energy produced over the course of $20$~ps is twice the initial hot-spot energy independent of whether radiation physics is included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Study of vacuum ultraviolet emission in helium and helium/nitrogen mixtures

Helium is frequently used as a working medium for the generation of plasmas and is capable of energetic photon emissions. These energetic photon emissions are often attributed to the formation of helium excimer and subsequent photon emission. When the plasma device is exposed to another gas, such as nitrogen, this energetic photon emission can cause photoionization and further ionization wave penetration into the additional gas. Often ignored are the helium resonance emissions that are assumed to be radiation trapped and therefore not pertinent to photoionization. Here, experimental evidence for the presence of helium atomic emission in a pulsed discharge at ten’s of Torr is shown. Simulations of a discharge in similar conditions agree with the experimental measurements. In this context, the role of atomic and molecular helium light emission on photoionization of molecular nitrogen in an ionization wave is studied using a kinetic modeling approach that accounts for radiation dynamics in a developing low-temperature plasma. Three different mixtures of helium at a total pressure of 250 Torr are studied in simulation. Photoionization of the nitrogen molecule by vacuum ultraviolet helium emission is used as the only seed source ahead of the ionization front. It is found that even though radiation trapped, the atomic helium emission lines are the significant source of photoionization of nitrogen. The significant effect of radiation trapped photon emission on ionization wave dynamics demonstrates the need to consider these radiation dynamics in plasma reactors where self-absorbed radiation is ignored.

42 ENGINEERING↗

Nexus of electrification and energy efficiency retrofit of commercial buildings at the district scale

Rapid electrification of buildings at the district scale is needed for cities to achieve climate change mitigation goals. However, most electrification studies focus on either the single building level or the city/region building stock level, and depend on the slow and uncertain process of requesting personally identifiable customer energy usage data from utilities. To answer a key question facing local policymakers: “Where can electrification proceed at scale without first upgrading the grid?” this study aims to quantify and inform building electrification impacts at the district scale using detailed building energy modeling and based on public records datasets. We explore how energy efficiency retrofits can help mitigate increased peak electric demand, and quantify impacts to energy use and carbon emissions. Building energy models of a baseline, and scenarios of simple electrification, energy retrofits, and electrification in combination with retrofits were created and simulated for 54 commercial buildings in two contiguous districts of San Francisco. A simple electrification scenario increased annual electricity consumption but reduced annual site energy usage by 15% to 17%, mainly due to replacing inefficient gas furnaces and boilers with more efficient heat pumps. Peak demand increased 7.4% for Fisherman's Wharf (e.g. within the capacity of the existing power grid), while the Design District showed a marginal decrease. Annual carbon emissions were reduced by 46% and 37%. Combining electrification with efficiency upgrades reduced peak demand by 26% and 40%, and annual carbon emissions by 63% and 64% for the two districts. Furthermore, these results indicate that impacts of electrification depend on the mix of building uses within a district, and coupling electrification with energy efficiency upgrades is an effective strategy to decarbonize buildings while maintaining or reducing the peak electric demand.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Fuel–shell mix and yield degradation in kinetic shock-driven inertial confinement fusion implosions

Fuel–shell mix in kinetic plasma conditions is probed using nuclear and x-ray self-emission in shock-driven, D 3 He-gas-filled inertial confinement fusion implosions. As initial gas fill density decreases, measured nuclear yields and ion temperatures are lower than expected as compared to radiation-hydrodynamic simulations. Spatially and temporally resolved x-ray emissions indicate significant mixing at the fuel–shell interface in implosions with low initial gas fill density. In conclusion, this observed fuel–shell mix introduces a substantial amount of shell ions into the center of the implosion prior to and during shock flash and is the key mechanism needed in the kinetic-ion simulations to match experimental nuclear yields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigating relationship between surface topography and emissivity of metallic additively manufactured parts

Due to the direct relationship between thermal history and mechanical behavior, in situ thermal monitoring is key in gauging quality of parts produced with additive manufacturing (AM). Accurate monitoring of temperatures in an AM process requires knowledge of environment and object parameters including object emissivity. The emissivity is dependent on several variables, including: wavelength, material composition, temperature, and surface topography. Researchers have been concerned with the thermal emissivity dependence on temperature since large ranges are seen in metal powder bed processes, but there is also an extensive range of surfaces produced by AM. This work focused on discovering what roughness characteristics control thermal emissivity through investigation of prototypic 316 stainless steel AM samples produced with a range of build conditions on a laser powder bed fusion machine. Through experimental measurements of emissivity using hemispherical directional reflectance (HDR), guided by simulations using a finite-difference time-domain (FDTD) Maxwell solver, it was found that combinations of existing roughness parameters describing both height and slope of the surface correlate well with emissivity changes. Furthermore, these parameters work well due to their apt description of surface features encouraging internal reflection, which is the phenomenon that increases emissivity when a surface falls under the geometric optical region conditions.

36 MATERIALS SCIENCE↗

Temporal evolution of the light emitted by a thin, laser-ionized plasma source

We present an experimental and simulation-based investigation of the temporal evolution of light emission from a thin, laser-ionized helium plasma source. We demonstrate an analytic model to calculate the approximate scaling of the time-integrated, on-axis light emission with the initial plasma density and temperature, supported by the experiment, which enhances the understanding of plasma light measurement for plasma wakefield accelerator (PWFA) plasma sources. Our model simulates the plasma density and temperature using a split-step Fourier code and a particle-in-cell code. A fluid simulation is then used to model the plasma and neutral density, and the electron temperature as a function of time and position. We then show the numerical results of the space-and-time-resolved light emission and that collisional excitation is the dominant source of light emission. We validate our model by measuring the light emitted by a laser-ionized plasma using a novel statistical method capable of resolving the nanosecond-scale temporal dynamics of the plasma light using a cost-effective camera with microsecond-scale timing jitter. This method is ideal for deployment in the high radiation environment of a particle accelerator that precludes the use of expensive nanosecond-gated cameras. Our results show that our models can effectively simulate the dynamics of a thin, laser-ionized plasma source. In addition, this work provides a detailed understanding of the plasma light measurement, which is one of the few diagnostic signals available for the direct measurement of PWFA plasma sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of RF Sheath-Driven Tungsten Erosion at RF Antenna in the WEST Tokamak

This study applies the newly developed STRIPE (Simulated Transport of RF Impurity Production and Emission) framework to analyze tungsten (W) erosion at RF antenna structures in the WEST tokamak. STRIPE integrates SolEdge3x for edge plasma backgrounds, COMSOL for 3D RF sheath potentials, RustBCA for sputtering yields, and GITR for impurity transport and ion energy–angle distributions. Building on prior work by Kumar et al. (2025) Nuclear Fusion, 65, 076039, which validated STRIPE for WEST ICRH discharge #57877, the present study provides a spatially resolved assessment of gross W erosion at both Q2 antenna limiters under ohmic and ICRH conditions. Simulations using 2D SolEdge3x profiles in COMSOL capture rectified sheath potentials exceeding 300 V, leading to strong upper-limiter localization. Both poloidal and toroidal asymmetries are observed and attributed to RF sheath effects, with modeled erosion patterns deviating from experiment—highlighting sensitivity to sheath geometry and plasma resolution. Erosion is driven primarily by high-charge-state oxygen ions (O6+–O8+), while D+ plays a negligible role. Assuming a plasma composition of 1% oxygen and 98% deuterium, STRIPE predicts a 30-fold increase in gross W erosion from ohmic to ICRH phases, consistent with a >25-fold rise in W-I (400.9 nm) brightness. Quantitative agreement is within 5% in the ohmic phase and 30% under ICRH, demonstrating predictive capability. Importantly, the study shows that the magnitude of ICRH-driven W erosion depends strongly on the concentration of light impurities (O, B, N, C), which drive sputtering through high charge states. Cleaner plasma conditions with reduced impurity content are therefore expected to substantially mitigate antenna W sources in WEST and other toroidal fusion devices. These findings establish STRIPE as a predictive framework for RF-induced plasma–material interactions and support its application to reactor-scale antenna design.

Kumar, Atul [ORNL] (ORCID:0000000277475210)↗

Mixing in ICF implosions on the National Ignition Facility caused by the fill-tube

The micrometer-scale tube that fills capsules with thermonuclear fuel in inertial confinement fusion experiments at the National Ignition Facility is also one of the implosion's main degradation sources. It seeds a perturbation that injects the ablator material into the center, radiating away some of the hot-spot energy. This paper discusses how the perturbation arises in experiments using high-density carbon ablators and how the ablator mix interacts once it enters the hot-spot. Both modeling and experiments show an in-flight areal-density perturbation and localized x-ray emission at stagnation from the fill-tube. Simulations suggest that the fill-tube is degrading an otherwise 1D implosion by ~2×, but when other degradation sources are present, the yield reduction is closer to 20%. Characteristics of the fill-tube assembly, such as the through-hole size and the glue mass, alter the dynamics and magnitude of the degradation. Overall, these aspects point the way toward improvements in the design, some of which (smaller diameter fill-tube) have already shown improvements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Inference of molecular divertor density from filtered camera analysis of molecularly induced Balmer line emission during detachment in JET L-mode plasmas

Here a previously presented Monte Carlo method for estimating local plasma conditions in 2D based on intensity ratios of deuterium Balmer D α , D γ and D ε lines was amended to consider also the D α and D γ emission contributions arising from molecular processes. The obtained estimates were used to infer the molecular divertor density with the help of the molecular databases of EIRENE. The method was benchmarked against EDGE2D-EIRENE simulations and observed to reproduce the molecularly induced emission fractions and the molecular divertor densities primarily within 25% of the references. Experimental analysis of a JET L-mode density scan suggested molecularly induced D α and D γ contributions of up to 60–70% and 20%, respectively, during the process of detachment. The independent estimates of the molecular divertor density inferred from the obtained molecularly induced D α and D γ intensities agree within uncertainties with each other. Both estimates show the molecular density increasing up to approximately 1.0–2.0 x 10 20 m -3 at the outer strike point in deep detachment with its ratio to the local electron density agreeing with EDGE2D-EIRENE predictions within the scatter of the experimental data.

interferometry↗

Experimental and numerical investigation of flame stabilization and pollutant formation in matrix stabilized ammonia-hydrogen combustion

Ammonia (NH 3 ) is a carbon-free fuel that offers an attractive alternative for reducing greenhouse gas emissions. However, the slow flame speed, low heating value, and emissions of nitrogen-containing pollutants present significant issues for practical combustion applications. To address these issues, we investigate the use of matrix stabilized combustion. In this type of burner, combustion is performed within an inert porous ceramic foam, heat is recirculated by solid conduction and radiation, which enhances flame speed and combustion stabilization, thereby permitting combustion over a wide range of equivalence ratio conditions. We present a new porous media burner (PMB) capable of stabilizing NH 3 /air flames at ambient conditions. An extensive experimental characterization of the stability of this burner is conducted with up to 30% by volume of hydrogen (H 2 ) in the fuel stream. A 15:1 turndown ratio is demonstrated, with a high thermal power density of 62 MW m -3 . Concentrations of NO, unburnt NH 3 , and H 2 in the exhaust stream are measured. Two regimes are identified for low NO operation: rich and very lean. For rich conditions, NO emissions decrease with increasing equivalence ratio and decreasing H 2 blending. Unburnt NH 3 emissions follow opposite trends. These measurements are complemented by simulations in which the burner is represented by a coupled solid-gas reactor network. In conclusion, this model captures the burner’s pollutant emissions to good accuracy and is used to analyze the mechanisms of pollutant formation.

42 ENGINEERING↗

Stratospheric aerosol forcing for CMIP7 – Part 1: optical properties for pre-industrial, historical, and scenario simulations

Stratospheric aerosols, most of which originate from explosive volcanic sulfur emissions into the stratosphere, are a key natural driver of climate variability. They are thus a forcing provided by the Coupled Model Intercomparison Project (CMIP) Climate Forcings Task Team to climate modelling groups participating in phase 7 of CMIP. For the historical period, we provide two datasets covering 1750–2023: (i) a volcanic upper tropospheric-stratospheric sulfur emission dataset, documented in a companion paper; and (ii) a stratospheric sulfate aerosol optical property dataset, documented in this paper. For the satellite era (from 1979 onwards), stratospheric aerosol optical properties are derived from the Global Space-based Stratospheric Aerosol Climatology (GloSSAC) dataset. For the pre-satellite era (1750–1978), optical properties are derived from our volcanic SO 2 emission dataset using a new version of the reduced-complexity volcanic aerosol model Easy Volcanic Aerosol (Height) (EVA_H v2). A background, non-volcanic stratospheric aerosol climatology is derived from the 1998–2001 period with a trend over 1850–1978 accounting for increasing anthropogenic aerosols. A monthly stratospheric aerosol climatology is derived from the 1850–2021 average for both pre-industrial and scenario (future) simulations, with a 9-year ramp over 2022–2030 for scenario simulations to ensure a smooth transition from the historical period. CMIP7 simulations should use version 2.2.1 of our pre-industrial and historical datasets, and version 2.2.2 of our scenario dataset. Our methodology to produce historical aerosol optical properties significantly differs from CMIP6 for the pre-satellite era, and the resulting forcings in turn largely differ. In particular, the CMIP6 dataset was mostly based on the sparse and uncertain pyrheliometer record, which resulted in strongly underrepresented emissions from small-to-moderate magnitude eruptions. The resulting bias is addressed in CMIP7, which is entirely emission-derived in the pre-satellite era and uses more recent ice-core-based volcanic sulfur emission inventories than CMIP6. Our approach results in an overall larger volcanic aerosol forcing for CMIP7, with the 1850–2014 mean mid-visible global mean stratospheric aerosol optical depth (SAOD) in CMIP7 (0.0138) being 29 % higher than in CMIP6 (0.0107). The pre-industrial mean of the same variable is 26 % higher in CMIP7 (0.0135, derived from the historical 1850–2021) than CMIP6 (0.0107, derived from the historical 1850–2014 mean). Using a reduced-complexity climate model, we simulate a global mean surface temperature that is 0.07 °C colder for 1850–1900 when using the CMIP7 dataset instead of CMIP6, whereas 2000–2014 is 0.03 °C warmer in CMIP7. Our dataset also exhibits lower forcing for 1960–1980, resulting in temperatures 0.06 °C warmer when averaged over 1960–1990, a period for which CMIP6 climate models exhibit a cold bias. Given the large uncertainties characterizing the dataset, in particular for the pre-satellite era, we advise against treating the CMIP7 or CMIP6 dataset as uniquely superior for any specific year and highlight the need for further evaluation. We conclude the study by discussing sources of uncertainty for the dataset, future research avenues to improve it, as well as requirements to operationalize the production of the dataset, i.e. extend it and update it on an annual basis instead of every 5–7 years following CMIP cycles.

Aubry, Thomas J. [Univ. of Exeter, Devon (United K↗

How to support EV adoption: Tradeoffs between charging infrastructure investments and vehicle subsidies in California

Supporting the adoption of zero-emission vehicle (ZEVs), including plug-in electric vehicles (EVs), has become a priority for governments due to their ability to reduce petroleum demand, improve air quality, and reduce carbon dioxide (CO2) emissions. Optimal strategies to accelerate EV adoption must weigh the relative value of alternative policy mechanisms to consumers, including public charging infrastructure and vehicle purchase subsidies. We use a historically validated light-duty vehicle consumer choice tool, the ADOPT model, to simulate personal light-duty vehicle adoption and related emissions in California. ADOPT is updated to incorporate a quantification of the tangible value of public charging infrastructure, allowing us to simulate the impact of investments in public charging infrastructure and vehicle purchase subsidies under different scenarios. We show that both policies result in increased EV adoption, with the most effective policy varying depending on vehicle technology assumptions. Under conservative technology improvement assumptions, infrastructure investments are most effective in promoting EV sales and reducing CO2 emissions, while under optimistic technology improvement assumptions a combination of infrastructure and subsidies best supports EV sales and CO2 emission reductions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Interactive Gas Chemistry for Enhanced Science Capabilities of the Energy Exascale Earth System Model Version 3

Atmospheric chemistry plays a crucial role in Earth system models (ESMs), controlling atmospheric composition and radiative balance; it is highly interactive with the physical climate, biogeochemical cycles, and human systems. However, it often imposes computational challenges in an ESM. Here we develop a full troposphere‐stratosphere interactive chemistry module for the US Department of Energy's Energy Exascale Earth System Model (E3SM). We intentionally build a streamlined module based on E3SM version 2 that interacts with other components and maintains all of major chemical and chemistry‐climate feedbacks. The module incorporates a new, highly efficient tracer advection scheme; linearization of stratospheric chemistry; and abridged tropospheric chemical mechanism with 28 reactive tracers. This new model, E3SM‐chem, can readily perform century‐long climate simulations of ozone, methane, and nitrous oxide based on emission scenarios as well as provide hourly budgets for the gas‐phase radicals that drive aerosol chemistry. We evaluate E3SM‐chem with an atmosphere‐only simulation as in the recent climate model intercomparison project (CMIP6) finding results similar to the other CMIP6 models. For the present‐day, E3SM‐chem matches the standard measurement metrics for stratospheric and tropospheric ozone, surface air quality, other key reactive gases like carbon monoxide, and the methane lifetime. Overall, E3SM‐chem maintains the climate fidelity of the baseline model while adding at most 20% to the computational cost of the atmosphere model. Hence, interactive chemistry can be a default configuration for long climate simulations at resolutions of 1° or finer, which is crucial for producing self‐consistent chemistry‐climate feedbacks that alter the climate system.

54 ENVIRONMENTAL SCIENCES↗

Significant impact of forcing uncertainty in a large ensemble of climate model simulations

Forcing due to solar and volcanic variability, on the natural side, and greenhouse gas and aerosol emissions, on the anthropogenic side, are the main inputs to climate models. Reliable climate model simulations of past and future climate change depend crucially upon them. Here we analyze large ensembles of simulations using a comprehensive Earth System Model to quantify uncertainties in global climate change attributable to differences in prescribed forcings. The different forcings considered here are those used in the two most recent phases of the Coupled Model Intercomparison Project (CMIP), namely CMIP5 and CMIP6. We show significant differences in simulated global surface air temperature due to volcanic aerosol forcing in the second half of the 19th century and in the early 21st century. The latter arise from small-to-moderate eruptions incorporated in CMIP6 simulations but not in CMIP5 simulations. We also find significant differences in global surface air temperature and Arctic sea ice area due to anthropogenic aerosol forcing in the second half of the 20th century and early 21st century. These differences are as large as those obtained in different versions of an Earth System Model employing identical forcings. In simulations from 2015 to 2100, we find significant differences in the rates of projected global warming arising from CMIP5 and CMIP6 concentration pathways that differ slightly but are equivalent in terms of their nominal radiative forcing levels in 2100. Our results highlight the influence of assumptions about natural and anthropogenic aerosol loadings on carbon budgets, the likelihood of meeting Paris targets, and the equivalence of future forcing scenarios.

54 ENVIRONMENTAL SCIENCES↗

Analysis of ducted fuel injection at high-pressure transcritical conditions using large-eddy simulations

Ducted fuel injection (DFI) is a proposed fuel injection concept for achieving substantial reductions in emissions. In this concept, the fuel is injected through a coannular duct, resulting in increased fuel-air mixing and minimized formation of soot and other unwanted combustion products. Apart from comprehensive experimental investigations on DFI, so far computational studies have been limited to single-point Reynolds-averaged Navier Stokes simulations. Therefore, the objective of this work is to complement these studies by performing large-eddy simulations using a diffuse-interface method to examine the physical mechanisms and combustion processes of DFI, specifically focusing on the mixing process and the effect of fuel-ducting on combustion and pollutant emissions. To this end, finite-rate chemistry simulations are performed of the DFI configuration corresponding to the Engine Combustion Network Spray A injector at transcritical conditions (n-dodecane fuel, 60 bar pressure and 1000 K temperature chamber conditions). A two-equation soot model is employed for the qualitative analysis of soot emissions. Direct comparisons of averaged and instantaneous flow field results with the Spray A configuration are performed to assess the effect of DFI on the first- and second-stage ignition and soot formation. Compared to the free-spray condition, the results show that the DFI case exhibits a combination of (i) increased mass flow rate and entrained air, (ii) larger pressure drop magnitude and flow velocity, and (iii) a closer-to-stoichiometric mixture composition (both globally and locally), each of which is conjectured to contribute toward reduced soot production.

33 ADVANCED PROPULSION SYSTEMS↗

r -process nucleosynthesis and kilonovae from hypermassive neutron star post-merger remnants

ABSTRACT We investigate r-process nucleosynthesis and kilonova emission resulting from binary neutron star (BNS) mergers based on a three-dimensional (3D) general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. The simulation includes a microphysical finite-temperature equation of state (EOS) and neutrino emission and absorption effects via a leakage scheme. We track the thermodynamic properties of the ejecta using Lagrangian tracer particles and determine its composition using the nuclear reaction network SkyNet. We investigate the impact of neutrinos on the nucleosynthetic yields by varying the neutrino luminosities during post-processing. The ejecta show a broad distribution with respect to their electron fraction Ye, peaking between ∼0.25–0.4 depending on the neutrino luminosity employed. We find that the resulting r-process abundance patterns differ from solar, with no significant production of material beyond the second r-process peak when using luminosities recorded by the tracer particles. We also map the HMNS outflows to the radiation hydrodynamics code SNEC and predict the evolution of the bolometric luminosity as well as broadband light curves of the kilonova. The bolometric light curve peaks on the timescale of a day and the brightest emission is seen in the infrared bands. This is the first direct calculation of the r-process yields and kilonova signal expected from HMNS winds based on 3D GRMHD simulations. For longer-lived remnants, these winds may be the dominant ejecta component producing the kilonova emission.

79 ASTRONOMY AND ASTROPHYSICS↗