The stellar-to-halo mass relation (SHMR) is one of the main sources of information we have on the connection between galaxies and their dark matter haloes. Here we analyse in detail two popular forms of the SHMR, M*/Mh vs Mh and M*/Mh vs M*, and compare them with another physically motivated scaling relation, M*/Mh vs GMh/j*sigmahat*. Although this relation cannot predict the halo mass explicitely, it connects the stellar mass fraction to fundamental galaxy properties such as specific angular momentum (j*) and velocity dispersion (sigmahat*) via disc gravitational instability. Our detailed comparative analysis is based on one of the largest sample of galaxies with both high-quality rotation curves and near-infrared surface photometry, and leads to the following results: (i) M*/Mh vs Mh and M*/Mh vs M* are not just two alternative parametrizations of the same relation, but two significantly different relations; (ii) M*/Mh vs GMh/j*sigmahat* outperforms the two popular relations in terms of tightness, correlation strength and significance; (iii) j* and sigmahat* play an equally important role in our scaling relation, and it is their interplay that constrains M*/Mh so tightly; (iv) the evolution of M*/Mh, j* and sigmahat* is regulated by disc gravitational instability: when M*/Mh varies, j* and sigmahat* also vary as predicted by our scaling relation, thus erasing the memory of such evolution. This implies that the process of disc gravitational instability is intriguingly uniform across disc galaxies of all morphological types: from lenticulars to blue compact dwarfs. In particular, the cosmic variance of Toomre's Q is 0.2 dex, a universal value for both stars and atomic gas.

Spectroscopic surveys of the Milky Way's stars have revealed spatial, chemical and kinematical structures that encode its history. In this work, we study their origins using a cosmological zoom simulation, VINTERGATAN, of a Milky Way-mass disc galaxy. We find that in connection to the last major merger at $z\sim 1.5$, cosmological accretion leads to the rapid formation of an outer, metal-poor, low-[$\alpha$/Fe] gas disc around the inner, metal-rich galaxy containing the old high-[$\alpha$/Fe] stars. This event leads to a bimodality in [$\alpha$/Fe] over a range of [Fe/H]. A detailed analysis of how the galaxy evolves since $z\sim 1$ is presented. We demonstrate the way in which inside-out growth shapes the radial surface density and metallicity profile and how radial migration preferentially relocates stars from the inner to the outer disc. Secular disc heating is found to give rise to increasing velocity dispersions and scaleheights with stellar age, which together with disc flaring explains several trends observed in the Milky Way, including shallower radial [Fe/H]-profiles above the midplane. We show how the galaxy formation scenario imprints non-trivial mappings between structural associations (i.e. thick and thin discs), velocity dispersions, $\alpha$-enhancements, and ages of stars, e.g. the most metal-poor stars in the low-[$\alpha$/Fe] sequence are found to have a scaleheight comparable to old high-[$\alpha$/Fe] stars. Finally, we illustrate how at low spatial resolution, comparable to the thickness of the galaxy, the proposed pathway to distinct sequences in [$\alpha$/Fe]-[Fe/H] cannot be captured.

Using the VINTERGATAN cosmological zoom simulation, we explore the contributions of the in situ and accreted material, and the effect of galaxy interactions and mergers in the assembly of a Milky Way-like galaxy. We find that the initial growth phase of galaxy evolution, dominated by repeated major mergers, provides the necessary physical conditions for the assembly of a thick, kinematically hot disk populated by high-[$\alpha$/Fe] stars, formed both in situ and in accreted satellite galaxies. We find that the diversity of evolutionary tracks followed by the simulated galaxy and its progenitors leads to very little overlap of the in situ and accreted populations for any given chemical composition. At a given age, the spread in [$\alpha$/Fe] abundance ratio results from the diversity of physical conditions in VINTERGATAN and its satellites, with an enhancement in [$\alpha$/Fe] found in stars formed during starburst episodes. Later, the cessation of the merger activity promotes the in situ formation of stars in the low-[$\alpha$/Fe] regime, in a radially extended, thin and overall kinematically colder disk, thus establishing chemically bimodal thin and thick disks, in line with observations. We draw links between notable features in the [Fe/H] - [$\alpha$/Fe] plane with their physical causes, and propose a comprehensive formation scenario explaining self-consistently, in the cosmological context, the main observed properties of the Milky Way.

Using the cosmological zoom simulation VINTERGATAN, we present a new scenario for the onset of star formation at the metal-poor end of the low-[$\alpha$/Fe] sequence in a Milky Way-like galaxy. In this scenario, the galaxy is fueled by two distinct gas flows. One is enriched by outflows from massive galaxies, but not the other. While the former feeds the inner galactic region, the latter fuels an outer gas disk, inclined with respect to the main galactic plane, and with a significantly poorer chemical content. The first passage of the last major merger galaxy triggers tidal compression in the outer disk, which increases the gas density and eventually leads to star formation, at a metallicity 0.75 dex lower than the inner galaxy. This forms the first stars of the low-[$\alpha$/Fe] sequence. These in situ stars have halo-like kinematics, similarly to what is observed in the Milky Way, due to the inclination of the outer disk which eventually aligns with the inner one via gravitational torques. We show that this tilting disk scenario is likely to be common in Milky-Way like galaxies. This process implies that the low-[$\alpha$/Fe] sequence is populated in situ, simultaneously from two formation channels, in the inner and the outer galaxy, with distinct metallicities. This contrasts with purely sequential scenarios for the assembly of the Milky Way disk and could be tested observationally.

Roughly ten per cent of OB stars are kicked out of their natal clusters before ending their life as supernovae. These so-called runaway stars can travel hundreds of parsecs into the low-density interstellar medium, where momentum and energy from stellar feedback is efficiently deposited. In this work we explore how this mechanism affects large scale properties of the galaxy, such as outflows. To do so we use a new model which treats OB stars and their associated feedback processes on a star-by-star basis. With this model we compare two hydrodynamical simulations of Milky Way-like galaxies, one where we include runaways, and one where we ignore them. Including runaway stars leads to twice as many supernovae explosions in regions with gas densities ranging from 1e-5 cm^-3 to 1e-3 cm^-3. This results in more efficient heating of the inter-arm regions, and drives strong galactic winds with mass loading factors boosted by up to one order of magnitude. These outflows produce a more massive and extended multi-phase circumgalactic medium, as well as a population of dense clouds in the halo. Conversely, since less energy and momentum is released in the dense star forming regions, the cold phase of the interstellar medium is less disturbed by feedback effects.

We study how star formation is regulated in low-mass field dwarf galaxies ($10^5 \leq M_{\star} \leq 10^6 \, \text{M}_{\odot}$), using cosmological high-resolution ($3 \, \text{pc}$) hydrodynamical simulations. Cosmic reionization quenches star formation in all our simulated dwarfs, but three galaxies with final dynamical masses of $3 \times 10^{9} \,\text{M}_{\odot}$ are subsequently able to replenish their interstellar medium by slowly accreting gas. Two of these galaxies re-ignite and sustain star formation until the present day at an average rate of $10^{-5} \, \text{M}_{\odot} \, \text{yr}^{-1}$, highly reminiscent of observed low-mass star-forming dwarf irregulars such as Leo T. The resumption of star formation is delayed by several billion years due to residual feedback from stellar winds and Type Ia supernovae; even at $z=0$, the third galaxy remains in a temporary equilibrium with a large gas content but without any ongoing star formation. Using the "genetic modification'' approach, we create an alternative mass growth history for this gas-rich quiescent dwarf and show how a small $(0.2\,\mathrm{dex})$ increase in dynamical mass can overcome residual stellar feedback, re-igniting star formation. The interaction between feedback and mass build-up produces a diversity in the stellar ages and gas content of low-mass dwarfs, which will be probed by combining next-generation HI and imaging surveys.

Galactic outflows driven by stellar feedback are crucial for explaining the inefficiency of galaxy formation. Although strong feedback can promote the formation of galactic discs by limiting star formation at early times and removing low angular momentum gas, it is not understood how the same feedback can result in diverse objects such as elliptical galaxies or razor thin spiral galaxies. In this work, we investigate this problem using cosmological zoom-in simulations of two galaxies forming within $10^{12}~\mathrm{M_\odot}$ halos with almost identical mass accretion histories. At $z>1.5$, the two galaxies feature a surface density of star formation $\Sigma_{\rm SFR}\simeq 10~\mathrm{M_\odot}~{\rm yr}^{-1}~{\rm kpc}^{-2}$, leading to strong outflows. After the last starburst episode, both galaxies feature a dramatic gaseous disc growth from 1 kpc to 5 kpc during 1 Gyr, a decisive epoch we dub "the Grand Twirl". After this episode, the evolutionary tracks diverge, with one galaxy ending up as a bulge-dominated galaxy, whereas the other ends up as a disc-dominated galaxy. The origins of this dichotomy are the angular momentum of the accreted gas, and whether it adds constructively to the initial disc angular momentum. The build-up of this extended disc leads to a rapid lowering of $\Sigma_{\rm SFR}$ by over two orders of magnitude with $\Sigma_{\rm SFR} \lesssim 0.1~\mathrm{M_\odot}~{\rm yr}^{-1}~{\rm kpc}^{-2}$, in remarkable agreement with what is derived from Milky Way stellar populations in the solar neighbourhood. As a consequence, supernovae explosions are spread out and cannot launch galactic outflows anymore, allowing for the persistence of a thin, gently star forming, extended disc.

We introduce the "Engineering Dwarfs at Galaxy formation's Edge" (EDGE) project to study the cosmological formation and evolution of the smallest galaxies in the Universe. In this first paper, we explore the effects of resolution and sub-grid physics on a single low mass halo ($M_{\rm halo}=10^{9}~M_\odot$), simulated to redshift $z=0$ at a mass and spatial resolution of $\sim 20~M_\odot$ and $\sim 3$ pc. We consider different star formation prescriptions, supernova feedback strengths and on-the-fly radiative transfer (RT). We show that RT changes the mode of galactic self-regulation at this halo mass, suppressing star formation by causing the interstellar and circumgalactic gas to remain predominantly warm ($\sim 10^4$ K) even before cosmic reionisation. By contrast, without RT, star formation regulation occurs only through starbursts and their associated vigorous galactic outflows. In spite of this difference, the entire simulation suite (with the exception of models without any feedback) matches observed dwarf galaxy sizes, velocity dispersions, $V$-band magnitudes and dynamical mass-to-light-ratios. This is because such structural scaling relations are predominantly set by the host dark matter halo, with the remaining model-to-model variation being smaller than the observational scatter. We find that only the stellar mass-metallicity relation differentiates the galaxy formation models. Explosive feedback ejects more metals from the dwarf, leading to a lower metallicity at a fixed stellar mass. We conclude that the stellar mass-metallicity relation of the very smallest galaxies provides a unique constraint on galaxy formation physics.

We demonstrate how the least luminous galaxies in the Universe, ultra-faint dwarf galaxies, are sensitive to their dynamical mass at the time of cosmic reionization. We select a low-mass ($\sim \text{1.5} \times 10^{9} \, \text{M}_{\odot}$) dark matter halo from a cosmological volume, and perform zoom hydrodynamical simulations with multiple alternative histories using "genetically modified" initial conditions. Earlier forming ultra-faints have higher stellar mass today, due to a longer period of star formation before their quenching by reionization. Our histories all converge to the same final dynamical mass, demonstrating the existence of extended scatter ($\geq$ 1 dex) in stellar masses at fixed halo mass due to the diversity of possible histories. One of our variants builds less than 2 % of its final dynamical mass before reionization, rapidly quenching in-situ star formation. The bulk of its final stellar mass is later grown by dry mergers, depositing stars in the galaxy's outskirts and hence expanding its effective radius. This mechanism constitutes a new formation scenario for highly diffuse ($\text{r}_{1 /2} \sim 820 \, \text{pc}$, $\sim 32 \, \text{mag arcsec}^2$), metal-poor ($\big[ \mathrm{Fe}\, / \mathrm{H} \big]= -2.9$), ultra-faint ($\mathcal{M}_V= -5.7$) dwarf galaxies within the reach of next-generation low surface brightness surveys.

We have investigated the possibilities to quantify how much stars move in the Milky Way stellar disk due to diffuse processes (i.e. so called blurring) and due to influences from spiral arms and the bar (i.e. so called churning). To this end we assume that it is possible to infer the formation radius of a star if we know their elemental abundances and age as well as the metallicity profile of the interstellar medium at the time of the formation of the star. Using this information, coupled with orbital information derived from Gaia DR2 data and radial velocities from large spectroscopic surveys, we show that it is possible to isolate stellar samples such that we can start to quantify how important the role of churning is. We use data from APOGEE DR14, parallaxes from Gaia and stellar ages based on C and N elemental abundances in the stars. In our sample, we find that about half of the stars have experienced some sort of radial migration (based solely on their orbital properties), 10 % have likely have suffered only from churning, whilst a modest 5-7 % of stars have never experienced either churning or blurring making them ideal tracers of the original properties of the cool stellar disk. Our investigation shows that it is possible to put up a framework where we can begin to quantify churning and blurring an important. Important aspects for future work would include to investigate how selection effects should be accounted for.

Observations find a median star formation efficiency per free-fall time in Milky Way Giant Molecular Clouds (GMCs) on the order of $\epsilon_{\rm ff}\sim 1\%$ with dispersions of $\sim0.5\,{\rm dex}$. The origin of this scatter in $\epsilon_{\rm ff}$ is still debated and difficult to reproduce with analytical models. We track the formation, evolution and destruction of GMCs in a hydrodynamical simulation of a Milky Way-like galaxy and by deriving cloud properties in an observationally motivated way, measure the distribution of star formation efficiencies which are in excellent agreement with observations. We find no significant link between $\epsilon_{\rm ff}$ and any measured global property of GMCs (e.g. gas mass, velocity dispersion). Instead, a wide range of efficiencies exist in the entire parameter space. From the cloud evolutionary tracks, we find that each cloud follow a \emph{unique} evolutionary path which gives rise to wide diversity in all properties. We argue that it is this diversity in cloud properties, above all else, that results in the dispersion of $\epsilon_{\rm ff}$.

The signatures of the formation and evolution of a galaxy are imprinted in its stars. Their velocities, ages, and chemical compositions present major constraints on models of galaxy formation, and on various processes such as the gas inflows and outflows, the accretion of cold gas, radial migration, and the variability of star formation activity. Understanding the evolution of the Milky Way requires large observational datasets of stars via which these quantities can be determined accurately. This is the science driver of the 4MOST MIlky way Disc And BuLgE High-Resolution (4MIDABLE-HR) survey: to obtain high-resolution spectra at $R \sim 20\,000$ and to provide detailed elemental abundances for large samples of stars in the Galactic disc and bulge. High data quality will allow us to provide accurate spectroscopic diagnostics of two million stellar spectra: precise radial velocities; rotation; abundances of many elements, including those that are currently only accessible in the optical, such as Li, s-, and r-process; and multi-epoch spectra for a sub-sample of stars. Synergies with complementary missions like Gaia and TESS will provide masses, stellar ages and multiplicity, forming a multi-dimensional dataset that will allow us to explore and constrain the origin and structure of the Milky Way.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

Modelling the propagation of supernova (SN) bubbles, in terms of energy, momentum and spatial extent, is critical for simulations of galaxy evolution which do not capture these scales. To date, small scale models of SN feedback predict that the evolution of above-mentioned quantities can be solely parameterised by average quantities of the surrounding gas, such as density. However, most of these studies neglect the turbulent motions of this medium. In this paper, we study the propagation and evolution of SNe in turbulent environments. We confirm that the time evolution of injected energy and momentum can be characterised by the average density. However, the details of the density structure of the interstellar medium play a crucial role in the spatial extent of the bubble, even at a given average density. We demonstrate that spherically symmetric models of SN bubbles do not model well their spatial extent, and therefore cannot not be used to design sub-grid models of SNe feedback at galactic and cosmological scales.

The physical origin of enhanced star formation activity in interacting galaxies remains an open question. Knowing whether starbursts are triggered by an increase of the quantity of dense gas or an increase of the star formation efficiency would improve our understanding of galaxy evolution and allow to transpose the results obtained in the local Universe to high redshift galaxies. In this paper, we analyze a parsec-resolution simulation of an Antennae-like model of interacting galaxies. We find that the interplay of physical processes has complex and important variations in time and space, through different combinations of mechanisms like tides, shear and turbulence. These can have similar imprints on observables like depletion time and CO emission. The densest gas only constitutes the tail of the density distribution of some clouds, but exists in large excess in others. The super-linearity of the star formation rate dependence on gas density implies that this excess translates into a reduction of depletion times, and thus leads to a deviation from the classical star formation regime, visible up to galactic scales. These clouds are found in all parts of the galaxies, but their number density varies from one region to the next, due to different cloud assembly mechanisms. Therefore, the dependence of cloud and star formation-related quantities (like CO flux and depletion time) on the scale at which they are measured also varies across the galaxies. We find that the $\alpha_{\rm CO}$ conversion factor between the CO luminosity and molecular gas mass has even stronger spatial than temporal variations in a system like the Antennae. These results raise a number of cautionary notes for the interpretation of observations of unresolved star-forming regions, but also predict that the diversity of environments for star formation will be better captured by the future generations of instruments.

We use state-of-the-art cosmological zoom simulations to explore the distribution of neutral gas in and around galaxies that gives rise to high column density \ion{H}{i} \mbox{Ly-$\alpha$} absorption (formally, sub-DLAs and DLAs) in the spectra of background quasars. Previous cosmological hydrodynamic simulations under-predict the mean projected separations $(b)$ of these absorbers relative to the host, and invoke selection effects to bridge the gap with observations. On the other hand, single lines of sight (LOS) in absorption cannot uniquely constrain the galactic origin. Our simulations match all observational data, with DLA and sub-DLA LOS existing over the entire probed parameter space ($-4\lesssim $[M/H]$\lesssim 0.5$, $b<50$ kpc) at all redshifts ($z\sim 0.4 - 3.0$). We demonstrate how the existence of DLA LOS at $b\gtrsim 20-30$ kpc from a massive host galaxy require high numerical resolution, and that these LOS are associated with dwarf satellites in the main halo, stripped metal-rich gas and outflows. Separating the galaxy into interstellar ("\ion{H}{i} disc") and circumgalactic ("halo") components, we find that both components significantly contribute to damped \ion{H}{i} absorption LOS. Above the sub-DLA (DLA) limits, the disc and halo contribute with $\sim 60 (80)$ and $\sim 40 (20)$ per cent, respectively. Our simulations confirm analytical model-predictions of the DLA-distribution at $z\lesssim 1$. At high redshift ($z\sim 2-3$) sub-DLA and DLAs occupy similar spatial scales, but on average separate by a factor of two by $z\sim 0.5$. On whether sub-DLA and DLA LOS sample different stellar-mass galaxies, such a correlation can be driven by a differential covering-fraction of sub-DLA to DLA LOS with stellar mass. This preferentially selects sub-DLA LOS in more massive galaxies in the low-$z$ universe.

We explore the multiphase structure of the circumgalactic medium (CGM) probed by synthetic spectra through a cosmological zoom-in galaxy formation simulation. We employ a Bayesian method for modelling a combination of absorption lines to derive physical properties of absorbers with a formal treatment of detections, including saturated systems, and non-detections in a uniform manner. We find that in the lines of sight passing through localized density structures, absorption lines of low, intermediate and high ions are present in the spectrum and overlap in velocity space. Low, intermediate and high ions can be combined to derive the mass-weighted properties of a density-varying peak, although the ions are not co-spatial within the structure. By contrast, lines of sight that go through the hot halo only exhibit detectable HI and high ions. In such lines of sight, the absorption lines are typically broad due to the complex velocity fields across the entire halo. We show that the derived gas density, temperature, and metallicity match closely the corresponding HI mass-weighted averages along the LOS. We also show that when the data quality allows, our Bayesian technique allows one to recover the underlying physical properties of LOS by incorporating both detections and non-detections. It is especially useful to include non-detections, of species such as NV or NeVIII, when the number of detections of strong absorbers, such as HI and OVI, is smaller than the number of model parameters (density, temperature, and metallicity).

Using hydrodynamical simulations of entire galactic discs similar to the Milky Way, reaching 4.6pc resolution, we study the origins of observed physical properties of giant molecular clouds (GMCs). We find that efficient stellar feedback is a necessary ingredient in order to develop a realistic interstellar medium (ISM), leading to molecular cloud masses, sizes, velocity dispersions and virial parameters in excellent agreement with Milky Way observations. GMC scaling relations observed in the Milky Way, such as the mass-size ($M$--$R$), velocity dispersion-size ($\sigma$--$R$), and the $\sigma$--$R\Sigma$ relations, are reproduced in a feedback driven ISM when observed in projection, with $M\propto R^{2.3}$ and $\sigma\propto R^{0.56}$. When analysed in 3D, GMC scaling relations steepen significantly, indicating potential limitations of our understanding of molecular cloud 3D structure from observations. Furthermore, we demonstrate how a GMC population's underlying distribution of virial parameters can strongly influence the scatter in derived scaling relations. Finally, we show that GMCs with nearly identical global properties exist in different evolutionary stages, where a majority of clouds being either gravitationally bound or expanding, but with a significant fraction being compressed by external ISM pressure, at all times.

We present a model for the concurrent formation of globular clusters (GCs) and supermassive stars (SMSs, $>10^3\,{\rm M}_\odot$) to address the origin of the HeCNONaMgAl abundance anomalies in GCs. GCs form in converging gas flows and accumulate low-angular momentum gas, which accretes onto protostars. This leads to an adiabatic contraction of the cluster and an increase of the stellar collision rate. A SMS can form via runaway collisions if the cluster reaches sufficiently high density before two-body relaxation halts the contraction. This condition is met if the number of stars $\gtrsim10^6$ and the gas accretion rate $\gtrsim10^5\,{\rm M}_\odot$/Myr, reminiscent of GC formation in high gas-density environments, such as -- but not restricted to -- the early Universe. The strong SMS wind mixes with the inflowing pristine gas, such that the protostars accrete diluted hot-hydrogen burning yields of the SMS. Because of continuous rejuvenation, the amount of processed material liberated by the SMS can be an order of magnitude higher than its maximum mass. This `conveyor-belt' production of hot-hydrogen burning products provides a solution to the mass budget problem that plagues other scenarios. Additionally, the liberated material is mildly enriched in helium and relatively rich in other hot-hydrogen burning products, in agreement with abundances of GCs today. Finally, we find a super-linear scaling between the amount of processed material and cluster mass, providing an explanation for the observed increase of the fraction of processed material with GC mass. We discuss open questions of this new GC enrichment scenario and propose observational tests.

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.

We discuss some of the key open questions regarding the formation and evolution of globular clusters (GCs) during galaxy formation and assembly within a cosmological framework. The current state-of-the-art for both observations and simulations is described, and we briefly mention directions for future research. The oldest GCs have ages $\ge$ 12.5 Gyr and formed around the time of reionisation. Resolved colour-magnitude diagrams of Milky Way GCs and direct imaging of lensed proto-GCs at z $\sim$ 6 with JWST promise further insight. Globular clusters are known to host multiple populations of stars with variations in their chemical abundances. Recently, such multiple populations have been detected in $\sim$2 Gyr old compact, massive star clusters. This suggests a common, single pathway for the formation of GCs at high and low redshift. The shape of the initial mass function for GCs remains unknown, however for massive galaxies a power-law mass function is favoured. Significant progress has been made recently modelling GC formation in the context of galaxy formation, with success in reproducing many of the observed GC-galaxy scaling relations.

We fit the rotation curves of isolated dwarf galaxies to directly measure the stellar mass-halo mass relation ($M_*-M_{200}$) over the mass range $5 \times 10^5 < M_{*}/{\rm M}_\odot < 10^{8}$. By accounting for cusp-core transformations due to stellar feedback, we find a monotonic relation with little scatter. Such monotonicity implies that abundance matching should yield a similar $M_*-M_{200}$ if the cosmological model is correct. Using the 'field galaxy' stellar mass function from the Sloan Digital Sky Survey (SDSS) and the halo mass function from the $\Lambda$ Cold Dark Matter Bolshoi simulation, we find remarkable agreement between the two. This holds down to $M_{200} \sim 5 \times 10^9$M$_\odot$, and to $M_{200} \sim 5 \times 10^8$M$_\odot$ if we assume a power law extrapolation of the SDSS stellar mass function below $M_* \sim 10^7$M$_\odot$. However, if instead of SDSS we use the stellar mass function of nearby galaxy groups, then the agreement is poor. This occurs because the group stellar mass function is shallower than that of the field below $M_* \sim 10^9$M$_\odot$, recovering the familiar 'missing satellites' and 'too big to fail' problems. Our result demonstrates that both problems are confined to group environments and must, therefore, owe to 'galaxy formation physics' rather than exotic cosmology. Finally, we repeat our analysis for a $\Lambda$ Warm Dark Matter cosmology, finding that it fails at 68% confidence for a thermal relic mass of $m_{\rm WDM} < 1.25$keV, and $m_{\rm WDM} < 2$keV if we use the power law extrapolation of SDSS. We conclude by making a number of predictions for future surveys based on these results.

We study the impact of stellar feedback in shaping the density and velocity structure of neutral hydrogen (HI) in disc galaxies. For our analysis, we carry out $\sim 4.6$pc resolution $N$-body+adaptive mesh refinement (AMR) hydrodynamic simulations of isolated galaxies, set up to mimic a Milky Way (MW), and a Large and Small Magellanic Cloud (LMC, SMC). We quantify the density and velocity structure of the interstellar medium using power spectra and compare the simulated galaxies to observed HI in local spiral galaxies from THINGS (The HI Nearby Galaxy Survey). Our models with stellar feedback give an excellent match to the observed THINGS HI density power spectra. We find that kinetic energy power spectra in feedback regulated galaxies, regardless of galaxy mass and size, show scalings in excellent agreement with super-sonic turbulence ($E(k)\propto k^{-2}$) on scales below the thickness of the HI layer. We show that feedback influences the gas density field, and drives gas turbulence, up to large (kpc) scales. This is in stark contrast to density fields generated by large scale gravity-only driven turbulence. We conclude that the neutral gas content of galaxies carries signatures of stellar feedback on all scales.

We present a cosmological zoom-in simulation of a Milky Way-like galaxy used to explore the formation and evolution of star clusters. We investigate in particular the origin of the bimodality observed in the colour and metallicity of globular clusters, and the environmental evolution through cosmic times in the form of tidal tensors. Our results self-consistently confirm previous findings that the blue, metal-poor clusters form in satellite galaxies which are accreted onto the Milky Way, while the red, metal-rich clusters form mostly in situ or, to a lower extent in massive, self-enriched galaxies merging with the Milky Way. By monitoring the tidal fields these populations experience, we find that clusters formed in situ (generally centrally concentrated) feel significantly stronger tides than the accreted ones, both in the present-day, and when averaged over their entire life. Furthermore, we note that the tidal field experienced by Milky Way clusters is significantly weaker in the past than at present-day, confirming that it is unlikely that a power-law cluster initial mass function like that of young massive clusters, is transformed into the observed peaked distribution in the Milky Way with relaxation-driven evaporation in a tidal field.

Galaxy mergers are believed to trigger strong starbursts. This is well assessed by observations in the local Universe. However the efficiency of this mechanism has poorly been tested so far for high redshift, actively star forming, galaxies. We present a suite of pc-resolution hydrodynamical numerical simulations to compare the star formation process along a merging sequence of high and low z galaxies, by varying the gas mass fraction between the two models. We show that, for the same orbit, high-redshift gas-rich mergers are less efficient than low-redshift ones at producing starbursts: the star formation rate excess induced by the merger and its duration are both around 10 times lower than in the low gas fraction case. The mechanisms that account for the star formation triggering at low redshift - the increased compressive turbulence, gas fragmentation, and central gas inflows - are only mildly, if not at all, enhanced for high gas fraction galaxy encounters. Furthermore, we show that the strong stellar feedback from the initially high star formation rate in high redshift galaxies does not prevent an increase of the star formation during the merger. Our results are consistent with the observed increase of the number of major mergers with increasing redshift being faster than the respective increase in the number of starburst galaxies.

Feedback from massive stars is thought to play an important role in the evolution of molecular clouds. In this work we analyse the effects of stellar winds and supernovae (SNe) in the evolution of two massive ($\sim 10^6\,M_\odot$) giant molecular clouds (GMCs): one gravitationally bound collapsing cloud and one unbound cloud undergoing disruption by galactic shear. These two clouds have been extracted from a large scale galaxy model and are re-simulated at a spatial resolution of $\sim 0.01$ pc, including feedback from winds, SNe, and the combined effect of both. We find that stellar winds stop accretion of gas onto sink particles, and can also trigger star formation in the shells formed by the winds, although the overall effect is to reduce the global star formation rate of both clouds. Furthermore, we observe that winds tend to escape through the corridors of diffuse gas. The effect of SNe is not so prominent and the star formation rate is similar to models neglecting stellar feedback. We find that most of the energy injected by the SNe is radiated away, but overdense areas are created by multiple and concurrent SN events especially in the most virialised cloud. Our results suggest that the impact of stellar feedback is sensitive to the morphology of star forming clouds, which is set by large scale galactic flows, being of greater importance in clouds undergoing gravitational collapse.

We present new Atacama Large Millimeter/Submillimeter Array (ALMA) 850um continuum observations of the original Lyman-alpha Blob (LAB) in the SSA22 field at z=3.1 (SSA22-LAB01). The ALMA map resolves the previously identified submillimeter source into three components with total flux density S_850 = 1.68+/-0.06 mJy, corresponding to a star formation rate of ~150 M_sun/yr. The submillimeter sources are associated with several faint (m~27 mag) rest-frame ultraviolet sources identified in Hubble Space Telescope Imaging Spectrograph (STIS) clear filter imaging (~5850A). One of these companions is spectroscopically confirmed with Keck MOSFIRE to lie within 20 projected kpc and 250 km/s of one of the ALMA components. We postulate that some of these STIS sources represent a population of low-mass star-forming satellites surrounding the central submillimeter sources, potentially contributing to their growth and activity through accretion. Using a high resolution cosmological zoom simulation of a 10^13 M_sun halo at z=3, including stellar, dust and Ly-alpha radiative transfer, we can model the ALMA+STIS observations and demonstrate that Ly-alpha photons escaping from the central submillimeter sources are expected to resonantly scatter in neutral hydrogen, the majority of which is predicted to be associated with halo substructure. We show how this process gives rise to extended Ly-alpha emission with similar surface brightness and morphology to observed giant LABs.

The shape and diversity of dwarf galaxy rotation curves is at apparent odds with dark matter halos in a $\Lambda$ Cold Dark Matter ($\Lambda$CDM) cosmology. We use mock data from isolated dwarf galaxy simulations to show that this owes to three main effects. Firstly, stellar feedback heats dark matter, leading to a 'coreNFW' dark matter density profile with a slowly rising rotation curve. Secondly, if close to a recent starburst, large HI bubbles push the rotation curve out of equilibrium, deforming the rotation curve shape. Thirdly, when galaxies are viewed near face-on, their best fit inclination is biased high. This can lead to a very shallow rotation curve that falsely implies a large dark matter core. All three problems can be avoided, however, by a combination of improved mass models and a careful selection of target galaxies. Fitting our coreNFW model to mock rotation curve data, we show that we can recover the rotation curve shape, dark matter halo mass $M_{200}$ and concentration parameter $c$ within our quoted uncertainties. We fit our coreNFW model to real data for four isolated dwarf irregulars, chosen to span a wide range of rotation curve shapes. We obtain an excellent fit for NGC 6822 and WLM, with tight constraints on $M_{200}$, and $c$ consistent with $\Lambda$CDM. However, IC 1613 and DDO 101 give a poor fit. For IC 1613, we show that this owes to disequilibria and its uncertain inclination $i$; for DDO 101, it owes to its uncertain distance $D$. If we assume $i_{\rm IC1613} \sim 15^\circ$ and $D_{\rm DDO101} \sim 12$ Mpc, consistent with current uncertainties, we are able to fit both galaxies very well. We conclude that $\Lambda$CDM appears to give an excellent match to dwarf galaxy rotation curves.

We use high resolution simulations of isolated dwarf galaxies to study the physics of dark matter cusp-core transformations at the edge of galaxy formation: M200 = 10^7 - 10^9 Msun. We work at a resolution (~4 pc minimum cell size; ~250 Msun per particle) at which the impact from individual supernovae explosions can be resolved, becoming insensitive to even large changes in our numerical 'sub-grid' parameters. We find that our dwarf galaxies give a remarkable match to the stellar light profile; star formation history; metallicity distribution function; and star/gas kinematics of isolated dwarf irregular galaxies. Our key result is that dark matter cores of size comparable to the stellar half mass radius (r_1/2) always form if star formation proceeds for long enough. Cores fully form in less than 4 Gyrs for the M200 = 10^8 Msun and 14 Gyrs for the 10^9 Msun dwarf. We provide a convenient two parameter 'coreNFW' fitting function that captures this dark matter core growth as a function of star formation time and the projected stellar half mass radius. Our results have several implications: (i) we make a strong prediction that if LCDM is correct, then 'pristine' dark matter cusps will be found either in systems that have truncated star formation and/or at radii r > r_1/2; (ii) complete core formation lowers the projected velocity dispersion at r_1/2 by a factor ~2, which is sufficient to fully explain the 'too big to fail problem'; and (iii) cored dwarfs will be much more susceptible to tides, leading to a dramatic scouring of the subhalo mass function inside galaxies and groups.

We analyze circumgalactic medium (CGM) in a suite of high-resolution cosmological re-simulations of a Milky-Way size galaxy and show that CGM properties are quite sensitive to details of star formation--feedback loop modelling. The simulation that produces a realistic late-type galaxy, fails to reproduce existing observations of the CGM. In contrast, simulation that does not produce a realistic galaxy has the predicted CGM in better agreement with observations. This illustrates that properties of galaxies and properties of their CGM provide strong ${\it complementary}$ constraints on the processes governing galaxy formation. Our simulations predict that column density profiles of ions are well described by an exponential function of projected distance $d$: $N \propto e^{-d/h_s}$. Simulations thus indicate that the sharp drop in absorber detections at larger distances in observations does not correspond to a "boundary" of an ion, but reflects the underlying steep exponential column density profile. Furthermore, we find that ionization energy of ions is tightly correlated with the scale height $h_s$: $h_s \propto E_{\rm ion}^{0.74}$. At $z \approx 0$, warm gas traced by low-ionization species (e.g., Mg II and C IV) has $ h_s \approx 0.03-0.07 R_{\rm vir}$, while higher ionization species (O VI and Ne VIII) have $h_s \approx 0.32-0.45R_{\rm vir}$. Finally, the scale heights of ions in our simulations evolve slower than the virial radius for $z\leq 2$, but similarly to the halo scale radius, $r_s$. Thus, we suggest that the column density profiles of galaxies at different redshifts should be scaled by $r_s$ rather than the halo virial radius.

The strong time-dependence of the dynamics of galactic bars yields a complex and rapidly evolving distribution of dense gas and star forming regions. Although bars mainly host regions void of any star formation activity, their extremities can gather the physical conditions for the formation of molecular complexes and mini-starbursts. Using a sub-parsec resolution hydrodynamical simulation of a Milky Way-like galaxy, we probe these conditions to explore how and where bar (hydro-)dynamics favours the formation or destruction of molecular clouds and stars. The interplay between the kpc-scale dynamics (gas flows, shear) and the parsec-scale (turbulence) is key to this problem. We find a strong dichotomy between the leading and trailing sides of the bar, in term of cloud fragmentation and in the age distribution of the young stars. After orbiting along the bar edge, these young structures slow down at the extremities of the bar, where orbital crowding increases the probability of cloud-cloud collision. We find that such events increase the Mach number of the cloud, leading to an enhanced star formation efficiency and finally the formation of massive stellar associations, in a fashion similar to galaxy-galaxy interactions. We highlight the role of bar dynamics in decoupling young stars from the clouds in which they form, and discuss the implications on the injection of feedback into the interstellar medium, in particular in the context of galaxy formation.

We use cosmological zoom-in simulations of galaxy formation in a Milky Way (MW)-sized halo started from identical initial conditions to investigate the evolution of galaxy sizes, baryon fractions, morphologies and angular momenta in runs with different parameters of the star formation--feedback cycle. Our fiducial model with a high local star formation efficiency, which results in efficient feedback, produces a realistic late-type galaxy that matches the evolution of basic properties of late-type galaxies: stellar mass, disk size, morphology dominated by a kinematically cold disk, stellar and gas surface density profiles, and specific angular momentum. We argue that feedback's role in this success is twofold: (1) removal of low-angular momentum gas and (2) maintaining a low disk-to-halo mass fraction which suppresses disk instabilities that lead to angular momentum redistribution and a central concentration of baryons. However, our model with a low local star formation efficiency, but large energy input per supernova, chosen to produce a galaxy with a similar star formation history as our fiducial model, leads to a highly irregular galaxy with no kinematically cold component, overly extended stellar distribution and low angular momentum. This indicates that only when feedback is allowed to become vigorous via locally efficient star formation in dense cold gas, resulting galaxy sizes, gas/stellar surface density profiles and stellar disk angular momenta agree with observed $z=0$ galaxies.

Radiation feedback is typically implemented using subgrid recipes in hydrodynamical simulations of galaxies. Very little work has so far been performed using radiation-hydrodynamics (RHD), and there is no consensus on the importance of radiation feedback in galaxy evolution. We present RHD simulations of isolated galaxy disks of different masses with a resolution of 18 pc. Besides accounting for supernova feedback, our simulations are the first galaxy-scale simulations to include RHD treatments of photo-ionisation heating and radiation pressure, from both direct optical/UV radiation and multi-scattered, re-processed infrared (IR) radiation. Photo-heating smooths and thickens the disks and suppresses star formation about as much as the inclusion of ("thermal dump") supernova feedback does. These effects decrease with galaxy mass and are mainly due to the prevention of the formation of dense clouds, as opposed to their destruction. Radiation pressure, whether from direct or IR radiation, has little effect, but for the IR radiation we show that its impact is limited by our inability to resolve the high optical depths for which multi-scattering becomes important. While artificially boosting the IR optical depths does reduce the star formation, it does so by smoothing the gas rather than by generating stronger outflows. We conclude that although higher-resolution simulations, and potentially also different supernova implementations, are needed for confirmation, our findings suggest that radiation feedback is more gentle and less effective than is often assumed in subgrid prescriptions.

Several scenarios have been proposed to explain the presence of multiple stellar populations in globular clusters. Many of them invoke multiple generations of stars to explain the observed chemical abundance anomalies, but it has also been suggested that self-enrichment could occur via accretion of ejecta from massive stars onto the circumstellar disc of low-mass pre-main sequence stars. These scenarios imply different initial conditions for the kinematics of the various stellar populations. Given some net angular momentum initially, models for which a second generation forms from gas that collects in a cooling flow into the core of the cluster predict an initially larger rotational amplitude for the polluted stars compared to the pristine stars. This is opposite to what is expected from the accretion model, where the polluted stars are the ones crossing the core and are on preferentially radial (low-angular momentum) orbits, such that their rotational amplitude is lower. Here we present the results of a suite of $N$-body simulations with initial conditions chosen to capture the distinct kinematic properties of these pollution scenarios. We show that initial differences in the kinematics of polluted and pristine stars can survive to the present epoch in the outer parts of a large fraction of Galactic globular clusters. The differential rotation of pristine and polluted stars is identified as a unique kinematic signature that could allow us to distinguish between various scenarios, while other kinematic imprints are generally very similar from one scenario to the other.

We describe a NOVel form of Adaptive softening (NovA) for collisionless $N$-body simulations, implemented in the Ramses adaptive mesh refinement code. We introduce a refinement criterion that the particle distribution within each cell be sufficiently isotropic, as measured by its moment of inertia tensor. In this way, collapse is only refined if it occurs along all three axes, ensuring that the softening $\epsilon$ is always of order twice the largest inter-particle spacing in a cell. This more conservative force softening criterion is designed to minimise spurious two-body effects, while maintaining high force resolution in collapsed regions of the flow. We test NovA using an antisymmetric perturbed plane wave collapse (`Valinia' test) before applying it to warm dark matter (WDM) simulations. For the Valinia test, we show that -- unlike the standard $N$-body method -- NovA produces no numerical fragmentation while still being able to correctly capture fine caustics and shells around the collapsing regions. For the WDM simulations, we find that NovA converges significantly more rapidly than standard $N$-body, producing little or no spurious halos on small scales. We show, however, that determining whether or not halos exist below the free streaming mass $M_{\rm fs}$ is complicated by the fact that our halo finder (AHF) likely incorrectly labels some caustics and criss-crossing filaments as halos, while one or two particularly massive filaments appear to fragment in any version of NovA where refinement is allowed. Such massive filaments may be physically unstable to collapse, as is the case for infinite, static, self-gravitating cylinders. We will use NovA in forthcoming papers to study the issue of halo formation below $M_{\rm fs}$; filament stability; and to obtain new constraints on the temperature of dark matter.

Gravitational instabilities play an important role in galaxy evolution and in shaping the interstellar medium (ISM). The ISM is observed to be highly turbulent, meaning that observables like the gas surface density and velocity dispersion depend on the size of the region over which they are measured. In this work we investigate, using simulations of Milky Way-like disc galaxies with a resolution of $\sim 9$ pc, the nature of turbulence in the ISM and how this affects the gravitational stability of galaxies. By accounting for the measured average turbulent scalings of the density and velocity fields in the stability analysis, we can more robustly characterize the average level of stability of the galaxies as a function of scale, and in a straightforward manner identify scales prone to fragmentation. Furthermore, we find that the stability of a disc with feedback-driven turbulence can be well described by a "Toomre-like" $Q$ stability criterion on all scales, whereas the classical $Q$ can formally lose its meaning on small scales if violent disc instabilities occur in models lacking pressure support from stellar feedback.

We investigate the star formation-feedback cycle in cosmological galaxy formation simulations, focusing on progenitors of Milky Way (MW)-sized galaxies. We find that in order to reproduce key properties of the MW progenitors, such as semi-empirically derived star formation histories and the shape of rotation curves, our implementation of star formation and stellar feedback requires 1) a combination of local early momentum feedback via radiation pressure and stellar winds and subsequent efficient supernovae feedback, and 2) efficacy of feedback that results in self-regulation of the global star formation rate on kiloparsec scales. We show that such feedback-driven self-regulation is achieved globally for a local star formation efficiency per free fall time of $\epsilon_{\rm ff}\approx 10\%$. Although this value is larger that the $\epsilon_{\rm ff}\sim 1\%$ value usually inferred from the Kennicutt-Schmidt (KS) relation, we show that it is consistent with direct observational estimates of $\epsilon_{\rm ff}$ in molecular clouds. Moreover, we show that simulations with local efficiency of $\epsilon_{\rm ff}\approx 10\%$ reproduce the global observed KS relation. Such simulations also reproduce the cosmic star formation history of the Milky Way sized galaxies and satisfy a number of other observational constraints. Conversely, we find that simulations that a priori assume an inefficient mode of star formation, instead of achieving it via stellar feedback regulation, fail to produce sufficiently vigorous outflows and do not reproduce observations. This illustrates the importance of understanding the complex interplay between star formation and feedback and the detailed processes that contribute to the feedback-regulated formation of galaxies.

Galaxy formation models and simulations rely on various feedback mechanisms to reproduce the observed baryonic scaling relations and galaxy morphologies. Although dwarf galaxy and giant elliptical properties can be explained using feedback from supernova and active galactic nuclei, Milky Way-sized galaxies still represent a challenge to current theories of galaxy formation. In this paper, we explore the possible role of feedback from stellar radiation in regulating the main properties of disk galaxies such as our own Milky Way. We have performed a suite of cosmological simulations of the same $\sim10^{12} {\rm M}_{\odot}$ halo selected based on its rather typical mass accretion history. We have implemented radiative feedback from young stars using a crude model of radiative transfer for ultraviolet (UV) and infrared (IR) radiation. However, the model is realistic enough such that the dust opacity plays a direct role in regulating the efficiency of our feedback mechanism. We have explored various models for the dust opacity, assuming different constant dust temperatures, as well as a varying dust temperature model. We find that while strong radiative feedback appears as a viable mechanism to regulate the stellar mass fraction in massive galaxies, it also prevents the formation of discs with reasonable morphologies. In models with strong stellar radiation feedback, stellar discs are systematically too thick while the gas disc morphology is completely destroyed due to the efficient mixing between the feedback-affected gas and its surroundings. At the resolution of our simulation suite, we find it impossible to preserve spiral disc morphology while at the same time expelling enough baryons to satisfy the abundance matching constraints.

Gravitational instabilities play a primary role in shaping the clumpy structure and powering the star formation activity of gas-rich high-redshift galaxies. Here we analyse the stability of such systems, focusing on the size and mass ranges of unstable regions in the disc. Our analysis takes into account the mass-size and linewidth-size scaling relations observed in molecular gas, originally discovered by Larson. We show that such relations can have a strong impact on the size and mass of star-forming clumps, as well as on the stability properties of the disc at all observable scales, making the classical Toomre parameter a highly unreliable indicator of gravitational instability. For instance, a disc with Q=1 can be far from marginal instability, while a disc with Q<<1 can be marginally unstable. Our work raises an important caveat: if clumpy discs at high redshift have scale-dependent surface densities and velocity dispersions, as implied by the observed clump scaling relations, then we cannot thoroughly understand their stability and star formation properties unless we perform multi-scale observations. This will soon be possible thanks to dedicated ALMA surveys, which will explore the physical properties of super-giant molecular clouds at the peak of cosmic star formation and beyond.

The current velocity of the Smith Cloud indicates that it has undergone at least one passage of the Galactic disc. Using hydrodynamic simulations we examine the present day structure of the Smith Cloud. We find that a dark matter supported cloud is able to reproduce the observed present day neutral hydrogen mass, column density distribution and morphology. In this case the dark matter halo becomes elongated, owing to the tidal interaction with the Galactic disc. Clouds in models neglecting dark matter confinement are destroyed upon disc passage, unless the initial cloud mass is well in excess of what is observed today. We then determine integrated flux upper limits to the gamma-ray emission around such a hypothesised dark matter core in the Smith Cloud. No statistically significant core or extended gamma-ray emission are detected down to a 95% confidence level upper limit of $1.4\times10^{-10}$ ph cm$^{-2}$ s$^{-1}$ in the 1-300 GeV energy range. For the derived distance of 12.4 kpc, the Fermi upper limits set the first tentative constraints on the dark matter cross sections annihilating into $\tau^{+}{\tau}^{-}$ and $b\bar{b}$ for a high-velocity cloud.

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

We present results from high-resolution hydrodynamic simulations of isolated SMC- and Milky Way-sized galaxies that include a model for feedback from galactic cosmic rays (CRs). We find that CRs are naturally able to drive winds with mass loading factors of up to ~10 in dwarf systems. The scaling of the mass loading factor with circular velocity between the two simulated systems is consistent with \propto v_c^{1-2} required to reproduce the faint end of the galaxy luminosity function. In addition, simulations with CR feedback reproduce both the normalization and the slope of the observed trend of wind velocity with galaxy circular velocity. We find that winds in simulations with CR feedback exhibit qualitatively different properties compared to SN driven winds, where most of the acceleration happens violently in situ near star forming sites. In contrast, the CR-driven winds are accelerated gently by the large-scale pressure gradient established by CRs diffusing from the star-forming galaxy disk out into the halo. The CR-driven winds also exhibit much cooler temperatures and, in the SMC-sized system, warm (T~10^4 K) gas dominates the outflow. The prevalence of warm gas in such outflows may provide a clue as to the origin of ubiquitous warm gas in the gaseous halos of galaxies detected via absorption lines in quasar spectra.

Stellar feedback plays a key role in galaxy formation by regulating star formation, driving interstellar turbulence and generating galactic scale outflows. Although modern simulations of galaxy formation can resolve scales of 10-100 pc, star formation and feedback operate on smaller, "subgrid" scales. Great care should therefore be taken in order to properly account for the effect of feedback on global galaxy evolution. We investigate the momentum and energy budget of feedback during different stages of stellar evolution, and study its impact on the interstellar medium using simulations of local star forming regions and galactic disks at the resolution affordable in modern cosmological zoom-in simulations. In particular, we present a novel subgrid model for the momentum injection due to radiation pressure and stellar winds from massive stars during early, pre-supernova evolutionary stages of young star clusters. Early injection of momentum acts to clear out dense gas in star forming regions, hence limiting star formation. The reduced gas density mitigates radiative losses of thermal feedback energy from subsequent supernova explosions, leading to an increased overall efficiency of stellar feedback. The detailed impact of stellar feedback depends sensitively on the implementation and choice of parameters. Somewhat encouragingly, we find that implementations in which feedback is efficient lead to approximate self-regulation of global star formation efficiency. We compare simulation results using our feedback implementation to other phenomenological feedback methods, where thermal feedback energy is allowed to dissipate over time scales longer than the formal gas cooling time. We find that simulations with maximal momentum injection suppress star formation to a similar degree as is found in simulations adopting adiabatic thermal feedback.

We study the formation of disc galaxies in a fully cosmological framework using adaptive mesh refinement simulations. We perform an extensive parameter study of the main subgrid processes that control how gas is converted into stars and the coupled effect of supernovae feedback. We argue that previous attempts to form disc galaxies have been unsuccessful because of the universal adoption of strong feedback combined with high star formation efficiencies. Unless extreme amounts of energy are injected into the interstellar medium during supernovae events, these star formation parameters result in bulge dominated S0/Sa galaxies as star formation is too efficient at z~3. We show that a low efficiency of star-formation more closely models the subparsec physical processes, especially at high redshift. We highlight the successful formation of extended disc galaxies with scale lengths r_d=4-5 kpc, flat rotation curves and bulge to disc ratios of B/D~1/4. Not only do we resolve the formation of a Milky Way-like spiral galaxy, we also observe the secular evolution of the disc as it forms a pseudo-bulge. The disc properties agree well with observations and are compatible with the photometric and baryonic Tully-Fisher relations, the Kennicutt-Schmidt relation and the observed angular momentum content of spiral galaxies. We conclude that underlying small-scale star formation physics plays a larger role than previously considered in simulations of galaxy formation.

A precise determination of the local dark matter density and an accurate control over the corresponding uncertainties are of paramount importance for Dark Matter (DM) searches. Using very recent high-resolution numerical simulations of a Milky Way like object, we study the systematic uncertainties that affect the determination of the local dark matter density based on dynamical measurements in the Galaxy. In particular, extracting from the simulation with baryons the orientation of the Galactic stellar disk with respect to the DM distribution, we study the DM density for an observer located at $\sim$8 kpc from the Galactic center {\it on the stellar disk}, $\rho_0$. This quantity is found to be always larger than the average density in a spherical shell of same radius $\bar{\rho}_0$, which is the quantity inferred from dynamical measurements in the Galaxy, and to vary in the range $\rho_0/\bar{\rho}_0=1.01-1.41$. This suggests that the actual dark matter density in the solar neighbourhood is on average 21\% larger than the value inferred from most dynamical measurements, and that the associated systematic errors are larger than the statistical errors recently discussed in the literature.

Standard formulations of smoothed particle hydrodynamics (SPH) are unable to resolve mixing at fluid boundaries. We use an error and stability analysis of the generalised SPH equations of motion to prove that this is due to two distinct problems. The first is a leading order error in the momentum equation. This should decrease with increasing neighbour number, but does not because numerical instabilities cause the kernel to be irregularly sampled. We identify two important instabilities: the clumping instability and the banding instability, and we show that both are cured by a suitable choice of kernel. The second problem is the local mixing instability (LMI). This occurs as particles attempt to mix on the kernel scale, but are unable to due to entropy conservation. The result is a pressure discontinuity at boundaries that pushes fluids of different entropy apart. We cure the LMI by using a weighted density estimate that ensures that pressures are single valued throughout the flow. This also gives a better volume estimate for the particles, reducing errors in the continuity and momentum equations. We demonstrate mixing in our new Optimised Smoothed Particle Hydrodynamics (OSPH) scheme using a Kelvin Helmholtz instability (KHI) test with density contrast 1:2, and the 'blob test' - a 1:10 density ratio gas sphere in a wind tunnel - finding excellent agreement between OSPH and Eulerian codes.

We explore the gravitational instability of clumpy and turbulent gas discs, taking into account the Larson-type scaling laws observed in giant molecular clouds (GMCs) and HI, as well as more general scaling relations. This degree of freedom is of special interest in view of the coming high-z ISM surveys, and is thus potentially important for understanding the dynamical effects of turbulence at all epochs of galaxy evolution. Our analysis shows that turbulence has a deep impact on the gravitational instability of the disc. It excites a rich variety of stability regimes, several of which have no classical counterpart. Among other diagnostics, we provide two useful tools for observers and simulators: (1) the stability map of turbulence, which illustrates our stability scenario and relates it to the phenomenology of interstellar turbulence: GMC/HI observations, simulations and models; (2) a Toomre-like stability criterion, $Q\geq\bar{Q}$, which applies to a large class of clumpy/turbulent discs. We make specific predictions about GMC and cold-HI turbulence, and point out the implications of our analysis for high-z galaxy surveys.

Observations of high redshift galaxies have revealed a multitude of large clumpy rapidly star-forming galaxies. Their formation scenario and their link to present day spirals is still unknown. In this Letter we perform adaptive mesh refinement simulations of disk formation in a cosmological context that are unrivalled in terms of mass and spatial resolution. We find that the so called "chain-galaxies" and "clump-clusters" are a natural outcome of early epochs of enhanced gas accretion from cold dense streams as well as tidally and ram-pressured stripped material from minor mergers and satellites. Through interaction with the hot halo gas, this freshly accreted cold gas settles into a large disk-like system, not necessarily aligned to an older stellar component, that undergoes fragmentation and subsequent star formation, forming large clumps in the mass range 10^7-10^9 M_sun. Galaxy formation is a complex process at this important epoch when most of the central baryons are being acquired through a range of different mechanisms - we highlight that a rapid mass loading epoch is required to fuel the fragmentation taking place in the massive arms in the outskirts of extended disks, an accretion mode that occurs naturally in the hierarchical assembly process at early epochs.

Predicting the local flux of dark matter particles is vital for dark matter direct detection experiments. To date, such predictions have been based on simulations that model the dark matter alone. Here we include the influence of the baryonic matter for the first time. We use two different approaches. Firstly, we use dark matter only simulations to estimate the expected merger history for a Milky Way mass galaxy, and then add a thin stellar disc to measure its effect. Secondly, we use three cosmological hydrodynamic simulations of Milky Way mass galaxies. In both cases, we find that a stellar/gas disc at high redshift (z~1) causes merging satellites to be preferentially dragged towards the disc plane. This results in an accreted dark matter disc that contributes ~0.25 - 1 times the non-rotating halo density at the solar position. An associated thick stellar disc forms with the dark disc and shares a similar velocity distribution. If these accreted stars can be separated from those that formed in situ, future astronomical surveys will be able to infer the properties of the dark disc from these stars. The dark disc, unlike dark matter streams, is an equilibrium structure that must exist in disc galaxies that form in a hierarchical cosmology. Its low rotation lag with respect to the Earth significantly boosts WIMP capture in the Earth and Sun, increases the likelihood of direct detection at low recoil energy, boosts the annual modulation signal, and leads to distinct variations in the flux as a function of recoil energy that allow the WIMP mass to be determined (see contribution from T. Bruch this volume).

Observations of turbulent velocity dispersions in the HI component of galactic disks show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high resolution adaptive mesh simulations of gaseous disks embedded within dark matter haloes to explore the roles of cooling, star-formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disk slowly cools until gravitational and thermal instabilities give rise to a multi-phase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large scale driving of global non-axisymmetric modes as well as cloud-cloud tidal interactions and merging. At low star-formation rates these processes alone can explain the observed HI velocity dispersion profiles and the characteristic value of ~10 km/s observed within a wide range of disk galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star-formation rate per unit area >10^-3 M_sun/yr/kpc^2.

The effects of particle discreteness in N-body simulations of Lambda Cold Dark Matter (LambdaCDM) are still an intensively debated issue. In this paper we explore such effects, taking into account the scatter caused by the randomness of the initial conditions, and focusing on the statistical properties of the cosmological density field. For this purpose, we run large sets of LambdaCDM simulations and analyse them using a large variety of diagnostics, including new and powerful wavelet statistics. Among other facts, we point out (1) that dynamical evolution does not propagate discreteness noise up from the small scales at which it is introduced, and (2) that one should aim to satisfy the condition epsilon ~ 2d, where "epsilon" is the force resolution and "d" is the interparticle distance. We clarify what such a condition means, and how to implement it in modern cosmological codes.

In a LCDM cosmology, the Milky Way accretes satellites into the stellar disc. We use cosmological simulations to assess the frequency of near disc plane and higher inclination accretion events, and collisionless simulations of satellite mergers to quantify the final state of the accreted material and the effect on the thin disc. On average, a Milky Way-sized galaxy has 3 subhalos with vmax>80km/s; 7 with vmax>60km/s; and 15 with vmax>40km/s merge at redshift z>1. Assuming isotropic accretion, a third of these merge at an impact angle <20 degrees and are dragged into the disc plane by dynamical friction. Their accreted stars and dark matter settle into a thick disc. The stellar thick disc qualitatively reproduces the observed thick disc at the solar neighbourhood, but is less massive by a factor ~2-10. The dark matter disc contributes 0.25-1 times the halo density at the solar position. Although not likely to be dynamically interesting, the dark disc has important implications for the direct detection of dark matter because of its low velocity with respect to the Earth. Higher inclination encounters (>20 degrees) are twice as likely as low inclination ones. These lead to structures that closely resemble the recently discovered inner/outer stellar halos. They also do more damage to the Milky Way stellar disc creating a more pronounced flare, and warp; both long-lived and consistent with current observations. The most massive mergers (vmax>80km/s) heat the thin disc enough to produce a thick disc. These heated thin disc stars are essential for obtaining a thick disc as massive as that seen in the Milky Way; they likely comprise some ~50-90% of the thick disc stars. The Milky Way thin disc must reform from fresh gas after z=1 [abridged].

Melott has made pioneering studies of the effects of particle discreteness in N-body simulations, a fundamental point that needs careful thought and analysis since all such simulations suffer from numerical noise arising from the use of finite-mass particles. Melott (arXiv:0804.0589) claims that the conclusions of our paper (arXiv:0804.0294) are essentially equivalent to those of his earlier work. Melott is wrong: he has jumped onto one of our conclusions and interpreted that in his own way. Here we point out the whys and the wherefores.

Since its discovery in 1996, the source of the bright H-alpha emission (up to 750 mR) along the Magellanic Stream has remained a mystery. There is no evidence of ionising stars within the HI stream, and the extended hot halo is far too tenuous to drive strong shocks into the clouds. We now present a hydrodynamical model that explains the known properties of the H-alpha emission and provides new insights on the lifetime of the Stream clouds. The upstream clouds are gradually disrupted due to their interaction with the hot halo gas. The clouds that follow plough into gas ablated from the upstream clouds, leading to shock ionisation at the leading edges of the downstream clouds. Since the following clouds also experience ablation, and weaker H-alpha (100-200 mR) is quite extensive, a disruptive cascade must be operating along much of the Stream. In our model, the clouds are evolving on timescales of 100-200 Myr, such that the Stream must be replenished by the Magellanic Clouds at a fairly constant rate. The ablated material falls onto the Galaxy as a warm drizzle which suggests that diffuse ionized gas at 10**4 K may be an important constituent of galactic accretion. The observed HI emission provides a new constraint on the rate of disruption of the Stream and, consequently, the infall rate of metal-poor gas onto the Galaxy. When the ionized component of the Stream is fully accounted for, the rate of gas accretion is 0.4 Msun/yr, roughly twice the rate deduced from HI observations alone.

N-body simulations sample their initial conditions on an initial particle distribution, which for cosmological simulations is usually a glass or grid, whilst a Poisson distribution is used for galaxy models, spherical collapse etc. These pre-initial conditions have inherent correlations, noise due to discreteness and preferential alignments, whilst the glass distribution is poorly defined and computationally expensive to construct. We present a novel particle distribution which can be useful as a pre-initial condition for N-body simulations, using a simple construction based on a ``quaquaversal'' tiling of space. This distribution has little preferred orientation (i.e. is statistically isotropic), has a rapidly vanishing large scale power-spectrum (P(k) ~ k^4), and is trivial to create. It should be particularly useful for warm dark matter and cold collapse simulations.

We have carried out a hydrodynamical code comparison study of interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of the SPH smoothing kernel over which information is not transferred.