Search

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 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., Comment: 17 pages, 10 figures, published by MNRAS

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., Comment: submitted to ApJ

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., Comment: Accepted for publication in MNRAS. 53 pages, 16 figures

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., Comment: 18 pages, 8 figures, MNRAS accepted

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., Comment: ApJ Letters in press

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., Comment: Accepted for publication in MNRAS. Changes (including a pamameter study in Appendix C) highlighted

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., Comment: 7 pages, 5 figures, ApJL, accepted

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., Comment: 26 pages, MNRAS, accepted

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., Comment: 57 pages, 16 figures, 3 tables. ApJ accepted

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., Comment: 13 pages, 9 figures, 3 tables

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., Comment: 33 pages, 18 figures, 5 tables. Accepted for publication in MNRAS. Full length tables available at http://www.purl.org/atlas3d

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., Comment: 26 pages, 11 figures, accepted for publication in MNRAS

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., Comment: MNRAS in press, On-line version of the catalog of galaxies with available deep images and figures (including navigation tools) are available from http://irfu.cea.fr/Projets/matlas/public/Atlas3D/atlas3d-pXXIX.html

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., Comment: 5 pages, 3 figures, accepted for publication in ApJL

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., Comment: 19 pages, 19 figures. Accepted for publication in MNRAS

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., Comment: 21 pages, 19 figures

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., Comment: 22 pages, submitted to MNRAS, comments welcome

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., Comment: 7 pages, 4 figures, ApJ, 789, L40

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., Comment: 19 pages, 13 figures, accepted for publication in MNRAS

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., Comment: 23 pages, 29 figures, accepted for publication in MNRAS

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., Comment: 16 (+4) pages, 10 figures, accepted for publication in MNRAS

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., Comment: 13 pages, 11 figures, accepted for the publication in MNRAS

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., Comment: Replaced with the version accepted by MNRAS

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., Comment: MNRAS accepted

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., Comment: 11 pages, 7 figures, accepted by the Astrophysical Journal

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., Comment: 5 pages, 5 figures, accepted in ApJL; corrected title

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^111 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., Comment: 23 pages, 16 figures, and 1 table. Accepted for publication in MNRAS

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., Comment: 16 pages, 10 figures, resubmitted to MNRAS

[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., Comment: 32 pages, 15 figures, accepted for publications on MNRAS. More information about the Atlas3D project and machine readable Table C1 see http://www-astro.physics.ox.ac.uk/atlas3d/. Corrected typos

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., Comment: 14 pages, 14 figures, 2 tables. Accepted for publication in MNRAS

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., Comment: 26 pages, 20 figures, accepted for publication in MNRAS

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

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], Comment: 16 pages, 10 figures, accepted for publications on MNRAS. More information about our Atlas3D project see http://www-astro.physics.ox.ac.uk/atlas3d/

[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., Comment: 28 pages, Accepted in MNRAS, 7 figures, 6 tables

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., Comment: 15 pages, 8 figures, accepted for publication in MNRAS

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., Comment: To appear in the Proceedings of IAU Symposium No. 295: "The intriguing life of massive galaxies", D. Thomas, A. Pasquali and I. Ferreras, eds

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., Comment: 13 pages, 9 figures, accepted by ApJ

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], Comment: 23 pages, 20 figures. Accepted for publication in MNRAS

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., Comment: Accepted by the Astrophysical Journal, 17 pages, 12 figures

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), Comment: 22 pages + appendix, 12 figures, accepted to MNRAS

[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., Comment: Accepted for publication in MNRAS. 29 pages 14 figures. More information about our Atlas3D project is available at http://purl.org/atlas3d

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., Comment: 52 pages, 5 figures (+ 33 appendix figures), MNRAS, accepted

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., Comment: MNRAS, accepted

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), Comment: 17 pages, 15 figures, 2 tables, accepted by MNRAS

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], Comment: 35 pages, 18 figures, LaTeX. Accepted for publication in MNRAS. Includes a new data table and figures of all dynamical models. More information about our Atlas3D project is available at http://purl.org/atlas3d

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], Comment: 35 pages, 16 figure, LaTeX. Accepted for publication in MNRAS. Includes a new data table. More information about our Atlas3D project is available at http://purl.org/atlas3d

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)., Comment: 7 pages, 3 figures, accepted by ApJ