We report the first detection of a correlation between gravitational lensing by large scale structure and the thermal Sunyaev-Zeldovich (tSZ) effect. Using the mass map from the Canada France Hawaii Telescope Lensing Survey (CFHTLenS) and a newly-con
structed tSZ map from Planck, we measure a non-zero correlation between the two maps out to one degree angular separation on the sky, with an overall significance of 6 sigma. The tSZ maps are formed in a manner that removes primary cosmic microwave background fluctuations and minimizes residual contamination by galactic and extragalactic dust emission, and by CO line emission. We perform numerous tests to show that our measurement is immune to these residual contaminants. The resulting correlation function is consistent with the existence of a warm baryonic gas tracing the large scale structure with a bias b_gas. Given the shape of the lensing kernel, our signal sensitivity peaks at a redshift z~0.4, where half a degree separation on the sky corresponds to a physical scale of ~10 Mpc. The amplitude of the signal constrains the product (b_gas/1)(T_e / 0.1 keV)(n_e / 1 m^-3)=2.01pm 0.31pm 0.21, at redshift zero. Our study suggests that a substantial fraction of the missing baryons in the universe may reside in a low density warm plasma that traces dark matter.
We consider the effects of different star formation criteria on galactic scales, in high-resolution simulations with explicitly resolved GMCs and stellar feedback. We compare: (1) a self-gravity criterion (based on the local virial parameter and the
assumption that self-gravitating gas collapses to high density in a free-fall time), (2) a fixed density threshold, (3) a molecular-gas law, (4) a temperature threshold, (5) a Jeans-instability requirement, (6) a criteria that cooling times be shorter than dynamical times, and (7) a convergent-flow criterion. We consider these both MW-like and high-density (starburst) galaxies. With feedback present, all models produce identical integrated star formation rates (SFRs), in agreement with the Kennicutt relation. Without feedback all produce orders-of-magnitude excessive SFRs. This is totally dependent on feedback and independent of the SF law. However, the spatial and density distribution of SF depend strongly on the SF criteria. Because cooling rates are generally fast and gas is turbulent, criteria (4)-(7) are weak and spread SF uniformly over the disk (above densities n~0.01-0.1 cm^-3). A molecular criterion (3) localizes to higher densities, but still a wide range; for Z Z_solar, it is similar to a density threshold at n~1 cm^-3 (well below mean densities in the MW center or starbursts). Fixed density thresholds (2) can always select the highest densities, but must be adjusted for simulation resolution and galaxy properties; the same threshold that works in a MW-like simulation will select nearly all gas in a starburst. Binding criteria (1) tend to adaptively select the largest over-densities, independent of galaxy model or resolution, and automatically predict clustered SF. We argue that this SF model is most physically-motivated and presents significant numerical advantages in large-dynamic range simulations.
Rapid accretion of cold gas plays a crucial role in getting gas into galaxies. It has been suggested that this accretion proceeds along narrow streams that might also directly drive the turbulence in galactic gas, dynamical disturbances, and bulge fo
rmation. In cosmological simulations, however, it is impossible to isolate and hence disentangle the effect of accretion from internal instabilities and mergers. Moreover, in most cosmological simulations, the phase structure and turbulence in the ISM arising from stellar feedback are treated in a sub-grid manner, so that feedback cannot generate ISM turbulence. In this paper we therefore test the effects of cold streams in extremely high-resolution simulations of otherwise isolated galaxy disks using detailed models for star formation and feedback; we then include or exclude mock cold flows falling onto the galaxies with accretion rates, velocities and geometry set to maximize their effect on the disk. We find: (1) Turbulent velocity dispersions in gas disks are identical with or without the cold flow; the energy injected by the flow is dissipated where it meets the disk. (2) In runs without stellar feedback, the presence of a cold flow has essentially no effect on runaway local collapse, resulting in star formation rates (SFRs) that are far too large. (3) Disks in runs with feedback and cold flows have higher SFRs, but only insofar as they have more gas. (4) Because flows are extended relative to the disk, they do not trigger strong resonant responses and so induce weak morphological perturbation (bulge formation via instabilities is not accelerated). (5) However, flows can thicken the disk by direct contribution of out-of-plane streams. We conclude that while inflows are critical over cosmological timescales to determine the supply and angular momentum of gas disks, they have weak instantaneous dynamical effects on galaxies.
We study galaxy super-winds driven in major mergers, using pc-resolution simulations with detailed models for stellar feedback that can self-consistently follow the formation/destruction of GMCs and generation of winds. The models include molecular c
ooling, star formation at high densities in GMCs, and gas recycling and feedback from SNe (I&II), stellar winds, and radiation pressure. We study mergers of systems from SMC-like dwarfs and Milky Way analogues to z~2 starburst disks. Multi-phase super-winds are generated in all passages, with outflow rates up to ~1000 M_sun/yr. However, the wind mass-loading efficiency (outflow rate divided by SFR) is similar to that in isolated galaxy counterparts of each merger: it depends more on global galaxy properties (mass, size, escape velocity) than on the dynamical state of the merger. Winds tend to be bi- or uni-polar, but multiple events build up complex morphologies with overlapping, differently-oriented bubbles/shells at a range of radii. The winds have complex velocity and phase structure, with material at a range of speeds up to ~1000 km/s, and a mix of molecular, ionized, and hot gas that depends on galaxy properties and different feedback mechanisms. These simulations resolve a problem in some sub-grid models, where simple wind prescriptions can dramatically suppress merger-induced starbursts. But despite large mass-loading factors (>~10) in the winds, the peak SFRs are comparable to those in no wind simulations. Wind acceleration does not act equally, so cold dense gas can still lose angular momentum and form stars, while blowing out gas that would not have participated in the starburst in the first place. Considerable wind material is not unbound, and falls back on the disk at later times post-merger, leading to higher post-starburst SFRs in the presence of stellar feedback. This may require AGN feedback to explain galaxy quenching.
We show that the mass fraction of GMC gas (n>100 cm^-3) in dense (n>>10^4 cm^-3) star-forming clumps, observable in dense molecular tracers (L_HCN/L_CO(1-0)), is a sensitive probe of the strength and mechanism(s) of stellar feedback. Using high-resol
ution galaxy-scale simulations with pc-scale resolution and explicit models for feedback from radiation pressure, photoionization heating, stellar winds, and supernovae (SNe), we make predictions for the dense molecular gas tracers as a function of GMC and galaxy properties and the efficiency of stellar feedback. In models with weak/no feedback, much of the mass in GMCs collapses into dense sub-units, predicting L_HCN/L_CO(1-0) ratios order-of-magnitude larger than observed. By contrast, models with feedback properties taken directly from stellar evolution calculations predict dense gas tracers in good agreement with observations. Changing the strength or timing of SNe tends to move systems along, rather than off, the L_HCN-L_CO relation (because SNe heat lower-density material, not the high-density gas). Changing the strength of radiation pressure (which acts efficiently in the highest density gas), however, has a much stronger effect on L_HCN than on L_CO. We predict that the fraction of dense gas (L_HCN/L_CO(1-0)) increases with increasing GMC surface density; this drives a trend in L_HCN/L_CO(1-0) with SFR and luminosity which has tentatively been observed. Our results make specific predictions for enhancements in the dense gas tracers in unusually dense environments such as ULIRGs and galactic nuclei (including the galactic center).
We use numerical simulations of isolated galaxies to study the effects of stellar feedback on the formation and evolution of giant star-forming gas clumps in high-redshift, gas-rich galaxies. Such galactic disks are unstable to the formation of bound
gas-rich clumps whose properties initially depend only on global disk properties, not the microphysics of feedback. In simulations without stellar feedback, clumps turn an order-unity fraction of their mass into stars and sink to the center, forming a large bulge and kicking most of the stars out into a much more extended stellar envelope. By contrast, strong radiative stellar feedback disrupts even the most massive clumps after they turn ~10-20% of their mass into stars, in a timescale of ~10-100 Myr, ejecting some material into a super-wind and recycling the rest of the gas into the diffuse ISM. This suppresses the bulge formation rate by direct clump coalescence by a factor of several. However, the galactic disks do undergo significant internal evolution in the absence of mergers: clumps form and disrupt continuously and torque gas to the galactic center. The resulting evolution is qualitatively similar to bar/spiral evolution in simulations with a more homogeneous ISM.
Feedback from massive stars is believed to play a critical role in shaping the galaxy mass function, the structure of the interstellar medium (ISM), and the low efficiency of star formation, but the exact form of the feedback is uncertain. In this pa
per, the first in a series, we present and test a novel numerical implementation of stellar feedback resulting from momentum imparted to the ISM by radiation, supernovae, and stellar winds. We employ a realistic cooling function, and find that a large fraction of the gas cools to <100K, so that the ISM becomes highly inhomogeneous. Despite this, our simulated galaxies reach an approximate steady state, in which gas gravitationally collapses to form giant molecular clouds (GMCs), dense clumps, and stars; subsequently, stellar feedback disperses the GMCs, repopulating the diffuse ISM. This collapse and dispersal cycle is seen in models of SMC-like dwarfs, the Milky-Way, and z~2 clumpy disk analogues. The simulated global star formation efficiencies are consistent with the observed Kennicutt-Schmidt relation. Moreover, the star formation rates are nearly independent of the numerically imposed high-density star formation efficiency, density threshold, and density scaling. This is a consequence of the fact that, in our simulations, star formation is regulated by stellar feedback limiting the amount of very dense gas available for forming stars. In contrast, in simulations without stellar feedback, i.e. under the action of only gravity and gravitationally-induced turbulence, the ISM experiences runaway collapse to very high densities. In these simulations without feedback, the global star formation rates exceed observed galactic star formation rates by 1-2 orders of magnitude, demonstrating that stellar feedback is crucial to the regulation of star formation in galaxies.
We compile observations of the surface mass density profiles of dense stellar systems, including globular clusters in the Milky Way and nearby galaxies, massive star clusters in nearby starbursts, nuclear star clusters in dwarf spheroidals and late-t
ype disks, ultra-compact dwarfs, and galaxy spheroids spanning the range from low-mass cusp bulges and ellipticals to massive core ellipticals. We show that in all cases the maximum stellar surface density attained in the central regions of these systems is similar, Sigma_max ~ 10^11 M_sun/kpc^2 (~20 g/cm^2), despite the fact that the systems span 7 orders of magnitude in total stellar mass M_star, 5 in effective radius R_e, and have a wide range in effective surface density M_star/R_e^2. The surface density limit is reached on a wide variety of physical scales in different systems and is thus not a limit on three-dimensional stellar density. Given the very different formation mechanisms involved in these different classes of objects, we argue that a single piece of physics likely determines Sigma_max. The radiation fields and winds produced by massive stars can have a significant influence on the formation of both star clusters and galaxies, while neither supernovae nor black hole accretion are important in star cluster formation. We thus conclude that feedback from massive stars likely accounts for the observed Sigma_max, plausibly because star formation reaches an Eddington-like flux that regulates the growth of these diverse systems. This suggests that current models of galaxy formation, which focus on feedback from supernovae and active galactic nuclei, are missing a crucial ingredient.
Supermassive black holes (BHs) obey tight scaling relations between their mass and their host galaxy properties such as total stellar mass, velocity dispersion, and potential well depth. This has led to the development of self-regulated models for BH
growth, in which feedback from the central BH halts its own growth upon reaching a critical threshold. However, models have also been proposed in which feedback plays no role: so long as a fixed fraction of the host gas supply is accreted, relations like those observed can be reproduced. Here, we argue that the scatter in the observed BH-host correlations, and its run with scale, presents a demanding constraint on any model for these correlations, and that it favors self-regulated models of BH growth. We show that the scatter in the stellar mass fraction within a radius R in observed ellipticals and spheroids increases strongly at small R. At fixed total stellar mass (or host velocity dispersion), on very small scales near the BH radius of influence, there is an order-of-magnitude scatter in the amount of gas that must have entered and formed stars. In short, the BH appears to know more about the global host galaxy potential on large scales than the stars and gas supply on small scales. This is predicted in self-regulated models; however, models where there is no feedback would generically predict order-of-magnitude scatter in the BH-host correlations. Likewise, models in which the BH feedback in the bright mode does not regulate the growth of the BH itself, but sets the stellar mass of the galaxy by inducing star formation or blowing out a mass in gas much larger than the galaxy stellar mass, are difficult to reconcile with the scatter on small scales.
