While many tensions between Local Group (LG) satellite galaxies and LCDM cosmology have been alleviated through recent cosmological simulations, the spatial distribution of satellites remains an important test of physical models and physical versus numerical disruption in simulations. Using the FIRE-2 cosmological zoom-in baryonic simulations, we examine the radial distributions of satellites with Mstar > 10^5 Msun around 8 isolated Milky Way- (MW) mass host galaxies and 4 hosts in LG-like pairs. We demonstrate that these simulations resolve the survival and physical destruction of satellites with Mstar >~ 10^5 Msun. The simulations broadly agree with LG observations, spanning the radial profiles around the MW and M31. This agreement does not depend strongly on satellite mass, even at distances <~ 100 kpc. Host-to-host variation dominates the scatter in satellite counts within 300 kpc of the hosts, while time variation dominates scatter within 50 kpc. More massive host galaxies within our sample have fewer satellites at small distances, likely because of enhanced tidal destruction of satellites via the baryonic disks of host galaxies. Furthermore, we quantify and provide fits to the tidal depletion of subhalos in baryonic relative to dark matter-only simulations as a function of distance. Our simulated profiles imply observational incompleteness in the LG even at Mstar >~ 10^5 Msun: we predict 2-10 such satellites to be discovered around the MW and possibly 6-9 around M31. To provide cosmological context, we compare our results with the radial profiles of satellites around MW analogs in the SAGA survey, finding that our simulations are broadly consistent with most SAGA systems.

Giant molecular clouds (GMCs) are well-studied in the local Universe, however, exactly how their properties vary during galaxy evolution is poorly understood due to challenging resolution requirements, both observational and computational. We present the first time-dependent analysis of giant molecular clouds in a Milky Way-like galaxy and an LMC-like dwarf galaxy of the FIRE-2 (Feedback In Realistic Environments) simulation suite, which have sufficient resolution to predict the bulk properties of GMCs in cosmological galaxy formation self-consistently. We show explicitly that the majority of star formation outside the galactic center occurs within self-gravitating gas structures that have properties consistent with observed bound GMCs. We find that the typical cloud bulk properties such as mass and surface density do not vary more than a factor of 2 in any systematic way after the first Gyr of cosmic evolution within a given galaxy from its progenitor. While the median properties are constant, the tails of the distributions can briefly undergo drastic changes, which can produce very massive and dense self-gravitating gas clouds. Once the galaxy forms, we identify only two systematic trends in bulk properties over cosmic time: a steady increase in metallicity produced by previous stellar populations and a weak decrease in bulk cloud temperatures. With the exception of metallicity we find no significant differences in cloud properties between the Milky Way-like and dwarf galaxies. These results have important implications for cosmological star and star cluster formation and put especially strong constraints on theories relating the stellar initial mass function to cloud properties.

We present models of CO(1-0) emission from Milky Way-mass galaxies at redshift zero in the FIRE-2 cosmological zoom-in simulations. We calculate the molecular abundances by post-processing the simulations with an equilibrium chemistry solver while accounting for the effects of local sources, and determine the emergent CO(1-0) emission using a line radiative transfer code. We find that the results depend strongly on the shielding length assumed, which in our models sets the attenuation of the incident UV radiation field. Commonly used choices for the shielding length, such as the Jeans length, result in CO abundances that are too high at a given H$_2$ abundance. We find that a model with a distribution of shielding lengths, which has a median shielding length of $\sim 3$ pc in cold gas ($T < 300$ K) for both CO and H$_{2}$, is able to reproduce both the observed CO(1-0) luminosity and inferred CO-to-H$_{2}$ conversion factor at a given star formation rate compared with observations. We suggest that this short shielding length can be thought of as a subgrid model which controls the amount of radiation that penetrates giant molecular clouds.

Surveys of the Milky Way (MW) and M31 enable detailed studies of stellar populations across ages and metallicities, with the goal of reconstructing formation histories across cosmic time. These surveys motivate key questions for galactic archaeology in a cosmological context: when did the main progenitor of a MW/M31-mass galaxy form, and what were the galactic building blocks that formed it? We investigate the formation times and progenitor galaxies of MW/M31-mass galaxies using the FIRE-2 cosmological simulations, including 6 isolated MW/M31-mass galaxies and 6 galaxies in Local Group (LG)-like pairs at z = 0. We examine main progenitor "formation" based on two metrics: (1) transition from primarily ex-situ to in-situ stellar mass growth and (2) mass dominance compared to other progenitors. We find that the main progenitor of a MW/M31-mass galaxy emerged typically at z ~ 3-4 (11.6-12.2 Gyr ago), while stars in the bulge region (inner 2 kpc) at z = 0 formed primarily in a single main progenitor at z < 5 (< 12.6 Gyr ago). Compared with isolated hosts, the main progenitors of LG-like paired hosts emerged significantly earlier (\Delta z ~ 2, \Delta t ~ 1.6 Gyr), with ~ 4x higher stellar mass at all z > 4 (> 12.2 Gyr ago). This highlights the importance of environment in MW/M31-mass galaxy formation, especially at early times. Overall, about 100 galaxies with M_star > 10^5 M_sun formed a typical MW/M31-mass system. Thus, surviving satellites represent a highly incomplete census (by ~ 5x) of the progenitor population.

Far-Ultraviolet (FUV) radiation greatly exceeds ultraviolet, supernovae and winds in the energy budget of young star clusters but is poorly modelled in galaxy simulations. We present results of the first full isolated galaxy disk simulations to include FUV radiation self-consistently. This is the first science application of the TREVR radiative transfer algorithm. We find that FUV radiation alone cannot regulate star formation. However, FUV radiation is crucial for producing warm neutral gas. FUV is also a long-range feedback and is more important in the outer disks of galaxies. We also use the super-bubble feedback model, which depends only on the supernova energy per stellar mass, is more physically realistic than common, parameter-driven alternatives and thus better constrains supernova feedback impacts. FUV and supernovae together can regulate star formation, reproducing the tight Kennicutt-Schmidt relation observed for an individual galaxy, without producing too much hot ionized medium and with less disruption to the ISM compared to supernovae alone.

We present the first measurement of the lifetimes of Giant Molecular Clouds (GMCs) in cosmological simulations at $z = 0$, using the Latte suite of FIRE-2 simulations of Milky Way-mass galaxies. We track GMCs with total gas mass $\gtrsim 10^5$ M$_\odot$ at high spatial ($\sim1$ pc), mass ($7100$ M$_{\odot}$), and temporal (1 Myr) resolution. Our simulated GMCs are consistent with the distribution of masses for massive GMCs in the Milky Way and nearby galaxies. We find GMC lifetimes of $5-7$ Myr, or 1-2 freefall times, on average, with less than 1\% of clouds living longer than 20 Myr. We find increasing GMC lifetimes with galactocentric radius, implying that environment affects the evolutionary cycle of GMCs. However, our GMC lifetimes show no systematic dependence on GMC mass or amount of star formation. These results are broadly consistent with inferences from the literature and provide an initial investigation into ultimately understanding the physical processes that govern GMC lifetimes in a cosmological setting.

(Abridged from Executive Summary) This white paper focuses on the interdisciplinary fields of astrostatistics and astroinformatics, in which modern statistical and computational methods are applied to and developed for astronomical data. Astrostatistics and astroinformatics have grown dramatically in the past ten years, with international organizations, societies, conferences, workshops, and summer schools becoming the norm. Canada's formal role in astrostatistics and astroinformatics has been relatively limited, but there is a great opportunity and necessity for growth in this area. We conducted a survey of astronomers in Canada to gain information on the training mechanisms through which we learn statistical methods and to identify areas for improvement. In general, the results of our survey indicate that while astronomers see statistical methods as critically important for their research, they lack focused training in this area and wish they had received more formal training during all stages of education and professional development. These findings inform our recommendations for the LRP2020 on how to increase interdisciplinary connections between astronomy and statistics at the institutional, national, and international levels over the next ten years. We recommend specific, actionable ways to increase these connections, and discuss how interdisciplinary work can benefit not only research but also astronomy's role in training Highly Qualified Personnel (HQP) in Canada.

We present a new method to study the characteristic scales of collapse and fragmentation in galactic disks. Clump formation is seeded in simulations via controlled perturbations with a specified wavelength and velocity. These are applied to otherwise quiet gas disks ranging from analogues of present day spirals to gas-rich, high-redshift galaxies. The results are compared to linear theory, turbulently perturbed disks and observations. The results reflect the expectations of linear, non-axisymmetric theory with a finite window for growth into a bound clump. We identify two new modes of clump formation: rotation-driven fission and fragmentation of tidal tails, though both are expected to rarely contribute to clump formation in observed disks. We find that bound clumps are generally much smaller than the so-called Toomre mass. The preferred scale for fragmentation increases with the disk gas mass but cannot produce bound objects larger than $\sim10^9$ M$_{\odot}$. The most likely bound clump mass increases from $3\times10^6$ in low mass disks up to $5\times10^8$ M$_{\odot}$. We conclude that observed massive stellar and gaseous clumps on 1 kpc scales at high redshift are most likely aggregates of many initially distinct bound clumps.

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

Tidal encounters are believed to be one of the key drivers of galactic spiral structure in the Universe. Such spirals are expected to produce different morphological and kinematic features compared to density wave and dynamic spiral arms. In this work we present high resolution simulations of a tidal encounter of a small mass companion with a disc galaxy. Included are the effects of gas cooling and heating, star formation and stellar feedback. The structure of the perturbed disc differs greatly from the isolated galaxy, showing clear spiral features that act as sites of new star formation, and displaying interarm spurs. The two arms of the galaxy, the bridge and tail, appear to behave differently; with different star formation histories and structure. Specific attention is focused on offsets between gas and stellar spiral features which can be directly compared to observations. We find some offsets do exist between different media, with gaseous arms appearing mostly on the convex side of the stellar arms, though the exact locations appear highly time dependent. These results further highlight the differences between tidal spirals and other theories of arm structure.

We explore the regulation of star formation in star-forming galaxies through a suite of high-resolution isolated galaxy simulations. We use the SPH code GASOLINE, including photoelectric heating and metal cooling, which produces a multi-phase interstellar medium. We show that representative star formation and feedback sub-grid models naturally lead to a weak, sub-linear dependence between the amount of star formation and changes to star formation parameters. We incorporate these sub-grid models into an equilibrium pressure-driven regulation framework. We show that the sub-linear scaling arises as a consequence of the non-linear relationship between scale height and the effective pressure generated by stellar feedback. Thus, simulated star-formation regulation is sensitive to how well vertical structure in the ISM is resolved. Full galaxy disks experience density waves which drive locally time-dependent star formation. We develop a simple time-dependent, pressure-driven model that reproduces the response extremely well.

With ALMA making it possible to resolve giant molecular clouds (GMCs) in other galaxies, it is becoming necessary to quantify the observational bias on measured GMC properties. Using a hydrodynamical simulation of a barred spiral galaxy, we compared the physical properties of GMCs formed in position-position-position space (PPP) to the observational position-position-velocity space (PPV). We assessed the effect of disc inclination: face-on (PPV_face) and edge-on (PPV_edge), and resolution: 1.5 pc versus 24 pc, on GMC properties and the further implications of using Larson's scaling relations for mass-radius and velocity dispersion-radius. The low-resolution PPV data are generated by simulating ALMA Cycle 3 observations using the CASA package. Results show that the median properties do not differ strongly between PPP and PPV_face under both resolutions, but PPV_edge clouds deviate from these two. The differences become magnified when switching to the lower, but more realistic resolution. The discrepancy can lead to opposite results for the virial parameter's measure of gravitational binding, and therefore the dynamical state of the clouds. The power-law indices for the two Larson's scaling relations decrease going from PPP, PPV_face to PPV_edge and decrease from high to low resolutions. We conclude that the relations are not entirely driven by the underlying physical origin and therefore have to be used with caution when considering the environmental dependence, dynamical state, and the extragalactic CO-to-H2 conversion factor of GMCs.

We explore the effect of galactic environment on properties of molecular clouds. Using clouds formed in a large-scale galactic disc simulation, we measure the observable properties from synthetic column density maps. We confirm that a significant fraction of unbound clouds forms naturally in a galactic disc environment and that a mixed population of bound and unbound clouds can match observed scaling relations and distributions for extragalactic molecular clouds. By dividing the clouds into inner and outer disc populations, we compare their distributions of properties and test whether there are statistically significant differences between them. We find that clouds in the outer disc have lower masses, sizes, and velocity dispersions as compared to those in the inner disc for reasonable choices of the inner/outer boundary. We attribute the differences to the strong impact of galactic shear on the disc stability at large galactocentric radii. In particular, our Toomre analysis of the disc shows a narrowing envelope of unstable masses as a function of radius, resulting in the formation of smaller, lower mass fragments in the outer disc. We also show that the star formation rate is affected by the environment of the parent cloud, and is particularly influenced by the underlying surface density profile of the gas throughout the disc. Our work highlights the strengths of using galaxy-scale simulations to understand the formation and evolution of cloud properties - and the star formation within them - in the context of their environment.

As both simulations and observations reach the resolution of the star-forming molecular clouds, it becomes important to clarify if these two techniques are discussing the same objects in galaxies. We compare clouds formed in a high resolution galaxy simulation identified as continuous structures within a contour, in the simulator's position-position-position (PPP) co-ordinate space and the observer's position-position-velocity space (PPV). Results indicate that the properties of the cloud populations are similar in both methods and up to 70% of clouds have a single counterpart in the opposite data structure. Comparing individual clouds in a one-to-one match reveals a scatter in properties mostly within a factor of two. However, the small variations in mass, radius and velocity dispersion produce significant differences in derived quantities such as the virial parameter. This makes it difficult to determine if a structure is truely gravitationally bound. The three cloud types originally found in the simulation in Fujimoto et al. (2014) are identified in both data sets, with around 80% of the clouds retaining their type between identification methods. We also compared our results when using a peak decomposition method to identify clouds in both PPP and PPV space. The number of clouds increased with this technique, but the overall cloud properties remained similar. However, the more crowded environment lowered the ability to match clouds between techniques to 40%. The three cloud types also became harder to separate, especially in the PPV data set. The method used for cloud identification therefore plays a critical role in determining cloud properties, but both PPP and PPV can potentially identify the same structures.

We present a new stellar feedback model that reproduces superbubbles. Superbubbles from clustered young stars evolve quite differently to individual supernovae and are substantially more efficient at generating gas motions. The essential new components of the model are thermal conduction, sub-grid evaporation and a sub-grid multi-phase treatment for cases where the simulation mass resolution is insufficient to model the early stages of the superbubble. The multi-phase stage is short compared to superbubble lifetimes. Thermal conduction physically regulates the hot gas mass without requiring a free parameter. Accurately following the hot component naturally avoids overcooling. Prior approaches tend to heat too much mass, leaving the hot ISM below $10^6$ K and susceptible to rapid cooling unless ad-hoc fixes were used. The hot phase also allows feedback energy to correctly accumulate from multiple, clustered sources, including stellar winds and supernovae. We employ high-resolution simulations of a single star cluster to show the model is insensitive to numerical resolution, unresolved ISM structure and suppression of conduction by magnetic fields. We also simulate a Milky Way analog and a dwarf galaxy. Both galaxies show regulated star formation and produce strong outflows.

We present the results of a study of simulated Giant Molecular Clouds (GMCs) formed in a Milky Way-type galactic disk with a flat rotation curve. This simulation, which does not include star formation or feedback, produces clouds with masses ranging between 10^4 Msun and 10^7 Msun. We compare our simulated cloud population to two observational surveys; The Boston University- Five College Radio Astronomy Observatory Galactic Ring Survey and the BIMA All-Disk Survey of M33. An analysis of the global cloud properties as well as a comparison of Larson's scaling relations is carried out. We find that simulated cloud properties agree well with the observed cloud properties, with the closest agreement occurring between the clouds at comparable resolution in M33. Our clouds are highly filamentary - a property that derives both from their formation due to gravitational instability in the sheared galactic environment, as well as to cloud- cloud gravitational encounters. We also find that the rate at which potentially star forming gas accumulates within dense regions - wherein n_{thresh} > 10^4 cm^{-3} - is 3% per 10 Myr, in clouds of roughly 10^6 Msun. This suggests that star formation rates in observed clouds are related to the rates at which gas can be accumulated into dense subregions within GMCs via filamentary flows. The most internally well-resolved clouds are chosen for listing in a catalogue of simulated GMCs; the first of its kind. The catalogued clouds are available as an extracted data set from the global simulation.