Search

Galaxy cluster outskirts mark the transition region from the mildly non-linear cosmic web to the highly non-linear, virialised, cluster interior. It is in this transition region that the intra-cluster medium (ICM) begins to influence the properties of accreting galaxies and groups, as ram pressure impacts a galaxy's cold gas content and subsequent star formation rate. Conversely, the thermodynamical properties of the ICM in this transition region should also feel the influence of accreting substructure (i.e. galaxies and groups), whose passage can drive shocks. In this paper, we use a suite of cosmological hydrodynamical zoom simulations of a single galaxy cluster, drawn from the nIFTy comparison project, to study how the dynamics of substructure accreted from the cosmic web influences the thermodynamical properties of the ICM in the cluster's outskirts. We demonstrate how features evident in radial profiles of the ICM (e.g. gas density and temperature) can be linked to strong shocks, transient and short-lived in nature, driven by the passage of substructure. The range of astrophysical codes and galaxy formation models in our comparison are broadly consistent in their predictions (e.g. agreeing when and where shocks occur, but differing in how strong shocks will be); this is as we would expect of a process driven by large-scale gravitational dynamics and strong, inefficiently radiating, shocks. This suggests that mapping such shock structures in the ICM in a cluster's outskirts (via e.g. radio synchrotron emission) could provide a complementary measure of its recent merger and accretion history., Comment: 15 pages, 13 figures, version accepted for publication in MNRAS

Galaxy cluster outskirts mark the transition region from the mildly non-linear cosmic web to the highly non-linear, virialised, cluster interior. It is in this transition region that the intra-cluster medium (ICM) begins to influence the properties of accreting galaxies and groups, as ram pressure impacts a galaxy's cold gas content and subsequent star formation rate. Conversely, the thermodynamical properties of the ICM in this transition region should also feel the influence of accreting substructure (i.e. galaxies and groups), whose passage can drive shocks. In this paper, we use a suite of cosmological hydrodynamical zoom simulations of a single galaxy cluster, drawn from the nIFTy comparison project, to study how the dynamics of substructure accreted from the cosmic web influences the thermodynamical properties of the ICM in the cluster's outskirts. We demonstrate how features evident in radial profiles of the ICM (e.g. gas density and temperature) can be linked to strong shocks, transient and short-lived in nature, driven by the passage of substructure. The range of astrophysical codes and galaxy formation models in our comparison are broadly consistent in their predictions (e.g. agreeing when and where shocks occur, but differing in how strong shocks will be); this is as we would expect of a process driven by large-scale gravitational dynamics and strong, inefficiently radiating, shocks. This suggests that mapping such shock structures in the ICM in a cluster's outskirts (via e.g. radio synchrotron emission) could provide a complementary measure of its recent merger and accretion history., Comment: 15 pages, 13 figures, version accepted for publication in MNRAS

There are many proposed mechanisms driving themorphological transformation of disc galaxies to elliptical galaxies. In this paper,we determine if the observed transformation in low-mass groups can be explained by the merger histories of galaxies. We measured the group mass- morphology relation for groups from the Galaxy and Mass Assembly group catalogue with masses from 1011 to 1015 M⊙. Contrary to previous studies, the fraction of elliptical galaxies in our more complete group sample increases significantly with group mass across the full range of group mass. The elliptical fraction increases at a rate of 0.163 ± 0.012 per dex of group mass for groups more massive than 1012.5 M⊙. If we allow for uncertainties in the observed group masses, our results are consistent with a continuous increase in elliptical fraction from group masses as low as 1011 M⊙. We tested if this observed relation is consistent with the merger activity using a GADGET-2 dark matter simulation of the galaxy groups. We specified that a simulated galaxy would be transformed to an elliptical morphology either if it experienced a major merger or if its cumulative mass gained from minor mergers exceeded 30 per cent of its final mass. We then calculated a group mass-morphology relation for the simulations. The position and slope of the simulated relation were consistent with the observational relation, with a gradient of 0.184 ± 0.010 per dex of group mass. These results demonstrate a strong correlation between the frequency of merger events and disc-to-elliptical galaxy transformation in galaxy group environments.

We present an analysis of the neutral hydrogen (H i) properties of a fully cosmological hydrodynamical dwarf galaxy, run with varying simulation parameters. As reported by Governato et al., the high-resolution, high star formation density threshold version of this galaxy is the first simulation to result in the successful reproduction of a (dwarf) spiral galaxy without any associated stellar bulge. We have set out to compare in detail the H i distribution and kinematics of this simulated bulgeless disc with what is observed in a sample of nearby dwarfs. To do so, we extracted the radial gas density profiles, velocity dispersion (e.g. velocity ellipsoid and turbulence) and the power spectrum of structure within the cold interstellar medium (ISM) from the simulations. The highest resolution dwarf, when using a high-density star formation threshold comparable to densities of giant molecular clouds, possesses bulk characteristics consistent with those observed in nature, though the cold gas is not as radially extended as that observed in nearby dwarfs, resulting in somewhat excessive surface densities. The lines-of-sight velocity dispersion radial profiles have values that are in good agreement with the observed dwarf galaxies, but due to the fact that only the streaming velocities of particles are tracked, a correction to include the thermal velocities can lead to profiles that are quite flat. The ISM power spectra of the simulations appear to possess more power on smaller spatial scales than that of the Small Magellanic Cloud. We conclude that unavoidable limitations remain due to the unresolved physics of star formation and feedback within parsec-scale molecular clouds

Theoretical galaxy formation models are an established and powerful tool for interpreting the astrophysical significance of observational data, particularly galaxy surveys. Such models have been utilised with great success by optical surveys such as 2dFGRS and SDSS, but their application to radio surveys of cold gas in galaxies has been limited. In this chapter we describe recent developments in the modelling of the cold gas properties in the models, and how these developments are essential if they are to be applied to cold gas surveys of the kind that will be carried out with the SKA. By linking explicitly a galaxy's star formation rate to the abundance of molecular hydrogen in the galaxy rather than cold gas abundance, as was assumed previously, the latest models reproduce naturally many of the global atomic and molecular hydrogen properties of observed galaxies. We review some of the key results of the latest models and highlight areas where further developments are necessary. We discuss also how model predictions can be most accurately compared with observational data, what challenges we expect when creating synthetic galaxy surveys in the SKA era, and how the SKA can be used to test models of dark matter., Comment: 20 pages, 9 figures. Contribution to the conference 'Advancing Astrophysics with the Square Kilometre Array', June 8-13, 2014, Giardini Naxos, Italy

Theoretical galaxy formation models are an established and powerful tool for interpreting the astrophysical significance of observational data, particularly galaxy surveys. Such models have been utilised with great success by optical surveys such as 2dFGRS and SDSS, but their application to radio surveys of cold gas in galaxies has been limited. In this chapter we describe recent developments in the modelling of the cold gas properties in the models, and how these developments are essential if they are to be applied to cold gas surveys of the kind that will be carried out with the SKA. By linking explicitly a galaxy's star formation rate to the abundance of molecular hydrogen in the galaxy rather than cold gas abundance, as was assumed previously, the latest models reproduce naturally many of the global atomic and molecular hydrogen properties of observed galaxies. We review some of the key results of the latest models and highlight areas where further developments are necessary. We discuss also how model predictions can be most accurately compared with observational data, what challenges we expect when creating synthetic galaxy surveys in the SKA era, and how the SKA can be used to test models of dark matter., Comment: 20 pages, 9 figures. Contribution to the conference 'Advancing Astrophysics with the Square Kilometre Array', June 8-13, 2014, Giardini Naxos, Italy

We use the Galaxies-Intergalactic Medium Interaction Calculation (gimic) cosmological hydrodynamic simulation at z = 0 to study the distribution and environmental dependence of neutral hydrogen (H i) gas in the outskirts of simulated galaxies. This gas can currently be probed directly in, for example, Ly alpha absorption via the observation of background quasars. Radio facilities, such as the Square Kilometre Array, will provide a complementary probe of the diffuse H i in emission and will constrain the physics underpinning the complex interplay between accretion and feedback mechanisms which affect the intergalactic medium. We extract a sample of 488 galaxies from a resimulation of the average cosmic density gimic region. We estimate the neutral hydrogen content of these galaxies and the surrounding intergalactic medium within which they reside. We investigate the average H i radial profiles by stacking the individual profiles according to both mass and environment. We find high H i column densities at large impact parameters in group environments and markedly lower H i densities for non-group galaxies. We suggest that these results likely arise from the combined effects of ram pressure stripping and tidal interactions present in group environments.

We use the Galaxies-Intergalactic Medium Interaction Calculation (gimic) cosmological hydrodynamic simulation at z = 0 to study the distribution and environmental dependence of neutral hydrogen (H i) gas in the outskirts of simulated galaxies. This gas can currently be probed directly in, for example, Ly alpha absorption via the observation of background quasars. Radio facilities, such as the Square Kilometre Array, will provide a complementary probe of the diffuse H i in emission and will constrain the physics underpinning the complex interplay between accretion and feedback mechanisms which affect the intergalactic medium. We extract a sample of 488 galaxies from a resimulation of the average cosmic density gimic region. We estimate the neutral hydrogen content of these galaxies and the surrounding intergalactic medium within which they reside. We investigate the average H i radial profiles by stacking the individual profiles according to both mass and environment. We find high H i column densities at large impact parameters in group environments and markedly lower H i densities for non-group galaxies. We suggest that these results likely arise from the combined effects of ram pressure stripping and tidal interactions present in group environments.

We use the Galaxies-Intergalactic Medium Interaction Calculation (gimic) cosmological hydrodynamic simulation at z = 0 to study the distribution and environmental dependence of neutral hydrogen (H i) gas in the outskirts of simulated galaxies. This gas can currently be probed directly in, for example, Ly alpha absorption via the observation of background quasars. Radio facilities, such as the Square Kilometre Array, will provide a complementary probe of the diffuse H i in emission and will constrain the physics underpinning the complex interplay between accretion and feedback mechanisms which affect the intergalactic medium. We extract a sample of 488 galaxies from a resimulation of the average cosmic density gimic region. We estimate the neutral hydrogen content of these galaxies and the surrounding intergalactic medium within which they reside. We investigate the average H i radial profiles by stacking the individual profiles according to both mass and environment. We find high H i column densities at large impact parameters in group environments and markedly lower H i densities for non-group galaxies. We suggest that these results likely arise from the combined effects of ram pressure stripping and tidal interactions present in group environments.

We use the Galaxies-Intergalactic Medium Interaction Calculation (gimic) cosmological hydrodynamic simulation at z = 0 to study the distribution and environmental dependence of neutral hydrogen (H i) gas in the outskirts of simulated galaxies. This gas can currently be probed directly in, for example, Ly alpha absorption via the observation of background quasars. Radio facilities, such as the Square Kilometre Array, will provide a complementary probe of the diffuse H i in emission and will constrain the physics underpinning the complex interplay between accretion and feedback mechanisms which affect the intergalactic medium. We extract a sample of 488 galaxies from a resimulation of the average cosmic density gimic region. We estimate the neutral hydrogen content of these galaxies and the surrounding intergalactic medium within which they reside. We investigate the average H i radial profiles by stacking the individual profiles according to both mass and environment. We find high H i column densities at large impact parameters in group environments and markedly lower H i densities for non-group galaxies. We suggest that these results likely arise from the combined effects of ram pressure stripping and tidal interactions present in group environments.

We use the Galaxies-Intergalactic Medium Interaction Calculation (gimic) cosmological hydrodynamic simulation at z = 0 to study the distribution and environmental dependence of neutral hydrogen (H i) gas in the outskirts of simulated galaxies. This gas can currently be probed directly in, for example, Ly alpha absorption via the observation of background quasars. Radio facilities, such as the Square Kilometre Array, will provide a complementary probe of the diffuse H i in emission and will constrain the physics underpinning the complex interplay between accretion and feedback mechanisms which affect the intergalactic medium. We extract a sample of 488 galaxies from a resimulation of the average cosmic density gimic region. We estimate the neutral hydrogen content of these galaxies and the surrounding intergalactic medium within which they reside. We investigate the average H i radial profiles by stacking the individual profiles according to both mass and environment. We find high H i column densities at large impact parameters in group environments and markedly lower H i densities for non-group galaxies. We suggest that these results likely arise from the combined effects of ram pressure stripping and tidal interactions present in group environments.

We use the Galaxies-Intergalactic Medium Interaction Calculation (gimic) cosmological hydrodynamic simulation at z = 0 to study the distribution and environmental dependence of neutral hydrogen (H i) gas in the outskirts of simulated galaxies. This gas can currently be probed directly in, for example, Ly alpha absorption via the observation of background quasars. Radio facilities, such as the Square Kilometre Array, will provide a complementary probe of the diffuse H i in emission and will constrain the physics underpinning the complex interplay between accretion and feedback mechanisms which affect the intergalactic medium. We extract a sample of 488 galaxies from a resimulation of the average cosmic density gimic region. We estimate the neutral hydrogen content of these galaxies and the surrounding intergalactic medium within which they reside. We investigate the average H i radial profiles by stacking the individual profiles according to both mass and environment. We find high H i column densities at large impact parameters in group environments and markedly lower H i densities for non-group galaxies. We suggest that these results likely arise from the combined effects of ram pressure stripping and tidal interactions present in group environments.

We use the Galaxies-Intergalactic Medium Interaction Calculation (gimic) cosmological hydrodynamic simulation at z = 0 to study the distribution and environmental dependence of neutral hydrogen (H i) gas in the outskirts of simulated galaxies. This gas can currently be probed directly in, for example, Ly alpha absorption via the observation of background quasars. Radio facilities, such as the Square Kilometre Array, will provide a complementary probe of the diffuse H i in emission and will constrain the physics underpinning the complex interplay between accretion and feedback mechanisms which affect the intergalactic medium. We extract a sample of 488 galaxies from a resimulation of the average cosmic density gimic region. We estimate the neutral hydrogen content of these galaxies and the surrounding intergalactic medium within which they reside. We investigate the average H i radial profiles by stacking the individual profiles according to both mass and environment. We find high H i column densities at large impact parameters in group environments and markedly lower H i densities for non-group galaxies. We suggest that these results likely arise from the combined effects of ram pressure stripping and tidal interactions present in group environments.

We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) cosmological hydrodynamic simulation at z=0 to study the distribution and environmental dependence of neutral hydrogen (HI) gas in the outskirts of simulated galaxies. This gas can currently be probed directly in, for example, Ly$\alpha$ absorption via the observation of background quasars. Radio facilities, such as the Square Kilometre Array, will provide a complementary probe of the diffuse HI in emission and will constrain the physics underpinning the complex interplay between accretion and feedback mechanisms which affect the intergalactic medium. We extract a sample of 488 galaxies from a re-simulation of the average cosmic density GIMIC region. We estimate the neutral hydrogen content of these galaxies and the surrounding intergalactic medium within which they reside. We investigate the average HI radial profiles by stacking the individual profiles according to both mass and environment. We find high HI column densities at large impact parameters in group environments and markedly lower HI densities for non-group galaxies. We suggest that these results likely arise from the combined effects of ram pressure stripping and tidal interactions present in group environments., Comment: 9 pages, 3 figures, MNRAS accepted

We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) cosmological hydrodynamic simulation at z=0 to study the distribution and environmental dependence of neutral hydrogen (HI) gas in the outskirts of simulated galaxies. This gas can currently be probed directly in, for example, Ly$\alpha$ absorption via the observation of background quasars. Radio facilities, such as the Square Kilometre Array, will provide a complementary probe of the diffuse HI in emission and will constrain the physics underpinning the complex interplay between accretion and feedback mechanisms which affect the intergalactic medium. We extract a sample of 488 galaxies from a re-simulation of the average cosmic density GIMIC region. We estimate the neutral hydrogen content of these galaxies and the surrounding intergalactic medium within which they reside. We investigate the average HI radial profiles by stacking the individual profiles according to both mass and environment. We find high HI column densities at large impact parameters in group environments and markedly lower HI densities for non-group galaxies. We suggest that these results likely arise from the combined effects of ram pressure stripping and tidal interactions present in group environments., Comment: 9 pages, 3 figures, MNRAS accepted

We highlight two research strands related to our ongoing chemodynamical Galactic Archaeology efforts: (i) the spatio-temporal infall rate of gas onto the disk, drawing analogies with the infall behaviour imposed by classical galactic chemical evolution models of inside-out disk growth; (ii) the radial age gradient predicted by spectrophometric models of disk galaxies. In relation to (i), at low-redshift, we find that half of the infall onto the disk is gas associated with the corona, while half can be associated with cooler gas streams; we also find that gas enters the disk preferentially orthogonal to the system, rather than in-plane. In relation to (ii), we recover age gradient troughs/inflections consistent with those observed in nature, without recourse to radial migrations., Comment: To appear in "From the First Structures to the Universe Today"; Proceedings of 2nd Workshop on Numerical and Observational Astrophysics; Buenos Aires; 14-18 Nov 2011; ed. M.W. De Rossi. S.E. Pedrosa & L.J. Pellizza; AAA Workshop Series (Argentine Astronomical Society)

We highlight two research strands related to our ongoing chemodynamical Galactic Archaeology efforts: (i) the spatio-temporal infall rate of gas onto the disk, drawing analogies with the infall behaviour imposed by classical galactic chemical evolution models of inside-out disk growth; (ii) the radial age gradient predicted by spectrophometric models of disk galaxies. In relation to (i), at low-redshift, we find that half of the infall onto the disk is gas associated with the corona, while half can be associated with cooler gas streams; we also find that gas enters the disk preferentially orthogonal to the system, rather than in-plane. In relation to (ii), we recover age gradient troughs/inflections consistent with those observed in nature, without recourse to radial migrations., Comment: To appear in "From the First Structures to the Universe Today"; Proceedings of 2nd Workshop on Numerical and Observational Astrophysics; Buenos Aires; 14-18 Nov 2011; ed. M.W. De Rossi. S.E. Pedrosa & L.J. Pellizza; AAA Workshop Series (Argentine Astronomical Society)

We present an analysis of the neutral hydrogen (HI) properties of a fully cosmological hydrodynamical dwarf galaxy, run with varying simulation parameters. As reported by Governato et al. (2010), the high resolution, high star formation density threshold version of this galaxy is the first simulation to result in the successful reproduction of a (dwarf) spiral galaxy without any associated stellar bulge. We have set out to compare in detail the HI distribution and kinematics of this simulated bulgeless disk with what is observed in a sample of nearby dwarfs. To do so, we extracted the radial gas density profiles, velocity dispersion (e.g., velocity ellipsoid, turbulence), and the power spectrum of structure within the cold interstellar medium from the simulations. The highest resolution dwarf, when using a high density star formation threshold comparable to densities of giant molecular clouds, possesses bulk characteristics consistent with those observed in nature, though the cold gas is not as radially extended as that observed in nearby dwarfs, resulting in somewhat excessive surface densities. The lines-of-sight velocity dispersion radial profiles have values that are in good agreement with observed dwarf galaxies, but due to the fact that only the streaming velocities of particles are tracked, a correction to include the thermal velocities can lead to profiles that are quite flat. The ISM power spectra of the simulations appear to possess more power on smaller spatial scales than that of the SMC. We conclude that unavoidable limitations remain due to the unresolved physics of star formation and feedback within pc-scale molecular clouds., Comment: Accepted for publication in MNRAS; 9 pages; 5 figures; Fig.4 de-graded to meet arXiv restrictions

We present an analysis of the neutral hydrogen (HI) properties of a fully cosmological hydrodynamical dwarf galaxy, run with varying simulation parameters. As reported by Governato et al. (2010), the high resolution, high star formation density threshold version of this galaxy is the first simulation to result in the successful reproduction of a (dwarf) spiral galaxy without any associated stellar bulge. We have set out to compare in detail the HI distribution and kinematics of this simulated bulgeless disk with what is observed in a sample of nearby dwarfs. To do so, we extracted the radial gas density profiles, velocity dispersion (e.g., velocity ellipsoid, turbulence), and the power spectrum of structure within the cold interstellar medium from the simulations. The highest resolution dwarf, when using a high density star formation threshold comparable to densities of giant molecular clouds, possesses bulk characteristics consistent with those observed in nature, though the cold gas is not as radially extended as that observed in nearby dwarfs, resulting in somewhat excessive surface densities. The lines-of-sight velocity dispersion radial profiles have values that are in good agreement with observed dwarf galaxies, but due to the fact that only the streaming velocities of particles are tracked, a correction to include the thermal velocities can lead to profiles that are quite flat. The ISM power spectra of the simulations appear to possess more power on smaller spatial scales than that of the SMC. We conclude that unavoidable limitations remain due to the unresolved physics of star formation and feedback within pc-scale molecular clouds., Comment: Accepted for publication in MNRAS; 9 pages; 5 figures; Fig.4 de-graded to meet arXiv restrictions

We present an analysis of the neutral hydrogen (H I) properties of a fully cosmological hydrodynamical dwarf galaxy, run with varying simulation parameters. As reported by Governato et al., the high-resolution, high star formation density threshold version of this galaxy is the first simulation to result in the successful reproduction of a (dwarf) spiral galaxy without any associated stellar bulge. We have set out to compare in detail the H I distribution and kinematics of this simulated bulgeless disc with what is observed in a sample of nearby dwarfs. To do so, we extracted the radial gas density profiles, velocity dispersion (e.g. velocity ellipsoid and turbulence) and the power spectrum of structure within the cold interstellar medium (ISM) from the simulations. The highest resolution dwarf, when using a high-density star formation threshold comparable to densities of giant molecular clouds, possesses bulk characteristics consistent with those observed in nature, though the cold gas is not as radially extended as that observed in nearby dwarfs, resulting in somewhat excessive surface densities. The lines-of-sight velocity dispersion radial profiles have values that are in good agreement with the observed dwarf galaxies, but due to the fact that only the streaming velocities of particles are tracked, a correction to include the thermal velocities can lead to profiles that are quite flat. The ISM power spectra of the simulations appear to possess more power on smaller spatial scales than that of the Small Magellanic Cloud. We conclude that unavoidable limitations remain due to the unresolved physics of star formation and feedback within parsec-scale molecular clouds.

We use the Millennium simulation to probe the correlation between cluster velocities and their shapes and the consequences for measurements of the kinetic Sunyaev-Zeldovich (kSZ) effect. Halos are generally prolate ellipsoids with orientations that are correlated with those of nearby halos. We measure the mean streaming velocities of halos along the lines that separate them, demonstrating that the peculiar velocities and the long axes of halos tend to be somewhat aligned, especially for the most massive halos. Since the kSZ effect is proportional to the line-of-sight velocity and the optical depth of the cluster, the alignment results in a strong enhancement of the kSZ signature in clusters moving along the line of sight. This effect has not been taken into account in many analyses of kSZ signatures., Comment: 5 pages, 5 figures, 1 table; accepted for publication in MNRAS

We use the Millennium simulation to probe the correlation between cluster velocities and their shapes and the consequences for measurements of the kinetic Sunyaev-Zeldovich (kSZ) effect. Halos are generally prolate ellipsoids with orientations that are correlated with those of nearby halos. We measure the mean streaming velocities of halos along the lines that separate them, demonstrating that the peculiar velocities and the long axes of halos tend to be somewhat aligned, especially for the most massive halos. Since the kSZ effect is proportional to the line-of-sight velocity and the optical depth of the cluster, the alignment results in a strong enhancement of the kSZ signature in clusters moving along the line of sight. This effect has not been taken into account in many analyses of kSZ signatures., Comment: 5 pages, 5 figures, 1 table; accepted for publication in MNRAS

We present an analysis of the neutral hydrogen (H i) properties of a fully cosmological hydrodynamical dwarf galaxy, run with varying simulation parameters. As reported by Governato et al., the high-resolution, high star formation density threshold version of this galaxy is the first simulation to result in the successful reproduction of a (dwarf) spiral galaxy without any associated stellar bulge. We have set out to compare in detail the H i distribution and kinematics of this simulated bulgeless disc with what is observed in a sample of nearby dwarfs. To do so, we extracted the radial gas density profiles, velocity dispersion (e.g. velocity ellipsoid and turbulence) and the power spectrum of structure within the cold interstellar medium (ISM) from the simulations. The highest resolution dwarf, when using a high-density star formation threshold comparable to densities of giant molecular clouds, possesses bulk characteristics consistent with those observed in nature, though the cold gas is not as radially extended as that observed in nearby dwarfs, resulting in somewhat excessive surface densities. The lines-of-sight velocity dispersion radial profiles have values that are in good agreement with the observed dwarf galaxies, but due to the fact that only the streaming velocities of particles are tracked, a correction to include the thermal velocities can lead to profiles that are quite flat. The ISM power spectra of the simulations appear to possess more power on smaller spatial scales than that of the Small Magellanic Cloud. We conclude that unavoidable limitations remain due to the unresolved physics of star formation and feedback within parsec-scale molecular clouds