Search NASA⌕ Search

Engineering topics

Klessen, Ralf S.

Publications and source records attributed to Klessen, Ralf S..

Metal Mixing in the r -process Enhanced Ultrafaint Dwarf Galaxy Reticulum II

The ultrafaint dwarf galaxy Reticulum II was enriched by a single rare and prolific r-process event. The r-process content of Reticulum II thus provides a unique opportunity to study metal mixing in a relic first galaxy. Using multi-object high-resolution spectroscopy with VLT/GIRAFFE and Magellan/M2FS, we identify 32 clear spectroscopic member stars and measure abundances of Mg, Ca, Fe, and Ba where possible. We find $72^{+10}_{-12}$ of the stars are r-process-enhanced, with a mean $\langle$[Ba / H]$\rangle$ = -1.68 ± 0.07 and unresolved intrinsic dispersion σ [Ba/H] <0.20. The homogeneous r-process abundances imply that Ret II’s metals are well mixed by the time the r enhanced stars form, which simulations have shown requires at least 100 Myr of metal mixing in between bursts of star formation to homogenize. This is the first direct evidence of bursty star formation in an ultrafaint dwarf galaxy. The homogeneous dilution prefers a prompt and high-yield r-process site, such as collapsar disk winds or prompt neutron star mergers. We also find evidence from [Ba/H] and [Mg/Ca] that the r-enhanced stars in Ret II formed in the absence of substantial pristine gas accretion, perhaps indicating that ≈70% of Ret II stars formed after reionization.

79 ASTRONOMY AND ASTROPHYSICS↗

The density distribution and physical origins of intermittency in supersonic, highly magnetized turbulence with diverse modes of driving

The probability density function (PDF) of the logarithmic density contrast, s = ln (ρ/ρ 0 ), with gas density ρ and mean density ρ 0 , for hydrodynamical supersonic turbulence is well known to have significant non-Gaussian (intermittent) features that monotonically increase with the turbulent Mach number, M. By studying the mass- and volume-weighted s-PDF for an ensemble of 36 sub-to-trans-Alfvenic mean-field, supersonic, isothermal turbulence simulations with different modes of driving, relevant to molecular gas in the cool interstellar medium, we show that a more intricate picture emerges for the non-Gaussian nature of s. Using four independent measures of the non-Gaussian components, we find hydrodynamical-like structure in the highly magnetized plasma for M≲4. However, for M≳4, the non-Gaussian signatures disappear, leaving approximately Gaussian s-statistics – exactly the opposite of hydrodynamical turbulence in the high-M limit. We also find that the non-Gaussian components of the PDF increase monotonically with more compressive driving modes. To understand the M≲4 non-Gaussian features, we use one-dimensional pencil beams to explore the dynamics along and across the large-scale magnetic field, B 0 . Here we discuss kinetic, density, and magnetic field fluctuations from the pencil beams, and identify physical sources of non-Gaussian components to the PDF as single, strong shocks coupled to fast magnetosonic compressions that form along B 0 . We discuss the Gaussianization of the M≳4 s-fields through the lens of two phenomenologies: the self-similarity of the s-field and homogenization of the dynamical time-scales between the over- and underdense regions in the compressible gas.

79 ASTRONOMY AND ASTROPHYSICS↗

Measuring Young Stars in Space and Time. I. The Photometric Catalog and Extinction Properties of N44

In order to better understand the role of high-mass stellar feedback in regulating star formation in giant molecular clouds, we carried out a Hubble Space Telescope (HST) Treasury Program Measuring Young Stars in Space and Time (MYSST) targeting the star-forming complex N44 in the Large Magellanic Cloud (LMC). Using the F555W and F814W broadband filters of both the ACS and WFC3/UVIS, we built a photometric catalog of 461,684 stars down to m {sub F555W} ≃ 29 mag and m {sub F814W} ≃ 28 mag, corresponding to the magnitude of an unreddened 1 Myr pre-main-sequence star of ≈ 0.09 M {sub ☉} at the LMC distance. In this first paper we describe the observing strategy of MYSST and the data reduction procedure and present the photometric catalog. We identify multiple young stellar populations tracing the gaseous rim of N44's superbubble, together with various contaminants belonging to the LMC field population. We also determine the reddening properties from the slope of the elongated red clump (RC) feature by applying the machine-learning algorithm RANSAC, and we select a set of upper-main-sequence stars as primary probes to build an extinction map, deriving a relatively modest median extinction A {sub F555W} ≃ 0.77 mag. The same procedure applied to the RC provides A {sub F555W} ≃ 0.68 mag.

47 OTHER INSTRUMENTATION↗

Measuring Young Stars in Space and Time. II. The Pre-main-sequence Stellar Content of N44

The Hubble Space Telescope survey Measuring Young Stars in Space and Time (MYSST) entails some of the deepest photometric observations of extragalactic star formation, capturing even the lowest-mass stars of the active star-forming complex N44 in the Large Magellanic Cloud. We employ the new MYSST stellar catalog to identify and characterize the content of young pre-main-sequence (PMS) stars across N44 and analyze the PMS clustering structure. To distinguish PMS stars from more evolved line of sight contaminants, a non-trivial task due to several effects that alter photometry, we utilize a machine-learning classification approach. This consists of training a support vector machine (SVM) and a random forest (RF) on a carefully selected subset of the MYSST data and categorize all observed stars as PMS or non-PMS. Combining SVM and RF predictions to retrieve the most robust set of PMS sources, we find ∼26,700 candidates with a PMS probability above 95% across N44. Employing a clustering approach based on a nearest neighbor surface density estimate, we identify 16 prominent PMS structures at 1σ significance above the mean density with sub-clusters persisting up to and beyond 3σ significance. The most active star-forming center, located at the western edge of N44's bubble, is a subcluster with an effective radius of ∼5.6 pc entailing more than 1100 PMS candidates. Furthermore, we confirm that almost all identified clusters coincide with known H ii regions and are close to or harbor massive young O stars or YSOs previously discovered by MUSE and Spitzer observations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

High-resolution CARMA Observation of Molecular Gas in the North America and Pelican Nebulae

In this paper, we present the first results from a CARMA high-resolution {sup 12}CO(1-0), {sup 13}CO(1-0), and C{sup 18}O(1-0) molecular line survey of the North America and Pelican (NAP) Nebulae. CARMA observations have been combined with single-dish data from the Purple Mountain 13.7 m telescope, to add short spacings and to produce high-dynamic-range images. We find that the molecular gas is predominantly shaped by the W80 H ii bubble, driven by an O star. Several bright rims noted in the observation are probably remnant molecular clouds, heated and stripped by the massive star. Matching these rims in molecular lines and optical images, we construct a model of the three-dimensional structure of the NAP complex. Two groups of molecular clumps/filaments are on the near side of the bubble: one is being pushed toward us, whereas the other is moving toward the bubble. Another group is on the far side of the bubble, and moving away. The young stellar objects in the Gulf region reside in three different clusters, each hosted by a cloud from one of the three molecular clump groups. Although all gas content in the NAP is impacted by feedback from the central O star, some regions show no signs of star formation, while other areas clearly exhibit star formation activity. Additional molecular gas being carved by feedback includes cometary structures in the Pelican Head region, and the boomerang features at the boundary of the Gulf region. The results show that the NAP complex is an ideal place for the study of feedback effects on star formation.

74 ATOMIC AND MOLECULAR PHYSICS↗

A multishock model for the density variance of anisotropic, highly magnetized, supersonic turbulence

ABSTRACT Shocks form the basis of our understanding for the density and velocity statistics of supersonic turbulent flows, such as those found in the cool interstellar medium (ISM). The variance of the density field, $\sigma ^2_{\rho /\rho _0}$, is of particular interest for molecular clouds (MCs), the birthplaces of stars in the Universe. The density variance may be used to infer underlying physical processes in an MC, and parametrizes the star formation (SF) rate of a cloud. However, models for $\sigma ^2_{\rho /\rho _0}$ all share a common feature – the variance is assumed to be isotropic. This assumption does not hold when a trans-/sub-Alfvénic mean magnetic field, ${B}_0$, is present in the cloud, which observations suggest is relevant for some MCs. We develop an anisotropic model for $\sigma _{\rho /\rho _0}^2$, using contributions from hydrodynamical and fast magnetosonic shocks that propagate orthogonal to each other. Our model predicts an upper bound for $\sigma _{\rho /\rho _0}^2$ in the high Mach number $(\mathcal {M})$ limit as small-scale density fluctuations become suppressed by the strong ${B}_0$. The model reduces to the isotropic $\sigma _{\rho /\rho _0}^2\!-\!\mathcal {M}$ relation in the hydrodynamical limit. To validate our model, we calculate $\sigma _{\rho /\rho _0}^2$ from 12 high-resolution, three-dimensional, supersonic, sub-Alfvénic magnetohydrodynamical (MHD) turbulence simulations and find good agreement with our theory. We discuss how the two MHD shocks may be the bimodally oriented overdensities observed in some MCs and the implications for SF theory in the presence of a sub-Alfvénic ${B}_0$. By creating an anisotropic, supersonic density fluctuation model, this study paves the way for SF theory in the highly anisotropic regime of interstellar turbulence.

Beattie, James R.↗

Applying the Tremaine–Weinberg Method to Nearby Galaxies: Stellar-mass-based Pattern Speeds and Comparisons with ISM Kinematics

We apply the Tremaine–Weinberg method to 19 nearby galaxies using stellar mass surface densities and velocities derived from the PHANGS-MUSE survey, to calculate (primarily bar) pattern speeds (Ω{sub P}). After quality checks, we find that around half (10) of these stellar-mass-based measurements are reliable. For those galaxies, we find good agreement between our results and previously published pattern speeds, and we use rotation curves to calculate major resonance locations (corotation radii and Lindblad resonances). We also compare these stellar-mass-derived pattern speeds with Hα (from MUSE) and CO(J = 2 − 1) emission from the PHANGS-ALMA survey. We find that in the case of these clumpy interstellar medium (ISM) tracers, this method erroneously gives a signal that is simply the angular frequency at a representative radius set by the distribution of these clumps (Ω{sub clump}), and that this Ω{sub clump} is significantly different from Ω{sub P} (∼20% in the case of Hα, and ∼50% in the case of CO). Thus, we conclude that it is inadvisable to use “pattern speeds” derived from ISM kinematics. Finally, we compare our derived pattern speeds and corotation radii, along with bar properties, to the global parameters of these galaxies. Consistent with previous studies, we find that galaxies with a later Hubble type have a larger ratio of corotation radius to bar length, more molecular-gas-rich galaxies have higher Ω{sub P}, and more bulge-dominated galaxies have lower Ω{sub P}. Unlike earlier works, however, there are no clear trends between the bar strength and Ω{sub P}, nor between the total stellar mass surface density and the pattern speed.

74 ATOMIC AND MOLECULAR PHYSICS↗

How the First Stars Regulated Star Formation. II. Enrichment by Nearby Supernovae

Metals from Population III (Pop III) supernovae led to the formation of less massive Pop II stars in the early universe, altering the course of evolution of primeval galaxies and cosmological reionization. There are a variety of scenarios in which heavy elements from the first supernovae were taken up into second-generation stars, but cosmological simulations only model them on the largest scales. We present small-scale, high-resolution simulations of the chemical enrichment of a primordial halo by a nearby supernova after partial evaporation by the progenitor star. We find that ejecta from the explosion crash into and mix violently with ablative flows driven off the halo by the star, creating dense, enriched clumps capable of collapsing into Pop II stars. Metals may mix less efficiently with the partially exposed core of the halo, so it might form either Pop III or Pop II stars. Both Pop II and III stars may thus form after the collision if the ejecta do not strip all the gas from the halo. The partial evaporation of the halo prior to the explosion is crucial to its later enrichment by the supernova.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Maximum Mass of Accreting Primordial Supermassive Stars

Supermassive primordial stars are suspected to be the progenitors of the most massive quasars at z ∼ 6. Previous studies of such stars were either unable to resolve hydrodynamical timescales or considered stars in isolation, not in the extreme accretion flows in which they actually form. Therefore, they could not self-consistently predict their final masses at collapse, or those of the resulting supermassive black hole seeds, but rather invoked comparison to simple polytropic models. Here, we systematically examine the birth, evolution, and collapse of accreting, non-rotating supermassive stars under accretion rates of 0.01–10 M {sub ⊙} yr{sup −1} using the stellar evolution code Kepler . Our approach includes post-Newtonian corrections to the stellar structure and an adaptive nuclear network and can transition to following the hydrodynamic evolution of supermassive stars after they encounter the general relativistic instability. We find that this instability triggers the collapse of the star at masses of 150,000–330,000 M {sub ⊙} for accretion rates of 0.1–10 M {sub ⊙} yr{sup −1}, and that the final mass of the star scales roughly logarithmically with the rate. The structure of the star, and thus its stability against collapse, is sensitive to the treatment of convection and the heat content of the outer accreted envelope. Comparison with other codes suggests differences here may lead to small deviations in the evolutionary state of the star as a function of time, that worsen with accretion rate. Since the general relativistic instability leads to the immediate death of these stars, our models place an upper limit on the masses of the first quasars at birth.

79 ASTRONOMY AND ASTROPHYSICS↗

GRAVITATIONAL CONTRACTION VERSUS SUPERNOVA DRIVING AND THE ORIGIN OF THE VELOCITY DISPERSION–SIZE RELATION IN MOLECULAR CLOUDS

Molecular cloud (MC) observations show that clouds have non-thermal velocity dispersions that scale with the cloud size as σ ∝ R {sup 1/2} at a constant surface density, and for varying surface density scale with both the cloud’s size and surface density, σ {sup 2} ∝ R Σ. The energy source driving these chaotic motions remains poorly understood. We describe the velocity dispersions observed in a cloud population formed in a numerical simulation of a magnetized, stratified, supernova (SN)-driven, interstellar medium, including diffuse heating and radiative cooling, before and after we include the effects of the self-gravity of the gas. We compare the relationships between velocity dispersion, size, and surface density measured in the simulated cloud population to those found in observations of Galactic MCs. Our simulations prior to the onset of self-gravity suggest that external SN explosions alone do not drive turbulent motions of the observed magnitudes within dense clouds. On the other hand, self-gravity induces non-thermal motions as gravitationally bound clouds begin to collapse in our model, approaching the observed relations between velocity dispersion, size, and surface density. Energy conservation suggests that the observed behavior is consistent with the kinetic energy being proportional to the gravitational energy. However, the clouds in our model show no sign of reaching a stable equilibrium state at any time, even for strongly magnetized clouds. We conclude that gravitationally bound MCs are always in a state of gravitational contraction and their properties are a natural result of this chaotic collapse. In order to agree with observed star formation efficiencies, this process must be terminated by the early destruction of the clouds, presumably from internal stellar feedback.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

Structure analysis of simulated molecular clouds with the Δ-variance

Here, we employ the Δ-variance analysis and study the turbulent gas dynamics of simulated molecular clouds (MCs). Our models account for a simplified treatment of time-dependent chemistry and the non-isothermal nature of the gas. We investigate simulations using three different initial mean number densities of n 0 = 30, 100 and 300 cm -3 that span the range of values typical for MCs in the solar neighbourhood. Furthermore, we model the CO line emission in a post-processing step using a radiative transfer code. We evaluate Δ-variance spectra for centroid velocity (CV) maps as well as for integrated intensity and column density maps for various chemical components: the total, H 2 and 12 CO number density and the integrated intensity of both the 12 CO and 13 CO (J = 1 → 0) lines. The spectral slopes of the Δ-variance computed on the CV maps for the total and H 2 number density are significantly steeper compared to the different CO tracers. We find slopes for the linewidth–size relation ranging from 0.4 to 0.7 for the total and H 2 density models, while the slopes for the various CO tracers range from 0.2 to 0.4 and underestimate the values for the total and H 2 density by a factor of 1.5–3.0. We demonstrate that optical depth effects can significantly alter the Δ-variance spectra. Furthermore, we report a critical density threshold of 100 cm -3 at which the Δ-variance slopes of the various CO tracers change sign. We thus conclude that carbon monoxide traces the total cloud structure well only if the average cloud density lies above this limit.

79 ASTRONOMY AND ASTROPHYSICS↗

A NEW APPROACH TO DETERMINE OPTICALLY THICK H{sub 2} COOLING AND ITS EFFECT ON PRIMORDIAL STAR FORMATION

We present a new method for estimating the H{sub 2} cooling rate in the optically thick regime in simulations of primordial star formation. Our new approach is based on the TreeCol algorithm, which projects matter distributions onto a spherical grid to create maps of column densities for each fluid element in the computational domain. We have improved this algorithm by using the relative gas velocities to weight the individual matter contributions with the relative spectral line overlaps, in order to properly account for the Doppler effect. We compare our new method to the widely used Sobolev approximation, which yields an estimate for the column density based on the local velocity gradient and the thermal velocity. This approach generally underestimates the photon escape probability because it neglects the density gradient and the actual shape of the cloud. We present a correction factor for the true line overlap in the Sobolev approximation and a new method based on local quantities, which fits the exact results reasonably well during the collapse of the cloud, with the error in the cooling rates always being less than 10%. Analytical fitting formulae fail at determining the photon escape probability after formation of the first protostar (error of ∼40%) because they are based on the assumption of spherical symmetry and therefore break down once a protostellar accretion disk has formed. Our method yields lower temperatures and hence promotes fragmentation for densities above ∼10{sup 10} cm{sup –3} at a distance of ∼200 AU from the first protostar. Since the overall accretion rates are hardly affected by the cooling implementation, we expect Pop III stars to have lower masses in our simulations, compared to the results of previous simulations that used the Sobolev approximation.

79 ASTRONOMY AND ASTROPHYSICS↗

Collective outflow from a small multiple stellar system

The formation of high-mass stars is usually accompanied by powerful protostellar outflows. Such high-mass outflows are not simply scaled-up versions of their lower-mass counterparts, since observations suggest that the collimation degree degrades with stellar mass. Theoretically, the origins of massive outflows remain open to question because radiative feedback and fragmentation of the accretion flow around the most massive stars, with M > 15 M {sub ☉}, may impede the driving of magnetic disk winds. We here present a three-dimensional simulation of the early stages of core fragmentation and massive star formation that includes a subgrid-scale model for protostellar outflows. We find that stars that form in a common accretion flow tend to have aligned outflow axes, so that the individual jets of multiple stars can combine to form a collective outflow. We compare our simulation to observations with synthetic H{sub 2} and CO observations and find that the morphology and kinematics of such a collective outflow resembles some observed massive outflows, such as Cepheus A and DR 21. We finally compare physical quantities derived from simulated observations of our models to the actual values in the models to examine the reliability of standard methods for deriving physical quantities, demonstrating that those methods indeed recover the actual values to within a factor of two to three.

79 ASTRONOMY AND ASTROPHYSICS↗