ABSTRACT
We study the time evolution of molecular clouds across three Milky Way-like isolated disc galaxy simulations at a temporal resolution of 1 Myr and at a range of spatial resolutions spanning ...two orders of magnitude in spatial scale from ∼10 pc up to ∼1 kpc. The cloud evolution networks generated at the highest spatial resolution contain a cumulative total of ∼80 000 separate molecular clouds in different galactic–dynamical environments. We find that clouds undergo mergers at a rate proportional to the crossing time between their centroids, but that their physical properties are largely insensitive to these interactions. Below the gas–disc scale height, the cloud lifetime τlife obeys a scaling relation of the form τlife∝ℓ−0.3 with the cloud size ℓ, consistent with over-densities that collapse, form stars, and are dispersed by stellar feedback. Above the disc scale height, these self-gravitating regions are no longer resolved, so the scaling relation flattens to a constant value of ∼13 Myr, consistent with the turbulent crossing time of the gas disc, as observed in nearby disc galaxies.
Long gamma-ray bursts (GRBs), among the most energetic events in the Universe, are explosions of massive and short-lived stars, so they pinpoint locations of recent star formation. However, several ...GRB host galaxies have recently been found to be deficient in molecular gas (H2), believed to be the fuel of star formation. Moreover, optical spectroscopy of GRB afterglows implies that the molecular phase constitutes only a small fraction of the gas along the GRB line of sight. Here we report the first ever 21 cm line observations of GRB host galaxies, using the AustraliaTelescope Compact Array, implying high levels of atomic hydrogen (H i), which suggests that the connection between atomic gas and star formation is stronger than previously thought. In this case, it is possible that star formation is directly fuelled by atomic gas (or that the H i-to-H2 conversion is very efficient, which rapidly exhaust molecular gas), as has been theoretically shown to be possible. This can happen in low-metallicity gas near the onset of star formation because cooling of gas (necessary for star formation) is faster than the H i-to-H2 conversion. Indeed, large atomic gas reservoirs, together with low molecular gas masses, stellar, and dust masses are consistent with GRB hosts being preferentially galaxies which have very recently started a star formation episode after accreting metal-poor gas from the intergalactic medium. This provides a natural route for forming GRBs in low-metallicity environments. The gas inflow scenario is also consistent with the existence of the companion H I object with no optical counterpart ~19 kpc from the GRB 060505 host, and with the fact that the H I centroids of the GRB 980425 and 060505 hosts do not coincide with optical centres of these galaxies, but are located close to the GRB positions.
We simulate the early stages of the evolution of turbulent, virialized, high-mass protostellar cores, with primary attention to how cores fragment and whether they form a small or large number of ...protostars. Our simulations use the Orion adaptive mesh refinement code to follow the collapse from 60.1 pc scales to 610 AU scales, for durations that cover the main fragmentation phase, using three-dimensional gravito-radiation hydrodynamics. We find that for a wide range of initial conditions radiation feedback from accreting protostars inhibits the formation of fragments, so that the vast majority of the collapsed mass accretes onto one or a few objects. Most of the fragmentation that does occur takes place in massive, self-shielding disks. These are driven to gravitational instability by rapid accretion, producing rapid mass and angular momentum transport that allows most of the gas to accrete onto the central star rather than forming fragments. In contrast, a control run using the same initial conditions but an isothermal equation of state produces much more fragmentation, both in and out of the disk. We conclude that massive cores with observed properties are not likely to fragment into many stars, so that, at least at high masses, the core mass function probably determines the stellar initial mass function. Our results also demonstrate that simulations of massive star-forming regions that do not include radiative transfer, and instead rely on a barotropic equation of state or optically thin heating and cooling curves, are likely to produce misleading results.
Abstract
Carbon monoxide (CO) provides crucial information about the molecular gas properties of galaxies. While 12CO has been targeted extensively, isotopologues such as 13CO have the advantage of ...being less optically thick and observations have recently become accessible across full galaxy discs. We present a comprehensive new data set of 13CO(1–0) observations with the IRAM 30-m telescope of the full discs of nine nearby spiral galaxies from the EMPIRE survey at a spatial resolution of ∼1.5 kpc. 13CO(1–0) is mapped out to 0.7 − 1 r25 and detected at high signal-to-noise ratio throughout our maps. We analyse the 12CO(1–0)-to-13CO(1–0) ratio (ℜ) as a function of galactocentric radius and other parameters such as the 12CO(2–1)-to-12CO(1–0) intensity ratio, the 70-to-160 μm flux density ratio, the star formation rate surface density, the star formation efficiency, and the CO-to-H2 conversion factor. We find that ℜ varies by a factor of 2 at most within and amongst galaxies, with a median value of 11 and larger variations in the galaxy centres than in the discs. We argue that optical depth effects, most likely due to changes in the mixture of diffuse/dense gas, are favoured explanations for the observed ℜ variations, while abundance changes may also be at play. We calculate a spatially resolved 13CO(1–0)-to-H2 conversion factor and find an average value of 1.0 × 1021 cm−2 (K km s−1)−1 over our sample with a standard deviation of a factor of 2. We find that 13CO(1–0) does not appear to be a good predictor of the bulk molecular gas mass in normal galaxy discs due to the presence of a large diffuse phase, but it may be a better tracer of the mass than 12CO(1–0) in the galaxy centres where the fraction of dense gas is larger.
We present a study of the dust-to-gas ratios in five nearby galaxies: NGC 628 (M74), NGC 6503, NGC 7793, UGC 5139 (Holmberg I), and UGC 4305 (Holmberg II). Using Hubble Space Telescope broadband ...WFC3/UVIS UV and optical images from the Treasury program Legacy ExtraGalactic UV Survey (LEGUS) combined with archival HST/Advanced Camera for Surveys data, we correct thousands of individual stars for extinction across these five galaxies using an isochrone-matching (reddening-free Q) method. We generate extinction maps for each galaxy from the individual stellar extinctions using both adaptive and fixed resolution techniques and correlate these maps with neutral H i and CO gas maps from the literature, including the H i Nearby Galaxy Survey and the HERA CO-Line Extragalactic Survey. We calculate dust-to-gas ratios and investigate variations in the dust-to-gas ratio with galaxy metallicity. We find a power-law relationship between dust-to-gas ratio and metallicity, consistent with other studies of dust-to-gas ratio compared to metallicity. We find a change in the relation when H2 is not included. This implies that underestimation of N H 2 in low-metallicity dwarfs from a too-low CO-to-H2 conversion factor XCO could have produced too low a slope in the derived relationship between dust-to-gas ratio and metallicity. We also compare our extinctions to those derived from fitting the spectral energy distribution (SED) using the Bayesian Extinction and Stellar Tool for NGC 7793 and find systematically lower extinctions from SED fitting as compared to isochrone matching.
Abstract
Turbulence is a key process in many fields of astrophysics. Advances in numerical simulations of fluids over the last several decades have revolutionized our understanding of turbulence and ...related processes such as star formation and cosmic ray propagation. However, data from numerical simulations of astrophysical turbulence are often not made public. We introduce a new simulation-oriented database for the astronomical community: the Catalogue for Astrophysical Turbulence Simulations (CATS), located at
www.mhdturbulence.com
. CATS includes magnetohydrodynamic (MHD) turbulent box simulation data products generated by the public codes
athena++
,
arepo
,
enzo,
and
flash
. CATS also includes several synthetic observational data sets, such as turbulent HI data cubes. We also include measured power spectra and three-point correlation functions from some of these data. We discuss the importance of open-source statistical and visualization tools for the analysis of turbulence simulations such as those found in CATS.
Abstract
We present the implementation of a Bayesian formalism within the Stochastically Lighting Up Galaxies (slug) stellar population synthesis code, which is designed to investigate variations in ...the initial mass function (IMF) of star clusters. By comparing observed cluster photometry to large libraries of clusters simulated with a continuously varying IMF, our formalism yields the posterior probability distribution function (PDF) of the cluster mass, age and extinction, jointly with the parameters describing the IMF. We apply this formalism to a sample of star clusters from the nearby galaxy NGC 628, for which broad-band photometry in five filters is available as part of the Legacy ExtraGalactic UV Survey (LEGUS). After allowing the upper-end slope of the IMF (α3) to vary, we recover PDFs for the mass, age and extinction that are broadly consistent with what is found when assuming an invariant Kroupa IMF. However, the posterior PDF for α3 is very broad due to a strong degeneracy with the cluster mass, and it is found to be sensitive to the choice of priors, particularly on the cluster mass. We find only a modest improvement in the constraining power of α3 when adding Hα photometry from the companion Hα-LEGUS survey. Conversely, Hα photometry significantly improves the age determination, reducing the frequency of multi-modal PDFs. With the aid of mock clusters, we quantify the degeneracy between physical parameters, showing how constraints on the cluster mass that are independent of photometry can be used to pin down the IMF properties of star clusters.
Summary Background During in-hospital cardiac arrests, how long resuscitation attempts should be continued before termination of efforts is unknown. We investigated whether duration of resuscitation ...attempts varies between hospitals and whether patients at hospitals that attempt resuscitation for longer have higher survival rates than do those at hospitals with shorter durations of resuscitation efforts. Methods Between 2000 and 2008, we identified 64 339 patients with cardiac arrests at 435 US hospitals within the Get With The Guidelines—Resuscitation registry. For each hospital, we calculated the median duration of resuscitation before termination of efforts in non-survivors as a measure of the hospital's overall tendency for longer attempts. We used multilevel regression models to assess the association between the length of resuscitation attempts and risk-adjusted survival. Our primary endpoints were immediate survival with return of spontaneous circulation during cardiac arrest and survival to hospital discharge. Findings 31 198 of 64 339 (48·5%) patients achieved return of spontaneous circulation and 9912 (15·4%) survived to discharge. For patients achieving return of spontaneous circulation, the median duration of resuscitation was 12 min (IQR 6–21) compared with 20 min (14–30) for non-survivors. Compared with patients at hospitals in the quartile with the shortest median resuscitation attempts in non-survivors (16 min IQR 15–17), those at hospitals in the quartile with the longest attempts (25 min 25–28) had a higher likelihood of return of spontaneous circulation (adjusted risk ratio 1·12, 95% CI 1·06–1·18; p<0·0001) and survival to discharge (1·12, 1·02–1·23; 0·021). Interpretation Duration of resuscitation attempts varies between hospitals. Although we cannot define an optimum duration for resuscitation attempts on the basis of these observational data, our findings suggest that efforts to systematically increase the duration of resuscitation could improve survival in this high-risk population. Funding American Heart Association, Robert Wood Johnson Foundation Clinical Scholars Program, and the National Institutes of Health.
We provide a model for how Kennicutt-Schmidt (KS) laws, which describe the correlation between star formation rate and gas surface or volume density, depend on the molecular line chosen to trace the ...gas. We show that, for lines that can be excited at low temperatures, the KS law depends on how the line critical density compares to the median density In a galaxy's star-forming molecular clouds. High critical density lines trace regions with similar physical properties across galaxy types, and this produces a linear correlation between line luminosity and star formation rate. Low critical density lines probe regions whose properties vary across galaxies, leading to a star formation rate that varies superlinearly with line luminosity. We show that a simple model In which molecular clouds are treated as isothermal and homogenous can quantitatively reproduce the observed correlations between galactic luminosities in far-infrared and in the CO(1 arrow right 0) and HCN(1 arrow right 0) lines, and naturally explains why these correlations have different slopes. We predict that IR similar to line luminosity correlations should change slope for galaxies in which the median density is close to the line critical density. This prediction may be tested by observations of lines such as HCO super(+) (1 arrow right 0) with intermediate critical densities, or by HCN (1 arrow right 0) observations of intensely star-forming high-redshift galaxies with very high densities. Recent observations by Gao et al. hint at just such a change in slope. We argue that deviations from linearity in the HCN (1 arrow right 0)-IR correlation at high luminosity are consistent with the assumption of a constant star formation efficiency.
We report the results of a series of adaptive mesh refinement radiation-hydrodynamic simulations of the collapse of massive star-forming clouds using the ORION code. These simulations are the first ...to include the feedback effects protostellar outflows, as well as protostellar radiative heating and radiation pressure exerted on the infalling, dusty gas. We find that outflows evacuate polar cavities of reduced optical depth through the ambient core. These enhance the radiative flux in the poleward direction so that it is 1.7-15 times larger than that in the midplane. As a result the radiative heating and outward radiation force exerted on the protostellar disk and infalling cloud gas in the equatorial direction are greatly diminished. This simultaneously reduces the Eddington radiation pressure barrier to high-mass star formation and increases the minimum threshold surface density for radiative heating to suppress fragmentation compared to models that do not include outflows. The strength of both these effects depends on the initial core surface density. Lower surface density cores have longer free-fall times and thus massive stars formed within them undergo more Kelvin contraction as the core collapses, leading to more powerful outflows. Furthermore, in lower surface density clouds the ratio of the time required for the outflow to break out of the core to the core free-fall time is smaller, so that these clouds are consequently influenced by outflows at earlier stages of the collapse. As a result, outflow effects are strongest in low surface density cores and weakest in high surface density ones. We also find that radiation focusing in the direction of outflow cavities is sufficient to prevent the formation of radiation pressure-supported circumstellar gas bubbles, in contrast to models which neglect protostellar outflow feedback.