Cold Gas Fractions Research Articles

STARBURST-DRIVEN GALACTIC SUPERBUBBLES RADIATING TO 10 K

ABSTRACT Our three-dimensional hydrodynamical simulations of starbursts examine the formation of superbubbles over a range of driving luminosities and mass loadings that determine superbubble growth and wind velocity. From this we determine the relationship between the velocity of a galactic wind (GW) and the power of the starburst. We find a threshold for the formation of a wind, above which the speed of the wind is not affected by grid resolution or the temperature floor of our radiative cooling. We investigate the effect that two different temperature floors in our radiative cooling prescription have on wind kinematics and content. We find that cooling to 10 K instead of to 104 K increases the mass fraction of cold neutral and hot X-ray gas in the GW, while halving that in warm Hα. Our simulations show that the mass of cold gas transported into the lower halo does not depend on the starburst strength. Optically bright filaments form at the edge of merging superbubbles, or where a cold dense cloud has been disrupted by the wind. Filaments formed by merging superbubbles will persist and grow to > 400 pc in length if anchored to a star forming complex. Filaments embedded in the hot GW contain warm and cold gas that moves 300−1200 km s−1 slower than the surrounding wind, with the coldest gas hardly moving with respect to the Galaxy. Warm and cold matter in the GW show asymmetric absorption profiles consistent with observations, with a thin tail up to the wind velocity.

Star formation and molecular hydrogen in dwarf galaxies: a non-equilibrium view

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

Properties of damped Ly α absorption systems and star-forming galaxies in semi-analytic models atz = 2

We investigate predictions from semi-analytic cosmological models of galaxy formation for the properties of star-forming galaxies (SFGs) and damped Ly α absorption systems (DLAS), and the relationship between these two populations. Our models reproduce fairly well the observed distributions of redshift, stellar mass, star formation rate (SFR), and dust extinction for z ˜ 2 SFGs. We predict that DLA hosts span a broad range of properties, with broad and relatively flat distributions of stellar and halo mass, SFR, and luminosity. The photometric colours of DLA host galaxies trace the colours of galaxies with similar luminosities, but the majority are much fainter than the limits of most existing surveys of SFGs. Generally, DLA host galaxies and SFGs at z = 2 follow similar trends between stellar mass, DLA cross-section, cold gas fraction, SFR, metallicity, and dust extinction as the global population of galaxies with the same stellar mass. Since DLAS select galaxies with larger cold gas masses, they tend to have larger cold gas fractions, lower metallicities, higher SFRs, and less dust extinction than galaxies at the same stellar mass. Our models reproduce the observed relations between impact parameter, column density, and metallicity, suggesting that the sizes of the gas discs giving rise to DLAS in our models are roughly correct. We find that molecular fractions and SFRs are in general significantly lower at the location of the DLA line of sight than the galaxy-averaged value.

The inferred evolution of the cold gas properties of CANDELS galaxies at 0.5 &lt;z&lt; 3.0

We derive the total cold gas, atomic hydrogen, and molecular gas masses of approximately 24 000 galaxies covering four decades in stellar mass at redshifts 0.5 < z < 3.0, taken from the CANDELS survey. Our inferences are based on the inversion of a molecular hydrogen based star formation law, coupled with a prescription to separate atomic and molecular gas. We find that: 1) there is an increasing trend between the inferred cold gas (HI and H2), HI, and H2 mass and the stellar mass of galaxies down to stellar masses of 10^8 Msun already in place at z = 3; 2) the molecular fractions of cold gas increase with increasing stellar mass and look-back time; 3) there is hardly any evolution in the mean HI content of galaxies at fixed stellar mass; 4) the cold gas fraction and relative amount of molecular hydrogen in galaxies decrease at a relatively constant rate with time, independent of stellar mass; 5) there is a large population of low-stellar mass galaxies dominated by atomic gas. These galaxies are very gas rich, but only a minor fraction of their gas is molecular; 6) the ratio between star-formation rate (SFR) and inferred total cold gas mass (HI + H2) of galaxies (i.e., star-formation efficiency; SFE) increases with star-formation at fixed stellar masses. Due to its simplicity, the presented approach is valuable to assess the impact of selection biases on small samples of directly-observed gas masses and to extend scaling relations down to stellar mass ranges and redshifts that are currently difficult to probe with direct measurements of gas content.

The impact of stellar feedback on high-z galaxy populations

AbstractOne major deficiency of state-of-the-art galaxy formation models consists in their inability of capturing the observed galaxy downsizing trend significantly over-estimating the number density of low-mass galaxies, in particular at high redshifts. Employing an enhanced galaxy formation model with a full chemical enrichment scheme (DeLucia et al., 2014), we present an improved model for stellar feedback (based on parametrizations from cosmological zoom simulations), in which strong gas outflows occur due to bursty star formation at high z, while star formation is mainly “quiescent” not causing any significant outflows anymore at low z. Due to the stronger gas outflows at high z, early star formation is strongly delayed towards later times. This helps to sufficiently detach the evolution of galaxy growth from the hiearchical dark matter assembly resulting in a fairly good agreement with the evolution of the observed stellar mass function (SMF, see Fig. 1). With our new feedback scheme, we can also successfully reproduce many other observational constraints, such as the metallicity content, the cold gas fractions or the quiescent galaxy fractions at both low and high redshifts. The resulting new-generation galaxy catalogues (Hirschmann et al., in prep) based on that model are expected to significantly contribute to the interpretation of current and up-coming large-scale surveys (HST, JWST, Euclid). This will, in turn, provide a rapid verification and refinement of our modeling.

AN ANALYTICAL MODEL FOR GALAXY METALLICITY: WHAT DO METALLICITY RELATIONS TELL US ABOUT STAR FORMATION AND OUTFLOW?

We develop a simple analytical model that tracks galactic metallicities governed by star formation and feedback to gain insight from the observed galaxy stellar mass-metallicity relations over a large range of stellar masses and redshifts. The model reveals the following implications of star formation and feedback processes in galaxy formation. First, the observed metallicity relations provide a stringent upper limit for the averaged outflow mass-loading factors of local galaxies, which is ~20 for M_*~10^9Msun galaxies and monotonically decreases to ~1 for M_*~10^{11}Msun galaxies. Second, the inferred upper-limit for the outflow mass-loading factor sensitively depends on whether the outflow is metal-enriched with respect to the ISM metallicity. If half of the metals ejected from SNe leave the galaxy in metal-enriched winds, the outflow mass-loading factor for galaxies at any mass can barely be higher than ~10, which puts strong constraints on galaxy formation models. Third, the relatively lower stellar-phase to gas-phase metallicity ratio for lower-mass galaxies indicate that low-mass galaxies are still rapidly enriching their metallicities in recent times, while high-mass galaxies are more settled, which seems to show a downsizing effect in the metallicity evolution of galaxies. The analysis presented in the paper demonstrates the importance of accurate measurements of galaxy metallicities and the cold gas fraction of galaxies at different redshifts for constraining star formation and feedback processes, and demonstrates the power of these relations in constraining the physics of galaxy formation.

THE 21-SPONGE H i ABSORPTION SURVEY. I. TECHNIQUES AND INITIAL RESULTS

We present methods and results from "21-cm Spectral Line Observations of Neutral Gas with the EVLA" (21-SPONGE), a large survey for Galactic neutral hydrogen (HI) absorption with the Karl G. Jansky Very Large Array (VLA). With the upgraded capabilities of the VLA, we reach median root-mean-square (RMS) noise in optical depth of $\sigma_{\tau}=9\times 10^{-4}$ per $0.42\rm\,km\,s^{-1}$ channel for the 31 sources presented here. Upon completion, 21-SPONGE will be the largest HI absorption survey with this high sensitivity. We discuss the observations and data reduction strategies, as well as line fitting techniques. We prove that the VLA bandpass is stable enough to detect broad, shallow lines associated with warm HI, and show that bandpass observations can be combined in time to reduce spectral noise. In combination with matching HI emission profiles from the Arecibo Observatory ($\sim3.5'$ angular resolution), we estimate excitation (or spin) temperatures ($\rm T_s$) and column densities for Gaussian components fitted to sightlines along which we detect HI absorption (30/31). We measure temperatures up to $\rm T_s\sim1500\rm\,K$ for individual lines, showing that we can probe the thermally unstable interstellar medium (ISM) directly. However, we detect fewer of these thermally unstable components than expected from previous observational studies. We probe a wide range in column density between $\sim10^{16}$ and $>10^{21}\rm\,cm^{-2}$ for individual HI clouds. In addition, we reproduce the trend between cold gas fraction and average $\rm T_s$ found by synthetic observations of a hydrodynamic ISM simulation by Kim et al. (2014). Finally, we investigate methods for estimating HI $\rm T_s$ and discuss their biases.

Evolution of the atomic and molecular gas content of galaxies in dark matter haloes

We present a semi-empirical model to infer the atomic and molecular hydrogen content of galaxies as a function of halo mass and time. Our model combines the SFR-halo mass-redshift relation (constrained by galaxy abundances) with inverted SFR-surface density relations to infer galaxy H I and H2 masses. We present gas scaling relations, gas fractions, and mass functions from z = 0 to z = 3 and the gas properties of galaxies as a function of their host halo masses. Predictions of our work include: 1) there is a ~ 0.2 dex decrease in the H I mass of galaxies as a function of their stellar mass since z = 1.5, whereas the H2 mass of galaxies decreases by > 1 dex over the same period. 2) galaxy cold gas fractions and H2 fractions decrease with increasing stellar mass and time. Galaxies with M* > 10^10 Msun are dominated by their stellar content at z < 1, whereas less-massive galaxies only reach these gas fractions at z = 0. We find the strongest evolution in relative gas content at z < 1.5. 3) the SFR to gas mass ratio decreases by an order of magnitude from z = 3 to z = 0. This is consistent with lower H2 fractions; these lower fractions in combination with smaller gas reservoirs correspond to decreased present-day galaxy SFRs. 4) an H2-based star- formation relation can simultaneously fuel the evolution of the cosmic star-formation and reproduce the observed weak evolution in the cosmic HI density. 5) galaxies residing in haloes with masses near 10^12 Msun are most efficient at obtaining large gas reservoirs and forming H2 at all redshifts. These two effects lie at the origin of the high star-formation efficiencies in haloes with the same mass.

The origin of the atomic and molecular gas contents of early-type galaxies – II. Misaligned gas accretion

We study the origin of the wide distribution of angles between the angular momenta of the stellar and gas components, αG, S, in early-type galaxies (ETGs). We use the GALFORM model of galaxy formation, set in the Λ cold dark matter framework, and coupled it with a Monte Carlo simulation to follow the angular momenta flips driven by matter accretion on to haloes and galaxies. We consider a gas disc to be misaligned with respect to the stellar body if αG,S > 30 deg. By assuming that the only sources of misalignments in galaxies are galaxy mergers, we place a lower limit of 2–5 per cent on the fraction of ETGs with misaligned gas/stellar components. These low fractions are inconsistent with the observed value of ≈42 ± 6 per cent in ATLAS3D. In the more general case, in which smooth gas accretion in addition to galaxy mergers can drive misalignments, our calculation predicts that ≈46 per cent of ETGs have αG, S > 30 deg. In this calculation, we find correlations between αG, S and stellar mass, cold gas fraction and star formation rate, such that ETGs with high masses, low cold gas fractions and low star formation rates are more likely to display aligned cold gas and stellar components. We confirm these trends observationally for the first time using ATLAS3D data. We argue that the high fraction of misaligned gas discs observed in ETGs is mostly due to smooth gas accretion (e.g. cooling from the hot halo of galaxies) which takes place after most of the stellar mass of the galaxy is in place and comes misaligned with respect to the stellar component. Galaxies that have accreted most of their cold gas content prior to the time where most of the stellar mass was in place show aligned components.

A PARAMETRIC STUDY OF POSSIBLE SOLUTIONS TO THE HIGH-REDSHIFT OVERPRODUCTION OF STARS IN MODELED DWARF GALAXIES

Both numerical hydrodynamic and semi-analytic cosmological models of galaxy formation struggle to match observed star formation histories of galaxies in low mass halos (M$_{\rm{H}} \lesssim 10^{11} M_\odot$), predicting more star formation at high redshift and less star formation at low redshift than observed. The fundamental problem is that galaxies' gas accretion and star formation rates are too closely coupled in the models: the accretion rate largely drives the star formation rate. Observations point to gas accretion rates that outpace star formation at high redshift, resulting in a buildup of gas and a delay in star formation until lower redshifts. We present three empirical adjustments of standard recipes in a semi-analytic model motivated by three physical scenarios that could cause this decoupling: 1) the mass-loading factors of outflows driven by stellar feedback may have a steeper dependence on halo mass at earlier times, 2) the efficiency of star formation may be lower in low mass halos at high redshift, and 3) gas may not be able to accrete efficiently onto the disk in low mass halos at high redshift. These new recipes, once tuned, better reproduce the evolution of $f_\star \equiv M_\star/M_{\rm{H}}$ as a function of halo mass as derived from abundance matching over redshifts $z=0$ to $3$, though they have different effects on cold gas fractions, star formation rates, and metallicities. Changes to gas accretion and stellar-driven winds are promising, while direct modification of the star formation timescale requires drastic measures that are not physically well-motivated.

Numerical hydrodynamic simulations based on semi-analytic galaxy merger trees: method and Milky Way-like galaxies

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

Cloud formation in colliding flows: influence of the choice of cooling function

We study the influence of the choice of cooling function on the formation of molecular clouds in high-resolution three-dimensional simulations of converging flows. We directly compare the results obtained using the simple, parametrized cooling function introduced by Koyama & Inutsuka (2002) and used by a number of converging flow studies with the results of the detailed calculation of the non-equilibrium chemistry and thermal balance of the gas. We find that a number of the cloud properties, such as the mass and volume filling fractions of cold gas, are relatively insensitive to the choice of cooling function. On the other hand, the cloud morphology and the large-scale velocity distribution of the gas do strongly depend on the cooling function. We show that the differences that we see can largely be explained by differences in the way that Lyman-alpha cooling is treated in the two complementary approaches, and that a proper non-equilibrium treatment of the ionisation and recombination of the gas is necessary in order to model the high-temperature cooling correctly. We also investigate the properties of the dense clumps formed within the cloud. In agreement with previous models, we find that the majority of these clumps are not self-gravitating, suggesting that some form of large-scale collapse of the cloud may be required in order to produce gravitationally unstable clumps and hence stars. Overall, the physical properties of the dense clumps are similar in both simulations, suggesting that they do not depend strongly on the choice of cooling function. However, we do find a systematic difference of around 10K in the mean temperatures of the clumps produced by the two models.

Origin of the antihierarchical growth of black holes

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

H i content, metallicities and spin temperatures of damped and sub-damped Lyα systems in the redshift desert (0.6 &lt; zabs &lt; 1.7)★

The HI 21cm absorption optical depth and the N(HI) derived from Lya absorption can be combined to yield the spin temperature (Ts) of DLAs. Although Ts measurements exist for samples of DLAs with z <0.6 and z >1.7, the intermediate redshift regime currently contains only 2 HI 21cm detections, leading to a `redshift desert' that spans 4 Gyrs of cosmic time. To connect the low and high z regimes, we present observations of the Lya line of six 0.6<z<1.7 HI 21cm absorbers. The dataset is complemented by both VLBA observations (to derive the absorber covering factor, f), and optical echelle spectra from which metal abundances are determined. Our dataset therefore not only offers the largest statistical study of HI 21cm absorbers to date, and bridges the redshift desert, but is also the first to use a fully f-corrected dataset to look for metallicity-based trends. In agreement with trends found in Galactic sightlines, we find that the lowest N(HI) absorbers tend to be dominated by warm gas. In the DLA regime, spin temperatures show a wider range of values than Galactic data, as may be expected in a heterogenous galactic population. However, we find that low metallicity DLAs are dominated by small cold gas fractions and only absorbers with relatively high metallicities exhibit significant fractions of cold gas. Using a compilation of HI 21cm absorbers which are selected to have f-corrected spin temperatures, we confirm an anti-correlation between metallicity and Ts at 3.4 sigma significance. Finally, one of the DLAs in our sample is a newly-discovered HI 21cm absorber (at z=0.602 towards J1431+3952), which we find to have the lowest f-corrected spin temperature yet reported in the literature: Ts=90+-23 K. The observed distribution of Ts and metallicities in DLAs and the implications for understanding the characteristics of the interstellar medium in high redshift galaxies are discussed.

The SWELLS survey - IV. Precision measurements of the stellar and dark matter distributions in a spiral lens galaxy

We construct a fully self-consistent mass model for the lens galaxy J2141 at z=0.14, and use it to improve on previous studies by modelling its gravitational lensing effect, gas rotation curve and stellar kinematics simultaneously. We adopt a very flexible axisymmetric mass model constituted by a generalized NFW dark matter halo and a stellar mass distribution obtained by deprojecting the MGE fit to the high-resolution K'-band LGSAO imaging data of the galaxy, with the (spatially constant) M/L ratio as a free parameter. We model the stellar kinematics by solving the anisotropic Jeans equations. We find that the inner logarithmic slope of the dark halo is weakly constrained (gamma = 0.82^{+0.65}_{-0.54}), and consistent with an unmodified NFW profile. We infer the galaxy to have (i) a dark matter fraction within 2.2 disk radii of 0.28^{+0.15}_{-0.10}, independent of the galaxy stellar population, implying a maximal disk for J2141; (ii) an apparently uncontracted dark matter halo, with concentration c_{-2} = 7.7_{-2.5}^{+4.2} and virial velocity v_{vir} = 242_{-39}^{+44} km/s, consistent with LCDM predictions; (iii) a slightly oblate halo (q_h = 0.75^{+0.27}_{-0.16}), consistent with predictions from baryon-affected models. Comparing the stellar mass inferred from the combined analysis (log_{10} Mstar/Msun = 11.12_{-0.09}^{+0.05}) with that inferred from SPS modelling of the galaxies colours, and accounting for a cold gas fraction of 20+/-10%, we determine a preference for a Chabrier IMF over Salpeter IMF by a Bayes factor of 5.7 (substantial evidence). We infer a value beta_{z} = 1 - sigma^2_{z}/sigma^2_{R} = 0.43_{-0.11}^{+0.08} for the orbital anisotropy parameter in the meridional plane, in agreement with most studies of local disk galaxies, and ruling out at 99% CL that the dynamics of this system can be described by a two-integral distribution function. [Abridged]

Constraints on star formation driven galaxy winds from the mass-metallicity relation at z= 0

We extend a chemical evolution model relating galaxy stellar mass and gas-phase oxygen abundance (the mass-metallicity relation) to explicitly consider the mass-dependence of galaxy gas fractions and outflows. Using empirically derived scalings of galaxy mass with halo virial velocity in conjunction with the most recent observations of z~0 total galaxy cold gas fractions and the mass-metallicity relation, we place stringent global constraints on the magnitude and scaling of the efficiency with which star forming galaxies expel metals. We demonstrate that under the assumptions that metal accretion is negligible and the stellar initial mass function does not vary, efficient outflows are required to reproduce the mass-metallicity relation; without winds, gas-to-stellar mass ratios >~ 0.3 dex higher than observed are needed. Moreover, z=0 gas fractions are low enough that while they have some effect on the magnitude of outflows required, the slope of the gas fraction--stellar mass relation does not strongly affect our conclusions on how the wind efficiencies must scale with galaxy mass. Despite systematic uncertainties in the normalization and slope of the mass-metallicity relation, we show that the metal expulsion efficiency zetaw=(Zw/Zg)etaw (where Zw is the wind metallicitiy and Zg is the interstellar medium metallicity) must be both high and scale steeply with mass. Specifically, we show that zetaw >> 1 and zetaw proportional to vvir^-3 or steeper. In contrast, momentum- or energy-driven outflow models suggest that etaw should scale as vvir^-1 or vvir^-2, respectively, implying that the Zw-Mstar relation should be shallower than the Zg-Mstar relation. [abridged]

HYDRODYNAMICS OF HIGH-REDSHIFT GALAXY COLLISIONS: FROM GAS-RICH DISKS TO DISPERSION-DOMINATED MERGERS AND COMPACT SPHEROIDS

Disk galaxies at high redshift (z~2) are characterized by high fractions of cold gas, strong turbulence, and giant star-forming clumps. Major mergers of disk galaxies at high redshift should then generally involve such turbulent clumpy disks. Merger simulations, however, model the ISM as a stable, homogeneous, and thermally pressurized medium. We present the first merger simulations with high fractions of cold, turbulent, and clumpy gas. We discuss the major new features of these models compared to models where the gas is artificially stabilized and warmed. Gas turbulence, which is already strong in high-redshift disks, is further enhanced in mergers. Some phases are dispersion-dominated, with most of the gas kinetic energy in the form of velocity dispersion and very chaotic velocity fields, unlike merger models using a thermally stabilized gas. These mergers can reach very high star formation rates, and have multi-component gas spectra consistent with SubMillimeter Galaxies. Major mergers with high fractions of cold turbulent gas are also characterized by highly dissipative gas collapse to the center of mass, with the stellar component following in a global contraction. The final galaxies are early-type with relatively small radii and high Sersic indices, like high-redshift compact spheroids. The mass fraction in a disk component that survives or re-forms after a merger is severely reduced compared to models with stabilized gas, and the formation of a massive disk component would require significant accretion of external baryons afterwards. Mergers thus appear to destroy extended disks even when the gas fraction is high, and this lends further support to smooth infall as the main formation mechanism for massive disk galaxies.

GROWING MASSIVE BLACK HOLE PAIRS IN MINOR MERGERS OF DISK GALAXIES

We perform a suite of high-resolution smoothed particle hydrodynamics simulations to investigate the orbital decay and mass evolution of massive black hole (MBH) pairs down to scales of ~30 pc during minor mergers of disk galaxies. Our simulation set includes star formation and accretion onto the MBHs, as well as feedback from both processes. We consider 1:10 merger events starting at z~3, with MBH masses in the sensitivity window of the Laser Interferometer Space Antenna, and we follow the coupling between the merger dynamics and the evolution of the MBH mass ratio until the satellite galaxy is tidally disrupted. While the more massive MBH accretes in most cases as if the galaxy were in isolation, the satellite MBH may undergo distinct episodes of enhanced accretion, owing to strong tidal torques acting on its host galaxy and to orbital circularization inside the disk of the primary galaxy. As a consequence, the initial 1:10 mass ratio of the MBHs changes by the time the satellite is disrupted. Depending on the initial fraction of cold gas in the galactic disks and the geometry of the encounter, the mass ratios of the MBH pairs at the time of satellite disruption can stay unchanged or become as large as 1:2. Remarkably, the efficiency of MBH orbital decay correlates with the final mass ratio of the pair itself: MBH pairs that increase significantly their mass ratio are also expected to inspiral more promptly down to nuclear-scale separations. These findings indicate that the mass ratios of MBH pairs in galactic nuclei do not necessarily trace the mass ratios of their merging host galaxies, but are determined by the complex interplay between gas accretion and merger dynamics.

ARE RADIO ACTIVE GALACTIC NUCLEI POWERED BY ACCRETION OR BLACK HOLE SPIN?

We compare accretion and black hole spin as potential energy sources for outbursts from AGN in brightest cluster galaxies (BCGs). Based on our adopted spin model, we find that the distribution of AGN power estimated from X-ray cavities is consistent with a broad range of both spin parameter and accretion rate. Sufficient quantities of molecular gas are available in most BCGs to power their AGN by accretion alone. However, we find no correlation between AGN power and molecular gas mass over the range of jet power considered here. For a given AGN power, the BCG's gas mass and accretion efficiency, defined as the fraction of the available cold molecular gas that is required to power the AGN, both vary by more than two orders of magnitude. Most of the molecular gas in BCGs is apparently consumed by star formation or is driven out of the nucleus by the AGN before it reaches the nuclear black hole. Bondi accretion from hot atmospheres is generally unable to fuel powerful AGN, unless their black holes are more massive than their bulge luminosities imply. We identify several powerful AGN that reside in relatively gas-poor galaxies, indicating an unusually efficient mode of accretion, or that their AGN are powered by another mechanism. If these systems are powered primarily by black hole spin, rather than by accretion, spin must also be tapped efficiently in some systems, i.e., $P_{\rm jet} > \dot Mc^2$, or their black hole masses must be substantially larger than the values implied by their bulge luminosities. We constrain the (model dependent) accretion rate at the transition from radiatively inefficient to radiatively efficient accretion flows to be a few percent of the Eddington rate, a value that is consistent with other estimates.

Forced turbulence in thermally bistable gas: a parameter study

Context: The thermal instability is one of the dynamical agents for turbulence in the diffuse interstellar medium, where both, turbulence and thermal instability interact in a highly non-linear manner. Aims: We study basic properties of turbulence in thermally bistable gas for variable simulation parameters. The resulting cold gas fractions can be applied as parameterisation in simulations on galactic scales. Methods: Turbulent flow is induced on large scales by means of compressive stochastic forcing in a periodic box. The compressible Euler equations with constant UV heating and a parameterised cooling function are solved on uniform grids. We investigate several values of the mean density of the gas and different magnitudes of the forcing. For comparison with other numerical studies, solenoidal forcing is applied as well. Results: After a transient phase, we observe that a state of statistically stationary turbulence is approached. Compressive forcing generally produces a two-phase medium, with a decreasing fraction of cold gas for increasing forcing strength. This behaviour can be explained on the basis of turbulent mixing. We also find indications for power-law tails of probability density functions of the gas density. Solenoidal forcing, on the other hand, appears to prevent the evolution into a two-phase-medium for certain parameter regions. Conclusions: The dynamics of thermally bistable turbulence shows a substantial sensitivity on the initial state and the forcing properties.