Engineering topics

Naab, Thorsten

Publications and source records attributed to Naab, Thorsten.

10 records

DOE OSTI2023

How do supernova remnants cool? – I. Morphology, optical emission lines, and shocks

ABSTRACT Supernovae (SNe) inject ∼1051 erg in the interstellar medium, thereby shocking and heating the gas. A substantial fraction of this energy is later lost via radiative cooling. We present a post-processing module for the flash code to calculate the cooling radiation from shock-heated gas using collisional excitation data from mappings v. When applying this tool to a simulated SN remnant (SNR), we find that most energy is emitted in the EUV. However, optical emission lines ([O iii], [N ii], [S ii], H α, H β) are usually best observable. Our shock detection scheme shows that [S ii] and [N ii] emissions arise from the thin shell surrounding the SNR, while [O iii], H $\rm \alpha$, and H $\rm \beta$ originate from the volume-filling hot gas inside the SNR bubble. We find that the optical emission lines are affected by the SNR’s complex structure and its projection on to the plane of the sky because the escaping line luminosity can be reduced by 10–80 per cent due to absorption along the line of sight. Additionally, the subtraction of contaminating background radiation is required for the correct classification of an SNR on the oxygen or sulphur BPT diagrams. The electron temperature and density obtained from our synthetic observations match well with the simulation but are very sensitive to the assumed metallicity.

Astronomy & Astrophysics↗

DOE OSTI2023

SILCC – VII. Gas kinematics and multiphase outflows of the simulated ISM at high gas surface densities

ABSTRACT We present magnetohydrodynamic (MHD) simulations of the star-forming multiphase interstellar medium (ISM) in stratified galactic patches with gas surface densities Σgas = 10, 30, 50, and 100 $\mathrm{M_\odot \, pc^{-2}}$. The silcc project simulation framework accounts for non-equilibrium thermal and chemical processes in the warm and cold ISM. The sink-based star formation and feedback model includes stellar winds, hydrogen-ionizing UV radiation, core-collapse supernovae, and cosmic ray (CR) injection and diffusion. The simulations follow the observed relation between Σgas and the star formation rate surface density ΣSFR. CRs qualitatively change the outflow phase structure. Without CRs, the outflows transition from a two-phase (warm and hot at 1 kpc) to a single-phase (hot at 2 kpc) structure. With CRs, the outflow always has three phases (cold, warm, and hot), dominated in mass by the warm phase. The impact of CRs on mass loading decreases for higher Σgas and the mass loading factors of the CR-supported outflows are of order unity independent of ΣSFR. Similar to observations, vertical velocity dispersions of the warm ionized medium (WIM) and the cold neutral medium (CNM) correlate with the star formation rate as $\sigma _\mathrm{z} \propto \Sigma _\mathrm{SFR}^a$, with a ∼ 0.20. In the absence of stellar feedback, we find no correlation. The velocity dispersion of the WIM is a factor ∼2.2 higher than that of the CNM, in agreement with local observations. For $\Sigma _\mathrm{SFR} \gtrsim 1.5 \times 10^{-2}\, \mathrm{M}_\odot \, \mathrm{yr}^{-1}\, \mathrm{kpc}^{-2}$ the WIM motions become supersonic.

Astronomy & Astrophysics↗

DOE OSTI2017

Falling Outer Rotation Curves of Star-forming Galaxies at 0.6 ≲ z ≲ 2.6 Probed with KMOS{sup 3D} and SINS/zC-SINF

We exploit the deep, resolved, Hα kinematic data from the KMOS{sup 3D} and SINS/zC-SINF surveys to examine the largely unexplored outer-disk kinematics of star-forming galaxies (SFGs), out to the peak of cosmic star formation. Our sample contains 101 SFGs, representative of the more massive (9.3≲logM{sub ∗}/M{sub ⊙}≲11.5) main sequence population at 0.6 ≤ z ≤ 2.6. Through a novel stacking approach, we are able to constrain a representative rotation curve extending out to ∼4 effective radii. This average rotation curve exhibits a significant drop in rotation velocity beyond the turnover, with a slope of ΔV/ΔR=−0.26{sub −0.09}{sup +0.10} in units of normalized coordinates V/V {sub max} and R/R {sub turn}. This result confirms that the fall-off seen in some individual galaxies is a common feature of our sample of high-z disks. The outer fall-off strikingly deviates from the flat or mildly rising rotation curves of local spiral galaxies that have similar masses. Through a comparison with models that include baryons and dark matter, we demonstrate that the falling stacked rotation curve is consistent with a high mass fraction of baryons, relative to the total dark matter halo (m {sub d} ≳ 0.05), in combination with a sizeable level of pressure support in the outer disk. These findings agree with recent studies demonstrating that high-z star-forming disks are strongly baryon-dominated within the disk scale, and furthermore suggest that pressure gradients caused by large, turbulent gas motions are present even in their outer disks. These results are largely independent of our model assumptions, such as the presence of stellar bulges, the effect of adiabatic contraction, and variations in halo concentration.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Post-Newtonian Dynamical Modeling of Supermassive Black Holes in Galactic-scale Simulations

We present KETJU, a new extension of the widely used smoothed particle hydrodynamics simulation code GADGET-3. The key feature of the code is the inclusion of algorithmically regularized regions around every supermassive black hole (SMBH). This allows for simultaneously following global galactic-scale dynamical and astrophysical processes, while solving the dynamics of SMBHs, SMBH binaries, and surrounding stellar systems at subparsec scales. The KETJU code includes post-Newtonian terms in the equations of motions of the SMBHs, which enables a new SMBH merger criterion based on the gravitational wave coalescence timescale, pushing the merger separation of SMBHs down to ∼0.005 pc. We test the performance of our code by comparison to NBODY7 and rVINE. We set up dynamically stable multicomponent merger progenitor galaxies to study the SMBH binary evolution during galaxy mergers. In our simulation sample the SMBH binaries do not suffer from the final-parsec problem, which we attribute to the nonspherical shape of the merger remnants. For bulge-only models, the hardening rate decreases with increasing resolution, whereas for models that in addition include massive dark matter halos, the SMBH binary hardening rate becomes practically independent of the mass resolution of the stellar bulge. The SMBHs coalesce on average 200 Myr after the formation of the SMBH binary. However, small differences in the initial SMBH binary eccentricities can result in large differences in the SMBH coalescence times. Finally, we discuss the future prospects of KETJU, which allows for a straightforward inclusion of gas physics in the simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Modeling for Stellar Feedback in Galaxy Formation Simulations

Various heuristic approaches to model unresolved supernova (SN) feedback in galaxy formation simulations exist to reproduce the formation of spiral galaxies and the overall inefficient conversion of gas into stars. Some models, however, require resolution-dependent scalings. We present a subresolution model representing the three major phases of supernova blast wave evolution—free expansion, energy-conserving Sedov–Taylor, and momentum-conserving snowplow—with energy scalings adopted from high-resolution interstellar-medium simulations in both uniform and multiphase media. We allow for the effects of significantly enhanced SN remnant propagation in a multiphase medium with the cooling radius scaling with the hot volume fraction, f{sub hot}, as (1−f{sub hot}){sup −4/5}. We also include winds from young massive stars and AGB stars, Strömgren sphere gas heating by massive stars, and a mechanism that limits gas cooling that is driven by radiative recombination of dense H ii regions. We present initial tests for isolated Milky Way-like systems simulated with the Gadget-based code SPHgal with improved SPH prescription. Compared to pure thermal SN input, the model significantly suppresses star formation at early epochs, with star formation extended both in time and space in better accord with observations. Compared to models with pure thermal SN feedback, the age at which half the stellar mass is assembled increases by a factor of 2.4, and the mass-loading parameter and gas outflow rate from the galactic disk increase by a factor of 2. Simulation results are converged for a variation of two orders of magnitude in particle mass in the range (1.3–130) × 10{sup 4} solar masses.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

KMOS3D: dynamical constraints on the mass budget in early star-forming disks

We exploit deep integral-field spectroscopic observations with KMOS/Very Large Telescope of 240 star-forming disks at 0.6 2. We conclude that baryons make up most of the mass within the disk regions of high-redshift star-forming disk galaxies, with typical disks at z>2 being strongly baryon-dominated within R {sub e}. Substantial object-to-object variations in both stellar and baryonic mass fractions are observed among the galaxies in our sample, larger than what can be accounted for by the formal uncertainties in their respective measurements. In both cases, the mass fractions correlate most strongly with measures of surface density. High-Σ{sub star} galaxies feature stellar mass fractions closer to unity, and systems with high inferred gas or baryonic surface densities leave less room for additional mass components other than stars and molecular gas. Our findings can be interpreted as more extended disks probing further (and more compact disks probing less far) into the dark matter halos that host them.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

The evolution of metallicity and metallicity gradients from z = 2.7 to 0.6 with KMOS{sup 3D}

We present measurements of the [N ii]/Hα ratio as a probe of gas-phase oxygen abundance for a sample of 419 star-forming galaxies at z = 0.6–2.7 from the KMOS{sup 3D} near-IR multi-integral field unit (IFU) survey. The mass–metallicity relation (MZR) is determined consistently with the same sample selection, metallicity tracer, and methodology over the wide redshift range probed by the survey. We find good agreement with long-slit surveys in the literature, except for the low-mass slope of the relation at z∼2.3, where this sample is less biased than previous samples based on optical spectroscopic redshifts. In this regime we measure a steeper slope than some literature results. Excluding the contribution from active galactic nuclei from the MZR reduces sensitivity at the high-mass end, but produces otherwise consistent results. There is no significant dependence of the [N ii]/Hα ratio on star formation rate at fixed redshift and stellar mass. The IFU data allow spatially resolved measurements of [N ii]/Hα, from which we can infer abundance gradients for 180 galaxies, thus tripling the current sample in the literature. The observed gradients are on average flat, with only 15 gradients statistically offset from zero at >3σ. We have modeled the effect of beam smearing, assuming a smooth intrinsic radial gradient and known seeing, inclination, and effective radius for each galaxy. Our seeing-limited observations can recover up to 70% of the intrinsic gradient for the largest, face-on disks, but only 30% for the smaller, more inclined galaxies. We do not find significant trends between observed or corrected gradients and any stellar population, dynamical, or structural galaxy parameters, mostly in agreement with existing studies with much smaller sample sizes. In cosmological simulations, strong feedback is generally required to produce flat gradients at high redshift.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Supernova feedback and the hot gas filling fraction of the interstellar medium

Supernovae (SNe), the most energetic stellar feedback mechanism, are crucial for regulating the interstellar medium (ISM) and launching galactic winds. We explore how supernova remnants (SNRs) create a multiphase medium by performing three-dimentional hydrodynamical simulations at various SN rates, S, and ISM average densities, n-bar . The evolution of an SNR in a self-consistently generated three-phase ISM is qualitatively different from that in a uniform or a two-phase warm/cold medium. By traveling faster and further in the low-density hot phase, the domain of an SNR increases by >10{sup 2.5}. Varying n-bar and S, we find that a steady state can only be achieved when the hot gas volume fraction f{sub V,hot}≲0.6±0.1. Above that level, overlapping SNRs render connecting topology of the hot gas, and the ISM is subjected to thermal runaway. Photoelectric heating (PEH) has a surprisingly strong impact on f{sub V,hot}. For n-bar ≳3 cm{sup −3}, a reasonable PEH rate is able to suppress the thermal runaway. Overall, we determine the critical SN rate for the onset of thermal runaway to be S{sub crit}=200( n-bar /1 cm{sup −3}){sup k}(E{sub SN}/10{sup 51} erg){sup −1} kpc{sup −3} Myr{sup −1}, where k = (1.2, 2.7) for n-bar ⩽1 and >1 cm{sup −3}, respectively. We present a fitting formula of the ISM pressure P( n-bar ,S), which can be used as an effective equation of state in cosmological simulations. Despite the five orders of magnitude span of ( n-bar ,S), the average Mach number varies little: M ≈ 0.5 ± 0.2, 1.2 ± 0.3, and 2.3 ± 0.9 for the hot, warm, and cold phases, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

IMPACT OF SUPERNOVA AND COSMIC-RAY DRIVING ON THE SURFACE BRIGHTNESS OF THE GALACTIC HALO IN SOFT X-RAYS

The halo of the Milky Way contains a hot plasma with a surface brightness in soft X-rays of the order 10{sup −12} erg cm{sup −2} s{sup −1} deg{sup −2}. The origin of this gas is unclear, but so far numerical models of galactic star formation have failed to reproduce such a large surface brightness by several orders of magnitude. In this paper, we analyze simulations of the turbulent, magnetized, multi-phase interstellar medium including thermal feedback by supernova explosions as well as cosmic-ray feedback. We include a time-dependent chemical network, self-shielding by gas and dust, and self-gravity. Pure thermal feedback alone is sufficient to produce the observed surface brightness, although it is very sensitive to the supernova rate. Cosmic rays suppress this sensitivity and reduce the surface brightness because they drive cooler outflows. Self-gravity has by far the largest effect because it accumulates the diffuse gas in the disk in dense clumps and filaments, so that supernovae exploding in voids can eject a large amount of hot gas into the halo. This can boost the surface brightness by several orders of magnitude. Although our simulations do not reach a steady state, all simulations produce surface brightness values of the same order of magnitude as the observations, with the exact value depending sensitively on the simulation parameters. We conclude that star formation feedback alone is sufficient to explain the origin of the hot halo gas, but measurements of the surface brightness alone do not provide useful diagnostics for the study of galactic star formation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

CONNECTION BETWEEN DYNAMICALLY DERIVED INITIAL MASS FUNCTION NORMALIZATION AND STELLAR POPULATION PARAMETERS

We report on empirical trends between the dynamically determined stellar initial mass function (IMF) and stellar population properties for a complete, volume-limited sample of 260 early-type galaxies from the ATLAS{sup 3D} project. We study trends between our dynamically derived IMF normalization α{sub dyn} ≡ (M/L){sub stars}/(M/L){sub Salp} and absorption line strengths, and interpret these via single stellar population-equivalent ages, abundance ratios (measured as [α/Fe]), and total metallicity, [Z/H]. We find that old and alpha-enhanced galaxies tend to have on average heavier (Salpeter-like) mass normalization of the IMF, but stellar population does not appear to be a good predictor of the IMF, with a large range of α{sub dyn} at a given population parameter. As a result, we find weak α{sub dyn}-[α/Fe] and α{sub dyn} –Age correlations and no significant α{sub dyn} –[Z/H] correlation. The observed trends appear significantly weaker than those reported in studies that measure the IMF normalization via the low-mass star demographics inferred through stellar spectral analysis.

79 ASTRONOMY AND ASTROPHYSICS↗