The James Webb Space Telescope (JWST) is expected to observe galaxies at $z>10$ that are presently inaccessible. Here, we use a self-consistent empirical model, the UniverseMachine, to generate mock galaxy catalogues and lightcones over the redshift range $z=0-15$. These data include realistic galaxy properties (stellar masses, star formation rates, and UV luminosities), galaxy-halo relationships, and galaxy-galaxy clustering. Mock observables are also provided for different model parameters spanning observational uncertainties at $z<10$. We predict that Cycle 1 JWST surveys will very likely detect galaxies with $M_*>10^7 M_\odot$ and/or $M_{1500}<-17$ out to at least $z\sim 13.5$. Number density uncertainties at $z>12$ expand dramatically, so efforts to detect $z>12$ galaxies will provide the most valuable constraints on galaxy formation models. The faint-end slopes of the stellar mass/luminosity functions at a given mass/luminosity threshold steepen as redshift increases. This is because observable galaxies are hosted by haloes in the exponentially falling regime of the halo mass function at high redshifts. Hence, these faint-end slopes are robustly predicted to become shallower below current observable limits ($M_\ast < 10^7M_\odot$ or $M_\mathrm{1500}>-17$). For reionization models, extrapolating luminosity functions with a constant faint-end slope from $M_{1500}=-17$ down to $M_{1500}=-12$ gives the most reasonable upper limit for the total UV luminosity and cosmic star formation rate up to $z\sim 12$. We compare to three other empirical models and one semi-analytic model, showing that the range of predicted observables from our approach encompasses predictions from other techniques. Public catalogues and lightcones for common fields are available online.

We present the novel wide & deep neural network GalaxyNet, which connects the properties of galaxies and dark matter haloes, and is directly trained on observed galaxy statistics using reinforcement learning. The most important halo properties to predict stellar mass and star formation rate (SFR) are halo mass, growth rate, and scale factor at the time the mass peaks, which results from a feature importance analysis with random forests. We train different models with supervised learning to find the optimal network architecture. GalaxyNet is then trained with a reinforcement learning approach: for a fixed set of weights and biases, we compute the galaxy properties for all haloes and then derive mock statistics (stellar mass functions, cosmic and specific SFRs, quenched fractions, and clustering). Comparing these statistics to observations we get the model loss, which is minimised with particle swarm optimisation. GalaxyNet reproduces the observed data very accurately ($\chi_\mathrm{red}=1.05$), and predicts a stellar-to-halo mass relation with a lower normalisation and shallower low-mass slope at high redshift than empirical models. We find that at low mass, the galaxies with the highest SFRs are satellites, although most satellites are quenched. The normalisation of the instantaneous conversion efficiency increases with redshift, but stays constant above $z\gtrsim0.7$. Finally, we use GalaxyNet to populate a cosmic volume of $(5.9~\mathrm{Gpc})^3$ with galaxies and predict the BAO signal, the bias, and the clustering of active and passive galaxies up to $z=4$, which can be tested with next-generation surveys, such as LSST and Euclid.

Supernovae (SN) generate hot gas in the interstellar medium (ISM), help setting the ISM structure and support the driving of outflows. It is important to resolve the hot gas generation for galaxy formation simulations at solar mass and sub-parsec resolution which realise individual supernova (SN) explosions with ambient densities varying by several orders of magnitude in a realistic multi-phase ISM. We test resolution requirements by simulating SN blast waves at three metallicities ($Z = 0.01, 0.1$ and $1 Z_{\odot}$), six densities and their respective equilibrium chemical compositions ($n=0.001$ cm$^{-3}$ - $100$ cm$^{-3}$), and four mass resolutions ($0.1$ - $100$ M$_{\odot}$), in three dimensions. We include non-equilibrium cooling and chemistry, a homogenous interstellar radiation field, and shielding with a modern pressure-energy smoothed particle hydrodynamics (SPH) method including isotropic thermal conduction and a meshless-finite-mass (MFM) solver. We find stronger resolution requirements for chemistry and hot phase generation than for momentum generation. While at $10$ M$_{\odot}$ the radial momenta at the end of the Sedov phase start converging, the hot phase generation and chemistry require higher resolutions to represent the neutral to ionised hydrogen fraction at the end of the Sedov phase correctly. Thermal conduction typically reduces the hot phase by $0.2$ dex and has little impact on the chemical composition. In general, our $1$, and $0.1$ M$_{\odot}$ results agree well with previous numerical and analytic estimates. We conclude that for the thermal energy injection SN model presented here resolutions higher than $10$ M$_{\odot}$ are required to model the chemistry, momentum and hot phase generation in a multi-phase ISM.

We describe a population of young star clusters (SCs) formed in a hydrodynamical simulation of a gas-rich dwarf galaxy merger resolved with individual massive stars at sub-parsec spatial resolution. The simulation is part of the GRIFFIN (Galaxy Realizations Including Feedback From INdividual massive stars) project. The star formation environment during the simulation spans seven orders of magnitude in gas surface density and thermal pressure, and the global star formation rate surface density ($\Sigma_\mathrm{SFR}$) varies by more than three orders of magnitude during the simulation. Young SCs more massive than $M_{\mathrm{*,cl}}\sim 10^{2.5}\,M_{\odot}$ form along a mass function with a power-law index $\alpha\sim-1.7$ ($\alpha\sim-2$ for $M_{\mathrm{*,cl}}\gtrsim10^{3}\,M_{\odot}$) at all merger phases, while the normalization and the highest SC masses (up to $\sim 10^6 M_{\odot}$) correlate with $\Sigma_\mathrm{SFR}$. The cluster formation efficiency varies from $\Gamma\sim20\%$ in early merger phases to $\Gamma\sim80\%$ at the peak of the starburst and is compared to observations and model predictions. The massive SCs ($\gtrsim10^4\,M_{\odot}$) have sizes and mean surface densities similar to observed young massive SCs. Simulated lower mass clusters appear slightly more concentrated than observed. All SCs form on timescales of a few Myr and lose their gas rapidly resulting in typical stellar age spreads between $\sigma\sim0.1-2$ Myr ($1\sigma$), consistent with observations. The age spreads increase with cluster mass, with the most massive cluster ($\sim10^6\, M_{\odot}$) reaching a spread of $5\, \mathrm{Myr}$ once its hierarchical formation finishes. Our study shows that it is now feasible to investigate the SC population of entire galaxies with novel high-resolution numerical simulations.

We explore the galaxy-galaxy merger rate with the empirical model for galaxy formation, EMERGE. On average, we find that between $2$ per cent and $6$ per cent of massive galaxies ($\log_{10}(m_{*}/M_{\odot}) \geq 10.3$) will experience a major merger per Gyr. Our model predicts galaxy merger rates that do not scale as a power-law with redshift, and exhibit a clear stellar mass and mass-ratio dependence, unlike the underlying halo-halo merger rate. Specifically, major mergers are more frequent at high masses and at low redshift. We show mergers are significant for the stellar mass growth of galaxies $\log_{10}(m_{*}/M_{\odot}) \gtrsim 11.0$. For the most massive galaxies major mergers dominate the accreted mass fraction, contributing as much as $90$ per cent of the total accreted stellar mass. We argue that these phenomena are a direct result of the stellar-to-halo mass relation, which results in massive galaxies having a higher likelihood of experiencing major mergers than low mass galaxies. Furthermore, we find no statistical evidence for mergers as a driving mechanism for quenching, with more than $58$ per cent of the most massive galaxies experiencing no major merger within $2$ Gyr of quenching. Our model produces a galaxy pair fraction consistent with resent observations, exhibiting a form best described by a power-law exponential function. Translating these pair fractions into merger rates results in an over prediction compared to the model intrinsic values. We find the pair fraction can be best mapped to the intrinsic merger rate by adopting a constant observation timescale $T_{obs}=2.56$ Gyr, assuming all observed pairs merge by $z=0$.

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized.

Observational systematics complicate comparisons with theoretical models limiting understanding of galaxy evolution. In particular, different empirical determinations of the stellar mass function imply distinct mappings between the galaxy and halo masses, leading to diverse galaxy evolutionary tracks. Using our state-of-the-art STatistical sEmi-Empirical modeL, STEEL, we show fully self-consistent models capable of generating galaxy growth histories that simultaneously and closely agree with the latest data on satellite richness and star-formation rates at multiple redshifts and environments. Central galaxy histories are generated using the central halo mass tracks from state-of-the-art statistical dark matter accretion histories coupled to abundance matching routines. We show that too flat high-mass slopes in the input stellar-mass-halo-mass relations as predicted by previous works, imply non-physical stellar mass growth histories weaker than those implied by satellite accretion alone. Our best-fit models reproduce the satellite distributions at the largest masses and highest redshifts probed, the latest data on star formation rates and its bi-modality in the local Universe, and the correct fraction of ellipticals. Our results are important to predict robust and self-consistent stellar-mass-halo-mass relations and to generate reliable galaxy mock catalogues for the next generations of extra-galactic surveys such as Euclid and LSST.

We present constraints on the formation and evolution of early-type galaxies (ETGs) with the empirical model EMERGE. The parameters of this model are adjusted so that it reproduces the evolution of stellar mass functions, specific star formation rates, and cosmic star formation rates since $z\approx10$ as well as 'quenched' galaxy fractions and correlation functions. We find that at fixed halo mass present-day ETGs are more massive than late-type galaxies, whereas at fixed stellar mass ETGs populate more massive halos in agreement with lensing results. This effect naturally results from the shape and scatter of the stellar-to-halo mass relation and the galaxy formation histories. The ETG stellar mass assembly is dominated by 'in-situ' star formation below a stellar mass of $3\times10^{11}\mathrm{M}_\odot$ and by merging and accretion of 'ex-situ' formed stars at higher mass. The mass dependence is in tension with current cosmological simulations. Lower mass ETGs show extended star formation towards low redshift in agreement with recent estimates from IFU surveys. All ETGs have main progenitors on the 'main sequence of star formation' with the 'red sequence' appearing at $z \approx 2$. Above this redshift, over 95 per cent of the ETG progenitors are star-forming. More than 90 per cent of $z \approx 2$ 'main sequence' galaxies with $m_* > 10^{10}\mathrm{M}_\odot$ evolve into present-day ETGs. Above redshift 6, more than 80 per cent of the observed stellar mass functions above $10^{9}\mathrm{M}_\odot$ can be accounted for by ETG progenitors with $m_* > 10^{10}\mathrm{M}_\odot$. This implies that current and future high redshift observations mainly probe the birth of present-day ETGs. The source code and documentation of EMERGE are available at github.com/bmoster/emerge.

We implement an optically thin approximation for the effects of the local radiation field from stars and hot gas on the gas heating and cooling in the N-body SPH code GASOLINE2. We resimulate three galaxies from the NIHAO project: one dwarf, one Milky Way-like and one massive spiral, and study what are the local radiation field effects on various galaxy properties. We also study the effects of varying the Ultra Violet Background (UVB) model, by running the same galaxies with two different UVBs. Galaxy properties at $z=0$ like stellar mass, stellar effective mass radius, HI mass, and radial extent of the HI disc, show significant changes between the models with and without the local radiation field, and smaller differences between the two UVB models. The intrinsic effect of the local radiation field through cosmic time is to increase the equilibrium temperature at the interface between the galaxies and their circumgalactic media (CGM), moving this boundary inwards, while leaving relatively unchanged the gas inflow rate. Consequently, the temperature of the inflow increases when considering the local radiation sources. This temperature increase is a function of total galaxy mass, with a median CGM temperature difference of one order of magnitude for the massive spiral. The local radiation field suppresses the stellar mass growth by $\sim$20 per cent by $z=0$ for all three galaxies, while the HI mass is roughly halfed. The differences in the gas phase diagrams, significantly impact the HI column densities, shifting their peaks in the distributions towards lower $N_{\rm HI}$.

We investigate the build-up of the galactic dynamo and subsequently the origin of a magnetic driven outflow. We use a setup of an isolated disc galaxy with a realistic circum-galactic medium (CGM). We find good agreement of the galactic dynamo with theoretical and observational predictions from the radial and toroidal components of the magnetic field as function of radius and disc scale height. We find several field reversals indicating dipole structure at early times and quadrupole structure at late times. Together with the magnetic pitch angle and the dynamo control parameters $R_{\alpha}$, $R_{\omega}$ and $D$ we present strong evidence for an $\alpha^2$-$\Omega$ dynamo. The formation of a bar in the centre leads to further amplification of the magnetic field via adiabatic compression which subsequently drives an outflow. Due to the Parker-Instability the magnetic field lines rise to the edge of the disc, break out and expand freely in the CGM driven by the magnetic pressure. Finally, we investigate the correlation between magnetic field and star formation rate. Globally, we find that the magnetic field is increasing as function of the star formation rate surface density with a slope between $0.3$ and $0.45$ in good agreement with predictions from theory and observations. Locally, we find that the magnetic field can decrease while star formation increases. We find that this effect is correlated with the diffusion of magnetic field from the spiral arms to the inter-arm regions which we explicitly include by solving the induction equation and accounting for non-linear terms.

Galaxy mergers are a fundamental part of galaxy evolution. To study the resulting mass distributions of different kinds of galaxy mergers, we present a simulation suite of 36 high-resolution isolated merger simulations, exploring a wide range of parameter space in terms of mass ratios (mu = 1:5, 1:10, 1:50, 1:100) and orbital parameters. We find that mini mergers deposit a higher fraction of their mass in the outer halo compared to minor mergers, while their contribution to the central mass distribution is highly dependent on the orbital impact parameter: for larger pericentric distances we find that the centre of the host galaxy is almost not contaminated by merger particles. We also find that the median of the resulting radial mass distribution for mini mergers differs significantly from the predictions of simple theoretical tidal-force models. Furthermore, we find that mini mergers can increase the size of the host disc significantly without changing the global shape of the galaxy, if the impact occurs in the disc plane, thus providing a possible explanation for extended low-surface brightness disks reported in observations. Finally, we find clear evidence that streams are a strong indication of nearly circular infall of a satellite (with large angular momentum), whereas the appearance of shells clearly points to (nearly) radial satellite infall.

We use 25 simulated galaxies from the NIHAO project to define and characterize a variety of kinematic stellar structures: thin and thick discs, large scale single discs, classical and pseudo bulges, spheroids, inner discs, and stellar haloes. These structures have masses, spins, shapes and rotational support in good agreement with theoretical expectations and observational data. Above a dark matter halo mass of $2.5\times10^{\rm~11}M_{\rm\odot}$, all galaxies have a classical bulge and 70\% have a thin and thick disc. The kinematic (thin) discs follow a power-law relation between angular momentum and stellar mass $J_{\rm *}=3.4M_{\rm *}^{\rm1.26\pm0.06}$, in very good agreement with the prediction based on the empirical stellar-to-halo mass relation in the same mass range, and show a strong correlation between maximum `observed' rotation velocity and dark matter halo circular velocity $v_{\rm c}=6.4v_{\rm max}^{0.64\pm0.04}$. Tracing back in time these structures' progenitors, we find all to lose a fraction $1-f_j$ of their maximum angular momentum. Thin discs are significantly better at retaining their high-redshift spins ($f_j\sim0.70$) than thick ones ($f_j\sim0.40$). Stellar haloes have their progenitor baryons assembled the latest ($z_{\rm~1/2}\sim1.1$) and over the longest timescales ($\tau\sim6.2$~Gyr), and have the smallest fraction of stars born in-situ ($f_{\rm in-situ}=0.35\pm0.14$). All other structures have $1.5\lesssim z_{\rm1/2}\lesssim3$, $\tau=4\pm2$~Gyr and $f_{\rm in-situ}\gtrsim0.9$.

Over the past decade, empirical constraints on the galaxy-dark matter halo connection have significantly advanced our understanding of galaxy evolution. Past techniques have focused on connections between halo properties and galaxy stellar mass and/or star formation rates. Empirical techniques in the next decade will link halo assembly histories with galaxies' circumgalactic media, supermassive black holes, morphologies, kinematics, sizes, colors, metallicities, and transient rates. Uncovering these links will resolve many critical uncertainties in galaxy formation and will enable much higher-fidelity mock catalogs essential for interpreting observations. Achieving these results will require broader and deeper spectroscopic coverage of galaxies and their circumgalactic media; survey teams will also need to meet several criteria (cross-comparisons, public access, and covariance matrices) to facilitate combining data across different surveys. Acting on these recommendations will continue enabling dramatic progress in both empirical modeling and galaxy evolution for the next decade.

We present simulations of isolated disc galaxies in a realistic environment performed with the Tree-SPMHD-Code Gadget-3. Our simulations include a spherical circum-galactic medium (CGM) surrounding the galactic disc, motivated by observations and the results of cosmological simulations. We present three galactic models with different halo masses between 10e10 Msol and 10e12 Msol, and for each we use two different approaches to seed the magnetic field, as well as a control simulation without a magnetic field. We find that the amplification of the magnetic field in the centre of the disc leads to a biconical magnetic outflow of gas that magnetizes the CGM. This biconical magnetic outflow reduces the star formation rate (SFR) of the galaxy by roughly 40 percent compared to the simulations without magnetic fields. As the key aspect of our simulations, we find that small scale turbulent motion of the gas in the disc leads to the amplification of the magnetic field up to tens of 10e-6 G, as long as the magnetic field strength is low. For stronger magnetic fields turbulent motion does not lead to significant amplification but is replaced by an alpha-omega dynamo. The occurance of a small scale turbulent dynamo becomes apparent through the magnetic power spectrum and analysis of the field lines' curvature. In accordance with recent observations we find an anti-correlation between the spiral structure in the gas density and in the magnetic field due to a diffusion term added to the induction equation.

We simulate recoiling black hole trajectories from $z=20$ to $z=0$ in dark matter halos, quantifying how parameter choices affect escape velocities. These choices include the strength of dynamical friction, the presence of stars and gas, the accelerating expansion of the universe (Hubble acceleration), host halo accretion and motion, and seed black hole mass. $\Lambda$CDM halo accretion increases escape velocities by up to 0.6 dex and significantly shortens return timescales compared to non-accreting cases. Other parameters change orbit damping rates but have subdominant effects on escape velocities; dynamical friction is weak at halo escape velocities, even for extreme parameter values. We present formulae for black hole escape velocities as a function of host halo mass and redshift. Finally, we discuss how these findings affect black hole mass assembly as well as minimum stellar and halo masses necessary to retain supermassive black holes.

We present the first results of applying Gaussian Mixture Models in the stellar kinematic space of normalized angular momentum and binding energy on NIHAO high resolution galaxies to separate the stars into multiple components. We exemplify this method using a simulated Milky Way analogue, whose stellar component hosts: thin and thick discs, classical and pseudo bulges, and a stellar halo. The properties of these stellar structures are in good agreement with observational expectations in terms of sizes, shapes and rotational support. Interestingly, the two kinematic discs show surface mass density profiles more centrally concentrated than exponentials, while the bulges and the stellar halo are purely exponential. We trace back in time the Lagrangian mass of each component separately to study their formation history. Between z~3 and the end of halo virialization, z~1.3, all components lose a fraction of their angular momentum. The classical bulge loses the most (~95%) and the thin disc the least (~60%). Both bulges formed their stars in-situ at high redshift, while the thin disc formed ~98% in-situ, but with a constant SFR~1.5M$_{\rm\odot}$/yr$^{\rm-1}$ over the last ~11 Gyr. Accreted stars (6% of total stellar mass) are mainly incorporated to the thick disc or the stellar halo, which formed ex-situ 8% and 45% of their respective masses. Our analysis pipeline is freely available at https://github.com/aobr/gsf.

We empirically constrain how galaxy size relates to halo virial radius using new measurements of the size- and stellar mass-dependent clustering of galaxies in the Sloan Digital Sky Survey. We find that small galaxies cluster much more strongly than large galaxies of the same stellar mass. The magnitude of this clustering difference increases on small scales, and decreases with increasing stellar mass. Using Halotools to forward model the observations, we test an empirical model in which present-day galaxy size is proportional to the size of the virial radius at the time the halo reached its maximum mass. This simple model reproduces the observed size-dependence of galaxy clustering in striking detail. The success of this model provides strong support for the conclusion that satellite galaxies have smaller sizes relative to central galaxies of the same halo mass. Our findings indicate that satellite size is set prior to the time of infall, and that a remarkably simple, linear size--virial radius relation emerges from the complex physics regulating galaxy size. We make quantitative predictions for future measurements of galaxy-galaxy lensing, including dependence upon size, scale, and stellar mass, and provide a scaling relation of the ratio of mean sizes of satellites and central galaxies as a function of their halo mass that can be used to calibrate hydrodynamical simulations and semi-analytic models.

We present EMERGE, an Empirical ModEl for the foRmation of GalaxiEs, describing the evolution of individual galaxies in large volumes from $z\sim10$ to the present day. We assign a star formation rate to each dark matter halo based on its growth rate, which specifies how much baryonic material becomes available, and the instantaneous baryon conversion efficiency, which determines how efficiently this material is converted to stars, thereby capturing the baryonic physics. Satellites are quenched following the delayed-then-rapid model, and they are tidally disrupted once their subhalo has lost a significant fraction of its mass. The model is constrained with observed data extending out to high redshift. The empirical relations are very flexible, and the model complexity is increased only if required by the data, assessed by several model selection statistics. We find that for the same final halo mass galaxies can have very different star formation histories. Nevertheless, the average star formation and accretion rates are in good agreement with models following an abundance matching strategy. Galaxies that are quenched at $z=0$ typically have a higher peak star formation rate compared to their star-forming counterparts. The accretion of stars can dominate the total mass of massive galaxies, but is insignificant for low-mass systems, independent of star-formation activity. EMERGE predicts stellar-to-halo mass ratios for individual galaxies and introduces scatter self-consistently. We find that at fixed halo mass, passive galaxies have a higher stellar mass on average. The intra-cluster-mass in massive haloes can be up to 8 times larger than the mass of the central galaxy. Clustering for star-forming and quenched galaxies is in good agreement with observational constraints, indicating a realistic assignment of galaxies to haloes.

Using the Illustris simulation, we follow thousands of elliptical galaxies back in time to identify how the dichotomy between fast and slow rotating ellipticals (FRs and SRs) develops. Comparing to the $\textrm{ATLAS}^\textrm{3D}$ survey, we show that Illustris reproduces similar elliptical galaxy rotation properties, quantified by the degree of ordered rotation, $\lambda_\textrm{R}$. There is a clear segregation between low-mass ($M_{\rm *} < 10^{11} M_{\rm \odot}$) ellipticals, which form a smooth distribution of FRs, and high-mass galaxies ($M_{\rm *} > 10^{11.5} M_{\rm \odot}$), which are mostly SRs, in agreement with observations. We find that SRs are very gas poor, metal rich and red in colour, while FRs are generally more gas rich and still star forming. We suggest that ellipticals begin naturally as FRs and, as they grow in mass, lose their spin and become SRs. While at $z = 1$, the progenitors of SRs and FRs are nearly indistinguishable, their merger and star formation histories differ thereafter. We find that major mergers tend to disrupt galaxy spin, though in rare cases can lead to a spin-up. No major difference is found between the effects of gas-rich and gas-poor mergers and the amount of minor mergers seem to have little correlation with galaxy spin. In between major mergers, lower-mass ellipticals, which are mostly gas-rich, tend to recover their spin by accreting gas and stars. For galaxies with $M_{\rm *}$ above $\sim 10^{11} M_{\rm \odot}$, this trend reverses; galaxies only retain or steadily lose their spin. More frequent mergers, accompanied by an inability to regain spin, lead massive ellipticals to lose most of ordered rotation and transition from FRs to SRs.

We study the effect of mergers on the morphology of galaxies by means of the simulated merger tree approach first proposed by Moster et al. This method combines N-body cosmological simulations and semi-analytic techniques to extract realistic initial conditions for galaxy mergers. These are then evolved using high resolution hydrodynamical simulations, which include dark matter, stars, cold gas in the disc and hot gas in the halo. We show that the satellite mass accretion is not as effective as previously thought, as there is substantial stellar stripping before the final merger. The fraction of stellar disc mass transferred to the bulge is quite low, even in the case of a major merger, mainly due to the dispersion of part of the stellar disc mass into the halo. We confirm the findings of Hopkins et al., that a gas rich disc is able to survive major mergers more efficiently. The enhanced star formation associated with the merger is not localised to the bulge of galaxy, but a substantial fraction takes place in the disc too. The inclusion of the hot gas reservoir in the galaxy model contributes to reducing the efficiency of bulge formation. Overall, our findings suggest that mergers are not as efficient as previously thought in transforming discs into bulges. This possibly alleviates some of the tensions between observations of bulgeless galaxies and the hierarchical scenario for structure formation.

We use semi-analytic models and cosmological merger trees to provide the initial conditions for multi-merger numerical hydrodynamic simulations, and exploit these simulations to explore the effect of galaxy interaction and merging on star formation (SF). We compute numerical realisations of twelve merger trees from z=1.5 to z=0. We include the effects of the large hot gaseous halo around all galaxies, following recent obervations and predictions of galaxy formation models. We find that including the hot gaseous halo has a number of important effects. Firstly, as expected, the star formation rate on long timescales is increased due to cooling of the hot halo and refuelling of the cold gas reservoir. Secondly, we find that interactions do not always increase the SF in the long term. This is partially due to the orbiting galaxies transferring gravitational energy to the hot gaseous haloes and raising their temperature. Finally we find that the relative size of the starburst, when including the hot halo, is much smaller than previous studies showed. Our simulations also show that the order and timing of interactions are important for the evolution of a galaxy. When multiple galaxies interact at the same time, the SF enhancement is less than when galaxies interact in series. All these effects show the importance of including hot gas and cosmologically motivated merger trees in galaxy evolution models.

We employ cosmological hydrodynamical simulations to investigate the effects of AGN feedback on the formation of massive galaxies with present-day stellar masses of $M_{stel} = 8.8 \times 10^{10} - 6.0 \times 10^{11} M_{sun}$. Using smoothed particle hydrodynamics simulations with a pressure-entropy formulation that allows an improved treatment of contact discontinuities and fluid mixing, we run three sets of simulations of 20 halos with different AGN feedback models: (1) no feedback, (2) thermal feedback, and (3) mechanical and radiation feedback. We assume that seed black holes are present at early cosmic epochs at the centre of emerging dark matter halos and trace their mass growth via gas accretion and mergers with other black holes. Both feedback models successfully recover the observed M_BH - sigma relation and black hole-to-stellar mass ratio for simulated central early-type galaxies. The baryonic conversion efficiencies are reduced by a factor of two compared to models without any AGN feedback at all halo masses. However, massive galaxies simulated with thermal AGN feedback show a factor of ~10-100 higher X-ray luminosities than observed. The mechanical/radiation feedback model reproduces the observed correlation between X-ray luminosities and velocity dispersion, e.g. for galaxies with sigma = 200 km/s, the X-ray luminosity is reduced from $10^{42}$ erg/s to $10^{40}$ erg/s. It also efficiently suppresses late time star formation, reducing the specific star formation rate from $10^{-10.5}$ $yr^{-1}$ to $10^{-14}$ $yr^{-1}$ on average and resulting in quiescent galaxies since z=2, whereas the thermal feedback model shows higher late time in-situ star formation rates than observed.

We present the formation of a Kinematically Decoupled Core (KDC) in an elliptical galaxy, resulting from a major merger simulation of two disk galaxies. We show that although the two progenitor galaxies are initially following a prograde orbit, strong reactive forces during the merger can cause a short-lived change of their orbital spin; the two progenitors follow a retrograde orbit right before their final coalescence. This results in a central kinematic decoupling and the formation of a large-scale (~2 kpc radius) counter-rotating core (CRC) at the center of the final elliptical-like merger remnant (M*=1.3x10^11 Msun), while its outer parts keep the rotation direction of the initial orbital spin. The stellar velocity dispersion distribution of the merger remnant galaxy exhibits two symmetrical off-centered peaks, comparable to the observed "2-sigma galaxies". The KDC/CRC consists mainly of old, pre-merger population stars (older than 5 Gyr), remaining prominent in the center of the galaxy for more than 2 Gyr after the coalescence of its progenitors. Its properties are consistent with KDCs observed in massive elliptical galaxies. This new channel for the formation of KDCs from prograde mergers is in addition to previously known formation scenarios from retrograde mergers and can help towards explaining the substantial fraction of KDCs observed in early-type galaxies.

We present the smoothed-particle hydrodynamics implementation SPHGal, which combines some recently proposed improvements in GADGET. This includes a pressure-entropy formulation with a Wendland kernel, a higher order estimate of velocity gradients, a modified artificial viscosity switch with a modified strong limiter, and artificial conduction of thermal energy. With a series of idealized hydrodynamic tests we show that the pressure-entropy formulation is ideal for resolving fluid mixing at contact discontinuities but performs conspicuously worse at strong shocks due to the large entropy discontinuities. Including artificial conduction at shocks greatly improves the results. In simulations of Milky Way like disk galaxies a feedback-induced instability develops if too much artificial viscosity is introduced. Our modified artificial viscosity scheme prevents this instability and shows efficient shock capturing capability. We also investigate the star formation rate and the galactic outflow. The star formation rates vary slightly for different SPH schemes while the mass loading is sensitive to the SPH scheme and significantly reduced in our favored implementation. We compare the accretion behavior of the hot halo gas. The formation of cold blobs, an artifact of simple SPH implementations, can be eliminated efficiently with proper fluid mixing, either by conduction and/or by using a pressure-entropy formulation.

We show that recently documented trends in galaxy sizes with mass and redshift can be understood in terms of the influence of underlying cosmic evolution; a holistic view which is complimentary to interpretations involving the accumulation of discreet evolutionary processes acting on individual objects. Using standard cosmology theory, supported with results from the Millennium simulations, we derive expected size trends for collapsed cosmic structures, emphasising the important distinction between these trends and the assembly paths of individual regions. We then argue that the observed variation in the stellar mass content of these structures can be understood to first order in terms of natural limitations of cooling and feedback. But whilst these relative masses vary by orders of magnitude, galaxy and host radii have been found to correlate linearly. We explain how these two aspects will lead to galaxy sizes that closely follow observed trends and their evolution, comparing directly with the COSMOS and SDSS surveys. Thus we conclude that the observed minimum radius for galaxies, the evolving trend in size as a function of mass for intermediate systems, and the observed increase in the sizes of massive galaxies, may all be considered an emergent consequence of the cosmic expansion.

We study the effect of AGN mechanical and radiation feedback on the formation of bulge dominated galaxies via mergers of disc galaxies. The merging galaxies have mass-ratios of 1:1 to 6:1 and include pre-existing hot gaseous halos to properly account for the global impact of AGN feedback. Using smoothed particle hydrodynamics simulation code (GADGET-3) we compare three models with different AGN feedback models: (1) no black hole and no AGN feedback; (2) thermal AGN feedback; and (3) mechanical and radiative AGN feedback. The last model is motivated by observations of broad line quasars which show winds with initial velocities of $v_w \ge$ 10,000 km/s and also heating associated with the central AGN X-ray radiation. The primary changes in gas properties due to mechanical AGN feedback are lower thermal X-ray luminosity from the final galaxy - in better agreement with observations - and galactic outflows with higher velocity $\sim 1000$ km/s similar to recent direct observations of nearby merger remnants. The kinetic energy of the outflowing gas is a factor of $\sim$ 20 higher than in the thermal feedback case. All merger remnants with momentum-based AGN feedback with $v_w \sim 10,000$ km/s and $\epsilon_w=2 \times 10^{-3}$, independent of their progenitor mass-ratios, reproduce the observed relations between stellar velocity dispersion and black hole mass ($M_{\rm bh} - \sigma$) as well as X-ray luminosity ($L_X - \sigma$) with $10^{37.5}$ erg/s $\lesssim L_X (0.3-8~{\rm keV}) \lesssim 10^{39.2}$ erg/s for velocity dispersions in the range of 120 km/s $\lesssim \sigma \lesssim$ 190 km/s. In addition, the mechanical feedback produces a much greater AGN variability. We also show that gas is more rapidly and impulsively stripped from the galactic centres driving a moderate increase in galaxy size and decrease in central density with the mechanical AGN feedback model.

We use the Millennium Simulation, and an abundance-matching framework, to investigate the dynamical behaviour of galaxy pairs embedded in a cosmological context. Our main galaxy-pair sample, selected to have separations under 250 kpc/h, consists of over 1.3 million pairs at redshift z = 0, with stellar masses greater than 10^9 Msun, probing mass ratios down to 1:1000. We use dark matter halo membership and energy to classify our galaxy pairs. In terms of halo membership, central-satellite pairs tend to be in isolation (in relation to external more massive galaxies), are energetically- bound to each other, and are also weakly-bound to a neighbouring massive galaxy. Satellite-satellite pairs, instead, inhabit regions in close proximity to a more massive galaxy, are energetically-unbound, and are often bound to that neighbour. We find that 60% of our paired galaxies are bound to both their companion and to a third external object. Moreover, only 9% of our pairs resemble the kind of systems described by idealised binary merger simulations in complete isolation. In sum, we demonstrate the importance of properly connecting galaxy pairs to the rest of the Universe.

We present a new approach to study galaxy evolution in a cosmological context. We combine cosmological merger trees and semi-analytic models of galaxy formation to provide the initial conditions for multi-merger hydrodynamic simulations. In this way we exploit the advantages of merger simulations (high resolution and inclusion of the gas physics) and semi-analytic models (cosmological background and low computational cost), and integrate them to create a novel tool. This approach allows us to study the evolution of various galaxy properties, including the treatment of the hot gaseous halo from which gas cools and accretes onto the central disc, which has been neglected in many previous studies. This method shows several advantages over other methods. As only the particles in the regions of interest are included, the run time is much shorter than in traditional cosmological simulations, leading to greater computational efficiency. Using cosmological simulations, we show that multiple mergers are expected to be more common than sequences of isolated mergers, and therefore studies of galaxy mergers should take this into account. In this pilot study, we present our method and illustrate the results of simulating ten Milky Way-like galaxies since z=1. We find good agreement with observations for the total stellar masses, star formation rates, cold gas fractions and disc scale length parameters. We expect that this novel numerical approach will be very useful for pursuing a number of questions pertaining to the transformation of galaxy internal structure through cosmic time.

We examine the present-day total stellar-to-halo mass (SHM) ratio as a function of halo mass for a new sample of simulated field galaxies using fully cosmological, LCDM, high resolution SPH + N-Body simulations.These simulations include an explicit treatment of metal line cooling, dust and self-shielding, H2 based star formation and supernova driven gas outflows. The 18 simulated halos have masses ranging from a few times 10^8 to nearly 10^12 solar masses. At z=0 our simulated galaxies have a baryon content and morphology typical of field galaxies. Over a stellar mass range of 2.2 x 10^3 to 4.5 x 10^10 solar masses, we find extremely good agreement between the SHM ratio in simulations and the present-day predictions from the statistical Abundance Matching Technique presented in Moster et al. (2012). This improvement over past simulations is due to a number systematic factors, each decreasing the SHM ratios: 1) gas outflows that reduce the overall SF efficiency but allow for the formation of a cold gas component 2) estimating the stellar masses of simulated galaxies using artificial observations and photometric techniques similar to those used in observations and 3) accounting for a systematic, up to 30 percent overestimate in total halo masses in DM-only simulations, due to the neglect of baryon loss over cosmic times. Our analysis suggests that stellar mass estimates based on photometric magnitudes can underestimate the contribution of old stellar populations to the total stellar mass, leading to stellar mass errors of up to 50 percent for individual galaxies. These results highlight the importance of using proper techniques to compare simulations with observations and reduce the perceived tension between the star formation efficiency in galaxy formation models and in real galaxies.

LHS1070 is a nearby multiple system of low mass stars. It is an important source of information for probing the low mass end of the main sequence, down to the hydrogen-burning limit. The primary of the system is a mid-M dwarf and two components are late-M to early L dwarfs, at the star-brown dwarf transition. Hence LHS1070 is a valuable object to understand the onset of dust formation in cool stellar atmospheres.This work aims at determining the fundamental stellar parameters of LHS1070 and to test recent model atmospheres: BT-Dusty,BT-Settl, DRIFT, and MARCS models.Unlike in previous studies, we have performed a chi^2-minimization comparing well calibrated optical and infrared spectra with recent cool star synthetic spectra leading to the determination of the physical stellar parameters Teff, radius, and log g for each of the three components of LHS1070.

We present a new statistical method to determine the relationship between the stellar masses of galaxies and the masses of their host dark matter haloes over the entire cosmic history from z~4 to the present. This multi-epoch abundance matching (MEAM) model self-consistently takes into account that satellite galaxies first become satellites at times earlier than they are observed. We employ a redshift-dependent parameterization of the stellar-to-halo mass relation to populate haloes and subhaloes in the Millennium simulations with galaxies, requiring that the observed stellar mass functions at several redshifts be reproduced simultaneously. Using merger trees extracted from the dark matter simulations in combination with MEAM, we predict the average assembly histories of galaxies, separating into star formation within the galaxies (in-situ) and accretion of stars (ex-situ). The peak star formation efficiency decreases with redshift from 23% at z=0 to 9% at z=4 while the corresponding halo mass increases from 10^11.8M\odot to 10^12.5M\odot. The star formation rate of central galaxies peaks at a redshift which depends on halo mass; for massive haloes this peak is at early cosmic times while for low-mass galaxies the peak has not been reached yet. In haloes similar to that of the Milky-Way about half of the central stellar mass is assembled after z=0.7. In low-mass haloes, the accretion of satellites contributes little to the assembly of their central galaxies, while in massive haloes more than half of the central stellar mass is formed ex-situ with significant accretion of satellites at z<2. We find that our method implies a cosmic star formation history and an evolution of specific star formation rates which are consistent with those inferred directly. We present convenient fitting functions for stellar masses, star formation rates, and accretion rates as functions of halo mass and redshift.

We investigate the idea that the interaction of Dark Matter (DM) sub-halos with the gaseous disk of galaxies can be the origin for the observed holes and shells found in their neutral hydrogen (HI) distributions. We use high resolution hydrodynamic simulations to show that pure DM sub-halos impacting a galactic disk are not able to produce holes; on the contrary, they result in high density regions in the disk. However, sub-halos containing a small amount of gas (a few percent of the total DM mass of the sub-halo) are able to displace the gas in the disk and form holes and shells. The size and lifetime of these holes depend on the sub-halo gas mass, density and impact velocity. A DM sub-halo of mass $10^8$ $M_{\odot}$ and a gas mass fraction of $\sim 3$%, is able to create a kpc scale hole, with a lifetime similar to those observed in nearby galaxies. We also register an increase in the star formation rate at the rim of the hole, again in agreement with observations. Even though the properties off these simulated structures resemble those found in observations we find that the number of predicted holes (based on mass and orbital distributions of DM halos derived from cosmological N-body simulations) falls short compared to the observations. Only a handful of holes are produced per Giga year. This leads us to conclude that DM halo impact is not the major channel through which these holes are formed.

We employ hydrodynamical simulations to study the effects of dissipational gas physics on the vertical heating and thickening of disc galaxies during minor mergers. For the first time we present a suite of simulations that includes a diffuse, rotating, cooling, hot gaseous halo, as predicted by cosmological hydrodynamical simulations as well as models of galaxy formation. We study the effect of this new gaseous component on the vertical structure of a Milky Way-like stellar disc during 1:10 and 1:5 mergers. For 1:10 mergers we find no increased final thin disc scale height compared to the isolated simulation, leading to the conclusion that thin discs can be present even after a 1:10 merger if a reasonable amount of hot gas is present. The reason for this is the accretion of new cold gas, leading to the formation of a massive new thin stellar disc that dominates the surface brightness profile. In a previous study, in which we included only cold gas in the disk, we showed that the presence of cold gas decreased the thickening by a minor merger relative to the no-gas case. Here, we show that the evolution of the scale height in the presences of a cooling hot halo is dominated by the formation of the new stellar disc. In this scenario, the thick disc is the old stellar disc that has been thickened in a minor merger at z>1, while the thin disc is the new stellar disc that reforms after this merger. In addition, we study the evolution of the scale height during a 1:5 merger and find that a thin disc can be present even after this merger, provided enough hot gas is available. The final scale height in our simulations depends on the mass of the hot gaseous halo, the efficiency of the winds and the merger mass ratio. We find post-merger values in the range 0.5<z0<1.0 kpc in good agreement with observational constraints by local galaxies.

Cosmological hydrodynamical simulations as well as observations indicate that spiral galaxies are comprised of five different components: dark matter halo, stellar disc, stellar bulge, gaseous disc and gaseous halo. While the first four components have been extensively considered in numerical simulations of binary galaxy mergers, the effect of a hot gaseous halo has usually been neglected even though it can contain up to 80% of the total gas within the galaxy virial radius. We present a series of hydrodynamic simulations of major mergers of disc galaxies, that for the first time include a diffuse, rotating, hot gaseous halo. Through cooling and accretion, the hot halo can dissipate and refuel the cold gas disc before and after a merger. This cold gas can subsequently form stars, thus impacting the morphology and kinematics of the remnant. Simulations of isolated systems with total mass M~10^12Msun show a nearly constant star formation rate of ~5Msun/yr if the hot gaseous halo is included, while the star formation rate declines exponentially if it is neglected. We conduct a detailed study of the star formation efficiency during mergers and find that the presence of a hot gaseous halo reduces the starburst efficiency (e=0.5) compared to simulations without a hot halo (e=0.68). Moreover we find cases where the stellar mass of the merger remnant is lower than the sum of the stellar mass of the two progenitor galaxies when evolved in isolation. This suggests a revision to semi-analytic galaxy formation models which assume that a merger always leads to enhanced star formation. We show that adding the hot gas component has a significant effect on the kinematics and internal structure of the merger remnants, like an increased abundance of fast rotators and an r^(1/4) surface brightness profile at small scales.

Deep pencil beam surveys (<1 deg^2) are of fundamental importance for studying the high-redshift universe. However, inferences about galaxy population properties are in practice limited by 'cosmic variance'. This is the uncertainty in observational estimates of the number density of galaxies arising from the underlying large-scale density fluctuations. This source of uncertainty can be significant, especially for surveys which cover only small areas and for massive high-redshift galaxies. Cosmic variance for a given galaxy population can be determined using predictions from cold dark matter theory and the galaxy bias. In this paper we provide tools for experiment design and interpretation. For a given survey geometry we present the cosmic variance of dark matter as a function of mean redshift z and redshift bin size Dz. Using a halo occupation model to predict galaxy clustering, we derive the galaxy bias as a function of mean redshift for galaxy samples of a given stellar mass range. In the linear regime, the cosmic variance of these galaxy samples is the product of the galaxy bias and the dark matter cosmic variance. We present a simple recipe using a fitting function to compute cosmic variance as a function of the angular dimensions of the field, z, Dz and stellar mass m*. We also provide tabulated values and a software tool. We find that for GOODS at z=2 and with Dz=0.5 the relative cosmic variance of galaxies with m*>10^11 Msun is ~38%, while it is ~27% for GEMS and ~12% for COSMOS. For galaxies of m*~10^10 Msun the relative cosmic variance is ~19% for GOODS, ~13% for GEMS and ~6% for COSMOS. This implies that cosmic variance is a significant source of uncertainty at z=2 for small fields and massive galaxies, while for larger fields and intermediate mass galaxies cosmic variance is less serious.

We study the effect of dissipational gas physics on the vertical heating and thickening of disc galaxies during minor mergers. We produce a suite of minor merger simulations for Milky Way-like galaxies. This suite consists of collisionless simulations as well as hydrodynamical runs including a gaseous component in the galactic disc. We find that in dissipationless simulations minor mergers cause the scale height of the disc to increase by up to a factor of ~2. When the presence of gas in the disc is taken into account this thickening is reduced by 25% (50%) for an initial disc gas fraction of 20% (40%), leading to a final scale height z0 between 0.6 and 0.7 kpc, regardless of the initial scale height. We argue that the presence of gas reduces disc heating via two mechanisms: absorption of kinetic impact energy by the gas and/or formation of a new thin stellar disc that can cause heated stars to recontract towards the disc plane. We show that in our simulations most of the gas is consumed during the merger and thus the regrowth of a new thin disc has a negligible impact on the z0 of the post merger galaxy. Final disc scale heights found in our simulations are in good agreement with studies of the vertical structure of spiral galaxies where the majority of the systems are found to have scale heights of 0.4 kpc < z0 < 0.8 kpc. We also found no tension between recent measurements of the scale height of the Milky Way thin disc and results coming from our hydrodynamical simulations. We conclude that the existence of a thin disc in the Milky Way and in external galaxies is not in obvious conflict with the predictions of the CDM model.

We use a statistical approach to determine the relationship between the stellar masses of galaxies and the masses of the dark matter halos in which they reside. We obtain a parameterized stellar-to-halo mass (SHM) relation by populating halos and subhalos in an N-body simulation with galaxies and requiring that the observed stellar mass function be reproduced. We find good agreement with constraints from galaxy-galaxy lensing and predictions of semi-analytic models. Using this mapping, and the positions of the halos and subhalos obtained from the simulation, we find that our model predictions for the galaxy two-point correlation function (CF) as a function of stellar mass are in excellent agreement with the observed clustering properties in the SDSS at z=0. We show that the clustering data do not provide additional strong constraints on the SHM function and conclude that our model can therefore predict clustering as a function of stellar mass. We compute the conditional mass function, which yields the average number of galaxies with stellar masses in the range [m, m+dm] that reside in a halo of mass M. We study the redshift dependence of the SHM relation and show that, for low mass halos, the SHM ratio is lower at higher redshift. The derived SHM relation is used to predict the stellar mass dependent galaxy CF and bias at high redshift. Our model predicts that not only are massive galaxies more biased than low mass ones at all redshifts, but the bias increases more rapidly with increasing redshift for massive galaxies than for low mass ones. We present convenient fitting functions for the SHM relation as a function of redshift, the conditional mass function, and the bias as a function of stellar mass and redshift.

To contribute to the knowledge of dynamics of interplanetary dust we are searching for structures in the spatial distribution of interplanetary dust near the orbit of Venus. To this end we study the radial gradient of zodiacal light brightness, as observed by the zodiacal light photometer on board the Helios space probes on several orbits from 1975 to 1979. The cleanest data result from Helios B (= Helios 2) launched in January 1976. With respect to the general increase of zodiacal light brightness towards the Sun, the data show an excess brightness of a few percent for positions of the Helios space probe just outside the orbit of Venus. We consider this as evidence for a dust ring associated with the orbit of Venus, somewhat similar to that found earlier along the Earth's orbit.