Source author record

S. C. O. Glover

S. C. O. Glover appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2026arXiv

PHANGS-JWST: the largest extragalactic molecular cloud catalog traced by polycyclic aromatic hydrocarbon emission

High-resolution JWST images of nearby spiral galaxies reveal polycyclic aromatic hydrocarbon (PAH) emission structures that trace molecular gas, including CO-dark regions. We identify ISM cloud structures in PHANGS-JWST 7.7 $μ$m PAH maps for 66 galaxies, smoothed to 30 pc and at native resolution, extracting 108,466 and 146,040 clouds, respectively. Molecular properties were inferred using a linear conversion from PAH to CO. Given the tendency for clouds in galaxy centers to overlap in velocity space, we opted to flag these and omit them from the analysis in this work. The remaining clouds correspond to giant molecular clouds, such as those detected in CO(2-1) emission by ALMA, or lower surface density clouds that either fall below the ALMA detection limits of existing maps or genuinely have no molecular counterpart. Cross-matching with ALMA CO maps at 90 pc in 27 galaxies shows that 41 % of PAH clouds have CO associations. The converted molecular properties vary little across environments, but the most massive clouds are preferentially found in spiral arms. Fitting lognormal mass distributions down to $2\times10^{3} M_{\odot}$ shows that spiral arms host the highest-mass clouds, consistent with enhanced formation in arm gravitational potentials. Cloud molecular surface densities decline by a factor of $\sim 1.5-2$ toward $2 - 3 R_{e}$. However, the trend largely varies in individual galaxies, with flat, decreasing, and even no trend as a function of galactocentric radius. Factors like large-scale processes and morphologies might influence the observed trends. We publish two catalogs online, one at the common resolution of 30 pc and another at the native resolution. We expect them to have broad utility for future PAH clouds, molecular clouds, and star formation studies.

preprint2022arXiv

Linking stellar populations to HII regions across nearby galaxies: I. Constraining pre-supernova feedback from young clusters in NGC1672

One of the fundamental factors regulating the evolution of galaxies is stellar feedback. However, we still do not have strong observational constraints on the relative importance of the different feedback mechanisms (e.g. radiation, ionised gas pressure, stellar winds) in driving HII region evolution and molecular cloud disruption. In this letter, we constrain the relative importance of the various feedback mechanisms from young massive star populations by resolving HII regions across the disk of the nearby star-forming galaxy NGC 1672. We combine measurements of ionised gas nebular lines obtained by PHANGS-MUSE, with high-resolution imaging from the HST in both the narrow-band Hα and broad-band filters. We identify a sample of 40 isolated, compact HII regions in the HST Hα image, for which we measure the sizes that were previously unresolved in seeing-limited ground-based observations. Additionally, we identify the ionisation source(s) for each HII region from catalogues produced as part of the PHANGS-HST survey. We find that the HII regions investigated are mildly dominated by thermal or wind pressure, yet their elevation above the radiation pressure is within the expected uncertainty range. We see that radiation pressure provides a substantially higher contribution to the total pressure than previously found in the literature over similar size scales. In general, we find higher pressures within more compact HII regions, which is driven by the inherent size scaling relations of each pressure term, albeit with significant scatter introduced by the variation in the stellar population properties (e.g. luminosity, mass, age, metallicity). For nearby galaxies, here we provide a promising approach that could yield the statistics required to map out how the importance of different stellar feedback mechanisms evolve over the lifetime of an HII region.

preprint2020arXiv

Cloud formation in the atomic and molecular phase: HI self absorption (HISA) towards a Giant Molecular Filament

Molecular clouds form from the atomic phase of the interstellar medium. However, characterizing the transition between the atomic and the molecular interstellar medium (ISM) is a difficult observational task. Here we address cloud formation processes by combining HSIA with molecular line data. One scenario proposed by numerical simulations is that the column density probability density functions (N-PDF) evolves from a log-normal shape at early times to a power-law-like shape at later times. In this paper, we study the cold atomic component of the giant molecular filament GMF38a (d=3.4 kpc, length$\sim230$ pc). We identify an extended HISA feature, which is partly correlated with the 13CO emission. The peak velocities of the HISA and 13CO observations agree well on the eastern side of the filament, whereas a velocity offset of approximately 4 km s$^{-1}$ is found on the western side. The sonic Mach number we derive from the linewidth measurements shows that a large fraction of the HISA, which is ascribed to the cold neutral medium (CNM), is at subsonic and transonic velocities. The column density of the CNM is on the order of 10$^{20}$ to 10$^{21}$ cm$^{-2}$. The column density of molecular hydrogen is an order of magnitude higher. The N-PDFs from HISA (CNM), HI emission (WNM+CNM), and 13CO (molecular component) are well described by log-normal functions, which is in agreement with turbulent motions being the main driver of cloud dynamics. The N-PDF of the molecular component also shows a power law in the high column-density region, indicating self-gravity. We suggest that we are witnessing two different evolutionary stages within the filament. The eastern subregion seems to be forming a molecular cloud out of the atomic gas, whereas the western subregion already shows high column density peaks, active star formation and evidence of related feedback processes.

preprint2020arXiv

Dynamical cloud formation traced by atomic and molecular gas

Context: Atomic and molecular cloud formation is a dynamical process. However, kinematic signatures of these processes are still observationally poorly constrained. Methods: Targeting the cloud-scale environment of the prototypical infrared dark cloud G28.3, we employ spectral line imaging observations of the two atomic lines HI and [CI] as well as molecular lines observations in 13CO in the 1--0 and 3--2 transitions. The analysis comprises investigations of the kinematic properties of the different tracers, estimates of the mass flow rates, velocity structure functions, a Histogram of Oriented Gradients (HOG) study as well as comparisons to simulations. Results: The central IRDC is embedded in a more diffuse envelope of cold neutral medium (CNM) traced by HI self-absorption (HISA) and molecular gas. The spectral line data as well as the HOG and structure function analysis indicate a possible kinematic decoupling of the HI from the other gas compounds. Spectral analysis and position-velocity diagrams reveal two velocity components that converge at the position of the IRDC. Estimated mass flow rates appear rather constant from the cloud edge toward the center. The velocity structure function analysis is consistent with gas flows being dominated by the formation of hierarchical structures. Conclusions: The observations and analysis are consistent with a picture where the IRDC G28 is formed at the center of two converging gas flows. While the approximately constant mass flow rates are consistent with a self-similar, gravitationally driven collapse of the cloud, external compression by, e.g., spiral arm shocks or supernovae explosions cannot be excluded yet. Future investigations should aim at differentiating the origin of such converging gas flows.

preprint2020arXiv

LEGO II: A 3 mm molecular line study covering 100 pc of one of the most actively star-forming portions within the Milky Way Disc

The current generation of (sub)mm-telescopes has allowed molecular line emission to become a major tool for studying the physical, kinematic, and chemical properties of extragalactic systems, yet exploiting these observations requires a detailed understanding of where emission lines originate within the Milky Way. In this paper, we present 60$^{\prime\prime}$ ($\sim$3pc) resolution observations of many 3mm-band molecular lines across a large map of the W49 massive star-forming region ($\sim$100$\times$100pc at 11kpc), which were taken as part of the "LEGO" IRAM-30m large project. We find that the spatial extent or brightness of the molecular line transitions are not well correlated with their critical densities, highlighting abundance and optical depth must be considered when estimating line emission characteristics. We explore how the total emission and emission efficiency (i.e. line brightness per H$_{2}$ column density) of the line emission vary as a function of molecular hydrogen column density and dust temperature. We find that there is not a single region of this parameter space responsible for the brightest and most efficiently emitting gas for all species. For example, we find that the HCN transition shows high emission efficiency at high column density ($10^{22}$cm$^{-2}$) and moderate temperatures (35K), whilst e.g. N$_2$H$^+$ emits most efficiently towards lower temperatures ($10^{22}$cm$^{-2}$; <20K). We determine $X_{\mathrm{CO} (1-0)} \sim 0.3 \times 10^{20} \mathrm{cm^{-2}(Kkms^{-1})^{-1}}$, and $α_{\mathrm{HCN} (1-0)} \sim 30\mathrm{M_\odot(Kkms^{-1}pc^2)^{-1}}$, which both differ significantly from the commonly adopted values. In all, these results suggest caution should be taken when interpreting molecular line emission.

preprint2020arXiv

Less than the sum of its parts: the dust-corrected H$α$ luminosity of star-forming galaxies explored at different spatial resolutions with MaNGA and MUSE

The H$α$ and H$β$ emission line luminosities measured in a single integrated spectrum are affected in non-trivial ways by point-to-point variations in dust attenuation in a galaxy. This work investigates the impact of this variation when estimating global H$α$ luminosities corrected for the presence of dust by a global Balmer decrement. Analytical arguments show that the dust-corrected H$α$ luminosity is always underestimated when using the global H$α$/H$β$ flux ratio to correct for dust attenuation. We measure this effect on 156 face-on star-forming galaxies from the Mapping Nearby Galaxies at APO (MaNGA) survey. At 1-2 kpc spatial resolution, the effect is small but systematic, with the integrated dust-corrected H$α$ luminosity underestimated by $2$-$4$ per cent (and typically not more than by $10$ per cent), and depends on the specific star formation rate of the galaxy. Given the spatial resolution of MaNGA, these are lower limits for the effect. From Multi Unit Spectroscopic Explorer (MUSE) observations of NGC 628 with a resolution of 36 pc we find the discrepancy between the globally and the point-by-point dust-corrected H$α$ luminosity to be $14 \pm 1$ per cent, which may still underestimate the true effect. We use toy models and simulations to show that the true difference depends strongly on the spatial variance of the H$α$/H$β$ flux ratio, and on the slope of the relation between H$α$ luminosity and dust attenuation within a galaxy. Larger samples of higher spatial resolution observations are required to quantify the dependence of this effect as a function of galaxy properties.

preprint2020arXiv

Ubiquitous velocity fluctuations throughout the molecular interstellar medium

The density structure of the interstellar medium (ISM) determines where stars form and release energy, momentum, and heavy elements, driving galaxy evolution. Density variations are seeded and amplified by gas motion, but the exact nature of this motion is unknown across spatial scale and galactic environment. Although dense star-forming gas likely emerges from a combination of instabilities, convergent flows, and turbulence, establishing the precise origin is challenging because it requires quantifying gas motion over many orders of magnitude in spatial scale. Here we measure the motion of molecular gas in the Milky Way and in nearby galaxy NGC 4321, assembling observations that span an unprecedented spatial dynamic range ($10^{-1}{-}10^3$ pc). We detect ubiquitous velocity fluctuations across all spatial scales and galactic environments. Statistical analysis of these fluctuations indicates how star-forming gas is assembled. We discover oscillatory gas flows with wavelengths ranging from $0.3{-}400$ pc. These flows are coupled to regularly-spaced density enhancements that likely form via gravitational instabilities. We also identify stochastic and scale-free velocity and density fluctuations, consistent with the structure generated in turbulent flows. Our results demonstrate that ISM structure cannot be considered in isolation. Instead, its formation and evolution is controlled by nested, interdependent flows of matter covering many orders of magnitude in spatial scale.

preprint2019arXiv

Feedback in W49A diagnosed with radio recombination lines and models

We present images of radio recombination lines (RRLs) at wavelengths around 17 cm from the star-forming region W49A to determine the kinematics of ionized gas in the THOR survey (The HI/OH/Recombination line survey of the inner Milky Way) at an angular resolution of 16.8"x13.8". The distribution of ionized gas appears to be affected by feedback processes from the star clusters in W49A. The velocity structure of the RRLs shows a complex behavior with respect to the molecular gas. We find a shell-like distribution of ionized gas as traced by RRL emission surrounding the central cluster of OB stars in W49A. We describe the evolution of the shell with the recent feedback model code WARPFIELD that includes the important physical processes and has previously been applied to the 30 Doradus region in the Large Magellanic Cloud. The cloud structure and dynamics of W49A are in agreement with a feedback-driven shell that is re-collapsing. The shell may have triggered star formation in other parts of W49A. We suggest that W49A is a potential candidate for star formation regulated by feedback-driven and re-collapsing shells.

preprint2019arXiv

Mapping metallicity variations across nearby galaxy disks

The distribution of metals within a galaxy traces the baryon cycle and the buildup of galactic disks, but the detailed gas phase metallicity distribution remains poorly sampled. We have determined the gas phase oxygen abundances for 7,138 HII regions across the disks of eight nearby galaxies using VLT/MUSE optical integral field spectroscopy as part of the PHANGS-MUSE survey. After removing the first order radial gradients present in each galaxy, we look at the statistics of the metallicity offset (Delta O/H) and explore azimuthal variations. Across each galaxy, we find low (sigma=0.03-0.05 dex) scatter at any given radius, indicative of efficient mixing. We compare physical parameters for those HII regions that are 1 sigma outliers towards both enhanced and reduced abundances. Regions with enhanced abundances have high ionization parameter, higher Halpha luminosity, lower Halpha velocity dispersion, younger star clusters and associated molecular gas clouds show higher molecular gas densities. This indicates recent star formation has locally enriched the material. Regions with reduced abundances show increased Halpha velocity dispersions, suggestive of mixing introducing more pristine material. We observe subtle azimuthal variations in half of the sample, but can not always cleanly associate this with the spiral pattern. Regions with enhanced and reduced abundances are found distributed throughout the disk, and in half of our galaxies we can identify subsections of spiral arms with clearly associated metallicity gradients. This suggests spiral arms play a role in organizing and mixing the ISM.

preprint2019arXiv

Strong excess Faraday rotation on the Inside of the Sagittarius spiral arm

We present first results for Faraday rotation of compact polarized sources (1 to 2 GHz continuum) in The HI/OH/Recombination line (THOR) survey of the inner Galaxy. In the Galactic longitude range 39 degr < l < 52 degr, we find rotation measures in the range -310 rad/m2 < RM < +4219 rad/m2, with the highest values concentrated within a degree of l = 48 degrees at the Sagittarius arm tangent. Most of the high RMs arise in diffuse plasma, along lines of sight that do not intersect HII regions. For l > 49 degr, RM drops off rapidly, while at l < 47 degr, the mean RM is higher with a larger standard deviation than at l > 49 degr. We attribute the RM structure to the compressed diffuse Warm Ionized Medium in the spiral arm, upstream of the major star formation regions. The Sagittarius arm acts as a significant Faraday screen inside the Galaxy. This has implications for models of the Galactic magnetic field and the expected amount of Faraday rotation of Fast Radio Bursts from their host galaxies. We emphasize the importance of sensitivity to high Faraday depth in future polarization surveys.

preprint2019arXiv

The HI/OH/Recombination line survey of the inner Milky Way (THOR): data release 2 and HI overview

With the $Karl~G.~Jansky$ Very Large Array (VLA) in C-configuration, we observed a large portion of the first Galactic quadrant ($l=14.0-67.4^\circ $ and $\lvert b \rvert \leq 1.25^\circ $) achieving an angular resolution of $\leq 40^{\prime\prime}$. At $L$ Band, the WIDAR correlator at the VLA was set to cover the 21~cm HI line, four OH transitions, a series of H$nα$ radio recombination lines (RRLs; $n=151$ to 186), and eight 128~MHz wide continuum spectral windows (SPWs) simultaneously. The HI emission shows clear filamentary substructures at negative velocities with low velocity crowding. The emission at positive velocities is more smeared-out likely due to higher spatial and velocity crowding of structures at the positive velocities. Comparing to the spiral arm model of the Milky Way, the atomic gas follows the Sagittarius and Perseus Arm well but with significant material in the inter-arm regions. With the C-configuration-only HI+continuum data, we produced a HI optical depth map of the THOR areal coverage from 228 absorption spectra with the nearest-neighbor method. With this $τ$ map, we corrected the HI emission for optical depth and the derived column density is 38% higher than the column density with optically thin assumption. The total HI mass with optical depth correction in the survey region is 4.7$\times10^8~M_\odot$, 31% more than the mass derived assuming the emission is optically thin. If we apply this 31% correction to the whole Milky Way, the total atomic gas mass would be 9.4-10.5$\times 10^9~M_\odot$. Comparing the HI with existing CO data, we find a significant increase in the atomic-to-molecular gas ratio from the spiral arms to the inter-arm regions.

preprint2019arXiv

WARPFIELD Population Synthesis: The physics of (extra-)Galactic star formation and feedback driven cloud structure and emission from sub-to-kpc scales

We present a novel method to model galactic scale star formation and the resulting emission from star clusters and the multi-phase interstellar medium. We combine global parameters, such as SFR and CMF, with {\sc warpfield} which provides a description of the feedback-driven evolution of individual star-forming regions. Our approach includes stellar evolution, stellar winds, radiation pressure, supernovae, all of which couple to the dynamical evolution of the parental cloud in a highly non-linear fashion. The heating of diffuse gas and dust is calculated self-consistently with the age, mass and density dependent escape fractions of the local star-forming regions. From this we construct the interstellar radiation field at any point in the galaxy, and we employ the multi-frequency Monte Carlo radiative transfer code {\sc polaris} to produce synthetic emission maps for the one-to-one comparison with observational data. We demonstrate the capabilities of our approach by applying the method to a Milky Way like galaxy built-up in a high-resolution cosmological MHD simulation. We give three examples. First, we compute the multi-scale distribution of electron $n_{e^-}$ and $T_{e^{-}}$ and synthesize the MW all-sky H$α$ emission. We use a multipole expansion method to show that the resulting maps are consistent with observations. Second, we predict the expected \SIII 9530~Å emission. This line is a key target of several planned large survey programs. It suffers less extinction than other diagnostic lines and provides information about star formation in very dense environments that are otherwise observationally not readily accessible. Third, we explore the effects of differential extinction as seen by an extra-galactic observer, and discuss the consequences for the correct interpretation of \Ha emission as a star-formation rate tracer at different viewing angles.(abridged)

preprint2016arXiv

Continuum sources from the THOR survey between 1 and 2 GHz

We carried out a large program with the Karl G. Jansky Very Large Array (VLA): "THOR: The HI, OH, Recombination line survey of the Milky Way". We observed a significant portion of the Galactic plane in the first quadrant of the Milky Way in the 21cm HI line, 4 OH transitions, 19 radio recombination lines, and continuum from 1 to 2 GHz. In this paper we present a catalog of the continuum sources in the first half of the survey (l=14.0-37.9deg and l=47.1-51.2deg, |b|<1.1deg) at a spatial resolution of 10-25", with a spatially varying noise level of ~0.3-1 mJy/beam. The catalog contains ~4400 sources. Around 1200 of these are spatially resolved, and ~1000 are possible artifacts, given their low signal-to-noise ratios. Since the spatial distribution of the unresolved objects is evenly distributed and not confined to the Galactic plane, most of them are extragalactic. Thanks to the broad bandwidth of the observations from 1 to 2 GHz, we are able to determine a reliable spectral index for ~1800 sources. The spectral index distribution reveals a double-peaked profile with maxima at spectral indices of alpha = -1 and alpha = 0 , corresponding to steep declining and flat spectra, respectively. This allows us to distinguish between thermal and non-thermal emission, which can be used to determine the nature of each source. We examine the spectral index of ~300 known HII regions, for which we find thermal emission with spectral indices around alpha = 0. In contrast, supernova remnants (SNR) show non-thermal emission with alpha = -0.5 and extragalactic objects generally have a steeper spectral index of alpha = -1. Using the spectral index information of the THOR survey, we investigate potential SNR candidates. We classify the radiation of four SNR candidates as non-thermal, and for the first time, we provide strong evidence for the SNR origin of these candidates.

preprint2015arXiv

Fervent: Chemistry-coupled, ionising and non-ionising radiative feedback in magnetohydrodynamical simulations

We introduce a radiative transfer code module for the magnetohydrodynamical adaptive mesh refinement code FLASH 4. It is coupled to an efficient chemical network which explicitly tracks the three hydrogen species H, H_2, H+ as well as C+ and CO. The module is geared towards modeling all relevant thermal feedback processes of massive stars, and is able to follow the non-equilibrium time-dependent thermal and chemical state of the present-day interstellar medium as well as that of dense molecular clouds. We describe in detail the implementation of all relevant thermal stellar feedback mechanisms, i.e. photoelectric, photoionization and H_2 dissociation heating as well as pumping of molecular hydrogen by UV photons. All included radiative feedback processes are extensively tested. We also compare our module to dedicated photon-dominated region (PDR) codes and find good agreement in our modeled hydrogen species once our radiative transfer solution reaches equilibrium. In addition, we show that the implemented radiative feedback physics is insensitive to the spatial resolution of the code and show under which conditions it is possible to obtain well-converged evolution in time. Finally, we briefly explore the robustness of our scheme for treating combined ionizing and non-ionizing radiation.

preprint2015arXiv

Modelling the supernova-driven ISM in different environments

We use hydrodynamical simulations in a $(256\;{\rm pc})^3$ periodic box to model the impact of supernova (SN) explosions on the multi-phase interstellar medium (ISM) for initial densities $n = 0.5-30$ cm$^{-3}$ and SN rates $1-720$ Myr$^{-1}$. We include radiative cooling, diffuse heating, and the formation of molecular gas using a chemical network. The SNe explode either at random positions, at density peaks, or both. We further present a model combining thermal energy for resolved and momentum input for unresolved SNe. Random driving at high SN rates results in hot gas ($T\gtrsim 10^6$ K) filling $> 90$% of the volume. This gas reaches high pressures ($10^4 < P/k_\mathrm{B} < 10^7$ K cm$^{-3}$) due to the combination of SN explosions in the hot, low density medium and confinement in the periodic box. These pressures move the gas from a two-phase equilibrium to the single-phase, cold branch of the cooling curve. The molecular hydrogen dominates the mass ($>50$%), residing in small, dense clumps. Such a model might resemble the dense ISM in high-redshift galaxies. Peak driving results in huge radiative losses, producing a filamentary ISM with virtually no hot gas, and a small molecular hydrogen mass fraction ($\ll 1$%). Varying the ratio of peak to random SNe yields ISM properties in between the two extremes, with a sharp transition for equal contributions. The velocity dispersion in HI remains $\lesssim 10$ km s$^{-1}$ in all cases. For peak driving the velocity dispersion in H$_α$ can be as high as $70$ km s$^{-1}$ due to the contribution from young, embedded SN remnants.

preprint2015arXiv

THOR - The HI, OH, Recombination Line Survey of the Milky Way - The pilot study: HI observations of the giant molecular cloud W43

To study the atomic, molecular and ionized emission of Giant Molecular Clouds (GMCs), we have initiated a Large Program with the VLA: 'THOR - The HI, OH, Recombination Line survey of the Milky Way'. We map the 21cm HI line, 4 OH lines, 19 H_alpha recombination lines and the continuum from 1 to 2 GHz of a significant fraction of the Milky Way (l=15-67deg, |b|<1deg) at ~20" resolution. In this paper, we focus on the HI emission from the W43 star-formation complex. Classically, the HI 21cm line is treated as optically thin with column densities calculated under this assumption. This might give reasonable results for regions of low-mass star-formation, however, it is not sufficient to describe GMCs. We analyzed strong continuum sources to measure the optical depth, and thus correct the HI 21cm emission for optical depth effects and weak diffuse continuum emission. Hence, we are able to measure the HI mass of W43 more accurately and our analysis reveals a lower limit of M=6.6x10^6 M_sun, which is a factor of 2.4 larger than the mass estimated with the assumption of optically thin emission. The HI column densities are as high as N(HI)~150 M_sun/pc^2 ~ 1.9x10^22 cm^-2, which is an order of magnitude higher than for low mass star formation regions. This result challenges theoretical models that predict a threshold for the HI column density of ~10 M_sun/pc^2, at which the formation of molecular hydrogen should set in. By assuming an elliptical layered structure for W43, we estimate the particle density profiles. While at the cloud edge atomic and molecular hydrogen are well mixed, the center of the cloud is strongly dominated by molecular hydrogen. We do not identify a sharp transition between hydrogen in atomic and molecular form. Our results are an important characterization of the atomic to molecular hydrogen transition in an extreme environment and challenges current theoretical models.

preprint2015arXiv

TOPoS: II. On the bimodality of carbon abundance in CEMP stars. Implications on the early chemical evolution of galaxies

In the course of the TOPoS (Turn Off Primordial Stars) survey, aimed at discovering the lowest metallicity stars, we have found several carbon-enhanced metal-poor (CEMP) stars. We here present our analysis of six CEMP stars. Calcium and carbon are the only elements that can be measured in all six stars. The range is -5.0<=[Ca/H]< -2.1 and 7.12<=A(C)<=8.65. For star SDSS J1742+2531 we were able to detect three FeI lines from which we deduced [Fe/H]=-4.80, from four CaII lines we derived [Ca/H]=-4.56, and from synthesis of the G-band we derived A(C)=7.26. For SDSS J1035+0641 we were not able to detect any iron lines, yet we could place a robust (3sigma) upper limit of [Fe/H]< -5.0 and measure the Ca abundance, with [Ca/H]=-5.0, and carbon, A(C)=6.90. No lithium is detected in the spectrum of SDSS J1742+2531 or SDSS J1035+0641, which implies a robust upper limit of A(Li)<1.8 for both stars. Our measured carbon abundances confirm the bimodal distribution of carbon in CEMP stars, identifying a high-carbon band and a low-carbon band. We propose an interpretation of this bimodality according to which the stars on the high-carbon band are the result of mass transfer from an AGB companion, while the stars on the low-carbon band are genuine fossil records of a gas cloud that has also been enriched by a faint supernova (SN) providing carbon and the lighter elements. (Abridged)

preprint2014arXiv

The identification of filaments on far infrared and submillimiter images. Morphology, physical conditions and relation with star formation of filamentary structure

Observations of molecular clouds reveal a complex structure, with gas and dust often arranged in filamentary rather than spherical geometries. The associations of pre- and proto- stellar cores with the filaments suggest a direct link with the process of star formation. Any study of the properties of such filaments requires a representative samples from different enviroments and so an unbiased detection method. We developed such an approach using the Hessian matrix of a surface-brightness distribution to identify filaments and determine their physical and morphological properties. After testing the method on simulated, but realistic filaments, we apply the algorithms to column-density maps computed from Herschel observations of the Galactic Plane obtained by the Hi-GAL project. We identified ~500 filaments, in the longitude range of l=216.5 to l=225.5, with lengths from ~1 pc up to ~30 pc and widths between 0.1 pc and 2.5 pc. Average column densities are between 10^20 cm^-2 and 10^22 cm^-2. Filaments include the majority of dense material with N_H_2 > 6x10^21cm^-2. We find that the pre- and proto-stellar compact sources already identified in the same region are mostly associated with filaments. However, surface densities in excess of the expected critical values for high-mass star formation are only found on the filaments, indicating that these structures are necessary to channel material into the clumps. Furthermore, we analyze the gravitational stability of filaments and discuss their relationship with star formation.

preprint2014arXiv

The SILCC (SImulating the LifeCycle of molecular Clouds) project: I. Chemical evolution of the supernova-driven ISM

The SILCC project (SImulating the Life-Cycle of molecular Clouds) aims at a more self-consistent understanding of the interstellar medium (ISM) on small scales and its link to galaxy evolution. We simulate the evolution of the multi-phase ISM in a 500 pc x 500 pc x 10 kpc region of a galactic disc, with a gas surface density of $Σ_{_{\rm GAS}} = 10 \;{\rm M}_\odot/{\rm pc}^2$. The Flash 4.1 simulations include an external potential, self-gravity, magnetic fields, heating and radiative cooling, time-dependent chemistry of H$_2$ and CO considering (self-) shielding, and supernova (SN) feedback. We explore SN explosions at different (fixed) rates in high-density regions (peak), in random locations (random), in a combination of both (mixed), or clustered in space and time (clustered). Only random or clustered models with self-gravity (which evolve similarly) are in agreement with observations. Molecular hydrogen forms in dense filaments and clumps and contributes 20% - 40% to the total mass, whereas most of the mass (55% - 75%) is in atomic hydrogen. The ionised gas contributes <10%. For high SN rates (0.5 dex above Kennicutt-Schmidt) as well as for peak and mixed driving the formation of H$_2$ is strongly suppressed. Also without self-gravity the H$_2$ fraction is significantly lower ($\sim$ 5%). Most of the volume is filled with hot gas ($\sim$90% within $\pm$2 kpc). Only for random or clustered driving, a vertically expanding warm component of atomic hydrogen indicates a fountain flow. Magnetic fields have little impact on the final disc structure. However, they affect dense gas ($n\gtrsim 10\;{\rm cm}^{-3}$) and delay H$_2$ formation. We highlight that individual chemical species, in particular atomic hydrogen, populate different ISM phases and cannot be accurately accounted for by simple temperature-/density-based phase cut-offs.

preprint2013arXiv

Molecular cooling in the diffuse interstellar medium

We use a simple one-zone model of the thermal and chemical evolution of interstellar gas to study whether molecular hydrogen (H2) is ever an important coolant of the warm, diffuse interstellar medium (ISM). We demonstrate that at solar metallicity, H2 cooling is unimportant and the thermal evolution of the ISM is dominated by metal line cooling. At metallicities below 0.1 Z_solar, however, metal line cooling of low density gas quickly becomes unimportant and H2 can become the dominant coolant, even though its abundance in the gas remains small. We investigate the conditions required in order for H2 to dominate, and show that it provides significant cooling only when the ratio of the interstellar radiation field strength to the gas density is small. Finally, we demonstrate that our results are insensitive to changes in the initial fractional ionization of the gas or to uncertainties in the nature of the dust present in the low-metallicity ISM.

preprint2012arXiv

The Density Variance--Mach Number Relation in Supersonic Turbulence: I. Isothermal, magnetised gas

It is widely accepted that supersonic, magnetised turbulence plays a fundamental role for star formation in molecular clouds. It produces the initial dense gas seeds out of which new stars can form. However, the exact relation between gas compression, turbulent Mach number, and magnetic field strength is still poorly understood. Here, we introduce and test an analytical prediction for the relation between the density variance and the root-mean-square Mach number in supersonic, isothermal, magnetised turbulent flows. We approximate the density and velocity structure of the interstellar medium as a superposition of shock waves. We obtain the density contrast considering the momentum equation for a single magnetised shock and extrapolate this result to the entire cloud. Depending on the field geometry, we then make three different assumptions based on observational and theoretical constraints: B independent of density, B proportional to the root square of the density and B proportional to the density. We test the analytically derived density variance--Mach number relation with numerical simulations, and find that for B proportional to the root square of the density, the variance in the logarithmic density contrast, $σ_{\ln ρ/ρ_0}^2=\ln[1+b^2\mathscr{M}^2β_0/(β_0+1)]$, fits very well to simulated data with turbulent forcing parameter b=0.4, when the gas is super-Alfvénic. However, this result breaks down when the turbulence becomes trans-Alfvénic or sub-Alfvénic, because in this regime the turbulence becomes highly anisotropic. Our density variance--Mach number relations simplify to the purely hydrodynamic relation as the ratio of thermal to magnetic pressure $β_0$ approaches infinite.

preprint2011arXiv

Approximations for modelling CO chemistry in GMCs: a comparison of approaches

We examine several different simplified approaches for modelling the chemistry of CO in three-dimensional numerical simulations of turbulent molecular clouds. We compare the different models both by looking at the behaviour of integrated quantities such as the mean CO fraction or the cloud-averaged CO-to-H2 conversion factor, and also by studying the detailed distribution of CO as a function of gas density and visual extinction. In addition, we examine the extent to which the density and temperature distributions depend on our choice of chemical model. We find that all of the models predict the same density PDF and also agree very well on the form of the temperature PDF for temperatures T > 30 K, although at lower temperatures, some differences become apparent. All of the models also predict the same CO-to-H2 conversion factor, to within a factor of a few. However, when we look more closely at the details of the CO distribution, we find larger differences. The more complex models tend to produce less CO and more atomic carbon than the simpler models, suggesting that the C/CO ratio may be a useful observational tool for determining which model best fits the observational data. Nevertheless, the fact that these chemical differences do not appear to have a strong effect on the density or temperature distributions of the gas suggests that the dynamical behaviour of the molecular clouds on large scales is not particularly sensitive to how accurately the small-scale chemistry is modelled.

preprint2011arXiv

Is molecular gas necessary for star formation?

On galactic scales, the surface density of star formation appears to be well correlated with the surface density of molecular gas. This has lead many authors to suggest that there exists a causal relationship between the chemical state of the gas and its ability to form stars -- in other words, the assumption that the gas must be molecular before star formation can occur. We test this hypothesis by modelling star formation within a dense cloud of gas with properties similar to a small molecular cloud using a series of different models of the chemistry, ranging from one in which the formation of molecules is not followed and the gas is assumed to remain atomic throughout, to one that tracks the formation of both H2 and CO. We find that presence of molecules in the gas has little effect on the ability of the gas to form stars: star formation can occur just as easily in atomic gas as in molecular gas. At low densities (< 10^4 cm^-3), the gas is able to cool via C+ fine-structure emission almost as efficiently as via CO rotational line emission, while at higher densities, the main cooling process involves the transfer of energy from gas to dust, meaning that the presence of molecules is again unimportant. Cooling by H2 is particularly inefficient, accounting for as little as 1 percent of the overall cooling in the cloud. Rather than the chemical makeup, we find that the most important factor controlling the rate of star formation is the ability of the gas to shield itself from the interstellar radiation field. As this is also a prerequisite for the survival of molecules within the gas, our results support a picture in which molecule formation and the formation of cold gas are both correlated with the column density of the cloud -- and thus its ability to shield itself -- rather than being directly correlated with each other.

preprint2011arXiv

The Chemistry of the Early Universe

The chemistry of the early Universe is a fascinating field of study. Even in the absence of any elements heavier than lithium, a surprising degree of chemical complexity proves to be possible, giving the topic considerable interest in its own right. In addition, the fact that molecular hydrogen plays a key role in the formation of the first stars and galaxies means that if we want to understand the formation of these objects, we must first develop a good understanding of the chemical evolution of the gas. In this review, I first give a brief introduction to the chemistry occurring in the gas prior to the formation of the first stars and galaxies, and then go on to discuss in more detail the main chemical processes occurring during the gravitational collapse of gas from intergalactic to protostellar densities, and how these processes influence the final outcome of the collapse.

preprint2011arXiv

Understanding the physics of the X-factor

We study the relationship between the H2 and CO abundances in simulated molecular clouds using a fully dynamical model of magnetized turbulence coupled to a detailed chemical network. We find that the CO-to-H2 conversion factor for a given molecular cloud, the so-called X-factor, is determined primarily by the mean extinction of the cloud, rather than by its metallicity. Our results explain the discrepancy observed in low metallicity systems between cloud masses derived from CO observations and other techniques such as infrared emission, and predict that CO-bright clouds in low metallicity systems should be systematically larger and/or denser than Milky Way clouds.

preprint2010arXiv

Effects of Varying the Three-Body Molecular Hydrogen Formation Rate in Primordial Star Formation

The transformation of atomic hydrogen to molecular hydrogen through three-body reactions is a crucial stage in the collapse of primordial, metal-free halos, where the first generation of stars (Population III stars) in the Universe are formed. However, in the published literature, the rate coefficient for this reaction is uncertain by nearly an order of magnitude. We report on the results of both adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) simulations of the collapse of metal-free halos as a function of the value of this rate coefficient. For each simulation method, we have simulated a single halo three times, using three different values of the rate coefficient. We find that while variation between halo realizations may be greater than that caused by the three-body rate coefficient being used, both the accretion physics onto Population III protostars as well as the long-term stability of the disk and any potential fragmentation may depend strongly on this rate coefficient.

preprint2010arXiv

Fragmentation in turbulent primordial gas

We report results from numerical simulations of star formation in the early universe that focus on the role of subsonic turbulence, and investigate whether it can induce fragmentation of the gas. We find that dense primordial gas is highly susceptible to fragmentation, even for rms turbulent velocity dispersions as low as 20% of the initial sound speed. The resulting fragments cover over two orders of magnitude in mass, ranging from 0.1 to 40 solar masses. However, our results suggest that the details of the fragmentation depend on the local properties of the turbulent velocity field and hence we expect considerable variations in the resulting stellar mass spectrum in different halos.

preprint2010arXiv

On the Relationship Between Molecular Hydrogen and Carbon Monoxide Abundances in Molecular Clouds

The most usual tracer of molecular gas is line emission from CO. However, the reliability of that tracer has long been questioned in environments different from the Milky Way. We study the relationship between H2 and CO abundances using a fully dynamical model of magnetized turbulence coupled to a chemical network simplified to follow only the dominant pathways for H2 and CO formation and destruction, and including photodissociation using a six-ray approximation. We find that the abundance of H2 is primarily determined by the amount of time available for its formation, which is proportional to the product of the density and the metallicity, but insensitive to photodissociation. Photodissociation only becomes important at extinctions under a few tenths of a visual magnitude, in agreement with both observational and prior theoretical work. On the other hand, CO forms quickly, within a dynamical time, but its abundance depends primarily on photodissociation, with only a weak secondary dependence on H2 abundance. As a result, there is a sharp cutoff in CO abundance at mean visual extinctions A_V < 3. At lower values of A_V we find that the ratio of H2 column density to CO emissivity X_CO is proportional to A_V^(-3.5). This explains the discrepancy observed in low metallicity systems between cloud masses derived from CO observations and other techniques such as infrared emission. Our work predicts that CO-bright clouds in low metallicity systems should be systematically larger or denser than Milky Way clouds, or both. Our results further explain the narrow range of observed molecular cloud column densities as a threshold effect, without requiring the assumption of virial equilibrium.

preprint2009arXiv

Modelling CO formation in the turbulent interstellar medium

We present results from high-resolution three-dimensional simulations of turbulent interstellar gas that self-consistently follow its coupled thermal, chemical and dynamical evolution, with a particular focus on the formation and destruction of H2 and CO. We quantify the formation timescales for H2 and CO in physical conditions corresponding to those found in nearby giant molecular clouds, and show that both species form rapidly, with chemical timescales that are comparable to the dynamical timescale of the gas. We also investigate the spatial distributions of H2 and CO, and how they relate to the underlying gas distribution. We show that H2 is a good tracer of the gas distribution, but that the relationship between CO abundance and gas density is more complex. The CO abundance is not well-correlated with either the gas number density n or the visual extinction A_V: both have a large influence on the CO abundance, but the inhomogeneous nature of the density field produced by the turbulence means that n and A_V are only poorly correlated. There is a large scatter in A_V, and hence CO abundance, for gas with any particular density, and similarly a large scatter in density and CO abundance for gas with any particular visual extinction. This will have important consequences for the interpretation of the CO emission observed from real molecular clouds. Finally, we also examine the temperature structure of the simulated gas. We show that the molecular gas is not isothermal. Most of it has a temperature in the range of 10--20 K, but there is also a significant fraction of warmer gas, located in low-extinction regions where photoelectric heating remains effective.

preprint2009arXiv

Star Formation at Very Low Metallicity. V. The Greater Importance of Initial Conditions Compared to Metallicity Thresholds

The formation of the first stars out of metal-free gas appears to result in stars at least an order of magnitude more massive than in the present-day case. We here consider what controls the transition from a primordial to a modern initial mass function. It has been proposed that this occurs when effective metal line cooling occurs at a metallicity threshold of Z/Z_sun > 10^{-3.5}. We study the influence of low levels of metal enrichment on the cooling and collapse of initially ionized gas in small protogalactic halos using three-dimensional, smoothed particle hydrodynamics simulations with particle splitting. Our initial conditions represent protogalaxies forming within a previously ionized H ii region that has not yet had time to cool and recombine. These differ considerably from those used in simulations predicting a metallicity threshold, where the gas was initially cold and only partially ionized. In the centrally condensed potential that we study here, a wide variety of initial conditions for the gas yield a monolithic central collapse. Our models show no fragmentation during collapse to number densities as high as 10^5 cm^{-3}, for metallicities reaching as high as 10^{-1} Z_sun in one rotating case, far above the threshold suggested by previous work. Rotation allows for the formation of gravitationally stable gas disks over large fractions of the local Hubble time. Turbulence slows the growth of the central density slightly, but both spherically symmetric and turbulent initial conditions collapse and form a single sink particle. We therefore argue that fragmentation at moderate density depends on the initial conditions for star formation more than on the metal abundances present.