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.

Globular clusters (GCs) are typically old, with most having formed at z >~ 2. This makes understanding their birth environments difficult, as they are typically too distant to observe with sufficient angular resolution to resolve GC birth sites. Using 25 cosmological zoom-in simulations of Milky Way-like galaxies from the E-MOSAICS project, with physically-motivated models for star formation, feedback, and the formation, evolution, and disruption of GCs, we identify the birth environments of present-day GCs. We find roughly half of GCs in these galaxies formed in-situ (52.0 +/- 1.0 per cent) between z ~ 2 - 4, in turbulent, high-pressure discs fed by gas that was accreted without ever being strongly heated through a virial shock or feedback. A minority of GCs form during mergers (12.6 +/- 0.6 per cent in major mergers, and 7.2 +/- 0.5 per cent in minor mergers), but we find that mergers are important for preserving the GCs seen today by ejecting them from their natal, high density interstellar medium (ISM), where proto-GCs are rapidly destroyed due to tidal shocks from ISM substructure. This chaotic history of hierarchical galaxy assembly acts to mix the spatial and kinematic distribution of GCs formed through different channels, making it difficult to use observable GC properties to distinguish GCs formed in mergers from ones formed by smooth accretion, and similarly GCs formed in-situ from those formed ex-situ. These results suggest a simple picture of GC formation, in which GCs are a natural outcome of normal star formation in the typical, gas-rich galaxies that are the progenitors of present-day galaxies.

Quenched galaxies are often observed to contain a strong bulge component. The key question is whether this reflects a causal connection - can star formation be quenched dynamically by bulges or the spheroids of early-type galaxies? We systematically investigate the impact of these morphological components on star formation, by performing a suite of hydrodynamical simulations of isolated galaxies containing a spheroid. We vary the bulge mass and scale radius, while the total initial stellar, halo and gas mass are kept constant, with a gas fraction of 5 per cent. In addition, we consider two different sub-grid star formation prescriptions. The first follows most simulations in the literature by assuming a constant star formation efficiency per free-fall time, whereas in the second model it depends on the gas virial parameter, following high-resolution simulations of turbulent fragmentation. Across all simulations, central spheroids increase the gas velocity dispersion towards the galactic centre. This increases the gravitational stability of the gas disc, suppresses fragmentation and star formation, and results in galaxies hosting extremely smooth and quiescent gas discs that fall below the galaxy main sequence. These effects amplify when using the more sophisticated, dynamics-dependent star formation model. Finally, we discover a pronounced relation between the central stellar surface density and star formation rate (SFR), such that the most bulge-dominated galaxies show the strongest deviation from the main sequence. We conclude that the SFR of galaxies is not only set by the balance between accretion and feedback, but carries a (sometimes dominant) dependence on the gravitational potential.

We present a detailed comparison of the Milky Way (MW) globular cluster (GC) kinematics with the 25 Milky Way-mass cosmological simulations from the E-MOSAICS project. While the MW falls within the kinematic distribution of GCs spanned by the simulations, the relative kinematics of its metal-rich ($[\rm{Fe/H}]>-1.2$) versus metal-poor ($[\rm{Fe/H}]<-1.2$), and inner ($r<8\rm{kpc}$) versus outer ($r>8\rm{kpc}$) populations are atypical for its mass. To understand the origins of these features, we perform a comprehensive statistical analysis of the simulations, and find 18 correlations describing the assembly of $L^*$ galaxies and their dark matter haloes based on their GC population kinematics. The correlations arise because the orbital distributions of accreted and in-situ GCs depend on the masses and accretion redshifts of accreted satellites, driven by the combined effects of dynamical fraction, tidal stripping, and dynamical heating. Because the kinematics of in-situ/accreted GCs are broadly traced by the metal-rich/metal-poor and inner/outer populations, the observed GC kinematics are a sensitive probe of galaxy assembly. We predict that relative to the population of $L^*$ galaxies, the MW assembled its dark matter and stellar mass rapidly through a combination of in-situ star formation, more than a dozen low-mass mergers, and $1.4\pm1.2$ early ($z=3.1\pm1.3$) major merger. The rapid assembly period ended early, limiting the fraction of accreted stars. We conclude by providing detailed quantitative predictions for the assembly history of the MW.

Feedback from core collapse supernovae (SNe), the final stage of evolution of the most massive stars, has long been a key element in simulations of galaxy formation. In this paper, we examine how simplifying assumptions made in approximating the SN rates along with underlying physical uncertainty in those rates can lead to large variations in the overall evolution of simulated Milky Way-like galaxies. We find that the clustering of star formation is strongly impacted by the delay between star formation and SN feedback. In addition, the choice to use a realistic delay time distribution or instantaneous injection for SN can have a significant impact on the galaxy. These effects appear even when identical sub-grid models are used for coupling SN energy and momentum, and the total SN energy budget is kept constant. In addition, we show that the uncertain minimum SN progenitor mass has a significant impact on the SN energy budget and injection timescale, and can completely change the overall evolution of the galaxy. These underlying uncertainties mean that despite advances in the sub-grid modelling of SN feedback, there are still serious difficulties in constraining the effects of SN feedback. This complicates the task of comparing different simulations to each other, as well as comparing simulations to observations. We conclude by providing practical limits on the parameters of subgrid models for SN feedback, which bound the uncertainty arising from SN progenitor physics for future predictions from galaxy simulations.

We present a new theoretical framework for using entropy to understand how outflows driven by supernovae are launched from disc galaxies: via continuous, buoyant acceleration through the circumgalactic medium (CGM). When young star clusters detonate supernovae in the interstellar medium (ISM) of a galaxy, they generate hot, diffuse bubbles that push on the surrounding ISM and evaporate that ISM into their interiors. As these bubbles reach the scale height of the ISM, they break out of the disc, rising into the CGM. Once these bubbles break out, if they have sufficiently high entropy, they will feel an upward acceleration, owing to a local buoyant force. This upward force will accelerate these bubbles, driving them to high galactocentric radii, keeping them in the CGM for $>\Gyr$, even if their initial velocity is much lower than the local escape velocity. We derive an equation of motion for these entropy-driven winds that connects the ISM properties, halo mass, and CGM profile of galaxies to the ultimate evolution of feedback-driven winds. We explore the parameter space of these equations, and show how this novel framework can explain both self-consistent simulations of star formation and galactic outflows as well as the new wealth of observations of CGM kinematics. We show that these entropy-driven winds can produce long wind recycling times, while still carrying a significant amount of mass. Comparisons to simulations and observations show entropy-driven winds convincingly explain the kinematics of galactic outflows.

We report on the development of a polarization-sensitive dichroic (150/220 GHz) detector array for the Cosmology Large Angular Scale Surveyor (CLASS) delivered to the telescope site in June 2019. In concert with existing 40 and 90 GHz telescopes, the 150/220 GHz telescope will make observations of the cosmic microwave background over large angular scales aimed at measuring the primordial B-mode signal, the optical depth to reionization, and other fundamental physics and cosmology. The 150/220 GHz focal plane array consists of three detector modules with 1020 transition edge sensor (TES) bolometers in total. Each dual-polarization pixel on the focal plane contains four bolometers to measure the two linear polarization states at 150 and 220 GHz. Light is coupled through a planar orthomode transducer (OMT) fed by a smooth-walled feedhorn array made from an aluminum-silicon alloy (CE7). In this work, we discuss the design, assembly, and in-lab characterization of the 150/220 GHz detector array. The detectors are photon-noise limited, and we estimate the total array noise-equivalent power (NEP) to be 2.5 and 4 aW$\sqrt{\mathrm{s}}$ for 150 and 220 GHz arrays, respectively.

As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release.

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.

The stellar birth environment can significantly shorten protoplanetary disc (PPD) lifetimes due to the influence of stellar feedback mechanisms. The degree to which these mechanisms suppress the time and mass available for planet formation is dependent on the local far-ultraviolet (FUV) field strength, stellar density, and ISM properties. In this work, we present the first theoretical framework quantifying the distribution of PPD dispersal time-scales as a function of parameters that describe the galactic environment. We calculate the probability density function for FUV flux and stellar density in the solar neighbourhood. In agreement with previous studies, we find that external photoevaporation is the dominant environment-related factor influencing local stellar populations after the embedded phase. Applying our general prescription to the Central Molecular Zone of the Milky Way (i.e. the central ~250 pc), we predict that 90% of PPDs in the region are destroyed within 1 Myr of the dispersal of the parent molecular cloud. Even in such dense environments, we find that external photoevaporation is the dominant disc depletion mechanism over dynamical encounters between stars. PPDs around low-mass stars are particularly sensitive to FUV-induced mass loss, due to a shallower gravitational potential. For stars of mass ~1 $M_\odot$, the solar neighbourhood lies at approximately the highest gas surface density for which PPD dispersal is still relatively unaffected by external FUV photons, with a median PPD dispersal timescale of ~4 Myr. We highlight the key questions to be addressed to further contextualise the significance of the local galactic environment for planet formation.

We use the NIHAO simulations to investigate the effects of baryonic physics on the time evolution of Dark Matter central density profiles. The sample is made of $\approx 70$ independent high resolution hydrodynamical simulations of galaxy formation and covers a wide mass range: 1e10< Mhalo <1e12, i.e., from dwarfs to L* . We confirm previous results on the dependence of the inner dark matter density slope, $\alpha$, on the ratio between stellar-to-halo mass. We show that this relation holds approximately at all redshifts (with an intrinsic scatter of ~0.18 in $\alpha$). This implies that in practically all haloes the shape of their inner density profile changes quite substantially over cosmic time, as they grow in stellar and total mass. Thus, depending on their final stellar-to-halo mass ratio, haloes can either form and keep a substantial density core (size~1 kpc), or form and then destroy the core and re-contract the halo, going back to a cuspy profile, which is even steeper than CDM predictions for massive galaxies (~1e12 Msun). We show that results from the NIHAO suite are in good agreement with recent observational measurements of $\alpha$ in dwarf galaxies. Overall our results suggest that the notion of a universal density profile for dark matter haloes is no longer valid in the presence of galaxy formation.

We present a study of environmental effects and preprocessing in a large galaxy group using a high-resolution, zoom-in simulation run with the GASOLINE2 hydrodynamics code. We categorize galaxies that were always in distinct haloes as unaccreted, galaxies that were distinct before accretion onto the main group as single, and galaxies that were in external sub-groups before accretion onto the main group as grouped. The unaccreted galaxy population experiences steady growth in dark matter, gas and stellar mass. Both single- and group-accreted galaxies begin to lose dark matter and gas after first accretion onto any host but continue to grow in stellar mass. Individual trajectories show that galaxies cease mass growth within roughly three virial radii of the main group. Single galaxies continue to form stars until the group virial radius is crossed, when they begin to lose both dark matter and gas. Grouped galaxies peak in mass when joining their external sub-group, indicating that they experience preprocessing. Most accreted galaxies retain their accumulated stellar mass. The total mass loss is dominated by tidal stripping, with evidence for additional gas stripping via ram pressure. Most accreted galaxies are quenched $\sim$(0.5-2.5) Gyr after accretion onto any group. These differing histories place unaccreted, single and grouped galaxies in distinct regions of the stellar mass-to-halo mass (SMHM) relation. This suggests that preprocessed galaxies are a key source of scatter in the SMHM relation for mixed galaxy populations.

The evolution of galaxies is governed by equations with chaotic solutions: gravity and compressible hydrodynamics. While this micro-scale chaos and stochasticity has been well studied, it is poorly understood how it couples to macro-scale properties examined in simulations of galaxy formation. In this paper, we show how perturbations introduced by floating-point roundoff, random number generators, and seemingly trivial differences in algorithmic behaviour can produce non-trivial differences in star formation histories, circumgalactic medium (CGM) properties, and the distribution of stellar mass. We examine the importance of stochasticity due to discreteness noise, variations in merger timings and how self-regulation moderates the effects of this stochasticity. We show that chaotic variations in stellar mass can grow until halted by feedback-driven self-regulation or gas exhaustion. We also find that galaxy mergers are critical points from which large (as much as a factor of 2) variations in quantities such as the galaxy stellar mass can grow. These variations can grow and persist for more than a Gyr before regressing towards the mean. These results show that detailed comparisons of simulations require serious consideration of the magnitude of effects compared to run-to-run chaotic variation, and may significantly complicate interpreting the impact of different physical models. Understanding the results of simulations requires us to understand that the process of simulation is not a mapping of an infinitesimal point in configuration space to another, final infinitesimal point. Instead, simulations map a point in a space of possible initial conditions points to a volume of possible final states.

In a series of three papers, Eadie et al. developed a hierarchical Bayesian method to estimate the Milky Way Galaxy's mass given a physical model for the potential, a measurement model, and kinematic data of test particles such as globular clusters (GCs) or halo stars in the Galaxy's halo. The Galaxy's virial mass was found to have a 95\% Bayesian credible region (c.r.) of $(0.67, 1.09) \times 10^{12} M_{\odot}$. In the present study, we test the hierarchical Bayesian method against simulated galaxies created in the McMaster Unbiased Galaxy Simulations 2 (MUGS2), for which the true mass is known. We estimate the masses of MUGS2 galaxies using GC analogs from the simulations as tracers. The analysis, completed as a blind test, recovers the true $M_{200}$ of the MUGS2 galaxies within 95\% Bayesian c.r. in 8 out of 18 cases. Of the 10 galaxy masses that were not recovered within the 95\% c.r., a large subset have posterior distributions that occupy extreme ends of the parameter space allowed by the priors. A few incorrect mass estimates are explained by the exceptional evolution history of the galaxies. We also find evidence that the model cannot describe both the galaxies' inner and outer structure simultaneously in some cases. After removing the GC analogs associated with the galactic disks, the true masses were found more reliably (13 out of 18 were predicted within the c.r.). Finally, we discuss how representative the GC analogs are of the real GC population in the Milky Way.

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.

The methods in the Gasoline2 Smoothed Particle Hydrodynamics (SPH) code are described and tested. Gasoline2 is the most recent version of the Gasoline code for parallel hydrodynamics and gravity with identical hydrodynamics to the Changa code. As with other Modern SPH codes, we prevent sharp jumps in time steps, use upgraded kernels and larger neighbour numbers and employ local viscosity limiters. Unique features in Gasoline2 include its Geometric-Density-Average Force expression, explicit Turbulent Diffusion terms and Gradient-Based shock detection to limit artificial viscosity. This last feature allows Gasoline2 to completely avoid artificial viscosity in non-shocking compressive flows. We present a suite of tests demonstrating the value of these features with the same code configuration and parameter choices used for production simulations.

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.

Recent analysis \citep{McGaugh2016} of the SPARC galaxy sample found a surprisingly tight relation between the radial acceleration inferred from the rotation curves, and the acceleration due to the baryonic components of the disc. It has been suggested that this relation may be evidence for new physics, beyond $\Lambda CDM$. In this letter we show that 32 galaxies from the MUGS2 match the SPARC acceleration relation. These cosmological simulations of star forming, rotationally supported discs were simulated with a {\sc WMAP3} $\Lambda CDM$ cosmology, and match the SPARC acceleration relation with less scatter than the observational data. These results show that this acceleration relation is a consequence of dissipative collapse of baryons, rather than being evidence for exotic dark-sector physics or new dynamical laws.

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.

We use the NIHAO (Numerical Investigation of Hundred Astrophysical Objects) cosmological simulations to study the effects of galaxy formation on key properties of dark matter (DM) haloes. NIHAO consists of $\simeq 90$ high-resolution SPH simulations that include (metal-line) cooling, star formation, and feedback from massive stars and SuperNovae, and cover a wide stellar and halo mass range: $10^6 < M_* / M_{\odot} < 10^{11}$ ( $10^{9.5} < M_{\rm halo} / M_{\odot} < 10^{12.5}$). When compared to DM-only simulations, the NIHAO haloes have similar shapes at the virial radius, R_{\rm vir}, but are substantially rounder inside $\simeq 0.1R_{\rm vir}$. In NIHAO simulations $c/a$ increases with halo mass and integrated star formation efficiency, reaching $\sim 0.8$ at the Milky Way mass (compared to 0.5 in DM-only), providing a plausible solution to the long-standing conflict between observations and DM-only simulations. The radial profile of the phase-space $Q$ parameter ($\rho/\sigma^3$) is best fit with a single power law in DM-only simulations, but shows a flattening within $\simeq 0.1R_{\rm vir}$ for NIHAO for total masses $M>10^{11} M_{\odot}$. Finally, the global velocity distribution of DM is similar in both DM-only and NIHAO simulations, but in the solar neighborhood, NIHAO galaxies deviate substantially from Maxwellian. The distribution is more symmetric, roughly Gaussian, with a peak that shifts to higher velocities for Milky Way mass haloes. We provide the distribution parameters which can be used for predictions for direct DM detection experiments. Our results underline the ability of the galaxy formation processes to modify the properties of dark matter haloes.

We study the dependence of fragmentation in massive gas-rich galaxy disks at $z > 1$ on feedback model and hydrodynamical method, employing the GASOLINE2 SPH code and the lagrangian mesh-less code GIZMO in finite mass mode. We compare non-cosmological galaxy disk runs with standard blastwave supernovae (SN)feedback, which introduces delayed cooling in order to drive winds, and runs with the new superbubble SN feedback, which produces winds naturally by modelling the detailed physics of SN-driven bubbles and leads to efficient self-regulation of star formation. We find that, with blastwave feedback, massive star forming clumps form in comparable number and with very similar masses in GASOLINE2 and GIZMO. The typical masses are in the range $10^7-10^8 M_{\odot}$, lower than in most previous works, while giant clumps with masses above $10^9 M_{\odot}$ are exceedingly rare. With superbubble feedback, instead, massive bound star forming clumps do not form because galaxies never undergo a phase of violent disk instability. Only sporadic, unbound star forming overdensities lasting only a few tens of Myr can arise that are triggered by perturbations of massive satellite companions. We conclude that there is a severe tension between explaining massive star forming clumps observed at $z > 1$ primarily as the result of disk fragmentation driven by gravitational instability and the prevailing view of feedback-regulated galaxy formation. The link between disk stability and star formation efficiency should thus be regarded as a key testing ground for galaxy formation theory.

We explore when supernovae can (and cannot) regulate the star formation and bulge growth in galaxies based on a sample of 18 simulated galaxies. The simulations include key physics such as evaporation and conduction, neglected in prior work, and required to correctly model superbubbles resulting from stellar feedback. We show that for galaxies with virial masses $>10^{12}\;M_\odot$, supernovae alone cannot prevent excessive star formation. This failure occurs due to a shutdown of galactic winds, with wind mass loadings falling from $\eta\sim10$ to $\eta<1$. In more massive systems, this transfer of baryons to the circumgalactic medium falters earlier on and the galaxies diverge significantly from observed galaxy scaling relations and morphologies. The decreasing efficiency is simply due to a deepening potential well preventing gas escape. This implies that non-supernova feedback mechanisms must become dominant for galaxies with stellar masses greater than $\sim4\times10^{10}\;M_\odot$. The runaway growth of the central stellar bulge, strongly linked to black hole growth, suggests that feedback from active galactic nuclei is the probable mechanism. Below this mass, supernovae alone are able to produce a realistic stellar mass fraction, star formation history and disc morphology.

We present the first cosmological galaxy evolved using the modern smoothed particle hydrodynamics (SPH) code GASOLINE2 with superbubble feedback. We show that superbubble-driven galactic outflows powered by Type II supernovae alone can produce $\rm{L^*}$ galaxies with flat rotation curves with circular velocities $\sim 200\; \rm{km/s}$, low bulge-to-disc ratios, and stellar mass fractions that match observed values from high redshift to the present. These features are made possible by the high mass loadings generated by the evaporative growth of superbubbles. Outflows are driven extremely effectively at high redshift, expelling gas at early times and preventing overproduction of stars before $z=2$. Centrally concentrated gas in previous simulations has often lead to unrealistically high bulge to total ratios and strongly peaked rotation curves. We show that supernova-powered superbubbles alone can produce galaxies that agree well with observed properties without the need for additional feedback mechanisms or increased feedback energy. We present additional results arising from properly modelled hot feedback.

We introduce project NIHAO (Numerical Investigation of a Hundred Astrophysical Objects), a set of 100 cosmological zoom-in hydrodynamical simulations performed using the gasoline code, with an improved implementation of the SPH algorithm. The haloes in our study range from dwarf to Milky Way masses, and represent an unbiased sampling of merger histories, concentrations and spin parameters. The particle masses and force softenings are chosen to resolve the mass profile to below 1% of the virial radius at all masses, ensuring that galaxy half-light radii are well resolved. Using the same treatment of star formation and stellar feedback for every object, the simulated galaxies reproduce the observed inefficiency of galaxy formation across cosmic time as expressed through the stellar mass vs halo mass relation, and the star formation rate vs stellar mass relation. We thus conclude that stellar feedback is the chief piece of physics required to limit the efficiency of star formation in galaxies less massive than the Milky Way.

We show that the cool gas masses of galactic discs reach a steady state that lasts many Gyr after their last major merger in cosmological hydrodynamic simulations. The mass of disc gas, M$_{\rm gas}$, depends upon a galaxy halo's spin and virial mass, but not upon stellar feedback. Halos with low spin have high star formation efficiency and lower disc gas mass. Similarly, lower stellar feedback leads to more star formation so the gas mass ends up nearly the same irregardless of stellar feedback strength. Even considering spin, the M$_{\rm gas}$ relation with halo mass, M$_{200}$ only shows a factor of 3 scatter. The M$_{\rm gas}$--M$_{200}$ relation show a break at M$_{200}$=$2\times10^{11}$ M$_\odot$ that corresponds to an observed break in the M$_{\rm gas}$--M$_\star$ relation. The constant disc mass stems from a shared halo gas density profile in all the simulated galaxies. In their outer regions, the profiles are isothermal. Where the profile rises above $n=10^{-3}$ cm$^{-3}$, the gas readily cools and the profile steepens. Inside the disc, rotation supports gas with a flatter density profile except where supernova explosions disrupt the disc. Energy injection from stellar feedback also provides pressure support to the halo gas to prevent runaway cooling flows. The resulting constant gas mass makes simpler models for galaxy formation possible, either using a "bathtub" model for star formation rates or when modeling chemical evolution.

Stacking polarized radio emission in SKA surveys provides statistical information on large samples that is not accessible otherwise due to limitations in sensitivity, source statistics in small fields, and averaging over frequency (including Faraday synthesis). Polarization is a special case because one obvious source of stacking targets is the Stokes I source catalog, possibly in combination with external catalogs, for example an SKA HI survey or a non-radio survey. We point out the significance of stacking sub-samples selected by additional observable parameters to investigate relations that reveal more about the physics of the source. Applications of stacking polarization include, but are not limited to, obtaining in a statistical sense polarization information to the detection limit in total intensity, depolarization as a function of cosmic time at consistent source-frame wavelengths, magnetic field properties in objects with a low radio luminosity such as dwarf and low-surface-brightness galaxies, and investigating potential correlations of observable parameters with the average magnetic field direction in a sample. We also point out the potential use of stacking in validating the polarization calibration of a survey. While stacking is flexible in terms of survey definition, we discuss optimal survey parameters for the science experiments presented, as well as computing and archiving requirements.

We use high resolution Hydro$+$N-Body cosmological simulations to compare the assembly and evolution of a small field dwarf (stellar mass ~ 10$^{6-7}$ M$\odot$, total mass 10$^{10}$ M$\odot$ in $\Lambda$ dominated CDM and 2keV WDM cosmologies. We find that star formation (SF) in the WDM model is reduced and delayed by 1-2 Gyr relative to the CDM model, independently of the details of SF and feedback. Independent of the DM model, but proportionally to the SF efficiency, gas outflows lower the central mass density through `dynamical heating', such that all realizations have circular velocities $<$ 20kms at 500$~$pc, in agreement with local kinematic constraints. As a result of dynamical heating, older stars are less centrally concentrated than younger stars, similar to stellar population gradients observed in nearby dwarf galaxies. Introducing an important diagnostic of SF and feedback models, we translate our simulations into artificial color-magnitude diagrams and star formation histories in order to directly compare to available observations. The simulated galaxies formed most of their stars in many $\sim$10 Myr long bursts. The CDM galaxy has a global SFH, HI abundance and Fe/H and alpha-elements distribution well matched to current observations of dwarf galaxies. These results highlight the importance of directly including `baryon physics' in simulations when 1) comparing predictions of galaxy formation models with the kinematics and number density of local dwarf galaxies and 2) differentiating between CDM and non-standard models with different DM or power spectra.

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 stacking polarized intensity as a means to study the polarization of sources that are too faint to be detected individually in surveys of polarized radio sources. Stacking offers not only high sensitivity to the median signal of a class of radio sources, but also avoids a detection threshold in polarized intensity, and therefore an arbitrary exclusion of source with a low percentage of polarization. Correction for polarization bias is done through a Monte Carlo analysis and tested on a simulated survey. We show that the non-linear relation between the real polarized signal and the detected signal requires knowledge of the shape of the distribution of fractional polarization, which we constrain using the ratio of the upper quartile to the lower quartile of the distribution of stacked polarized intensities. Stacking polarized intensity for NVSS sources down to the detection limit in Stokes I, we find a gradual increase in median fractional polarization that is consistent with a trend that was noticed before for bright NVSS sources, but is much more gradual than found by previous deep surveys of radio polarization. Consequently, the polarized radio source counts derived from our stacking experiment predict fewer polarized radio sources for future surveys with the Square Kilometre Array and its pathfinders.

Detection thresholds in polarized intensity and polarization bias correction are investigated for surveys where the polarization information is obtained from RM synthesis. Considering unresolved sources with a single rotation measure, a detection threshold of $8 \sigma_{QU}$ applied to the Faraday spectrum will retrieve the RM with a false detection rate less than $10^{-4}$, but polarized intensity is more strongly biased than Ricean statistics suggest. For a detection threshold of $5 \sigma_{QU}$, the false detection rate increases to ~4%, depending also on $\lambda^2$ coverage and the extent of the Faraday spectrum. Non-Gaussian noise in Stokes Q and U due to imperfect imaging and calibration can be represented by a distribution that is the sum of a Gaussian and an exponential. The non-Gaussian wings of the noise distribution increase the false detection rate in polarized intensity by orders of magnitude. Monte-Carlo simulations assuming non-Gaussian noise in Q and U, give false detection rates at $8 \sigma_{QU}$ similar to Ricean false detection rates at $4.9 \sigma_{QU}$.