Researcher profile

Simon C. O. Glover

Simon C. O. Glover contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
42works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

42 published item(s)

preprint2025arXiv

Evolved Supergiants in PHANGS I: Red Supergiants in 19 Galaxies between 5-20 Mpc with HST and JWST

Red supergiants (RSGs) are important for our understanding of supernova progenitors, stellar populations, stellar evolution, mass loss and dust production. Extragalactic surveys of RSGs have a long history in the Local Group, but few studies exist beyond that due to the limited resolution and sensitivity of ground-based and previous space-based infrared observatories. Here we demonstrate the combined power of HST and JWST to push systematic searches of RSGs out to $\sim$20 Mpc. We introduce a catalog of 97057 RSGs -- the largest single-survey release of RSGs -- with masses $\gtrsim$10 M$_{\odot}$ in 19 galaxies from the PHANGS HST+JWST Treasury program. We use HST F814W and JWST F200W photometry to select stars as RSGs based on predicted colors and magnitudes from PARSEC isochrones. The spatial distribution of our recovered RSGs follow the familiar pattern of mostly being concentrated in active star-forming regions such as spiral arms and central starburst rings. The RSG number density on kpc-scales is strongly correlated ($r_s$$\sim$0.82) with local star-formation rate density ($Σ_{SFR}$) traced by extinction-corrected far-ultraviolet (FUV) from GALEX+WISE, and weakly correlated ($r_s$$\sim$0.57) with the total stellar mass density ($Σ_*$), traced by near-infrared emission from WISE+Spitzer. The number of RSGs per mass of stellar populations with ages 6-30 Myr (the likely age range of RSGs $>$10 M$_{\odot}$) is $\sim$1 per 10$^{3.77\pm0.27}$ M$_{\odot}$, assuming constant star-formation rates from FUV+W4. Our sample will be a useful resource for tracking progenitors and feedback sites of future supernovae in PHANGS, age-dating stellar populations, and more.

preprint2025arXiv

The Noctua Suite of Simulations -- The Difficulty of Growing Massive Black Holes in Low-Mass Dwarf Galaxies

Aims. We study the individual and cumulative impact of stellar feedback processes on massive black hole (MBH) growth in a simulated low-mass dwarf galaxy. Methods. A suite of high-resolution radiation-hydrodynamic simulations called Noctua is performed, using the ArepoNoctua numerical framework for BHs in galaxy simulations. The chemical evolution of the gas is explicitly modelled in a time-dependent non-equilibrium way. Two types of stellar feedback are considered: individually-traced type II supernova (SNII) explosions, and radiatively transferred (on-the-fly) ionising stellar radiation (ISR) from OB stars. As part of the numerical framework, we develop and apply a novel physically-motivated model for MBH gas accretion, taking into account the angular momentum of the gas in the radiatively efficient regime, to estimate the gas accretion rate from the sub-grid accretion disc. Results. Without any stellar feedback, an initial $10^4~\mathrm{M}_\odot$ MBH is able to steadily grow over time, roughly doubling its mass after 800 Myr. Surprisingly, the growth of the MBH is more than doubled when only ISR feedback is considered, compared to the no stellar feedback run. This is due to the star formation rate (SFR) being highly suppressed (to a similar level or slightly above that when SNII feedback is considered), enabling a higher cumulative net gas inflow onto the MBH from not only the cold neutral- and molecular medium phases, but also the unstable- and warm neutral medium phases. With SNII feedback included, the gas accretion onto the MBH is episodic over time, and is suppressed by more than an order of magnitude already during the first 150 Myr. When combining SNII with ISR feedback, the growth of the MBH remains suppressed due to SNII feedback, but to a lesser extent compared to the SNII-only feedback run, due to a slightly lower SFR, and hence a reduced number of SNII events.

preprint2023arXiv

Magnetic fields do not suppress global star formation in low metallicity dwarf galaxies

Many studies concluded that magnetic fields suppress star formation in molecular clouds and Milky Way like galaxies. However, most of these studies are based on fully developed fields that have reached the saturation level, with little work on investigating how an initial weak primordial field affects star formation in low metallicity environments. In this paper, we investigate the impact of a weak initial field on low metallicity dwarf galaxies. We perform high-resolution AREPO simulations of five isolated dwarf galaxies. Two models are hydrodynamical, two start with a primordial magnetic field of 10$^{-6} μ$G and different sub-solar metallicities, and one starts with a saturated field of 10$^{-2} μ$G. All models include a non-equilibrium, time-dependent chemical network that includes the effects of gas shielding from the ambient ultraviolet field. Sink particles form directly from the gravitational collapse of gas and are treated as star-forming clumps that can accrete gas. We vary the ambient uniform far ultraviolet field, and cosmic ray ionization rate between 1\% and 10\% of solar values. We find that the magnetic field has little impact on the global star formation rate, which is in tension with some previously published results. We further find that the initial field strength has little impact on the global star formation rate. We show that an increase in the mass fractions of both molecular hydrogen and cold gas, along with changes in the perpendicular gas velocity dispersion and the magnetic field acting in the weak-field model, overcome the expected suppression in star formation.

preprint2023arXiv

Towards the impact of GMC collisions on the star formation rate

Collisions between giant molecular clouds (GMCs) are one of the pathways for massive star formation, due to the high densities created. However the enhancement of the star formation rate (SFR) is not well constrained. In this study we perform a parameter study of cloud-cloud collisions, and investigate how the resulting SFR depends on the details of set-up. Our parameter study explores variations in: collision speed; magnetic field inclination (with respect to the collisional axis); and resolution, as defined by the number of cells per Jeans length. In all our collision simulations we find a factor of 2-3 increase in the SFR compared to our no collision simulation, with star formation beginning sooner with a) high collisional velocities, b) parallel orientation between the magnetic field and collision axis, c) and lower resolution. The mean virial parameter of high density (and thus possible star-forming) gas increases with collisional velocity, but has little variation with magnetic field inclination. The alignment of the velocity and magnetic field remains uniform in low density environments but becomes more perpendicular with increasing density, indicating the compression of the magnetic field by collapsing gas. Comparing the trends in the SFR with other GMC collision studies, we find good agreement with studies that account for the gravitational boundedness of the gas in their star formation algorithm, but not with those that simply form stars above a prescribed density threshold. This suggests that the latter approach should be used with caution when modelling star formation on resolved cloud scales.

preprint2022arXiv

A 2-3 mm high-resolution molecular line survey towards the centre of the nearby spiral galaxy NGC 6946

The complex physical, kinematic, and chemical properties of galaxy centres make them interesting environments to examine with molecular line emission. We present new $2-4$" (${\sim}75{-}150$ pc at $7.7$ Mpc) observations at 2 and 3 mm covering the central $50$" (${\sim}1.9$ kpc) of the nearby double-barred spiral galaxy NGC 6946 obtained with the IRAM Plateau de Bure Interferometer. We detect spectral lines from ten molecules: CO, HCN, HCO$^+$, HNC, CS, HC$_3$N, N$_2$H$^+$, C$_2$H, CH$_3$OH, and H$_2$CO. We complemented these with published 1mm CO observations and 33 GHz continuum observations to explore the star formation rate surface density ${Σ_{\mathrm{SFR}}}$ on 150 pc scales. In this paper, we analyse regions associated with the inner bar of NGC 6946 $-$ the nuclear region (NUC), the northern (NBE), and southern inner bar end (SBE) and we focus on short-spacing corrected bulk (CO) and dense gas tracers (HCN, HCO$^+$, and HNC). We find that HCO$^+$ correlates best with ${Σ_{\mathrm{SFR}}}$, but the dense gas fraction ($f_{\mathrm{dense}}$) and star formation efficiency of the dense gas (${\mathrm{SFE_{dense}}}$) fits show different behaviours than expected from large-scale disc observations.The SBE has a higher ${Σ_{\mathrm{SFR}}}$, $f_{\mathrm{dense}}$, and shocked gas fraction than the NBE. We examine line ratio diagnostics and find a higher CO(2-1)/CO(1-0) ratio towards NBE than for the NUC. Moreover, comparison with existing extragalactic datasets suggests that using the HCN/HNC ratio to probe kinetic temperatures is not suitable on kiloparsec and sub-kiloparsec scales in extragalactic regions. Lastly, our study shows that the HCO$^+$/HCN ratio might not be a unique indicator to diagnose AGN activity in galaxies.

preprint2022arXiv

A CO isotopologue Line Atlas within the Whirlpool galaxy Survey (CLAWS)

We present the CO isotopologue Line Atlas within the Whirpool galaxy Survey (CLAWS) based on an IRAM 30-m large programme which provides a benchmark study of numerous, faint CO isotopologues in the mm-wavelength regime across the full disc of M51 (NGC 5194). The survey's core goal is to use the low-J CO isotopologue lines to constrain CO excitation and chemistry, and therefrom the local physical conditions of the gas. In this survey paper, we describe the CLAWS observing and data reduction strategies. We map the J=1-0 and 2-1 transitions of the CO isotopologues $^{12}$CO,$^{13}$CO, C$^{18}$O and C$^{17}$O, as well as several supplementary lines within the 1 mm and 3 mm window (CN(1-0), CS(2-1), CH$_3$OH(2-1), N$_2$H$^+$(1-0), HC$_3$N(10-9)) at ~1 kpc resolution. A total observation time of 149 h offers unprecedented sensitivity. We use these data to explore several CO isotopologue line ratios in detail, study their radial (and azimuthal) trends and investigate whether changes in line ratios stem from changes in ISM properties such as gas temperatures, densities or chemical abundances. For example, we find negative radial trends for the $^{13}$CO}/$^{12}$CO, C$^{18}$O/$^{12}$CO and C$^{18}$O/$^{13}$CO line ratios in their J=1-0 transitions. We also find variations with local environment, such as higher $^{12}$CO(2-1)/(1-0) or $^{13}$CO/$^{12}$CO(1-0) line ratios in interarm regions compared to spiral arm regions. We propose that these aforementioned variations of CO line ratios are most likely due to a variation of the optical depth, while abundance variations due to selective nucleosynthesis on a galaxy-wide scale could also play a role. We also study the CO spectral line energy distribution (SLED) using archival JCMT $^{12}$CO(3-2) data and find a variation of the SLED shape with local environmental parameters further underlying changes in optical depth, gas temperatures or densities.

preprint2022arXiv

A tale of two DIGs: The relative role of HII regions and low-mass hot evolved stars in powering the diffuse ionised gas (DIG) in PHANGS-MUSE galaxies

We use integral field spectroscopy from the PHANGS-MUSE survey, which resolves the ionised interstellar medium at ${\sim}50$ pc resolution in 19 nearby spiral galaxies, to study the origin of the diffuse ionised gas (DIG). We examine the physical conditions of the diffuse gas by first removing morphologically-defined HII regions and then binning the low-surface-brightness areas to achieve significant detections of the key nebular lines. A simple model for the leakage and propagation of ionising radiation from HII regions is able to reproduce the observed distribution of H$α$ in the DIG. Leaking radiation from HII regions also explains the observed decrease in line ratios of low-ionisation species ([SII]/H$α$, [NII]/H$α$ and [OI]/H$α$) with increasing H$α$ surface brightness ($Σ_{Hα}$). Emission from hot low-mass evolved stars, however, is required to explain: (1) the enhanced low-ionisation line ratios observed in the central regions of some galaxies; (2) the observed trends of a flat or decreasing [OIII]/H$β$ with $Σ_{Hα}$; and (3) the offset of some DIG regions from the locus of HII regions in the Baldwin-Phillips-Terlevich (BPT) diagram, extending into the area of low-ionisation (nuclear) emission-line regions (LI[N]ERs). Hot low-mass evolved stars make a small contribution to the energy budget of the DIG (2% of the galaxy-integrated H$α$ emission), but their harder spectra make them fundamental contributors to [OIII] emission. The DIG might result from a superposition of two components, an energetically dominant contribution from young stars and a more diffuse background of harder ionising photons from old stars. This unified framework bridges observations of the Milky Way DIG with LI(N)ER-like emission observed in nearby galaxy bulges.

preprint2022arXiv

Emission-line diagnostics of HII regions using conditional Invertible Neural Networks

Young massive stars play an important role in the evolution of the interstellar medium (ISM) and the self-regulation of star formation in giant molecular clouds (GMCs) by injecting energy, momentum, and radiation (stellar feedback) into surrounding environments, disrupting the parental clouds, and regulating further star formation. Information of the stellar feedback inheres in the emission we observe, however inferring the physical properties from photometric and spectroscopic measurements is difficult, because stellar feedback is a highly complex and non-linear process, so that the observational data are highly degenerate. On this account, we introduce a novel method that couples a conditional invertible neural network (cINN) with the WARPFIELD-emission predictor (WARPFIELD-EMP) to estimate the physical properties of star-forming regions from spectral observations. We present a cINN that predicts the posterior distribution of seven physical parameters (cloud mass, star formation efficiency, cloud density, cloud age which means age of the first generation stars, age of the youngest cluster, the number of clusters, and the evolutionary phase of the cloud) from the luminosity of 12 optical emission lines, and test our network with synthetic models that are not used during training. Our network is a powerful and time-efficient tool that can accurately predict each parameter, although degeneracy sometimes remains in the posterior estimates of the number of clusters. We validate the posteriors estimated by the network and confirm that they are consistent with the input observations. We also evaluate the influence of observational uncertainties on the network performance.

preprint2022arXiv

Environmental dependence of the molecular cloud lifecycle in 54 main sequence galaxies

The processes of star formation and feedback, regulating the cycle of matter between gas and stars on the scales of giant molecular clouds (GMCs; $\sim$100pc), play a major role in governing galaxy evolution. Measuring the time-scales of GMC evolution is important to identify and characterise the specific physical mechanisms that drive this transition. By applying a robust statistical method to high-resolution CO and narrow-band H$α$ imaging from the PHANGS survey, we systematically measure the evolutionary timeline from molecular clouds to exposed young stellar regions on GMC scales, across the discs of an unprecedented sample of 54 star-forming main-sequence galaxies (excluding their unresolved centres). We find that clouds live for about $1{-}3$ GMC turbulence crossing times ($5{-}30$Myr) and are efficiently dispersed by stellar feedback within $1{-}5$Myr once the star-forming region becomes partially exposed, resulting in integrated star formation efficiencies of $1{-}8$%. These ranges reflect physical galaxy-to-galaxy variation. In order to evaluate whether galactic environment influences GMC evolution, we correlate our measurements with average properties of the GMCs and their local galactic environment. We find several strong correlations that can be physically understood, revealing a quantitative link between galactic-scale environmental properties and the small-scale GMC evolution. Notably, the measured CO-visible cloud lifetimes become shorter with decreasing galaxy mass, mostly due to the increasing presence of CO-dark molecular gas in such environment. Our results represent a first step towards a comprehensive picture of cloud assembly and dispersal, which requires further extension and refinement with tracers of the atomic gas, dust, and deeply-embedded stars.

preprint2022arXiv

Metal Mixing in Minihalos: The Descendants of Pair-Instability Supernovae

The lack of observations of abundance patterns originating in pair-instability supernovae has been a long-standing problem in relation to the first stars. This class of supernovae is expected to have an abundance pattern with a strong odd-even effect, making it substantially different from present-day supernovae. In this study, we use a cosmological radiation hydrodynamics simulation to model such supernovae and the subsequent formation of the second generation of stars. We incorporate streaming velocities for the first time. There are 14 star-forming minihalos in our $1\,\mathrm{Mpc}\,h^{-1}$ box, leading to 14 supernovae occurring before redshift $z=19.5$, where we start reducing the complexity of the simulation. Following the explosions, extremely metal-poor stars form in 10 halos via internal and external enrichment, which makes it the most common outcome. Only one halo does not recollapse during the simulations. This result is at tension with the current (lack of) observations of metal-poor stars with pair instability supernova abundance patterns, suggesting that these very massive stars might be rare even in the early Universe. The results from this simulation also give us insights into what drives different modes of recollapse and what determines the mixing behavior of metals after very energetic supernovae.

preprint2022arXiv

Planetary Nebula Luminosity Function distances for 19 galaxies observed by PHANGS-MUSE

We provide new planetary nebula luminosity function (PNLF) distances to 19 nearby spiral galaxies that were observed with VLT/MUSE by the PHANGS collaboration. Emission line ratios are used to separate planetary nebulae (PNe) from other bright [OIII] emitting sources like compact supernovae remnants (SNRs) or HII regions. While many studies have used narrowband imaging for this purpose, the detailed spectral line information provided by integral field unit (IFU) spectroscopy grants a more robust way of categorising different [OIII] emitters. We investigate the effects of SNR contamination on the PNLF and find that we would fail to classify all objects correctly, when limited to the same data narrowband imaging provides. However, the few misclassified objects usually do not fall on the bright end of the luminosity function, and only in three cases does the distance change by more than $1σ$. We find generally good agreement with literature values from other methods. Using metallicity constraints that have also been derived from the same IFU data, we revisit the PNLF zero point calibration. Over a range of $8.34 < 12 + \log(\mathrm{O}/\mathrm{H}) < 8.59$, our sample is consistent with a constant zero point and yields $M^*=-4.542^{+0.103}_{-0.059}\, \mathrm{mag}$, within $1σ$ of other literature values. MUSE pushes the limits of PNLF studies and makes galaxies beyond $20\, \mathrm{Mpc}$ accessible for this kind of analysis. This approach to the PNLF shows great promise for leveraging existing archival IFU data on nearby galaxies.

preprint2022arXiv

Public Release of A-SLOTH: Ancient Stars and Local Observables by Tracing Halos

The semi-analytical model A-SLOTH (Ancient Stars and Local Observables by Tracing Halos) is the first public code that connects the formation of the first stars and galaxies to observables. After several successful projects with this model, we publish the source code and describe the public version in this paper. The model is based on dark matter merger trees that can either be generated based on Extended Press-Schechter theory or that can be imported from dark matter simulations. On top of these merger trees, A-SLOTH applies analytical recipes for baryonic physics to model the formation of both metal-free and metal-poor stars and the transition between them with unprecedented precision and fidelity. A-SLOTH samples individual stars and includes radiative, chemical, and mechanical feedback. It is calibrated based on six observables, such as the optical depth to Thomson scattering, the stellar mass of the Milky Way and its satellite galaxies, the number of extremely-metal poor stars, and the cosmic star formation rate density at high redshift. A-SLOTH has versatile applications with moderate computational requirements. It can be used to constrain the properties of the first stars and high-z galaxies based on local observables, predicts properties of the oldest and most metal-poor stars in the Milky Way, can serve as a subgrid model for larger cosmological simulations, and predicts next-generation observables of the early Universe, such as supernova rates or gravitational wave events.

preprint2022arXiv

The Galactic dynamics revealed by the filamentary structure in atomic hydrogen emission

We present a study of the filamentary structure in the atomic hydrogen (HI) emission at the 21 cm wavelength toward the Galactic plane using the observations in the HI4PI survey. Using the Hessian matrix method across radial velocity channels, we identified the filamentary structures and quantified their orientations using circular statistics. We found that the regions of the Milky Way&#39;s disk beyond 10 kpc and up to roughly 18 kpc from the Galactic center display HI filamentary structures predominantly parallel to the Galactic plane. For regions at lower Galactocentric radii, we found that the HI filaments are mostly perpendicular or do not have a preferred orientation with respect to the Galactic plane. We interpret these results as the imprint of supernova feedback in the inner Galaxy and Galactic rotation in the outer Milky Way. We found that the HI filamentary structures follow the Galactic warp and that they highlight some of the variations interpreted as the effect of the gravitational interaction with satellite galaxies. In addition, the mean scale height of the filamentary structures is lower than that sampled by the bulk of the HI emission, thus indicating that the cold and warm atomic hydrogen phases have different scale heights in the outer galaxy. Finally, we found that the fraction of the column density in HI filaments is almost constant up to approximately 18 kpc from the Galactic center. This is possibly a result of the roughly constant ratio between the cold and warm atomic hydrogen phases inferred from the HI absorption studies. Our results indicate that the HI filamentary structures provide insight into the dynamical processes shaping the Galactic disk. Their orientations record how and where the stellar energy input, the Galactic fountain process, the cosmic ray diffusion, and the gas accretion have molded the diffuse interstellar medium in the Galactic plane.

preprint2022arXiv

The Gas-Star Formation Cycle in Nearby Star-forming Galaxies II. Resolved Distributions of CO and H$α$ Emission for 49 PHANGS Galaxies

The relative distribution of molecular gas and star formation in galaxies gives insight into the physical processes and timescales of the cycle between gas and stars. In this work, we track the relative spatial configuration of CO and H$α$ emission at high resolution in each of our galaxy targets, and use these measurements to quantify the distributions of regions in different evolutionary stages of star formation: from molecular gas without star formation traced by H$α$ to star-forming gas, and to HII regions. The large sample, drawn from the Physics at High Angular resolution in Nearby GalaxieS ALMA and narrowband H$α$ (PHANGS-ALMA and PHANGS-H$α$) surveys, spans a wide range of stellar mass and morphological types, allowing us to investigate the dependencies of the gas-star formation cycle on global galaxy properties. At a resolution of 150 pc, the incidence of regions in different stages shows a dependence on stellar mass and Hubble type of galaxies over the radial range probed. Massive and/or earlier-type galaxies exhibit a significant reservoir of molecular gas without star formation traced by H$α$, while lower-mass galaxies harbor substantial HII regions that may have dispersed their birth clouds or formed from low-mass, more isolated clouds. Galactic structures add a further layer of complexity to relative distribution of CO and H$α$ emission. Trends between galaxy properties and distributions of gas traced by CO and H$α$ are visible only when the observed spatial scale is $\ll$ 500 pc, reflecting the critical resolution requirement to distinguish stages of star formation process.

preprint2022arXiv

The PHANGS-MUSE survey -- Probing the chemo-dynamical evolution of disc galaxies

We present the PHANGS-MUSE survey, a programme using the MUSE IFS at the ESO VLT to map 19 massive $(9.4 < \log(M_{*}/M_\odot) < 11.0)$ nearby (D < 20 Mpc) star-forming disc galaxies. The survey consists of 168 MUSE pointings (1&#39;x1&#39; each), a total of nearly 15 Million spectra, covering ~1.5 Million independent spectra. PHANGS-MUSE provides the first IFS view of star formation across different local environments (including galaxy centres, bars, spiral arms) in external galaxies at a median resolution of 50~pc, better than the mean inter-cloud distance in the ionised interstellar medium. This `cloud-scale&#39; resolution allows detailed demographics and characterisations of HII regions and other ionised nebulae. PHANGS-MUSE further delivers a unique view on the associated gas and stellar kinematics, and provides constraints on the star formation history. The PHANGS-MUSE survey is complemented by dedicated ALMA CO(2-1) and multi-band HST observations, therefore allowing us to probe the key stages of the star formation process from molecular clouds to HII regions and star clusters. This paper describes the scientific motivation, sample selection, observational strategy, data reduction and analysis process of the PHANGS-MUSE survey. We present our bespoke automated data-reduction framework, which is built on the reduction recipes provided by ESO, but additionally allows for mosaicking and homogenisation of the point spread function. We further present a detailed quality assessment and a brief illustration of the potential scientific applications of the large set of PHANGS-MUSE data products generated by our data analysis framework. The data cubes and analysis data products described in this paper represent the basis for the first PHANGS-MUSE public data release and are available in the ESO archive and via the Canadian Astronomy Data Centre.

preprint2022arXiv

Tracing stars in Milky Way satellites with A-SLOTH

We study the stellar mass-to-halo mass relation at $z=0$ in 30 Milky Way-like systems down to the ultra-faint ($M_* < 10^5 M_\odot$) regime using the semi-analytic model A-SLOTH. A new model allows us to follow star formation and the stochastic stellar feedback from individually sampled Pop II stars. Our fiducial model produces consistent results with the stellar mass-to-halo mass relation derived from abundance matching and the observed cumulative stellar mass function above the observational completeness. We find a plateau in the stellar mass-to-halo mass relation in the ultra-faint regime. The stellar mass of this plateau tells us how many stars formed before supernovae occur and regulate further star formation, which is determined by the Pop~II star formation efficiency. We also find that the number of luminous satellites increases rapidly as $M_*$ decreases until $M_* \approx 10^4 M_\odot$. Finally, we find that the relative streaming velocity between baryons and dark matter at high redshift is important in determining the number of ultra-faint dwarf galaxies at $z=0$. The new model in A-SLOTH provides a framework to study the stellar properties and the formation history of metal-poor stars in Milky Way and its satellites.

preprint2022arXiv

Trapping of HII regions in Population III star formation

Radiative feedback from massive Population III (Pop III) stars in the form of ionising and photodissociating photons is widely believed to play a central role in shutting off accretion onto these stars. Understanding whether and how this occurs is vital for predicting the final masses reached by these stars and the form of the Pop III stellar initial mass function. To help us better understand the impact of UV radiation from massive Pop III stars on the gas surrounding them, we carry out high resolution simulations of the formation and early evolution of these stars, using the AREPO moving-mesh code coupled with the innovative radiative transfer module SPRAI. Contrary to most previous results, we find that the ionising radiation from these stars is trapped in the dense accretion disk surrounding them. Consequently, the inclusion of radiative feedback has no significant impact on either the number or the total mass of protostars formed during the 20 kyr period that we simulate. We show that the reason that we obtain qualitatively different results from previous studies of Pop III stellar feedback lies in how the radiation is injected into the simulation. HII region trapping only occurs if the photons are injected on scales smaller than the local scale height of the accretion disk, a criterion not fulfilled in previous 3D simulations of this process. Finally, we speculate as to whether outflows driven by the magnetic field or by Lyman-alpha radiation pressure may be able to clear enough gas away from the star to allow the HII region to escape from the disk.

preprint2021arXiv

Effect of the cosmological transition to metal-enriched star-formation on the hydrogen 21-cm signal

Mapping Cosmic Dawn with 21-cm tomography offers an exciting new window into the era of primordial star formation. However, self-consistent implementation of both the process of star formation and the related 21-cm signal is challenging, due to the multi-scale nature of the problem. In this study, we develop a flexible semi-analytical model to follow the formation of the first stars and the process of gradual transition from primordial to metal-enriched star formation. For this transition we use different in scenarios with varying time-delays (or recovery times) between the first supernovae and the formation of the second generation of stars. We use recovery times between 10 and 100\,Myr and find that these delays have a strong impact on the redshift at which the transition to metal-enriched star formation occurs. We then explore the effect of this transition on the 21-cm signal and find that the recovery time has a distinctive imprint in the signal. Together with an improved understanding of how this time-delay relates to the properties of Population~III stars, future 21-cm observations can give independent constraints on the earliest epoch of star formation.

preprint2021arXiv

Fragmentation induced starvation in Population III star formation: a resolution study

The Population III initial mass function (IMF) is currently unknown, but recent studies agree that fragmentation of primordial gas gives a broader IMF than the initially suggested singular star per halo. In this study we introduce sink particle mergers into Arepo, to perform the first resolution study for primordial star formation simulations and present the first Population III simulations to run up to densities of 10-6g cm-3 for hundreds of years after the formation of sink particles. The total number of sinks formed increases with increasing sink particle creation density, without achieving numerical convergence. The total mass in sinks remains invariant to the maximum resolution and is safe to estimate using low resolution studies. This results in an IMF that shifts towards lower masses with increasing resolution. Greater numbers of sinks cause increased fragmentation-induced starvation of the most massive sink, yielding lower accretion rates, masses and ionising photons emitted per second. The lack of convergence up to densities 2 orders of magnitudes higher than all relevant chemical reactions suggests that the number of sinks will continue to grow with increasing resolution until H2 is fully dissociated and the collapse becomes almost adiabatic at 10-4g cm-3. These results imply that many Population III studies utilising sink particles have produced IMFs which have overestimated the masses of primordial stars, and underestimated the number of stars formed. In the highest resolution runs, sinks with masses capable of surviving until the present day had an ejection fraction of 0.21.

preprint2021arXiv

Giant Molecular Cloud Catalogues for PHANGS-ALMA: Methods and Initial Results

We present improved methods for segmenting CO emission from galaxies into individual molecular clouds, providing an update to the CPROPS algorithms presented by Rosolowsky & Leroy (2006; arXiv:astro-ph/0601706 ). The new code enables both homogenization of the noise and spatial resolution among data, which allows for rigorous comparative analysis. The code also models the completeness of the data via false source injection and includes an updated segmentation approach to better deal with blended emission. These improved algorithms are implemented in a publicly available python package, PYCPROPS. We apply these methods to ten of the nearest galaxies in the PHANGS-ALMA survey, cataloguing CO emission at a common 90 pc resolution and a matched noise level. We measure the properties of 4986 individual clouds identified in these targets. We investigate the scaling relations among cloud properties and the cloud mass distributions in each galaxy. The physical properties of clouds vary among galaxies, both as a function of galactocentric radius and as a function of dynamical environment. Overall, the clouds in our target galaxies are well-described by approximate energy equipartition, although clouds in stellar bars and galaxy centres show elevated line widths and virial parameters. The mass distribution of clouds in spiral arms has a typical mass scale that is 2.5x larger than interarm clouds and spiral arms clouds show slightly lower median virial parameters compared to interarm clouds (1.2 versus 1.4).

preprint2021arXiv

Low-J CO Line Ratios From Single Dish CO Mapping Surveys and PHANGS-ALMA

We measure the low-J CO line ratio R21=CO(2-1)/CO(1-0), R32=CO(3-2)/CO(2-1), and R31 = CO(3-2)/CO(1-0) using whole-disk CO maps of nearby galaxies. We draw CO(2-1) from PHANGS--ALMA, HERACLES, and follow-up IRAM surveys; CO(1-0) from COMING and the Nobeyama CO Atlas of Nearby Spiral Galaxies; and CO(3-2) from the JCMT NGLS and APEX LASMA mapping. Altogether this yields 76, 47, and 29 maps of R21, R32, and R31 at 20&#34; \sim 1.3 kpc resolution, covering 43, 34, and 20 galaxies. Disk galaxies with high stellar mass, log10 M_* [Msun]=10.25-11 and star formation rate, SFR=1-5 Msun/yr, dominate the sample. We find galaxy-integrated mean values and 16%-84% range of R21 = 0.65 (0.50-0.83), R32=0.50 (0.23-0.59), and R31=0.31 (0.20-0.42). We identify weak trends relating galaxy-integrated line ratios to properties expected to correlate with excitation, including SFR/M_* and SFR/L_CO. Within galaxies, we measure central enhancements with respect to the galaxy-averaged value of \sim 0.18^{+0.09}_{-0.14} dex for R21, 0.27^{+0.13}_{-0.15} dex for R31, and 0.08^{+0.11}_{-0.09} dex for R32. All three line ratios anti-correlate with galactocentric radius and positively correlate with the local star formation rate surface density and specific star formation rate, and we provide approximate fits to these relations. The observed ratios can be reasonably reproduced by models with low temperature, moderate opacity, and moderate densities, in good agreement with expectations for the cold ISM. Because the line ratios are expected to anti-correlate with the CO(1-0)-to-H_2 conversion factor, alphaCO^(1-0), these results have general implications for the interpretation of CO emission from galaxies.

preprint2021arXiv

PHANGS-MUSE: the HII region luminosity function of local star-forming galaxies

We use an unprecedented sample of about 23,000 HII regions detected at an average physical resolution of 67pc in the PHANGS-MUSE sample to study the extragalactic HII region Ha luminosity function (LF). Our observations probe the star-forming disk of 19 nearby spiral galaxies with low inclination and located close to the star formation main sequence at z=0. The mean LF slope $α$ in our sample is =1.73 with a $σ$ of 0.15. We find that $α$ decreases with the galaxy&#39;s star formation rate surface density and argue that this is driven by an enhanced clustering of young stars at high gas surface densities. Looking at the HII regions within single galaxies we find that no significant variations occur between the LF of the inner and outer part of the star-forming disk, whereas the LF in the spiral arm areas is shallower than in the inter-arm areas for six out of the 13 galaxies with clearly visible spiral arms. We attribute these variations to the spiral arms increasing the molecular clouds arm--inter-arm mass contrast and find suggestive evidence that they are more evident for galaxies with stronger spiral arms. Furthermore, we find systematic variations in $α$ between samples of HII regions with high and low ionization parameter q and argue that they are driven by the aging of HII regions.

preprint2021arXiv

The PHANGS-HST Survey: Physics at High Angular resolution in Nearby GalaxieS with the Hubble Space Telescope

The PHANGS program is building the first dataset to enable the multi-phase, multi-scale study of star formation across the nearby spiral galaxy population. This effort is enabled by large survey programs with ALMA, VLT/MUSE, and HST, with which we have obtained CO(2-1) imaging, optical spectroscopic mapping, and high resolution UV-optical imaging, respectively. Here, we present PHANGS-HST, which has obtained five band NUV-U-B-V-I imaging of the disks of 38 spiral galaxies at distances of 4-23 Mpc, and parallel V and I band imaging of their halos, to provide a census of tens of thousands of compact star clusters, and multi-scale stellar associations. The combination of HST, ALMA, and VLT/MUSE observations will yield an unprecedented joint catalog of the observed and physical properties of ~100,000 star clusters, associations, HII regions, and molecular clouds. With these basic units of star formation, PHANGS will systematically chart the evolutionary cycling between gas and stars, across a diversity of galactic environments found in nearby galaxies. We discuss the design of the PHANGS-HST survey, and provide an overview of the HST data processing pipeline and first results, highlighting new methods for selecting star cluster candidates, morphological classification of candidates with convolutional neural networks, and identification of stellar associations over a range of physical scales with a watershed algorithm. We describe the cross-observatory imaging, catalogs, and software products to be released. These high-level science products will seed a broad range of investigations, in particular, the study of embedded stellar populations and dust with JWST, for which a PHANGS Cycle 1 Treasury program to obtain eight band 2-21 $μ$m imaging has been approved.

preprint2020arXiv

A Minimum Dilution Scenario for Supernovae and Consequences for Extremely Metal-Poor Stars

To date no metal-free stars have been identified by direct observations. The most common method of constraining their properties is searching the spectra of the most metal-poor stars for the chemical elements created in the first stars and their supernova. In this approach, modelled supernova yields are compared to the observed abundance patterns in extremely metal-poor stars. The method typically only uses the abundance ratios, i.e., the yields are diluted to the observed level. Following the usual assumption of spherical symmetry we compute a simple lower limit of the mass a supernova can mix with and find that it is consistent with all published simulations of early chemical enrichment in the interstellar medium. For three different cases, we demonstrate that this dilution limit can change the conclusions from the abundance fitting. There is a large discrepancy between the dilution found in simulations of SN explosions in minihaloes and the dilution assumed in many abundance fits. Limiting the dilution can significantly alter the likelihood of which supernovae are possible progenitors of observed CEMP-no stars. In particular, some of the faint, very low-yield SNe, which have been suggested as models for the abundance pattern of SMSS0313-6708, cannot explain the measured metal abundances, as their predicted metal yields are too small by two orders of magnitude. Altogether, the new dilution model presented here emphasizes the need to better understand the mixing and dilution behaviour of aspherical SNe.

preprint2020arXiv

A model for the onset of self-gravitation and star formation in molecular gas governed by galactic forces: II. the bottleneck to collapse set by cloud-environment decoupling

In Meidt et al. (2018), we showed that gas kinematics on the scale of individual molecular clouds are not dominated by self-gravity but also track a component that originates with orbital motion in the potential of the host galaxy. This agrees with observed cloud line widths, which show systematic variations from virial motions with environment, pointing at the influence of the galaxy potential. In this paper, we hypothesize that these motions act to slow down the collapse of gas and so help regulate star formation. Extending the results of Meidt et al. (2018), we derive a dynamical collapse timescale that approaches the free-fall time only once the gas has fully decoupled from the galactic potential. Using this timescale we make predictions for how the fraction of free-falling, strongly self-gravitating gas varies throughout the disks of star-forming galaxies. We also use this collapse timescale to predict variations in the molecular gas star formation efficiency, which is lowered from a maximum, feedback-regulated level in the presence of strong coupling to the galactic potential. Our model implies that gas can only decouple from the galaxy to collapse and efficiently form stars deep within clouds. We show that this naturally explains the observed drop in star formation rate per unit gas mass in the Milky Way&#39;s CMZ and other galaxy centers. The model for a galactic bottleneck to star formation also agrees well with resolved observations of dense gas and star formation in galaxy disks and the properties of local clouds.

preprint2020arXiv

A SOFIA Survey of [CII] in the galaxy M51 II. [CII] and CO kinematics across spiral arms

We present the first complete, velocity-resolved [CII] 158um image of the M51 grand-design spiral galaxy, observed with the upGREAT instrument on SOFIA. [CII] is an important tracer of various phases of the interstellar medium (ISM), including ionized gas, neutral atomic, and diffuse molecular regions. We combine the [CII] data with HI, CO, 24um dust continuum, FUV, and near-infrared K-band observations to study the evolution of the ISM across M51&#39;s spiral arms in both position-position, and position-velocity space. Our data show strong velocity gradients in HI, 12CO, and [CII] at the locations of stellar arms (traced by K--band data) with a clear offset in position-velocity space between upstream molecular gas (traced by 12CO) and downstream star formation (traced by [CII]). We compare the observed position--velocity maps across spiral arms with synthetic observations from numerical simulations of galaxies with both dynamical and quasi-stationary steady spiral arms that predict both tangential and radial velocities at the location of spiral arms. We find that our observations, based on the observed velocity gradients and associated offset between CO and [CII], are consistent with the presence of shocks in spiral arms in the inner parts of M51 and in the arm connecting the companion galaxy, M51b, in the outer parts of M51.

preprint2020arXiv

Distances to PHANGS Galaxies: New Tip of the Red Giant Branch Measurements and Adopted Distances

PHANGS-HST is an ultraviolet-optical imaging survey of 38 spiral galaxies within ~20 Mpc. Combined with the PHANGS-ALMA, PHANGS-MUSE surveys and other multiwavelength data, the dataset will provide an unprecedented look into the connections between young stars, HII regions, and cold molecular gas in these nearby star-forming galaxies. Accurate distances are needed to transform measured observables into physical parameters (e.g., brightness to luminosity, angular to physical sizes of molecular clouds, star clusters and associations). PHANGS-HST has obtained parallel ACS imaging of the galaxy halos in the F606W and F814W bands. Where possible, we use these parallel fields to derive tip of the red giant branch (TRGB) distances to these galaxies. In this paper, we present TRGB distances for 11 galaxies from ~4 to ~15 Mpc, based on the first year of PHANGS-HST observations. Five of these represent the first published TRGB distance measurements (IC 5332, NGC 2835, NGC 4298, NGC 4321, and NGC 4328), and eight of which are the best available distances to these targets. We also provide a compilation of distances for the 118 galaxies in the full PHANGS sample, which have been adopted for the first PHANGS-ALMA public data release.

preprint2020arXiv

Dynamical Equilibrium in the Molecular ISM in 28 Nearby Star-Forming Galaxies

We compare the observed turbulent pressure in molecular gas, $P_\mathrm{turb}$, to the required pressure for the interstellar gas to stay in equilibrium in the gravitational potential of a galaxy, $P_\mathrm{DE}$. To do this, we combine arcsecond resolution CO data from PHANGS-ALMA with multi-wavelength data that traces the atomic gas, stellar structure, and star formation rate (SFR) for 28 nearby star-forming galaxies. We find that $P_\mathrm{turb}$ correlates with, but almost always exceeds the estimated $P_\mathrm{DE}$ on kiloparsec scales. This indicates that the molecular gas is over-pressurized relative to the large-scale environment. We show that this over-pressurization can be explained by the clumpy nature of molecular gas; a revised estimate of $P_\mathrm{DE}$ on cloud scales, which accounts for molecular gas self-gravity, external gravity, and ambient pressure, agrees well with the observed $P_\mathrm{turb}$ in galaxy disks. We also find that molecular gas with cloud-scale ${P_\mathrm{turb}}\approx{P_\mathrm{DE}}\gtrsim{10^5\,k_\mathrm{B}\,\mathrm{K\,cm^{-3}}}$ in our sample is more likely to be self-gravitating, whereas gas at lower pressure appears more influenced by ambient pressure and/or external gravity. Furthermore, we show that the ratio between $P_\mathrm{turb}$ and the observed SFR surface density, $Σ_\mathrm{SFR}$, is compatible with stellar feedback-driven momentum injection in most cases, while a subset of the regions may show evidence of turbulence driven by additional sources. The correlation between $Σ_\mathrm{SFR}$ and kpc-scale $P_\mathrm{DE}$ in galaxy disks is consistent with the expectation from self-regulated star formation models. Finally, we confirm the empirical correlation between molecular-to-atomic gas ratio and kpc-scale $P_\mathrm{DE}$ reported in previous works.

preprint2020arXiv

Formation sites of Population III star formation: The effects of different levels of rotation and turbulence on the fragmentation behavior of primordial gas

We use the moving-mesh code AREPO to investigate the effects of different levels of rotation and turbulence on the fragmentation of primordial gas and the formation of Population III stars. We consider 9 different combinations of turbulence and rotation and carry out 5 different realizations of each setup, yielding one of the largest sets of simulations of Population III star formation ever performed. We find that fragmentation in Population III star-forming systems is a highly chaotic process and show that the outcomes of individual realizations of the same initial conditions often vary significantly. However, some general trends are apparent. Increasing the turbulent energy promotes fragmentation, while increasing the rotational energy inhibits fragmentation. Within the 1000 yr period that we simulate, runs including turbulence yield flat protostellar mass functions while purely rotational runs show a more top-heavy distribution. The masses of the individual protostars are distributed over a wide range from a few $10^{-3} \, {\rm M_{\odot}}$ to several tens of ${\rm M_\odot}$. The total mass growth rate of the stellar systems remains high throughout the simulations and depends only weakly on the degree of rotation and turbulence. Mergers between protostars are common, but predictions of the merger fraction are highly sensitive to the criterion used to decide whether two protostars should merge. Previous studies of Population III star formation have often considered only one realization per set of initial conditions. However, our results demonstrate that robust trends can only be reliably identified by considering averages over a larger sample of runs.

preprint2020arXiv

PHANGS CO kinematics: disk orientations and rotation curves at 150 pc resolution

We present kinematic orientations and high resolution (150 pc) rotation curves for 67 main sequence star-forming galaxies surveyed in CO (2-1) emission by PHANGS-ALMA. Our measurements are based on the application of a new fitting method tailored to CO velocity fields. Our approach identifies an optimal global orientation as a way to reduce the impact of non-axisymmetric (bar and spiral) features and the uneven spatial sampling characteristic of CO emission in the inner regions of nearby galaxies. The method performs especially well when applied to the large number of independent lines-of-sight contained in the PHANGS CO velocity fields mapped at 1&#39;&#39; resolution. The high resolution rotation curves fitted to these data are sensitive probes of mass distribution in the inner regions of these galaxies. We use the inner slope as well as the amplitude of our fitted rotation curves to demonstrate that CO is a reliable global dynamical mass tracer. From the consistency between photometric orientations from the literature and kinematic orientations determined with our method, we infer that the shapes of stellar disks in the mass range of log($\rm M_{\star}(M_{\odot})$)=9.0-10.9 probed by our sample are very close to circular and have uniform thickness.

preprint2020arXiv

Physical Processes in Star Formation

Star formation is a complex multi-scale phenomenon that is of significant importance for astrophysics in general. Stars and star formation are key pillars in observational astronomy from local star forming regions in the Milky Way up to high-redshift galaxies. From a theoretical perspective, star formation and feedback processes (radiation, winds, and supernovae) play a pivotal role in advancing our understanding of the physical processes at work, both individually and of their interactions. In this review we will give an overview of the main processes that are important for the understanding of star formation. We start with an observationally motivated view on star formation from a global perspective and outline the general paradigm of the life-cycle of molecular clouds, in which star formation is the key process to close the cycle. After that we focus on the thermal and chemical aspects in star forming regions, discuss turbulence and magnetic fields as well as gravitational forces. Finally, we review the most important stellar feedback mechanisms.

preprint2020arXiv

The role of galactic dynamics in shaping the physical properties of giant molecular clouds in Milky Way-like galaxies

We examine the role of the large-scale galactic-dynamical environment in setting the properties of giant molecular clouds in Milky Way-like galaxies. We perform three high-resolution simulations of Milky Way-like discs with the moving-mesh hydrodynamics code Arepo, yielding a statistical sample of $\sim 80,000$ giant molecular clouds and $\sim 55,000$ HI clouds. We account for the self-gravity of the gas, momentum and thermal energy injection from supernovae and HII regions, mass injection from stellar winds, and the non-equilibrium chemistry of hydrogen, carbon and oxygen. By varying the external gravitational potential, we probe galactic-dynamical environments spanning an order of magnitude in the orbital angular velocity, gravitational stability, mid-plane pressure and the gradient of the galactic rotation curve. The simulated molecular clouds are highly overdense ($\sim 100\times$) and over-pressured ($\sim 25\times$) relative to the ambient interstellar medium. Their gravo-turbulent and star-forming properties are decoupled from the dynamics of the galactic mid-plane, so that the kpc-scale star formation rate surface density is related only to the number of molecular clouds per unit area of the galactic mid-plane. Despite this, the clouds display clear, statistically-significant correlations of their rotational properties with the rates of galactic shearing and gravitational free-fall. We find that galactic rotation and gravitational instability can influence their elongation, angular momenta, and tangential velocity dispersions. The lower pressures and densities of the HI clouds allow for a greater range of significant dynamical correlations, mirroring the rotational properties of the molecular clouds, while also displaying a coupling of their gravitational and turbulent properties to the galactic-dynamical environment.

preprint2020arXiv

When Gas Dynamics Decouples from Galactic Rotation: Characterizing ISM Circulation in Disk Galaxies

In galactic disks, galactic rotation sets the bulk motion of gas, and its energy and momentum can be transferred toward small scales. Additionally, in the interstellar medium, random and noncircular motions arise from stellar feedback, cloud-cloud interactions, and instabilities, among other processes. Our aim is to comprehend to which extent small-scale gas dynamics is decoupled from galactic rotation. We study the relative contributions of galactic rotation and local noncircular motions to the circulation of gas, $Γ$, a macroscopic measure of local rotation, defined as the line integral of the velocity field around a closed path. We measure the circulation distribution as a function of spatial scale in a set of simulated disk galaxies and we model the velocity field as the sum of galactic rotation and a Gaussian random field. The random field is parameterized by a broken power law in Fourier space, with a break at the scale $λ_c$. We define the spatial scale $λ_{\rm eq}$ at which galactic rotation and non-circular motions contribute equally to $Γ$. For our simulated galaxies, the gas dynamics at the scale of molecular clouds is usually dominated by noncircular motions, but in the center of galactic disks galactic rotation is still relevant. Our model shows that the transfer of rotation from large scales breaks at the scale $λ_c$ and this transition is necessary to reproduce the circulation distribution. We find that $λ_{\rm eq}$, and therefore the structure of the gas velocity field, is set by the local conditions of gravitational stability and stellar feedback.

preprint2019arXiv

Binary outflows from young stars: interaction of co-orbital jet and wind

Jets from young stellar objects provide insight into the workings of the beating heart at the centre of star forming cores. In some cases, multiple pulsed outflows are detected such as the atomic and molecular jets from a proposed binary system in the T\,Tauri star HH30. We investigate here the development and propagation of duelling atomic and molecular outflows stemming from the two stars in co-orbit. We perform a series of numerical experiments with the {\small ZEUS-MP} code with enhanced cooling and chemistry modules. The aim of this work is to identify signatures on scales of order 100 AU. The jet sources are off the grid domain and so it is the propagation and interaction from ~ 20AU out to 100\,AU simulated here. We find that the molecular flow from the orbiting source significantly disturbs the atomic jet, deflecting and twisting the jet and disrupting the jet knots. Regions of high ionisation are generated as the atomic jet rams through the dense molecular outflow. Synthetic images in atomic and molecular lines are presented which demonstrate identifying signatures. In particular, the structure within the atomic jet is lost and H-alpha may trace the walls of the present CO cavity or where the walls have been recently. These results provide a framework for the interpretation of upcoming high resolution observations.

preprint2019arXiv

Mapping electron temperature variations across a spiral arm in NGC 1672

We report one of the first extragalactic observations of electron temperature variations across a spiral arm. Using MUSE mosaic observations of the nearby galaxy NGC 1672, we measure the [N II]5755 auroral line in a sample of 80 HII regions in the eastern spiral arm of NGC1672. We discover systematic temperature variations as a function of distance perpendicular to the spiral arm. The electron temperature is lowest on the spiral arm itself and highest on the downstream side. Photoionization models of different metallicity, pressure, and age of the ionizing source are explored to understand what properties of the interstellar medium drive the observed temperature variations. An azimuthally varying metallicity appears to be the most likely cause of the temperature variations. The electron temperature measurements solidify recent discoveries of azimuthal variations of oxygen abundance based on strong lines, and rule out the possibility that the abundance variations are artefacts of the strong-line calibrations.

preprint2019arXiv

Relations Between Molecular Cloud Structure Sizes and Line Widths in the Large Magellanic Cloud

We present a comparative study of the size-line width relation for substructures within six molecular clouds in the Large Magellanic Cloud (LMC) mapped with the Atacama Large Millimeter/submillimeter Array (ALMA). Our sample extends our previous study, which compared a Planck detected cold cloud in the outskirts of the LMC with the 30 Doradus molecular cloud and found the typical line width for 1 pc radius structures to be 5 times larger in 30 Doradus. By observing clouds with intermediate levels of star formation activity, we find evidence that line width at a given size increases with increasing local and cloud-scale 8$μ$m intensity. At the same time, line width at a given size appears to independently correlate with measures of mass surface density. Our results suggest that both virial-like motions due to gravity and local energy injection by star formation feedback play important roles in determining intracloud dynamics.

preprint2019arXiv

Simulations of the star-forming molecular gas in an interacting M51-like galaxy

We present here the first of a series of papers aimed at better understanding the evolution and properties of giant molecular clouds (GMCs) in a galactic context. We perform high resolution, three-dimensional {\sc arepo} simulations of an interacting galaxy inspired by the well-observed M51 galaxy. Our fiducial simulations include a non-equilibrium, time-dependent, chemical network that follows the evolution of atomic and molecular hydrogen as well as carbon and oxygen self-consistently. Our calculations also treat gas self-gravity and subsequent star formation (described by sink particles), and coupled supernova feedback. In the densest parts of the simulated interstellar medium (ISM) we reach sub-parsec resolution, granting us the ability to resolve individual GMCs and their formation and destruction self-consistently throughout the galaxy. In this initial work we focus on the general properties of the ISM with a particular focus on the cold star-forming gas. We discuss the role of the interaction with the companion galaxy in generating cold molecular gas and controlling stellar birth. We find that while the interaction drives large-scale gas flows and induces spiral arms in the galaxy, it is of secondary importance in determining gas fractions in the different ISM phases and the overall star-formation rate. The behaviour of the gas on small GMC scales instead is mostly controlled by the self-regulating property of the ISM driven by coupled feedback.

preprint2019arXiv

The Cloud Factory I: Generating resolved filamentary molecular clouds from galactic-scale forces

We introduce a new suite of simulations, &#34;The Cloud Factory&#34;, which self-consistently forms molecular cloud complexes at high enough resolution to resolve internal substructure (up to 0.25 Msol in mass) all while including galactic-scale forces. We use a version of the Arepo code modified to include a detailed treatment of the physics of the cold molecular ISM, and an analytical galactic gravitational potential for computational efficiency. The simulations have nested levels of resolution, with the lowest layer tied to tracer particles injected into individual cloud complexes. These tracer refinement regions are embedded in the larger simulation so continue to experience forces from outside the cloud. This allows the simulations to act as a laboratory for testing the effect of galactic environment on star formation. Here we introduce our method and investigate the effect of galactic environment on filamentary clouds. We find that cloud complexes formed after a clustered burst of feedback, have shorter lengths and are less likely to fragment compared to quiescent clouds (e.g. the Musca filament) or those dominated by the galactic potential (e.g. Nessie). Spiral arms and differential rotation preferentially align filaments, but strong feedback randomises them. Long filaments formed within the cloud complexes are necessarily coherent with low internal velocity gradients, which has implications for the formation of filamentary star-clusters. Cloud complexes formed in regions dominated by supernova feedback have fewer star-forming cores, and these are more widely distributed. These differences show galactic-scale forces can have a significant impact on star formation within molecular clouds.

preprint2019arXiv

The headlight cloud in NGC 628: An extreme giant molecular cloud in a typical galaxy disk

Cloud-scale surveys of molecular gas reveal the link between molecular clouds properties and star formation (SF) across a range of galactic environments. Cloud populations in galaxy disks are considered to be representative of the `normal&#39; SF. At high resolution, however, clouds with exceptional gas properties and SF activity may also be observed in normal disk environments. In this paper, we study the brightest cloud traced in CO emission in the disk of NGC628. The cloud is spatially coincident with an extremely bright HII region. We characterize its molecular gas properties and investigate how feedback and large-scale processes influence the properties of the molecular gas. High resolution CO ALMA observations are used to characterize its mass and dynamical state, which are compared to other clouds in NGC628. A LVG analysis is used to constrain the beam-diluted density and temperature of the molecular gas. We analyze the MUSE spectrum using Starburst99 to characterize the young stellar population associated with the HII region. The cloud is massive ($1-2\times10^7$M$_{\odot}$), with a beam-diluted density of $n_{\rm H_2}=5\times10^4$ cm$^{-3}$. It has a low virial parameter, suggesting that its CO emission may be overluminous due to heating by the HII region. A young ($2-4$ Myr), massive $3\times10^{5}$ M$_{\odot}$ stellar population is associated. We argue that the cloud is currently being destroyed by feedback from young massive stars. Due to the cloud&#39;s large mass, this phase of the cloud&#39;s evolution is long enough for the impact of feedback on the excitation of the gas to be observed. Its high mass may be related to its location at a spiral co-rotation radius, where gas experiences reduced galactic shear compared to other regions of the disk, and receives a sustained inflow of gas that can promote the cloud&#39;s mass growth.

preprint2019arXiv

The lifecycle of molecular clouds in nearby star-forming disc galaxies

It remains a major challenge to derive a theory of cloud-scale ($\lesssim100$ pc) star formation and feedback, describing how galaxies convert gas into stars as a function of the galactic environment. Progress has been hampered by a lack of robust empirical constraints on the giant molecular cloud (GMC) lifecycle. We address this problem by systematically applying a new statistical method for measuring the evolutionary timeline of the GMC lifecycle, star formation, and feedback to a sample of nine nearby disc galaxies, observed as part of the PHANGS-ALMA survey. We measure the spatially-resolved ($\sim100$ pc) CO-to-H$α$ flux ratio and find a universal de-correlation between molecular gas and young stars on GMC scales, allowing us to quantify the underlying evolutionary timeline. GMC lifetimes are short, typically 10-30 Myr, and exhibit environmental variation, between and within galaxies. At kpc-scale molecular gas surface densities $Σ_{\rm H_2}\geqslant8$M$_{\odot}$pc$^{-2}$, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at $Σ_{\rm H_2}\leqslant8$M$_{\odot}$pc$^{-2}$ GMCs decouple from galactic dynamics and live for an internal dynamical time-scale. After a long inert phase without massive star formation traced by H$α$ (75-90% of the cloud lifetime), GMCs disperse within just 1-5 Myr once massive stars emerge. The dispersal is most likely due to early stellar feedback, causing GMCs to achieve integrated star formation efficiencies of 4-10% These results show that galactic star formation is governed by cloud-scale, environmentally-dependent, dynamical processes driving rapid evolutionary cycling. GMCs and HII regions are the fundamental units undergoing these lifecycles, with mean separations of 100-300 pc in star-forming discs. Future work should characterise the multi-scale physics and mass flows driving these lifecycles.

preprint2019arXiv

WARPFIELD-EMP: The Self-Consistent Prediction of Emission Lines from Evolving HII Regions in Dense Molecular Clouds

We present the {\sc warpfield} emission predictor, {\sc warpfield-emp}, which couples the 1D stellar feedback code {\sc warpfield} with the {\sc cloudy} \hii region/PDR code and the {\sc polaris} radiative transfer code, in order to make detailed predictions for the time-dependent line and continuum emission arising from the H{\sc ii} region and PDR surrounding an evolving star cluster. {\sc warpfield-emp} accounts for a wide range of physical processes (stellar winds, supernovae, radiation pressure, gravity, thermal conduction, radiative cooling, dust extinction etc.) and yet runs quickly enough to allow us to explore broad ranges of different cloud parameters. We compare the results of an extensive set of models with SITELLE observations of a large sample of \hii regions in NGC~628 and find very good agreement, particularly for the highest signal-to-noise observations. We show that our approach of modeling individual clouds from first principles (instead of in terms of dimensionless quantities such as the ionization parameter) allows us to avoid long-standing degeneracies in the interpretation of \hii region diagnostics and enables us to relate these diagnostics to important physical parameters such as cloud mass or cluster age. Finally, we explore the implications of our models regarding the reliability of simple metallicity diagnostics, the properties of long-lived embedded clusters, and the role played by winds and supernovae in regulating \hii region and PDR line emission.

preprint2018arXiv

Efficacy of early stellar feedback in low gas surface density environments

We present a suite of high resolution radiation hydrodynamic simulations of a small patch ($1 \ {\rm kpc}^2$) of the inter-stellar medium (ISM) performed with Arepo-RT, with the aim to quantify the efficacy of various feedback processes like supernovae explosions (SNe), photoheating and radiation pressure in low gas surface density galaxies ($Σ_{\rm gas} \simeq 10 \ {\rm M}_\odot \ {\rm pc}^{-2}$). We show that radiation fields decrease the star formation rate and therefore the total stellar mass formed by a factor of $\sim 2$. This increases the gas depletion timescale and brings the simulated Kennicutt-Schmidt relation closer to the observational estimates. Radiation feedback coupled with SNe is more efficient at driving outflows with the mass and energy loading increasing by a factor of $\sim 10$. This increase is mainly driven by the additional entrainment of medium density ($10^{-2} \leq n< 1 \ {\rm cm}^{-3}$), warm ($300 \ {\rm K}\leq T<8000 \ {\rm K}$) material. Therefore including radiation fields tends to launch colder, denser and higher mass and energy loaded outflows. This is because photoheating of the high density gas around a newly formed star over-pressurises the region, causing it to expand. This reduces the ambient density in which the SNe explode by a factor of $10-100$ which in turn increases their momentum output by a factor of $\sim 1.5-2.5$. Finally, we note that in these low gas surface density environments, radiation fields primarily impact the ISM via photoheating and radiation pressure has only a minimal role in regulating star formation.