BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Sadegh Khochfar

Sadegh Khochfar appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-ph.HEastro-phastro-ph.SRgr-qc

Catalog footprint

What is connected

82works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

82 published item(s)

preprint2022arXiv

FOREVER22: galaxy formation in protocluster regions

We present results from a new cosmological hydrodynamics simulation campaign of protocluster (PC) regions, FOREVER22: FORmation and EVolution of galaxies in Extremely-overdense Regions motivated by SSA22. The simulations cover a wide range of cosmological scales using three different zoom set-ups in a parent volume of $(714.2~\rm cMpc)^{3}$: PCR (Proto-Cluster Region; $V= (28.6~{\rm cMpc})^{3} $, SPH particle mass, $m_{\rm{SPH}} = 4.1 \times 10^{6}~\rm M_{\odot}$ and final redshift, $z_{\rm end}=2.0$), BCG (Brightest proto-Cluster Galaxy; $V \sim (10~{\rm cMpc})^{3} $, $m_{\rm SPH} = 5.0\times10^{5}~\rm M_{\odot}$ and $z_{\rm end}=4.0$ ), and First ( $V \sim (3~{\rm cMpc})^{3} $, $m_{\rm SPH} = 7.9 \times 10^{3}~\rm M_{\odot}$ and $z_{\rm end}=9.5$) runs, that allow to focus on different aspects of galaxy formation. In the PCR runs, we follow 10 PCs, each harbouring 1 - 4 SMBHs with $M_{\rm BH} \ge 10^{9}~\rm M_{\odot}$. One of the PC cores shows a spatially close arrangement of seven starburst galaxies with ${\rm SFR} \gtrsim 100~\rm M_{\odot}~yr^{-1}$ each, that are dust-obscured and would appear as submillimeter galaxies with flux $\gtrsim 1~$ mJy at $1.1~ \rm mm$ in observations. The BCG runs show that the total SFRs of haloes hosting BCGs are affected by AGN feedback, but exceed $1000~\rm M_{\odot}~yr^{-1}$ at $z \lesssim 6$. The First runs resolve mini-haloes hosting population (Pop) III stars and we show that, in PC regions, the dominant stellar population changes from Pop III to Pop II at $z \gtrsim 20$, and the first galaxies with ${\rm SFR} \gtrsim 18~\rm M_{\odot}~yr^{-1}$ form at $z \sim 10$. These can be prime targets for future observations with the James Webb Space Telescope. Our simulations successfully reproduce the global star formation activities in observed PCs and suggest that PCs can kickstart cosmic reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multi-Epoch Machine Learning 1: Unravelling Nature vs Nurture for Galaxy Formation

We present a novel machine learning method for predicting the baryonic properties of dark matter only subhalos from N-body simulations. Our model is built using the extremely randomized tree (ERT) algorithm and takes subhalo properties over a wide range of redshifts as its input features. We train our model using the IllustrisTNG simulations to predict blackhole mass, gas mass, magnitudes, star formation rate, stellar mass, and metallicity. We compare the results of our method with a baseline model from previous works, and against a model that only considers the mass history of the subhalo. We find that our new model significantly outperforms both of the other models. We then investigate the predictive power of each input by looking at feature importance scores from the ERT algorithm. We produce feature importance plots for each baryonic property, and find that they differ significantly. We identify low redshifts as being most important for predicting star formation rate and gas mass, with high redshifts being most important for predicting stellar mass and metallicity, and consider what this implies for nature vs nurture. We find that the physical properties of galaxies investigated in this study are all driven by nurture and not nature. The only property showing a somewhat stronger impact of nature is the present-day star formation rate of galaxies. Finally we verify that the feature importance plots are discovering physical patterns, and that the trends shown are not an artefact of the ERT algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Cosmic Mach Number as an Environment Measure for the Underlying Dark Matter Density Field

Using cosmological dark matter only simulations of a $(1.6$ Gpc$/h)^3$ volume from the Legacy simulation project, we calculate Cosmic Mach Numbers (CMN) and perform a theoretical investigation of their relation with halo properties and features of the density field to gauge their use as an measure of the environment. CMNs calculated on individual spheres show correlations with both the overdensity in a region and the density gradient in the direction of the bulk flow around that region. To reduce the scatter around the median of these correlations, we introduce a new measure, the rank ordered Cosmic Mach number ($\hat{\mathcal{M}}_g$), which shows a tight correlations with the overdensity $δ=\frac{ρ-\barρ}{\barρ}$. Measures of the large scale density gradient as well as other average properties of the halo population in a region show tight correlations with $\hat{\mathcal{M}}_g$ as well. Our results in this first empirical study suggest that $\hat{\mathcal{M}}_g$ is an excellent proxy for the underlying density field and hence environment that can circumvent reliance on number density counts in a region. For scales between $10$ and $100 Mpc$/h, Mach numbers calculated using dark matter halos $(> 10^{12}$ M$_{\odot})$ that would typically host massive galaxies are consistent with theoretical predictions of the linear matter power spectrum at a level of $10\%$ due to non-linear effects of gravity. At redshifts $z\geq 3$, these deviations disappear. We also quantify errors due to missing large scale modes in simulations. Simulations of box size $\leq 1 $ Gpc/$h$ typically predict CMNs 10-30\% too small on scales of$\sim 100$ Mpc$/h$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Turbulent Cold Flows Gave Birth to the First Quasars

How quasars powered by supermassive black holes (SMBHs) formed less than a billion years after the Big Bang is still one of the outstanding problems in astrophysics 20 years after their discovery$^{1-4}$. Cosmological simulations suggest that rare cold flows converging on primordial haloes in low-shear environments could have created these quasars if they were 10$^4$ - 10$^5$ M$_{\odot}$ at birth but could not resolve their formation$^{5-8}$. Semianalytical studies of the progenitor halo of a primordial quasar found that it favours the formation of such seeds but could not verify if one actually appeared$^9$. Here we show that a halo at the rare convergence of strong, cold accretion flows creates massive BH seeds without the need for UV backgrounds, supersonic streaming motions, or even atomic cooling. Cold flows drive violent, supersonic turbulence in the halo that prevents star formation until it reaches a mass that triggers sudden, catastrophic baryon collapse that forms 31,000 and 40,000 M$_{\odot}$ stars. This simple, robust process ensures that haloes capable of forming quasars by z $>$ 6 produce massive seeds. The first quasars were thus a natural consequence of structure formation in cold dark matter cosmologies, not exotic, finely-tuned environments as previously thought$^{10-14}$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Mapping the inhomogeneous Universe with Standard Sirens: Degeneracy between inhomogeneity and modified gravity theories

The detection of gravitational waves (GWs) and an accompanying electromagnetic (E/M) counterpart have been suggested as a future probe for cosmology and theories of gravity. In this paper, we present calculations of the luminosity distance of sources taking into account inhomogeneities in the matter distribution that are predicted in numerical simulations of structure formation. In addition, we show that inhomogeneities resulting from clustering of matter can mimic certain classes of modified gravity theories, or other effects that dampen GW amplitudes, and deviations larger than $δν\sim \mathcal{O}(0.1)\ (99\%\ \rm{C.L.})$ to the extra friction term $ν$, from zero, would be necessary to distinguish them. For these, we assume mock GWs sources, with known redshift, based on binary population synthesis models, between redshifts $z=0$ and $z=5$. We show that future GW detectors, like Einstein Telescope or Cosmic Explorer, will be needed for strong constraints on the inhomogeneity parameters and breaking the degeneracy between modified gravity effects and matter anisotropies by measuring $ν$ at $5 \%$ and $1 \%$ level with $100$ and $350$ events respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Birth Mass Function of Pop III Stars

Population III (Pop III) stars ended the cosmic Dark Ages and began early cosmological reionization and chemical enrichment. However, in spite of their importance to the evolution of the early Universe, their properties remain uncertain because of limitations to previous numerical simulations and the lack of any observational constraints. Here we investigate Pop III star formation in five primordial halos with 3D radiation-hydrodynamical cosmological simulations. We find that multiple stars form in each minihalo and that their numbers increase over time, with up to 23 stars forming in one of the halos. Radiative feedback from the stars generates strong outflows, deforms the surrounding protostellar disk, and delays star formation for a few thousand years. Star formation rates vary with halo and depend on mass accretion onto the disk, halo spin number, and the fraction of massive stars in the halo. Stellar masses in our models range from 0.1-37 $\rm M_{\odot}$, and of the 55 stars that form in our models twelve are $\rm > 10~ M_{\odot}$ and most of the others are 1-10 $\rm M_{\odot}$. Our simulations thus suggest that Pop III stars have characteristic masses of 1-10 $\rm M_{\odot}$ and a top-heavy IMF with dN/dM $\propto M_*^{-1.18}$. Up to 70\% of the stars are ejected from their disks by three-body interactions which, along with ionizing UV feedback, limits their final masses.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Calibration of a star formation and feedback model for cosmological simulations with Enzo

We present results from seventy-one zoom simulations of a Milky Way-sized (MW) halo, exploring the parameter space for a widely-used star formation and feedback model in the {\tt Enzo} simulation code. We propose a novel way to match observations, using functional fits to the observed baryon makeup over a wide range of halo masses. The model MW galaxy is calibrated using three parameters: the star formation efficiency $\left(f_*\right)$, the efficiency of thermal energy from stellar feedback $\left(ε\right)$ and the region into which feedback is injected $\left(r\ {\rm and}\ s\right)$. We find that changing the amount of feedback energy affects the baryon content most significantly. We then identify two sets of feedback parameter values that are both able to reproduce the baryonic properties for haloes between $10^{10}\,\mathrm{M_\odot}$ and $10^{12}\,\mathrm{M_\odot}$. We can potentially improve the agreement by incorporating more parameters or physics. If we choose to focus on one property at a time, we can obtain a more realistic halo baryon makeup. We show that the employed feedback prescription is insensitive to dark matter mass resolution between $10^5\,{\rm M_\odot}$ and $10^7\,{\rm M_\odot}$. Contrasting both star formation criteria and the corresponding combination of optimal feedback parameters, we also highlight that feedback is self-consistent: to match the same baryonic properties, with a relatively higher gas to stars conversion efficiency, the feedback strength required is lower, and vice versa. Lastly, we demonstrate that chaotic variance in the code can cause deviations of approximately 10\% and 25\% in the stellar and baryon mass in simulations evolved from identical initial conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inception of a first quasar at cosmic dawn

Earliest quasars at the cosmic dawn are powered by mass accretion onto supermassive black holes of a billion solar masses. Massive black hole seeds forming through the direct collapse mechanism are considered the most promising candidates but how do they grow and co-evolve with their host galaxies at early cosmic times remains unknown. We here present results from a cosmological radiation hydrodynamical simulation including self-consistent modeling of both Pop III and Pop II star formation, their radiative and supernova feedback in the host galaxy along with X-ray feedback from an accreting massive black hole (MBH) of $\rm 10^5 ~M_{\odot}$ in a halo of $\rm 2 \times 10^9~M_{\odot}$ from $z=26$ down to $z=16$. Our results show that energy deposition from X-rays in the proximity of MBH suppresses Pop III star formation for about 12 Myr while at the same time these X-rays catalyze $\rm H_2$ formation which leads to the formation of a Pop III star cluster of 500 $\rm M_{\odot}$ in the close vicinity of the MBH. We find that mode of star formation for Pop III is episodic and bursty due to the clumpy accretion while for Pop II it is continuous. The stellar mass of the host galaxy at $z \sim 16$ is $\rm 2 \times 10^7~M_{\odot}$ with a star formation rate (SFR) of $\rm \sim 0.1-1~M_{\odot}/yr$. In total, the MBH accretes $\rm 1.5 \times 10^6~M_{\odot}$ during 120 Myr with the mean accretion rate of $\rm \sim 0.01~M_{\odot}/yr$ corresponding to an average Eddington fraction of 50\%.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Birth of Binary Direct-Collapse Black Holes

Supermassive primordial stars forming during catastrophic baryon collapse in atomically-cooling halos at $z \sim$ 15 - 20 may be the origin of the first quasars in the universe. However, no simulation to date has followed the evolution of these halos at resolutions that are high enough or for times that are long enough to determine if collapse actually produces SMSs. Here we report new cosmological simulations of baryon collapse in atomically-cooled halos for times that are long enough for SMSs to form and die as direct-collapse black holes (DCBHs). We find that the high infall rates required to build up such stars do persist until the end of their lives and could fuel the rapid growth of their BHs thereafter. Our simulations also demonstrate that binary and even small multiples of SMSs can form in low-spin and high-spin halos, respectively. This discovery raises the exciting prospect of detecting gravitational waves from DCBH mergers with LISA and tidal disruption events in the near infrared with the {\em James Webb Space Telescope} and ground-based telescopes in the coming decade.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The First Billion Years Project: Finding Infant Globular Clusters at z=6

We explored a suite of high-resolution cosmological simulations from the First Billion Years Project (FiBY) at $z \geq 6$. All substructures within the simulations have been identified with the SUBFIND algorithm. From our analysis, two distinct groups of objects emerge. We hypothesise that the substructures in the first group, which appear to have a high baryon fraction ($f_{\rm b} \geq 0.95$), are possible infant GC candidates. Objects belonging to the second group have a high stellar fraction ($f_{\rm star} \geq 0.95$) and show a potential resemblance to infant ultra-faint dwarf galaxies. The high baryon fraction objects identified in this study are characterised by a stellar content similar to the one observed in present-day GCs, but they still contain a high gas fraction ($f_{\rm gas} \sim 0.95$) and a relatively low amount of dark matter. They are compact, dense systems. Their sizes are consistent with recent estimates based on the first observations of possible proto-GCs at high redshifts. These types of infant GC candidates appear to be more massive and more abundant in massive host galaxies, indicating that the assembly of galaxies via mergers may play an important role in building several GC-host scaling relations. Specifically, we express the relation between the mass of the most massive infant GC and its host stellar mass as $\log(M_{\rm cl}) = (0.31\pm0.15)\log(M_{\rm *,gal} + (4.17\pm1.06)$. We also report a new relation between the most massive infant GC and the parent specific star formation rate of the form $\log(M_{\rm cl}) = (0.85\pm0.30)\log(sSFR) + α$ that describes the data at both low and high redshift. Finally, we assess the present-day GC mass (GC number) -- halo mass relation offers a satisfactory description of the behaviour of our infant GC candidates at high redshift, suggesting that such a relation may be set at formation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Galaxy evolution and radiative properties in the early Universe: multi-wavelength analysis in cosmological simulations

Recent observations have successfully detected UV or infrared flux from galaxies at the epoch of reionization. However, the origin of their radiative properties has not been fully understood yet. Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present theoretical predictions of multi-wavelength radiative properties of the first galaxies at z=6-15. We find that most of the gas and dust are ejected from star-forming regions due to supernova (SN) feedback, which allows UV photons to escape. We show that the peak of SED rapidly shifts between UV and infrared wavelengths on a timescale of 100 Myr due to intermittent star formation and feedback. When dusty gas covers the star-forming regions, the galaxies become bright in the observed-frame sub-millimeter wavelengths. In addition, we find that the escape fraction of ionizing photons also changes between 1-40% at z>10. The mass fraction of HII region changes with the star formation history, resulting in the fluctuations of metal lines and Lyman-alpha line luminosities. In the starbursting phase of galaxies with the halo mass $\sim 10^{11}\,{\rm M_{\odot}}$ ($10^{12}\,{\rm M_{\odot}}$), the simulated galaxy has $L_{\rm [O\,III]} \sim 10^{42}\, (10^{43})\, {\rm erg\,s^{-1}}$, which is consistent with the observed star-forming galaxies at z>7. Our simulations suggest that deep [C II] observation with ALMA can trace the distribution of neutral gas extending over $\sim 20$ physical kpc. We also find that the luminosity ratio $L_{\rm [O\,III]}/L_{\rm [C\,II]}$ decreases with bolometric luminosity due to metal enrichment. Our simulations show that the combination of multi-wavelength observations by ALMA and JWST will be able to reveal the multi-phase ISM structure and the transition from starbursting to outflowing phases of high-z galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Detecting Direct Collapse Black Holes: making the case for CR7

We propose that one of the sources in the recently detected system CR7 by Sobral et al. (2015) through spectro-photometric measurements at $z = 6.6$ harbors a direct collapse blackhole (DCBH). We argue that the LW radiation field required for direct collapse in source A is provided by sources B and C. By tracing the LW production history and star formation rate over cosmic time for the halo hosting CR7 in a $Λ$CDM universe, we demonstrate that a DCBH could have formed at $z\sim 20$. The spectrum of source A is well fit by nebular emission from primordial gas around a BH with MBH $\sim 4.4 \times 10^6 \ \rm M_{\odot}$ accreting at a 40 % of the Eddington rate, which strongly supports our interpretation of the data. Combining these lines of evidence, we argue that CR7 might well be the first DCBH candidate.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

New constraints on direct collapse black hole formation in the early Universe

Direct collapse black holes (DCBH) have been proposed as a solution to the challenge of assembling supermassive black holes by $z>6$ to explain the bright quasars observed at this epoch. The formation of a DCBH seed with $\rm M_{BH}\sim10^{4-5}\ \rm M_{\odot}$ requires a pristine atomic-cooling halo to be illuminated by an external radiation field that is sufficiently strong to entirely suppress H$_{2}$ cooling in the halo. Many previous studies have attempted to constrain the critical specific intensity that is likely required to suppress H$_{2}$ cooling, denoted as $J_{\rm crit}$. However, these studies have typically assumed that the incident external radiation field can be modeled with a black-body spectrum. Under this assumption, it is possible to derive a {unique} value for $J_{\rm crit}$ that depends only on the temperature of the black-body. In this study we consider a more realistic spectral energy distribution (SED) for the external source of radiation that depends entirely on its star formation history and age. The rate of destruction of the species responsible for suppressing molecular hydrogen cooling depends on the detailed shape of the SED. Therefore the value of $J_{\rm crit}$ is tied to the shape of the incident SED of the neighbouring galaxy. We fit a parametric form to the rates of destruction of H$_2$ and H$^-$ that permit direct collapse. Owing to this, we find that $J_{\rm crit}$ is not a fixed threshold but can lie anywhere in the range $J_{\rm crit} \sim 0.5$--$10^{3}$, depending on the details of the source stellar population.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nuclear discs as clocks for the assembly history of early-type galaxies: the case of NGC4458

Approximately 20% of early-type galaxies host small nuclear stellar discs that are tens to a few hundred parsecs in size. Such discs are expected to be easily disrupted during major galactic encounters, hence their age serve to constrain their assembly history. We use VIMOS integral-field spectroscopic observations for the intermediate-mass E0 galaxy NGC4458 and age-date its nuclear disc via high-resolution fitting of various model spectra. We find that the nuclear disc is at least 6 Gyr old. A clue to gain narrow limits to the stellar age is our knowledge of the nuclear disc contribution to the central surface brightness. The presence of an old nuclear disk, or the absence of disruptive encounters since z~0.6, for a small galaxy such as NGC4458 which belongs to the Virgo cluster, may be consistent with a hierarchical picture for galaxy formation where the smallest galaxies assembles earlier and the crowded galactic environments reduce the incidence of galaxy mergers. On the other hand, NGC4458 displays little or no bulk rotation except for a central kpc-scale kinematically-decoupled core. Slow rotation and decoupled core are usually explained in terms of mergers. The presence and age of the nuclear disc constraint these mergers to have happened at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Atlas3D project -- XXXI. Nuclear radio emission in nearby early-type galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The role of stellar relaxation in the formation and evolution of the first massive black holes

We present calculations on the formation of massive black holes with 10^5 Msun at z > 6 that can be the seeds of supermassive black holes at z > 6. Under the assumption of compact star cluster formation in merging galaxies, star clusters in haloes of 10^8 ~ 10^9 Msun can undergo rapid core-collapse leading to the formation of very massive stars (VMSs) with ~1000 Msun which directly collapse into black holes with similar masses. Star clusters in halos of > 10^9 Msun experience type-II supernovae before the formation of VMSs due to long core-collapse time scales. We also model the subsequent growth of black holes via accretion of residual stars in clusters. 2-body relaxation efficiently re-fills the loss cones of stellar orbits at larger radii and resonant relaxation at small radii is the main driver for accretion of stars onto black holes. As a result, more than ninety percent of stars in the initial cluster are swallowed by the central black holes before z=6. Using dark matter merger trees we derive black hole mass functions at z=6-20. The mass function ranges from 10^3 to 10^5 Msun at z <~ 15. Major merging of galaxies of >~ 4*10^8 Msun at z ~ 20 successfully leads to the formation of >~ 10^5 Msun BHs by z >~ 10 which can be the potential seeds of supermassive black holes seen today.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Can the 21 cm signal probe Population III and II star formation?

Using varying models for the star formation rate (SFR) of Population (Pop) III and II stars at z>6 we derive the expected redshift history of the global 21 cm signal from the inter-galactic medium (IGM). To recover the observed Thomson scattering optical depth of the cosmic microwave background (CMB) requires SFRs at the level of ~ 10^-3 Msun/yr/Mpc^3 at z ~ 15 from Pop III stars, or ~ 10^-1Msun/yr/Mpc^3 at z ~ 7 from Pop II stars. In the case the SFR is dominated by Pop III stars, the IGM quickly heats above the CMB at z > 12 due to heating from supernovae. In addition, Lyman-alpha photons from haloes hosting Pop III stars couple the spin temperature to that of the gas, resulting in a deep absorption signal. If the SFR is dominated by Pop II stars, the IGM slowly heats and exceeds the CMB temperature at z ~ 10. However, the larger and varying fraction of Pop III stars is able to break this degeneracy. We find that the impact of the initial mass function (IMF) of Pop III stars on the 21 cm signal results in an earlier change to a positive signal if the IMF slope is ~ -1.2. Measuring the 21 cm signal at z > 10 with next generation radio telescopes such as the Square Kilometre Array will be able to investigate the contribution from Pop III and Pop II stars to the global star formation rate.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Atlas3D Project - XXX. Star formation histories and stellar population scaling relations of early-type galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The First Billion Years Project: The escape fraction of ionizing photons in the epoch of reionization

Proto-galaxies forming in low-mass dark matter haloes are thought to provide the majority of ionizing photons needed to reionize the Universe, due to their high escape fractions of ionizing photons. We study how the escape fraction in high-redshift galaxies relates to the physical properties of the halo in which the galaxies form, by computing escape fractions in more than 75000 haloes between redshifts 27 and 6 that were extracted from the First Billion Years project, high-resolution cosmological hydrodynamics simulations of galaxy formation. We find that the main constraint on the escape fraction is the gas column density in a radius of 10 pc around the stellar populations, causing a strong mass dependence of the escape fraction. The lower potential well in haloes with virial mass below 1e8 solar mass results in low column densities that can be penetrated by radiation from young stars (age < 5 Myr). In haloes with virial mass above 1e8 solar mass supernova feedback is important, but only 30% of the haloes in this mass range have an escape fraction higher than 1%. We find a large range of escape fractions in haloes with similar properties, caused by different distributions of the dense gas in the halo. This makes it very hard to predict the escape fraction on the basis of halo properties and results in a highly anisotropic escape fraction. The strong mass dependence, the large spread and the large anisotropy of the escape fraction may strongly affect the topology of reionization and is something current models of cosmic reionization should strive to take into account.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The First Population II Stars Formed in Externally Enriched Mini-halos

We present a simulation of the formation of the earliest Population II stars, starting from cosmological initial conditions and ending when metals created in the first supernovae are incorporated into a collapsing gas-cloud. This occurs after a supernova blast-wave collides with a nearby mini-halo, inducing further turbulence that efficiently mixes metals into the dense gas in the center of the halo. The gas that first collapses has been enriched to a metallicity of Z ~ 2e-5 Zsun. Due to the extremely low metallicity, collapse proceeds similarly to metal-free gas until dust cooling becomes efficient at high densities, causing the cloud to fragment into a large number of low mass objects. This external enrichment mechanism provides a plausible origin for the most metal-poor stars observed, such as SMSS J031300.36-670839.3, that appear to have formed out of gas enriched by a single supernova. This mechanism operates on shorter timescales than the time for low-mass mini-halos (M < 5e5 Msun) to recover their gas after experiencing a supernova. As such, metal-enriched stars will likely form first via this channel if the conditions are right for it to occur. We identify a number of other externally enriched halos that may form stars in this manner. These halos have metallicities as high as 0.01 Zsun, suggesting that some members of the first generation of metal-enriched stars may be hiding in plain sight in current stellar surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The impact of reionization on the formation of supermassive black hole seeds

Direct collapse black holes (DCBHs) formed from the collapse of atomically-cooled primordial gas in the early Universe are strong candidates for the seeds of supermassive BHs. DCBHs are thought to form in atomic cooling haloes in the presence of a strong molecule-dissociating, Lyman-Werner (LW) radiation field. Given that star forming galaxies are likely to be the source of the LW radiation in this scenario, ionizing radiation from these galaxies may accompany the LW radiation. We present cosmological simulations resolving the collapse of primordial gas into an atomic cooling halo, including the effects of both LW and ionizing radiation. We find that in cases where the gas is not self-shielded from the ionizing radiation, the collapse can be delayed by ~ 25 Myr. When the ionized gas does collapse, the free electrons that are present catalyze H2 formation. In turn, H2 cooling becomes efficient in the center of the halo, and DCBH formation is prevented. We emphasize, however, that in many cases the gas collapsing into atomic cooling haloes at high redshift is self-shielding to ionizing radiation. Therefore, it is only in a fraction of such haloes in which DCBH formation is prevented due to reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Angular momentum loss of primordial gas in Lyman-alpha radiation field

We present results on the radiation drag exerted by an isotropic and homogenous background of Lya photons on neutral gas clouds orbiting within HII regions around Population III stars of different masses. The Doppler shift causes a frequency difference between photons moving in the direction of the cloud and opposite to it resulting in a net momentum loss of the cloud in the direction of motion. We find that half of the angular momentum of gas with v_theta <~ 20 km/s near (r <~ 3 kpc) a Population III star of 120 Msun at z=20 is lost within ~ 10^6 yr. The radiation drag is a strong function of cloud velocity that peaks at v ~ 20 km/s reflecting the frequency dependence of the photon cross section. Clouds moving with velocities larger than ~ 100 km/s loose their angular momentum on time scales of ~ 10^8 yr. At lower redshifts radiation drag becomes inefficient as the Lya photon density in HII regions decreases by a factor (1+z)^3 and angular momentum is lost on time scales > 10^8 yr even for low velocity clouds. Our results suggest that a sweet spot exists for the loss of angular momentum by radiation drag for gas clouds at z > 10 and with v ~ 20 km/s. Comparison to dynamical friction forces acting on typical gas clouds suggest that radiation drag is the dominant effect impacting the orbit. We propose that this effect can suppress the formation of extended gas discs in the first galaxies and help gas accretion near galactic centres and central black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Connection between dynamically derived IMF normalisation and stellar population parameters

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Revised rate coefficients for H$_2$ and H$^-$ destruction by realistic stellar spectra

Understanding the processes that can destroy H$_2$ and H$^-$ species is quintessential in governing the formation of the first stars, black holes and galaxies. In this study we compute the reaction rate coefficients for H$_2$ photo--dissociation by Lyman--Werner photons ($11.2 - 13.6$ eV), and H$^-$ photo--detachment by 0.76 eV photons emanating from self-consistent stellar populations that we model using publicly available stellar synthesis codes. So far studies that include chemical networks for the formation of molecular hydrogen take these processes into account by assuming that the source spectra can be approximated by a power-law dependency or a black-body spectrum at 10$^4$ or $10^5$ K. We show that using spectra generated from realistic stellar population models can alter the reaction rates for photo-dissociation, $\rm k_{\rm{di}}$, and photo-detachment, $\rm k_{\rm{de}}$, significantly. In particular, $\rm k_{\rm{de}}$ can be up to $\sim 2-4$ orders of magnitude lower in the case of realistic stellar spectra suggesting that previous calculations have over-estimated the impact that radiation has on lowering H$_2$ abundances. In contrast to burst modes of star formation, we find that models with continuous star formation predict increasing $\rm k_{\rm{de}}$ and $\rm k_{\rm{di}}$, which makes it necessary to include the star formation history of sources to derive self-consistent reaction rates, and that it is not enough to just calculate J$_{21}$ for the background. For models with constant star formation rate the change in shape of the spectral energy distribution leads to a non-negligible late-time contribution to $\rm k_{\rm{de}}$ and $\rm k_{\rm{di}}$, and we present self-consistently derived cosmological reaction rates based on star formation rates consistent with observations of the high redshift Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Atlas3D project - XXIV. The intrinsic shape distribution of early-type galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Atlas3D project - XXVI. HI discs in real and simulated fast and slow rotators

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The ATLAS3D Project - XXVIII. Dynamically-driven star formation suppression in early-type galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Atlas3D project -- XXIX. The new look of early-type galaxies and surrounding fields disclosed by extremely deep optical images

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The First Billion Years project: dark matter haloes going from contraction to expansion and back again

We study the effect of baryons on the inner dark matter profile of the first galaxies using the First Billion Years simulation between z=16-6 before secular evolution sets in. Using a large statistical sample from two simulations of the same volume and cosmological initial conditions, one with and one without baryons, we are able to directly compare haloes with their baryon-free counterparts, allowing a detailed study of the modifications to the dark matter density profile due to the presence of baryons during the first billion years of galaxy formation. For each of the ~ 5000 haloes in our sample we quantify the impact of the baryons using eta, defined as the ratio of dark matter mass enclosed in 100 pc in the baryonic run to its counterpart without baryons. During this epoch of rapid growth of galaxies, we find that many haloes of these first galaxies show an enhancement of dark matter in the halo centre compared to the baryon-free simulation, while many others show a deficit. We find that the mean value of eta is close to unity, but there is a large dispersion, with a standard deviation of 0.677. The enhancement is cyclical in time and tracks the star formation cycle of the galaxy; as gas falls to the centre and forms stars, the dark matter moves in as well. Supernova feedback then removes the gas, and the dark matter again responds to the changing potential. We study three physical models relating the motion of baryons to that of the dark matter: adiabatic contraction, dynamical friction, and rapid outflows. Abridged, see text for full abstract

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The HST/ACS Coma Cluster Survey - VII. Structure and Assembly of Massive Galaxies in the Center of the Coma Cluster

We constrain the assembly history of galaxies in the projected central 0.5 Mpc of the Coma cluster by performing structural decomposition on 69 massive (M_star >= 10^9 M_sun) galaxies using high-resolution F814W images from the HST Treasury Survey of Coma. Each galaxy is modeled with up to three Sersic components having a free Sersic index n. After excluding the two cDs in the projected central 0.5 Mpc of Coma, 57% of the galactic stellar mass in the projected central 0.5 Mpc of Coma resides in classical bulges/ellipticals while 43% resides in cold disk-dominated structures. Most of the stellar mass in Coma may have been assembled through major (and possibly minor) mergers. Hubble types are assigned based on the decompositions, and we find a strong morphology-density relation; the ratio of (E+S0):spirals is (91.0%):9.0%. In agreement with earlier work, the size of outer disks in Coma S0s/spirals is smaller compared with lower-density environments captured with SDSS (Data Release 2). Among similar-mass clusters from a hierarchical semi-analytic model, no single cluster can simultaneously match all the global properties of the Coma cluster. The model strongly overpredicts the mass of cold gas and underpredicts the mean fraction of stellar mass locked in hot components over a wide range of galaxy masses. We suggest that these disagreements with the model result from missing cluster physics (e.g., ram-pressure stripping), and certain bulge assembly modes (e.g., mergers of clumps). Overall, our study of Coma underscores that galaxy evolution is not solely a function of stellar mass, but also of environment.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

NGC1266 as a Local Candidate for Rapid Cessation of Star Formation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Atlas3D project - XIX. The hot-gas content of early-type galaxies: fast versus slow rotators

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]

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Atlas3D project - XV. Benchmark for early-type galaxies scaling relations from 260 dynamical models: mass-to-light ratio, dark matter, Fundamental Plane and Mass Plane

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]

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Atlas3D project - XX. Mass-size and mass-sigma distributions of early-type galaxies: bulge fraction drives kinematics, mass-to-light ratio, molecular gas fraction and stellar initial mass function

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]

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The ATLAS3D Project -- XXIII. Angular momentum and nuclear surface brightness profiles

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Unravelling obese black holes in the first galaxies

We predict the existence and observational signatures of a new class of objects that assembled early, during the first billion years of cosmic time: Obese Black-hole Galaxies (OBGs). OBGs are objects in which the mass of the central black hole initially exceeds that of the stellar component of the host galaxy, and the luminosity from black-hole accretion dominates the starlight. From a cosmological simulation, we demonstrate that there are sites where star formation is initially inhibited and direct-collapse black holes (DCBHs) form due to the photo-dissociating effect of Lyman-Werner radiation on molecular hydrogen. We show that the formation of OBGs is inevitable, because the probability of finding the required extra-galactic environment and the right physical conditions in a halo conducive to DCBH formation is quite high in the early universe. We estimate an OBG number density of 0.009/Mpc^3 at z~8 and 0.03/Mpc^3 at z~6. Extrapolating from our simulation volume, we infer that the most luminous quasars detected at z~6 likely transited through an earlier OBG phase. We find that these primordial galaxies start off with an over-massive BH and acquire their stellar component from subsequent merging as well as in-situ star formation. In doing so, they inevitably go through an OBG phase dominated by the accretion luminosity at the Eddington rate or below, released from the growing BH. The OBG phase is characterised by an ultra-violet (UV) spectrum with slope of beta ~ -2.3 and the absence of a Balmer Break. OBGs should also be spatially unresolved, and are expected to be brighter than the majority of known high-redshift galaxies. OBGs could potentially be revealed via HST follow-up imaging of samples of brighter Lyman-break galaxies provided by wide-area ground-based surveys such as UltraVISTA, and should be easily uncovered and studied with instruments aboard JWST...(abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A systematic variation of the stellar initial mass function in early-type galaxies

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Discovery of a giant HI tail in the galaxy group HCG 44

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gemini GMOS and WHT SAURON integral-field spectrograph observations of the AGN driven outflow in NGC 1266

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Hydrodynamical simulations and semi-analytic models of galaxy formation: two sides of the same coin

In this work we develop a new method to turn a state-of-the-art hydrodynamical cosmological simulation of galaxy formation (HYD) into a simple semi-analytic model (SAM). This is achieved by summarizing the efficiencies of accretion, cooling, star formation, and feedback given by the HYD, as functions of the halo mass and redshift. Surprisingly, by turning the HYD into a SAM, we conserve the mass of individual galaxies, with deviations at the level of 0.1 dex, on an object-by-object basis. This is true for all redshifts, and for the mass of stars and gas components, although the agreement reaches 0.2 dex for satellite galaxies at low redshift. We show that the same level of accuracy is obtained even in case the SAM uses only one phase of gas within each galaxy. Moreover, we demonstrate that the formation history of one massive galaxy provides sufficient information for the SAM to reproduce the population of galaxies within the entire cosmological box. The reasons for the small scatter between the HYD and SAM galaxies are: a) The efficiencies are matched as functions of the halo mass and redshift, meaning that the evolution within merger-trees agrees on average. b) For a given galaxy, efficiencies fluctuate around the mean value on time scales of 0.2-2 Gyr. c) The various mass components of galaxies are obtained by integrating the efficiencies over time, averaging out these fluctuations. We compare the efficiencies found here to standard SAM recipes and find that they often deviate significantly. For example, here the HYD shows smooth accretion that is less effective for low mass haloes, and is always composed of hot or dilute gas; cooling is less effective at high redshift; and star formation changes only mildly with cosmic time. The method developed here can be applied in general to any HYD, and can thus serve as a common language for both HYDs and SAM (Abridged).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Jet interactions with a giant molecular cloud in the Galactic centre and ejection of hypervelocity stars

The hypervelocity OB stars in the Milky Way Galaxy were ejected from the central regions some 10-100 million years ago. We argue that these stars, {as well as many more abundant bound OB stars in the innermost few parsecs,} were generated by the interactions of an AGN jet from the central black hole with a dense molecular cloud. Considerations of the associated energy and momentum injection have broader implications for the possible origin of the Fermi bubbles and for the enrichment of the intergalactic medium.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Jet-induced star formation in gas-rich galaxies

Feedback from active galactic nuclei (AGN) has become a major component in simulations of galaxy evolution, in particular for massive galaxies. AGN jets have been shown to provide a large amount of energy and are capable of quenching cooling flows. Their impact on the host galaxy, however, is still not understood. Subgrid models of AGN activity in a galaxy evolution context so far have been mostly focused on the quenching of star formation. To shed more light on the actual physics of the "radio mode" part of AGN activity, we have performed simulations of the interaction of a powerful AGN jet with the massive gaseous disc (10^11 solar masses) of a high-redshift galaxy. We spatially resolve both the jet and the clumpy, multi-phase interstellar medium (ISM) and include an explicit star formation model in the simulation. Following the system over more than 10^7 years, we find that the jet activity excavates the central region, but overall causes a significant change to the shape of the density probability distribution function and hence the star formation rate due to the formation of a blast wave with strong compression and cooling in the ISM. This results in a ring- or disc-shaped population of young stars. At later times, the increase in star formation rate also occurs in the disc regions further out since the jet cocoon pressurizes the ISM. The total mass of the additionally formed stars may be up to 10^10 solar masses for one duty cycle. We discuss the details of this jet-induced star formation (positive feedback) and its potential consequences for galaxy evolution and observable signatures.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Looking at the relation between haloes and galaxies under the lens

We use an empirical approach to model the stellar mass of galaxies according to their host dark-matter haloes and subhaloes ('HASH'), where each galaxy resides in a subhalo taken from a large N-body cosmological simulation. This approach allows us to study the mass relation between subhaloes and galaxies (MR) using: the weak lensing signal (WL), the two point auto-correlation function (CF), and the stellar mass function of galaxies (SMF). Our method is based on modeling the lensing signal directly from the cosmological simulation, and should thus be more accurate than other methods based on analytic halo models. We find that the WL does not provide a strong constraint on the MR, allowing for more than a factor of 10 freedom in the subhalo mass of central galaxies, for a given stellar mass. The freedom for satellite galaxies is much larger, providing a very poor constraint on the fraction of satellite galaxies (0.05-0.8). These results are not significantly modified when using both the SMF and WL as constraints. We show that for the most massive galaxies, observational constraints based on the CF with 0.1 dex errors, are equivalent to 0.05 dex error in the WL. For intermediate and low mass galaxies the WL and CF constrain the MR in a different way. Although the WL is currently not adding much information at these masses, it has the potential of being important using future, more accurate measurements. The models found here do not match simultaneously the observed CF and WL signals, and show a limited ability to match the WL & SMF. We suspect that this is partially due to the cosmological model assumed here and we therefore adopt a mock WL signal through most of this work. In comparison to previous models in the literature, the method presented here is probably more general, as it reveals a larger range of solutions for a given set of observational constraints [abridged].

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulations of binary galaxy mergers and the link with Fast Rotators, Slow Rotators, and Kinematically Distinct Cores

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Atlas3D project - XII. Recovery of the mass-to-light ratio of simulated early-type barred galaxies with axisymmetric dynamical models

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Atlas3D project - XIV. The extent and kinematics of molecular gas in early-type galaxies

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)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The ATLAS3D project - XVII. Linking photometric and kinematic signatures of stellar discs in early-type galaxies

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The ATLAS3D project - XVIII. CARMA CO imaging survey of early-type galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The ATLAS3D Project - XXI. Correlations between gradients of local escape velocity and stellar populations in early-type galaxies

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]

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The First Billion Years project - III: The impact of stellar radiation on the coevolution of Populations II and III

With the first metal enrichment by Population (Pop) III supernovae (SNe), the formation of the first metal-enriched, Pop II stars becomes possible. In turn, Pop III star formation and early metal enrichment are slowed by the high energy radiation emitted by Pop II stars. Thus, through the SNe and radiation they produce, Populations II and III coevolve in the early Universe, one regulated by the other. We present large (4 Mpc)^3, high resolution cosmological simulations in which we self-consistently model early metal enrichment and the stellar radiation responsible for the destruction of the coolants (H2 and HD) required for Pop III star formation. We find that the molecule-dissociating stellar radiation produced both locally and over cosmological distances reduces the Pop III star formation rate at z > 10 by up to an order of magnitude compared to the case in which this radiation is not included. However, we find that the effect of LW feedback is to enhance the amount of Pop II star formation. We attribute this to the reduced rate at which gas is blown out of dark matter haloes by SNe in the simulation with LW feedback, which results in larger reservoirs for metal-enriched star formation. Even accounting for metal enrichment, molecule-dissociating radiation and the strong suppression of low-mass galaxy formation due to reionization at z < 10, we find that Pop III stars are still formed at a rate of ~ 10^-5 M_sun yr^-1 Mpc^-3 down to z ~ 6. This suggests that the majority of primordial pair-instability SNe that may be uncovered in future surveys will be found at z < 10. We also find that the molecule-dissociating radiation emitted from Pop II stars may destroy H2 molecules at a high enough rate to suppress gas cooling and allow for the formation of supermassive primordial stars which collapse to form ~ 100,000 solar mass black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The First Billion Years project - IV: Proto-galaxies reionising the Universe

The contribution of stars in galaxies to cosmic reionisation depends on the star formation history in the Universe, the abundance of galaxies during reionisation, the escape fraction of ionising photons and the clumping factor of the inter-galactic medium (IGM). We compute the star formation rate and clumping factor during reionisation in a cosmological volume using a high-resolution hydrodynamical simulation. We post-process the output with detailed radiative transfer simulations to compute the escape fraction of ionising photons. Together, this gives us the opportunity to assess the contribution of galaxies to reionisation self-consistently. The strong mass and redshift dependence of the escape fraction indicates that reionisation occurred between z=15 and z=10 and was mainly driven by proto-galaxies forming in dark-matter haloes with masses between 1e7 and 1e8 solar mass. More massive galaxies that are rare at these redshifts and have significantly lower escape fractions contribute less photons to the reionisation process than the more-abundant low-mass galaxies. Star formation in the low-mass haloes is suppressed by radiative feedback from reionisation, therefore these proto-galaxies only contribute when the part of the Universe they live in is still neutral. After z~10, massive galaxies become more abundant and provide most of the ionising photons. In addition, we find that Population (Pop) III stars are too short-lived and not frequent enough to have a major contribution to reionisation. Although the stellar component of the proto-galaxies that produce the bulk of ionising photons during reionisation is too faint to be detected by the James Webb Space Telescope (JWST), these sources are brightest in the H-alpha and Ly-alpha recombination lines, which will likely be detected by JWST in deep surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Growth of the Stellar Seeds of Supermassive Black Holes

One of the most promising explanations for the origin of the billion solar mass black holes (BHs) inferred to power quasars at redshifts z > 6 is that supermassive stars (SMSs) with masses > 10,000 solar masses collapse to form the seed BHs from which they grow. Here we review recent theoretical advances which provide support for this scenario. Firstly, given sufficiently high accretion rates of gas into the cores of primordial protogalaxies, it appears that neither the high energy radiation emitted from the stellar surface nor the limited lifetime of SMSs can prevent their growth to masses of up to > 100,000 solar masses. Secondly, recent cosmological simulations suggest that the high fluxes of molecule-dissociating radiation which may be required in order to achieve such high accretion rates may be more common in the early universe than previously thought. We conclude that the majority of supermassive BHs may originate from SMSs at high redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The road to the red sequence: A detailed view of the formation of a massive galaxy at z~2

(Abridged) We present here a detailed analysis of the star formation history (SFH) of FW4871, a massive galaxy at z=1.893+-0.002. We compare rest-frame optical and NUV slitless grism spectra from the Hubble Space Telescope with a large set of composite stellar populations to constrain the underlying star formation history. Even though the morphology features prominent tidal tails, indicative of a recent merger, there is no sign of on-going star formation within an aperture encircling one effective radius, which corresponds to a physical extent of 2.6 kpc. A model assuming truncation of an otherwise constant SFH gives a formation epoch zF~10, with a truncation after 2.7 Gyr, giving a mass-weighted age of 1.5 Gyr and a stellar mass of 0.8-3E11Msun, implying star formation rates of 30-110 Msun/yr. A more complex model including a recent burst of star formation places the age of the youngest component at 145 Myr, with a mass contribution lower than 20%, and a maximum amount of dust reddening of E(B-V)<0.4 mag (95% confidence levels). This low level of dust reddening is consistent with the low emission observed at 24 micron, corresponding to rest-frame 8 micron, where PAH emission should contribute significantly if a strong formation episode were present. The colour profile of FW4871 does not suggest a significant radial trend in the properties of the stellar populations out to 3Re. We suggest that the recent merger that formed FW4871 is responsible for the quenching of its star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ubiquitous seeding of supermassive black holes by direct collapse

We study for the first time the environment of massive black hole (BH) seeds (~10^4-5 Msun) formed via the direct collapse of pristine gas clouds in massive haloes (>10^7 Msun) at z>6. Our model is based on the evolution of dark matter haloes within a cosmological N-body simulation, combined with prescriptions for the formation of BH along with both Pop III and Pop II stars. We calculate the spatially-varying intensity of Lyman Werner (LW) radiation from stars and identify the massive pristine haloes in which it is high enough to shut down molecular hydrogen cooling. In contrast to previous BH seeding models with a spatially constant LW background, we find that the intensity of LW radiation due to local sources, J_local, can be up to 10^6 times the spatially averaged background in the simulated volume and exceeds the critical value, J_crit, for the complete suppression of molecular cooling, in some cases by 4 orders of magnitude. Even after accounting for possible metal pollution in a halo from previous episodes of star formation, we find a steady rise in the formation rate of direct collapse (DC) BHs with decreasing redshift from 10^{-3}/Mpc^3/z at z=12 to 10^{-2}/Mpc^3/z at z=6. The onset of Pop II star formation at z~16 simultaneously marks the onset of the epoch of DCBH formation, as the increased level of LW radiation from Pop II stars is able to elevate the local levels of the LW intensity to J_local > J_crit while Pop III stars fail to do so at any time. The number density of DCBHs is sensitive to the number of LW photons and can vary by an order of magnitude at z=6 after accounting for reionisation feedback. Haloes hosting DCBHs are more clustered than similar massive counterparts that do not host DCBHs, especially at redshifts z>10. We also show that planned surveys with JWST should be able to detect the supermassive stellar precursors of DCBHs.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A tale of two populations: the stellar mass of central and satellite galaxies

We develop a new empirical methodology to study the relation between the stellar mass of galaxies and the mass of their host subhaloes. Our approach is similar to abundance matching, and is based on assigning a stellar mass to each subhalo within a large cosmological N-body simulation. However, we significantly extend the method to include a different treatment for central and satellite galaxies, allowing the stellar mass of satellite galaxies to depend on both the host subhalo mass, and on the halo mass. Unlike in previous studies, our models are constrained by two observations: the stellar mass function of galaxies, and their auto-correlation functions (CFs). We present results for ~10^6 different successful models, showing that the uncertainty in the host subhalo mass reaches a factor of ~10 for a given stellar mass. Our results thus indicate that the host subhalo mass of central and satellite galaxies is poorly constrained by using information coming solely from the abundance and CFs of galaxies. In addition, we demonstrate that the fraction of stellar mass relative to the universal baryon fraction can vary between different models, and can reach ~0.6 for a specific population of low mass galaxies. We conclude that additional observational constraints are needed, in order to better constrain the mass relation between haloes and galaxies. These might include weak lensing, satellite kinematics, or measures of clustering other than the CFs. Moreover, CFs based on larger surveys with a better estimate of the systematic uncertainties are needed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Discovery of an AGN-Driven Molecular Outflow in the Local Early-Type Galaxy NGC 1266

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gravity at Work: How the Build-Up of Environments Shape Galaxy Properties

We present results on the heating of the inter-cluster medium (ICM) by gravitational potential energy from in-falling satellites. We calculate the available excess energy of baryons once they are stripped from their satellite and added to the ICM of the hosting environment. this excess energy is a strong function of environment and we find that it can exceed the contribution from AGNs or supernovae (SN) by up to two orders of magnitude in the densest environments/haloes. Cooling by radiative losses is in general fully compensated by gravitational heating in massive groups and clusters with hot gas temperature > 1 keV. The reason for the strong environment dependence is the continued infall of substructure onto dense environments during their formation in contrast to field-like environments. We show that gravitational heating is able to reduce the number of too luminous galaxies in models and to produce model luminosity functions in agreement with observations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Atlas-3D project - IX. The merger origin of a fast and a slow rotating Early-Type Galaxy revealed with deep optical imaging: first results

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The ATLAS3D project - II. Morphologies, kinemetric features and alignment between photometric and kinematic axes of early-type galaxies

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The ATLAS3D project - III. A census of the stellar angular momentum within the effective radius of early-type galaxies: unveiling the distribution of Fast and Slow Rotators

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Atlas3D project - V. The CO Tully-Fisher relation of early-type galaxies

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Atlas3D Project - VI. Simulations of binary galaxy mergers and the link with Fast Rotators, Slow Rotators, and Kinematically Distinct Cores

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Atlas3D project - VII. A new look at the morphology of nearby galaxies: the kinematic morphology-density relation

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Atlas3D project - X. On the origin of the molecular and ionised gas in early-type galaxies

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Atlas3D project -- IV: the molecular gas content of early-type galaxies

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Atlas3D Project -- XI. Dense molecular gas properties of CO-luminous early-type galaxies

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Atlas3D project -- XIII. Mass and morphology of HI in early-type galaxies as a function of environment

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The ATLAS3D Project-- VIII: Modelling the Formation and Evolution of Fast and Slow Rotator Early-Type Galaxies within $Λ$CDM

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Contribution of Supernovae to Cosmic Reionization

While stars are widely discussed as the source of the high energy photons which reionized the universe, an additional source of ionizing photons that must also contribute to reionization in this scenario is the supernovae (SNe) which mark the end of the life of massive stars. Here we estimate the relative contributions of SNe and stars to reionization. While the rate at which ionizing photons are produced in SN shocks is well below that at which they are produced by stars, the harder spectra of radiation emitted from SNe leads to an enhanced escape fraction of SN generated photons relative to that of stellar photons. In particular, along a given line of sight out of a galaxy, we find that for neutral hydrogen column densities N_H >~10^18 cm^-2 the contribution to reionization from SNe is greater than that from stars. Drawing on the results of simulations presented in the literature, we find that the overall (line of sight-averaged) SNe shock-generated ionizing photon escape fraction is larger than the stellar photon escape fraction by a factor of \simeq 4 to \simeq 7, depending on the metallicity of the stellar population. Overall, our results suggest that the effect of SNe is an enhancement of up to ~ 10 percent in the fraction of hydrogen reionized by stellar sources. We briefly discuss the implications of our results for the population of galaxies responsible for reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The imprint of cosmological non-Gaussianities on primordial structure formation

We study via numerical N-body/SPH chemistry simulations the effects of primordial non-Gaussianities on the formation of the first stars and galaxies, and investigate the impact of supernova feedback in cosmologies with different fnl. Density distributions are biased to higher values, so star formation and the consequent feedback processes take place earlier in high-fnl models and later in low-fnl ones. Mechanical feedback is responsible for shocking and evacuating the gas from star forming sites earlier in the highly non-Gaussian cases, because of the larger bias at high densities. Chemical feedback translates into high-redshift metal filling factors that are larger by some orders of magnitude for larger fnl, but that converge within one Gyr, for both population III and population II-I stellar regimes. The efficient enrichment process, though, leads to metallicities > 0.01 Zsun by redshift ~9, almost independently from fnl. The impact of non-Gaussianities on the formation of dark-matter haloes at high redshift is directly reflected in the properties of the gas in these haloes, as models with larger fnl show more concentrated gas profiles at early times. Non-Gaussian signatures in the gas behaviour are lost after the first feedback takes place and introduces a significant degree of turbulence and chaotic motions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The interplay between chemical and mechanical feedback from the first generation of stars

We study cosmological simulations of early structure formation, including non-equilibrium molecular chemistry, metal pollution from stellar evolution, transition from population III (popIII) to population II (popII) star formation, regulated by a given critical metallicity, and feedback effects. We investigate the properties of early metal spreading from the different stellar populations and its interplay with primordial molecular gas. We find that, independently of the details about popIII modeling, after the onset of star formation, regions enriched below the critical level are mostly found in isolated environments, while popII star formation regions are much more clumped. Typical star forming haloes show average SN driven outflow rates of up to 10^{-4} Msun/yr in enriched gas, initially leaving the original star formation regions almost devoid of metals. The polluted material, which is gravitationally incorporated in over-dense environments on timescales of 10^7 yr, is mostly coming from external, nearby star forming sites ("gravitational enrichment"). In parallel, the pristine-gas inflow rates are between 10^{-3} - 10^{-1} Msun/yr. However, thermal feedback from SN generates turbulence and destroys molecules within the pristine gas, and only the polluted material, incorporated via gravitational enrichment, can continue to cool by atomic metal fine-structure transitions on time scales short enough to end the initial popIII regime within less than 10^8 yr.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Asymmetries in Extragalactic Double Radio Sources: Clues from 3D Simulations of Jet - Disc Interaction

Observational and theoretical studies of extragalactic radio sources have suggested that an inhomogeneous environment may be responsible for observed arm length asymmetries of jets and the properties of extended emission line regions in high redshift radio galaxies. We perform 3D hydrodynamic simulations of the interaction of a powerful extragalactic bipolar jet with a disc-shaped clumpy interstellar medium of log-normal density distribution and analyze the asymmetry. Furthermore, we compute the relation between jet asymmetry and the ISM properties by means of Monte Carlo simulations based on a 1D propagation model for the jet through the dense medium. We find that the properties of the ISM can be related to a probability distribution of jet arm length asymmetries: Disc density and height are found to have the largest effect on the asymmetry for realistic parameter ranges, while the Fourier energy spectrum of the ISM and turbulent Mach number only have a smaller effect. The hydrodynamic simulations show that asymmetries generally may be even larger than expected from the 1D model due to the complex interaction of the jet and its bow shock with gaseous clumps, which goes much beyond simple energy disposal. From our results, observed asymmetries of medium-sized local radio galaxies may be explained by gas masses of 10^9 to 10^10 solar masses in massive elliptical galaxies. Furthermore, the simulations provide a theoretical basis for the observed correlation that emission line nebulae are generally found to be brighter on the side of the shorter lobe in high redshift radio galaxies (McCarthy et al. 1991). This interaction of jets with the cold gas phase suggests that star formation in evolving high redshift galaxies may be affected considerably by jet activity.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Further constraining galaxy evolution models through the Size Function of SDSS Early-type galaxies

We discuss how the effective radius Phi(Re) function (ERF) recently worked out by Bernardi et al. (2009) represents a new testbed to improve the current understanding of Semi-analytic Models of Galaxy formation. In particular, we here show that a detailed hierarchical model of structure formation can broadly reproduce the correct peak in the size distribution of local early-type galaxies, although it significantly overpredicts the number of very compact and very large galaxies. This in turn is reflected in the predicted size-mass relation, much flatter than the observed one, due to too large (~3 kpc) low-mass galaxies (<10^11 \msun), and to a non-negligible fraction of compact (< 0.5-1 kpc) and massive galaxies (> 10^11 \msun). We also find that the latter discrepancy is smaller than previously claimed, and limited to only ultracompact (Re < 0.5 kpc) galaxies when considering elliptical-dominated samples. We explore several causes behind these effects. We conclude that the former problem might be linked to the initial conditions, given that large and low-mass galaxies are present at all epochs in the model. The survival of compact and massive galaxies might instead be linked to their very old ages and peculiar merger histories. Overall, knowledge of the galactic stellar mass {\em and} size distributions allows a better understanding of where and how to improve models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Mergers and Bulge Formation in Lambda-CDM: Which Mergers Matter?

We use a suite of semi-empirical models to predict galaxy merger rates and contributions to bulge growth as functions of merger mass, redshift, and mass ratio. The models use empirical halo occupation constraints to identify mergers, together with high-resolution simulations to quantify how mergers with different properties contribute to the bulge population. We find good agreement with a variety of observational constraints, and provide fitting functions for merger rates and contributions to bulge growth. We identify several robust conclusions. (1) Major mergers dominate formation and assembly of L* bulges and the spheroid mass density, minor mergers contribute ~30%. (2) This is mass-dependent: bulge formation is dominated by more minor mergers in lower-mass systems. At higher masses, bulges form in major mergers near L*, but subsequently assemble in minor mergers. (3) The minor/major contribution is also morphology-dependent: higher B/T systems form in more major mergers, lower B/T systems form in situ from minor mergers. (4) Low-mass galaxies, being gas-rich, require more major mergers to reach the same B/T as high-mass systems. (5) Absolute merger rates increase with galaxy mass. (6) Predicted rates agree well with observations, but suggest that some morphology-selected samples include contamination from minor mergers. (7) Predicted rates agree with integrated growth in bulge mass with cosmic time, but with factor ~2 uncertainty - half the bulge mass density could come from non-mergers. We consider ~1000 model variations and quantify resulting uncertainties. Conclusions regarding the major/minor contribution to bulge growth are very robust, absolute merger rates have systematic factor ~2 uncertainties.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Mergers in Lambda-CDM: Uncertainties in Theoretical Predictions and Interpretations of the Merger Rate

Different methodologies lead to order-of-magnitude variations in predicted galaxy merger rates. We examine and quantify the dominant uncertainties. Different halo merger rates and subhalo 'destruction' rates agree to within a factor ~2 given proper care in definitions. If however (sub)halo masses are not appropriately defined or are under-resolved, the major merger rate can be dramatically suppressed. The dominant differences in galaxy merger rates owe to baryonic physics. Hydrodynamic simulations without feedback and older models that do not agree with the observed galaxy mass function propagate factor ~5 bias in the resulting merger rates. However, if the model matches the galaxy mass function, properties of central galaxies are sufficiently converged to give small differences in merger rates. But variations in baryonic physics of satellites also have dramatic effects. The known problem of satellite 'over-quenching' in most semi-analytic models (SAMs), whereby SAM satellites are too efficiently stripped of gas, could lead to order-of-magnitude under-estimates of merger rates for low-mass, gas-rich galaxies. Fixing the satellite properties to observations tends to predict higher merger rates, but with factor ~2 empirical uncertainties. Choice of mass ratio definition matters: at low masses, most true major mergers (in baryonic/dynamical galaxy mass) will appear to be minor mergers in their stellar or luminosity mass ratio. Observations and models using these criteria may underestimate major merger rates by factors ~5. Orbital parameters and gas fractions also introduce factor ~3 differences in amount of bulge formed by mergers, even for fixed mass ratio encounters.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Suppression of accretion onto low-mass Population III stars

Motivated by recent theoretical work suggesting that a substantial fraction of Population (Pop) III stars may have had masses low enough for them to survive to the present day, we consider the role that the accretion of metal-enriched gas may have had in altering their surface composition, thereby disguising them as Pop II stars. We demonstrate that if weak, Solar-like winds are launched from low-mass Pop III stars formed in the progenitors of the dark matter halo of the Galaxy, then such stars are likely to avoid significant enrichment via accretion of material from the interstellar medium. We find that at early times accretion is easily prevented if the stars are ejected from the central regions of the haloes in which they form, either by dynamical interactions with more massive Pop III stars, or by violent relaxation during halo mergers. While accretion may still take place during passage through sufficiently dense molecular clouds at later times, we find that the probability of such a passage is generally low (< 0.1), assuming that stars have velocities of order the maximum circular velocity of their host haloes and accounting for the orbital decay of merging haloes. In turn, due to the higher gas density required for accretion onto stars with higher velocities, we find an even lower probability of accretion (~ 0.01) for the subset of Pop III stars formed at z > 10, which are more quickly incorporated into massive haloes than stars formed at lower redshift. While there is no a priori reason to assume that low-mass Pop III stars do not have Solar-like winds, without them surface enrichment via accretion is likely to be inevitable. We briefly discuss the implications that our results hold for stellar archaeology.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Atlas3D project -- I. A volume-limited sample of 260 nearby early-type galaxies: science goals and selection criteria

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Impact of Gas Stripping and Stellar Mass Loss on Satellite Galaxy Evolution

Current semi-analytic models (SAMs) of galaxy formation over-predict the fraction of passive small late-type satellite galaxies in dense environments by a factor of two to three. We hypothesize that this is due to inaccurate prescriptions on cold gas evolution. In the hope of solving this problem we apply detailed prescriptions on the evolution of diffuse hot gases in satellites and on stellar mass loss, both of which are critical to model cold gas evolution. We replace the conventional shock-heating motivated instant stripping with a realistic gradual prescription based on ram pressure and tidal stripping. We also carefully consider stellar mass loss in our model. When both mechanisms are included, the fraction of passive late types matches the data much more closely. The satellite over-quenching problem is still present in small galaxies in massive haloes, however. In terms of the detectable residual star formation rates, gradual diffuse gas stripping appears to be much more important than stellar mass loss in our model. The implications of these results and other possibilities, such as redshift-dependent merging geometry and tidal disruption, are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Bulge n and B/T in High Mass Galaxies: Constraints on the Origin of Bulges in Hierarchical Models

We use the bulge Sersic index n and bulge-to-total ratio (B/T) to explore the fundamental question of how bulges form. We perform 2D bulge-disk-bar decomposition on H-band images of 143 bright, high stellar mass (>1.0e10 solar masses) low-to-moderately inclined (i<70 degrees) spirals. Our results are: (1) Our H-band bar fraction (~58%) is consistent with that from ellipse fits. (2) 70% of the stellar mass is in disks, 10% in bars, and 20% in bulges. (3) A large fraction (~69%) of bright spirals have B/T<0.2, and ~76% have low n<2 bulges. These bulges exist in barred and unbarred galaxies across a wide range of Hubble types. (4) About 65% (68%) of bright spirals with n<2 (B/T<0.2) bulges host bars, suggesting a possible link between bars and bulges. (5) We compare the results with predictions from a set of LCDM models. In the models, a high mass spiral can have a bulge with a present-day low B/T<0.2 only if it did not undergo a major merger since z<2. The predicted fraction (~1.6%) of high mass spirals, which have undergone a major merger since z<4 and host a bulge with a present-day low B/T<0.2, is a factor of over thirty smaller than the observed fraction (~66%) of high mass spirals with B/T<0.2. Thus, contrary to common perception, bulges built via major mergers since z<4 seriously fail to account for the bulges present in ~66% of high mass spirals. Most of these present-day low B/T<0.2 bulges are likely to have been built by a combination of minor mergers and/or secular processes since z<4.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Galaxy Zoo: A sample of blue early-type galaxies at low redshift

We report the discovery of a population of nearby, blue early-type galaxies with high star formation rates (0.5 < SFR < 50 Msun/yr). They are identified by their visual morphology as provided by Galaxy Zoo for SDSS DR6 and their u-r colour. We select a volume-limited sample in the redshift range 0.02 < z < 0.05, corresponding to luminosities of approximately L* and above, and with u-r colours significantly bluer than the red sequence. We confirm the early-type morphology of the objects in this sample and investigate their environmental dependence and star formation properties. Blue early-type galaxies tend to live in lower-density environments than `normal' red sequence early-types and make up 5.7 +/-0.4% of the low-redshift early-type galaxy population. We find that such blue early-type galaxies are virtually absent at high velocity dispersions above 200 km/s. Our analysis uses emission line diganostic diagrams and we find that ~25% of them are actively starforming, while another ~25% host both star formation and an AGN. Another ~12% are AGN. The remaining 38% show no strong emission lines. When present and uncontaminated by an AGN contribution, the star formation is generally intense. We consider star formation rates derived from Halpha, u-band and infrared luminosities, and radial colour profiles, and conclude that the star formation is spatially extended. Of those objects that are not currently undergoing star formation must have ceased doing so recently in order to account for their blue optical colours. The gas phase metallicity of the actively starforming blue early-types galaxies is supersolar in all cases. We discuss the place of these objects in the context of galaxy formation. A catalogue of all 204 blue early-type galaxies in our sample, including star formation rates and emission line classification, is provided.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the formation of massive galaxies: A simultaneous study of number density, size and intrinsic colour evolution in GOODS

The evolution of number density, size and intrinsic colour is determined for a volume-limited sample of visually classified early-type galaxies selected from the HST/ACS images of the GOODS North and South fields (version 2). The sample comprises 457 galaxies over 320 arcmin2 with stellar masses above 3E10 Msun in the redshift range 0.4<z<1.2. Our data allow a simultaneous study of number density, intrinsic colour distribution and size. We find that the most massive systems (>3E11 Msun) do not show any appreciable change in comoving number density or size in our data. Furthermore, when including the results from 2dFGRS, we find that the number density of massive early-type galaxies is consistent with no evolution between z=1.2 and 0, i.e. over an epoch spanning more than half of the current age of the Universe. Massive galaxies show very homogeneous **intrinsic** colour distributions, featuring red cores with small scatter. The distribution of half-light radii -- when compared to z=0 and z>1 samples -- is compatible with the predictions of semi-analytic models relating size evolution to the amount of dissipation during major mergers. However, in a more speculative fashion, the observations can also be interpreted as weak or even no evolution in comoving number density **and size** between 0.4<z<1.2, thus pushing major mergers of the most massive galaxies towards lower redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

A simple model for the size-evolution of elliptical galaxies

We use semi-analytical modelling of galaxy formation to predict the redshift-size-evolution of elliptical galaxies. Using a simple model in which relative sizes of elliptical galaxies of a given mass correlate with the fraction of stars formed in a star burst during a major merger event, we are able to reproduce the observed size redshif t evolution. The size evolution is a result of the amount of cold gas available during the major merger. Mergers at high redshifts are gas-rich and produce ellipticals with smaller sizes. In particular we find a power-law relation between the sizes at different redshifts, with the power-law index giving a measure of the relative amount of dissipation during the mergers that lead to the formation of an elliptical. The size evolution is found to be stronger for more massive galaxies as they involve more gas at high redshifts when they form compared to less massive ellipticals. Local ellipticals more massive than $5 \times 10^{11}$ M$_{\odot}$ will be approximately 4 times larger than their counterparts at $z=2$. Our results indicate that the scatter in the size of similar massive present day elliptical galaxies is a result of their formation epoch, with smaller ellipticals being formed earlier.

0 citations0 reviewsNo code linkedOpen access