Researcher profile

T. Naab

T. Naab contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
46works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

46 published item(s)

preprint2022arXiv

CO and [CII] line emission of molecular clouds -- the impact of stellar feedback and non-equilibrium chemistry

We analyse synthetic $^{12}$CO, $^{13}$CO, and [CII] emission maps of simulated molecular clouds of the SILCC-Zoom project, which include an on-the-fly evolution of H$_2$, CO, and C$^+$. We use simulations of hydrodynamical and magnetohydrodynamical clouds, both with and without stellar feedback. We introduce a novel post-processing of the C$^+$ abundance using CLOUDY, to account for further ionization states of carbon due to stellar radiation. We report the first self-consistent synthetic emission maps of [CII] in feedback bubbles, largely devoid of emission inside them, as recently found in observations. The C$^+$ mass is only poorly affected by stellar feedback but the [CII] luminosity increases by $50 - 85$ per cent compared to runs without feedback. Furthermore, we investigate the capability of the CO/[CII] line ratio as a tracer of the amount of H$_2$ in the clouds and their evolutionary stage. We obtain, for both $^{12}$CO and $^{13}$CO, no clear trend of the luminosity ratio, $L_\mathrm{CO}/L_\mathrm{[CII]}$. It can therefore \textit{not} be used as a reliable measure of the H$_2$ mass fraction. We note a monotonic relation between $L_\mathrm{CO}/L_\mathrm{[CII]}$ and the H$_2$ fraction when considering the ratio for individual pixels of our synthetic maps, but with large scatter. Moreover, we show that assuming chemical equilibrium results in an overestimation of H$_2$ and CO masses of up to 110 and 30 per cent, respectively, and in an underestimation of H and C$^+$ masses of 65 and 7 per cent, respectively. In consequence, $L_\mathrm{CO}$ would be overestimated by up to 50 per cent, and $L_\mathrm{C[II]}$ be underestimated by up to 35 per cent. Hence, the assumption of chemical equilibrium in molecular cloud simulations introduces intrinsic errors of a factor of up to $\sim2$ in chemical abundances, luminosities and luminosity ratios.

preprint2020arXiv

From Nuclear to Circumgalactic: Zooming in on AGN-Driven Outflows at z~2.2 with SINFONI

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

preprint2016arXiv

The angular momentum distribution and baryon content of star forming galaxies at z~1-3

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$α$ 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 <$λ$> ~ 0.037 and its dispersion ($σ_{log(λ)}$~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 $λ$ x (j_d/j_DM) and the nuclear stellar density $Σ_{*}$(1kpc) favors a scenario where disk-internal angular momentum redistribution leads to &#34;compaction&#34; 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.

preprint2016arXiv

The low dark matter content of the lenticular galaxy NGC 3998

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&#39;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.

preprint2015arXiv

Comparison of simple mass estimators for slowly rotating elliptical galaxies

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&#39;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 $σ_p(R_2)$ as $M_{vir} \left[M_{\odot}h^{-1}\right] \approx 6\cdot 10^{12} \left( \frac{σ_p(R_2)}{200\, \rm km\, s^{-1}} \right)^{4}$ with RMS scatter $\approx 40 \%$.

preprint2015arXiv

Dynamical evolution of massive black holes in galactic-scale N-body simulations - introducing the regularized tree code &#34;rVINE&#34;

We present a hybrid code combining the OpenMP-parallel tree code VINE with an algorithmic chain regularization scheme. The new code, called &#34;rVINE&#34;, 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&#39;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.

preprint2015arXiv

Modelling the supernova-driven ISM in different environments

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$_α$ can be as high as $70$ km s$^{-1}$ due to the contribution from young, embedded SN remnants.

preprint2014arXiv

Combined CO & Dust Scaling Relations of Depletion Time and Molecular Gas Fractions with Cosmic Time, Specific Star Formation Rate and Stellar Mass

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.

preprint2014arXiv

Evidence for Wide-Spread AGN Driven Outflows in the Most Massive z~1-2 Star Forming Galaxies

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.

preprint2014arXiv

The energy and momentum input of supernova explosions in structured and ionised molecular clouds

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.

preprint2014arXiv

The SILCC (SImulating the LifeCycle of molecular Clouds) project: I. Chemical evolution of the supernova-driven ISM

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 $Σ_{_{\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.

preprint2014arXiv

The SINS/zC-SINF survey of z~2 galaxy kinematics: Evidence for powerful AGN-driven nuclear outflows in massive star-forming galaxies

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.

preprint2013arXiv

Constrained simulations of the Antennae Galaxies: Comparison with Herschel-PACS observations

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.

preprint2013arXiv

Nebular Excitation in z~2 Star-forming Galaxies from the SINS and LUCI Surveys: The Influence of Shocks and AGN

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&#39; 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.

preprint2013arXiv

PHIBSS: molecular gas, extinction, star formation and kinematics in the z=1.5 star forming galaxy EGS13011166

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 &#34;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 &#39;mixed&#39; 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.

preprint2013arXiv

The ATLAS$^{\rm{3D}}$ project - XXV: Two-dimensional kinematic analysis of simulated galaxies and the cosmological origin of fast and slow rotators

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 $λ_{\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/σ$. 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/σ$. 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)

preprint2013arXiv

The Atlas3D project -- XXVII. Cold Gas and the Colours and Ages of Early-type Galaxies

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

preprint2013arXiv

The SINS/zC-SINF survey of z~2 galaxy kinematics: evidence for gravitational quenching

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.

preprint2013arXiv

Towards a resolved Kennicutt-Schmidt law at high redshift

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&#34;. 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.

preprint2012arXiv

AGN Feedback Driven Molecular Outflow in NGC 1266

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.

preprint2012arXiv

PHIBSS: molecular gas content and scaling relations in z~1-3 normal star forming galaxies

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.

preprint2012arXiv

Quenching of Star Formation in Molecular Outflow Host NGC 1266

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β 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).

preprint2012arXiv

Testing a simple recipe for estimating galaxy masses from minimal observational data

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 ($σ\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 ($σ< 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.

preprint2011arXiv

Galactic ménage à trois: Simulating magnetic fields in colliding galaxies

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.

preprint2011arXiv

Galaxy formation in semi-analytic models and cosmological hydrodynamic zoom simulations

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 &#34;hot&#34; relative to &#34;cold&#34; 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.

preprint2011arXiv

Investigating the merger origin of Early-Type Galaxies using ultra-deep optical images

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.

preprint2011arXiv

Molecular Gas and Star Formation in Local Early-Type Galaxies

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.

preprint2011arXiv

The metallicity dependence of the CO {\rightarrow} H_2 conversion factor in z>1 star forming galaxies

We use the first systematic samples of CO millimeter emission in z>1 &#39;main-sequence&#39; star forming galaxies (SFGs) to study the metallicity dependence of the conversion factor α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 &#39;CO-dark&#39; 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(α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.

preprint2011arXiv

The SINS survey of z~2 galaxy kinematics: properties of the giant star forming clumps

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α/[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.

preprint2010arXiv

A Study of the Gas-Star Formation Relation over Cosmic Time

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 &#39;main-sequence&#39;, 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.

preprint2010arXiv

Formation of slowly rotating early-type galaxies via major mergers: a Resolution Study

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.

preprint2010arXiv

High molecular gas fractions in normal massive star forming galaxies in the young Universe

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&#39;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.

preprint2010arXiv

High-Redshift Star-Forming Galaxies: Angular Momentum and Baryon Fraction, Turbulent Pressure Effects and the Origin of Turbulence

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.

preprint2010arXiv

NGC6240: Merger-Induced Star Formation & Gas Dynamics

We present spatially resolved integral field spectroscopic K-band data at a resolution of 0.13&#34; (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&#39;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&#39; bulges.

preprint2010arXiv

On the evolution of the intrinsic scatter in black hole versus galaxy mass relations

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 σ is well approximated by σmerg(m) = σ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.

preprint2010arXiv

Simulating magnetic fields in the Antennae galaxies

We present self-consistent high-resolution simulations of NGC4038/4039 (the &#34;Antennae galaxies&#34;) 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.

preprint2010arXiv

The Effects of X-Ray and UV Background Radiation on the Low-Mass Slope of the Galaxy Mass Function

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.

preprint2009arXiv

Magnetic field structure due to the global velocity field in spiral 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.

preprint2009arXiv

Protostellar discs formed from turbulent cores

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.

preprint2009arXiv

Protostellar disks formed from rigidly rotating cores

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 Ω_O. Protostellar discs forming from cores with Ω_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 Ω_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 Ω_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.

preprint2009arXiv

The halo merger rate in the Millennium Simulation and implications for observed galaxy merger fractions

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 &#34;knee&#34; 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/.

preprint2009arXiv

The SINS survey: modeling the dynamics of z~2 galaxies and the high-z Tully-Fisher relation

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.

preprint2008arXiv

Mergers and Mass Accretion Rates in Galaxy Assembly: The Millennium Simulation Compared to Observations of z~2 Galaxies

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.

preprint2008arXiv

The flattening and the orbital structure of early-type galaxies and collisionless N-body binary disk mergers

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.

preprint2006arXiv

Triggered Star Formation in the Environment of Young Massive Stars

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.