We report high quality, H{\alpha} or CO rotation curves (RCs) to several Re for 41 large, massive, star-forming disk galaxies (SFGs), across the peak of cosmic galaxy evolution (z~0.67-2.45), taken with the ESO-VLT, the LBT and IRAM-NOEMA. Most RC41 SFGs have reflection symmetric RCs plausibly described by equilibrium dynamics. We fit the major axis position-velocity cuts with beam-convolved, forward modeling with a bulge, a turbulent rotating disk, and a dark matter (DM) halo. We include priors for stellar and molecular gas masses, optical light effective radii and inclinations, and DM masses from abundance matching scaling relations. Two-thirds or more of the z>1.2 SFGs are baryon dominated within a few Re of typically 5.5 kpc, and have DM fractions less than maximal disks (<fDM (Re)>=0.12). At lower redshift (z<1.2) that fraction is less than one-third. DM fractions correlate inversely with the baryonic angular momentum parameter, baryonic surface density and bulge mass. Inferred low DM fractions cannot apply to the entire disk & halo but more plausibly reflect a flattened, or cored, inner DM density distribution. The typical central 'DM deficit' in these cores relative to NFW distributions is ~30% of the bulge mass. The observations are consistent with rapid radial transport of baryons in the first generation massive gas rich halos forming globally gravitationally unstable disks, and leading to efficient build-up of massive bulges and central black holes. A combination of heating due to dynamical friction and AGN feedback may drive DM out of the initial cusps.

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.

We present an analysis of the origin and properties of the circum-galactic medium (CGM) in a suite of 11 cosmological zoom simulations resembling present day spiral galaxies. On average the galaxies retain about 50\% of the cosmic fraction in baryons, almost equally divided into disc (interstellar medium) gas, cool CGM gas and warm-hot CGM gas. At radii smaller than 50 kpc the CGM is dominated by recycled warm-hot gas injected from the central galaxy, while at larger radii it is dominated by cool gas accreted onto the halo. The recycled gas typically accounts for one-third of the CGM mass. We introduce the novel publicly available analysis tool \textsc{pygad} to compute ion abundances and mock absorption spectra. For Lyman-${\alpha}$ absorption we find good agreement of the simulated equivalent width (EW) distribution and observations out to large radii. Disc galaxies with quiescent assembly histories show significantly more absorption along the disc major axis. By comparing the EW and HI column densities we find that CGM Lyman-${\alpha}$ absorbers are best represented by an effective line-width $b\approx 50 - 70$ km s$^{-1}$ that increases mildly with halo mass, larger than typically assumed.

We use deep adaptive optics assisted integral field spectroscopy from SINFONI on the VLT to study the spatially resolved properties of ionized gas outflows driven by active galactic nuclei (AGN) in three galaxies at z~2.2 -- K20-ID5, COS4-11337 and J0901+1814. These systems probe AGN feedback from nuclear to circumgalactic scales, and provide unique insights into the different mechanisms by which AGN-driven outflows interact with their host galaxies. K20-ID5 and COS4-11337 are compact star forming galaxies with powerful $\sim$1500 km s$^{-1}$ AGN-driven outflows that dominate their nuclear H$\alpha$ emission. The outflows do not appear to have any impact on the instantaneous star formation activity of the host galaxies, but they carry a significant amount of kinetic energy which could heat the halo gas and potentially lead to a reduction in the rate of cold gas accretion onto the galaxies. The outflow from COS4-11337 is propagating directly towards its companion galaxy COS4-11363, at a projected separation of 5.4 kpc. COS4-11363 shows signs of shock excitation and recent truncation of star formation activity, which could plausibly have been induced by the outflow from COS4-11337. J0901+1814 is gravitationally lensed, giving us a unique view of a compact (R = 470 $\pm$ 70 pc), relatively low velocity ($\sim$650 km s$^{-1}$) AGN-driven outflow. J0901+1814 has a similar AGN luminosity to COS4-11337, suggesting that the difference in outflow properties is not related to the current AGN luminosity, and may instead reflect a difference in the evolutionary stage of the outflow and/or the coupling efficiency between the AGN ionizing radiation field and the gas in the nuclear regions.

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 study the evidence for a diversity of formation processes in early-type galaxies by presenting the first complete volume-limited sample of slow rotators with both integral-field kinematics from the ATLAS3D Project and high spatial resolution photometry from the Hubble Space Telescope. Analysing the nuclear surface brightness profiles of 12 newly imaged slow rotators, we classify their light profiles as core-less, and place an upper limit to the core size of ~10 pc. Considering the full magnitude and volume-limited ATLAS3D sample, we correlate the presence or lack of cores with stellar kinematics, including the proxy for the stellar angular momentum and the velocity dispersion within one half-light radius, stellar mass, stellar age, $\alpha$-element abundance, and age and metallicity gradients. More than half of the slow rotators have core-less light profiles, and they are all less massive than $10^{11}$ Msun. Core-less slow rotators show evidence for counter-rotating flattened structures, have steeper metallicity gradients, and a larger dispersion of gradient values than core slow rotators. Our results suggest that core and core-less slow rotators have different assembly processes, where the former are the relics of massive dissipation-less merging in the presence of central supermassive black holes. Formation processes of core-less slow rotators are consistent with accretion of counter-rotating gas or gas-rich mergers of special orbital configurations, which lower the final net angular momentum of stars, but support star formation. We also highlight core fast rotators as galaxies that share properties of core slow rotators and core-less slow rotators. Formation processes similar to those for core-less slow rotators can be invoked to explain the assembly of core fast rotators, with the distinction that these processes form or preserve cores.[Abridged]

We present the novel algorithmically regularised integration method MSTAR for high accuracy ($|\Delta E/E| \gtrsim 10^{-14}$) integrations of N-body systems using minimum spanning tree coordinates. The two-fold parallelisation of the $\mathcal{O}(N_\mathrm{part}^2)$ force loops and the substep divisions of the extrapolation method allows for a parallel scaling up to $N_\mathrm{CPU} = 0.2 \times N_\mathrm{part}$. The efficient parallel scaling of MSTAR makes the accurate integration of much larger particle numbers possible compared to the traditional algorithmic regularisation chain (AR-CHAIN) methods, e.g. $N_\mathrm{part} = 5000$ particles on $400$ CPUs for $1$ Gyr in a few weeks of wall-clock time. We present applications of MSTAR on few particle systems, studying the Kozai mechanism and N-body systems like star clusters with up to $N_\mathrm{part} =10^4$ particles. Combined with a tree or a fast multipole based integrator the high performance of MSTAR removes a major computational bottleneck in simulations with regularised subsystems. It will enable the next generation galactic-scale simulations with up to $10^9$ stellar particles (e.g. $m_\star = 100 M_\odot$ for a $M_\star = 10^{11} M_\odot$ galaxy) including accurate collisional dynamics in the vicinity of nuclear supermassive black holes.

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$.

We study the orbital evolution and gravitational wave (GW) emission of supermassive black hole (SMBH) binaries formed in gas-free mergers of massive early-type galaxies using the hybrid tree-regularized N-body code KETJU. The evolution of the SMBHs and the surrounding galaxies is followed self-consistently from the large-scale merger down to the final few orbits before the black holes coalesce. Post-Newtonian corrections are included up to PN3.5-level for the binary dynamics, and the GW calculations include the corresponding corrections up to PN1.0-level. We analyze the significance of the stellar environment on the evolution of the binary and the emitted GW signal during the final GW emission dominated phase of the binary hardening and inspiral. Our simulations are compared to semi-analytic models that have often been used for making predictions for the stochastic GW background emitted by SMBHs. We find that the commonly used semi-analytic parameter values produce large differences in merger timescales and eccentricity evolution, but result in only $\sim 10\%$ differences in the GW spectrum emitted by a single binary at frequencies $f\gtrsim 10^{-1} \, \rm yr^{-1}$, which are accessible by current pulsar timing arrays. These differences are in part caused by the strong effects of the SMBH binaries on the surrounding stellar population, which are not included in the semi-analytic models.

Cosmic ray (CR) protons are an important component in many astrophysical systems. Processes like CR injection, cooling, adiabatic changes as well as active CR transport through the medium strongly modify the CR momentum distribution and have to be taken into account in hydrodynamical simulations. We present an efficient novel numerical scheme to accurately compute the evolution of the particle distribution function by solving the Fokker-Planck equation with a low number of spectral bins (10 - 20), which is required to include a full spectrum for every computational fluid element. The distribution function is represented by piecewise power laws and is not forced to be continuous, which enables an optimal representation of the spectrum. The Fokker-Planck equation is solved with a two-moment approach evolving the CR number and energy density. The low numerical diffusion of the scheme reduces the numerical errors by orders of magnitude in comparison to classical schemes with piecewise constant spectral representations. With this method not only the spectral evolution of CRs can be computed accurately in magnetohydrodynamic simulations but also their dynamical impact as well as CR ionisation. This allows for more accurate models for astrophysical plasmas, like the interstellar medium, and direct comparisons with observations.

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 explore the effect of active galactic nucleus (AGN) feedback from central galaxies on their satellites by comparing two sets of cosmological zoom-in runs of 27 halos with masses ranging from $10^{12}$ to $10^{13.4}$ solar masses at z=0, with (wAGN) and without (noAGN) AGN feedback. Both simulations include stellar feedback from multiple processes, including powerful winds from supernovae, stellar winds from young massive stars, AGB stars, radiative heating within Str\"omgren spheres and photoelectric heating. Our wAGN model is identical to the noAGN model except that it also includes a model for black hole seeding and accretion, as well as AGN feedback via high-velocity broad absorption line winds and Compton/photoionization heating. We show that the inclusion of AGN feedback from the central galaxy significantly affects the star formation history and the gas content of the satellite galaxies. AGN feedback starts to affect the gas content and the star formation of the satellites as early as z=2. The mean gas rich fraction of satellites at z=0 decreases from 15% in the noAGN simulation to 5% in the wAGN simulation. The difference between the two sets extends as far out as five times the virial radius of the central galaxy at z=1. We investigate the quenching mechanism by studying the physical conditions in the surroundings of pairs of satellites matched across the wAGN and noAGN simulations and find an increase in the temperature and relative velocity of the intergalactic gas.

We present a hydro-dynamical simulation at sub-parsec and few solar mass resolution of a merger between two gas-rich dwarf galaxies. Our simulation includes a detailed model for the multi-phase interstellar medium (ISM) and is able to follow the entire formation history of spatially resolved star clusters, including feedback from individual massive stars. Shortly after the merger we find a population of $\sim 900$ stellar clusters with masses above $10^{2.5}\; \rm{M_\odot}$ and a cluster mass function (CMF), which is well fitted with a power-law with a slope of $\alpha=-1.70\pm0.08$. We describe here in detail the formation of the three most massive clusters (M$_{*} \gtrsim 10^5$ M$_\odot$), which populate the high-mass end of the CMF. The simulated clusters form rapidly on a timescale of $6$-$8$ Myr in converging flows of dense gas. The embedded merger phase has extremely high star formation rate surface densities of $\Sigma_\mathrm{SFR}>10\; \mathrm{M}_\odot\; \mathrm{yr}^{-1}\; \mathrm{kpc}^{-2}$ and thermal gas pressures in excess of $P_{\rm th}\sim10^7 \; \mathrm{k}_{\rm B}\;(\rm K\;\mathrm{cm}^{-3})^{-1}$. The formation process is terminated by rapid gas expulsion driven by the first generation of supernovae, after which the cluster centers relax and both their structure and kinematics become indistinguishable from observed local globular clusters. The simulation presented here provides a general model for the formation of metal-poor globular clusters in chemically unevolved starbursting environments of low-mass dwarf galaxies, which are common at high redshifts.

We present the $0.6<z<2.6$ evolution of the ionized gas velocity dispersion in 175 star-forming disk galaxies based on data from the full KMOS$^{\rm 3D}$ integral field spectroscopic survey. In a forward-modelling Bayesian framework including instrumental effects and beam-smearing, we fit simultaneously the observed galaxy velocity and velocity dispersion along the kinematic major axis to derive the intrinsic velocity dispersion $\sigma_0$. We find a reduction of the average intrinsic velocity dispersion of disk galaxies as a function of cosmic time, from $\sigma_0\sim45$ km s$^{-1}$ at $z\sim2.3$ to $\sigma_0\sim30$ km s$^{-1}$ at $z\sim0.9$. There is substantial intrinsic scatter ($\sigma_{\sigma_0, {\rm int}}\approx10$ km s$^{-1}$) around the best-fit $\sigma_0-z$-relation beyond what can be accounted for from the typical measurement uncertainties ($\delta\sigma_0\approx12$ km s$^{-1}$), independent of other identifiable galaxy parameters. This potentially suggests a dynamic mechanism such as minor mergers or variation in accretion being responsible for the scatter. Putting our data into the broader literature context, we find that ionized and atomic+molecular velocity dispersions evolve similarly with redshift, with the ionized gas dispersion being $\sim10-15$ km s$^{-1}$ higher on average. We investigate the physical driver of the on average elevated velocity dispersions at higher redshift, and find that our galaxies are at most marginally Toomre-stable, suggesting that their turbulent velocities are powered by gravitational instabilities, while stellar feedback as a driver alone is insufficient. This picture is supported through comparison with a state-of-the-art analytical model of galaxy evolution.

We study the destruction of interstellar dust via sputtering in supernova (SN) shocks using three-dimensional hydrodynamical simulations. With a novel numerical framework, we follow both sputtering and dust dynamics governed by direct collisions, plasma drag and betatron acceleration. Grain-grain collisions are not included and the grain-size distribution is assumed to be fixed. The amount of dust destroyed per SN is quantified for a broad range of ambient densities and fitting formulae are provided. Integrated over the grain-size distribution, nonthermal (inertial) sputtering dominates over thermal sputtering for typical ambient densities. We present the first simulations that explicitly follow dust sputtering within a turbulent multiphase interstellar medium. We find that the dust destruction timescales $\tau$ are 0.35 Gyr for silicate dust and 0.44 Gyr for carbon dust in solar neighborhood conditions. The SN environment has an important impact on $\tau$. SNe that occur in preexisting bubbles destroy less dust as the destruction is limited by the amount of dust in the shocked gas. This makes $\tau$ about 2.5 times longer than the estimate based on results from a single SN explosion. We investigate the evolution of the dust-to-gas mass ratio (DGR), and find that a spatial inhomogeneity of $\sim$ 14\% develops for scales below 10 pc. It locally correlates positively with gas density but negatively with gas temperature even in the exterior of the bubbles due to incomplete gas mixing. This leads to a $\sim$ 30\% lower DGR in the volume filling warm gas compared to that in the dense clouds.

We present a census of ionized gas outflows in 599 normal galaxies at redshift 0.6<z<2.7, mostly based on integral field spectroscopy of Ha, [NII], and [SII] line emission. The sample fairly homogeneously covers the main sequence of star-forming galaxies with masses 9.0<log(M*/Msun)<11.7, and probes into the regimes of quiescent galaxies and starburst outliers. About 1/3 exhibits the high-velocity component indicative of outflows, roughly equally split into winds driven by star formation (SF) and active galactic nuclei (AGN). The incidence of SF-driven winds correlates mainly with star formation properties. These outflows have typical velocities of ~450 km/s, local electron densities of n_e~380 cm^-3, modest mass loading factors of ~0.1-0.2 at all galaxy masses, and energetics compatible with momentum driving by young stellar populations. The SF-driven winds may escape from log(M*/Msun)<10.3 galaxies but substantial mass, momentum, and energy in hotter and colder outflow phases seem required to account for low galaxy formation efficiencies in the low-mass regime. Faster AGN-driven outflows (~1000-2000 km/s) are commonly detected above log(M*/Msun)~10.7, in up to ~75% of log(M*/Msun)>11.2 galaxies. The incidence, strength, and velocity of AGN-driven winds strongly correlates with stellar mass and central concentration. Their outflowing ionized gas appears denser (n_e~1000 cm^-3), and possibly compressed and shock-excited. These winds have comparable mass loading factors as the SF-driven winds but carry ~10 (~50) times more momentum (energy). The results confirm our previous findings of high duty cycle, energy-driven outflows powered by AGN above the Schechter mass, which may contribute to star formation quenching.

We investigate the relationship between star formation activity and outflow properties on kiloparsec scales in a sample of 28 star forming galaxies at $z\sim$ 2-2.6, using adaptive optics assisted integral field observations from SINFONI on the VLT. The narrow and broad components of the H$\alpha$ emission are used to simultaneously determine the local star formation rate surface density ($\Sigma_{\rm SFR}$), and the outflow velocity $v_{\rm out}$ and mass outflow rate $\dot{M}_{\rm out}$, respectively. We find clear evidence for faster outflows with larger mass loading factors at higher $\Sigma_{\rm SFR}$. The outflow velocities scale as $v_{\rm out}$ $\propto$ $\Sigma_{\rm SFR}^{0.34 \pm 0.10}$, which suggests that the outflows may be driven by a combination of mechanical energy released by supernova explosions and stellar winds, as well as radiation pressure acting on dust grains. The majority of the outflowing material does not have sufficient velocity to escape from the galaxy halos, but will likely be re-accreted and contribute to the chemical enrichment of the galaxies. In the highest $\Sigma_{\rm SFR}$ regions the outflow component contains an average of $\sim$45% of the H$\alpha$ flux, while in the lower $\Sigma_{\rm SFR}$ regions only $\sim$10% of the H$\alpha$ flux is associated with outflows. The mass loading factor, $\eta$ = $\dot{M}_{\rm out}$/SFR, is positively correlated with $\Sigma_{\rm SFR}$ but is relatively low even at the highest $\Sigma_{\rm SFR}$: $\eta \lesssim$ 0.5 $\times$ (380 cm$^{-3}$/n$_e$). This may be in tension with the $\eta$ $\gtrsim$ 1 required by cosmological simulations, unless a significant fraction of the outflowing mass is in other gas phases and has sufficient velocity to escape the galaxy halos.

We study the impact of merging supermassive black holes (SMBHs) on the central regions of massive early-type galaxies (ETGs) using a series of merger simulations with varying mass ratios. The ETG models include realistic stellar and dark matter components and are evolved with the GADGET based regularized tree code KETJU. We show that observed key properties of the nuclear stellar populations of massive ETGs, namely flat stellar density distributions (cores), tangentially biased velocity distributions and kinematically decoupled (counter-)rotation can naturally result from a single process $-$ the scouring by SMBHs. Major mergers with mass ratios of $q>1/3$ produce flat, tangentially biased cores with kinematically distinct components. These features originate from spin reversals during the decay of the SMBH orbits caused by gravitational torques from tidally expelled nuclear material. Minor mergers ($q\lesssim1/3$) on the other hand, form non-rotating cores and the tidal reversal becomes less important. Low-density stellar cores scoured in (multiple) minor mergers are less tangentially biased. This implies that the nuclear stellar properties of massive ETGs can be solely explained by stellar dynamical processes during their final assembly without any need for `feedback' from accreting black holes. We predict a strong correlation between decoupled cores, central anisotropy and merger history: decoupled cores form in binary mergers and we predict them to occur in elliptical galaxies with the strongest central anisotropy. Measurements of the central orbital structure are the key to understanding the number of mergers a given galaxy has experienced.

Following the success of the Plateau de Bure high-z Blue Sequence Survey (PHIBSS), we present the PHIBSS2 legacy program, a survey of the molecular gas properties of star-forming galaxies on and around the star formation main sequence (MS) at different redshifts using NOEMA. This survey significantly extends the existing sample of star-forming galaxies with CO molecular gas measurements, probing the peak epoch of star formation (z=1-1.6) as well as its building-up (z=2-3) and winding-down (z=0.5-0.8) phases. The targets are drawn from the GOODS, COSMOS, and AEGIS deep fields and uniformly sample the MS in the stellar mass (M*) - star formation rate (SFR) plane with log(M*/Msun) = 10-11.8. We describe the survey strategy and sample selection before focusing on the results obtained at z=0.5-0.8, where we report 60 CO(2-1) detections out of 61 targets. We determine their molecular gas masses and separately obtain disc sizes and bulge-to-total (B/T) luminosity ratios from HST I-band images. The median molecular gas-to-stellar mass ratio, gas fraction, and depletion time as well as their dependence with M* and offset from the MS follow published scaling relations for a much larger sample of galaxies spanning a wider range of redshifts. The galaxy-averaged Kennicutt-Schmidt (KS) relation between molecular gas and SFR surface densities is strikingly linear, pointing towards similar star formation timescales within galaxies at any given epoch. In terms of morphology, the molecular gas content, the SFR, the disc stellar mass, and the disc molecular gas fraction do not seem to correlate with B/T and the stellar surface density, which suggests an ongoing supply of fresh molecular gas to compensate for the build-up of the bulge. Our measurements do not yield any significant variation of the depletion time with B/T and hence no strong evidence for morphological quenching within the scatter of the MS.

Galactic outflows are critical to our understanding of galaxy formation and evolution. However the details of the underlying feedback process remain unclear. We compare Ly$\alpha$ observations of the circumgalactic medium (CGM) of Lyman Break Galaxies (LBGs) with mock observations of their simulated CGM. We use cosmological hydrodynamical `zoom-in' simulations of an LBG which contains strong, momentum-driven galactic outflows. Simulation snapshots at $z=2.2$ and $z=2.65$ are used, corresponding to the available observational data. The simulation is post-processed with the radiative transfer code \textsc{crash} to account for the impact of ionising photons on hydrogen gas surrounding the simulated LBG. We generate mock absorption line maps for comparison with data derived from observed close galaxy-galaxy pairs. We perform calculations of Ly$\alpha$ photons scattering through the CGM with our newly developed Monte-Carlo code \textsc{slaf}, and compare to observations of diffuse Ly$\alpha$ halos around LBGs. Our fiducial galactic outflow model comes closer to reproducing currently observed characteristics of the CGM in Ly$\alpha$ than a reference inefficient feedback model used for comparison. Nevertheless, our fiducial model still struggles to reproduce the observed data of the inner CGM (at impact parameter $b<30$kpc). Our results suggest that galactic outflows affect Ly$\alpha$ absorption and emission around galaxies mostly at impact parameters $b<50$ kpc, while cold accretion flows dominate at larger distances. We discuss the implications of this result, and underline the potential constraining power of CGM observations - in emission and absorption - on galactic outflow models.

As part of the SILCC-ZOOM project we present our first sub-parsec resolution radiation-hydrodynamic simulations of two molecular clouds self-consistently forming from a turbulent, multi-phase ISM. The clouds have similar initial masses of few 10$^4$ M$_{\odot}$, escape velocities of ~5 km s$^{-1}$, and a similar initial energy budget. We follow the formation of star clusters with a sink based model and the impact of radiation from individual massive stars with the tree-based radiation transfer module TreeRay. Photo-ionizing radiation is coupled to a chemical network to follow gas heating, cooling and molecule formation and dissociation. For the first 3 Myr of cloud evolution we find that the overall star formation effciency is considerably reduced by a factor of ~4 to global cloud values of < 10 % as the mass accretion of sinks that host massive stars is terminated after <1 Myr. Despite the low effciency, star formation is triggered across the clouds. Therefore, a much larger region of the cloud is affected by radiation and the clouds begin to disperse. The time scale on which the clouds are dispersed sensitively depends on the cloud substructure and in particular on the amount of gas at high visual extinction. The damage of radiation done to the highly shielded cloud (MC1) is delayed. We also show that the radiation input can sustain the thermal and kinetic energy of the clouds at a constant level. Our results strongly support the importance of ionizing radiation from massive stars for explaining the low observed star formation effciency of molecular clouds

We investigate the connection between the HI~content, SFR and environment of galaxies using a hydrodynamic simulation that incorporates scaling relations for galactic wind and a heuristic halo mass-based quenching prescription. We run two zoom-in simulations of galaxy groups with $M_{halo}>10^{13}M_\odot$ at z=0, selected to have quiet merger histories. We track galaxies as they become satellites, and compute the delay time $\tau_{d}$ during which the satellites are similar to central galaxies at a given stellar mass, and a fading time $\tau_{f}$ during which satellites go from gas-rich and star-forming to gas-poor and quiescent. We find $0.7< \tau_{d}< 3$ Gyr at $z=0$, and depends inversely on the satellite halo mass at infall. At z ~ 1 we find ~$0.3< \tau_{d} < 2$ Gyr, broadly consistent with a positive correlation with the Hubble time. For a given halo mass, lower stellar mass galaxies at infall time have higher $\tau_{d}$. We generally find $\tau_{f}\ll \tau_{d}$, ranging between ~ 150 Myr at z~0 and ~ 80 Myr at z ~ 1 based on linear interpolation, with some uncertainty because they are smaller than our simulation output frequency ($200-300$ Myr). $\tau_{f}$ has no obvious dependency on infall halo mass. Both timescales show little difference between HI depletion and SF quenching, indicating that using up the gas reservoir by star formation without refilling is the main mechanism to transform satellite galaxies at these halo masses. At a given physical distance from the center of the main halo of interest, higher redshift galaxies have on average higher cold gas content, but the ratio of gas (HI or H$_2$) to star formation rate is similar, indicating that star formation is consistently fed through reservoirs of HI then H$_2$. For a given amount of HI, galaxies have shorter consumption times in more massive halo structures.

We present a series of 20 cosmological zoom simulations of the formation of massive galaxies with and without a model for AGN feedback. Differences in stellar population and kinematic properties are evaluated by constructing mock integral field unit (IFU) maps. The impact of the AGN is weak at high redshift when all systems are mostly fast-rotating and disc-like. After $z \sim 1$ the AGN simulations result in lower mass, older, less metal rich and slower rotating systems with less disky isophotes - in general agreement with observations. Two-dimensional kinematic maps of in-situ and accreted stars show that these differences result from reduced in-situ star formation due to AGN feedback. A full analysis of stellar orbits indicates that galaxies simulated with AGN are typically more triaxial and have higher fractions of x-tubes and box orbits and lower fractions of z-tubes. This trend can also be explained by reduced late in-situ star formation. We introduce a global parameter, $\xi_3$ , to characterise the anti-correlation between the third-order kinematic moment $h_3$ and the line-of-sight velocity ($v_{los}/{\sigma}$), and compare to ATLAS$^{3D}$ observations. The kinematic asymmetry parameter $\xi_3$ might be a useful diagnostic for large integral field surveys as it is a kinematic indicator for intrinsic shape and orbital content.

We compute synthetic optical and ultraviolet (UV) emission-line properties of galaxies in a full cosmological framework by coupling, in post-processing, new-generation nebular-emission models with high-resolution, cosmological zoom-in simulations of massive galaxies. Our self-consistent modelling accounts for nebular emission from young stars and accreting black holes (BHs). We investigate which optical- and UV-line diagnostic diagrams can best help to discern between the main ionizing sources, as traced by the ratio of BH accretion to star formation rates in model galaxies, over a wide range of redshifts. At low redshift, simulated star-forming galaxies, galaxies dominated by active galactic nuclei and composite galaxies are appropriately differentiated by standard selection criteria in the classical [OIII]$\lambda$5007/H$\beta$ versus [NII]$\lambda$6584/H$\alpha$ diagram. At redshifts $z \gt 1$, however, this optical diagram fails to discriminate between active and inactive galaxies at metallicities below $0.5\ Z_\odot$. To robustly classify the ionizing radiation of such metal-poor galaxies, which dominate in the early Universe, we confirm 3 previous, and propose 11 novel diagnostic diagrams based on equivalent widths and luminosity ratios of UV emission lines, such as EW(OIII]$\lambda$1663) versus OIII]$\lambda$1663/HeII$\lambda$1640, CIII]$\lambda$1908/HeII$\lambda$1640 versus OIII]$\lambda$1663/HeII$\lambda$1640, and CIV$\lambda$1550/CIII]$\lambda$1908 versus CIII]$\lambda$1908/CII$\lambda$2326. We formulate associated UV selection criteria and discuss some caveats of our results (e.g., uncertainties in the modelling of the HeII$\lambda$1640 line). These UV diagnostic diagrams are potentially important for the interpretation of high-quality spectra of very distant galaxies to be gathered by next-generation telescopes, such as the James Webb Space Telescope.

The C$^{+}$ ion is an important coolant of interstellar gas, and so the [CII] fine structure line is frequently observed in the interstellar medium. However, the physical and chemical properties of the [CII]-emitting gas are still unclear. We carry out non-LTE radiative transfer simulations with RADMC-3D to study the [CII] line emission from a young, turbulent molecular cloud before the onset of star formation, using data from the SILCC-Zoom project. The [CII] emission is optically thick over 40% of the observable area with $I_{[\textrm{CII}]} > 0.5$ K km s$^{-1}$. To determine the physical properties of the [CII] emitting gas, we treat the [CII] emission as optically thin. We find that the [CII] emission originates primarily from cold, moderate density gas ($40 \lesssim T \lesssim 65$ K and $50 \lesssim n \lesssim 440$ cm$^{-3}$), composed mainly of atomic hydrogen and with an effective visual extinction between $\sim 0.50$ and $\sim 0.91$. Gas dominated by molecular hydrogen contributes only $\lesssim$20% of the total [CII] line emission. Thus, [CII] is not a good tracer for CO-dark H$_2$ at this early phase in the cloud's lifetime. We also find that the total gas, H and C$^+$ column densities are all correlated with the integrated [CII] line emission, with power law slopes ranging from 0.5 to 0.7. Further, the median ratio between the total column density and the [CII] line emission is $Y_{{\rm CII}}\approx 1.1 \times 10^{21}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$, and $Y_{{\rm CII}}$ scales with $I_{[\textrm{CII}]}^{-0.3}$. We expect $Y_{{\rm CII}}$ to change in environments with a lower or higher radiation field than simulated here.

We use cosmological hydrodynamical simulations to investigate the role of feedback from accreting black holes on the evolution of sizes, compactness, stellar core density and specific star-formation of massive galaxies with stellar masses of $ M_{star} > 10^{10.9} M_{\odot}$. We perform two sets of cosmological zoom-in simulations of 30 halos to z=0: (1) without black holes and Active Galactic Nucleus (AGN) feedback and (2) with AGN feedback arising from winds and X-ray radiation. We find that AGN feedback can alter the stellar density distribution, reduce the core density within the central 1 kpc by 0.3 dex from z=1, and enhance the size growth of massive galaxies. We also find that galaxies simulated with AGN feedback evolve along similar tracks to those characterized by observations in specific star formation versus compactness. We confirm that AGN feedback plays an important role in transforming galaxies from blue compact galaxies into red extended galaxies in two ways: (1) it effectively quenches the star formation, transforming blue compact galaxies into compact quiescent galaxies and (2) it also removes and prevents new accretion of cold gas, shutting down in-situ star formation and causing subsequent mergers to be gas-poor or mixed. Gas poor minor mergers then build up an extended stellar envelope. AGN feedback also puffs up the central region through the fast AGN driven winds as well as the slow expulsion of gas while the black hole is quiescent. Without AGN feedback, large amounts of gas accumulate in the central region, triggering star formation and leading to overly massive blue galaxies with dense stellar cores.

Given its velocity dispersion, the early-type galaxy NGC 1600 has an unusually massive ($M_\bullet = 1.7 \times 10^{10} M_\odot$) central supermassive black hole (SMBH), surrounded by a large core ($r_\mathrm{b} = 0.7$ kpc) with a tangentially biased stellar distribution. We present high-resolution equal-mass merger simulations including SMBHs to study the formation of such systems. The structural parameters of the progenitor ellipticals were chosen to produce merger remnants resembling NGC 1600. We test initial stellar density slopes of $\rho \propto r^{-1}$ and $\rho \propto r^{-3/2}$ and vary the initial SMBH masses from $8.5 \times 10^8$ to $8.5 \times 10^9$ $M_\odot$. With increasing SMBH mass the merger remnants show a systematic decrease in central surface brightness, an increasing core size, and an increasingly tangentially biased central velocity anisotropy. Two-dimensional kinematic maps reveal decoupled, rotating core regions for the most massive SMBHs. The stellar cores form rapidly as the SMBHs become bound, while the velocity anisotropy develops more slowly after the SMBH binaries become hard. The simulated merger remnants follow distinct relations between the core radius and the sphere-of-influence, and the SMBH mass, similar to observed systems. We find a systematic change in the relations as a function of the progenitor density slope, and present a simple scouring model reproducing this behavior. Finally, we find the best agreement with NGC 1600 using SMBH masses totaling the observed value of $M_\bullet = 1.7 \times 10^{10} M_\odot$. In general, density slopes of $\rho \propto r^{-3/2}$ for the progenitor galaxies are strongly favored for the equal-mass merger scenario.

Magnetic fields are ubiquitously observed in the interstellar medium (ISM) of present-day star-forming galaxies with dynamically relevant energy densities. Using three-dimensional magneto-hydrodynamic (MHD) simulations of the supernova (SN) driven ISM in the flux-freezing approximation (ideal MHD) we investigate the impact of the magnetic field on the chemical and dynamical evolution of the gas, fragmentation and the formation of molecular clouds. We follow the chemistry with a network of six species (H$^{+}$, H, H$_2$, C$^+$, CO, free electrons) including local shielding effects. We find that magnetic fields thicken the disc by a factor of a few to a scale height of $\sim100\,\mathrm{pc}$, delay the formation of dense (and molecular) gas by $\sim25\,\mathrm{Myr}$ and result in differently shaped gas structures. The magnetised gas fragments into fewer clumps, which are initially at subcritical mass-to-flux ratios, $M/\Phi\approx0.3(M/\Phi)_\mathrm{crit}$, and accrete gas preferentially parallel to the magnetic field lines until supercritial mass-to-flux ratios of up to order 10 are reached. The accretion rates onto molecular clouds scale with $\dot{M}\propto M^{1.5}$. The median of the inter-cloud velocity dispersion is $\sim2-5\,\mathrm{km\,s}^{-1}$ and lower than the internal velocity dispersion in the clouds ($\sim3-7\,\mathrm{km\,s}^{-1}$). However, individual cloud-cloud collisions occur at speeds of a few $10\,\mathrm{km\,s}^{-1}$.

Nuclear outflows driven by accreting massive black holes are one of the main feedback mechanisms invoked at high-z to reproduce the distinct separation between star-forming, disk galaxies and quiescent spheroidal systems. Yet, our knowledge of feedback at high-z remains limited by the lack of observations of the multiple gas phases in galaxy outflows. In this work we use new deep, high-spatial resolution ALMA CO(3-2) and archival VLT/SINFONI H$\alpha$ observations to study the molecular and ionized components of the AGN-driven outflow in zC400528 ---a massive, main sequence galaxy at z=2.3 in the process of quenching. We detect a powerful molecular outflow that shows a positive velocity gradient and extends for at least ~10 kpc from the nuclear region, about three times the projected size of the ionized wind. The molecular gas in the outflow does not reach velocities high enough to escape the galaxy and is therefore expected to be reaccreted. Keeping in mind the various assumptions involved in the analysis, we find that the mass and energetics of the outflow are dominated by the molecular phase. The AGN-driven outflow in zC400528 is powerful enough to deplete the molecular gas reservoir on a timescale at least twice shorter than that needed to exhaust it by star formation. This suggests that the nuclear outflow is one of the main quenching engines at work in the observed suppression of the central star-formation activity in zC400528.

Recent magneto-hydrodynamic simulations of the star-forming interstellar medium (ISM) with parsec scale resolution indicate that relativistic cosmic rays support the launching of galactic outflows on scales of a few kpc. If these fundamental constituents of the ISM are injected at the sites of supernova (SN) explosions, the outflows are smoother, colder, and denser than the highly structured, hot-phase driven outflows forming, e.g., by thermal SN energy injection alone. In this Letter we present computations of resonant Lyman-$\alpha$ (Ly$\alpha$) radiation transfer through snapshots of a suite of stratified disk simulations from the SILCC project. For a range of thermal, radiative, and kinetic feedback models only simulations including non-thermal cosmic rays produce Ly$\alpha$ spectra with enhanced red peaks and strong absorption at line center -- similar to observed systems. The absence of cosmic ray feedback leads to spectra incompatible with observations. We attribute this to the smoother neutral gas distribution of cosmic ray supported outflows within a few kpc from the disk midplane.

We investigate the impact of cosmic rays (CRs) on galactic outflows from a multi-phase interstellar medium with solar neighbourhood conditions. The three-dimensional magneto-hydrodynamical simulations include CRs as a relativistic fluid in the advection-diffusion approximation. The thermal and chemical state of the ISM is computed with a non-equilibrium chemical network. We find that CRs (injected with 10 \% of the supernova energy) efficiently support the launching of outflows and strongly affect their phase structure. Outflows leaving the midplane are denser ($\rho \sim 10^{-26}\,\mathrm{g\,cm}^{-3}$), colder ($\sim 10^4\,\mathrm{K}$), and slower ($\sim 30\,\mathrm{km\,s}^{-1}$) if CRs are considered in addition to thermal SNe. The CR supported outflows are also smoother, in particular at larger heights ($> 1\,\mathrm{kpc}$ above the midplane) without the direct impact of SN explosions. Approximately $5\% - 25\%$ of the injected CR energy is lost via hadronic cooling. Smaller diffusion coefficients lead to slightly larger hadronic losses but allow for steeper CR pressure gradients, stronger outflows and larger accelerations. Up to a height of $z \sim1\,\mathrm{kpc}$ there are large volumes in approximate pressure equilibrium between thermal and CR component. At larger altitudes the CR pressure is $10-100$ times as large as the thermal counterpart. More than $\sim 1\,\mathrm{kpc}$ away from the midplane, CRs provide the dominant gas acceleration mechanism.

We explore the effect of momentum-driven winds representing radiation pressure driven outflows from accretion onto supermassive black holes in a set of numerical hydrodynamical simulations. We explore two matched sets of cosmological zoom-in runs of 24 halos with masses ~$10^{12.0}-10^{13.4}$ M_sun run with two different feedback models. Our `NoAGN' model includes stellar feedback via UV heating, stellar winds and supernovae, photoelectric heating and cosmic X-ray background heating from a meta-galactic background. Our fiducial `MrAGN' model is identical except that it also includes a model for black hole seeding and accretion, as well as heating and momentum injection associated with the radiation from black hole accretion. Our MrAGN model launches galactic outflows which result in both `ejective' feedback - the outflows themselves which drive gas out of galaxies - and `preventative' feedback, which suppresses the inflow of new and recycling gas. As much as 80 % of outflowing galactic gas can be expelled, and accretion can be suppressed by as much as a factor of 30 in the MrAGN runs when compared with the NoAGN runs. The histories of NoAGN galaxies are recycling-dominated, with ~70% of material that leaves the galaxy eventually returning, and the majority of outflowing gas re-accretes on 1 Gyr timescales without AGN feedback. Outflowing gas in the MrAGN runs has higher characteristic velocity (500 - 1,000 km/s versus 100-300 km/s for outflowing NoAGN gas) and travels as far as a few Mpcs. Only ~10% of ejected material is re-accreted in the MrAGN galaxies.

Photoionizing radiation and stellar winds from massive stars deposit energy and momentum into the interstellar medium (ISM). They might disperse the local ISM, change its turbulent multi-phase structure, and even regulate star formation. Ionizing radiation dominates the massive stars' energy output, but the relative effect of winds might change with stellar mass and the properties of the ambient ISM. We present simulations of the interaction of stellar winds and ionizing radiation of 12, 23, and 60 M$_{\odot}$ stars within a cold neutral (CNM, $n_{0}$ = 100 cm$^{-3}$), warm neutral (WNM, $n_{0}$ = 1, 10 cm$^{-3}$) or warm ionized (WIM, $n_{0}$ = 0.1 cm$^{-3}$) medium. The FLASH simulations adopt the novel tree-based radiation transfer algorithm TreeRay. With the On-the-Spot approximation and a temperature-dependent recombination coefficient, it is coupled to a chemical network with radiative heating and cooling. In the homogeneous CNM, the total momentum injection ranges from 1.6$\times$10$^{4}$ to 4$\times$10$^{5}$ M$_{\odot}$ km s$^{-1}$ and is always dominated by the expansion of the ionized H$_{\text{II}}$ region. In the WIM, stellar winds dominate (2$\times$10$^{2}$ to 5$\times$10$^{3}$ M$_{\odot}$ km s$^{-1}$), while the input from radiation is small ($\sim$ 10$^{2}$ M$_{\odot}$ km s$^{-1}$). The WNM ($n_{0}$ = 1 cm$^{-3}$) is a transition regime. Energetically, stellar winds couple more efficiently to the ISM ($\sim$ 0.1 percent of wind luminosity) than radiation ($<$ 0.001 percent of ionizing luminosity). For estimating the impact of massive stars, the strongly mass-dependent ratios of wind to ionizing luminosity and the properties of the ambient medium have to be considered.

We present high-resolution ($\sim$ 0.1 pc), hydrodynamical and magnetohydrodynamical simulations to investigate whether the observed level of molecular cloud (MC) turbulence can be generated and maintained by external supernova (SN) explosions. The MCs are formed self-consistently within their large-scale galactic environment following the non-equilibrium formation of H$_2$ and CO including (self-) shielding and important heating and cooling processes. The MCs inherit their initial level of turbulence from the diffuse ISM, where turbulence is injected by SN explosions. However, by systematically exploring the effect of individual SNe going off outside the clouds, we show that at later stages the importance of SN driven turbulence is decreased significantly. This holds for different MC masses as well as MCs with and without magnetic fields. The SN impact also decreases rapidly with larger distances. Nearby SNe ($d$ $\sim$ 25 pc) boost the turbulent velocity dispersions of the MC by up to 70 per cent (up to a few km s$^{-1}$). For $d$ $>$ 50 pc, however, their impact decreases fast with increasing $d$ and is almost negligible. For all probed distances the gain in velocity dispersion decays rapidly within a few 100 kyr. This is significantly shorter than the average timescale for a MC to be hit by a nearby SN under solar neighbourhood conditions ($\sim$ 2 Myr). Hence, at these conditions SNe are not able to sustain the observed level of MC turbulence. However, in environments with high gas surface densities and SN rates like the Central Molecular Zone, observed elevated MC dispersions could be triggered by external SNe.

We present the "SINS/zC-SINF AO survey" of 35 star-forming galaxies, the largest sample with deep adaptive optics-assisted (AO) near-infrared integral field spectroscopy at z~2. The observations, taken with SINFONI at the Very Large Telescope, resolve the Ha and [NII] line emission and kinematics on scales of ~1.5 kpc. In stellar mass, star formation rate, rest-optical colors and size, the AO sample is representative of its parent seeing-limited sample and probes the massive (M* ~ 2x10^9 - 3x10^11 Msun), actively star-forming (SFR ~ 10-600 Msun/yr) part of the z~2 galaxy population over a wide range in colors ((U-V)_rest ~ 0.15-1.5 mag) and half-light radii (R_e,H ~ 1-8.5 kpc). The sample overlaps largely with the "main sequence" of star-forming galaxies in the same redshift range to a similar K_AB = 23 magnitude limit; it has ~0.3 dex higher median specific SFR, ~0.1 mag bluer median (U-V)_rest color, and ~10% larger median rest-optical size. We describe the observations, data reduction, and extraction of basic flux and kinematic properties. With typically 3-4 times higher resolution and 4-5 times longer integrations (up to 23hr) than the seeing-limited datasets of the same objects, the AO data reveal much more detail in morphology and kinematics. The now complete AO observations confirm the majority of kinematically-classified disks and the typically elevated disk velocity dispersions previously reported based on subsets of the data. We derive typically flat or slightly negative radial [NII]/Ha gradients, with no significant trend with global galaxy properties, kinematic nature, or the presence of an AGN. Azimuthal variations in [NII]/Ha are seen in several sources and are associated with ionized gas outflows, and possible more metal-poor star-forming clumps or small companions. [Abridged]

We discuss constraints on the mass density distribution (parameterized as $\rho\propto r^{-\gamma}$) in early-type galaxies provided by strong lensing and stellar kinematics data. The constraints come from mass measurements at two `pinch' radii. One `pinch' radius $r_1=2.2 R_{Einst}$ is defined such that the Einstein (i.e. aperture) mass can be converted to the spherical mass almost independently of the mass-model. Another `pinch' radius $r_2=R_{opt}$ is chosen so that the dynamical mass, derived from the line-of-sight velocity dispersion, is least sensitive to the anisotropy of stellar orbits. We verified the performance of this approach on a sample of simulated elliptical galaxies and on a sample of 15 SLACS lens galaxies at $0.01 \leq z \leq 0.35$, which have already been analysed in Barnabe et al. (2011) by the self-consistent joint lensing and kinematic code. For massive simulated galaxies the density slope $\gamma$ is recovered with an accuracy of $\sim 13\%$, unless $r_1$ and $r_2$ happen to be close to each other. For SLACS galaxies, we found good overall agreement with the results of Barnabe et al. (2011) with a sample-averaged slope $\gamma=2.1\pm0.05$. While the two-pinch-radii approach has larger statistical uncertainties, it is much simpler and uses only few arithmetic operations with directly observable quantities.

This paper provides an update of our previous scaling relations (Genzel et al.2015) between galaxy integrated molecular gas masses, stellar masses and star formation rates, in the framework of the star formation main-sequence (MS), with the main goal to test for possible systematic effects. For this purpose our new study combines three independent methods of determining molecular gas masses from CO line fluxes, far-infrared dust spectral energy distributions, and ~1mm dust photometry, in a large sample of 1444 star forming galaxies (SFGs) between z=0 and 4. The sample covers the stellar mass range log(M*/M_solar)=9.0-11.8, and star formation rates relative to that on the MS, delta_MS=SFR/SFR(MS), from 10^{-1.3} to 10^{2.2}. Our most important finding is that all data sets, despite the different techniques and analysis methods used, follow the same scaling trends, once method-to-method zero point offsets are minimized and uncertainties are properly taken into account. The molecular gas depletion time t_depl, defined as the ratio of molecular gas mass to star formation rate, scales as (1+z)^{-0.6}x(delta_MS)^{-0.44}, and is only weakly dependent on stellar mass. The ratio of molecular-to-stellar mass mu_gas depends on (1+z)^{2.5}x (delta_MS)^{0.52}x(M*)^{-0.36}, which tracks the evolution of the specific star formation rate. The redshift dependence of mu_gas requires a curvature term, as may the mass-dependences of t_depl and mu_gas. We find no or only weak correlations of t_depl and mu_gas with optical size R or surface density once one removes the above scalings, but we caution that optical sizes may not be appropriate for the high gas and dust columns at high-z.

We present a high-resolution smoothed particle hydrodynamics simulation of the Antennae galaxies (NGC 4038/4039) and follow the evolution $3$ Gyrs beyond the final coalescence. The simulation includes metallicity dependent cooling, star formation, and both stellar feedback and chemical enrichment. The simulated best-match Antennae reproduces well both the observed morphology and the off-nuclear starburst. We also produce for the first time a simulated two-dimensional metallicity map of the Antennae and find good agreement with the observed metallicity of off-nuclear stellar clusters, however the nuclear metallicities are overproduced by $\sim 0.5$ dex. Using the radiative transfer code SKIRT we produce multi-wavelength observations of both the Antennae and the merger remnant. The $1$ Gyr old remnant is well fitted with a S\'ersic profile of $n=7.07$, and with an $r$-band effective radius of $r_{\mathrm{e}}= 1.6$ kpc and velocity dispersion of $\sigma_{\mathrm{e}}=180$ km$/$s the remnant is located on the fundamental plane of early-type galaxies (ETGs). The initially blue Antennae remnant evolves onto the red sequence after $\sim 2.5$ Gyr of secular evolution. The remnant would be classified as a fast rotator, as the specific angular momentum evolves from $\lambda_{\mathrm{Re}}\approx0.11$ to $\lambda_{\mathrm{Re}}\approx0.14$ during its evolution. The remnant shows ordered rotation and a double peaked maximum in the mean 2D line-of-sight velocity. These kinematical features are relatively common among local ETGs and we specifically identify three local ETGs (NGC 3226, NGC 3379 and NGC 4494) in the ATLAS$^\mathrm{3D}$ sample, whose photometric and kinematic properties most resemble the Antennae remnant.

We present 3D "zoom-in" simulations of the formation of two molecular clouds out of the galactic interstellar medium. We model the clouds - identified from the SILCC simulations - with a resolution of up to 0.06 pc using adaptive mesh refinement in combination with a chemical network to follow heating, cooling, and the formation of H$_2$ and CO including (self-) shielding. The two clouds are assembled within a few million years with mass growth rates of up to $\sim$ 10$^{-2}$ M$_\mathrm{sun}$ yr$^{-1}$ and final masses of $\sim$ 50 000 M$_\mathrm{sun}$. A spatial resolution of $\lesssim$ 0.1 pc is required for convergence with respect to the mass, velocity dispersion, and chemical abundances of the clouds, although these properties also depend on the cloud definition such as based on density thresholds, H$_2$ or CO mass fraction. To avoid grid artefacts, the progressive increase of resolution has to occur within the free-fall time of the densest structures (1 - 1.5 Myr) and $\gtrsim$ 200 time steps should be spent on each refinement level before the resolution is progressively increased further. This avoids the formation of spurious, large-scale, rotating clumps from unresolved turbulent flows. While CO is a good tracer for the evolution of dense gas with number densities $n \geq$ 300 cm$^{-3}$, H$_2$ is also found for $n \lesssim 30$ cm$^{-3}$ due to turbulent mixing and becomes dominant at column densities around 30 - 50 M$_\mathrm{sun}$ pc$^{-2}$. The CO-to-H$_2$ ratio steadily increases within the first 2 Myr whereas $X_\mathrm{CO}\simeq$ 1 - 4 $\times$ 10$^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$ is approximately constant since the CO(1-0) line quickly becomes optically thick.

We present evidence from cosmological hydrodynamical simulations for a co-evolution of the slope of the total (dark and stellar) mass density profile, gamma_tot, and the dark matter fraction within the half-mass radius, f_DM, in early-type galaxies. The relation can be described as gamma_tot = A f_DM + B for all systems at all redshifts. The trend is set by the decreasing importance of gas dissipation towards lower redshifts and for more massive systems. Early-type galaxies are smaller, more concentrated, have lower f_DM and steeper gamma_tot at high redshifts and at lower masses for a given redshift; f_DM and gamma_tot are good indicators for growth by "dry" merging. The values for A and B change distinctively for different feedback models, and this relation can be used as a test for such models. A similar correlation exists between gamma_tot and the stellar mass surface density Sigma_*. A model with weak stellar feedback and feedback from black holes is in best agreement with observations. All simulations, independent of the assumed feedback model, predict steeper gamma_tot and lower f_DM at higher redshifts. While the latter is in agreement with the observed trends, the former is in conflict with lensing observations, which indicate constant or decreasing gamma_tot. This discrepancy is shown to be artificial: the observed trends can be reproduced from the simulations using observational methodology to calculate the total density slopes.

Galaxies occupy different regions of the [OIII]$\lambda5007$/H$\beta$-versus-[NII]$\lambda6584$/H$\alpha$ emission-line ratio diagram in the distant and local Universe. We investigate the origin of this intriguing result by modelling self-consistently, for the first time, nebular emission from young stars, accreting black holes (BHs) and older, post-asymptotic-giant-branch (post-AGB) stellar populations in galaxy formation simulations in a full cosmological context. In post-processing, we couple new-generation nebular-emission models with high-resolution, cosmological zoom-in simulations of massive galaxies to explore which galaxy physical properties drive the cosmic evolution of the optical-line ratios [OIII]$\lambda5007$/H$\beta$, [NII]$\lambda6584$/H$\alpha$, [SII]$\lambda\lambda6717,6731$/H$\alpha$ and [OI]$\lambda6300$/H$\alpha$. The line ratios of simulated galaxies agree well with observations of both star-forming and active local SDSS galaxies. Towards higher redshifts, at fixed galaxy stellar mass, the average [OIII]/H$\beta$ increases and [NII]/H$\alpha$, [SII]/H$\alpha$ and [OI]/H$\alpha$ decrease -- widely consistent with observations. At fixed stellar mass, we identify star formation history, which controls nebular emission from young stars via the ionization parameter, as the primary driver of the cosmic evolution of [OIII]/H$\beta$ and [NII]/H$\alpha$. For [SII]/H$\alpha$ and [OI]/H$\alpha$, this applies only to redshifts above $z=1.5$, the evolution at lower redshift being driven in roughly equal parts by nebular emission from AGN and post-AGB stars. Instead, changes in the hardness of ionizing radiation, ionized-gas density, the prevalence of BH accretion relative to star formation and the dust-to-metal mass ratio (whose impact on the gas-phase N/O ratio we model at fixed O/H) play at most a minor role in the cosmic evolution of simulated galaxy line ratios.

We observed twelve nearby HI -detected early-type galaxies (ETGs) of stellar mass $\sim 10^{10}M\odot \leq M_* \leq \sim 10^{11}M\odot$ with the Mitchell Integral-Field Spectrograph, reaching approximately three half-light radii in most cases. We extracted line-of-sight velocity distributions for the stellar and gaseous components. We find little evidence of transitions in the stellar kinematics of the galaxies in our sample beyond the central effective radius, with centrally fast-rotating galaxies remaining fast-rotating and centrally slow-rotating galaxies likewise remaining slow-rotating. This is consistent with these galaxies having not experienced late dry major mergers; however, several of our objects have ionised gas that is misaligned with respect to their stars, suggesting some kind of past interaction. We extract Lick index measurements of the commonly-used H$\beta$, Fe5015, Mg\, b, Fe5270 and Fe5335 absorption features, and we find most galaxies to have flat H$\beta$ gradients and negative Mg, b gradients. We measure gradients of age, metallicity and abundance ratio for our galaxies using spectral fitting, and for the majority of our galaxies find negative age and metallicity gradients. We also find the stellar mass-to-light ratios to decrease with radius for most of the galaxies in our sample. Our results are consistent with a view in which intermediate-mass ETGs experience mostly quiet evolutionary histories, but in which many have experienced some kind of gaseous interaction in recent times.

We ask how the inclusion of various physical heating processes due to the metal content of gas affect the evolution of massive galaxies and compute a suite of cosmological hydrodynamical simulations that follow these systems and their supermassive black holes. We use a smoothed particle hydrodynamics code with a pressure-entropy formulation and a more accurate treatment of the metal production, turbulent diffusion and cooling rate based on individual element abundances. The feedback models include (1) AGN feedback via high velocity BAL winds and Compton/photoionization heating, (2) explicit stellar feedback from multiple processes including powerful winds from supernova events, stellar winds from young massive stars and AGB stars as well as radiative heating within Stromgren spheres around massive stars, and (3) additional heating effects due to the presence of metals including grain photoelectric heating, metallicity dependent X-ray heating by nearby accreting black holes and from the cosmic X-ray background, which are the major improvement in our feedback model. With a suite of zoom-in simulations of 30 halos with $M_{vir} \sim 10^{12-13.4}$, we show that energy and momentum budget from all feedback effects generate realistic galaxy properties. We explore the detailed role of each feedback model with three additional sets of simulations with varying input physics. We show that the metal induced heating mechanisms reduce the fraction of accreted stellar material by mainly suppressing the growth of diffuse small stellar systems at high redshift but overall have a relatively minor effect on the final stellar and gas properties of massive galaxies. The inclusion of AGN feedback significantly improves the ability of our cosmological simulations to yield realistic gas and stellar properties of massive galaxies with reasonable fraction of the final stellar mass accreted from other galaxies.

We investigate the stellar mass and baryonic mass Tully-Fisher relations (TFRs) of massive star-forming disk galaxies at redshift z~2.3 and z~0.9 as part of the KMOS^3D integral field spectroscopy survey. Our spatially resolved data allow reliable modelling of individual galaxies, including the effect of pressure support on the inferred gravitational potential. At fixed circular velocity, we find higher baryonic masses and similar stellar masses at z~2.3 as compared to z~0.9. Together with the decreasing gas-to-stellar mass ratios with decreasing redshift, this implies that the contribution of dark matter to the dynamical mass at the galaxy scale increases towards lower redshift. A comparison to local relations reveals a negative evolution of the stellar and baryonic TFR zero-points from z=0 to z~0.9, no evolution of the stellar TFR zero-point from z~0.9 to z~2.3, and a positive evolution of the baryonic TFR zero-point from z~0.9 to z~2.3. We discuss a toy model of disk galaxy evolution to explain the observed, non-monotonic TFR evolution, taking into account the empirically motivated redshift dependencies of galactic gas fractions, and of the relative amount of baryons to dark matter on galaxy and halo scales.

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.

We present KETJU, a new extension of the widely-used smoothed particle hydrodynamics simulation code GADGET-3. The key feature of the code is the inclusion of algorithmically regularized regions around every supermassive black hole (SMBH). This allows for simultaneously following global galactic-scale dynamical and astrophysical processes, while solving the dynamics of SMBHs, SMBH binaries and surrounding stellar systems at sub-parsec scales. The KETJU code includes Post-Newtonian terms in the equations of motions of the SMBHs which enables a new SMBH merger criterion based on the gravitational wave coalescence timescale pushing the merger separation of SMBHs down to $\sim 0.005$ pc. We test the performance of our code by comparison to NBODY7 and rVINE. We set up dynamically stable multi-component merger progenitor galaxies to study the SMBH binary evolution during galaxy mergers. In our simulation sample the SMBH binaries do not suffer from the final-parsec problem, which we attribute to the non-spherical shape of the merger remnants. For bulge-only models, the hardening rate decreases with increasing resolution, whereas for models which in addition include massive dark matter halos the SMBH binary hardening rate becomes practically independent of the mass resolution of the stellar bulge. The SMBHs coalesce on average 200 Myr after the formation of the SMBH binary. However, small differences in the initial SMBH binary eccentricities can result in large differences in the SMBH coalescence times. Finally, we discuss the future prospects of KETJU, which allows for a straightforward inclusion of gas physics in the simulations.

We present three-dimensional hydrodynamical simulations showing the effect of kinetic and radiative AGN feedback on a model galaxy representing a massive quiescent low-redshift early-type galaxy of $M_* = 8.41\times 10^{10} M_\odot$, harbouring a $M_\mathrm{BH} = 4\times 10^8 M_\odot $ black hole surrounded by a cooling gaseous halo. We show that, for a total baryon fraction of $\sim 20\%$ of the cosmological value, feedback from the AGN can keep the galaxy quiescent for about 4.35 Gyr and with properties consistent with black hole mass and X-ray luminosity scaling relations. However, this can only be achieved if the AGN feedback model includes both kinetic and radiative feedback modes. The simulation with only kinetic feedback fails to keep the model galaxy fully quiescent, while one with only radiative feedback leads to excessive black-hole growth. For higher baryon fractions (e.g. 50\% of the cosmological value), the X-ray luminosities exceed observed values by at least one order of magnitude, and rapid cooling results in a star-forming galaxy. The AGN plays a major role in keeping the circumgalactic gas at observed metallicities of $Z/Z_\odot \gtrsim 0.3 $ within the central $\sim 30$ kpc by venting nuclear gas enriched with metals from residual star formation activity. As indicated by previous cosmological simulations, our results are consistent with a model for which the black hole mass and the total baryon fraction are set at higher redshifts $z > 1$ and the AGN alone can keep the model galaxy on observed scaling relations. Models without AGN feedback violate both the quiescence criterion as well as CGM metallicity constraints.

We exploit the deep resolved Halpha kinematic data from the KMOS^3D and SINS/zC-SINF surveys to examine the largely unexplored outer disk kinematics of star-forming galaxies (SFGs) out to the peak of cosmic star formation. Our sample contains 101 SFGs representative of the more massive (9.3 < log(M*/Msun) < 11.5) main sequence population at 0.6<z<2.6. Through a novel stacking approach we are able to constrain a representative rotation curve extending out to ~4 effective radii. This average rotation curve exhibits a significant drop in rotation velocity beyond the turnover, with a slope of Delta(V)/Delta(R) = $-0.26^{+0.10}_{-0.09}$ in units of normalized coordinates V/V_max and R/R_turn. This result confirms that the fall-off seen previously in some individual galaxies is a common feature of our sample of high-z disks. We show that this outer fall-off strikingly deviates from the flat or mildly rising rotation curves of local spiral galaxies of similar masses. We furthermore compare our data with models including baryons and dark matter demonstrating that the falling stacked rotation curve can be explained by a high mass fraction of baryons relative to the total dark matter halo (m_d>~0.05) in combination with a sizeable level of pressure support in the outer disk. These findings are in agreement with recent studies demonstrating that star-forming disks at high redshift are strongly baryon dominated within the disk scale, and furthermore suggest that pressure gradients caused by large turbulent gas motions are present even in their outer disks. We demonstrate that these results are largely independent of our model assumptions such as the presence of a central stellar bulge, the effect of adiabatic contraction at fixed m_d, and variations in the concentration parameter.

In cold dark matter cosmology, the baryonic components of galaxies are thought to be mixed with and embedded in non-baryonic and non-relativistic dark matter, which dominates the total mass of the galaxy and its dark matter halo. In the local Universe, the mass of dark matter within a galactic disk increases with disk radius, becoming appreciable and then dominant in the outer, baryonic regions of the disks of star-forming galaxies. This results in rotation velocities of the visible matter within the disk that are constant or increasing with disk radius. Comparison between the dynamical mass and the sum of stellar and cold gas mass at the peak epoch of galaxy formation, inferred from ancillary data, suggest high baryon factions in the inner, star-forming regions of the disks. Although this implied baryon fraction may be larger than in the local Universe, the systematic uncertainties (stellar initial mass function, calibration of gas masses) render such comparisons inconclusive in terms of the mass of dark matter. Here we report rotation curves for the outer disks of six massive star-forming galaxies, and find that the rotation velocities are not constant, but decrease with radius. We propose that this trend arises because of two main factors: first, a large fraction of the massive, high-redshift galaxy population was strongly baryon dominated, with dark matter playing a smaller part than in the local Universe; and second, the large velocity dispersion in high-redshift disks introduces a substantial pressure term that leads to a decrease in rotation velocity with increasing radius. The effect of both factors appears to increase with redshift. Qualitatively, the observations suggest that baryons in the early Universe efficiently condensed at the centres of dark matter halos when gas fractions were high, and dark matter was less concentrated. [Abridged]

We present high-resolution hydrodynamical simulations of isolated dwarf galaxies including self-gravity, non-equilibrium cooling and chemistry, interstellar radiation fields (IRSF) and shielding, star formation, and stellar feedback. This includes spatially and temporally varying photoelectric (PE) heating, photoionization, resolved supernova (SN) blast waves and metal enrichment. A new flexible method to sample the stellar initial mass function allows us to follow the contribution to the ISRF, the metal output and the SN delay times of individual massive stars. We find that SNe play the dominant role in regulating the global star formation rate, shaping the multi-phase interstellar medium (ISM) and driving galactic outflows. Outflow rates (with mass-loading factors of a few) and hot gas fractions of the ISM increase with the number of SNe exploding in low-density environments where radiative energy losses are low. While PE heating alone can suppress star formation slightly more (a factor of a few) than SNe alone can do, it is unable to drive outflows and reproduce the multi-phase ISM that emerges naturally when SNe are included. These results are in conflict with recent results of Forbes et al. who concluded that PE heating is the dominant process suppressing star formation in dwarfs, about an order of magnitude more efficient than SNe. Potential origins for this discrepancy are discussed. In the absence of SNe and photoionization (mechanisms to disperse dense clouds), the impact of PE heating is highly overestimated owing to the (unrealistic) proximity of dense gas to the radiation sources. This leads to a substantial boost of the infrared continuum emission from the UV-irradiated dust and a far infrared line-to-continuum ratio too low compared to observations. Though sub-dominant in regulating star formation, the ISRF controls the abundance of molecular hydrogen via photodissociation.

Various heuristic approaches to model unresolved supernova (SN) feedback in galaxy formation simulations exist to reproduce the formation of spiral galaxies and the overall inefficient conversion of gas into stars. Some models, however, require resolution dependent scalings. We present a sub-resolution model representing the three major phases of supernova blast wave evolution $-$free expansion, energy conserving Sedov-Taylor, and momentum conserving snowplow$-$ with energy scalings adopted from high-resolution interstellar-medium simulations in both uniform and multiphase media. We allow for the effects of significantly enhanced SN remnant propagation in a multiphase medium with the cooling radius scaling with the hot volume fraction, $f_{\mathrm{hot}}$, as $(1 - f_{\mathrm{hot}})^{-4/5}$. We also include winds from young massive stars and AGB stars, Str\"omgren sphere gas heating by massive stars, and a gas cooling limiting mechanism driven by radiative recombination of dense HII regions. We present initial tests for isolated Milky-Way like systems simulated with the GADGET based code SPHgal with improved SPH prescription. Compared to pure thermal SN input, the model significantly suppresses star formation at early epochs, with star formation extended both in time and space in better accord with observations. Compared to models with pure thermal SN feedback, the age at which half the stellar mass is assembled increases by a factor of 2.4, and the mass loading parameter and gas outflow rate from the galactic disk increase by a factor of 2. Simulation results are converged for a two order of magnitude variation in particle mass in the range (1.3$-$130)$\times 10^4$ solar masses.

Numerical simulations have become a major tool for understanding galaxy formation and evolution. Over the decades the field has made significant progress. It is now possible to simulate the formation of individual galaxies and galaxy populations from well defined initial conditions with realistic abundances and global properties. An essential component of the calculation is to correctly estimate the inflow to and outflow from forming galaxies since observations indicating low formation efficiency and strong circum-glactic presence of gas are persuasive. Energetic 'feedback' from massive stars and accreting super-massive black holes - generally unresolved in cosmological simulations - plays a major role for driving galactic outflows, which have been shown to regulate many aspects of galaxy evolution. A surprisingly large variety of plausible sub-resolution models succeeds in this exercise. They capture the essential characteristics of the problem, i.e. outflows regulating galactic gas flows, but their predictive power is limited. In this review we focus on one major challenge for galaxy formation theory: to understand the underlying physical processes that regulate the structure of the interstellar medium, star formation and the driving of galactic outflows. This requires accurate physical models and numerical simulations, which can precisely describe the multi-phase structure of the interstellar medium on the currently unresolved few hundred parsecs scales of large scale cosmological simulations. Such models ultimately require the full accounting for the dominant cooling and heating processes, the radiation and winds from massive stars and accreting black holes, an accurate treatment of supernova explosions as well as the non-thermal components of the interstellar medium like magnetic fields and cosmic rays.

We present three-dimensional magneto-hydrodynamical simulations of the self-gravitating interstellar medium (ISM) in a periodic (256 pc)$^3$ box with a mean number density of 0.5 cm$^{-3}$. At a fixed supernova rate we investigate the multi-phase ISM structure, H$_{2}$ molecule formation and density-magnetic field scaling for varying initial magnetic field strengths (0, $6\times 10^{-3}$, 0.3, 3 $\mu$G). All magnetic runs saturate at mass weighted field strengths of $\sim$ 1 $-$ 3 $\mu$G but the ISM structure is notably different. With increasing initial field strengths (from $6\times 10^{-3}$ to 3 $\mu$G) the simulations develop an ISM with a more homogeneous density and temperature structure, with increasing mass (from 5% to 85%) and volume filling fractions (from 4% to 85%) of warm (300 K $<$ T $<$ 8000 K) gas, with decreasing volume filling fractions (VFF) from $\sim$ 35% to $\sim$ 12% of hot gas (T $> 10^5$ K) and with a decreasing H$_{2}$ mass fraction (from 70% to $<$ 1%). Meanwhile the mass fraction of gas in which the magnetic pressure dominates over the thermal pressure increases by a factor of 10, from 0.07 for an initial field of $6\times 10^{-3}$ $\mu$G to 0.7 for a 3 $\mu$G initial field. In all but the simulations with the highest initial field strength self-gravity promotes the formation of dense gas and H$_{2}$, but does not change any other trends. We conclude that magnetic fields have a significant impact on the multi-phase, chemical and thermal structure of the ISM and discuss potential implications and limitations of the model.

Dust grains are an important component of the interstellar medium (ISM) of galaxies. We present the first direct measurement of the residence times of interstellar dust in the different ISM phases, and of the transition rates between these phases, in realistic hydrodynamical simulations of the multi-phase ISM. Our simulations include a time-dependent chemical network that follows the abundances of H^+, H, H_2, C^+ and CO and take into account self-shielding by gas and dust using a tree-based radiation transfer method. Supernova explosions are injected either at random locations, at density peaks, or as a mixture of the two. For each simulation, we investigate how matter circulates between the ISM phases and find more sizeable transitions than considered in simple mass exchange schemes in the literature. The derived residence times in the ISM phases are characterised by broad distributions, in particular for the molecular, warm and hot medium. The most realistic simulations with random and mixed driving have median residence times in the molecular, cold, warm and hot phase around 17, 7, 44 and 1 Myr, respectively. The transition rates measured in the random driving run are in good agreement with observations of Ti gas-phase depletion in the warm and cold phases in a simple depletion model, although the depletion in the molecular phase is under-predicted. ISM phase definitions based on chemical abundance rather than temperature cuts are physically more meaningful, but lead to significantly different transition rates and residence times because there is no direct correspondence between the two definitions.

We present three-dimensional radiation-hydrodynamical simulations of the impact of stellar winds, photoelectric heating, photodissociating and photoionising radiation, and supernovae on the chemical composition and star formation in a stratified disc model. This is followed with a sink-based model for star clusters with populations of individual massive stars. Stellar winds and ionising radiation regulate the star formation rate at a factor of ~10 below the simulation with only supernova feedback due to their immediate impact on the ambient interstellar medium after star formation. Ionising radiation (with winds and supernovae) significantly reduces the ambient densities for most supernova explosions to rho < 10^-25 g cm^-3, compared to 10^-23 g cm^-3 for the model with only winds and supernovae. Radiation from massive stars reduces the amount of molecular hydrogen and increases the neutral hydrogen mass and volume filling fraction. Only this model results in a molecular gas depletion time scale of 2 Gyr and shows the best agreement with observations. In the radiative models, the Halpha emission is dominated by radiative recombination as opposed to collisional excitation (the dominant emission in non-radiative models), which only contributes ~1-10 % to the total Halpha emission. Individual massive stars (M >= 30 M_sun) with short lifetimes are responsible for significant fluctuations in the Halpha luminosities. The corresponding inferred star formation rates can underestimate the true instantaneous star formation rate by factors of ~10.

We exploit deep integral-field spectroscopic observations with KMOS/VLT of 240 star-forming disks at 0.6 < z < 2.6 to dynamically constrain their mass budget. Our sample consists of massive ($\gtrsim 10^{9.8} M_\odot$) galaxies with sizes $R_e \gtrsim 2$ kpc. By contrasting the observed velocity and dispersion profiles to dynamical models, we find that on average the stellar content contributes $32^{+8}_{-7}\%$ of the total dynamical mass, with a significant spread among galaxies (68th percentile range f_star ~ 18 - 62%). Including molecular gas as inferred from CO- and dust-based scaling relations, the estimated baryonic mass adds up to $56^{+17}_{-12}\%$ of total for the typical galaxy in our sample, reaching ~ 90% at z > 2. We conclude that baryons make up most of the mass within the disk regions of high-redshift star-forming disk galaxies, with typical disks at z > 2 being strongly baryon-dominated within $R_e$. Substantial object-to-object variations in both stellar and baryonic mass fractions are observed among the galaxies in our sample, larger than what can be accounted for by the formal uncertainties in their respective measurements. In both cases, the mass fractions correlate most strongly with measures of surface density. High $\Sigma_{star}$ galaxies feature stellar mass fractions closer to unity, and systems with high inferred gas or baryonic surface densities leave less room for additional mass components other than stars and molecular gas. Our findings can be interpreted as more extended disks probing further (and more compact disks probing less far) into the dark matter halos that host them.

We study the impact of stellar winds and supernovae on the multi-phase interstellar medium using three-dimensional hydrodynamical simulations carried out with FLASH. The selected galactic disc region has a size of (500 pc)$^2$ x $\pm$ 5 kpc and a gas surface density of 10 M$_{\odot}$/pc$^2$. The simulations include an external stellar potential and gas self-gravity, radiative cooling and diffuse heating, sink particles representing star clusters, stellar winds from these clusters which combine the winds from indi- vidual massive stars by following their evolution tracks, and subsequent supernova explosions. Dust and gas (self-)shielding is followed to compute the chemical state of the gas with a chemical network. We find that stellar winds can regulate star (cluster) formation. Since the winds suppress the accretion of fresh gas soon after the cluster has formed, they lead to clusters which have lower average masses (10$^2$ - 10$^{4.3}$ M$_{\odot}$) and form on shorter timescales (10$^{-3}$ - 10 Myr). In particular we find an anti-correlation of cluster mass and accretion time scale. Without winds the star clusters easily grow to larger masses for ~5 Myr until the first supernova explodes. Overall the most massive stars provide the most wind energy input, while objects beginning their evolution as B-type stars contribute most of the supernova energy input. A significant outflow from the disk (mass loading $\gtrsim$ 1 at 1 kpc) can be launched by thermal gas pressure if more than 50% of the volume near the disc mid-plane can be heated to T > 3x10$^5$ K. Stellar winds alone cannot create a hot volume-filling phase. The models which are in best agreement with observed star formation rates drive either no outflows or weak outflows.

We analyze the angular momenta of massive star forming galaxies (SFGs) at the peak of the cosmic star formation epoch (z~0.8-2.6). Our sample of ~360 log(M*/Msun) ~ 9.3-11.8 SFGs is mainly based on the KMOS3D and SINS/zC-SINF surveys of H$\alpha$ kinematics, and collectively provides a representative subset of the massive star forming population. The inferred halo scale angular momentum distribution is broadly consistent with that theoretically predicted for their dark matter halos, in terms of mean spin parameter <$\lambda$> ~ 0.037 and its dispersion ($\sigma_{log(\lambda)}$~0.2). Spin parameters correlate with the disk radial scale, and with their stellar surface density, but do not depend significantly on halo mass, stellar mass, or redshift. Our data thus support the long-standing assumption that on average, even at high redshifts, the specific angular momentum of disk galaxies reflects that of their dark matter halos (j_d = j_DM). The lack of correlation between $\lambda$ x (j_d/j_DM) and the nuclear stellar density $\Sigma_{*}$(1kpc) favors a scenario where disk-internal angular momentum redistribution leads to "compaction" inside massive high-redshift disks. For our sample, the inferred average stellar-to-dark matter mass ratio is ~2%, consistent with abundance matching results. Including the molecular gas, the total baryonic disk-to-dark matter mass ratio is ~5% for halos near $10^{12}$ Msun, which corresponds to 31% of the cosmologically available baryons, implying that high-redshift disks are strongly baryon dominated.

We observed the lenticular galaxy NGC 3998 with the Mitchell Integral-Field Spectrograph and extracted line-of-sight velocity distributions out to 3 half-light radii. We constructed collisionless orbit models in order to constrain NGC 3998's dark and visible structure, using kinematics from both the Mitchell and SAURON instruments. We find NGC 3998 to be almost axisymmetric, seen nearly face on with a flattened intrinsic shape - i.e., a face-on fast-rotator. We find an I-band mass-to-light ratio of $4.7_{-0.45}^{+0.32}$ in good agreement with previous spectral fitting results for this galaxy. Our best-fit orbit model shows a both a bulge and a disc component, with a non-negligible counter-rotating component also evident. We find that relatively little dark matter is needed to model this galaxy, with an inferred dark mass fraction of just $(7.1^{+8.1}_{-7.1})\%$ within one half-light radius.

Supernova (SN) blast waves inject energy and momentum into the interstellar medium (ISM), control its turbulent multiphase structure and the launching of galactic outflows. Accurate modelling of the blast wave evolution is therefore essential for ISM and galaxy formation simulations. We present an efficient method to compute the input of momentum, thermal energy, and the velocity distribution of the shock-accelerated gas for ambient media with uniform (and with stellar wind blown bubbles), power-law, and turbulent density distributions. Assuming solar metallicity cooling, the blast wave evolution is followed to the beginning of the momentum conserving snowplough phase. The model recovers previous results for uniform ambient media. The momentum injection in wind-blown bubbles depend on the swept-up mass and the efficiency of cooling, when the blast wave hits the wind shell. For power-law density distributions with $n(r) \sim$ $r^{-2}$ (for $n(r) > n_{_{\rm floor}}$) the amount of momentum injection is solely regulated by the background density $n_{_{\rm floor}}$ and compares to $n_{_{\rm uni}}$ = $n_{_{\rm floor}}$. However, in turbulent ambient media with log-normal density distributions the momentum input can increase by a factor of 2 (compared to the homogeneous case) for high Mach numbers. The average momentum boost can be approximated as $p_{_{\rm turb}}/\mathrm{p_{_{0}}}\ =23.07\, \left(\frac{n_{_{0,\rm turb}}}{1\,{\rm cm}^{-3}}\right)^{-0.12} + 0.82 (\ln(1+b^{2}\mathcal{M}^{2}))^{1.49}\left(\frac{n_{_{0,\rm turb}}}{1\,{\rm cm}^{-3}}\right)^{-1.6}$. The velocity distributions are broad as gas can be accelerated to high velocities in low-density channels. The model values agree with results from recent, computationally expensive, three-dimensional simulations of SN explosions in turbulent media.

We study the connection of star formation to atomic (HI) and molecular hydrogen (H$_2$) in isolated, low metallicity dwarf galaxies with high-resolution ($m_{\rm gas}$ = 4 M$_\odot$, $N_{\rm ngb}$ = 100) SPH simulations. The model includes self-gravity, non-equilibrium cooling, shielding from an interstellar radiation field, the chemistry of H$_2$ formation, H$_2$-independent star formation, supernova feedback and metal enrichment. We find that the H$_2$ mass fraction is sensitive to the adopted dust-to-gas ratio and the strength of the interstellar radiation field, while the star formation rate is not. Star formation is regulated by stellar feedback, keeping the gas out of thermal equilibrium for densities $n <$ 1 cm$^{-3}$. Because of the long chemical timescales, the H$_2$ mass remains out of chemical equilibrium throughout the simulation. Star formation is well-correlated with cold ( T $\leqslant$ 100 K ) gas, but this dense and cold gas - the reservoir for star formation - is dominated by HI, not H$_2$. In addition, a significant fraction of H$_2$ resides in a diffuse, warm phase, which is not star-forming. The ISM is dominated by warm gas (100 K $<$ T $\leqslant 3\times 10^4$ K) both in mass and in volume. The scale height of the gaseous disc increases with radius while the cold gas is always confined to a thin layer in the mid-plane. The cold gas fraction is regulated by feedback at small radii and by the assumed radiation field at large radii. The decreasing cold gas fractions result in a rapid increase in depletion time (up to 100 Gyrs) for total gas surface densities $\Sigma_{\rm HI+H_2} \lesssim$ 10 M$_\odot$pc$^{-2}$, in agreement with observations of dwarf galaxies in the Kennicutt-Schmidt plane.

We present measurements of the [NII]/Ha ratio as a probe of gas-phase oxygen abundance for a sample of 419 star-forming galaxies at z=0.6-2.7 from the KMOS3D near-IR multi-IFU survey. The mass-metallicity relation (MZR) is determined consistently with the same sample selection, metallicity tracer, and methodology over the wide redshift range probed by the survey. We find good agreement with long-slit surveys in the literature, except for the low-mass slope of the relation at z~2.3, where this sample is less biased than previous samples based on optical spectroscopic redshifts. In this regime we measure a steeper slope than some literature results. Excluding the AGN contribution from the MZR reduces sensitivity at the high mass end, but produces otherwise consistent results. There is no significant dependence of the [NII]/Ha ratio on SFR or environment at fixed redshift and stellar mass. The IFU data allow spatially resolved measurements of [NII]/Ha, from which we can infer abundance gradients for 180 galaxies, thus tripling the current sample in the literature. The observed gradients are on average flat, with only 15 gradients statistically offset from zero at >3sigma. We have modelled the effect of beam-smearing, assuming a smooth intrinsic radial gradient and known seeing, inclination and effective radius for each galaxy. Our seeing-limited observations can recover up to 70% of the intrinsic gradient for the largest, face-on disks, but only 30% for the smaller, more inclined galaxies. We do not find significant trends between observed or corrected gradients and any stellar population, dynamical or structural galaxy parameters, mostly in agreement with existing studies with much smaller sample sizes. In cosmological simulations, strong feedback is generally required to produce flat gradients at high redshift.

We present the results of a high-resolution, 5 GHz, Karl G. Jansky Very Large Array study of the nuclear radio emission in a representative subset of the Atlas3D survey of early-type galaxies (ETGs). We find that 51 +/- 4% of the ETGs in our sample contain nuclear radio emission with luminosities as low as 10^18 W/Hz. Most of the nuclear radio sources have compact (< 25-110 pc) morphologies, although < 10% display multi-component core+jet or extended jet/lobe structures. Based on the radio continuum properties, as well as optical emission line diagnostics and the nuclear X-ray properties, we conclude that the majority of the central 5 GHz sources detected in the Atlas3D galaxies are associated with the presence of an active galactic nucleus (AGN). However, even at sub-arcsecond spatial resolution, the nuclear radio emission in some cases appears to arise from low-level nuclear star formation rather than an AGN, particularly when molecular gas and a young central stellar population is present. This is in contrast to popular assumptions in the literature that the presence of a compact, unresolved, nuclear radio continuum source universally signifies the presence of an AGN. Additionally, we examine the relationships between the 5 GHz luminosity and various galaxy properties including the molecular gas mass and - for the first time - the global kinematic state. We discuss implications for the growth, triggering, and fueling of radio AGNs, as well as AGN-driven feedback in the continued evolution of nearby ETGs.

Introducing the DRAGON simulation project, we present direct $N$-body simulations of four massive globular clusters (GCs) with $10^6$ stars and 5$\%$ primordial binaries at a high level of accuracy and realism. The GC evolution is computed with NBODY6++GPU and follows the dynamical and stellar evolution of individual stars and binaries, kicks of neutron stars and black holes, and the effect of a tidal field. We investigate the evolution of the luminous (stellar) and dark (faint stars and stellar remnants) GC components and create mock observations of the simulations (i.e. photometry, color-magnitude diagrams, surface brightness and velocity dispersion profiles). By connecting internal processes to observable features we highlight the formation of a long-lived 'dark' nuclear subsystem made of black holes (BHs), which results in a two-component structure. The inner core is dominated by the BH subsystem and experiences a core collapse phase within the first Gyr. It can be detected in the stellar (luminous) line-of-sight velocity dispersion profiles. The outer extended core - commonly observed in the (luminous) surface brightness profiles - shows no collapse features and is continuously expanding. We demonstrate how a King (1966) model fit to observed clusters might help identify the presence of post core-collapse BH subsystems. For global observables like core and half-mass radii the direct simulations agree well with Monte-Carlo models. Variations in the initial mass function can result in significantly different GC properties (e.g. density distributions) driven by varying amounts of early mass loss and the number of forming BHs.

We present a hydrodynamical simulation of the turbulent, magnetized, supernova (SN)-driven interstellar medium (ISM) in a stratified box that dynamically couples the injection and evolution of cosmic rays (CRs) and a self-consistent evolution of the chemical composition. CRs are treated as a relativistic fluid in the advection-diffusion approximation. The thermodynamic evolution of the gas is computed using a chemical network that follows the abundances of H+, H, H2, CO, C+, and free electrons and includes (self-)shielding of the gas and dust. We find that CRs perceptibly thicken the disk with the heights of 90% (70%) enclosed mass reaching ~1.5 kpc (~0.2 kpc). The simulations indicate that CRs alone can launch and sustain strong outflows of atomic and ionized gas with mass loading factors of order unity, even in solar neighborhood conditions and with a CR energy injection per SN of 10^50 erg, 10% of the fiducial thermal energy of an SN. The CR-driven outflows have moderate launching velocities close to the midplane (~100 km/s) and are denser (\rho~1e-24 - 1e-26 g/cm^3), smoother, and colder than the (thermal) SN-driven winds. The simulations support the importance of CRs for setting the vertical structure of the disk as well as the driving of winds.

The SILCC project (SImulating the Life-Cycle of molecular Clouds) aims at a more self-consistent understanding of the interstellar medium (ISM) on small scales and its link to galaxy evolution. We present three-dimensional (magneto)hydrodynamic simulations of the ISM in a vertically stratified box including self-gravity, an external potential due to the stellar component of the galactic disc, and stellar feedback in the form of an interstellar radiation field and supernovae (SNe). The cooling of the gas is based on a chemical network that follows the abundances of H+, H, H2, C+, and CO and takes shielding into account consistently. We vary the SN feedback by comparing different SN rates, clustering and different positioning, in particular SNe in density peaks and at random positions, which has a major impact on the dynamics. Only for random SN positions the energy is injected in sufficiently low-density environments to reduce energy losses and enhance the effective kinetic coupling of the SNe with the gas. This leads to more realistic velocity dispersions (\sigma_HI ~ 0.8\sigma_(300-8000K) ~ 10-20km/s, \sigma_H\alpha ~ 0.6\sigma_(8000-3e5K) ~ 20-30km/s), and strong outflows with mass loading factors of up to 10 even for solar neighbourhood conditions. Clustered SNe abet the onset of outflows compared to individual SNe but do not influence the net outflow rate. The outflows do not contain any molecular gas and are mainly composed of atomic hydrogen. The bulk of the outflowing mass is dense (\rho ~ 1e-25-1e-24g/cc) and slow (v ~ 20-40km/s) but there is a high-velocity tail of up to v ~ 500km/s with \rho ~ 1e-28-1e-27g/cc.

The halo of the Milky Way contains a hot plasma with a surface brightness in soft X-rays of the order $10^{-12}$erg cm$^{-2}$ s$^{-1}$ deg$^{-2}$. The origin of this gas is unclear, but so far numerical models of galactic star formation have failed to reproduce such a large surface brightness by several orders of magnitude. In this paper, we analyze simulations of the turbulent, magnetized, multi-phase interstellar medium including thermal feedback by supernova explosions as well as cosmic-ray feedback. We include a time-dependent chemical network, self-shielding by gas and dust, and self-gravity. Pure thermal feedback alone is sufficient to produce the observed surface brightness, although it is very sensitive to the supernova rate. Cosmic rays suppress this sensitivity and reduce the surface brightness because they drive cooler outflows. Self-gravity has by far the largest effect because it accumulates the diffuse gas in the disk in dense clumps and filaments, so that supernovae exploding in voids can eject a large amount of hot gas into the halo. This can boost the surface brightness by several orders of magnitude. Although our simulations do not reach a steady state, all simulations produce surface brightness values of the same order of magnitude as the observations, with the exact value depending sensitively on the simulation parameters. We conclude that star formation feedback alone is sufficient to explain the origin of the hot halo gas, but measurements of the surface brightness alone do not provide useful diagnostics for the study of galactic star formation.

Supernovae (SN), the most energetic stellar feedback mechanism, are crucial for regulating the interstellar medium (ISM) and launching galactic winds. We explore how supernova remnants (SNRs) create a multiphase medium by performing 3D hydrodynamical simulations at various SN rates, $S$, and ISM average densities, $\bar{n}$. The evolution of a SNR in a self-consistently generated three-phase ISM is qualitatively different from that in a uniform or a two-phase warm/cold medium. By travelling faster and further in the low-density hot phase, the domain of a SNR increases by $>10^{2.5}$. Varying $\bar{n}$ and $S$, we find that a steady state can only be achieved when the hot gas volume fraction $f_{\rm{V,hot}}\lesssim 0.6 \pm 0.1 $. Above that level, overlapping SNRs render connecting topology of the hot gas, and the ISM is subjected to thermal runaway. Photoelectric heating (PEH) has a surprisingly strong impact on $f_{\rm{V,hot}}$. For $\bar{n}\gtrsim 3 \cm-3 $, a reasonable PEH rate is able to suppress the thermal runaway. Overall, we determine the critical SN rate for the onset of thermal runaway to be $S_{\rm{crit}} = 200 (\bar{n}/1\cm-3)^k (E_{\rm{SN}}/10^{51}\erg)^{-1} \kpc^{-3} \myr-1$, where $k = (1.2,2.7)$ for $\bar{n} \leq 1$ and $> 1\cm-3 $, respectively. We present a fitting formula of the ISM pressure $P(\bar{n}$, $S$), which can be used as an effective equation of state in cosmological simulations. Despite the 5 orders of magnitude span of $(\bar{n},S)$, the average Mach number varies little: $\mathcal{M} \approx \ 0.5\pm 0.2, \ 1.2\pm 0.3,\ 2.3\pm 0.9$ for the hot, warm and cold phases, respectively.

We study the HI K-band Tully-Fisher relation and the baryonic Tully-Fisher relation for a sample of 16 early-type galaxies, taken from the ATLAS3D sample, which all have very regular HI disks extending well beyond the optical body (> 5 R_eff). We use the kinematics of these disks to estimate the circular velocity at large radii for these galaxies. We find that the Tully-Fisher relation for our early-type galaxies is offset by about 0.5-0.7 magnitudes from the relation for spiral galaxies. The residuals with respect to the spiral Tully-Fisher relation correlate with estimates of the stellar mass-to-light ratio, suggesting that the offset between the relations is mainly driven by differences in stellar populations. We also observe a small offset between our Tully-Fisher relation with the relation derived for the ATLAS3D sample based on CO data representing the galaxies' inner regions (< 1 R_eff). This indicates that the circular velocities at large radii are systematically 10% lower than those near 0.5-1 R_eff, in line with recent determinations of the shape of the mass profile of early-type galaxies. The baryonic Tully-Fisher relation of our sample is distinctly tighter than the standard one, in particular when using mass-to-light ratios based on dynamical models of the stellar kinematics. We find that the early-type galaxies fall on the spiral baryonic Tully-Fisher relation if one assumes M/L_K = 0.54 M_sun/L_sun for the stellar populations of the spirals, a value similar to that found by recent studies of the dynamics of spiral galaxies. Such a mass-to-light ratio for spiral galaxies would imply that their disks are 60-70% of maximal. Our analysis increases the range of galaxy morphologies for which the baryonic Tully-Fisher relations holds, strengthening previous claims that it is a more fundamental scaling relation than the classical Tully-Fisher relation.

We present a hybrid code combining the OpenMP-parallel tree code VINE with an algorithmic chain regularization scheme. The new code, called "rVINE", aims to significantly improve the accuracy of close encounters of massive bodies with supermassive black holes in galaxy-scale numerical simulations. We demonstrate the capabilities of the code by studying two test problems, the sinking of a single massive black hole to the centre of a gas-free galaxy due to dynamical friction and the hardening of a supermassive black hole binary due to close stellar encounters. We show that results obtained with rVINE compare well with NBODY7 for problems with particle numbers that can be simulated with NBODY7. In particular, in both NBODY7 and rVINE we find a clear N-dependence of the binary hardening rate, a low binary eccentricity and moderate eccentricity evolution, as well as the conversion of the galaxy's inner density profile from a cusp to a a core via the ejection of stars at high velocity. The much larger number of particles that can be handled by rVINE will open up exciting opportunities to model stellar dynamics close to SMBHs much more accurately in a realistic galactic context. This will help to remedy the inherent limitations of commonly used tree solvers to follow the correct dynamical evolution of black holes in galaxy scale simulations.

Accurate direct $N$-body simulations help to obtain detailed information about the dynamical evolution of star clusters. They also enable comparisons with analytical models and Fokker-Planck or Monte-Carlo methods. NBODY6 is a well-known direct $N$-body code for star clusters, and NBODY6++ is the extended version designed for large particle number simulations by supercomputers. We present NBODY6++GPU, an optimized version of NBODY6++ with hybrid parallelization methods (MPI, GPU, OpenMP, and AVX/SSE) to accelerate large direct $N$-body simulations, and in particular to solve the million-body problem. We discuss the new features of the NBODY6++GPU code, benchmarks, as well as the first results from a simulation of a realistic globular cluster initially containing a million particles. For million-body simulations, NBODY6++GPU is $400-2000$ times faster than NBODY6 with 320 CPU cores and 32 NVIDIA K20X GPUs. With this computing cluster specification, the simulations of million-body globular clusters including $5\%$ primordial binaries require about an hour per half-mass crossing time.

We compare the performance of mass estimators for elliptical galaxies that rely on the directly observable surface brightness and velocity dispersion profiles, without invoking computationally expensive detailed modeling. These methods recover the mass at a specific radius where the mass estimate is expected to be least sensitive to the anisotropy of stellar orbits. One method (Wolf et al. 2010) uses the total luminosity-weighted velocity dispersion and evaluates the mass at a 3D half-light radius $r_{1/2}$, i.e., it depends on the GLOBAL galaxy properties. Another approach (Churazov et al. 2010) estimates the mass from the velocity dispersion at a radius $R_2$ where the surface brightness declines as $R^{-2}$, i.e., it depends on the LOCAL properties. We evaluate the accuracy of the two methods for analytical models, simulated galaxies and real elliptical galaxies that have already been modeled by the Schwarzschild's orbit-superposition technique. Both estimators recover an almost unbiased circular speed estimate with a modest RMS scatter ($\lesssim 10 \%$). Tests on analytical models and simulated galaxies indicate that the local estimator has a smaller RMS scatter than the global one. We show by examination of simulated galaxies that the projected velocity dispersion at $R_2$ could serve as a good proxy for the virial galaxy mass. For simulated galaxies the total halo mass scales with $\sigma_p(R_2)$ as $M_{vir} \left[M_{\odot}h^{-1}\right] \approx 6\cdot 10^{12} \left( \frac{\sigma_p(R_2)}{200\, \rm km\, s^{-1}} \right)^{4}$ with RMS scatter $\approx 40 \%$.

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 use hydrodynamical simulations in a $(256\;{\rm pc})^3$ periodic box to model the impact of supernova (SN) explosions on the multi-phase interstellar medium (ISM) for initial densities $n = 0.5-30$ cm$^{-3}$ and SN rates $1-720$ Myr$^{-1}$. We include radiative cooling, diffuse heating, and the formation of molecular gas using a chemical network. The SNe explode either at random positions, at density peaks, or both. We further present a model combining thermal energy for resolved and momentum input for unresolved SNe. Random driving at high SN rates results in hot gas ($T\gtrsim 10^6$ K) filling $> 90$% of the volume. This gas reaches high pressures ($10^4 < P/k_\mathrm{B} < 10^7$ K cm$^{-3}$) due to the combination of SN explosions in the hot, low density medium and confinement in the periodic box. These pressures move the gas from a two-phase equilibrium to the single-phase, cold branch of the cooling curve. The molecular hydrogen dominates the mass ($>50$%), residing in small, dense clumps. Such a model might resemble the dense ISM in high-redshift galaxies. Peak driving results in huge radiative losses, producing a filamentary ISM with virtually no hot gas, and a small molecular hydrogen mass fraction ($\ll 1$%). Varying the ratio of peak to random SNe yields ISM properties in between the two extremes, with a sharp transition for equal contributions. The velocity dispersion in HI remains $\lesssim 10$ km s$^{-1}$ in all cases. For peak driving the velocity dispersion in H$_\alpha$ can be as high as $70$ km s$^{-1}$ due to the contribution from young, embedded SN remnants.

We investigate the evolution of stellar population gradients from $z=2$ to $z=0$ in massive galaxies at large radii ($r > 2R_{\mathrm{eff}}$) using ten cosmological zoom simulations of halos with $6 \times 10^{12} M_{\odot} < M_{\mathrm{halo}} < 2 \times 10^{13}M_{\odot}$. The simulations follow metal cooling and enrichment from SNII, SNIa and AGB winds. We explore the differential impact of an empirical model for galactic winds that reproduces the mass-metallicity relation and its evolution with redshift. At larger radii the galaxies, for both models, become more dominated by stars accreted from satellite galaxies in major and minor mergers. In the wind model, fewer stars are accreted, but they are significantly more metal poor resulting in steep global metallicity ($\langle \nabla Z_{\mathrm{stars}} \rangle= -0.35$ dex/dex) and color (e.g. $\langle \nabla g-r \rangle = -0.13$ dex/dex) gradients in agreement with observations. In contrast, colour and metallicity gradients of the models without winds are inconsistent with observations. Age gradients are in general mildly positive at $z=0$ ($\langle \nabla Age_{\mathrm{stars}} \rangle= 0.04$ dex/dex) with significant differences between the models at higher redshift. We demonstrate that for the wind model, stellar accretion is steepening existing in-situ metallicity gradients by about 0.2 dex by the present day and helps to match observed gradients of massive early-type galaxies at large radii. Colour and metallicity gradients are significantly steeper for systems which have accreted stars in minor mergers, while galaxies with major mergers have relatively flat gradients, confirming previous results. This study highlights the importance of stellar accretion for stellar population properties of massive galaxies at large radii, which can provide important constraints for formation models.

We present the stellar population content of early-type galaxies from the Atlas3D survey. Using spectra integrated within apertures covering up to one effective radius, we apply two methods: one based on measuring line-strength indices and applying single stellar population (SSP) models to derive SSP-equivalent values of stellar age, metallicity, and alpha enhancement; and one based on spectral fitting to derive non-parametric star-formation histories, mass-weighted average values of age, metallicity, and half-mass formation timescales. Using homogeneously derived effective radii and dynamically-determined galaxy masses, we present the distribution of stellar population parameters on the Mass Plane (M_JAM, Sigma_e, R_maj), showing that at fixed mass, compact early-type galaxies are on average older, more metal-rich, and more alpha-enhanced than their larger counterparts. From non-parametric star-formation histories, we find that the duration of star formation is systematically more extended in lower mass objects. Assuming that our sample represents most of the stellar content of today's local Universe, approximately 50% of all stars formed within the first 2 Gyr following the big bang. Most of these stars reside today in the most massive galaxies (>10^10.5 M_sun), which themselves formed 90% of their stars by z~2. The lower-mass objects, in contrast, have formed barely half their stars in this time interval. Stellar population properties are independent of environment over two orders of magnitude in local density, varying only with galaxy mass. In the highest-density regions of our volume (dominated by the Virgo cluster), galaxies are older, alpha-enhanced and have shorter star-formation histories with respect to lower density regions.

We analyze orbits of stars and dark matter out to three effective radii for 42 galaxies formed in cosmological zoom simulations. Box orbits always dominate at the centers and $z$-tubes become important at larger radii. We connect the orbital structure to the formation histories and specific features (e.g. disk, counter-rotating core, minor axis rotation) in two-dimensional kinematic maps. Globally, fast rotating galaxies with significant recent in situ star formation are dominated by $z$-tubes. Slow rotators with recent mergers have significant box orbit and $x$-tube components. Rotation, quantified by the $\lambda_R$-parameter often originates from streaming motion of stars on $z$-tubes but sometimes from figure rotation. The observed anti-correlation of $h_3$ and $V_0 / \sigma$ in rotating galaxies can be connected to a dissipative formation history leading to high $z$-tube fractions. For galaxies with recent mergers in situ formed stars, accreted stars and dark matter particles populate similar orbits. Dark matter particles have isotropic velocity dispersions. Accreted stars are typically radially biased ($\beta \approx 0.2 - 0.4$). In situ stars become tangentially biased (as low as $\beta \approx -1.0$) if dissipation was relevant during the late assembly of the galaxy. We discuss the relevance of our analysis for integral field surveys and for constraining galaxy formation models.

The SILCC project (SImulating the Life-Cycle of molecular Clouds) aims at a more self-consistent understanding of the interstellar medium (ISM) on small scales and its link to galaxy evolution. We simulate the evolution of the multi-phase ISM in a 500 pc x 500 pc x 10 kpc region of a galactic disc, with a gas surface density of $\Sigma_{_{\rm GAS}} = 10 \;{\rm M}_\odot/{\rm pc}^2$. The Flash 4.1 simulations include an external potential, self-gravity, magnetic fields, heating and radiative cooling, time-dependent chemistry of H$_2$ and CO considering (self-) shielding, and supernova (SN) feedback. We explore SN explosions at different (fixed) rates in high-density regions (peak), in random locations (random), in a combination of both (mixed), or clustered in space and time (clustered). Only random or clustered models with self-gravity (which evolve similarly) are in agreement with observations. Molecular hydrogen forms in dense filaments and clumps and contributes 20% - 40% to the total mass, whereas most of the mass (55% - 75%) is in atomic hydrogen. The ionised gas contributes <10%. For high SN rates (0.5 dex above Kennicutt-Schmidt) as well as for peak and mixed driving the formation of H$_2$ is strongly suppressed. Also without self-gravity the H$_2$ fraction is significantly lower ($\sim$ 5%). Most of the volume is filled with hot gas ($\sim$90% within $\pm$2 kpc). Only for random or clustered driving, a vertically expanding warm component of atomic hydrogen indicates a fountain flow. Magnetic fields have little impact on the final disc structure. However, they affect dense gas ($n\gtrsim 10\;{\rm cm}^{-3}$) and delay H$_2$ formation. We highlight that individual chemical species, in particular atomic hydrogen, populate different ISM phases and cannot be accurately accounted for by simple temperature-/density-based phase cut-offs.

We combine molecular gas masses inferred from CO emission in 500 star forming galaxies (SFGs) between z=0 and 3, from the IRAM-COLDGASS, PHIBSS1/2 and other surveys, with gas masses derived from Herschel far-IR dust measurements in 512 galaxy stacks over the same stellar mass/redshift range. We constrain the scaling relations of molecular gas depletion time scale (tdepl) and gas to stellar mass ratio (Mmolgas/M*) of SFGs near the star formation main-sequence with redshift, specific star formation rate (sSFR) and stellar mass (M*). The CO- and dust-based scaling relations agree remarkably well. This suggests that the CO-H2 mass conversion factor varies little within 0.6dex of the main sequence (sSFR(ms,z,M*)), and less than 0.3dex throughout this redshift range. This study builds on and strengthens the results of earlier work. We find that tdepl scales as (1+z)^-0.3 *(sSFR/sSFR(ms,z,M*))^-0.5, with little dependence on M*. The resulting steep redshift dependence of Mmolgas/M* ~(1+z)^3 mirrors that of the sSFR and probably reflects the gas supply rate. The decreasing gas fractions at high M* are driven by the flattening of the SFR-M* relation. Throughout the redshift range probed a larger sSFR at constant M* is due to a combination of an increasing gas fraction and a decreasing depletion time scale. As a result galaxy integrated samples of the Mmolgas-SFR rate relation exhibit a super-linear slope, which increases with the range of sSFR. With these new relations it is now possible to determine Mmolgas with an accuracy of 0.1dex in relative terms, and 0.2dex including systematic uncertainties.

Galactic archeology based on star counts is instrumental to reconstruct the past mass assembly of Local Group galaxies. The development of new observing techniques and data-reduction, coupled with the use of sensitive large field of view cameras, now allows us to pursue this technique in more distant galaxies exploiting their diffuse low surface brightness (LSB) light. As part of the Atlas3D project, we have obtained with the MegaCam camera at the Canada-France Hawaii Telescope extremely deep, multi--band, images of nearby early-type galaxies. We present here a catalog of 92 galaxies from the Atlas3D sample, that are located in low to medium density environments. The observing strategy and data reduction pipeline, that achieve a gain of several magnitudes in the limiting surface brightness with respect to classical imaging surveys, are presented. The size and depth of the survey is compared to other recent deep imaging projects. The paper highlights the capability of LSB--optimized surveys at detecting new prominent structures that change the apparent morphology of galaxies. The intrinsic limitations of deep imaging observations are also discussed, among those, the contamination of the stellar halos of galaxies by extended ghost reflections, and the cirrus emission from Galactic dust. The detection and systematic census of fine structures that trace the present and past mass assembly of ETGs is one of the prime goals of the project. We provide specific examples of each type of observed structures -- tidal tails, stellar streams and shells --, and explain how they were identified and classified. We give an overview of the initial results. The detailed statistical analysis will be presented in future papers.

We investigate the early impact of single and binary supernova (SN) explosions on dense gas clouds with three-dimensional, high-resolution, hydrodynamic simulations. The effect of cloud structure, radiative cooling, and ionising radiation from the progenitor stars on the net input of kinetic energy, f_kin=E_kin/E_SN, thermal energy, f_therm=E_therm/E_SN, and gas momentum f_P=P/P_SN to the interstellar medium (ISM) is tested. For clouds with n=100 cm^{-3}, the momentum generating Sedov and pressure-driven snowplough phases are terminated early (~ 0.01 Myr) and radiative cooling limits the coupling to f_therm ~ 0.01, f_kin ~ 0.05, and f_P ~ 9, significantly lower than for the case without cooling. For pre-ionised clouds these numbers are only increased by ~ 50%, independent of the cloud structure. This only suffices to accelerate ~ 5% of the cloud to radial velocities >30km/s. A second SN might further enhance the coupling efficiencies if delayed past the Sedov phase of the first explosion. Such very low coupling efficiencies cast doubts on many galaxy-scale sub-resolution models for supernova feedback, most of which are validated a posteriori by qualitative agreement of galaxy properties with observations. Ionising radiation appears not to significantly enhance the immediate coupling of SNe to the surrounding gas as it drives the ISM into inert dense shells and cold clumps, a process which is unresolved in galaxy scale simulations. Our results support previous conclusions that supernovae might only drive a wind if a significant fraction explodes in low-density environments or if they are supported by processes other than ionising radiation.

In this paper we follow up on our previous detection of nuclear ionized outflows in the most massive (log(M*/Msun) >= 10.9) z~1-3 star-forming galaxies (Forster Schreiber et al.), by increasing the sample size by a factor of six (to 44 galaxies above log(M*/Msun) >= 10.9) from a combination of the SINS/zC-SINF, LUCI, GNIRS, and KMOS^3D spectroscopic surveys. We find a fairly sharp onset of the incidence of broad nuclear emission (FWHM in the Ha, [NII], and [SII] lines ~ 450-5300 km/s), with large [NII]/Ha ratios, above log(M*/Msun) ~ 10.9, with about two thirds of the galaxies in this mass range exhibiting this component. Broad nuclear components near and above the Schechter mass are similarly prevalent above and below the main sequence of star-forming galaxies, and at z~1 and ~2. The line ratios of the nuclear component are fit by excitation from active galactic nuclei (AGN), or by a combination of shocks and photoionization. The incidence of the most massive galaxies with broad nuclear components is at least as large as that of AGNs identified by X-ray, optical, infrared or radio indicators. The mass loading of the nuclear outflows is near unity. Our findings provide compelling evidence for powerful, high-duty cycle, AGN-driven outflows near the Schechter mass, and acting across the peak of cosmic galaxy formation.

We report on empirical trends between the dynamically determined stellar initial mass function (IMF) and stellar population properties for a complete, volume-limited sample of 260 early-type galaxies from the Atlas3D project. We study trends between our dynamically-derived IMF normalisation and absorption line strengths, and interpret these via single stellar population- (SSP-) equivalent ages, abundance ratios (measured as [alpha/Fe]), and total metallicity, [Z/H]. We find that old and alpha-enhanced galaxies tend to have on average heavier (Salpeter-like) mass normalisation of the IMF, but stellar population does not appear to be a good predictor of the IMF, with a large range of normalisation at a given population parameter. As a result, we find weak IMF-[alpha/Fe] and IMF-age correlations, and no significant IMF-[Z/H] correlation. The observed trends appear significantly weaker than those reported in studies that measure the IMF normalisation via low-mass star demographics inferred through stellar spectral analysis.

We use the Atlas3D sample to perform a study of the intrinsic shapes of early-type galaxies, taking advantage of the available combined photometric and kinematic data. Based on our ellipticity measurements from the Sloan Digital Sky Survey Data Release 7, and additional imaging from the Isaac Newton Telescope, we first invert the shape distribution of fast and slow rotators under the assumption of axisymmetry. The so-obtained intrinsic shape distribution for the fast rotators can be described with a Gaussian with a mean flattening of q=0.25 and standard deviation sigma_q = 0.14, and an additional tail towards rounder shapes. The slow rotators are much rounder, and are well described with a Gaussian with mean q = 0.63 and sigma_q =0.09. We then checked that our results were consistent when applying a different and independent method to obtain intrinsic shape distributions, by fitting the observed ellipticity distributions directly using Gaussian parametrisations for the intrinsic axis ratios. Although both fast and slow rotators are identified as early-type galaxies in morphological studies, and in many previous shape studies are therefore grouped together, their shape distributions are significantly different, hinting at different formation scenarios. The intrinsic shape distribution of the fast rotators shows similarities with the spiral galaxy population. Including the observed kinematic misalignment in our intrinsic shape study shows that the fast rotators are predominantly axisymmetric, with only very little room for triaxiality. For the slow rotators though there are very strong indications that they are (mildly) triaxial.

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 study how feedback influences baryon infall onto galaxies using cosmological, zoom-in simulations of haloes with present mass $M_{vir}=6.9\times10^{11} M_{\odot}$ to $1.7\times10^{12} M_{\odot}$. Starting at z=4 from identical initial conditions, implementations of weak and strong stellar feedback produce bulge- and disc-dominated galaxies, respectively. Strong feedback favours disc formation: (1) because conversion of gas into stars is suppressed at early times, as required by abundance matching arguments, resulting in flat star formation histories and higher gas fractions; (2) because 50% of the stars form in situ from recycled disc gas with angular momentum only weakly related to that of the z=0 dark halo; (3) because late-time gas accretion is typically an order of magnitude stronger and has higher specific angular momentum, with recycled gas dominating over primordial infall; (4) because 25-30% of the total accreted gas is ejected entirely before z~1, removing primarily low angular momentum material which enriches the nearby inter-galactic medium. Most recycled gas roughly conserves its angular momentum, but material ejected for long times and to large radii can gain significant angular momentum before re-accretion. These processes lower galaxy formation efficiency in addition to promoting disc formation.

We present the correlations between stellar mass, star formation rate (SFR) and [NII]/Ha flux ratio as indicator of gas-phase metallicity for a sample of 222 galaxies at 0.8 < z < 2.6 and log(M*/Msun)=9.0-11.5 from the LUCI, SINS/zC-SINF and KMOS3D surveys. This sample provides a unique analysis of the mass-metallicity relation (MZR) over an extended redshift range using consistent data analysis techniques and strong-line metallicity indicator. We find a constant slope at the low-mass end of the relation and can fully describe its redshift evolution through the evolution of the characteristic turnover mass where the relation begins to flatten at the asymptotic metallicity. At fixed mass and redshift, our data do not show a correlation between the [NII]/Ha ratio and SFR, which disagrees with the 0.2-0.3dex offset in [NII]/Ha predicted by the "fundamental relation" between stellar mass, SFR and metallicity discussed in recent literature. However, the overall evolution towards lower [NII]/Ha at earlier times does broadly agree with these predictions.

We present a novel implementation of cosmic rays (CR) in the magneto-hydrodynamic code FLASH. CRs are described as separate fluids with different energies. CR advection, energy dependent anisotropic diffusion with respect to the magnetic field and adiabatic losses to follow the evolution of spectra are taken into account. We present a first study of the transport and immediate (~150 kyr) dynamical impact of CRs on the turbulent magnetised interstellar medium around supernova remnants on scales up to 80 pc. CR diffusion quickly leads to an efficient acceleration of low-density gas (mainly perpendicular to the magnetic field) with accelerations up to two orders of magnitude above the thermal values. Peaked (at 1 GeV) CR injection spectra have a stronger impact on the dynamics than power-law spectra. For realistic magnetic field configurations low energy CRs (with smaller diffusion coefficients) distribute anisotropically with large spatial variations of a factor of ten and more. Adiabatic losses can change the local spectra perceptibly but do not have an integral effect on the dynamics at the spatial and temporal scales considered here. We discuss the potential global impact of CRs and anisotropic transport on the dynamical structure of the ISM and also detail the limitations of the model.

We analyze the formation histories of 19 galaxies from cosmological smoothed particle hydrodynamics zoom-in resimulations. We construct mock three-colour images and show that the models reproduce observed trends in the evolution of galaxy colours and morphologies. However, only a small fraction of galaxies contains bars. Many galaxies go through phases of central mass growth by in-situ star formation driven by gas-rich mergers or misaligned gas infall. These events lead to accretion of low-angular momentum gas to the centres and leave imprints on the distributions of z=0 stellar circularities, radii and metallicities as functions of age. Observations of the evolution of structural properties of samples of disc galaxies at z=2.5-0.0 infer continuous mass assembly at all radii. Our simulations can only explain this if there is a significant contribution from mergers or misaligned infall, as expected in a LambdaCDM universe. Quiescent merger histories lead to high kinematic disc fractions and inside-out growth, but show little central growth after the last `destructive' merger at z>1.5. For sufficiently strong feedback, as assumed in our models, a moderate amount of merging does not seem to be a problem for the z=0 disc galaxy population, but may rather be a requirement. The average profiles of simulated disc galaxies agree with observations at z>=1.5. At z<=1, there is too much growth in size and too little growth in central mass, possibly due to the under-abundance of bars. The discrepancies may partly be caused by differences between the star formation histories of the simulations and those assumed for observations.

We report the detection of ubiquitous powerful nuclear outflows in massive (> 10^11 Msun) z~2 star-forming galaxies (SFGs), which are plausibly driven by an Active Galactic Nucleus (AGN). The sample consists of the eight most massive SFGs from our SINS/zC-SINF survey of galaxy kinematics with the imaging spectrometer SINFONI, six of which have sensitive high-resolution adaptive optics (AO) assisted observations. All of the objects are disks hosting a significant stellar bulge. The spectra in their central regions exhibit a broad component in Halpha and forbidden [NII] and [SII] line emission, with typical velocity FWHM ~ 1500 km/s, [NII]/Halpha ratio ~ 0.6, and intrinsic extent of 2 - 3 kpc. These properties are consistent with warm ionized gas outflows associated with Type 2 AGN, the presence of which is confirmed via independent diagnostics in half the galaxies. The data imply a median ionized gas mass outflow rate of ~ 60 Msun/yr and mass loading of ~ 3. At larger radii, a weaker broad component is detected but with lower FWHM ~ 485 km/s and [NII]/Halpha ~ 0.35, characteristic for star formation-driven outflows as found in the lower-mass SINS/zC-SINF galaxies. The high inferred mass outflow rates and frequent occurrence suggest the nuclear outflows efficiently expel gas out of the centers of the galaxies with high duty cycles, and may thus contribute to the process of star formation quenching in massive galaxies. Larger samples at high masses will be crucial to confirm the importance and energetics of the nuclear outflow phenomenon, and its connection to AGN activity and bulge growth.

We present measurements of the star formation rate (SFR) in the early-type galaxies (ETGs) of the ATLAS3D sample, based on Wide-field Infrared Survey Explorer (WISE) 22um and Galaxy Evolution Explorer (GALEX) far-ultraviolet emission. We combine these with gas masses estimated from 12CO and HI data in order to investigate the star formation efficiency (SFE) in a larger sample of ETGs than previously available. We first recalibrate (based on WISE data) the relation between old stellar populations (traced at Ks-band) and 22um luminosity, allowing us to remove the contribution of 22um emission from circumstellar dust. We then go on to investigate the position of ETGs on the Kennicutt-Schmidt (KS) relation. Molecular gas-rich ETGs have comparable star formation surface densities to normal spiral galaxy centres, but they lie systematically offset from the KS relation, having lower star formation efficiencies by a factor of ~2.5 (in agreement with other authors). This effect is driven by galaxies where a substantial fraction of the molecular material is in the rising part of the rotation curve, and shear is high. We show here for the first time that although the number of stars formed per unit gas mass per unit time is lower in ETGs, it seems that the amount of stars formed per free-fall time is approximately constant. The scatter around this dynamical relation still correlates with galaxy properties such as the shape of the potential in the inner regions. This leads us to suggest that dynamical properties (such as shear or the global stability of the gas) may be important second parameters that regulate star formation and cause much of the scatter around star-formation relations.

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.

One quarter of all nearby early-type galaxies (ETGs) outside Virgo host a disc/ring of HI with size from a few to tens of kpc and mass up to ~1e+9 solar masses. Here we investigate whether this HI is related to the presence of a stellar disc within the host making use of the classification of ETGs in fast and slow rotators (FR/SR). We find a large diversity of HI masses and morphologies within both families. Surprisingly, SRs are detected as often, host as much HI and have a similar rate of HI discs/rings as FRs. Accretion of HI is therefore not always linked to the growth of an inner stellar disc. The weak relation between HI and stellar disc is confirmed by their frequent kinematical misalignment in FRs, including cases of polar and counterrotating gas. In SRs the HI is usually polar. This complex picture highlights a diversity of ETG formation histories which may be lost in the relative simplicity of their inner structure and emerges when studying their outer regions. We find that LCDM hydrodynamical simulations have difficulties reproducing the HI properties of ETGs. The gas discs formed in simulations are either too massive or too small depending on the star formation feedback implementation. Kinematical misalignments match the observations only qualitatively. The main point of conflict is that nearly all simulated FRs and a large fraction of all simulated SRs host corotating HI. This establishes the HI properties of ETGs as a novel challenge to simulations.

We present a study of the cold gas contents of the Atlas3D early-type galaxies, in the context of their optical colours, near-UV colours, and H\beta\ absorption line strengths. Early-type (elliptical and lenticular) galaxies are not as gas-poor as previously thought, and at least 40% of local early-type galaxies are now known to contain molecular and/or atomic gas. This cold gas offers the opportunity to study recent galaxy evolution through the processes of cold gas acquisition, consumption (star formation), and removal. Molecular and atomic gas detection rates range from 10% to 34% in red sequence early-type galaxies, depending on how the red sequence is defined, and from 50% to 70% in blue early-type galaxies. Notably, massive red sequence early-type galaxies (stellar masses > 5e10 Msun, derived from dynamical models) are found to have HI masses up to M(HI)/Mstar ~ 0.06 and H_2 masses up to M(H$_2$)/Mstar ~ 0.01. Some 20% of all massive early-type galaxies may have retained atomic and/or molecular gas through their transition to the red sequence. However, kinematic and metallicity signatures of external gas accretion (either from satellite galaxies or the intergalactic medium) are also common, particularly at stellar masses <= 5e10 Msun, where such signatures are found in ~ 50% of H$_2$-rich early-type galaxies. Our data are thus consistent with a scenario in which fast rotator early-type galaxies are quenched former spiral galaxies which have undergone some bulge growth processes, and in addition, some of them also experience cold gas accretion which can initiate a period of modest star formation activity. We discuss implications for the interpretation of colour-magnitude diagrams.

We study the dark and luminous mass distributions, circular velocity curves (CVC), line-of-sight kinematics, and angular momenta for a sample of 42 cosmological zoom simulations of massive galaxies. Using a temporal smoothing technique, we are able to reach large radii. We find that: (i)The dark matter halo density profiles outside a few kpc follow simple power-law models, with flat dark matter CVCs for lower-mass systems, and rising CVCs for high-mass haloes. The projected stellar density distributions at large radii can be fitted by Sersic functions with n>10, larger than for typical ETGs. (ii)The massive systems have nearly flat total CVCs at large radii, while the less massive systems have mildly decreasing CVCs. The slope of the CVC at large radii correlates with v_circ itself. (iii)The dark matter fractions within Re are in the range 15-30% and increase to 40-65% at 5Re. Larger and more massive galaxies have higher dark matter fractions. (iv)The short axes of simulated galaxies and their host dark matter haloes are well aligned and their short-to-long axis ratios are correlated. (v)The stellar vrms(R) profiles are slowly declining, in agreement with planetary nebulae observations in the outer haloes of most ETGs. (vi)The line-of-sight velocity fields v show that rotation properties at small and large radii are correlated. Most radial profiles for the cumulative specific angular momentum parameter lambda(R) are nearly flat or slightly rising, with values in [0.06,0.75] from 2Re to 5Re. (vii)Stellar mass, ellipticity at 5Re, and lambda(5Re) are correlated: the more massive systems have less angular momentum and are rounder, as for observed ETGs. (viii)More massive galaxies with a large fraction of accreted stars have radially anisotropic velocity distributions outside Re. Tangential anisotropy is seen only for galaxies with high fraction of in-situ stars. (Full abstract in PDF)

We present new Spectrographic Areal Unit for Research on Optical Nebulae (SAURON) integral-field spectroscopy and Swift Ultraviolet Optical Telescope (UVOT) observations of molecular outflow host galaxy NGC 1266 that indicate NGC 1266 has experienced a rapid cessation of star formation. Both the SAURON maps of stellar population age and the Swift UVOT observations demonstrate the presence of young ($< 1$ Gyr) stellar populations within the central 1 kpc, while existing Combined Array for Research in Millimeter-wave Astronomy (CARMA) CO(1--0) maps indicate that the sites of current star formation are constrained to the inner few hundred parsecs of the galaxy only. The optical spectrum of NGC 1266 from Moustakas & Kennicutt (2006) reveal a characteristic post-starburst (K+A) stellar population and Davis et al. (2012) confirm that ionized gas emission in the system originate from a shock. Galaxies with K+A spectra and shock-like ionized gas line ratios may comprise an important, overlooked segment of the post-starburst population, containing exactly those objects in which the AGN is actively expelling the star-forming material. While AGN activity is not the likely driver of the post-starburst event that occurred 500 Myr ago, the faint spiral structure seen in the Hubble Space Telescope (HST) Wide-field Camera 3 (WFC3) Y-, J- and H-band imaging seems to point to the possibility of gravitational torques being the culprit. If the molecular gas were driven into the center at the same time as the larger scale galaxy disk underwent quenching, the AGN might be able to sustain the presence of molecular gas for $\gtrsim 1$ Gyr by cyclically injecting turbulent kinetic energy into the dense molecular gas via a radio jet, inhibiting star formation.

We present a detailed two-dimensional stellar dynamical analysis of a sample of 44 cosmological hydrodynamical simulations of individual central galaxies and their satellites. Kinematic maps of the stellar line-of-sight velocity, velocity dispersion, and higher-order Gauss-Hermite moments $h_3$ and $h_4$ are constructed for each central galaxy and for the most massive satellites. The amount of rotation is quantified using the $\lambda_{\mathrm{R}}$-parameter. The velocity, velocity dispersion, $h_3$, and $h_4$ fields of the simulated galaxies show a diversity similar to observed kinematic maps of early-type galaxies in the ATLAS$^{\rm{3D}}$ survey. This includes fast (regular), slow, and misaligned rotation, hot spheroids with embedded cold disk components as well as galaxies with counter-rotating cores or central depressions in the velocity dispersion. We link the present day kinematic properties to the individual cosmological formation histories of the galaxies. In general, major galaxy mergers have a significant influence on the rotation properties resulting in both a spin-down as well as a spin-up of the merger remnant. Lower mass galaxies with significant in-situ formation of stars, or with additional gas-rich major mergers - resulting in a spin-up - in their formation history, form elongated fast rotators with a clear anti-correlation of $h_3$ and $v/\sigma$. An additional formation path for fast rotators includes gas poor major mergers leading to a spin-up of the remnants. This formation path does not result in anti-correlated $h_3$ and $v/\sigma$. The galaxies most consistent with the rare class of non-rotating round early-type galaxies grow by gas-poor minor mergers alone. In general, more massive galaxies have less in-situ star formation since $z \sim 2$, rotate slower and have older stellar populations. (shortened)

We present simulations of the magnetized interstellar medium (ISM) in models of massive star forming (40 Msun / yr) disk galaxies with high gas surface densities (~100 Msun / pc^2) similar to observed star forming high-redshift disks. We assume that type II supernovae deposit 10 per cent of their energy into the ISM as cosmic rays and neglect the additional deposition of thermal energy or momentum. With a typical Galactic diffusion coefficient for CRs (3e28 cm^2 / s) we demonstrate that this process alone can trigger the local formation of a strong low density galactic wind maintaining vertically open field lines. Driven by the additional pressure gradient of the relativistic fluid the wind speed can exceed 1000 km/s, much higher than the escape velocity of the galaxy. The global mass loading, i.e. the ratio of the gas mass leaving the galactic disk in a wind to the star formation rate becomes of order unity once the system has settled into an equilibrium. We conclude that relativistic particles accelerated in supernova remnants alone provide a natural and efficient mechanism to trigger winds similar to observed mass-loaded galactic winds in high-redshift galaxies. These winds also help explaining the low efficiencies for the conversion of gas into stars in galaxies as well as the early enrichment of the intergalactic medium with metals. This mechanism can be at least of similar importance than the traditionally considered momentum feedback from massive stars and thermal and kinetic feedback from supernova explosions.

As part of the SINS/zC-SINF surveys of high-z galaxy kinematics, we derive the radial distributions of H-alpha surface brightness, stellar mass surface density, and dynamical mass at ~2 kpc resolution in 19 z~2 star-forming disks with deep SINFONI AO spectroscopy at the ESO VLT. From these data we infer the radial distribution of the Toomre Q-parameter for these main-sequence star forming galaxies (SFGs), covering almost two decades of stellar mass (10^9.6 to 10^11.5 solar masses). In more than half of our SFGs, the H-alpha distributions cannot be fit by a centrally peaked distribution, such as an exponential, but are better described by a ring, or the combination of a ring and an exponential. At the same time the kinematic data indicate the presence of a mass distribution more centrally concentrated than a single exponential distribution for 5 of the 19 galaxies. The resulting Q-distributions are centrally peaked for all, and significantly exceed unity there for three quarters of the SFGs. The occurrence of H-alpha rings and of large nuclear Q-values is strongly correlated, and is more common for the more massive SFGs. While our sample is small and there remain substantial uncertainties and caveats, our observations are consistent with a scenario in which cloud fragmentation and global star formation are secularly suppressed in gas rich high-z disks from the inside out, as the central stellar mass density of the disks grows.

We investigate the differential effects of metal cooling and galactic stellar winds on the cosmological formation of individual galaxies with three sets of cosmological, hydrodynamical zoom simulations of 45 halos in the mass range 10^11<M_halo<10^13M_sun. Models including both galactic winds and metal cooling (i) suppress early star formation at z>1 and predict reasonable star formation histories, (ii) produce galaxies with high cold gas fractions (30-60 per cent) at high redshift, (iii) significantly reduce the galaxy formation efficiencies for halos (M_halo<10^12M_sun) at all redshifts in agreement with observational and abundance matching constraints, (iv) result in high-redshift galaxies with reduced circular velocities matching the observed Tully-Fisher relation at z~2, and (v) significantly increase the sizes of low-mass galaxies (M_stellar<3x10^10M_sun) at high redshift resulting in a weak size evolution - a trend in agreement with observations. However, the low redshift (z<0.5) star formation rates of massive galaxies are higher than observed (up to ten times). No tested model predicts the observed size evolution for low-mass and high-mass galaxies simultaneously. Due to the delayed onset of star formation in the wind models, the metal enrichment of gas and stars is delayed and agrees well with observational constraints. Metal cooling and stellar winds are both found to increase the ratio of in situ formed to accreted stars - the relative importance of dissipative vs. dissipationless assembly. For halo masses below ~10^12M_sun, this is mainly caused by less stellar accretion and compares well to predictions from semi-analytical models but still differs from abundance matching models. For higher masses, the fraction of in situ stars is over-predicted due to the unrealistically high star formation rates at low redshifts.

We present an update to the multiphase SPH galaxy formation code by Scannapieco et al. We include a more elaborate treatment of the production of metals, cooling rates based on individual element abundances, and a scheme for the turbulent diffusion of metals. Our SN feedback model now transfers energy to the ISM in kinetic and thermal form, and we include a prescription for the effects of radiation pressure from massive young stars on the ISM. We calibrate our new code on the well studied Aquarius haloes and then use it to simulate a sample of 16 galaxies with halo masses between 1x10^11 and 3x10^12 M_sun. In general, the stellar masses of the sample agree well with the stellar mass to halo mass relation inferred from abundance matching techniques for redshifts z=0-4. There is however a tendency to overproduce stars at z>4 and to underproduce them at z<0.5 in the least massive haloes. Overly high SFRs at z<1 for the most massive haloes are likely connected to the lack of AGN feedback in our model. The simulated sample also shows reasonable agreement with observed star formation rates, sizes, gas fractions and gas-phase metallicities at z=0-3. Remaining discrepancies can be connected to deviations from predictions for star formation histories from abundance matching. At z=0, the model galaxies show realistic morphologies, stellar surface density profiles, circular velocity curves and stellar metallicities, but overly flat metallicity gradients. 15 out of 16 of our galaxies contain disk components with kinematic disk fraction ranging between 15 and 65 %. The disk fraction depends on the time of the last destructive merger or misaligned infall event. Considering the remaining shortcomings of our simulations we conclude that even higher kinematic disk fractions may be possible for LambdaCDM haloes with quiet merger histories, such as the Aquarius haloes.

Based on high-resolution, spatially resolved data of 10 z~2 star-forming galaxies from the SINS/zC-SINF survey and LUCI data for 12 additional galaxies, we probe the excitation properties of high-z galaxies and the impact of active galactic nuclei (AGN), shocks and photoionization. We explore how these spatially-resolved line ratios can inform our interpretation of integrated emission line ratios obtained at high redshift. Many of our galaxies fall in the `composite' region of the z~0 [NII]/Ha versus [OIII]/Hb diagnostic (BPT) diagram, between star-forming galaxies and those with AGN. Based on our resolved measurements, we find that some of these galaxies likely host an AGN, while others appear to be affected by the presence of shocks possibly caused by an outflow or from enhanced ionization parameter as compared with HII regions in normal local star-forming galaxies. We find that the Mass-Excitation (MEx) diagnostic, which separates purely star-forming and AGN hosting local galaxies in the [OIII]/Hb versus stellar mass plane, does not properly separate z~2 galaxies classified according to the BPT diagram. However, if we shift the galaxies based on the offset between the local and z~2 mass-metallicity relation (i.e. to the mass they would have at z~0 with the same metallicity), we find better agreement between the MEx and BPT diagnostics. Finally, we find that metallicity calibrations based on [NII]/Ha are more biased by shocks and AGN at high-z than the [OIII]/Hb/[NII]/Ha calibration.

Stellar bars are the most common non-axisymmetric structures in galaxies and their impact on the evolution of disc galaxies at all cosmological times can be significant. Classical theory predicts that stellar discs are stabilized against bar formation if embedded in massive spheroidal dark matter halos. However, dark matter halos have been shown to facilitate the growth of bars through resonant gravitational interaction. Still, it remains unclear why some galaxies are barred and some are not. In this study, we demonstrate that co-rotating (i.e., in the same sense as the disc rotating) dark matter halos with spin parameters in the range of $0 \le \lambda_{\mathrm{dm}} \le 0.07$ - which are a definite prediction of modern cosmological models - promote the formation of bars and boxy bulges and therefore can play an important role in the formation of pseudobulges in a kinematically hot dark matter dominated disc galaxies. We find continuous trends for models with higher halo spins: bars form more rapidly, the forming slow bars are stronger, and the final bars are longer. After 2 Gyrs of evolution, the amplitude of the bar mode in a model with $\lambda_{\mathrm{dm}} = 0.05$ is a factor of ~6 times higher, A_2/A_0 = 0.23, than in the non-rotating halo model. After 5 Gyrs, the bar is ~ 2.5 times longer. The origin of this trend is that more rapidly spinning (co-rotating) halos provide a larger fraction of trailing dark matter particles that lag behind the disc bar and help growing the bar by taking away its angular momentum by resonant interactions. A counter-rotating halo suppresses the formation of a bar in our models. We discuss potential consequences for forming galaxies at high-redshift and present day low mass galaxies which have converted only a small fraction of their baryons into stars.

[Abridged and Edited] We investigate nuclear light profiles in 135 ATLAS3D galaxies for which the Hubble Space Telescope (HST) imaging is available and compare them to the large scale kinematics obtained with the SAURON integral-field spectrograph. Specific angular momentum, lambda_R, correlates with the shape of nuclear light profiles, where cores are typically found in slow rotators and core-less galaxies are fast rotators. Cores are found only in massive galaxies and only in systems with the stellar mass M>8x10^10 Msun. Based on our sample, we, however, see no evidence for a bimodal distribution of nuclear slopes. The best predictor for finding a core is based on the stellar velocity dispersion within an effective radius, sigma_e, and specific angular momentum, where cores are found for lambda_R<0.25 and sigma_e>160 km/s. We estimate that only about 10% of nearby early-type galaxies contain cores. Furthermore, we show that there is a genuine population of fast rotators with cores. We also show that core fast rotators are morphologically, kinematically and dynamically different from core slow rotators. The cores of fast rotators could harbour black holes of similar masses to those in core slow rotators, but typically more massive than those found in core-less fast rotators. Core-less galaxies, and especially core-less fast rotators, are under-luminous in the diffuse X-ray emission, but the presence of a core does not imply high X-ray luminosities. We postulate a possible population of core-less galaxies among slow rotators, which can not be explained as face-on discs, but comprise a genuine sub-population of slow rotators. These galaxies are typically less massive and flatter than core slow rotators, and show evidence for dynamical cold structures and exponential photometric components. We discuss possible processes for formation of cores and their subsequent preservation.

We present a set of hydro-dynamical numerical simulations of the Antennae galaxies in order to understand the origin of the central overlap starburst. Our dynamical model provides a good match to the observed nuclear and overlap star formation, especially when using a range of rather inefficient stellar feedback efficiencies (0.01 < q_EoS < 0.1). In this case a simple conversion of local star formation to molecular hydrogen surface density motivated by observations accounts well for the observed distribution of CO. Using radiative transfer post-processing we model synthetic far-infrared spectral energy distributions (SEDs) and two-dimensional emission maps for direct comparison with Herschel-PACS observations. For a gas-to-dust ratio of 62:1 and the best matching range of stellar feedback efficiencies the synthetic far-infrared SEDs of the central star forming region peak at values of ~65 - 81 Jy at 99 - 116 um, similar to a three-component modified black body fit to infrared observations. Also the spatial distribution of the far-infrared emission at 70 um, 100 um, and 160 um compares well with the observations: >50% (> 35%) of the emission in each band is concentrated in the overlap region while only < 30% (< 15%) is distributed to the combined emission from the two galactic nuclei in the simulations (observations). As a proof of principle we show that parameter variations in the feedback model result in unambiguous changes both in the global and in the spatially resolved observable far-infrared properties of Antennae galaxy models. Our results strengthen the importance of direct, spatially resolved comparative studies of matched galaxy merger simulations as a valuable tool to constrain the fundamental star formation and feedback physics.

Massive galaxies in the distant Universe form stars at much higher rates than today. Although direct resolution of the star forming regions of these galaxies is still a challenge, recent molecular gas observations at the IRAM Plateau de Bure interferometer enable us to study the star formation efficiency on subgalactic scales around redshift z = 1.2. We present a method for obtaining the gas and star formation rate (SFR) surface densities of ensembles of clumps composing galaxies at this redshift, even though the corresponding scales are not resolved. This method is based on identifying these structures in position-velocity diagrams corresponding to slices within the galaxies. We use unique IRAM observations of the CO(3-2) rotational line and DEEP2 spectra of four massive star forming distant galaxies - EGS13003805, EGS13004291, EGS12007881, and EGS13019128 in the AEGIS terminology - to determine the gas and SFR surface densities of the identifiable ensembles of clumps that constitute them. The integrated CO line luminosity is assumed to be directly proportional to the total gas mass, and the SFR is deduced from the [OII] line. We identify the ensembles of clumps with the angular resolution available in both CO and [OII] spectroscopy; i.e., 1-1.5". SFR and gas surface densities are averaged in areas of this size, which is also the thickness of the DEEP2 slits and of the extracted IRAM slices, and we derive a spatially resolved Kennicutt-Schmidt (KS) relation on a scale of ~8 kpc. The data generally indicates an average depletion time of 1.9 Gyr, but with significant variations from point to point within the galaxies.

We study the (M,sigma) and (M,Re) projections of the thin Mass Plane (MP) (M,sigma,Re) which describes the distribution of the galaxy population. The distribution of galaxy properties on the MP is characterized by: (i) a zone of exclusion described by two power-laws joined by a break at M 3e10 Msun. This results in a break in the mean M-sigma relation with Msigma^2.3 and M sigma^4.7 at small and large sigma respectively; (ii) a mass M 2e11 Msun which separates a population dominated by flat fast rotator with disks and spiral galaxies at lower masses, from one dominated by quite round slow rotators at larger masses; (iii) below that mass the distribution of ETGs properties tends to be constant along lines of constant sigma; (iv) it forms a parallel sequence with the distribution of spiral galaxies; (v) at even lower masses, the distribution of fast rotator ETGs and late spirals naturally extends to that of dwarf ETGs (Sph) and dwarf irregulars (Im) respectively. We show via dynamical models that sigma traces the bulge fraction, which drives the observed trends in M/L, Hbeta, colour, IMF and molecular gas fraction. We interpret this as due to a combination of two main effects: (i) an increase of the bulge fraction which increases sigma and greatly enhances the likelihood for a galaxy to have its star formation quenched, and (ii) dry merging, increasing galaxy mass along lines of nearly constant sigma, while leaving the population unchanged. [Abriged]

We report matched resolution, imaging spectroscopy of the CO J=3-2 line (with the IRAM Plateau de Bure millimeter interferometer) and of the H-alpha line (with LUCI at the Large Binocular Telescope)in the massive z=1.53 main-sequence galaxy EGS 13011166, as part of the "Plateau de Bure high-z, blue sequence survey (PHIBSS). We combine these data with HST V-J-J-H-band maps to derive spatially resolved distributions of stellar surface density, star formation rate, molecular gas surface density, optical extinction and gas kinematics. The spatial distribution and kinematics of the ionized and molecular gas are remarkably similar and are well modeled by a turbulent, globally Toomre unstable rotating disk. The stellar surface density distribution is smoother than the clumpy rest-frame UV/optical light distribution, and peaks in an obscured, star forming massive bulge near the dynamical center. The molecular gas surface density and the effective optical screen extinction track each other and are well modeled by a 'mixed' extinction model. The inferred slope of the spatially resolved molecular gas to star formation rate relation depends strongly on the adopted extinction model and can vary from 0.8 to 1.7. For the preferred mixed dust-gas model we find a near linear slope.

We analyze the spectra, spatial distributions and kinematics of Ha, [NII] and [SII] emission in a sample of 42, z~2.2 UV/optically selected star forming galaxies (SFGs) from the SINS & zC-SINF surveys, 35 of which were observed in the adaptive optics mode of SINFONI. This is supplemented by kinematic data from 48 z~1-2.5 galaxies from the literature. We find that the kinematic classification of the high-z SFGs as `dispersion dominated' or `rotation dominated' correlates most strongly with their intrinsic sizes. Smaller galaxies are more likely `dispersion-dominated' for two main reasons: 1) The rotation velocity scales linearly with galaxy size but intrinsic velocity dispersion does not depend on size, and as such, their ratio is systematically lower for smaller galaxies, and 2) Beam smearing strongly decreases large-scale velocity gradients and increases observed dispersion much more for galaxies with sizes at or below the resolution. Dispersion dominated SFGs may thus have intrinsic properties similar to `rotation dominated' SFGs, but are primarily more compact, lower mass, less metal enriched and may have higher gas fractions, plausibly because they represent an earlier evolutionary state.

We study the volume-limited and nearly mass selected (stellar mass M>6*10^9 Msun) Atlas3D sample of 260 early-type galaxies. We construct detailed axisymmetric dynamical models (JAM), which allow for orbital anisotropy, include a dark matter halo, and reproduce in detail both the galaxy images and the high-quality integral-field stellar kinematics. We derive accurate total M/L and dark matter fractions f_DM, within a sphere of radius r=Re. We also measure the stellar M/L and derive a median dark matter fraction f_DM=13%. We find that the thin two-dimensional subset spanned by galaxies in the (M_JAM,sigma_e,R_e) coordinates system, which we call the Mass Plane (MP) has an observed rms scatter of 19% and an intrinsic one of 11%. The MP satisfies the scalar virial relation M_JAM sigma_e^2 R_e within our tight errors. However, the details of how both Re and sigma_e are determined are critical in defining the precise deviation from the virial exponents. We revisit the (M/L)-sigma_e relation, which describes most of the deviations between the MP and the FP. The best-fitting relation is (M/L) sigma_e^0.72 (r-band). It provides an upper limit to any systematic increase of the IMF mass normalization with sigma_e. We study of the link between sigma_e and the galaxies circular velocity V_circ within 1Re (where stars dominate) and find the relation max(V_circ)~1.76*sigma_e, which has an observed scatter of 7%. The accurate parameters described in this paper are used in the companion Paper XX of this series to explore the variation of global galaxy properties, including the IMF, on the projections of the MP. [Abridged]

Dynamical modeling and strong lensing data indicate that the total density profiles of early-type galaxies are close to isothermal, i.e., rho_tot ~ r^gamma with gamma approx -2. To understand the origin of this universal slope we study a set of simulated spheroids formed in isolated binary mergers as well as the formation within the cosmological framework. The total stellar plus dark matter density profiles can always be described by a power law with an index of gamma approx -2.1 with a tendency toward steeper slopes for more compact, lower-mass ellipticals. In the binary mergers the amount of gas involved in the merger determines the precise steepness of the slope. This agrees with results from the cosmological simulations where ellipticals with steeper slopes have a higher fraction of stars formed in situ. Each gas-poor merger event evolves the slope toward gamma ~ -2, once this slope is reached further merger events do not change it anymore. All our ellipticals have flat intrinsic combined stellar and dark matter velocity dispersion profiles. We conclude that flat velocity dispersion profiles and total density distributions with a slope of gamma ~ -2 for the combined system of stars and dark matter act as a natural attractor. The variety of complex formation histories as present in cosmological simulations, including major as well as minor merger events, is essential to generate the full range of observed density slopes seen for present-day elliptical galaxies.

We study the formation and evolution of Brightest Cluster Galaxies starting from a $z=2$ population of quiescent ellipticals and following them to $z=0$. To this end, we use a suite of nine high-resolution dark matter-only simulations of galaxy clusters in a $\Lambda$CDM universe. We develop a scheme in which simulation particles are weighted to generate realistic and dynamically stable stellar density profiles at $z=2$. Our initial conditions assign a stellar mass to every identified dark halo as expected from abundance matching; assuming there exists a one-to-one relation between the visible properties of galaxies and their host haloes. We set the sizes of the luminous components according to the observed relations for $z\sim2$ massive quiescent galaxies. We study the evolution of the mass-size relation, the fate of satellite galaxies and the mass aggregation of the cluster central. From $z=2$, these galaxies grow on average in size by a factor 5 to 10 of and in mass by 2 to 3. The stellar mass growth rate of the simulated BCGs in our sample is of 1.9 in the range $0.3<z<1.0$ consistent with observations, and of 1.5 in the range $0.0<z<0.3$. Furthermore the satellite galaxies evolve to the present day mass-size relation by $z=0$. Assuming passively evolving stellar populations, we present surface brightness profiles for our cluster centrals which resemble those observed for the cDs in similar mass clusters both at $z=0$ and at $z=1$. This demonstrates that the $\Lambda$CDM cosmology does indeed predict minor and major mergers to occur in galaxy clusters with the frequency and mass ratio distribution required to explain the observed growth in size of passive galaxies since $z=2$. Our experiment shows that Brightest Cluster Galaxies can form through dissipationless mergers of quiescent massive $z=2$ galaxies, without substantial additional star formation.

For early-type galaxies, the ability to sustain a corona of hot, X-ray emitting gas could have played a key role in quenching their star-formation history. Yet, it is still unclear what drives the precise amount of hot gas around these galaxies. By combining photometric and spectroscopic measurements for the early-type galaxies observed during the Atlas3D integral-field survey with measurements of their X-ray luminosity based on X-ray data of both low and high spatial resolution we conclude that the hot-gas content of early-type galaxies can depend on their dynamical structure. Specifically, whereas slow rotators generally have X-ray halos with luminosity L_X,gas and temperature T values that are in line with what is expected if the hot-gas emission is sustained by the thermalisaton of the kinetic energy carried by the stellar-mass loss material, fast rotators tend to display L_X,gas values that fall consistently below the prediction of this model, with similar T values that do not scale with the stellar kinetic energy as observed in the case of slow rotators. Considering that fast rotators are likely to be intrinsically flatter than slow rotators, and that the few L_X,gas-deficient slow rotators also happen to be relatively flat, the observed L_X,gas deficiency in these objects would support the hypothesis whereby flatter galaxies have a harder time in retaining their hot gas. We discuss the implications that a different hot-gas content could have on the fate of both acquired and internally-produced gaseous material, considering in particular how the L_X,gas deficiency of fast rotators would make them more capable to recycle the stellar-mass loss material into new stars than slow rotators. This is consistent with the finding that molecular gas and young stars are detected only in fast rotators in the Atlas3D sample, and that fast rotators tend to dustier than slow rotators. [Abridged]

[Abridged] We present a detailed study of the physical properties of the molecular gas in a sample of 18 molecular gas-rich early-type galaxies (ETGs) from the ATLAS$ 3D sample. Our goal is to better understand the star formation processes occurring in those galaxies, starting here with the dense star-forming gas. We use existing integrated $^{12}$CO(1-0, 2-1), $^{13}$CO(1-0, 2-1), HCN(1-0) and HCO$^{+}$(1-0) observations and present new $^{12}$CO(3-2) single-dish data. From these, we derive for the first time the average kinetic temperature, H$_{2}$ volume density and column density of the emitting gas, this using a non-LTE theoretical model. Since the CO lines trace different physical conditions than of those the HCN and HCO$^{+}$ lines, the two sets of lines are treated separately. We also compare for the first time the predicted CO spectral line energy distributions (SLEDs) and gas properties of our molecular gas-rich ETGs with those of a sample of nearby well-studied disc galaxies. The gas excitation conditions in 13 of our 18 ETGs appear analogous to those in the centre of the Milky Way. Such results have never been obtained before for ETGs and open a new window to explore further star-formation processes in the Universe. The conclusions drawn should nevertheless be considered carefully, as they are based on a limited number of observations and on a simple model. In the near future, with higher CO transition observations, it should be possible to better identify the various gas components present in ETGs, as well as more precisely determine their associated physical conditions. To achieve these goals, we show here from our theoretical study, that mid-J CO lines (such as the $^{12}$CO(6-5) line) are particularly useful.

We study the global efficiency of star formation in high resolution hydrodynamical simulations of gas discs embedded in isolated early-type and spiral galaxies. Despite using a universal local law to form stars in the simulations, we find that the early-type galaxies are offset from the spirals on the large-scale Kennicutt relation, and form stars 2 to 5 times less efficiently. This offset is in agreement with previous results on morphological quenching: gas discs are more stable against star formation when embedded in early-type galaxies due to the lower disc self-gravity and increased shear. As a result, these gas discs do not fragment into dense clumps and do not reach as high densities as in the spiral galaxies. Even if some molecular gas is present, the fraction of very dense gas (above 10^4 cm-3) is significantly reduced, which explains the overall lower star formation efficiency. We also analyse a sample of local early-type and spiral galaxies, measuring their CO and HI surface densities and their star formation rates as determined by their non-stellar 8um emission. As predicted by the simulations, we find that the early-type galaxies are offset from the Kennicutt relation compared to the spirals, with a twice lower efficiency. Finally, we validate our approach by performing a direct comparison between models and observations. We run a simulation designed to mimic the stellar and gaseous properties of NGC524, a lenticular galaxy, and find a gas disc structure and global star formation rate in good agreement with the observations. Morphological quenching thus seems to be a robust mechanism, and is also consistent with other observations of a reduced star formation efficiency in early-type galaxies in the COLD GASS survey. This lower efficiency of star formation is not enough to explain the formation of the whole Red Sequence, but can contribute to the reddening of some galaxies.

There is observational evidence for inside-out growth of elliptical galaxies since $z \gtrsim 2-3$, which is not driven by in-situ star formation. Many systems at high redshift have small sizes $\sim 1kpc$ and surface brightness profiles with low Sersic indices n. The most likely descendants have, on average, grown by a factor of two in mass and a factor of four in size, indicating $r \propto M^{\alpha}$ with $\alpha \gtrsim 2$. They also have surface brightness profiles with $n \gtrsim 5$. This evolution can be qualitatively explained on the basis of two assumptions: compact ellipticals predominantly grow by collisionless minor or intermediate 'dry' mergers, and they are embedded in massive dark matter halos. We draw these conclusions from idealized collisionless mergers spheroidal galaxies - with and without dark matter - with mass ratios of 1:1, 1:5, and 1:10. The sizes evolve as $r \propto M^{\alpha}$ with $\alpha < 2$ for mass-ratios of 1:1. For minor mergers of galaxies embedded in dark matter halos, the sizes grow significantly faster and the profile shapes change more rapidly. Mergers with moderate mass-ratios of 1:5 give $\alpha \sim 2.3$ and a final Sersic index of $n = 9.5$ after doubling the stellar mass. This is accompanied by a significant increase of the dark matter fraction within the stellar half-mass radius, driven by the strong size increase probing larger, dark matter dominated regions. Only a few intermediate mass-ratio mergers of galaxies embedded in massive dark matter halos can result in the observed concurrent inside-out growth and the rapid evolution in profile shapes. Apart from negative stellar metallicity gradients such a 'minor' merger scenario also predicts significantly lower dark matter fractions for $z \sim 2$ compact quiescent galaxies and their rare present day analogues (abbreviated).

The discovery of a population of massive, compact and quiescent early-type galaxies has changed the view on plausible formation scenarios for the present day population of elliptical galaxies. Traditionally assumed formation histories dominated by 'single events' like early collapse or major mergers appear to be incomplete and have to be embedded in the context of hierarchical cosmological models with continuous gas accretion and the merging of small stellar systems (minor mergers). Once these processes are consistently taken into account the hierarchical models favor a two-phase assembly process and are in much better shape to capture the observed trends. We review some aspects of recent progress in the field.

Based on SINFONI Ha, [NII] and [SII] AO data of 30 z \sim 2 star-forming galaxies (SFGs) from the SINS and zcSINF surveys, we find a strong correlation of the Ha broad flux fraction with the star formation surface density of the galaxy, with an apparent threshold for strong outflows occurring at 1 Msun yr^-1 kpc^-2. Above this threshold, we find that SFGs with logm_\ast>10 have similar or perhaps greater wind mass loading factors (eta = Mdotout/SFR) and faster outflow velocities than lower mass SFGs. This trend suggests that the majority of outflowing gas at z \sim 2 may derive from high-mass SFGs, and that the z \sim 2 mass-metallicity relation is driven more by dilution of enriched gas in the galaxy gas reservoir than by the efficiency of outflows. The mass loading factor is also correlated with the SFR and inclination, such that more star-forming and face-on galaxies launch more powerful outflows. For galaxies that have evidence for strong outflows, we find that the broad emission is spatially extended to at least the half-light radius (\sim a few kpc). We propose that the observed threshold for strong outflows and the observed mass loading of these winds can be explained by a simple model wherein break-out of winds is governed by pressure balance in the disk. Using the ratio of the [SII] doublet in a broad and narrow component, we find that outflowing gas has a density of \sim10-100 cm^-3, significantly less than that of the star forming gas (600 cm^-3).

We present PHIBSS, the IRAM Plateau de Bure high-z blue sequence CO 3-2 survey of the molecular gas properties in normal star forming galaxies (SFGs) near the cosmic star formation peak. PHIBSS provides 52 CO detections in two redshift slices at z~1.2 and 2.2, with log(M*(M_solar))>10.4 and log(SFR(M_solar/yr))>1.5. Including a correction for the incomplete coverage of the M*-SFR plane, we infer average gas fractions of ~0.33 at z~1.2 and ~0.47 at z~2.2. Gas fractions drop with stellar mass, in agreement with cosmological simulations including strong star formation feedback. Most of the z~1-3 SFGs are rotationally supported turbulent disks. The sizes of CO and UV/optical emission are comparable. The molecular gas - star formation relation for the z=1-3 SFGs is near-linear, with a ~0.7 Gyrs gas depletion timescale; changes in depletion time are only a secondary effect. Since this timescale is much less than the Hubble time in all SFGs between z~0 and 2, fresh gas must be supplied with a fairly high duty cycle over several billion years. At given z and M*, gas fractions correlate strongly with the specific star formation rate. The variation of specific star formation rate between z~0 and 3 is mainly controlled by the fraction of baryonic mass that resides in cold gas.

We explore the connection between the local escape velocity, V_esc, and the stellar population properties in the ATLAS3D survey, a complete, volume-limited sample of nearby early-type galaxies. We make use of ugriz photometry to construct Multi-Gaussian Expansion models of the surface brightnesses of our galaxies. We are able to fit the full range of surface brightness profiles found in our sample, and in addition we reproduce the results of state-of-the-art photometry in the literature with residuals of 0.04 mags. We utilise these photometric models and SAURON integral-field spectroscopy, combined with Jeans dynamical modelling, to determine the local V_esc derived from the surface brightness. We find that the local V_esc is tightly correlated with the Mgb and Fe5015 linestrengths and optical colours, and anti-correlated with the Hbeta linestrength. In the case of the Mgb and Colour - V_esc relations we find that the relation within individual galaxies follows the global relation between different galaxies. We intentionally ignored any uncertain contribution due to dark matter since we are seeking an empirical description of stellar population gradients in early-type galaxies that is ideal for quantitative comparison with model predictions. We identify a population of galaxies that occur only at low V_esc that exhibit negative gradients in the Mgb - and Colour -V_esc relations. These galaxies typically have young central stellar populations and contain significant amounts of molecular gas and dust. Combining these results with N-body simulations of binary mergers we use the Mgb-V_esc relation to constrain the possible number of dry mergers experienced by the local early-type galaxy population - a typical massive ETG can have experienced only ~1.5 major mergers before becoming a significant outlier in the Mgb-V_esc relation. [Abridged]

We detail the rich molecular story of NGC 1266, its serendipitous discovery within the ATLAS3D survey (Cappellari et al. 2011) and how it plays host to an AGN-driven molecular outflow, potentially quenching all of its star formation (SF) within the next 100 Myr. While major mergers appear to play a role in instigating outflows in other systems, deep imaging of NGC 1266 as well as stellar kinematic observations from SAURON, have failed to provide evidence that NGC 1266 has recently been involved in a major interaction. The molecular gas and the instantaneous SF tracers indicate that the current sites of star formation are located in a hypercompact disk within 200 pc of the nucleus (Fig. 1; SF rate ~ 2 Msuns/yr). On the other hand, tracers of recent star formation, such as the H{\beta} absorption map from SAURON and stellar population analysis show that the young stars are distributed throughout a larger area of the galaxy than current star formation. As the AGN at the center of NGC 1266 continues to drive cold gas out of the galaxy, we expect star formation rates to decline as the star formation is ultimately quenched. Thus, NGC 1266 is in the midst of a key portion of its evolution and continued studies of this unique galaxy may help improve our understanding of how galaxies transition from the blue to the red sequence (Alatalo et al. 2011).

NGC 1266 is a nearby field galaxy observed as part of the ATLAS3D survey (Cappellari et al. 2011). NGC 1266 has been shown to host a compact (< 200 pc) molecular disk and a mass-loaded molecular outflow driven by the AGN (Alatalo et al. 2011). Very Long Basline Array (VLBA) observations at 1.65 GHz revealed a compact (diameter < 1.2 pc), high bright- ness temperature continuum source most consistent with a low-level AGN origin. The VLBA continuum source is positioned at the center of the molecular disk and may be responsible for the expulsion of molecular gas in NGC 1266. Thus, the candidate AGN-driven molecular outflow in NGC 1266 supports the picture in which AGNs do play a significant role in the quenching of star formation and ultimately the evolution of the red sequence of galaxies.

[Abridged] We analyse the morphological structures in galaxies of the ATLAS3D sample by fitting a single Sersic profile and decomposing all non-barred objects (180 of 260 objects) in two components parameterised by an exponential and a general Sersic function. The aim of this analysis is to look for signatures of discs in light distributions of nearby early-type galaxies and compare them to kinematic properties. Using Sersic index from single component fits for a distinction between slow and fast rotators, or even late- and early-type galaxies, is not recommended. Assuming that objects with n>3 are slow rotators (or ellipticals), there is only a 22 per cent probability to correctly classify objects as slow rotators (or 37 per cent of previously classified as ellipticals). We show that exponential sub-components, as well as light profiles fitted with only a single component of a low Sersic index, can be linked with the kinematic evidence for discs in early-type galaxies. The median disk-to-total light ratio for fast and slow rotators is 0.41 and 0.0, respectively. Similarly, the median Sersic indices of the bulge (general Sersic component) are 1.7 and 4.8 for fast and slow rotators, respectively. Overall, discs or disc-like structures, are present in 83 per cent of early-type galaxies which do not have bars, and they show a full range of disk-to-total light ratios. Discs in early-type galaxies contribute with about 40 per cent to the total mass of the analysed (non-barred) objects. The decomposition into discs and bulges can be used as a rough approximation for the separation between fast and slow rotators, but it is not a substitute, as there is only a 59 per cent probability to correctly recognise slow rotators. Kinematics (i.e. projected angular momentum) remains the best approach to mitigate the influence of the inclination effects.

We use interferometric CO observations to compare the extent, surface brightness profiles and kinematics of the molecular gas in CO-rich Atlas3D early-type galaxies (ETGs) and spiral galaxies. We find that the molecular gas extent is smaller in absolute terms in ETGs than in late-type galaxies, but that the size distributions are similar once scaled by the galaxies optical/stellar characteristic scale-lengths. Virgo cluster ETGs have less extended molecular gas reservoirs than field counterparts. Approximately half of ETGs have molecular gas surface brightness profiles that follow the stellar light profile. These systems often have relaxed gas out to large radii, suggesting they are unlikely to have had recent merger/accretion events. A third of the sample galaxies show molecular gas surface brightness profiles that fall off slower than the light, and sometimes show a truncation. We suggest that ram pressure stripping and/or the presence of hot gas has compressed/truncated the gas in these systems. The remaining galaxies have rings, or composite profiles, that we argue can be caused by the effects of bars. We investigated the kinematics of the molecular gas using position-velocity diagrams, and compared the observed kinematics with dynamical model predictions, and the observed stellar and ionised gas velocities. We confirm that the molecular gas reaches beyond the turnover of the circular velocity curve in ~70% of our CO-rich Atlas3D ETGs. In general we find that in most galaxies the molecular gas is relaxed and dynamically cold. Molecular gas is a better direct tracer of the circular velocity than the ionised gas, justifying its use as a kinematic tracer for Tully-Fisher and similar analyses. (abridged)

We present the Combined Array for Research in Millimeter Astronomy (CARMA) ATLAS3D molecular gas imaging survey, a systematic study of the distribution and kinematics of molecular gas in CO-rich early-type galaxies. Our full sample of 40 galaxies (30 newly mapped and 10 taken from the literature) is complete to a 12CO(1-0) integrated flux of 18.5 Jy km/s, and it represents the largest, best-studied sample of its type to date. A comparison of the CO distribution of each galaxy to the g-r color image (representing dust) shows that the molecular gas and dust distributions are in good agreement and trace the same underlying interstellar medium. The galaxies exhibit a variety of CO morphologies, including discs (50%), rings (15%), bars+rings (10%), spiral arms (5%), and mildly (12.5%) and strongly (7.5%) disrupted morphologies. There appear to be weak trends between galaxy mass and CO morphology, whereby the most massive galaxies in the sample tend to have molecular gas in a disc morphology. We derive a lower limit to the total accreted molecular gas mass across the sample of 2.48x10^10 Msuns, or approximately 8.3x10^8 Msuns per minor merger within the sample, consistent with minor merger stellar mass ratios.

Observational studies have revealed a "downsizing" trend in black hole (BH) growth: the number densities of luminous AGN peak at higher redshifts than those of faint AGN. This would seem to imply that massive black holes formed before low mass black holes, in apparent contradiction to hierarchical clustering scenarios. We investigate whether this observed "downsizing" in BH growth is reproduced in a semi-analytic model for the formation and evolution of galaxies and black holes, set within the hierarchical paradigm for structure formation (Somerville et al. 2008; S08). In this model, black holes evolve from light seeds (\sim100M\odot) and their growth is merger-driven. The original S08 model (baseline model) reproduces the number density of AGN at intermediate redshifts and luminosities, but underproduces luminous AGN at very high redshift (z > 3) and overproduces them at low redshift (z < 1). In addition, the baseline model underproduces low-luminosity AGN at low redshift (z < 1). To solve these problems we consider several modifications to the physical processes in the model: (1) a 'heavy' black hole seeding scenario (2) a sub-Eddington accretion rate ceiling that depends on the cold gas fraction, and (3) an additional black hole accretion mode due to disk instabilities. With these three modifications, the models can explain the observed downsizing, successfully reproduce the bolometric AGN luminosity function and simultaneously reproduce galaxy and black hole properties in the local Universe. We also perform a comparison with the observed soft and hard X-ray luminosity functions of AGN, including an empirical correction for torus-level obscuration, and reach similar conclusions. Our best-fit model suggests a scenario in which disk instabilities are the main driver for moderately luminous Seyfert galaxies at low redshift, while major mergers are the main trigger for luminous AGN.

We report the discovery of a giant HI tail in the intra-group medium of HCG 44 as part of the Atlas3D survey. The tail is ~300 kpc long in projection and contains ~5x10^8 M_sun of HI. We detect no diffuse stellar light at the location of the tail down to ~28.5 mag/arcsec^2 in g band. We speculate that the tail might have formed as gas was stripped from the outer regions of NGC 3187 (a member of HCG 44) by the group tidal field. In this case, a simple model indicates that about 1/3 of the galaxy's HI was stripped during a time interval of <1 Gyr. Alternatively, the tail may be the remnant of an interaction between HCG 44 and NGC 3162, a spiral galaxy now ~650 kpc away from the group. Regardless of the precise formation mechanism, the detected HI tail shows for the first time direct evidence of gas stripping in HCG 44. It also highlights that deep HI observations over a large field are needed to gather a complete census of this kind of events in the local Universe.

We use the SAURON and GMOS integral field spectrographs to observe the active galactic nucleus (AGN) powered outflow in NGC 1266. This unusual galaxy is relatively nearby (D=30 Mpc), allowing us to investigate the process of AGN feedback in action. We present maps of the kinematics and line strengths of the ionised gas emission lines Halpha, Hbeta, [OIII], [OI], [NII] and [SII], and report on the detection of Sodium D absorption. We use these tracers to explore the structure of the source, derive the ionised and atomic gas kinematics and investigate the gas excitation and physical conditions. NGC 1266 contains two ionised gas components along most lines of sight, tracing the ongoing outflow and a component closer to the galaxy systemic, the origin of which is unclear. This gas appears to be disturbed by a nascent AGN jet. We confirm that the outflow in NGC 1266 is truly multiphase, containing radio plasma, atomic, molecular and ionised gas and X-ray emitting plasma. The outflow has velocities up to \pm900 km/s away from the systemic velocity, and is very likely to be removing significant amounts of cold gas from the galaxy. The LINER-like line-emission in NGC 1266 is extended, and likely arises from fast shocks caused by the interaction of the radio jet with the ISM. These shocks have velocities of up to 800 km/s, which match well with the observed velocity of the outflow. Sodium D equivalent width profiles are used to set constraints on the size and orientation of the outflow. The ionised gas morphology correlates with the nascent radio jets observed in 1.4 GHz and 5 GHz continuum emission, supporting the suggestion that an AGN jet is providing the energy required to drive the outflow.

We study the growth of black holes (BHs) in galaxies using three-dimensional smoothed particle hydrodynamic simulations with new implementations of the momentum mechanical feedback, and restriction of accreted elements to those that are gravitationally bound to the BH. We also include the feedback from the X-ray radiation emitted by the BH, which heats the surrounding gas in the host galaxies, and adds radial momentum to the fluid. We perform simulations of isolated galaxies and merging galaxies and test various feedback models with the new treatment of the Bondi radius criterion. We find that overall the BH growth is similar to what has been obtained by earlier workers using the Springel, Di Matteo, & Hernquist algorithms. However, the outflowing wind velocities and mechanical energy emitted by winds are considerably higher (v_w ~ 1000-3000 km/s) compared to the standard thermal feedback model (v_w ~ 50-100 km/s). While the thermal feedback model emits only 0.1 % of BH released energy in winds, the momentum feedback model emits more than 30 % of the total energy released by the BH in winds. In the momentum feedback model, the degree of fluctuation in both radiant and wind output is considerably larger than in the standard treatments. We check that the new model of the BH mass accretion agrees with analytic results for the standard Bondi problem.

We revisit collisionless major and minor mergers of spheroidal galaxies in the context of the size evolution of elliptical galaxies. The simulations are performed as a series of mergers with mass-ratios of 1:1 and 1:10 for models representing pure bulges as well as bulges embedded in dark matter halos. For major and minor mergers, respectively, we identify and analyze two different processes, violent relaxation and stripping, leading to size evolution and a change of the dark matter fraction within the observable effective radius. Violent relaxation - which is the dominant mixing process for major mergers but less important for minor mergers - scatters relatively more dark matter particles than bulge particles to small. Stripping in minor mergers assembles stellar satellite particles at large radii in halo dominated regions of the massive host. This strongly increases the size of the bulge into regions with higher dark matter fractions leaving the inner host structure almost unchanged. A factor of two mass increase by minor mergers increases the dark matter fraction by 20 per cent. We present analytic corrections to simple one-component virial estimates for the evolution of the gravitational radii. If such a two-component system grows by minor mergers alone its size growth, $r_{\mathrm{e}} \propto M^\alpha$, reaches values of $\alpha \approx 2.4$, significantly exceeding the simple theoretical limit of $\alpha = 2$. For major mergers the sizes grow with $\alpha \lesssim 1$. Our results indicate that minor mergers of galaxies embedded in massive dark matter halos provide a potential mechanism for explaining the rapid size growth and the build-up of massive elliptical systems predicting significant dark matter fractions and radially biased velocity dispersions at large radii (abbreviated)

We present a sample of nine high resolution cosmological simulations in the mass range of M_vir=7x10^11-4x10^12 M_sun starting from LambdaCDM initial conditions. Our simulations include primordial radiative cooling, photoionization, star formation, supernova II feedback, but exclude supernova driven winds and AGN feedback. The simulated galaxies assemble in two phases, with the initial growth dominated by compact (r<r_eff) in situ star formation fueled by cold, low entropy gas streams resulting in a very similar mean assembly redshift of z_{f,ins}~2.5 for the in situ stellar component in all galaxies. The late growth is dominated by accretion of old stars formed in subunits outside the main galaxy (r>r_eff) resulting in an assembly redshift of z_{f,acc}~0.5-1.5 with much larger scatter. We find a positive correlation between the fraction of accreted stars and the final mass of our galaxies. We show that gravitational feedback strongly suppresses late star formation in massive galaxies contributing to the observed galaxy color bimodality. The accretion of stellar material is also responsible for the observed size growth of early-type galaxies. In addition, we find that the dark matter fractions within the stellar half-mass radii continuously increase towards lower redshift from about f_DM~0.05 at z~3 to f_DM~0.1-0.3 at z=0. Furthermore, the logarithmic slope of the total density profile is nearly isothermal at the present-day (gamma'~1.9-2.2). Finally, the input of gravitational heating lowers the central dark matter densities in the galaxies, with the effect being smaller compared to simulations without supernova feedback.

We investigate the accuracy in the recovery of the stellar dynamics of barred galaxies when using axisymmetric dynamical models. We do this by trying to recover the mass-to-light ratio (M/L) and the anisotropy of realistic galaxy simulations using the Jeans Anisotropic Multi-Gaussian Expansion (JAM) method. However, given that the biases we find are mostly due to an application of an axisymmetric modeling algorithm to a non-axisymmetric system and in particular to inaccuracies in the de-projected mass model, our results are relevant for general axisymmetric modelling methods. We run N-body collisionless simulations to build a library with various luminosity distribution, constructed to mimic real individual galaxies, with realistic anisotropy. The final result of our evolved library of simulations contains both barred and unbarred galaxies. The JAM method assumes an axisymmetric mass distribution, and we adopt a spatially constant M/L and anisotropy beta_z=1-sigma_z^2/sigma_R^2 distributions. The models are fitted to two-dimensional maps of the second velocity moments V_rms=sqrt(V^2+sigma^2) of the simulations for various viewing angles (position angle of the bar and inclination of the galaxy). We find that the inclination is generally well recovered by the JAM models, for both barred and unbarred simulations. For unbarred simulations the M/L is also accurately recovered, with negligible median bias and with a maximum one of just Delta(M/L)<1.5% when the galaxy is not too close to face on. At very low inclinations (i<30 deg) the M/L can be significantly overestimated (9% in our tests, but errors can be larger for very face-on views). For barred simulations the M/L is on average (when PA=45 deg) essentially unbiased, but we measure an over/under estimation of up to Delta(M/L)=15% in our tests. The sign of the M/L bias depends on the position angle of the bar as expected. [Abridged]

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 have obtained high-resolution data of the z 2 ring-like, clumpy star-forming galaxy (SFG) ZC406690 using the VLT/SINFONI with AO (in K-band) and in seeing-limited mode (in H- and J-band). Our data includes all of the main strong optical emission lines: [OII], [OIII], Ha, Hb, [NII] and [SII]. We find broad, blueshifted Ha and [OIII] emission line wings in the spectra of the galaxy's massive, star-forming clumps (sigma \sim 85 km s^-1) and even broader wings (up to 70% of the total Ha flux, with sigma \sim 290 km s^-1) in regions spatially offset from the clumps by \sim 2 kpc. The broad emission likely originates from large-scale outflows with mass outflow rates from individual clumps that are 1-8x the SFR of the clumps. Based on emission line ratio diagnostics ([NII]/Ha and [SII]/Ha) and photoionization and shock models, we find that the emission from the clumps is due to a combination of photoionization from the star-forming regions and shocks generated in the outflowing component, with 5-30% of the emission deriving from shocks. In terms of the ionization parameter (6x10^7-10^8 cm/s, based on both the SFR and the O32 ratio), density (local electron densities of 300-1800 cm^-3 in and around the clumps, and ionized gas column densities of 1200-8000 Msol/pc^2), and SFR (10-40 Msol/yr), these clumps more closely resemble nuclear starburst regions of local ULIRGs and dwarf irregulars than HII regions in local galaxies. However, the star-forming clumps are not located in the nucleus as in local starburst galaxies but instead are situated in a ring several kpc from the center of their high-redshift host galaxy, and have an overall disk-like morphology. The two brightest clumps are quite different in terms of their internal properties, energetics and relative ages, and thus we are given a glimpse at two different stages in the formation and evolution of rapidly star-forming giant clumps at high-z.

The accuracy and robustness of a simple method to estimate the total mass profile of a galaxy is tested using a sample of 65 cosmological zoom-simulations of individual galaxies. The method only requires information on the optical surface brightness and the projected velocity dispersion profiles and therefore can be applied even in case of poor observational data. In the simulated sample massive galaxies ($\sigma \simeq 200-400$ $\kms$) at redshift $z=0$ have almost isothermal rotation curves for broad range of radii (RMS $\simeq 5%$ for the circular speed deviations from a constant value over $0.5R_{\rm eff} < r < 3R_{\rm eff}$). For such galaxies the method recovers the unbiased value of the circular speed. The sample averaged deviation from the true circular speed is less than $\sim 1%$ with the scatter of $\simeq 5-8%$ (RMS) up to $R \simeq 5R_{\rm eff}$. Circular speed estimates of massive non-rotating simulated galaxies at higher redshifts ($z=1$ and $z=2$) are also almost unbiased and with the same scatter. For the least massive galaxies in the sample ($\sigma < 150$ $\kms$) at $z=0$ the RMS deviation is $\simeq 7-9%$ and the mean deviation is biased low by about $1-2%$. We also derive the circular velocity profile from the hydrostatic equilibrium (HE) equation for hot gas in the simulated galaxies. The accuracy of this estimate is about RMS $\simeq 4-5%$ for massive objects ($M > 6.5\times 10^{12} M_\odot$) and the HE estimate is biased low by $\simeq 3-4%$, which can be traced to the presence of gas motions. This implies that the simple mass estimate can be used to determine the mass of observed massive elliptical galaxies to an accuracy of $5-8 %$ and can be very useful for galaxy surveys.

Much of our knowledge of galaxies comes from analysing the radiation emitted by their stars. It depends on the stellar initial mass function (IMF) describing the distribution of stellar masses when the population formed. Consequently knowledge of the IMF is critical to virtually every aspect of galaxy evolution. More than half a century after the first IMF determination, no consensus has emerged on whether it is universal in different galaxies. Previous studies indicated that the IMF and the dark matter fraction in galaxy centres cannot be both universal, but they could not break the degeneracy between the two effects. Only recently indications were found that massive elliptical galaxies may not have the same IMF as our Milky Way. Here we report unambiguous evidence for a strong systematic variation of the IMF in early-type galaxies as a function of their stellar mass-to-light ratio, producing differences up to a factor of three in mass. This was inferred from detailed dynamical models of the two-dimensional stellar kinematics for the large Atlas3D representative sample of nearby early-type galaxies spanning two orders of magnitude in stellar mass. Our finding indicates that the IMF depends intimately on a galaxy's formation history.

We study the evolution of the stellar and dark matter components in a galaxy cluster of $10^{15} \, \rm{M_{\odot}}$ from $z=3$ to the present epoch using the high-resolution collisionless simulations of Ruszkowski & Springel (2009). At $z=3$ the dominant progenitor halos were populated with spherical model galaxies with and without accounting for adiabatic contraction. We apply a weighting scheme which allows us to change the relative amount of dark and stellar material assigned to each simulation particle in order to produce luminous properties which agree better with abundance matching arguments and observed bulge sizes at $z=3$. This permits the study of the effect of initial compactness on the evolution of the mass-size relation. We find that for more compact initial stellar distributions the size of the final Brightest Cluster Galaxy grows with mass according to $r\propto M^{2}$, whereas for more extended initial distributions, $r\propto M$. Our results show that collisionless mergers in a cosmological context can reduce the strength of inner dark matter cusps with changes in logarithmic slope of 0.3 to 0.5 at fixed radius. Shallow cusps such as those found recently in several strong lensing clusters thus do not necessarily conflict with CDM, but may rather reflect on the initial structure of the progenitor galaxies, which was shaped at high redshift by their formation process.

We study the formation of early-type galaxies (ETGs) through mergers with a sample of 70 high-resolution numerical simulations of binary mergers of disc galaxies. These simulations encompass various mass ratios, initial conditions and orbital parameters. We find that binary mergers of disc galaxies with mass ratios of 3:1 and 6:1 are nearly always classified as Fast Rotators according to the Atlas3D criterion: they preserve the structure of the input fast rotating spiral progenitors. Major disc mergers (mass ratios of 2:1 and 1:1) lead to both Fast and Slow Rotators. Most of the Slow Rotators hold a stellar Kinematically Distinct Core (KDC) in their 1-3 central kilo-parsec: these KDCs are built from the stellar components of the progenitors. The mass ratio of the progenitors is a fundamental parameter for the formation of Slow Rotators in binary mergers, but it also requires a retrograde spin for the progenitor galaxies with respect to the orbital angular momentum. The importance of the initial spin of the progenitors is also investigated in the library of galaxy mergers of the GalMer project.

Surveying eighteen 12CO-bright galaxies from the ATLAS3D early-type galaxy sample with the Institut de Radio Astronomie Millim\'etrique (IRAM) 30m telescope, we detect 13CO(1-0) and 13CO(2-1) in all eighteen galaxies, HCN(1-0) in 12/18 and HCO+(1-0) in 10/18. We find that the line ratios 12CO(1-0)/13CO(1-0) and 12CO(1-0)/HCN(1-0) are clearly correlated with several galaxy properties: total stellar mass, luminosity-weighted mean stellar age, molecular to atomic gas ratio, dust temperature and dust morphology. We suggest that these correlations are primarily governed by the optical depth in the 12CO lines; interacting, accreting and/or starbursting early-type galaxies have more optically thin molecular gas while those with settled dust and gas discs host optically thick molecular gas. The ranges of the integrated line intensity ratios generally overlap with those of spirals, although we note some outliers in the 12CO(1- 0)/13CO(1-0), 12CO(2-1)/13CO(2-1) and HCN/HCO+(1-0) ratios. In particular, three galaxies are found to have very low 12CO(1-0)/13CO(1-0) and 12CO(2-1)/13CO(2-1) ratios. Such low ratios may signal particularly stable molecular gas which creates stars less efficiently than 'normal' (i.e. below Schmidt-Kennicutt prediction), consistent with the low dust temperatures seen in these galaxies.

We present the Atlas3D HI survey of 166 nearby early-type galaxies (ETGs) down to M(HI)~10^7 M_sun. We detect HI in ~40% of all ETGs outside the Virgo cluster and in ~10% of all ETGs inside it. This demonstrates that it is common for non-cluster ETGs to host HI. The HI morphology varies from regular discs/rings (the majority of the detections) to unsettled gas distributions. The former are either small discs (M(HI)<10^8 M_sun) confined within the stellar body and sharing the same kinematics of the stars, or large discs/rings (M(HI) up to 5x10^9 M_sun) extending to tens of kpc from the host galaxy and frequently kinematically decoupled from the stars. Neutral hydrogen provides material for star formation in ETGs. Galaxies with central HI exhibit signatures of star formation in ~70% of the cases, ~5 times more frequently than galaxies without central HI. The central ISM is dominated by molecular gas. In ETGs with a small gas disc the conversion of HI into H_2 is as efficient as in spirals. The ETG HI mass function has M*~2x10^9 M_sun and slope=-0.7. ETGs host much less HI than spirals as a family. However, a significant fraction of them is as HI-rich as spirals. The main difference between ETGs and spirals is that the former lack the high-column-density HI typical of the bright stellar disc of the latter. We find an envelope of decreasing M(HI) with increasing environment density. The gas-richest ETGs live in the poorest environments (where star-formation is more common), galaxies in the centre of Virgo have the lowest HI content, and the cluster outskirts are a transition region. We find an HI morphology-density relation. At low environment density HI is mostly distributed on large discs/rings. More disturbed HI morphologies dominate environment densities typical of rich groups, confirming the importance of processes occurring on a galaxy-group scale for the evolution of ETGs.

We use the first systematic samples of CO millimeter emission in z>1 'main-sequence' star forming galaxies (SFGs) to study the metallicity dependence of the conversion factor {\alpha}CO, from CO line luminosity to molecular gas mass. The molecular gas depletion rate inferred from the ratio of the star formation rate (SFR) to CO luminosity, is ~1 Gyr-1 for near-solar metallicity galaxies with stellar masses above M_S~1e11 M_sun. In this regime the depletion rate does not vary more than a factor of two to three as a function of molecular gas surface density, or redshift between z~0 and 2. Below M_S the depletion rate increases rapidly with decreasing metallicity. We argue that this trend is not caused by starburst events, by changes in the physical parameters of the molecular clouds, or by the impact of the fundamental metallicity-SFR-stellar mass relation. A more probable explanation is that the conversion factor is metallicity dependent and that star formation can occur in 'CO-dark' gas. The trend is also expected theoretically from the effect of enhanced photodissociation of CO by ultraviolet radiation at low metallicity. From the available z~0 and z~1-3 samples we constrain the slope of the log({\alpha}CO) -log (metallicity) relation to range between -1 and -2, fairly insensitive to the assumed slope of the gas-star formation rate relation. Because of the lower metallicities near the peak of the galaxy formation activity at z~1-2 compared to z~0, we suggest that molecular gas masses estimated from CO luminosities have to be substantially corrected upward for galaxies below M_S.

We analyze 40 cosmological re-simulations of individual massive galaxies with present-day stellar masses of $M_{*} > 6.3 \times 10^{10} M_{\odot}$ in order to investigate the physical origin of the observed strong increase in galaxy sizes and the decrease of the stellar velocity dispersions since redshift $z \approx 2$. At present 25 out of 40 galaxies are quiescent with structural parameters (sizes and velocity dispersions) in agreement with local early type galaxies. At z=2 all simulated galaxies with $M_* \gtrsim 10^{11}M_{\odot}$ (11 out of 40) at z=2 are compact with projected half-mass radii of $\approx$ 0.77 ($\pm$0.24) kpc and line-of-sight velocity dispersions within the projected half-mass radius of $\approx$ 262 ($\pm$28) kms$^{-1}$ (3 out of 11 are already quiescent). Similar to observed compact early-type galaxies at high redshift the simulated galaxies are clearly offset from the local mass-size and mass-velocity dispersion relations. Towards redshift zero the sizes increase by a factor of $\sim 5-6$, following $R_{1/2} \propto (1+z)^{\alpha}$ with $\alpha = -1.44$ for quiescent galaxies ($\alpha = -1.12$ for all galaxies). The velocity dispersions drop by about one-third since $z \approx 2$, following $\sigma_{1/2} \propto (1+z)^{\beta}$ with $\beta = 0.44$ for the quiescent galaxies ($\beta = 0.37$ for all galaxies). The simulated size and dispersion evolution is in good agreement with observations and results from the subsequent accretion and merging of stellar systems at $z\lesssim 2$ which is a natural consequence of the hierarchical structure formation. A significant number of the simulated massive galaxies (7 out of 40) experience no merger more massive than 1:4 (usually considered as major mergers). On average, the dominant accretion mode is stellar minor mergers with a mass-weighted mass-ratio of 1:5. (abridged)

Many observed massive star-forming z\approx2 galaxies are large disks that exhibit irregular morphologies, with \sim1kpc, \sim10^(8-10)Msun clumps. We present the largest sample to date of high-resolution cosmological SPH simulations that zoom-in on the formation of individual M*\sim10^(10.5)Msun galaxies in \sim10^(12)Msun halos at z\approx2. Our code includes strong stellar feedback parameterized as momentum-driven galactic winds. This model reproduces many characteristic features of this observed class of galaxies, such as their clumpy morphologies, smooth and monotonic velocity gradients, high gas fractions (f_g\sim50%) and high specific star-formation rates (\gtrsim1Gyr^(-1)). In accord with recent models, giant clumps (Mclump\sim(5x10^8-10^9)Msun) form in-situ via gravitational instabilities. However, the galactic winds are critical for their subsequent evolution. The giant clumps we obtain are short-lived and are disrupted by wind-driven mass loss. They do not virialise or migrate to the galaxy centers as suggested in recent work neglecting strong winds. By phenomenologically implementing the winds that are observed from high-redshift galaxies and in particular from individual clumps, our simulations reproduce well new observational constraints on clump kinematics and clump ages. In particular, the observation that older clumps appear closer to their galaxy centers is reproduced in our simulations, as a result of inside-out formation of the disks rather than inward clump migration.

We present a detailed comparison between numerical cosmological hydrodynamic zoom simulations and semi-analytic models (SAMs) run within merger trees extracted from the simulations. The high-resolution simulations represent 48 individual halos with virial masses in the range 2.4*10^11M_sun < M_Halo < 3.3*10^13M_sun. They include radiative H & He cooling, photo-ionization, star formation and thermal SN feedback. We compare with different SAM versions including only this complement of physical processes, and also ones including supernova driven winds, metal cooling, and feedback from AGN. Our analysis is focused on the cosmic evolution of the baryon content in galaxies and its division into various components (stars, cold gas, and hot gas). Both the SAMs and simulations are compared with observational relations between halo mass and stellar mass, and between stellar mass and star formation rate, at low and high redshift. The simulations turn out to have much higher star formation efficiencies (by about a factor of ten) than the SAMs. Therefore the cold gas is consumed much more rapidly in the simulations and stars form much earlier. Also, simulations show a transition between stellar mass growth that is dominated by in situ formation of stars to growth that is predominantly through accretion of stars formed in external galaxies. In SAMs, stellar growth is always dominated by in situ star formation. In addition, SAMs overpredict the overall gas accretion rates relative to the simulations, and overestimate the fraction of "hot" relative to "cold" accretion. We discuss the reasons for these discrepancies, and identify several physical processes that are missing in our SAM. We also highlight physical processes that are neglected in the simulations studied here, but which appear to be crucial in order to understand the properties of real galaxies.

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.

We propose a simple model for the origin of fast and slow rotator early-type galaxies (ETG) within the hierarchical $\Lambda$CDM scenario, that is based on the assumption that the mass fraction of stellar discs in ETGs is a proxy for the specific angular momentum expressed via $\lambda_R$. Within our model we reproduce the fraction of fast and slow rotators as a function of magnitude in the \atl survey, assuming that fast rotating ETGs have at least 10% of their total stellar mass in a disc component. In agreement with \atl observations we find that slow rotators are predominantly galaxies with $ M_* > 10^{10.5}$ M$_{\odot}$ contributing $\sim 20%$ to the overall ETG population. We show in detail that the growth histories of fast and slow rotators are different, supporting the classification of ETGs into these two categories. Slow rotators accrete between $\sim 50% -90%$ of their stellar mass from satellites and their most massive progenitors have on average up to 3 major mergers during their evolution. Fast rotators in contrast, accrete less than 50% and have on average less than one major merger in their past. We find that the underlying physical reason for the different growth histories is the slowing down and ultimately complete shut-down of gas cooling in massive galaxies. Once cooling and associated star formation in disc stops, galaxies grow via infall from satellites. Frequent minor mergers thereby, destroy existing stellar discs via violent relaxation and also tend to lower the specific angular momentum of the main stellar body, lowering $\lambda_R$ into the slow rotator regime. Abridged...

We make use of interferometric CO and HI observations, and optical integral-field spectroscopy to probe the origin of the molecular and ionised interstellar medium (ISM) in local early-type galaxies (ETGs). We find that 36\pm5% of our sample of fast rotating ETGs have their ionised gas kinematically misaligned with respect to the stars, setting a strong lower limit on the importance of externally acquired gas (e.g. from mergers and cold accretion). Slow rotators have a flat distribution of misalignments, indicating that the dominant source of gas is external. The molecular, ionised and atomic gas in all the detected galaxies are always kinematically aligned, even when they are misaligned from the stars, suggesting that all these three phases of the ISM share a common origin. In addition, we find that the origin of the cold and warm gas in fast-rotating ETGs is strongly affected by environment, despite the molecular gas detection rate and mass fractions being fairly independent of group/cluster membership. Galaxies in dense groups and the Virgo cluster nearly always have their molecular gas kinematically aligned with the stellar kinematics, consistent with a purely internal origin. In the field, however, kinematic misalignments between the stellar and gaseous components indicate that >46% of local fast-rotating ETGs have their gas supplied from external sources. We discuss several scenarios which could explain the environmental dichotomy, but find it difficult to simultaneously explain the kinematic misalignment difference and the constant detection rate. Furthermore, our results suggest that galaxy mass may be an important independent factor associated with the origin of the gas, with the most massive fast-rotating galaxies in our sample (M_K<-24 mag; stellar mass of >8x10^10 Msun) always having kinematically aligned gas. (abridged)

Hydrodynamic simulations of galaxies with active galactic nuclei (AGN) have typically employed feedback that is purely local: i.e., an injection of energy to the immediate neighborhood of the black hole. We perform GADGET-2 simulations of massive elliptical galaxies with an additional feedback component: an observationally calibrated X-ray radiation field which emanates from the black hole and heats gas out to large radii from the galaxy center. We find that including the heating and radiation pressure associated with this X-ray flux in our simulations enhances the effects which are commonly reported from AGN feedback. This new feedback model is twice as effective as traditional feedback at suppressing star formation, produces 3 times less star formation in the last 6 Gyr, and modestly lowers the final BH mass (30%). It is also significantly more effective than an X-ray background in reducing the number of satellite galaxies.

The mass assembly of galaxies leaves imprints in their outskirts, such as shells and tidal tails. The frequency and properties of such fine structures depend on the main acting mechanisms - secular evolution, minor or major mergers - and on the age of the last substantial accretion event. We use this to constrain the mass assembly history of two apparently relaxed nearby Early-Type Galaxies (ETGs) selected from the Atlas-3D sample, NGC 680 and NGC 5557. Our ultra deep optical images obtained with MegaCam on the Canada-France-Hawaii Telescope reach 29 mag/arcsec^2 in the g-band. They reveal very low-surface brightness (LSB) filamentary structures around these ellipticals. Among them, a gigantic 160 kpc long tail East of NGC 5557 hosts gas-rich star-forming objects. NGC 680 exhibits two major diffuse plumes apparently connected to extended HI tails, as well as a series of arcs and shells. Comparing the outer stellar and gaseous morphology of the two ellipticals with that predicted from models of colliding galaxies, we argue that the LSB features are tidal debris, the star-forming objects near NGC 5557, long lived Tidal Dwarf Galaxies and that each of the two ETGs was assembled during a relatively recent, major wet merger, which likely occurred at a redshift below z = 0.5. The inner kinematics of NGC 680 is typical for fast rotators which make the bulk of nearby ETGs in the Atlas-3D sample. NGC 5557 belongs to the poorly populated class of massive, round, slow rotators that are predicted by semi-analytic models and cosmological simulations to be the end-product of a complex mass accretion history, involving ancient major mergers and more recent minor mergers. Our observations suggest that under specific circumstances a single binary merger may dominate the formation history of such objects and thus that at least some massive ETGs may form at relatively low redshift (abridged).

We study the formation of early-type galaxies through mergers with a sample of 70 high-resolution (softening length < 60 pc and 12*10^6 particles) numerical simulations of binary mergers of disc galaxies and 16 simulations of ETG remergers. These simulations, designed to accompany observations and models conducted within the Atlas3D project, encompass various mass ratios (from 1:1 to 6:1), initial conditions and orbital parameters. The progenitor disc galaxies are spiral-like with bulge to disc ratios typical of Sb and Sc galaxies. We find that binary mergers of disc galaxies with mass ratios of 3:1 and 6:1 are nearly always classified as Fast Rotators according to the Atlas3D criterion (based on the lambda_R parameter): they preserve the structure of the input fast rotating spiral progenitors. They have intrinsic ellipticities larger than 0.5, cover intrinsic lambda_R values between 0.2 and 0.6, within the range of observed Fast Rotators. Major disc mergers (mass ratios of 2:1 and 1:1) lead to both Fast and Slow Rotators. Most of the Fast Rotators produced in major mergers have intermediate flattening, with ellipticities between 0.4 and 0.6. Most Slow Rotators formed in these binary disc mergers hold a stellar Kinematically Distinct Core (KDC) in their 1-3 central kilo-parsec: these KDCs are built from the stellar components of the progenitors. Besides a handful of specific observed systems -- the counter-rotating discs (2-sigma galaxies) -- these therefore cannot reproduce the observed population of Slow Rotators in the nearby Universe. The mass ratio of the progenitors is a fundamental parameter for the formation of Slow Rotators in these binary mergers, but it also requires a retrograde spin for the earlier-type (Sb) progenitor galaxy with respect to the orbital angular momentum. (Abridged)

We look at the morphology of fast and slow rotator early-type galaxies. Edge-on fast rotators are lenticular galaxies. They appear like spiral galaxies with the gas and dust removed, and in some cases are flat ellipticals with disky isophotes. Fast rotators are often barred and span the same full range of bulge fractions as spiral galaxies. The slow rotators are rounder and are generally consistent with being genuine, namely spheroidal-like, elliptical galaxies. We propose a revision to the tuning-fork diagram by Hubble as it gives a misleading description of ETGs. We study for the first time the kinematic morphology-density T-Sigma relation using fast and slow rotators to replace lenticulars and ellipticals. We find that our relation is cleaner than using classic morphology. Slow rotators are nearly absent at the lowest density environments [f(SR)<2%] and generally constitute a small fraction [f(SR)~4%] of the total galaxy population in the environments explored by our survey, with the exception of the densest core of the Virgo cluster [f(SR)~20%]. We find a clean log-linear relation between the fraction f(Sp) of spiral galaxies and the local galaxy surface density. The existence of a smooth kinematic T-Sigma relation in the field excludes processes related to the cluster environment as main contributors to the apparent conversion of spirals into fast-rotators in low-density environments. It shows that the segregation is driven by local effects at the small-group scale. Only at the largest densities in the Virgo core does the f(Sp) relation break down and steepens sharply, while the fraction of slow-rotators starts to significantly increase. This suggests that a different mechanism is at work there. (Abridged)

We report the discovery of a powerful molecular wind from the nucleus of the non-interacting nearby S0 field galaxy NGC 1266. The single-dish CO profile exhibits emission to +/- 400 km/s and requires a nested Gaussian fit to be properly described. Interferometric observations reveal a massive, centrally- concentrated molecular component with a mass of 1.1x10^9 Msuns and a molecular outflow with a molecular mass of 2.4x10^7 Msuns. The molecular gas close to the systemic velocity consists of a rotating, compact nucleus with a mass of about 4.1x108 Msuns within a radius of approximately 60 pc. This compact molecular nucleus has a surface density of \approx 2.7 \times 10^4 Msuns/pc^2, more than two orders of magnitude larger than that of giant molecular clouds in the disk of the Milky Way, and it appears to sit on the Kennicutt-Schmidt relation despite its extreme kinematics and energetic activity. We interpret this nucleus as a disk that confines the outflowing wind. A mass outflow rate of 13 Msuns/yr leads to a depletion timescale of about 85 Myr. The star formation in NGC 1266 is insufficient to drive the outflow, and thus it is likely driven by the active galactic nucleus (AGN). The concentration of the majority of the molecular gas in the central 100 pc requires an extraordinary loss of angular momentum, but no obvious companion or interacting galaxy is present to enable the transfer. NGC 1266 is the first known outflowing molecular system that does not show any evidence of a recent interaction.

We present high resolution simulations of a multiple merger of three disk galaxies including the evolution of magnetic fields performed with the N-body/SPH code Gadget. For the first time, we embed the galaxies in a magnetized, low-density medium, thus modeling an ambient IGM. The simulations include radiative cooling and a model for star formation and supernova feedback. The progenitor disks have initial magnetic seed fields in the range of 10e-9 to 10e-6 G and the IGM has initial fields of 10e-12 to 10e-9 G. The simulations are compared to a run excluding magnetic fields. We show that the propagation of interaction-driven shocks depends significantly on the initial magnetic field strength. The shocks propagate faster in simulations with stronger initial field, suggesting that the shocks are supported by magnetic pressure. The Mach numbers of the shocks range from approximately M=1.5 for the non-magnetized case up to M=6 for the highest initial magnetization, resulting in higher temperatures of the shock heated IGM gas. The magnetic field in the system saturates rapidly after the mergers at ~ 10e-6 G within the galaxies and ~ 10e-8 G in the IGM independent of the initial value. These field strengths agree with observed values and correspond to the equipartition value of the magnetic pressure with the turbulent pressure in the system. We also present synthetic radio and polarization maps for different phases of the evolution showing that shocks driven by the interaction produce a high amount of polarized emission. These idealized simulations indicate that magnetic fields play an important role for the hydrodynamics of the IGM during galactic interactions. We also show that even weak seed fields are efficiently strengthened during multiple galactic mergers. This interaction driven amplification might have been a key process for the magnetization of the Universe.

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.

We reexamine the age distribution of star clusters in the Antennae in the context of N-body+hydrodynamical simulations of these interacting galaxies. All of the simulations that account for the observed morphology and other properties of the Antennae have star formation rates that vary relatively slowly with time, by factors of only 1.3 - 2.5 in the past 10^8 yr. In contrast, the observed age distribution of the clusters declines approximately as a power law, dN/dt \propto t^{gamma} with gamma = -1.0, for ages 10^6 yr \la t \la 10^9 yr. These two facts can only be reconciled if the clusters are disrupted progressively for at least 10^8 yr and possibly 10^9 yr. When we combine the simulated formation rates with a power-law model, f_surv \propto t^{delta}, for the fraction of clusters that survive to each age t, we match the observed age distribution with exponents in the range -0.9 \la delta \la -0.6 (with a slightly different delta for each simulation). The similarity between delta and gamma indicates that dN/dt is shaped mainly by the disruption of clusters rather than variations in their formation rate. Thus, the situation in the interacting Antennae resembles that in relatively quiescent galaxies such as the Milky Way and the Magellanic Clouds.

The fragmentation of shocked flows in a thermally bistable medium provides a natural mechanism to form turbulent cold clouds as precursors to molecular clouds. Yet because of the large density and temperature differences and the range of dynamical scales involved, following this process with numerical simulations is challenging. We compare two-dimensional simulations of flow-driven cloud formation without self-gravity, using the Lagrangian Smoothed Particle Hydrodynamics (SPH) code VINE and the Eulerian grid code Proteus. Results are qualitatively similar for both methods, yet the variable spatial resolution of the SPH method leads to smaller fragments and thinner filaments, rendering the overall morphologies different. Thermal and hydro-dynamical instabilities lead to rapid cooling and fragmentation into cold clumps with temperatures below 300K. For clumps more massive than 1 Msun/pc, the clump mass function has an average slope of -0.8. The internal velocity dispersion of the clumps is nearly an order of magnitude smaller than their relative motion, rendering it subsonic with respect to the internal sound speed of the clumps, but supersonic as seen by an external observer. For the SPH simulations most of the cold gas resides at temperatures below 100K, while the grid-based models show an additional, substantial component between 100 and 300K. Independently of the numerical method our models confirm that converging flows of warm neutral gas fragment rapidly and form high-density, low-temperature clumps as possible seeds for star formation.

[Abridged] We use the ATLAS3D sample of 260 early-type galaxies (ETGs) to study the apparent kinematic misalignment angle, Psi, defined as the angle between the photometric and kinematic major axis. We find that 71% of nearby ETGs are strictly aligned systems (Psi > 5 deg), an additional 14% have 5 < Psi < 10 deg and 90% of galaxies have Psi < 15 deg. Taking into account measurement uncertainties, 90% of galaxies can be considered aligned to better than 5 deg, suggesting that only a small fraction of early-type galaxies (~10%) are not consistent with axisymmetry within the projected half-light radius. We identify morphological features such as bars and rings (30%), dust structures (16%), blue nuclear colours (6%) and evidence of interactions (8%). We use kinemetry to analyse the mean velocity maps and separate galaxies in two broad types of regular and non-regular rotators. We find 82% of regular rotators and 17% non-regular rotators, with 2 galaxies that we were not able to classify due to data quality. The non-regular rotators are typically found in dense regions and are massive. The majority of galaxies do not have any specific kinemetric features, but we highlight here the frequency of the kinematically distinct cores (7%) and the aligned double peaks in the velocity dispersion maps (4%). Most of the galaxies that are misaligned have complex kinematics and are non-regular rotators. While the trends are weak, there is a tendency that large values of Psi are found in galaxies at intermediate environmental densities and among the most massive galaxies in the sample. We suggest that the most common formation mechanism for ETGs preserves the axisymmetry of the disk progenitors. Less commonly, the formation process results in a triaxial galaxy with much lower net angular momentum.

We have carried out a survey for 12CO J=1-0 and J=2-1 emission in the 260 early-type galaxies of the volume-limited Atlas3D sample, with the goal of connecting their star formation and assembly histories to their cold gas content. This is the largest volume-limited CO survey of its kind and is the first to include many Virgo Cluster members. Sample members are dynamically hot galaxies with a median stellar mass 3\times 10^{10} Msun; they are selected by morphology rather than colour, and the bulk of them lie on the red sequence. The overall CO detection rate is 56/259 = 0.22 \error 0.03, with no dependence on K luminosity and only a modest dependence on dynamical mass. There are a dozen CO detections among the Virgo Cluster members; statistical analysis of their H_2 mass distributions and their dynamical status within the cluster shows that the cluster's influence on their molecular masses is subtle at best, even though (unlike spirals) they seem to be virialized within the cluster. We suggest that the cluster members have retained their molecular gas through several Gyr residences in the cluster. There are also a few extremely CO-rich early-type galaxies with H_2 masses >= 10^9 Msun, and these are in low density environments. We do find a significant trend between molecular content and the stellar specific angular momentum. The galaxies of low angular momentum also have low CO detection rates, suggesting that their formation processes were more effective at destroying molecular gas or preventing its re-accretion. We speculate on the implications of these data for the formation of various sub-classes of early-type galaxies.

We demonstrate here using both single-dish and interferometric observations that CO molecules are an excellent kinematic tracer, even in high-mass galaxies, allowing us to investigate for the first time the CO Tully-Fisher relation of early-type galaxies. We compare the Tully-Fisher relations produced using both single-dish and interferometric data and various inclination estimation methods, and evaluate the use of the velocity profile shape as a criterion for selecting galaxies in which the molecular gas extends beyond the peak of the rotation curve. We show that the gradient and zero-point of the best-fit relations are robust, independent of the velocity measure and inclination used, and agree with those of relations derived using stellar kinematics. We also show that the early-type CO Tully-Fisher (CO-TF) relation is offset from the CO-TF of spirals by 0.98 \pm 0.22 magnitude at Ks-band, in line with other results. The intrinsic scatter of the relation is found to be \approx 0.4 magnitudes, similar to the level found in the spiral galaxy population. Next generation facilities such as the Large Millimeter Telescope (LMT) and the Atacama Large Millimeter/Sub-millimeter Array (ALMA) should allow this technique to be used in higher-redshift systems, providing a simple new tool to trace the mass-to-light ratio evolution of the most massive galaxies over cosmic time.

We provide a census of the apparent stellar angular momentum within 1 Re of a volume-limited sample of 260 early-type galaxies (ETGs) in the nearby Universe, using integral-field spectroscopy obtained in the course of the ATLAS3D project. We exploit the LambdaR parameter to characterise the existence of two families of ETGs: Slow Rotators which exhibit complex stellar velocity fields and often include stellar kinematically Distinct Cores (KDCs), and Fast Rotators which have regular velocity fields. Our complete sample of 260 ETGs leads to a new criterion to disentangle Fast and Slow Rotators which now includes a dependency on the apparent ellipticity (Epsilon). It separates the two classes significantly better than the previous prescription, and than a criterion based on V/Sigma: Slow Rotators and Fast Rotators have LambdaR lower and larger than kFSxSQRT(Epsilon), respectively, where kFS=0.31 for measurements made within 1 Re. We show that the vast majority of early-type galaxies are Fast Rotators: these have regular stellar rotation, with aligned photometric and kinematic axes (Paper II, Krajnovic et al. 2011}, include discs and often bars and represent 86% (224/260) of all early-type galaxies in the volume-limited ATLAS3D sample. Fast Rotators span the full range of apparent ellipticities from 0 to 0.85, and we suggest that they cover intrinsic ellipticities from about 0.35 to 0.85, the most flattened having morphologies consistent with spiral galaxies. Only a small fraction of ETGs are Slow Rotators representing 14% (36/260) of the ATLAS3D sample of ETGs. Of all Slow Rotators, 11% (4/36) exhibit two counter-rotating stellar disc-like components and are rather low mass objects (Mdyn<10^10.5 M_Sun). All other Slow Rotators (32/36) appear relatively round on the sky (Epsilon_e<0.4), tend to be massive (Mdyn>10^10.5 M_Sun), and often (17/32) exhibit KDCs.

We have studied the properties of giant star forming clumps in five z~2 star-forming disks with deep SINFONI AO spectroscopy at the ESO VLT. The clumps reside in disk regions where the Toomre Q-parameter is below unity, consistent with their being bound and having formed from gravitational instability. Broad H{\alpha}/[NII] line wings demonstrate that the clumps are launching sites of powerful outflows. The inferred outflow rates are comparable to or exceed the star formation rates, in one case by a factor of eight. Typical clumps may lose a fraction of their original gas by feedback in a few hundred million years, allowing them to migrate into the center. The most active clumps may lose much of their mass and disrupt in the disk. The clumps leave a modest imprint on the gas kinematics. Velocity gradients across the clumps are 10-40 km/s/kpc, similar to the galactic rotation gradients. Given beam smearing and clump sizes, these gradients may be consistent with significant rotational support in typical clumps. Extreme clumps may not be rotationally supported; either they are not virialized, or they are predominantly pressure supported. The velocity dispersion is spatially rather constant and increases only weakly with star formation surface density. The large velocity dispersions may be driven by the release of gravitational energy, either at the outer disk/accreting streams interface, and/or by the clump migration within the disk. Spatial variations in the inferred gas phase oxygen abundance are broadly consistent with inside-out growing disks, and/or with inward migration of the clumps.

The molecular gas content of local early-type galaxies is constrained and discussed in relation to their evolution. First, as part of the Atlas3D survey, we present the first complete, large (260 objects), volume-limited single-dish survey of CO in normal local early-type galaxies. We find a surprisingly high detection rate of 22%, independent of luminosity and at best weakly dependent on environment. Second, the extent of the molecular gas is constrained with CO synthesis imaging, and a variety of morphologies is revealed. The kinematics of the molecular gas and stars are often misaligned, implying an external gas origin in over a third of the systems, although this behaviour is drastically diffferent between field and cluster environments. Third, many objects appear to be in the process of forming regular kpc-size decoupled disks, and a star formation sequence can be sketched by piecing together multi-wavelength information on the molecular gas, current star formation, and young stars. Last, early-type galaxies do not seem to systematically obey all our usual prejudices regarding star formation, following the standard Schmidt-Kennicutt law but not the far infrared-radio correlation. This may suggest a greater diversity in star formation processes than observed in disk galaxies. Using multiple molecular tracers, we are thus starting to probe the physical conditions of the cold gas in early-types.

The mass assembly of galaxies leaves various imprints on their surroundings, such as shells, streams and tidal tails. The frequency and properties of these fine structures depend on the mechanism driving the mass assembly: e.g. a monolithic collapse, rapid cold-gas accretion followed by violent disk instabilities, minor mergers or major dry / wet mergers. Therefore, by studying the outskirts of galaxies, one can learn about their main formation mechanism. I present here our on-going work to characterize the outskirts of Early-Type Galaxies (ETGs), which are powerful probes at low redshift of the hierarchical mass assembly of galaxies. This work relies on ultra-deep optical images obtained at CFHT with the wide-field of view MegaCam camera of field and cluster ETGs obtained as part of the Atlas-3D and NGVS projects. State of the art numerical simulations are used to interpret the data. The images reveal a wealth of unknown faint structures at levels as faint as 29 mag arcsec-2 in the g-band. Initial results for two galaxies are presented here.

Cosmological simulations of galaxy formation appear to show a two-phase character with a rapid early phase at z>2 during which in-situ stars are formed within the galaxy from infalling cold gas followed by an extended phase since z<3 during which ex-situ stars are primarily accreted. In the latter phase massive systems grow considerably in mass and radius by accretion of smaller satellite stellar systems formed at quite early times (z>3) outside of the virial radius of the forming central galaxy. These tentative conclusions are obtained from high resolution re-simulations of 39 individual galaxies in a full cosmological context with present-day virial halo masses ranging from 7e11 M_sun h^-1 < M_vir < 2.7e13 M_sun h^-1 and central galaxy masses between 4.5e10 M_sun h^-1 < M_* < 3.6e11 M_sun h^-1. The simulations include the effects of a uniform UV background, radiative cooling, star formation and energetic feedback from SNII. The importance of stellar accretion increases with galaxy mass and towards lower redshift. In our simulations lower mass galaxies ($M_* < 9e10 M_sun h^-1) accrete about 60 per cent of their present-day stellar mass. High mass galaxy ($M_* > 1.7e11 M_sun h^-1) assembly is dominated by accretion and merging with about 80 per cent of the stars added by the present-day. In general the simulated galaxies approximately double their mass since z=1. For massive systems this mass growth is not accompanied by significant star formation. The majority of the in-situ created stars is formed at z>2, primarily out of cold gas flows. We recover the observational result of archaeological downsizing, where the most massive galaxies harbor the oldest stars. We find that this is not in contradiction with hierarchical structure formation. Most stars in the massive galaxies are formed early on in smaller structures, the galaxies themselves are assembled late.

The Atlas3D project is a multi-wavelength survey combined with a theoretical modeling effort. The observations span from the radio to the millimeter and optical, and provide multi-colour imaging, two-dimensional kinematics of the atomic (HI), molecular (CO) and ionized gas (Hbeta, [OIII] and [NI]), together with the kinematics and population of the stars (Hbeta, Fe5015 and Mgb), for a carefully selected, volume-limited (1.16*10^5 Mpc^3) sample of 260 early-type (elliptical E and lenticular S0) galaxies (ETGs). The models include semi-analytic, N-body binary mergers and cosmological simulations of galaxy formation. Here we present the science goals for the project and introduce the galaxy sample and the selection criteria. The sample consists of nearby (D<42 Mpc) morphologically-selected ETGs extracted from a parent sample of 871 galaxies (8% E, 22% S0 and 70% spirals) brighter than M_K<-21.5 mag (stellar mass M_Star>6*10^9 M_Sun). We analyze possible selection biases and we conclude that the parent sample is essentially complete and statistically representative of the nearby galaxy population. We present the size-luminosity relation for the spirals and ETGs and show that the ETGs in the Atlas3D sample define a tight red sequence in a colour-magnitude diagram, with few objects in the transition from the blue cloud. We describe the strategy of the SAURON integral-field observations and the extraction of the stellar kinematics with the pPXF method. We give an overview of the characteristics of the other main datasets already available for our sample and of the ongoing modelling projects.

The structure of a sample of high-redshift (z=2), rotating galaxies with high star formation rates and turbulent gas velocities of sigma=40-80 km/s is investigated. Fitting the observed disk rotational velocities and radii with a Mo, Mao, White (1998) (MMW) model requires unusually large disk spin parameters lambda_d>0.1 and disk-to-dark halo mass fraction m_d=0.2, close to the cosmic baryon fraction. The galaxies segregate into dispersion-dominated systems with 1<vmax/sigma<3, maximum rotational velocities vmax<200 km/s and disk half-light radii rd=1-3 kpc and rotation-dominated systems with vmax>200 km/s, vmax/sigma>3 and rd=4-8 kpc. For the dispersion-dominated sample, radial pressure gradients partly compensate the gravitational force, reducing the rotational velocities. Including this pressure effect in the MMW model, dispersion-dominated galaxies can be fitted well with spin parameters lf lambda_d=0.03-0.05 for high disk mass fractions of m_d=0.2 and with lambda_d=0.01-0.03 for m_d=0.05. These values are in good agreement with cosmological expectations. For the rotation-dominated sample however pressure effects are small and better agreement with theoretically expected disk spin parameters can only be achieved if the dark halo mass contribution in the visible disk regime (2-3*rd) is smaller than predicted by the MMW model. We argue that these galaxies can still be embedded in standard cold dark matter halos if the halos did not contract adiabatically in response to disk formation. It is shown that the observed high turbulent gas motions of the galaxies are consistent with a Toomre instability parameter Q=1 which is equal to the critical value, expected for gravitational disk instability to be the major driver of turbulence. The dominant energy source of turbulence is then the potential energy of the gas in the disk.

Even though the dark-matter power spectrum in the absence of biasing predicts a number density of halos n(M) ~ M^-2 (i.e. a Schechter alpha value of -2) at the low-mass end (M < 10^10 M_solar), hydrodynamic simulations have typically produced values for stellar systems in good agreement with the observed value alpha ~ -1. We explain this with a simple physical argument and show that an efficient external gas-heating mechanism (such as the UV background included in all hydro codes) will produce a critical halo mass below which halos cannot retain their gas and form stars. We test this conclusion with GADGET-2-based simulations using various UV backgrounds, and for the first time we also investigate the effect of an X-ray background. We show that at the present epoch alpha is depends primarily on the mean gas temperature at the star-formation epoch for low-mass systems (z <~ 3): with no background we find alpha ~ -1.5, with UV only alpha ~ -1.0, and with UV and X-rays alpha ~ -0.75. We find the critical final halo mass for star formation to be ~4x10^8 M_solar with a UV background and ~7x10^8 M_solar with UV and X-rays.

We present spatially resolved integral field spectroscopic K-band data at a resolution of 0.13" (60pc) and interferometric CO(2-1) line observations of the prototypical merging system NGC6240. Despite the clear rotational signature, the stellar kinematics in the two nuclei are dominated by dispersion. We use Jeans modelling to derive the masses and the mass-to-light ratios of the nuclei. Combining the luminosities with the spatially resolved Br-gamma equivalent width shows that only 1/3 of the K-band continuum from the nuclei is associated with the most recent star forming episode; and that less than 30% of the system's bolometric luminosity and only 9% of its stellar mass is due to this starburst. The star formation properties, calculated from typical merger star formation histories, demonstrate the impact of different assumptions about the star formation history. The properties of the nuclei, and the existence of a prominent old stellar population, indicate that the nuclei are remnants of the progenitor galaxies' bulges.

We study the structural evolution of turbulent molecular clouds under the influence of ionizing radiation emitted from a nearby massive star by performing a high resolution parameter study with the iVINE code. The temperature is taken to be 10K or 100K, the mean number density is either 100cm^3 or 300cm^3. Besides, the turbulence is varied between Mach 1.5 and Mach 12.5 and the main driving scale between 1pc and 8pc. We vary the ionizing flux by an order of magnitude. In our simulations the ionizing radiation enhances the initial turbulent density distribution and thus leads to the formation of pillar-like structures observed adjacent to HII regions in a natural way. Gravitational collapse occurs regularly at the tips of the structures. We find a clear correlation between the initial state of the turbulent cold cloud and the final morphology and physical properties of the structures formed. The most favorable regime for the formation of pillars is Mach 4-10. Structures and therefore stars only form if the initial density contrast between the high density unionized gas and the gas that is going to be ionized is lower than the temperature contrast between the hot and the cold gas. The density of the resulting pillars is determined by a pressure equilibrium between the hot and the cold gas. A thorough analysis of the simulations shows that the complex kinematical and geometrical structure of the formed elongated filaments reflects that of observed pillars to an impressive level of detail. In addition, we find that the observed line-of sight velocities allow for a distinct determination of different formation mechanisms. Comparing the current simulations to previous results and recent observations we conclude that e.g. the pillars of creation in M16 formed by the mechanism proposed here and not by the radiation driven implosion of pre-existing clumps.

We study the growth of dark matter halos in the concordance LCDM cosmology using several N-body simulations of large cosmological volumes. We build merger trees from the Millennium and Millennium-II simulations, covering a range 10^9-10^15 Msun in halo mass and 1-10^5 in merger mass ratio. Our algorithm takes special care of halo fragmentation and ensures that the mass contribution of each merger to halo growth is only counted once. This way the integrated merger rate converges and we can consistently determine the contribution of mergers of different mass ratios to halo growth. We find that all resolved mergers, up to mass ratios of 10^5, contribute only ~60% of the total halo mass growth, while major mergers are subdominant, e.g. mergers with mass ratios smaller than 3:1 (10:1) contribute only ~20% (~30%). This is verified with an analysis of two additional simulation boxes, where we follow all particles individually throughout cosmic time. Our results are also robust against using several halo definitions. Under the assumption that the power-law behaviour of the merger rate at large mass ratios can be extrapolated to arbitrarily large mass ratios, it is found that, independently of halo mass, ~40% of the mass in halos comes from genuinely smooth accretion of dark matter that was never bound in smaller halos. We discuss possible implications of our findings for galaxy formation. One implication, assuming as is standard that the pristine intergalactic medium is heated and photoionized by UV photons, is that all halos accrete >40% of their baryons in smooth "cold" T>~10^4K gas, rather than as warm, enriched or clumpy gas or as stars.

We present self-consistent high-resolution simulations of NGC4038/4039 (the "Antennae galaxies") including star formation, supernova feedback and magnetic fields performed with the N-body/SPH code Gadget, in which magnetohydrodynamics are followed with the SPH method. We vary the initial magnetic field in the progenitor disks from 1 nG to 100 muG. At the time of the best match with the central region of the Antennae system the magnetic field has been amplified by compression and shear flows to an equilibrium field of approximately 10 muG, independent of the initial seed field. These simulations are a proof of the principle that galaxy mergers are efficient drivers for the cosmic evolution of magnetic fields. We present a detailed analysis of the magnetic field structure in the central overlap region. Simulated radio and polarization maps are in good morphological and quantitative agreement with the observations. In particular, the two cores with the highest synchrotron intensity and ridges of regular magnetic fields between the cores and at the root of the southern tidal arm develop naturally in our simulations. This indicates that the simulations are capable of realistically following the evolution of the magnetic fields in a highly non-linear environment. We also discuss the relevance of the amplification effect for present day magnetic fields in the context of hierarchical structure formation.

We present results on the evolution of the intrinsic scatter of black hole masses considering different implementations of a model in which black holes only grow via mergers. We demonstrate how merger driven growth affects the correlations between black hole mass and host bulge mass. The simple case of an initially log-normal distributed scatter in black hole and bulge masses combined with random merging within the galaxy population results in a decreasing scatter with merging generation/number as predicted by the Central-limit theorem. In general we find that the decrease in scatter {\sigma} is well approximated by {\sigma}merg(m) = {\sigma}ini \times (m + 1)^(-a/2) with a = 0.42 for a range of mean number of mergers m < 50. For a large mean number of mergers (m > 100) we find a convergence to a = 0.61. This is valid for a wide range of different initial distributions, refill-scenarios or merger mass-ratios. Growth scenarios based on halo merger trees of a (100 Mpc)^3 dark matter LambdaCDM-simulation show a similar behaviour with a scatter decrease of a = 0.30 with typical number of mergers m < 50 consistent with random merging (best matching model: a = 0.34). Assuming a present day scatter of 0.3 dex in black hole mass and a mean number of mergers not exceeding m = 50 our results imply a scatter of 0.6 dex at z = 3 and thus a possible scenario in which overmassive (and undermassive) black holes at high redshift are a consequence of a larger intrinsic scatter in black hole mass. A simple toy model connecting the growth of black holes to the growth of LambdaCDM dark matter halos via mergers, neglecting any contribution from accretion, yields a consistent M\cdot -MBulge relation at z = 0 - if we assume the correct initial relation.

We use the first systematic data sets of CO molecular line emission in z~1-3 normal star forming galaxies for a comparison of the dependence of galaxy-averaged star formation rates on molecular gas masses at low and high redshifts, and in different galactic environments. Although the current high-z samples are still small and biased toward the luminous and massive tail of the actively star-forming 'main-sequence', a fairly clear picture is emerging. Independent of whether galaxy integrated quantities or surface densities are considered, low- and high-z SFG galaxy populations appear to follow similar molecular gas-star formation relations with slopes 1.1 to 1.2. The gas-depletion time scale in these SFGs grows from 0.5 Gyrs at z~2 to 1.5 Gyrs at z~0. Because star formation depletion times are significantly smaller than the Hubble time at all redshifts sampled, star formation rates and gas fractions are set by the balance between gas accretion from the halo and stellar feedback. In contrast, very luminous gas rich major mergers at both low-z and high-z produce on average 4 to10 times more far-infrared luminosity per unit gas mass. Only some fraction of this difference can be explained by uncertainties in gas-mass or luminosity estimators; much of it must be intrinsic. The most likely interpretation is that the star formation relation is driven by global dynamical effects. For a given mass, the more compact merger systems produce stars more rapidly because their gas clouds are more compressed with shorter dynamical times, so that they churn more quickly through the available gas reservoir than the typical normal disk galaxies. When the dependence on galactic dynamical time scale is explicitly included, disk galaxies and mergers appear to follow similar gas to star-formation relations. The mergers may be forming stars at slightly higher efficiencies than the disks.

We study resolution effects in numerical simulations of gas-rich and gas-poor major mergers, and show that the formation of slowly-rotating elliptical galaxies often requires a resolution that is beyond the present-day standards to be properly modelled. Our sample of equal-mass merger models encompasses various masses and spatial resolutions, ranging from about 200pc and 10^5 particles per component, typical of some recently published major merger simulations, to up to 32pc and 10^3 M_sun in simulations using 2.4 x 10^7 collisionless particles and 1.2 x 10^7 gas particles, among the highest resolutions reached so far for gas-rich major merger of massive disc galaxies. We find that the formation of fast-rotating early-type galaxies, that are flattened by a significant residual rotation, is overall correctly reproduced at all such resolutions. However, the formation of slow-rotating early-type galaxies, which have a low residual angular momentum and are supported mostly by anisotropic velocity dispersions, is strongly resolution-dependent. The evacuation of angular momentum from the main stellar body is largely missed at standard resolution, and systems that should be slow rotators are then found to be fast rotators. The effect is most important for gas-rich mergers, but is also witnessed in mergers with an absent or modest gas component. The effect is robust with respect to our initial conditions and interaction orbits, and originates in the physical treatment of the relaxation process during the coalescence of the galaxies. Our findings show that a high-enough resolution is required to accurately model the global properties of merger remnants and the evolution of their angular momentum. The role of gas-rich mergers of spiral galaxies in the formation of slow-rotating ellipticals may therefore have been underestimated.

We present a new high-resolution N-body/SPH simulation of an encounter of two gas-rich disk galaxies which closely matches the morphology and kinematics of the interacting Antennae galaxies (NGC 4038/39). The simulation includes radiative cooling, star formation and feedback from SNII. The large-scale morphology and kinematics are determined by the internal structure and the orbit of the progenitor disks. The properties of the central region, in particular the starburst in the overlap region, only match the observations for a very short time interval of ~20 Myr after the second encounter. This indicates that the Antennae galaxies are in a special phase only about 40 Myr after the second encounter and 50 Myr before their final collision. This is the only phase in the simulation when a gas-rich overlap region between the nuclei is forming accompanied by enhanced star formation. The star formation rate as well as the recent star formation history in the central region agree well with observational estimates. For the first time this new model explains the distributed extra-nuclear star formation in the Antennae galaxies as a consequence of the recent second encounter. The proposed model predicts that the Antennae are in a later merger stage than the Mice (NGC 4676) and would therefore lose their first place in the classical Toomre sequence.

Stars form from cold molecular interstellar gas. Since this is relatively rare in the local Universe, galaxies like the Milky Way form only a few new stars per year. Typical massive galaxies in the distant Universe formed stars an order of magnitude more rapidly. Unless star formation was significantly more efficient, this difference suggests that young galaxies were much more gas rich. Molecular gas observations in the distant Universe have so far been largely restricted to very luminous, rare objects, including mergers and quasars. Here we report the results of a systematic survey of molecular gas in samples of typical massive star forming galaxies at <z>~1.2 and 2.3, when the Universe was 40% and 24% of its current age. Our measurements provide empirical evidence that distant star forming galaxies indeed were gas rich, and that the star formation efficiency is not strongly dependent on cosmic epoch. The average fraction of cold gas relative to total galaxy baryonic mass at z= 2.3 and z=1.2 is ~44% and 34%, three to ten times higher than in today's massive spiral galaxies. The slow decrease between z~2 and 1 probably requires a mechanism of semi-continuous replenishment of fresh gas to the young galaxies.

We present a preliminary analysis of the HI properties of early-type galaxies in the ATLAS3D sample. Using WSRT data for ~100 galaxies outside the Virgo cluster and data from the Alfalfa project for galaxies inside Virgo, we discuss the dependence of HI properties on environment. We detect HI in about half of the galaxies outside Virgo. For these systems, the HI morphology and kinematics change as a function of environment, going from regular, rotating systems around isolated galaxies to progressively more disturbed structures for galaxies with neighbours or in groups. In denser environment, inside Virgo, nearly none of the galaxies contains HI. We discuss future work in this field which will be enabled by next-generation, pre-SKA radio instruments. We present a simulated Apertif HI observation of an ATLAS3D early-type galaxy, showing how its appearance and detection level vary as a function of redshift.

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.

We study the evolution of black holes (BHs) on the M_BH-sigma and M_BH-M_bulge planes as a function of time in disk galaxies undergoing mergers. We begin the simulations with the progenitor black hole masses being initially below (Delta log M_BH=-2), on (Delta log M_BH=0) and above (Delta log M_BH=0.5) the observed local relations. The final relations are rapidly established after the final coalescense of the galaxies and their BHs. Progenitors with low initial gas fractions (f_gas=0.2) starting below the relations evolve onto the relations (Delta log M_BH=-0.18), progenitors on the relations stay there (Delta log M_BH=0) and finally progenitors above the relations evolve towards the relations, but still remaining above them (Delta log M_BH=0.35). Mergers in which the progenitors have high initial gas fractions (f_gas=0.8) evolve above the relations in all cases (Delta log M_BH=0.5). We find that the initial gas fraction is the prime source of scatter in the observed relations, dominating over the scatter arising from the evolutionary stage of the merger remnants. The fact that BHs starting above the relations do not evolve onto the relations, indicates that our simulations rule out the scenario in which overmassive BHs evolve onto the relations through gas-rich mergers. By implication our simulations thus disfavor the picture in which supermassive BHs develop significantly before their parent bulges.

The disruptive effect of galactic tides is a textbook example of gravitational dynamics. However, depending on the shape of the potential, tides can also become fully compressive. When that is the case, they might trigger or strengthen the formation of galactic substructures (star clusters, tidal dwarf galaxies), instead of destroying them. We perform N-body simulations of interacting galaxies to quantify this effect. We demonstrate that tidal compression occurs repeatedly during a galaxy merger, independently of the specific choice of parameterization. With a model tailored to the Antennae galaxies, we show that the distribution of compressive tides matches the locations and timescales of observed substructures. After extending our study to a broad range of parameters, we conclude that neither the importance of the compressive tides (~15% of the stellar mass) nor their duration (~ 10 Myr) are strongly affected by changes in the progenitors' configurations and orbits. Moreover, we show that individual clumps of matter can enter compressive regions several times in the course of a simulation. We speculate that this may spawn multiple star formation episodes in some star clusters, through e.g., enhanced gas retention.

We find that the amount and nature of the assumed ionizing background can strongly affect galaxy formation and evolution. Galaxy evolution simulations typically incorporate an ultraviolet background which falls off rapidly above z=3; e.g., that of Haardt & Madau (1996). However, this decline may be too steep to fit the WMAP constraints on electron scattering optical depth or observations of intermediate redshift (z ~ 2-4) Ly-alpha forest transmission. As an alternative, we present simulations of the cosmological formation of individual galaxies with UV backgrounds that decline more slowly at high redshift: both a simple intensity rescaling and the background recently derived by Faucher-Giguere (2009), which softens the spectrum at higher redshifts. We also test an approximation of the X-ray background with a similar z-dependence. We find for the test galaxies that an increase in either the intensity or hardness of ionizing radiation generically pushes star formation towards lower redshifts: although overall star formation in the simulation boxes is reduced by 10-25%, the galaxies show a factor of ~2 increase in the fraction of stars within a 30 kpc radius that are formed after z=1. Other effects include late gas inflows enhanced up to 30 times, stellar half-mass radii decreased by up to 30%, central velocity dispersions increased up to 40%, and a strong reduction in substructure. The magnitude of the effects depends on the environmental/accretion properties of the particular galaxy.

Physics lectures always refer to the tides as a disruptive effect. However, tides can also be compressive. When the potential of two galaxies overlap, as happens during a merger, fully compressive tides can develop and have a strong impact on the dynamics of substructures such as star clusters or tidal dwarf galaxies. Using N-body simulations of a large set of mergers, we noticed the importance of these tidal modes at cluster scale. With a model of the Antennae galaxies, we conclude that the positions and timescales of these tidal modes match the actual distribution of young clusters. A detailed study of the statistics of the compressive tides shows a stunning correlation between this purely gravitational effect and the observed properties of the star clusters. In this contribution, we introduce the concept of compressive tide and show its relevance in the special case of the Antennae galaxies. We extend our conclusions to a broad range of parameters and discuss their implications on several critical points such as the infant mortality, multiple star formation epochs in clusters or the age distribution.

We investigate the collapse and fragmentation of low-mass, trans-sonically turbulent prestellar cores, using SPH simulations. The initial conditions are slightly supercritical Bonnor-Ebert spheres, all with the same density profile, the same mass (M_O=6.1 Msun) and the same radius (R_O=17,000 AU), but having different initial turbulent velocity fields. Four hundred turbulent velocity fields have been generated, all scaled so that the mean Mach number is M=1. Then a subset of these, having a range of net angular momenta, j, has been evolved. The evolution of these turbulent cores is not strongly correlated with j. Instead it is moderated by the formation of filamentary structures due to converging turbulent flows. A high fraction (~ 82%) of the protostars forming from turbulent cores are attended by protostellar accretion discs, but only a very small fraction (~16%) of these discs is sufficiently cool and extended to develop non-linear gravitational instabilities and fragment.

We study the origin of the scaling relations of early-type galaxies (ETGs) by constructing detailed models of the stellar dynamics for the K-band selected, volume-limited ATLAS3D sample of 263 nearby ETGs, spanning a large range of masses and stellar velocity dispersions (60 < sigma < 350 km/s).

Abridged: We use three-dimensional SPH simulations to investigate the collapse of low-mass prestellar cores and the formation and early evolution of protostellar discs. The initial conditions are slightly supercritical Bonnor-Ebert spheres in rigid rotation. The core mass and initial radius are held fixed at M_O=6.1 M_sun and R_O=17,000 AU, and the only parameter that we vary is the initial angular speed \Omega_O. Protostellar discs forming from cores with \Omega_O<1.35 10d-13 1/s have radii between 100 and 300 AU and are quite centrally concentrated; due to heating by gas infall onto the disc and accretion onto the central object, they are also quite warm, T>100 K, and therefore stable against gravitational fragmentation. In contrast, more rapidly rotating cores form discs which are less concentrated and cooler, and have radii between 400 and 1000 AU; as a consequence they are prone to gravitational fragmentation and the formation of multiple systems. We derive a criterion that predicts whether a rigidly rotating core having given M_O, R_O and \Omega_O will produce a protostellar disc which fragments whilst material is still infalling from the core envelope. We then apply this criterion to core samples for which M_O, R_O and \Omega_O have been estimated observationally. We conclude that the observed cores are stable against fragmentation at this stage, due to their low angular speeds and the heat delivered at the accretion shock where the infalling material hits the disc.

We continue our presentation of VINE. We begin with a description of relevant architectural properties of the serial and shared memory parallel computers on which VINE is intended to run, and describe their influences on the design of the code itself. We continue with a detailed description of a number of optimizations made to the layout of the particle data in memory and to our implementation of a binary tree used to access that data for use in gravitational force calculations and searches for SPH neighbor particles. We describe modifications to the code necessary to obtain forces efficiently from special purpose `GRAPE' hardware. We conclude with an extensive series of performance tests, which demonstrate that the code can be run efficiently and without modification in serial on small workstations or in parallel using OpenMP compiler directives on large scale, shared memory parallel machines. We analyze the effects of the code optimizations and estimate that they improve its overall performance by more than an order of magnitude over that obtained by many other tree codes. Scaled parallel performance of the gravity and SPH calculations, together the most costly components of most simulations, is nearly linear up to maximum machine sizes available to us (120 processors on an Origin~3000). At similar accuracy, performance of VINE, used in GRAPE-tree mode, is approximately a factor two slower than that of VINE, used in host-only mode. Optimizations of the GRAPE/host communications could improve the speed by as much as a factor of three, but have not yet been implemented in VINE.

We present a Fortran 95 code for simulating the evolution of astrophysical systems using particles to represent the underlying fluid flow. The code is designed to be versatile, flexible and extensible, with modular options that can be selected either at compile time or at run time. We include a number of general purpose modules describing a variety of physical processes commonly required in the astrophysical community. The code can be used as an N-body code to evolve a set of particles in two or three dimensions using either a Leapfrog or Runge-Kutta-Fehlberg integrator, with or without individual timesteps for each particle. Particles may interact gravitationally as $N$-body particles, and all or any subset may also interact hydrodynamically, using the Smoothed Particle Hydrodynamic (SPH) method. Massive point particles (`stars') which may accrete nearby SPH or $N$-body particles may also be included. The default free boundary conditions can be replaced by a module to include periodic boundaries. Cosmological expansion may also be included. An interface with special purpose `GRAPE' hardware may also be selected. If available, forces obtained from the GRAPE coprocessors may be transparently substituted for those obtained from the default tree based calculation. The code may be run without modification on single processors or in parallel using OpenMP compiler directives on large scale, shared memory parallel machines. In comparison to the Gadget-2 code of Springel 2005, the gravitational force calculation is $\approx 3.5 - 4.8$ times faster with VINE when run on 8 Itanium~2 processors in an SGI Altix, while producing nearly identical outcomes in our test problems. We present simulations of several test problems, including a merger simulation of two elliptical galaxies with 800000 particles.

We developed a new method to extract halo merger rates from the Millennium Simulation. First, by removing superfluous mergers that are artifacts of the FOF halo finder, we find a lower merger rate compared to previous work. The reductions are more significant at lower redshifts, lower halo masses, and minor mergers. Our approach agrees better with predictions from the EPS model. Second, we find that the FOF halo finder overestimates the halo mass by up to 50% for halos that are about to merge, which leads to an additional ~20% overestimate of the merger rate. Therefore, we define halo masses by including only gravitationally bound particles. We provide new best-fitting parameters for a global formula to account for these improvements. In addition, we extract the merger rate per progenitor halo, as well as per descendant halo. The former is the quantity that is related to observed galaxy merger fractions when they are measured via pair counting. At low mass/redshift the merger rate increases with mass and redshift. At high mass/redshift (for halos with masses a few times the "knee" of the mass function) these trends break down, and the merger rate per progenitor halo decreases with mass and increases only moderately with redshift. Defining the merger rate per progenitor halo also allows us to quantify the rate at which halos are being accreted onto larger halos. We provide an analytic formula to convert given merger rates per descendant halo into merger rates per progenitor halo. Finally, we perform a comparison between observed merger fractions and the fraction of halos that have undergone a major merger during the recent dynamical friction time, and find a fair agreement, within the large uncertainties of the observations. Our new halo merger trees are available at http://www.mpe.mpg.de/ir/MillenniumMergerTrees/.

Using a high resolution hydrodynamical cosmological simulation of the formation of a massive spheroidal galaxy we show that elliptical galaxies can be very compact and massive at high redshift in agreement with recent observations. Accretion of stripped in-falling stellar material increases the size of the system with time and the central concentration is reduced by dynamical friction of the surviving stellar cores. In a specific case of a spheroidal galaxy with a final stellar mass of $1.5 \times 10^{11} M_{\odot}$ we find that the effective radius $r_e$ increases from $0.7 \pm 0.2 \rm kpc$ at z = 3 to $r_e = 2.4 \pm 0.4 \rm kpc$ at z = 0 with a concomitant decrease in the effective density of an order of magnitude and a decrease of the central velocity dispersion by approximately 20% over this time interval. A simple argument based on the virial theorem shows that during the accretion of weakly bound material (minor mergers) the radius can increase as the square of the mass in contrast to the usual linear rate of increase for major mergers. By undergoing minor mergers compact high redshift spheroids can evolve into present-day systems with sizes and concentrations similar to observed local ellipticals. This indicates that minor mergers may be the main driver for the late evolution of sizes and densities of early-type galaxies.

We present a set of global, self-consistent N-body/SPH simulations of the dynamic evolution of galactic discs with gas and including magnetic fields. We have implemented a description to follow the evolution of magnetic fields with the ideal induction equation in the SPH part of the Vine code. Results from a direct implementation of the field equations are compared to a representation by Euler potentials, which pose a div(B)-free description, an constraint not fulfilled for the direct implementation. All simulations are compared to an implementation of magnetic fields in the Gadget code which includes also cleaning methods for div(B). Starting with a homogeneous seed field we find that by differential rotation and spiral structure formation of the disc the field is amplified by one order of magnitude within five rotation periods of the disc. The amplification is stronger for higher numerical resolution. Moreover, we find a tight connection of the magnetic field structure to the density pattern of the galaxy in our simulations, with the magnetic field lines being aligned with the developing spiral pattern of the gas. Our simulations clearly show the importance of non-axisymmetry for the evolution of the magnetic field.

We study the thermal formation history of four simulated galaxies that were shown in Naab et al. (2007) to reproduce a number of observed properties of elliptical galaxies. The temperature of the gas in the galaxies is steadily increasing with decreasing redshift, although much of the gas has a cooling time shorter than the Hubble time. The gas is being heated and kept hot by gravitational heating processes through the release of potential energy from infalling stellar clumps. The energy is dissipated in supersonic collisions of infalling gas lumps with the ambient gas and through the dynamical capturing of satellite systems causing gravitational wakes that transfer energy to the surrounding gas. Furthermore dynamical friction from the infalling clumps pushes out dark matter, lowering the central dark matter density by up to a factor of two from z=3 to z=0. In galaxies in which the late formation history (z<2) is dominated by minor merging and accretion the energy released (E~5x10^{59} ergs) from gravitational feedback is sufficient to form red and dead elliptical galaxies by z~1 even in the absence of supernova and AGN feedback.

We present the modeling of SINFONI integral field dynamics of 18 star forming galaxies at z ~ 2 from Halpha line emission. The galaxies are selected from the larger sample of the SINS survey, based on the prominence of ordered rotational motions with respect to more complex merger induced dynamics. The quality of the data allows us to carefully select systems with kinematics dominated by rotation, and to model the gas dynamics across the whole galaxy using suitable exponential disk models. We obtain a good correlation between the dynamical mass and the stellar mass, finding that large gas fractions Mgas~M*) are required to explain the difference between the two quantities. We use the derived stellar mass and maximum rotational velocity Vmax from the modeling to construct for the first time the stellar mass Tully-Fisher relation at z ~ 2.2. The relation obtained shows a slope similar to what is observed at lower redshift, but we detect an evolution of the zero point. We find that at z ~ 2.2 there is an offset in log(M*) for a given rotational velocity of 0.41+-0.11 with respect to the local Universe. This result is consistent with the predictions of the latest N-body/hydrodynamical simulations of disk formation and evolution, which invoke gas accretion onto the forming disk in filaments and cooling flows. This scenario is in agreement with other dynamical evidence from SINS, where gas accretion from the halo is required to reproduce the observed properties of a large fraction of the z ~ 2 galaxies.

We present high resolution simulations on the impact of ionizing radiation of massive O-stars on the surrounding turbulent interstellar medium (ISM). The simulations are performed with the newly developed software iVINE which combines ionization with smoothed particle hydrodynamics (SPH) and gravitational forces. We show that radiation from hot stars penetrates the ISM, efficiently heats cold low density gas and amplifies over-densities seeded by the initial turbulence. The formation of observed pillar-like structures in star forming regions (e.g. in M16) can be explained by this scenario. At the tip of the pillars gravitational collapse can be induced, eventually leading to the formation of low mass stars. Detailed analysis of the evolution of the turbulent spectra shows that UV-radiation of O-stars indeed provides an excellent mechanism to sustain and even drive turbulence in the parental molecular cloud.

We use oblate axisymmetric dynamical models including dark halos to determine the orbital structure of intermediate mass to massive Coma early-type galaxies. We find a large variety of orbital compositions. Averaged over all sample galaxies the unordered stellar kinetic energy in the azimuthal and the radial direction are of the same order, but they can differ by up to 40 percent in individual systems. In contrast, both for rotating and non-rotating galaxies the vertical kinetic energy is on average smaller than in the other two directions. This implies that even most of the rotating ellipticals are flattened by an anisotropy in the stellar velocity dispersions. Using three-integral axisymmetric toy models we show that flattening by stellar anisotropy maximises the entropy for a given density distribution. Collisionless disk merger remnants are radially anisotropic. The apparent lack of strong radial anisotropy in observed early-type galaxies implies that they may not have formed from mergers of disks unless the influence of dissipational processes was significant.

We present a three-dimensional, fully parallelized, efficient implementation of ionizing UV radiation for smoothed particle hydrodynamics (SPH) including self-gravity. Our method is based on the SPH/tree code VINE. We therefore call it iVINE (for Ionization + VINE). This approach allows detailed high-resolution studies of the effects of ionizing radiation from e.g. young massive stars on their turbulent parental molecular clouds. In this paper we describe the concept and the numerical implementation of the radiative transfer for a plain-parallel geometry and we discuss several test cases demonstrating the efficiency and accuracy of the new method. As a first application, we study the radiatively driven implosion of marginally stable molecular clouds at various distances of a strong UV source and show that they are driven into gravitational collapse. The resulting cores are very compact and dense exactly as it is observed in clustered environments. Our simulations indicate that the time of triggered collapse depends on the distance of the core from the UV source. Clouds closer to the source collapse several $10^5$ years earlier than more distant clouds. This effect can explain the observed age spread in OB associations where stars closer to the source are found to be younger. We discuss possible uncertainties in the observational derivation of shock front velocities due to early stripping of proto-stellar envelopes by ionizing radiation.

A popular formation scenario for giant elliptical galaxies proposes that they might have formed from binary mergers of disk galaxies. Difficulties with the scenario that emerged from earlier studies included providing the necessary stellar mass and metallicity, maintaining the tight color-magnitude relation and avoiding phase space limits. In this paper we revisit the issue and put constraints on the binary disc merger scenario based on the stellar populations of disc galaxies. We draw the following conclusions: Low redshift collisionless or gaseous mergers of present day Milky Way like disc galaxies do not form present day elliptical galaxies. Binary mergers of the progenitors of present day Milky Way like disc galaxies can have evolved into intermediate mass elliptical galaxies ($M < M_*$) if they have merged earlier than $\approx$ 3-4 Gyrs ago. Assuming that most present day disk galaxies formed in a similar way to the Milky Way model presented here, more massive giant ellipticals in general can not have formed from binary mergers of the progenitors of present day disc galaxies. A major reason for these conclusions is that the mass in metals of typical disk galaxy is approximately a factor of 4-8 smaller than the mass in metals of a typical early-type galaxy and this ratio grows to larger values with increasing redshift.

We use two-dimensional kinematic maps of simulated binary disc mergers to investigate the $\lambda_R$-parameter, which is a luminosity weighted measure of projected angular momentum per unit mass. This parameter was introduced to subdivide the SAURON sample of early type galaxies in so called fast $\lambda_R > 0.1$ and slow rotators $\lambda_R < 0.1$. Tests on merger remnants reveal that $\lambda_R$ is a robust indicator of the true angular momentum content in elliptical galaxies. We find the same range of $\lambda_R$ values in our merger remnants as in the SAURON galaxies. The merger mass ratio is decisive in creating a slow or a fast rotator in a single binary merger, the former being created mostly in an equal mass merger. Slow rotators have a $\lambda_R$ which does not vary with projection. The confusion rate with face-on fast rotators is very small. Merger with low gas fractions form slow rotators with smaller ellipticities and are in much better agreement with the SAURON slow rotators. Remergers of merger remnants are slow rotators but tend to have too high ellipticities. Fast rotators maintain the angular momentum content from the progenitor disc galaxy if merger mass ratio is high. Some SAURON galaxies have values of $\lambda_R$ as high as our progenitor disc galaxies.

In the framework of hierarchical structure formation ellipticals can form from merging of smaller disk galaxies. The nearby interacting 'Antennae' galaxy pair (NGC 4038/39) is one of the best-studied local systems of merging spirals, thus presenting us with an ideal laboratory for the study of galaxy evolution models. The Antennae are believed to be in a state prior to their final encounter with rapid subsequent merging, which puts them in the first position of the Toomre (1977) merger sequence. Here we present first numerical high-resolution, self-consistent, smoothed particle hydrodynamics (SPH) simulations of the Antennae system, including star formation and stellar feedback, and compare our results to VLA HI observations by Hibbard et al. (2001). We are able to obtain a close, but not yet perfect match to the observed morphology and kinematics of the system.

We present binary galaxy merger simulations of gas-rich disks (Sp-Sp), of early-type galaxies and disks (E-Sp, mixed mergers), and mergers of early-type galaxies (E-E, dry mergers) with varying mass ratios and different progenitor morphologies. The simulations include radiative cooling, star formation and black hole (BH) accretion and the associated feedback processes. We find for Sp-Sp mergers, that the peak star formation rate and BH accretion rate decrease and the growth timescales of the central black holes and newly formed stars increase with higher progenitor mass ratios. The termination of star formation by BH feedback in disk mergers is significantly less important for higher progenitor mass ratios (e.g. 3:1 and higher). In addition, the inclusion of BH feedback suppresses efficiently star formation in dry E-E mergers and mixed E-Sp mergers.

Gravitational tides are widely understood to strip and destroy galactic substructures. In the course of a galaxy merger, however, transient totally compressive tides may develop and prevent star forming regions from dissolving, after they condensed to form clusters of stars. We study the statistics of such compressive modes in an N-body model of the galaxy merger NGC 4038/39 (the Antennae) and show that ~15% of the disc material undergoes compressive tides at pericentre. The spatial distribution of observed young clusters in the overlap and nuclear regions of the Antennae matches surprisingly well the location of compressive tides obtained from simulation data. Furthermore, the statistics of time intervals spent by individual particles embedded in a compressive tide yields a log-normal distribution of characteristic time ~10 Myr, comparable to star cluster formation timescales. We argue that this generic process is operative in galaxy mergers at all redshifts and possibly enhances the formation of star clusters. We show with a model calculation that this process will prevent the dissolution of a star cluster during the formation phase, even for a star formation efficiency as low as ~10%. The transient nature of compressive tides implies that clusters may dissolve rapidly once the tidal field switches to the usual disruptive mode.

We present binary galaxy merger simulations with varying mass ratios and different progenitor morphologies. The simulations include mergers of gas-rich disks (Sp-Sp), of early-type galaxies and disks (E-Sp, mixed mergers), and mergers of early-type galaxies (E-E, dry mergers). We follow the dynamics of gas, stars and dark matter, and include radiative cooling, star formation and black hole (BH) accretion. For Sp-Sp mergers, the peak star formation rate and BH accretion rate decrease and the growth timescales of the central black holes and newly formed stars increase with higher progenitor mass ratios. The peak BH accretion rate typically occurs shortly after the time of BH merging for low progenitor mass ratios (e.g. 3:1 and lower), whereas for higher progenitor mass ratios there is no clear correlation between the peak BH accretion rate and BH merging time. The termination of star formation by BH feedback in disk mergers is significantly less important for higher progenitor mass ratios (e.g. 3:1 and higher). In addition, the inclusion of BH feedback suppresses efficiently star formation in dry E-E mergers and mixed E-Sp mergers. All merger remnants, independent of their progenitors, follow the observed relations between the central BH mass and the stellar velocity dispersion M_BH-sigma, the bulge mass M_BH-M_* and the bulge binding energy M_BH-M_{*}sigma^2, with the dominant source of scatter arising from variations in the initial gas mass fraction. The normalizations for all relations and the simulated slope of the M_BH-sigma and M_BH-M_{*}sigma^2 relations are in good agreement with the observations, whereas the simulated slope of the M_BH-M_* relation is slightly steeper compared to the observations. (abridged)

We examine the origin of photometrical properties of N-body merger remnants from the perpective of their orbital content. We show that disc mergers alone are unlikely to form a boxy-discy dichotomy in isophotal shape, because of the survival of disc-like orbits even in violent mergers. The shape of the different orbit families can vary strongly, depending on how violent the merger was. Minor axis tubes can become boxy and box orbits can be round. However, for edge-on projections isophotal shape, orbital content and line-of-sight velocity are well connected.

Recent observations of UV-/optically selected, massive star forming galaxies at z~2 indicate that the baryonic mass assembly and star formation history is dominated by continuous rapid accretion of gas and internal secular evolution, rather than by major mergers. We use the Millennium Simulation to build new halo merger trees, and extract halo merger fractions and mass accretion rates. We find that even for halos not undergoing major mergers the mass accretion rates are plausibly sufficient to account for the high star formation rates observed in z~2 disks. On the other hand, the fraction of major mergers in the Millennium Simulation is sufficient to account for the number counts of submillimeter galaxies (SMGs), in support of observational evidence that these are major mergers. When following the fate of these two populations in the Millennium Simulation to z=0, we find that subsequent mergers are not frequent enough to convert all z~2 turbulent disks into elliptical galaxies at z=0. Similarly, mergers cannot transform the compact SMGs/red sequence galaxies at z~2 into observed massive cluster ellipticals at z=0. We argue therefore, that secular and internal evolution must play an important role in the evolution of a significant fraction of z~2 UV-/optically and submillimeter selected galaxy populations.

We present Ha integral field spectroscopy of well resolved, UV/optically selected z~2 star-forming galaxies as part of the SINS survey with SINFONI on the ESO VLT. Our laser guide star adaptive optics and good seeing data show the presence of turbulent rotating star forming rings/disks, plus central bulge/inner disk components, whose mass fractions relative to total dynamical mass appears to scale with [NII]/Ha flux ratio and star formation age. We propose that the buildup of the central disks and bulges of massive galaxies at z~2 can be driven by the early secular evolution of gas-rich proto-disks. High redshift disks exhibit large random motions. This turbulence may in part be stirred up by the release of gravitational energy in the rapid cold accretion flows along the filaments of the cosmic web. As a result dynamical friction and viscous processes proceed on a time scale of <1 Gyr, at least an order of magnitude faster than in z~0 disk galaxies. Early secular evolution thus drives gas and stars into the central regions and can build up exponential disks and massive bulges, even without major mergers. Secular evolution along with increased efficiency of star formation at high surface densities may also help to account for the short time scales of the stellar buildup observed in massive galaxies at z~2.

We derive new constraints on the mass of the Milky Way's dark matter halo, based on a set of halo stars from SDSS as kinematic tracers. Our sample comprises 2401 rigorously selected Blue Horizontal-Branch (BHB) halo stars drawn from SDSS DR-6. To interpret these distributions, we compare them to matched mock observations drawn from two different cosmological galaxy formation simulations designed to resemble the Milky Way, which we presume to have an appropriate orbital distribution of halo stars. We then determine which value of $\rm V_{cir}(r)$ brings the observed distribution into agreement with the corresponding distributions from the simulations. This procedure results in an estimate of the Milky Way's circular velocity curve to $\sim 60$ kpc, which is found to be slightly falling from the adopted value of $\rm 220 km s^{-1}$ at the Sun's location, and implies M$(<60 \rm kpc) = 4.0\pm 0.7\times 10^{11}$M$_\odot$. The radial dependence of $\rm V_{cir}(r)$, derived in statistically independent bins, is found to be consistent with the expectations from an NFW dark matter halo with the established stellar mass components at its center. If we assume an NFW halo profile of characteristic concentration holds, we can use the observations to estimate the virial mass of the Milky Way's dark matter halo, M$_{\rm vir}=1.0^{+0.3}_{-0.2} \times 10^{12}$M$_\odot$, which is lower than many previous estimates. This estimate implies that nearly 40% of the baryons within the virial radius of the Milky Way's dark matter halo reside in the stellar components of our Galaxy. A value for M$_{\rm vir}$ of only $\sim 1\times10^{12}$M$_\odot$ also (re-)opens the question of whether all of the Milky Way's satellite galaxies are on bound orbits.

We present a simple set of kinematic criteria that can distinguish between galaxies dominated by ordered rotational motion and those involved in major merger events. Our criteria are based on the dynamics of the warm ionized gas (as traced by H-alpha) within galaxies, making this analysis accessible to high-redshift systems, whose kinematics are primarily traceable through emission features. Using the method of kinemetry (developed by Krajnovic and co-workers), we quantify asymmetries in both the velocity and velocity dispersion maps of the warm gas, and the resulting criteria enable us to empirically differentiate between non-merging and merging systems at high redshift. We apply these criteria to 11 of our best-studied rest-frame UV/optical-selected z~2 galaxies for which we have near infrared integral field spectroscopic data from SINFONI on the VLT. Of these 11 systems, we find that >50% have kinematics consistent with a single rotating disk interpretation, while the remaining systems are more likely undergoing major mergers. This result, combined with the short formation timescales of these systems, provides evidence that rapid, smooth accretion of gas plays a significant role in galaxy formation at high redshift.

The intrinsic anisotropy delta and flattening epsilon of simulated merger remnants is compared with elliptical galaxies that have been observed by the SAURON collaboration, and that were analysed using axisymmetric Schwarzschild models. Collisionless binary mergers of stellar disks and disk mergers with an additional isothermal gas component, neglecting star formation cannot reproduce the observed trend delta = 0.55 epsilon (SAURON relationship). An excellent fit of the SAURON relationship for flattened ellipticals with epsilon >= 0.25 is however found for merger simulations of disks with gas fractions >= 20%, including star formation and stellar energy feedback. Massive black hole feedback does not strongly affect this result. Subsequent dry merging of merger remnants however does not generate the slowly-rotating SAURON ellipticals which are characterized by low ellipticities epsilon < 0.25 and low anisotropies. This indicates that at least some ellipticals on the red galaxy sequence did not form by binary mergers of disks or early-type galaxies. We show that stellar spheroids resulting from multiple, hierarchical mergers of star-bursting subunits in a cosmological context are in excellent agreement with the low ellipticities and anisotropies of the slowly rotating SAURON ellipticals and their observed trend of delta with epsilon. The numerical simulations indicate that the SAURON relation might be a result of strong violent relaxation and phase mixing of multiple, kinematically cold stellar subunits with the angular momentum of the system determining its location on the relation.

Galactic spheroids can form as a result of galaxy interactions and mergers of disks. Detailed analyses of the photometric properties, the intrinsic orbital structure, the line-of-sight velocity distributions and the kinemetry of simulated merger remnants, which depend critically on the geometry and the gas content of the interacting progenitors, indicate that low and intermediate mass rotating ellipticals can form from mergers of disks. The masses and metallicities of all massive ellipticals and the kinematics of some massive non-rotating ellipticals cannot be explained by binary mergers. Thus these galaxies might have formed in a different way.

We model mock observations of collisionless N-body disc-disc mergers with an axisymmetric orbit superposition program that has been used to model Coma ellipticals. The remnants sample representatively the shapes of disc-disc mergers including prolate, triaxial and oblate objects. The aim is to better understand how the assumption of axial symmetry affects reconstructed masses and stellar motions of systems which are intrinsically not axisymmetric, whether it leads to a bias and how such a potential bias can be recognised in models of real galaxies. The mass recovery at the half-light radius depends on viewing-angle and intrinsic shape: edge-on views allow to reconstruct total masses with an accuracy between 20% (triaxial/prolate remnants) and 3% (oblate remnant). Masses of highly flattened, face-on systems are underestimated by up to 50%. Deviations in local mass densities can be larger where remnants are strongly triaxial or prolate. Luminous M/L are sensitive to box orbits in the remnants. Box orbits cause the central kinematics to vary with viewing-angle. Reconstructed luminous M/L and central masses follow this variation. Luminous M/L are always underestimated (up to a factor of 2.5). Respective dark halos in the models can be overestimated by about the same amount, depending again on viewing angle. Reconstructed velocity anisotropies depend on viewing angle and on the orbital composition of the remnant. We construct N-body realisations of the Schwarzschild models to discuss chaotic orbits and the virial equilibrium in our models. Apparently flattened, rotating ellipticals of intermediate mass are likely close to both, axial symmetry and edge-on orientation. Our results imply that Schwarzschild models allow a reconstruction of their masses and stellar anisotropies with high accuracy. (abridged)

We present a two-dimensional kinematic analysis for a sample of simulated binary disc merger remnants with mass ratios 1:1 and 3:1. For the progenitor discs we used pure stellar models as well as models with 10% of their mass in gas. A multitude of phenomena also observed in real galaxies are found in the simulations. These include misaligned rotation, embedded discs, gas rings, counter-rotating cores and kinematic misaligned discs. Using the 2D maps we illustrate projection effects and the change in properties of a merger remnant when gas is included in the merger. We find that kinematic peculiar subsystems are preferably formed in equal mass mergers. Equal-mass collisionless remnants can show almost no rotation, regular rotation or strong kinematic misalignment. The inclusion of gas makes the remnants appear more round(1:1) and axisymmetric(3:1). Counter-Rotating Cores (CRCs) are almost exclusively formed in equal-mass mergers with a dissipational component. 3:1 remnants show a much more regular structure. We quantify these properties by applying the kinemetric methods recently developed by Krajnovi\'c et al. This work will help to understand observations of elliptical galaxies with 2D field spectrographs, like SAURON.

The formation of tidal dwarf galaxies (TDG) inside tidal arms of interacting disk galaxies has been studied with N-body and N-body/SPH simulations at different resolutions. In pure N-body simulations no bound objects are formed at high resolution. At low resolution bound objects can form in tidal tails in agreement with previous work. We conclude that tidal dwarf galaxies are not likely to form by pure collisionless collapse in tidal tails. However, the presence of a sufficiently massive and extended gas component in the progenitor disk supports the formation of bound stellar objects in the tidal arms. Our results clearly favor a dissipation supported scenario in which the formation of TDGs is induced by the local collapse of gas which then triggers the collapse of the stellar component.

We describe high resolution Smoothed Particle Hydrodynamics (SPH) simulations of three approximately $M_*$ field galaxies starting from \LCDM initial conditions. The simulations are made intentionally simple, and include photoionization, cooling of the intergalactic medium, and star formation but not feedback from AGN or supernovae. All of the galaxies undergo an initial burst of star formation at $z \approx 5$, accompanied by the formation of a bubble of heated gas. Two out of three galaxies show early-type properties at present whereas only one of them experienced a major merger. Heating from shocks and -PdV work dominates over cooling so that for most of the gas the temperature is an increasing function of time. By $z \approx 1$ a significant fraction of the final stellar mass is in place and the spectral energy distribution resembles those of observed massive red galaxies. The galaxies have grown from $z=1 \to 0$ on average by 25% in mass and in size by gas poor (dry) stellar mergers. By the present day, the simulated galaxies are old ($\approx 10 {\rm Gyrs}$), kinematically hot stellar systems surrounded by hot gaseous haloes. Stars dominate the mass of the galaxies up to $\approx 4$ effective radii ($\approx 10$ kpc). Kinematic and most photometric properties are in good agreement with those of observed elliptical galaxies. The galaxy with a major merger develops a counter-rotating core. Our simulations show that realistic intermediate mass giant elliptical galaxies with plausible formation histories can be formed from \LCDM initial conditions even without requiring recent major mergers or feedback from supernovae or AGN.

Recent observations with the Spitzer Space Telescope show clear evidence that star formation takes place in the surrounding of young massive O-type stars, which are shaping their environment due to their powerful radiation and stellar winds. In this work we investigate the effect of ionising radiation of massive stars on the ambient interstellar medium (ISM): In particular we want to examine whether the UV-radiation of O-type stars can lead to the observed pillar-like structures and can trigger star formation. We developed a new implementation, based on a parallel Smooth Particle Hydrodynamics code (called IVINE), that allows an efficient treatment of the effect of ionising radiation from massive stars on their turbulent gaseous environment. Here we present first results at very high resolution. We show that ionising radiation can trigger the collapse of an otherwise stable molecular cloud. The arising structures resemble observed structures (e.g. the pillars of creation in the Eagle Nebula (M16) or the Horsehead Nebula B33). Including the effect of gravitation we find small regions that can be identified as formation places of individual stars. We conclude that ionising radiation from massive stars alone can trigger substantial star formation in molecular clouds.

We present a large set of merger simulations of early-type disc galaxies with mass ratios of 1:1 and 3:1 and 10% of the total disc mass in gas. In contrast to the collisionless case equal-mass mergers with gas do not result in very boxy remnants which is caused by the suppression of box orbits and the change of the projected shape of minor-axis tube orbits in the more axisymmetric remnants. The isophotal shape of 3:1 remnants and the global kinematic properties of 1:1 and 3:1 remnants are only weakly affected by the presence of gas. 1:1 remnants are slowly rotating whereas 3:1 remnants are fast rotating and discy. The shape of the stellar LOSVD is strongly influenced by gas. The LOSVDs of collisionless remnants have broad leading wings while their gaseous counterparts show steep leading wings, more consistent with observations of elliptical galaxies. We show that this change is also caused by the suppressed populating of box orbits and it is amplified by the formation of extended gas discs in the merger remnants. If elliptical galaxies have formed from mergers our results indicate that massive, slowly rotating boxy elliptical galaxies can not have formed from dissipative mergers of discs. Pure stellar (dry) mergers are the more likely candidates. On the other hand lower mass, fast rotating and discy ellipticals can have formed from dissipative (wet) mergers of early-type discs. So far, only unequal-mass disc mergers with gas can successfully explain their observed substructure. This is consistent with the revised morphological classification scheme of increasing importance of gas dissipation when moving from boxy ellipticals to discy ellipticals and then to spiral galaxies, proposed by Kormendy & Bender (abbreviated).

We present results from our Very Large Telescope large program to study the dynamical evolution of local Ultraluminous Infrared Galaxies (ULIRGs) and QSOs. This paper is the second in a series presenting the stellar kinematics of 54 ULIRGs, derived from high resolution, long-slit H- and K-band spectroscopy. The data presented here, including observations of 17 new targets, are mainly focused on sources that have coalesced into a single nucleus. The stellar kinematics, extracted from the CO ro-vibrational bandheads in our spectra, indicate that ULIRG remnants are dynamically heated systems with a mean dispersion of 161 km/s. The combination of kinematic, structural, and photometric properties of the remnants indicate that they mostly originate from major mergers and that they result in the formation of systems supported by random motions, therefore, elliptical galaxies. The peak of the velocity dispersion distribution and the locus of ULIRGs on the fundamental plane of early-type galaxies indicate that the end products of ultraluminous mergers are typically moderate-mass ellipticals (of stellar mass ~10^10 - 10^11 M_sun). Converting the host dispersion into black hole mass with the aid of the M_BH-sigma relation yields black hole mass estimates of the order 10^7 - 10^8 M_sun and high accretion rates with Eddington efficiencies often >0.5.

We present a detailed surface density analysis of a large sample of simulated collisionless mergers of disc galaxies with bulges (mass ratios 1:1, 2:1, 3:1, 4:1, and 6:1) and without bulges (mass ratios 1:1 and 3:1). They were classified, according to their bulge--to--total ($B/T$) ratio, either as a one-component system or as a two-component systems. In general projection effects change the classification of a remnant. Only merger remnants of discs with bulges show properties similar to observed early--type galaxies. Their B/T ratios are in the range $0.2<B/T<1$. Surprisingly, the initial mass ratio has a weak influence on the distributions of $B/T$, effective radius and S\'ersic index $n$. For all one-component projections ($\approx 60%$ of all projections) the S\'ersic index distribution peaks at $3<n<4$. However, the mass ratio is tightly linked to the properties of the outer exponential components which resemble pressure supported, spheroidal halos for 1:1 and 2:1 remnants and elongated heated discs for 6:1 remnants. We found distinct correlations between the fitting parameters which are very similar to observed relations (e.g. larger bulges have lower effective surface densities). No indications for a correlation between the surface density profiles and other global parameters are found. The remnants have properties similar to giant elliptical galaxies in the intermediate mass regime. A binary disc merger origin for all early-type galaxies, especially the most massive ones, is unlikely. Observed nearby merger remnants have properties similar to the simulated remnants. They can have formed from binary disc mergers and might evolve into early--type galaxies within a few Gyrs (abbreviated).

We present a study of the dynamical evolution of Ultraluminous Infrared Galaxies (ULIRGs), merging galaxies of infrared luminosity >10^12 L_sun. During our Very Large Telescope large program, we have obtained ISAAC near-infrared, high-resolution spectra of 54 ULIRGs (at several merger phases) and 12 local Palomar-Green QSOs to investigate whether ULIRGs go through a QSO phase during their evolution. One possible evolutionary scenario is that after nuclear coalescence, the black hole radiates close to Eddington to produce QSO luminosities. The mean stellar velocity dispersion that we measure from our spectra is similar (~160 km/s) for 30 post-coalescence ULIRGs and 7 IR-bright QSOs. The black holes in both populations have masses of order 10^7-10^8 M_sun (calculated from the relation to the host dispersion) and accrete at rates >0.5 Eddington. Placing ULIRGs and IR-bright QSOs on the fundamental plane of early-type galaxies shows that they are located on a similar region (that of moderate-mass ellipticals), in contrast to giant ellipticals and radio-loud QSOs. While this preliminary comparison of the ULIRG and QSO host kinematical properties indicates that (some) ULIRGs may undergo a QSO phase in their evolutionary history before they settle down as ellipticals, further data on non-IR excess QSOs are necessary to test this scenario.

The luminosity dependence of kinematical and isophotal properties of elliptical galaxies is investigated using numerical simulations of galaxy merging, combined with semi-analytical models of hierarchical structure formation. Mergers of spiral galaxies as the only formation mechanism for elliptical galaxies can neither reproduce the kinematical and photometric properties of very massive elliptical galaxies nor the change from rotationally flattened disky to anisotropic boxy systems with increasing luminosity. We present numerical simulations showing that binary mergers of early-type galaxies open an additional channel for the formation of anisotropic, slowly rotating and boxy ellipticals. Including this channel in a semi-analytical model we can successfully reproduce the observed trend that more luminous giant ellipticals are more boxy and less flattened by rotation. This trend can be strengthened by suppressing residual gas infall and star formation for galaxies with stellar bulge masses $M_* \geq 3 \times 10^{10} M_{\odot}$. Hence we propose that mergers of early-type galaxies play an important role for the assembly of massive elliptical galaxies.

We have used the 28'x 28' HST image mosaic from the GEMS (Galaxy Evolution from Morphology and SEDs) survey in conjunction with the COMBO-17 photometric redshift survey to constrain the incidence of major mergers between spheroid-dominated galaxies with little cold gas (dry mergers) since z = 0.7. A set of N-body merger simulations was used to explore the morphological signatures of such interactions: they are recognizable either as < 5kpc separation close pairs or because of broad, low surface brightness tidal features and asymmetries. Data with the depth and resolution of GEMS are sensitive to dry mergers between galaxies with M_V < -20.5 for z < 0.7; dry mergers at higher redshifts are not easily recovered in single-orbit HST imaging. Six dry mergers (12 galaxies) with luminosity ratios between 1:1 and 4:1 were found from a sample of 379 red early-type galaxies with M_V < -20.5 and 0.1 < z < 0.7. The simulations suggest that the morphological signatures of dry merging are visible for ~250Myr and we use this timescale to convert the observed merger incidence into a rate. On this basis we find that present day spheroidal galaxies with M_V < -20.5 on average have undergone between 0.5 and 2 major dry mergers since z ~ 0.7. We have compared this result with the predictions of a Cold Dark Matter based semi-analytic galaxy formation model. The model reproduces the observed declining major merger fraction of bright galaxies and the space density of luminous early-type galaxies reasonably well. The predicted dry merger fraction is consistent with our observational result. Hence, hierarchical models predict and observations now show that major dry mergers are an important driver of the evolution of massive early-type galaxies in recent epochs.

A simple model for the evolution of disc galaxies is presented. We adopt three numbers from observations of the Milky Way disc, the local surface mass density, the stellar scale length (of the assumedly exponential disc) and the amplitude of the (assumedly flat) rotation curve, and physically, the (local) dynamical Kennicutt star formation prescription, standard chemical evolution equations assuming and a model for spectral evolution of stellar populations. We can determine the detailed evolution of the model with only the addition of standard cosmological scalings with time of the dimensional parameters. A surprising wealth of detailed specifications follows from this prescription including the gaseous infall rate as a function of radius and time, the distribution of stellar ages and metallicities with time and radius, surface brightness profiles at different wavelengths, colours etc. At the solar neighbourhood stars start to form $\approx 10 Gyrs$ ago at an increasing rate peaking 4 billion years ago and then slowly declining in good agreement with observations. The mean age of long lived stars at the solar neighbourhood is about $4 Gyrs$. The local surface density of the stars and gas are 35 and $15 M_{\odot}pc^{-2}$, respectively. The metallicity distribution of the stars at the solar radius is narrow with a peak at $[Z/Z_{\odot}] = -0.1$.Both a Salpeter IMF and a Chabrier IMF are consistent with observations. Comparisons with the current and local fossil evidence provides support for the model which can then be used to assess other local disc galaxies, the evolution of disc galaxies in deep optical surveys and also for theoretical investigations such as simulations of merging disc galaxies (abbreviated).

We present first results from our Very Large Telescope large program to study the dynamical evolution of Ultraluminous Infrared Galaxies (ULIRGs), which are the products of mergers of gas-rich galaxies. The full data set consists of high resolution, long-slit, H- and K-band spectra of 38 ULIRGs and 12 QSOs (between 0.042<z<0.268). In this paper, we present the sources that have not fully coalesced, and therefore have two distinct nuclei. This sub-sample consists of 21 ULIRGs, the nuclear separation of which varies between 1.6 and 23.3 kpc. From the CO bandheads that appear in our spectra, we extract the stellar velocity dispersion, sigma, and the rotational velocity, V_rot. The stellar dispersion equals 142 km/s on average, while V_rot is often of the same order. We combine our spectroscopic results with high-resolution infrared (IR) imaging data to study the conditions for ULIRG activity in interacting pairs. We find that the majority of ULIRGs are triggered by almost equal-mass major mergers of 1.5:1 average ratio. Less frequently, 3:1 encounters are also observed in our sample. However, less violent mergers of mass ratio >3:1 typically do not force enough gas into the center to generate ULIRG luminosities.

We present results from our VLT Large Program to study the dynamical evolution of Ultraluminous Infrared Galaxies (ULIRGs). We have so far obtained near-infrared high-resolution ISAAC spectra of 53 local ULIRGs at several merger timescales and 12 Palomar-Green QSOs (more than half of which are IR-bright sources). We have extracted the stellar velocity dispersion and rotational velocity along our slits to derive the kinematics of the merging galaxies. These quantities enable us to answer the following questions about the evolution of ULIRGs: 1) What are the progenitor mass ratios?, 2) How do the stellar kinematics evolve?, and, 3) Is there a connection between ULIRGs and QSOs? We find that the Ultraluminous phase is mainly triggered by mergers of approximately equal mass galaxies, however, less violent minor mergers (of progenitor mass ratio ~3:1) also exist in our sample. Dynamical heating of the merging hosts is observed as the stellar systematic rotation decreases with time in favour of the increase of random motions. The merger remnants, being dispersion-dominated systems with non-negligible rotation, resemble elliptical galaxies. Placing ULIRGs on the fundamental plane of early-type galaxies shows that they resemble intermediate-mass ellipticals. After the nuclear coalescence, the black hole masses of ULIRGs, calculated from their relation to the host dispersions, are of the order 10^7-10^8 M_sun. To investigate whether ULIRGs go through a QSO phase during their evolution, we perform a similar (preliminary) analysis of the IR-bright-QSO kinematics. We find that the average dispersion of the IR-bright QSOs is similar to that of the ULIRG remnants, indicating that evolutionary links between the two populations may exist.

We present a detailed investigation of the relation between the orbital content of merger remnants and observable properties of elliptical and S0 galaxies. Our analysis is based on the statistical sample of collisionless mergers of disk galaxies with different mass ratios and orbital parameters, published by Naab & Burkert. We use the spectral method by Carpintero & Aguilar to determine the orbital content of every remnant and correlate it with its intrinsic shape, and its projected kinematic and photometric properties. We discuss the influence of the bulge component and varying pericenter distances. The two most abundant orbit classes are the minor axis tubes and the box orbits. Their ratio seems to determine the basic properties of a remnant. On average, the fraction of minor axis tubes increases by a factor of two from a merger mass ratio of 1:1 to 4:1, whereas the fraction of box orbits decreases by 10%. At a given mass the central velocity dispersion of a remnant scales with the ratio of minor axis tubes to box orbits. Interestingly, the division line between rotational supported systems and pressure supported objects, $(v_{maj}/\sigma_0)^*=0.7$, turns out to coincide with a box to minor axis tube ratio of unity. The observed $h_3$-$v/\sigma$ anti-correlation for ellipticals can not be reproduced by collisionless merger remnants. We propose that this can only be reconciled by an additional physical process that significantly reduces the box orbit content. Remnants which are dominated by minor axis tube orbits have predominantly disky projections. Boxy remnants have always a box to minor axis tube ratio larger than one. This study will enable to identify observed ellipticals that could have formed, in the collisionless limit, by gas-poor disk mergers.

The projected kinematical properties of unequal-mass merger remnants of disk galaxies are analysed and shown to agree well with observations of disky, fast rotating elliptical galaxies. This supports the major merger hypothesis of early-type galaxy formation. However, in contrast to previous claims, the merger remnants are very anisotropic with values of the anisotropy parameter that are similar to equal-mass merger remnants that form boxy, slowly rotating ellipticals. Including gas in the simulations does not change this result although the line-of-sight velocity profile and the intrinsic orbital structure are strongly affected by the presenceof gas. The kinematical difference between boxy and disky ellipticals appears not to be the amount of anisotropy but rather rotation and the shape of the velocity dispersion tensor. The apparent isotropy of observed disky ellipticals is shown to result from inclination effects. Even small inclination angles strongly reduce the measured anisotropy of fast rotating systems, seen in projection. A second problem is the limited amount of information that is available when measuring only the central velocity dispersion and a characteristic rotation and ellipticity. Methods are investigated that allow a better determination of the intrinsic anisotropy of fast rotating early-type galaxies with known inclination angles.

We perform a large parameter survey of 112 collisionless N-body simulations of binary mergers of disk-galaxies with mass ratios of 1:1, 2:1, 3:1, and 4:1 using the special purpose hardware GRAPE. The basic properties of the remnants correlate with the mass ratio of the progenitor disks. We find that about 80% of the equal-mass merger simulations lead to slowly rotating remnants with significant minor-axis rotation. One half of all projected 1:1 remnants shows boxy and the other half shows disky isophotes. Remnants with mass ratios of 3:1 and 4:1 all rotate fast and show a small amount of minor-axis rotation. About 90% of the projected 3:1 and 4:1 remnants show disky isophotes. 2:1 remnants show intermediate properties. Projection effects lead to a large spread in the data in good agreement with observations. They do not change the fundamental kinematical differences between equal-mass and unequal-mass merger remnants. The results of this study weaken the disk merger scenario as the possible formation mechanism of massive boxy giant ellipticals. In general, the properties of equal-mass (and 2:1) merger remnants are consistent with the observed population of giant ellipticals in the intermediate mass regime between low mass fast rotating disky and bright massive boxy giant ellipticals. 3:1 and 4:1 merger remnants, however, are in very good agreement with the class of low luminosity, fast rotating giant elliptical galaxies. Binary mergers of disk galaxies are therefore still very good candidates for being the main formation mechanism for intermediate and low mass giant ellipticals. The homogeneous class of massive boxy ellipticals most likely formed by a different process.

We present the results of a series of gas dynamical cosmological simulations of the formation of individual massive field galaxies in the standard concordance \LCDM and in a \LWDM cosmology. The galaxies form in an initial burst of star formation followed by accretion of small satellites. The simulated galaxies are old ($\approx$ 10 $Gyrs$) hot stellar systems with masses of $\approx 1.7 \times 10^{11} M_{\odot}$. Baryonic matter dominates the mass in the luminous part of the galaxies up to $\approx 5$ effective radii. The \LCDM galaxy is a slowly rotating spheroidal stellar system (E2) with predominantly disky isophotes. The line-of-sight velocity distributions deviate from Gaussian shape and $h_3$ is anticorrelated with $v_{los}$. The corresponding \LWDM galaxy is more elongated (E3 - E4) and rotates faster. The anisotropy parameter $(v/\sigma)^*$ is close to unity. Within the uncertainties of $M/L$ both galaxies follow the Fundamental Plane. We conclude that the properties of the two galaxies simulated in the \LCDM and \LWDM cosmology are in good agreement with observations of intermediate mass elliptical or S0 galaxies (abstract abbreviated).

We summarize current models of the formation of spheroidal stellar systems. Whereas globular clusters form in an efficient mode of star formation inside turbulent molecular clouds, the origin of galactic spheroids, that is bulges, dwarf ellipticals, and giant ellipticals, is directly coupled with structure formation and merging of structures in the Universe. Disks are the fundamental building blocks of galaxies and the progenitors of galactic spheroids. The origin of the various types of spheroids and their global properties can be understood as a result of disk heating by external perturbations, internal disk instabilities, or minor and major mergers.

The formation of elliptical galaxies as a result of the merging of spiral galaxies is discussed. We analyse a large set of numerical N-Body merger simulations which show that major mergers can in principle explain the observed isophotal fine structure of ellipticals and its correlation with kinematical properties. Equal-mass mergers lead to boxy, slowly rotating systems, unequal-mass mergers produce fast rotating and disky ellipticals. However, several problems remain. Anisotropic equal mass mergers appear disky under certain projections which is not observed. The intrinsic ellipticities of remnants are often larger than observed. Finally, although unequal-mass mergers produce fast rotating ellipticals, the remnants are in general more anisotropic than expected from observations. Additional processes seem to play an important role which are not included in dissipationless mergers. Resolving these problems might provide interesting new information on the structure and gas content of the progenitors of early-type galaxies.

We use high resolution numerical simulations to investigate the adiabatic contraction of dark matter halos with a Hernquist density profile. We test the response of the halos to the growth of additional axisymmetric disk potentials with various central concentrations and the spherically symmetric potential of a softened point mass. Adding the potentials on timescales that are long compared to the dynamical time scale of the halo, the contracted halos have density profiles that are in excellent agreement with analytical predictions based on the conservation of the adiabatic invariant $M(r)r$. This is surprising as this quantity is strictly conserved only for particles on circular orbits and in spherically symmetric potentials. If the same potentials are added on timescales that are short compared to the dynamical timescale, the result depends strongly on the adopted potential. The adiabatic approximation still works for disk potentials. It does, however, fail for the central potential.

Challenging stellar dynamical problems, such as the study of gravothermal oscillations in star clusters, have in the past initiated the very successful building of GRAPE special purpose computers. It is discussed, that present day tasks such as the formation and evolution of galactic nuclei with one or more massive black holes and the coupled stellar and gas dynamical processes in the formation of nuclei and star clusters, demand a new kind of hybrid architecture, using both GRAPE and a reconfigurable logics board called RACE. For such a system we have developed first implementations and floating point performance studies in the case of the SPH algorithm (smoothed particle hydrodynamics), which will be of great advantage for SPH modelling and also for direct $N$-body simulations using the more efficient Ahmad-Cohen neighbour schemes.

We performed N-body/SPH simulations of merging gas-rich disk galaxies with mass ratios of 1:1 and 3:1. A stellar disk and bulge component and a dark halo was realized with collisionless particles, the gas was represented by SPH particles. Since we did not include star formation we focused on the gas dynamics and its influence on the structure of merger remnants. We find that equal-mass mergers are in general more effective in driving gas to the center. Around 50% of the gas resides in the central regions after the merger is complete. The gas in unequal-mass mergers keeps a large amount of its initial angular momentum and settles in a large scale disk while only 20% -30% of the gas is driven to the center. This general result is almost independent of the merger geometry. The gas in the outer regions accumulates in dense knots within tidal tails which could lead to the formation of open clusters or dwarf satellites. Later on, the gas knots loose angular momentum by dynamical friction and successively sink to the center of the remnant thereby increasing the total gas content of the disk. Due to the influence of the gas the simulated merger remnant becomes more oblate than its pure stellar counterpart. The shape of the stellar LOSVD changes and the third-order Gauss-Hermite coefficient $h_3$ is in good agreement with observations.

We investigate detailed kinematical properties of simulated collisionless merger remnants of disk galaxies with mass ratios of 1:1 and 3:1. The simulations were performed by direct summation using the new special hardware device GRAPE-5. In agreement with observations, the shape of the line-of-sight velocity distribution (LOSVD) is Gaussian with small deviations. For most cases we find that the retrograde wings of the LOSVD are steeper than the prograde ones. This is in contradiction with observations which show broad retrograde and steep prograde wings. This serious problem in the collisionless formation scenario of massive elliptical galaxies can be solved if all rotating ellipticals, even boxy ones, contain an additional stellar disk component with $\propto 15%$ of the total stellar mass and a scale length of order the effective radius of the spheroid. We propose that the progenitor galaxies of massive ellipticals must have contained a significant amount of gas that did not condense into stars during the merger process itself but formed an extended gaseous disk before the star formation epoch. The heating source that prevented the gas from forming stars early and the origin of the large specific angular momentum required for the gas component to form an extended disk are still unsolved problems.

We test the formation of low-luminosity elliptical galaxies through collisionless mergers of unequal-mass disk galaxies. The kinematic properties of a small survey of simulated merger end-products with initial disk mass ratios of 3:1 is compared to a sample of seven low-luminosity galaxies observed by Rix et al. that were chosen photometrically to be "ellipticals". In this paper, we go beyond a comparison in terms of global properties (using one number to characterize a model or a galaxy e.g. <a4>, ellipticity at some fixed radius or central velocity dispersion) and examine the detailed kinematics as function of galactocentric distance. The merger remnants are ``observed'' through a slit along the major and minor axis, using a pixel binning and slit width similar to the one used during the spectroscopic observations. We compare the rotational support (v/sigma) of the merger remnants, the normalized rotation on the minor axis and the major axis Gauss-Hermite moments h3, to that of the observed galaxies. The N-body remnants are very flattened when viewed edge-on (<epsilon> = 0.6) and should be inclined before making a fair comparison with the Rix et al. data set (which has <epsilon> = 0.3). When the merger remnants are appropriately inclined, their v/sigma profiles rise slower than the observed ones. Note that even when the remnants are viewed edge-on, the v/sigma profiles do not match the observations. The detailed comparison of the observations with purely collisionless 3:1-merger remnants suggests that this kind of dissipationless merger is not likely to be the dominant formation channel for low-luminosity elliptical galaxies. (abridged)

We performed N-body/SPH simulations of merging gas rich disk-galaxies with a mass ratio of 3:1. A stellar disk and bulge component and a dark halo is realized with collisionless particles, the gas is simulated using SPH particles. In this paper we focus on the gas dynamics in mergers and its influence on the structure of the final merger remnant, neglecting star formation. During the merger around $8\times 10^9 M_{\odot}$ of gas accumulate in the central region of the merger remnant (inside 300 pc). Part of the gas falls to the center, channeled by a tidally induced bar in the more massive galaxy. A peak accretion rate of 150 solar masses/year is reached when the two galaxy centers merge. It is likely that the gas will experience a central starburst and/or fuel a central black hole. The gas in the outer regions accumulates in dense knots within tidal tails which could lead to the formation of open clusters or dwarf satellites. Later on, the gas knots loose angular momentum by dynamical friction and also sink to the center of the remnant, thereby increasing the central gas content. In the end, the gas settles in a central power-law disk with surface density $\Sigma \propto r^{-2}$ surrounded by dense, orbiting gas clumps. The stellar profile is of de Vaucouleurs-type, like in simulations of pure stellar mergers

Collisionless N-body simulations of merging disk-galaxies with mass ratios ($\eta$) of 1:1, 2:1, 3:1, and 4:1 have been performed using direct summation with the special purpose hardware GRAPE. The simulations are used to examine whether the formation of elliptical galaxies can be explained in the context of the merger scenario. The photometric, kinematical and isophotal properties of the merger remnants are investigated and turn out to be in very good agreement with observations of giant elliptical galaxies. We conclude that equal mass mergers lead to slowly rotating, anisotropic remnants having predominantly boxy isophotes. Mergers with a mass ratio of 3:1 and 4:1, on the other hand, are fast isotropic rotators with disky isophotes. Projection effects can explain the observed scatter in the kinematical and isophotal properties of elliptical galaxies.

High-resolution numerical simulations of galaxy mergers are analysed. The global structure and isophotal shapes of the merger remnants are in good agreement with the observations. Whereas equal-mass mergers lead to anisotropic, boxy ellipticals, unequal-mass mergers result in disky and isotropic systems. The line-of-sight velocity distributions show small deviations from a Gaussian distribution. In all cases we find that the retrograde wings are steeper and that the prograde wings are broader than a Gaussian distribution. This is in contradiction with the observations which show broader retrograde and steeper prograde wings in all ellipticals. This fundamental difference between observation and theory can be explained if all ellipticals, even anistropic boxy ones, contain extended stellar disk components with luminosities of order 10 to 20 percent the total luminosity and scale radii of order the effective radii of the spheroids.

N-body simulations of major spiral-spiral mergers with different mass ratios q <= 4:1 are presented. All mergers lead to spheroidal stellar systems with de Vaucouleurs-like surface density profiles, resembling elliptical galaxies. The internal dynamical structure of the remnants does however depend critically on the adopted mass ratio. Equal-mass mergers with q <= 2:1 form radially anisotropic, slowly rotating systems with preferentially boxy isophotes. Unequal-mass mergers with q > 2:1 result in the formation of rotationally supported ellipticals with disky isophotes and small minor axis rotation. Projection effects lead in general to a large scatter in the kinematical and isophotal properties. If the observed dichotomy between boxy and disky ellipticals is a result of different mass ratios of the merger components massive ellipticals must have formed preferentially from equal-mass mergers in contrast to low-mass ellipticals which had to form from unequal-mass mergers.

The origin of boxy and disky elliptical galaxies is investigated. The results of two collisionless N-body simulations of spiral-spiral mergers with mass ratios of 1:1 and 3:1 are discussed and the projected properties of the merger remnants are investigated. It is shown that the equal-mass merger leads to an anisotropic, slowly rotating system with preferentially boxy isophotes and significant minor axis rotation. The unequal-mass merger results in the formation of a rotationally supported elliptical with disky isophotes and small minor axis rotation. The observed scatter in the kinematical and isophotal properties of both classes of elliptical galaxies can be explained by projection effects.

We study a possible formation mechanism for elliptical/S0 galaxies using N-body simulations with GRAPE. A close galaxy encounter which does not lead to a merger can induce a strong bar in an initially axisymmetric disk. This bar is unstable against bending oscillations and loses angular momentum to the dark halo. After 6 Gyrs the remnant is spheroidal with a de Vaucouleurs-type surface profile and has ellipticities and rotation properties comparable to observed elliptical/S0 galaxies: the system rotates fast (v/sigma \approx 1) and has disky isophotes.