Source author record

Norman Murray

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

ResearcherUnclaimed source record

astro-ph.GA astro-ph.CO astro-ph.SR astro-ph.HE astro-ph.EP astro-ph astro-ph.IM gr-qc physics.ao-ph

Catalog footprint

What is connected

58works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

COMAP Early Science: IV. Power Spectrum Methodology and Results

We present the power spectrum methodology used for the first-season COMAP analysis, and assess the quality of the current data set. The main results are derived through the Feed-feed Pseudo-Cross-Spectrum (FPXS) method, which is a robust estimator with respect to both noise modeling errors and experimental systematics. We use effective transfer functions to take into account the effects of instrumental beam smoothing and various filter operations applied during the low-level data processing. The power spectra estimated in this way have allowed us to identify a systematic error associated with one of our two scanning strategies, believed to be due to residual ground or atmospheric contamination. We omit these data from our analysis and no longer use this scanning technique for observations. We present the power spectra from our first season of observing and demonstrate that the uncertainties are integrating as expected for uncorrelated noise, with any residual systematics suppressed to a level below the noise. Using the FPXS method, and combining data on scales $k=0.051-0.62 \,\mathrm{Mpc}^{-1}$ we estimate $P_\mathrm{CO}(k) = -2.7 \pm 1.7 \times 10^4μ\textrm{K}^2\mathrm{Mpc}^3$, the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum in the literature.

preprint2022arXiv

COMAP Early Science: V. Constraints and Forecasts at $z \sim 3$

We present the current state of models for the $z\sim3$ carbon monoxide (CO) line-intensity signal targeted by the CO Mapping Array Project (COMAP) Pathfinder in the context of its early science results. Our fiducial model, relating dark matter halo properties to CO luminosities, informs parameter priors with empirical models of the galaxy-halo connection and previous CO(1-0) observations. The Pathfinder early science data spanning wavenumbers $k=0.051$-$0.62\,$Mpc$^{-1}$ represent the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum. Our 95% upper limit on the redshift-space clustering amplitude $A_{\rm clust}\lesssim70\,μ$K$^2$ greatly improves on the indirect upper limit of $420\,μ$K$^2$ reported from the CO Power Spectrum Survey (COPSS) measurement at $k\sim1\,$Mpc$^{-1}$. The COMAP limit excludes a subset of models from previous literature, and constrains interpretation of the COPSS results, demonstrating the complementary nature of COMAP and interferometric CO surveys. Using line bias expectations from our priors, we also constrain the squared mean line intensity-bias product, $\langle{Tb}\rangle^2\lesssim50\,μ$K$^2$, and the cosmic molecular gas density, $ρ_\text{H2}<2.5\times10^8\,M_\odot\,$Mpc$^{-3}$ (95% upper limits). Based on early instrument performance and our current CO signal estimates, we forecast that the five-year Pathfinder campaign will detect the CO power spectrum with overall signal-to-noise of 9-17. Between then and now, we also expect to detect the CO-galaxy cross-spectrum using overlapping galaxy survey data, enabling enhanced inferences of cosmic star-formation and galaxy-evolution history.

preprint2022arXiv

Gas infall and radial transport in cosmological simulations of Milky Way-mass disks

Observations indicate that a continuous supply of gas is needed to maintain observed star formation rates in large, disky galaxies. To fuel star formation, gas must reach the inner regions of such galaxies. Despite its crucial importance for galaxy evolution, how and where gas joins galaxies is poorly constrained observationally and is rarely explored in fully cosmological simulations. To investigate gas accretion in the vicinity of galaxies, we analyze the FIRE-2 cosmological zoom-in simulations for 4 Milky Way mass galaxies (M_halo ~ 10E12 solar masses), focusing on simulations with cosmic ray physics. We find that at z~0, gas approaches the disk with angular momentum similar to the gaseous disk edge and low radial velocities, piling-up near the edge and settling into full rotational support. Accreting gas moves predominantly parallel to the disk with small but nonzero vertical velocity components, and joins the disk largely in the outskirts as opposed to "raining" down onto the disk. Once in the disk, gas trajectories are complex, being dominated by spiral arm induced oscillations and feedback. However, time and azimuthal averages show clear but slow net radial infall with transport speeds of 1-3 km/s and net mass fluxes through the disk on the order of one solar mass per year, comparable to the star formation rates of the galaxies and decreasing towards galactic center as gas is sunk into star formation. These rates are slightly higher in simulations without cosmic rays (1-7 km/s, ~4-5 solar masses per year). We find overall consistency of our results with observational constraints and discuss prospects of future observations of gas flows in and around galaxies.

preprint2022arXiv

Kinematics and Feedback in H II regions in the Dwarf Starburst Galaxy IC 10

We present a survey of the central region of the nearest starburst galaxy, IC 10, using the W. M. Keck Observatory Keck Cosmic Web Imager (KCWI) at high spectral and spatial resolution. We map the central starburst of IC 10 to sample the kinematic and ionization properties of the individual star-forming regions. Using the low spectral resolution mode of KCWI we map the oxygen abundance and with the high spectral resolution mode we identify 46 individual H II regions. These H II regions have an average radius of 4.0 pc, star formation rate $\sim1.3\times10^{-4}$ M$_\odot$ yr$^{-1}$, and velocity dispersion $\sim$16 km s$^{-1}$. None of the H II regions appear to be virialized ($\rm α_{vir}>>1$), and, on average, they show evidence of ongoing expansion. IC 10's H II regions are offset from the star forming region size-luminosity scaling relationships, as well as Larson's Law that relates size and velocity dispersion. We investigate the balance of inward and outward pressure, $\rm P_{in}$ and $\rm P_{out}$, finding $\rm P_{out}>P_{in}$ in 89% of H II regions, indicating feedback driven expansion even in these low mass H II regions. We find warm gas pressure ($\rm P_{gas}$) provides the dominant contribution to the outward pressure ($\rm P_{out}$). This counteracts the inward pressure which is dominated by turbulence in the surrounding gas rather than self-gravity. Five H II regions show evidence of outflows which are most likely supported by either stellar winds (2 regions) or champagne flows (3 regions). These observations provide new insights into the state of the star-forming regions in IC 10 and negative feedback from low mass clusters.

preprint2021arXiv

COMAP Early Science: I. Overview

The CO Mapping Array Project (COMAP) aims to use line intensity mapping of carbon monoxide (CO) to trace the distribution and global properties of galaxies over cosmic time, back to the Epoch of Reionization (EoR). To validate the technologies and techniques needed for this goal, a Pathfinder instrument has been constructed and fielded. Sensitive to CO(1-0) emission from $z=2.4$-$3.4$ and a fainter contribution from CO(2-1) at $z=6$-8, the Pathfinder is surveying $12$ deg$^2$ in a 5-year observing campaign to detect the CO signal from $z\sim3$. Using data from the first 13 months of observing, we estimate $P_\mathrm{CO}(k) = -2.7 \pm 1.7 \times 10^4μ\mathrm{K}^2 \mathrm{Mpc}^3$ on scales $k=0.051-0.62 \mathrm{Mpc}^{-1}$ - the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum. Based on these observations alone, we obtain a constraint on the amplitude of the clustering component (the squared mean CO line temperature-bias product) of $\langle Tb\rangle^2<49$ $μ$K$^2$ - nearly an order-of-magnitude improvement on the previous best measurement. These constraints allow us to rule out two models from the literature. We forecast a detection of the power spectrum after 5 years with signal-to-noise ratio (S/N) 9-17. Cross-correlation with an overlapping galaxy survey will yield a detection of the CO-galaxy power spectrum with S/N of 19. We are also conducting a 30 GHz survey of the Galactic plane and present a preliminary map. Looking to the future of COMAP, we examine the prospects for future phases of the experiment to detect and characterize the CO signal from the EoR.

preprint2021arXiv

COMAP Early Science: III. CO Data Processing

We describe the first season COMAP analysis pipeline that converts raw detector readouts to calibrated sky maps. This pipeline implements four main steps: gain calibration, filtering, data selection, and map-making. Absolute gain calibration relies on a combination of instrumental and astrophysical sources, while relative gain calibration exploits real-time total-power variations. High efficiency filtering is achieved through spectroscopic common-mode rejection within and across receivers, resulting in nearly uncorrelated white noise within single-frequency channels. Consequently, near-optimal but biased maps are produced by binning the filtered time stream into pixelized maps; the corresponding signal bias transfer function is estimated through simulations. Data selection is performed automatically through a series of goodness-of-fit statistics, including $χ^2$ and multi-scale correlation tests. Applying this pipeline to the first-season COMAP data, we produce a dataset with very low levels of correlated noise. We find that one of our two scanning strategies (the Lissajous type) is sensitive to residual instrumental systematics. As a result, we no longer use this type of scan and exclude data taken this way from our Season 1 power spectrum estimates. We perform a careful analysis of our data processing and observing efficiencies and take account of planned improvements to estimate our future performance. Power spectrum results derived from the first-season COMAP maps are presented and discussed in companion papers.

preprint2021arXiv

COMAP Early Science: VII. Prospects for CO Intensity Mapping at Reionization

We introduce COMAP-EoR, the next generation of the Carbon Monoxide Mapping Array Project aimed at extending CO intensity mapping to the Epoch of Reionization. COMAP-EoR supplements the existing 30 GHz COMAP Pathfinder with two additional 30 GHz instruments and a new 16 GHz receiver. This combination of frequencies will be able to simultaneously map CO(1--0) and CO(2--1) at reionization redshifts ($z\sim5-8$) in addition to providing a significant boost to the $z\sim3$ sensitivity of the Pathfinder. We examine a set of existing models of the EoR CO signal, and find power spectra spanning several orders of magnitude, highlighting our extreme ignorance about this period of cosmic history and the value of the COMAP-EoR measurement. We carry out the most detailed forecast to date of an intensity mapping cross-correlation, and find that five out of the six models we consider yield signal to noise ratios (S/N) $\gtrsim20$ for COMAP-EoR, with the brightest reaching a S/N above 400. We show that, for these models, COMAP-EoR can make a detailed measurement of the cosmic molecular gas history from $z\sim2-8$, as well as probe the population of faint, star-forming galaxies predicted by these models to be undetectable by traditional surveys. We show that, for the single model that does not predict numerous faint emitters, a COMAP-EoR-type measurement is required to rule out their existence. We briefly explore prospects for a third-generation Expanded Reionization Array (COMAP-ERA) capable of detecting the faintest models and characterizing the brightest signals in extreme detail.

preprint2021arXiv

The Mass of the Milky Way from the H3 Survey

The mass of the Milky Way is a critical quantity which, despite decades of research, remains uncertain within a factor of two. Until recently, most studies have used dynamical tracers in the inner regions of the halo, relying on extrapolations to estimate the mass of the Milky Way. In this paper, we extend the hierarchical Bayesian model applied in Eadie & Jurić (2019) to study the mass distribution of the Milky Way halo; the new model allows for the use of all available 6D phase-space measurements. We use kinematic data of halo stars out to $142~{\rm kpc}$, obtained from the H3 Survey and $\textit{Gaia}$ EDR3, to infer the mass of the Galaxy. Inference is carried out with the No-U-Turn sampler, a fast and scalable extension of Hamiltonian Monte Carlo. We report a median mass enclosed within $100~{\rm kpc}$ of $\rm M(<100 \; kpc) = 0.69_{-0.04}^{+0.05} \times 10^{12} \; M_\odot$ (68% Bayesian credible interval), or a virial mass of $\rm M_{200} = M(<216.2_{-7.5}^{+7.5} \; kpc) = 1.08_{-0.11}^{+0.12} \times 10^{12} \; M_\odot$, in good agreement with other recent estimates. We analyze our results using posterior predictive checks and find limitations in the model's ability to describe the data. In particular, we find sensitivity with respect to substructure in the halo, which limits the precision of our mass estimates to $\sim 15\%$.

preprint2020arXiv

Properties of the Circumgalactic Medium in Cosmic Ray-Dominated Galaxy Halos

We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass halos hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, collisional and streaming losses, with constant parallel diffusivity $κ\sim3\times10^{29}\,\mathrm{cm^2\ s^{-1}}$ chosen to match $γ$-ray observations. With this, CRs become more important at larger halo masses and lower redshifts, and dominate the pressure in the CGM in MW-mass halos at $z\lesssim 1-2$. The gas in these "CR-dominated" halos differs significantly from runs without CRs: the gas is primarily cool (a few $\sim10^{4}\,$K), and the cool phase is volume-filling and has a thermal pressure below that needed for virial or local thermal pressure balance. Ionization of the "low" and "mid" ions in this diffuse cool gas is dominated by photo-ionization, with O VI columns $\gtrsim 10^{14.5}\,\mathrm{cm^{-2}}$ at distances $\gtrsim 150\,\mathrm{kpc}$. CR and thermal gas pressure are locally anti-correlated, maintaining total pressure balance, and the CGM gas density profile is determined by the balance of CR pressure gradients and gravity. Neglecting CRs, the same halos are primarily warm/hot ($T\gtrsim 10^{5}\,$K) with thermal pressure balancing gravity, collisional ionization dominates, O VI columns are lower and Ne VIII higher, and the cool phase is confined to dense filaments in local thermal pressure equilibrium with the hot phase.

preprint2020arXiv

Radiative Stellar Feedback in Galaxy Formation: Methods and Physics

Radiative feedback (RFB) from stars plays a key role in galaxies, but remains poorly-understood. We explore this using high-resolution, multi-frequency radiation-hydrodynamics (RHD) simulations from the Feedback In Realistic Environments (FIRE) project. We study ultra-faint dwarf through Milky Way mass scales, including H+He photo-ionization; photo-electric, Lyman Werner, Compton, and dust heating; and single+multiple scattering radiation pressure (RP). We compare distinct numerical algorithms: ray-based LEBRON (exact when optically-thin) and moments-based M1 (exact when optically-thick). The most important RFB channels on galaxy scales are photo-ionization heating and single-scattering RP: in all galaxies, most ionizing/far-UV luminosity (~1/2 of lifetime-integrated bolometric) is absorbed. In dwarfs, the most important effect is photo-ionization heating from the UV background suppressing accretion. In MW-mass galaxies, meta-galactic backgrounds have negligible effects; but local photo-ionization and single-scattering RP contribute to regulating the galactic star formation efficiency and lowering central densities. Without some RFB (or other 'rapid' FB), resolved GMCs convert too-efficiently into stars, making galaxies dominated by hyper-dense, bound star clusters. This makes star formation more violent and 'bursty' when SNe explode in these hyper-clustered objects: thus, including RFB 'smoothes' SFHs. These conclusions are robust to RHD methods, but M1 produces somewhat stronger effects. Like in previous FIRE simulations, IR multiple-scattering is rare (negligible in dwarfs, ~10% of RP in massive galaxies): absorption occurs primarily in 'normal' GMCs with A_v~1.

preprint2020arXiv

Self-consistent proto-globular cluster formation in cosmological simulations of high-redshift galaxies

We report the formation of bound star clusters in a sample of high-resolution cosmological zoom-in simulations of z>5 galaxies from the FIRE project. We find that bound clusters preferentially form in high-pressure clouds with gas surface densities over 10^4 Msun pc^-2, where the cloud-scale star formation efficiency is near unity and young stars born in these regions are gravitationally bound at birth. These high-pressure clouds are compressed by feedback-driven winds and/or collisions of smaller clouds/gas streams in highly gas-rich, turbulent environments. The newly formed clusters follow a power-law mass function of dN/dM~M^-2. The cluster formation efficiency is similar across galaxies with stellar masses of ~10^7-10^10 Msun at z>5. The age spread of cluster stars is typically a few Myrs and increases with cluster mass. The metallicity dispersion of cluster members is ~0.08 dex in [Z/H] and does not depend on cluster mass significantly. Our findings support the scenario that present-day old globular clusters (GCs) were formed during relatively normal star formation in high-redshift galaxies. Simulations with a stricter/looser star formation model form a factor of a few more/fewer bound clusters per stellar mass formed, while the shape of the mass function is unchanged. Simulations with a lower local star formation efficiency form more stars in bound clusters. The simulated clusters are larger than observed GCs due to finite resolution. Our simulations are among the first cosmological simulations that form bound clusters self-consistently in a wide range of high-redshift galaxies.

preprint2019arXiv

Swirls of FIRE: Spatially Resolved Gas Velocity Dispersions and Star Formation Rates in FIRE-2 Disk Environments

We study the spatially resolved (sub-kpc) gas velocity dispersion ($σ$)--star formation rate (SFR) relation in the FIRE-2 (Feedback in Realistic Environments) cosmological simulations. We specifically focus on Milky Way mass disk galaxies at late times. In agreement with observations, we find a relatively flat relationship, with $σ\approx 15-30$ km/s in neutral gas across 3 dex in SFRs. We show that higher dense gas fractions (ratios of dense gas to neutral gas) and SFRs are correlated at constant $σ$. Similarly, lower gas fractions (ratios of gas to stellar mass) are correlated with higher $σ$ at constant SFR. The limits of the $σ$-$Σ_{\rm SFR}$ relation correspond to the onset of strong outflows. We see evidence of "on-off" cycles of star formation in the simulations, corresponding to feedback injection timescales of 10-100 Myr, where SFRs oscillate about equilibrium SFR predictions. Finally, SFRs and velocity dispersions in the simulations agree well with feedback-regulated and marginally stable gas disk (Toomre's $Q =1$) model predictions, and the data effectively rule out models assuming that gas turns into stars at (low) constant efficiency (i.e., ${\rm 1\%}$ per free-fall time). And although the simulation data do not entirely exclude gas accretion/gravitationally powered turbulence as a driver of $σ$, it appears to be strongly subdominant to stellar feedback in the simulated galaxy disks.

preprint2016arXiv

Binary Stars Can Provide the "Missing Photons" Needed for Reionization

Empirical constraints on reionization require galactic ionizing photon escape fractions fesc>20%, but recent high-resolution radiation-hydrodynamic calculations have consistently found much lower values ~1-5%. While these models include strong stellar feedback and additional processes such as runaway stars, they almost exclusively consider stellar evolution models based on single (isolated) stars, despite the fact that most massive stars are in binaries. We re-visit these calculations, combining radiative transfer and high-resolution cosmological simulations from the Feedback in Realistic Environments (FIRE) project. For the first time, we use a stellar evolution model that includes a physically and observationally motivated treatment of binaries (the BPASS model). Binary mass transfer and mergers enhance the population of massive stars at late times (>3 Myr) after star formation, which in turn strongly enhances the late-time ionizing photon production (especially at low metallicities). These photons are produced after feedback from massive stars has carved escape channels in the ISM, and so efficiently leak out of galaxies. As a result, the time-averaged "effective" escape fraction (ratio of escaped ionizing photons to observed 1500 A photons) increases by factors 4-10, sufficient to explain reionization. While important uncertainties remain, we conclude that binary evolution may be critical for understanding the ionization of the Universe.

preprint2016arXiv

Detection of lensing substructure using ALMA observations of the dusty galaxy SDP.81

We study the abundance of substructure in the matter density near galaxies using ALMA Science Verification observations of the strong lensing system SDP.81. We present a method to measure the abundance of subhalos around galaxies using interferometric observations of gravitational lenses. Using simulated ALMA observations, we explore the effects of various systematics, including antenna phase errors and source priors, and show how such errors may be measured or marginalized. We apply our formalism to ALMA observations of SDP.81. We find evidence for the presence of a $M=10^{8.96\pm 0.12} M_{\odot}$ subhalo near one of the images, with a significance of $6.9σ$ in a joint fit to data from bands 6 and 7; the effect of the subhalo is also detected in both bands individually. We also derive constraints on the abundance of dark matter subhalos down to $M\sim 2\times 10^7 M_{\odot}$, pushing down to the mass regime of the smallest detected satellites in the Local Group, where there are significant discrepancies between the observed population of luminous galaxies and predicted dark matter subhalos. We find hints of additional substructure, warranting further study using the full SDP.81 dataset (including, for example, the spectroscopic imaging of the lensed carbon monoxide emission). We compare the results of this search to the predictions of $Λ$CDM halos, and find that given current uncertainties in the host halo properties of SDP.81, our measurements of substructure are consistent with theoretical expectations. Observations of larger samples of gravitational lenses with ALMA should be able to improve the constraints on the abundance of galactic substructure.

preprint2016arXiv

Giant clumps in the FIRE simulations: a case study of a massive high-redshift galaxy

The morphology of massive star-forming galaxies at high redshift is often dominated by giant clumps of mass ~10^8-10^9 Msun and size ~100-1000 pc. Previous studies have proposed that giant clumps might have an important role in the evolution of their host galaxy, particularly in building the central bulge. However, this depends on whether clumps live long enough to migrate from their original location in the disc or whether they get disrupted by their own stellar feedback before reaching the centre of the galaxy. We use cosmological hydrodynamical simulations from the FIRE (Feedback in Realistic Environments) project that implement explicit treatments of stellar feedback and ISM physics to study the properties of these clumps. We follow the evolution of giant clumps in a massive (stellar mass ~10^10.8 Msun at z=1), discy, gas-rich galaxy from redshift z>2 to z=1. Even though the clumpy phase of this galaxy lasts over a gigayear, individual gas clumps are short-lived, with mean lifetime of massive clumps of ~20 Myr. During that time, they turn between 0.1% and 20% of their gas into stars before being disrupted, similar to local GMCs. Clumps with M>10^7 Msun account for ~20% of the total star formation in the galaxy during the clumpy phase, producing ~10^10 Msun of stars. We do not find evidence for net inward migration of clumps within the galaxy. The number of giant clumps and their mass decrease at lower redshifts, following the decrease in the overall gas fraction and star-formation rate.

preprint2016arXiv

Observational Evidence of Dynamic Star Formation Rate in Milky Way Giant Molecular Clouds

Star formation on galactic scales is known to be a slow process, but whether it is slow on smaller scales is uncertain. We cross-correlate 5469 giant molecular clouds (GMCs) from a new all-sky catalog with 256 star forming complexes (SFCs) to build a sample of 191 SFC-GMC complexes---collections of multiple clouds each matched to 191 SFCs. The total mass in stars harbored by these clouds is inferred from WMAP free-free fluxes. We measure the GMC mass, the virial parameter, the star formation efficiency $ε$ and the star formation rate per free-fall time $ε_{\rm ff}$. Both $ε$ and $ε_{\rm ff}$ range over 3--4 orders of magnitude. We find that 68.3% of the clouds fall within $σ_{\logε}=0.79\pm0.22\,{\rm dex}$ and $σ_{\logε_{\rm ff}}=0.91\pm0.22\,{\rm dex}$ about the median. Compared to these observed scatters, a simple model with a time independent $ε_{\rm ff}$ that depends on the host GMC properties predicts $σ_{\logε_{\rm ff}}=$0.12-0.24. Allowing for a time-variable $ε_{\rm ff}$, we can recover the large dispersion in the rate of star formation. This strongly suggests that star formation in the Milky Way is a dynamic process on GMC scales. We also show that the surface star formation rate profile of the Milky Way correlates well with the molecular gas surface density profile.

preprint2016arXiv

Stellar and Quasar Feedback in Concert: Effects on AGN Accretion, Obscuration, and Outflows

We study the interaction of feedback from active galactic nuclei (AGN) and a multi-phase interstellar medium (ISM), in simulations including explicit stellar feedback, multi-phase cooling, accretion-disk winds, and Compton heating. We examine radii ~0.1-100 pc around a black hole (BH), where the accretion rate onto the BH is determined and where AGN-powered winds and radiation couple to the ISM. We conclude: (1) The BH accretion rate is determined by exchange of angular momentum between gas and stars in gravitational instabilities. This produces accretion rates ~0.03-1 Msun/yr, sufficient to power luminous AGN. (2) The gas disk in the galactic nucleus undergoes an initial burst of star formation followed by several Myrs where stellar feedback suppresses the star formation rate (SFR). (3) AGN winds injected at small radii with momentum fluxes ~L/c couple efficiently to the ISM and have dramatic effects on ISM properties within ~100 pc. AGN winds suppress the nuclear SFR by factors ~10-30 and BH accretion rate by factors ~3-30. They increase the outflow rate from the nucleus by factors ~10, consistent with observational evidence for galaxy-scale AGN-driven outflows. (4) With AGN feedback, the predicted column density distribution to the BH is consistent with observations. Absent AGN feedback, the BH is isotropically obscured and there are not enough optically-thin sightlines to explain Type-I AGN. A 'torus-like' geometry arises self-consistently as AGN feedback evacuates gas in polar regions.

preprint2015arXiv

Asynchronous rotation of Earth-mass planets in the habitable zone of lower-mass stars

Planets in the habitable zone of lower-mass stars are often assumed to be in a state of tidally synchronized rotation, which would considerably affect their putative habitability. Although thermal tides cause Venus to rotate retrogradely, simple scaling arguments tend to attribute this peculiarity to the massive Venusian atmosphere. Using a global climate model, we show that even a relatively thin atmosphere can drive terrestrial planets' rotation away from synchronicity. We derive a more realistic atmospheric tide model that predicts four asynchronous equilibrium spin states, two being stable, when the amplitude of the thermal tide exceeds a threshold that is met for habitable Earth-like planets with a 1-bar atmosphere around stars more massive than 0.5-0.7Msun. Thus, many recently discovered terrestrial planets could exhibit asynchronous spin-orbit rotation, even with a thin atmosphere.

preprint2015arXiv

Dynamics of Dusty Radiation Pressure Driven Shells and Clouds: Fast Outflows from Galaxies, Star Clusters, Massive Stars, and AGN

It is typically assumed that radiation pressure driven winds are accelerated to an asymptotic velocity of V ~ v_esc, where v_esc is the escape velocity from the central source. We note that this is not the case for dusty shells and clouds. Instead, if the shell or cloud is initially optically-thick to the UV emission from the source of luminosity L, then there is a significant boost in V that reflects the integral of the momentum absorbed as it is accelerated. For shells reaching a generalized Eddington limit, we show that V ~ (4R_UV L/M_sh c)^1/2, in both point-mass and isothermal-sphere potentials, where R_UV is the radius where the shell becomes optically-thin to UV photons, and M_sh is the mass of the shell. The asymptotic velocity significantly exceeds v_esc for typical parameters, and can explain the ~1000-2000km/s outflows observed from rapidly star-forming galaxies and active galactic nuclei if the surrounding halo has low gas density. Similarly fast outflows from massive stars can be accelerated on few - 10^3 yr timescales. These results carry over to clouds that subtend only a small fraction of the solid angle from the source of radiation and that expand as a consequence of their internal sound speed. We further consider the dynamics of shells that sweep up a dense circumstellar or circumgalactic medium. We calculate the "momentum ratio" Mdot v/(L/c) in the shell limit and show that it can only significantly exceed ~2 if the effective optical depth of the shell to re-radiated FIR photons is much larger than unity. We discuss simple prescriptions for the properties of galactic outflows for use in large-scale cosmological simulations. We also briefly discuss applications to the dusty ejection episodes of massive stars, the disruption of giant molecular clouds, and AGN.

preprint2015arXiv

Forged in FIRE: cusps, cores, and baryons in low-mass dwarf galaxies

We present ultra-high resolution cosmological hydrodynamic simulations of $M_*\simeq10^{4-6.3}M_{\odot}$ dwarf galaxies that form within $M_{v}=10^{9.5-10}M_{\odot}$ dark matter halos. Our simulations rely on the FIRE implementation of star formation feedback and were run with high enough force and mass resolution to directly resolve stellar and dark matter structure on the ~200 pc scales of interest for classical and ultra-faint dwarfs in the Local Group. The resultant galaxies sit on the $M_*$ vs. $M_{v}$ relation required to match the Local Group stellar mass function. They have bursty star formation histories and also form with half-light radii and metallicities that broadly match those observed for local dwarfs at the same stellar mass. We demonstrate that it is possible to create a large (~1 kpc) dark matter core in a cosmological simulation of an $M_*\simeq10^{6.5}M_{\odot}$ dwarf galaxy that resides within an $M_{v}=10^{10}M_{\odot}$ halo -- precisely the scale of interest for resolving the Too Big to Fail problem. However, these large cores are not ubiquitous and appear to correlate closely with the star formation histories of the dwarfs: dark matter cores are largest in systems that form their stars late ($z\lesssim2$), after the early epoch of cusp building mergers has ended. Our $M_*\simeq10^4M_{\odot}$ dwarf retains a cuspy dark matter halo density profile that matches almost identically that of a dark-matter only run of the same system. Despite forming in a field environment, this very low mass dwarf has observable properties that match closely to those of ultra-faint satellite galaxies of the Milky Way, including a uniformly old stellar population (>10 Gyr). Though ancient, most of the stars in our ultra-faint form after reionization; the UV field acts mainly to suppress fresh gas accretion, not to boil away gas that is already present in the proto-dwarf.

preprint2015arXiv

Gusty, gaseous flows of FIRE: galactic winds in cosmological simulations with explicit stellar feedback

We present an analysis of the galaxy-scale gaseous outflows from the FIRE (Feedback in Realistic Environments) simulations. This suite of hydrodynamic cosmological zoom simulations resolves formation of star-forming giant molecular clouds to $z=0$, and features an explicit stellar feedback model on small scales. Our simulations reveal that high redshift galaxies undergo bursts of star formation followed by powerful gusts of galactic outflows that eject much of the ISM and temporarily suppress star formation. At low redshift, however, sufficiently massive galaxies corresponding to L*-progenitors develop stable disks and switch into a continuous and quiescent mode of star formation that does not drive outflows far into the halo. Mass-loading factors for winds in L*-progenitors are $η\approx 10$ at high redshift, but decrease to $η\ll 1$ at low redshift. Although lower values of $η$ are expected as halos grow in mass over time, we show that the strong suppression of outflows with decreasing redshift cannot be explained by mass evolution alone. Circumgalactic outflow velocities are variable and broadly distributed, but typically range between one and three times the circular velocity of the halo. Much of the ejected material builds a reservoir of enriched gas within the circumgalactic medium, some of which could be later recycled to fuel further star formation. However, a fraction of the gas that leaves the virial radius through galactic winds is never regained, causing most halos with mass $M_h \le 10^{12} M_{\odot}$ to be deficient in baryons compared to the cosmic mean by $z=0$.

preprint2015arXiv

Insights into Tidal Disruption of Stars from PS1-10jh

Was PS1-10jh (Gezari et al. 2012), an optical/UV transient discovered by the Pan-STARRS Medium Deep Survey, the tidal disruption of a star by a massive black hole? We address two aspects of the problem: the composition of the putative disrupted object (using the spectroscopic data), and the energetics of the observed gas and radiation (using the photometric data). We perform photoionization calculations and compare with the observed lower limit of the line ratio L(He II 4686)/L(Halpha) > 5 to argue that this event was not the disruption of a solar-type star, and instead was likely the disruption of a helium core (as first proposed by Gezari et al. 2012). Disruption of such a dense object requires a relatively small central BH, M_BH <~ 2 x 10^5 M_sun. We use the photometric data to infer that PS1-10jh comprised an outflow of ~ 0.01 M_sun of gas, escaping from the BH at ~1000 km/s, and we propose that this outflow was driven primarily by radiation pressure trapped by Thomson and resonance line scattering. The large ratio of radiated energy to kinetic energy, E_rad/E_K ~ 10^4, together with the large value of E_rad ~ 2 x 10^(51) erg, suggests that the outflow was shocked at large radius (perhaps similar to super-luminous supernovae or the internal shock model for gamma-ray bursts). We describe puzzles in the physics of PS1-10jh, and discuss how this event may help us understand future tidal disruptions and super-Eddington accretion events as well.

preprint2015arXiv

The Difficulty of Getting High Escape Fractions of Ionizing Photons from High-redshift Galaxies: a View from the FIRE Cosmological Simulations

We present a series of high-resolution (20-2000 Msun, 0.1-4 pc) cosmological zoom-in simulations at z~6 from the Feedback In Realistic Environment (FIRE) project. These simulations cover halo masses 10^9-10^11 Msun and rest-frame ultraviolet magnitude Muv = -9 to -19. These simulations include explicit models of the multi-phase ISM, star formation, and stellar feedback, which produce reasonable galaxy properties at z = 0-6. We post-process the snapshots with a radiative transfer code to evaluate the escape fraction (fesc) of hydrogen ionizing photons. We find that the instantaneous fesc has large time variability (0.01%-20%), while the time-averaged fesc over long time-scales generally remains ~5%, considerably lower than the estimate in many reionization models. We find no strong dependence of fesc on galaxy mass or redshift. In our simulations, the intrinsic ionizing photon budgets are dominated by stellar populations younger than 3 Myr, which tend to be buried in dense birth clouds. The escaping photons mostly come from populations between 3-10 Myr, whose birth clouds have been largely cleared by stellar feedback. However, these populations only contribute a small fraction of intrinsic ionizing photon budgets according to standard stellar population models. We show that fesc can be boosted to high values, if stellar populations older than 3 Myr produce more ionizing photons than standard stellar population models (as motivated by, e.g., models including binaries). By contrast, runaway stars with velocities suggested by observations can enhance fesc by only a small fraction. We show that "sub-grid" star formation models, which do not explicitly resolve star formation in dense clouds with n >> 1 cm^-3, will dramatically over-predict fesc.

preprint2015arXiv

Time Varying Dynamical Star Formation Rate

We present numerical evidence of dynamic star formation in which the accreted stellar mass grows superlinearly with time, roughly as $t^2$. We perform simulations of star formation in self-gravitating hydrodynamic and magneto-hydrodynamic turbulence that is continuously driven. By turning the self-gravity of the gas in the simulations on or off, we demonstrate that self-gravity is the dominant physical effect setting the mass accretion rate at early times before feedback effects take over, contrary to theories of turbulence-regulated star formation. We find that gravitational collapse steepens the density profile around stars, generating the power-law tail on what is otherwise a lognormal density probability distribution function. Furthermore, we find turbulent velocity profiles to flatten inside collapsing regions, altering the size-linewidth relation. This local flattening reflects enhancements of turbulent velocity on small scales, as verified by changes to the velocity power spectra. Our results indicate that gas self-gravity dynamically alters both density and velocity structures in clouds, giving rise to a time-varying star formation rate. We find that a substantial fraction of the gas that forms stars arrives via low density flows, as opposed to accreting through high density filaments.

preprint2014arXiv

Detection of warm and diffuse baryons in large scale structure from the cross-correlation of gravitational lensing and the thermal Sunyaev-Zeldovich effect

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-constructed 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.01\pm 0.31\pm 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.

preprint2014arXiv

Hot Galactic Winds Constrained by the X-Ray Luminosities of Galaxies

Galactic superwinds may be driven by very hot outflows generated by overlapping supernovae within the host galaxy. We use the Chevalier & Clegg (CC85) wind model and the observed correlation between X-ray luminosities of galaxies and their SFRs to constrain the mass loss rates (\dot{M}_hot) across a wide range of star formation rates (SFRs), from dwarf starbursts to ultra-luminous infrared galaxies. We show that for fixed thermalization efficiency and mass loading rate, the X-ray luminosity of the hot wind scales as L_X ~ SFR^2, significantly steeper than is observed for star-forming galaxies: L_X ~ SFR. Using this difference we constrain the mass-loading and thermalization efficiency of hot galactic winds. For reasonable values of the thermalization efficiency (<~ 1) and for SFR >~ 10 M_sun/yr we find that \dot{M}_hot/SFR <~ 1, significantly lower than required by integrated constraints on the efficiency of stellar feedback in galaxies, and potentially too low to explain observations of winds from rapidly star-forming galaxies. In addition, we highlight the fact that heavily mass-loaded winds cannot be described by the adiabatic CC85 model because they become strongly radiative.

preprint2014arXiv

On the Dynamics of Ultra Compact X-ray Binaries: 4U 1850-087, 4U 0513-40 and M15 X-2

In this work we extend our dynamical study of Ultra Compact X-ray Binaries (UCXB) 4U 1820-30 from Prodan and Murray 2012 to three more UCXBs in globular clusters: 4U 1850-087, 4U 0513-40 and M15 X-2. These three UCXBs have orbital periods < 20 mins. Two of them, 4U 1850-087 and 4U 0513-40, have suspected luminosity variations of order of ~ 1yr. There is insufficient observational data to make any statements regarding the long periodicity in the light curve of M15 X-2 at this point. The properties of these three systems are quite similar to 4U 1820-30, which prompt us to model their dynamics in the same manner. As in the case of 4U 1820-30, we interpret the suspected long periods as the period of small oscillations around a stable fixed point in the Kozai resonance. We provide a lower limit on the tidal dissipation factor Q which is in agreement with results obtained for the case of 4U 1820-30.

preprint2014arXiv

Radiation Feedback in ULIRGS: Are Photons Movers and Shakers?

We use our variable Eddington tensor (VET) radiation hydrodynamics code to perform two-dimensional simulations to study the impact of radiation forces on atmospheres composed of dust and gas. Our setup closely follows that of Krumholz & Thompson, assuming that dust and gas are well-coupled and that the radiation field is characterized by blackbodies with temperatures >~ 80 K, as might be found in ultraluminous infrared galaxies. In agreement with previous work, we find that Rayleigh-Taylor instabilities develop in radiation supported atmospheres, leading to inhomogeneities that limit momentum exchange between radiation and dusty gas, and eventually providing a near balance of the radiation and gravitational forces. However, the evolution of the velocity and spatial distributions of the gas differs significantly from previous work, which utilized a less accurate flux-limited diffusion (FLD) method. Our VET simulations show continuous net acceleration of the gas, with no steady-state reached by the end of the simulation. In contrast, FLD results show little net acceleration of the gas and settle in to a quasi-steady, turbulent state with low velocity dispersion. The discrepancies result primarily from the inability of FLD to properly model the variation of the radiation field around structures that are less than a few optical depths across. We conclude that radiation feedback remains a viable mechanism for driving high-Mach number turbulence. We discuss implications for observed systems and global numerical simulations of feedback, but more realistic setups are needed to make robust observational predictions and assess the prospect of launching outflows with radiation.

preprint2014arXiv

Secular Effects of Tidal Damping in Compact Planetary Systems

We describe the long-term evolution of compact systems of terrestrial planets, using a set of simulations that match the statistical properties of the observed exoplanet distribution. The evolution is driven by tidal dissipation in the planetary interiors, but the systems evolve as a whole due to secular gravitational interactions. We find that, for Earth-like dissipation levels, planetary orbits can be circularised out to periods of order 100 days, an order of magnitude larger than is possible for single planets. The resulting distribution of eccentricities is a qualitative match to that inferred from transit timing variations, with a minority of non-zero eccentricities maintained by particular secular configurations. The coupling of the tidal and secular processes enhance the inward migration of the innermost planets in these systems, and can drive them to short orbital periods. Resonant interactions of both the mean motion and secular variety are observed, although the interactions are not strong enough to drive systemic instability in most cases. However, we demonstrate that these systems can easily be driven unstable if coupled to giant planets on longer period orbits.

preprint2014arXiv

What controls star formation in the central 500 pc of the Galaxy?

The star formation rate (SFR) in the Central Molecular Zone (CMZ, i.e. the central 500 pc) of the Milky Way is lower by a factor of >10 than expected for the substantial amount of dense gas it contains, which challenges current star formation theories. In this paper, we quantify which physical mechanisms could be responsible. On scales larger than the disc scale height, the low SFR is found to be consistent with episodic star formation due to secular instabilities or possibly variations of the gas inflow along the Galactic bar. The CMZ is marginally Toomre-stable when including gas and stars, but highly Toomre-stable when only accounting for the gas, indicating a low condensation rate of self-gravitating clouds. On small scales, we find that the SFR in the CMZ may be caused by an elevated critical density for star formation due to the high turbulent pressure. The existence of a universal density threshold for star formation is ruled out. The HI-H$_2$ phase transition of hydrogen, the tidal field, a possible underproduction of massive stars due to a bottom-heavy initial mass function, magnetic fields, and cosmic ray or radiation pressure feedback also cannot individually explain the low SFR. We propose a self-consistent cycle of star formation in the CMZ, in which the effects of several different processes combine to inhibit star formation. The rate-limiting factor is the slow evolution of the gas towards collapse - once star formation is initiated it proceeds at a normal rate. The ubiquity of star formation inhibitors suggests that a lowered central SFR should be a common phenomenon in other galaxies. We discuss the implications for galactic-scale star formation and supermassive black hole growth, and relate our results to the star formation conditions in other extreme environments.

preprint2013arXiv

Accretion Does Not Drive the Turbulence in Galactic Disks

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

preprint2013arXiv

Black hole triple dynamics: breakdown of the orbit average approximation and implications for gravitational wave detections

Coalescing black hole (BH) binaries forming in the dense core of globular clusters (GCs) are expected to be one the brightest sources of gravitational wave (GW) radiation for the next generation of ground-based laser interferometers. Favorable conditions for merger are initiated by the Kozai resonance in which the gravitational interaction with a third distant object, typically another BH, induces quasi-periodic variations of the inner BH binary eccentricity. In this paper we perform high precision N-body simulations of the long term evolution of hierarchical BH triples and investigate the conditions that lead to the merging of the BH binary and the way it might become an observable source of GW radiation. We find that the secular orbit average treatment, adopted in previous works, does not reliably describe the dynamics of these systems if the binary is orbited by the outer BH on a highly inclined orbit at a moderate distance. During the high eccentricity phase of a Kozai cycle the torque due to the outer BH can drive the binary to extremely large eccentricities in a fraction of the binary's orbital period. This occurs before relativistic terms become important to the evolution and allows the binary GW signal to reach large GW frequencies (>~10 Hz) at high eccentricities. We show that ~50 % of coalescing BH binaries driven by the Kozai mechanism in GCs will have eccentricities larger than 0.1, with 10 % of them being extremely eccentric, (1-e)<~10^-5, when they first chirp in the frequency band of ground based laser interferometers. This implies that a large fraction of such GW sources could be missed if conventional quasi-circular templates are used for analysis of GW detectors data. The efficient detection of all coalescing BH binaries in GCs will therefore require template banks of eccentric inspiral waveforms for matched-filtering and dedicated search strategies.

preprint2013arXiv

Dense Molecular Gas: A Sensitive Probe of Stellar Feedback Models

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

preprint2013arXiv

On WD-WD Mergers in Triple Systems: The Role of Kozai Resonance with Tidal Friction

White dwarf-white dwarf (WD-WD) mergers may lead to type Ia supernovae events. Thompson (2011) suggested that many such binaries are produced in hierarchical triple systems. The tertiary induces eccentricity oscillations in the inner binary via the Kozai-Lidov mechanism, driving the binary to high eccentricities, and significantly reducing the gravitational wave merger timescale (T_GW) over a broad range of parameter space. Here, we investigate the role of tidal forces in these systems. We show that tidal effects are important in the regime of moderately high initial relative inclination between the inner binary and the outer tertiary. For 85 < i_0 < 90 degrees (prograde) and 97 < i_0 < 102 degrees (retrograde), tides combine with GW radiation to dramatically decrease T_GW. In the regime of high inclinations between 91 < i_0 < 96 degrees, the inner binary likely suffers a direct collision, as in the work of Katz & Dong (2012) and tidal effects do not play an important role.

preprint2013arXiv

Resolving The Generation of Starburst Winds in Galaxy Mergers

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 cooling, 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.

preprint2013arXiv

Scattered Emission from z~1 Galactic Outflows

Mapping Mg II resonance emission scattered by galactic winds offers a means to determine the spatial extent and density of the warm outflow. Using Keck/LRIS spectroscopy, we have resolved scattered Mg II emission to the east of 32016857, a star-forming galaxy at z =0.9392 with an outflow. The Mg II emission from this galaxy exhibits a P-Cygni profile, extends further than both the continuum and [O II] emission along the eastern side of the slit, and has a constant Doppler shift along the slit which does not follow the velocity gradient of the nebular [O II] emission. Using the Sobolev approximation, we derive the density of Mg+ ions at a radius of 12 to 18 kpc in the outflow. We model the ionization correction and find that much of the outflowing Mg is in Mg++. We estimate that the total mass flux could be as large as 330 - 500 solar masses per year, with the largest uncertainties coming from the depletion of Mg onto grains and the clumpiness of the warm outflow. We show that confining the warm clouds with a hot wind reduces the estimated mass flux of the warm outflow and indicates amass-loading factor near unity in the warm phase alone. Based on the high blue luminosities that distinguish 32016857 and TKRS 4389, described by Rubin et al. 2011, from other galaxies with P-Cygni emission, we suggest that, as sensitivity to diffuse emission improves, scattering halos may prove to be a generic property of star-forming galaxies at intermediate redshifts.

preprint2013arXiv

Star Formation in Galaxy Mergers with Realistic Models of Stellar Feedback & the Interstellar Medium

We use simulations with realistic models for stellar feedback to study galaxy mergers. These high resolution (1 pc) simulations follow formation and destruction of individual GMCs and star clusters. The final starburst is dominated by in situ star formation, fueled by gas which flows inwards due to global torques. The resulting high gas density results in rapid star formation. The gas is self gravitating, and forms massive (~10^10 M_sun) GMCs and subsequent super-starclusters (masses up to 10^8 M_sun). However, in contrast to some recent simulations, the bulk of new stars which eventually form the central bulge are not born in superclusters which then sink to the center of the galaxy, because feedback efficiently disperses GMCs after they turn several percent of their mass into stars. Most of the mass that reaches the nucleus does so in the form of gas. The Kennicutt-Schmidt law emerges naturally as a consequence of feedback balancing gravitational collapse, independent of the small-scale star formation microphysics. The same mechanisms that drive this relation in isolated galaxies, in particular radiation pressure from IR photons, extend over seven decades in SFR to regulate star formation in the most extreme starbursts (densities >10^4 M_sun/pc^2). Feedback also drives super-winds with large mass loss rates; but a significant fraction of the wind material falls back onto the disks at later times, leading to higher post-starburst SFRs in the presence of stellar feedback. Strong AGN feedback is required to explain sharp cutoffs in star formation rate. We compare the predicted relic structure, mass profile, morphology, and efficiency of disk survival to simulations which do not explicitly resolve GMCs or feedback. Global galaxy properties are similar, but sub-galactic properties and star formation rates can differ significantly.

preprint2013arXiv

The Meaning and Consequences of Star Formation Criteria in Galaxy Models with Resolved Stellar Feedback

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.

preprint2012arXiv

A Giant Metrewave Radio Telescope/Chandra view of IRAS 09104+4109: A type 2 QSO in a cooling flow

IRAS 09104+4109 is a rare example of a dust enshrouded type 2 QSO in the centre of a cool-core galaxy cluster. Previous observations of this z=0.44 system showed that as well as powering the hyper-luminous infrared emission of the cluster-central galaxy, the QSO is associated with a double-lobed radio source. However, the steep radio spectral index and misalignment between the jets and ionised optical emission suggested that the orientation of the QSO had recently changed. We use a combination of new, multi-band Giant Metrewave Radio Telescope observations and archival radio data to confirm that the jets are no longer powered by the QSO, and estimate their age to be 120-160 Myr. This is in agreement with the ~70-200 Myr age previously estimated for star-formation in the galaxy. Previously unpublished Very Long Baseline Array data reveal a 200 pc scale double radio source in the galaxy core which is more closely aligned with the current QSO axis and may represent a more recent period of jet activity. These results suggest that the realignment of the QSO, the cessation of jet activity, and the onset of rapid star-formation may have been caused by a gas-rich galaxy merger. A Chandra X-ray observation confirms the presence of cavities associated with the radio jets, and we estimate the energy required to inflate them to be ~7.7x10^60 erg. The mechanical power of the jets is sufficient to balance radiative cooling in the cluster, provided they are efficiently coupled to the intra-cluster medium (ICM). We find no evidence of direct radiative heating and conclude that the QSO either lacks the radiative luminosity to heat the ICM, or that it requires longer than 100-200 Myr to significantly impact its environment. [Abridged]

preprint2012arXiv

Dark Matter Substructure Detection Using Spatially Resolved Spectroscopy of Lensed Dusty Galaxies

We investigate how strong lensing of dusty, star-forming galaxies by foreground galaxies can be used as a probe of dark matter halo substructure. We find that spatially resolved spectroscopy of lensed sources allows dramatic improvements to measurements of lens parameters. In particular we find that modeling of the full, three-dimensional (angular position and radial velocity) data can significantly facilitate substructure detection, increasing the sensitivity of observables to lower mass subhalos. We carry out simulations of lensed dusty sources observed by early ALMA (Cycle 1) and use a Fisher matrix analysis to study the parameter degeneracies and mass detection limits of this method. We find that, even with conservative assumptions, it is possible to detect galactic dark matter subhalos of ~ 10^8 M_{\odot} with high significance in most lensed DSFGs. Specifically, we find that in typical DSFG lenses, there is a ~ 55 % probability of detecting a substructure with M>10^8 M_{\odot} with more than 5 sigma detection significance in each lens, if the abundance of substructure is consistent with previous lensing results. The full ALMA array, with its significantly enhanced sensitivity and resolution, should improve these estimates considerably. Given the sample of ~100 lenses provided by surveys like the South Pole Telescope, our understanding of dark matter substructure in typical galaxy halos is poised to improve dramatically over the next few years.

preprint2012arXiv

Stellar Feedback & Bulge Formation in Clumpy Disks

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.

preprint2012arXiv

Stellar Feedback in Galaxies and the Origin of Galaxy-scale Winds

Feedback from massive stars is believed to play a critical role in driving galactic super-winds that enrich the IGM and shape the galaxy mass function and mass-metallicity relation. In previous papers, we introduced new numerical methods for implementing stellar feedback on sub-GMC through galactic scales in galaxy simulations. This includes radiation pressure (UV through IR), SNe (Type-I & II), stellar winds ('fast' O-star through 'slow' AGB winds), and HII photoionization. Here, we show that these feedback mechanisms drive galactic winds with outflow rates as high as ~10-20 times the galaxy SFR. The mass-loading efficiency (wind mass loss rate divided by SFR) scales inversely with circular velocity, consistent with momentum-conservation expectations. We study the contributions of each feedback mechanism to galactic winds in a range of galaxy models, from SMC-like dwarfs & MW-analogues to z~2 clumpy disks. In massive, gas-rich systems (local starbursts and high-z galaxies), radiation pressure dominates the wind generation. For MW-like spirals and dwarf galaxies the gas densities are much lower, and shock-heated gas from SNe and stellar winds dominates production of large-scale outflows. In all models, however, winds have a multi-phase structure that depends on interactions between multiple feedback mechanisms operating on different spatial & time scales: any single mechanism fails to reproduce the winds observed. We provide fitting functions for wind mass-loading and velocities as a function of galaxy properties, for use in cosmological simulations and semi-analytic models. These differ from typically-adopted formulae with explicit dependence on gas surface density that can be very important in both low-density dwarf galaxies and high-density gas-rich galaxies.

preprint2012arXiv

The Structure of the Interstellar Medium of Star Forming Galaxies

We present numerical methods for including stellar feedback in galaxy-scale simulations. We include heating by SNe (I & II), gas recycling and shock-heating from O-star & AGB winds, HII photoionization, and radiation pressure from stellar photons. The energetics and time-dependence are taken directly from stellar evolution models. We implement these in simulations with pc-scale resolution, modeling galaxies from SMC-like dwarfs and MW analogues to massive z~2 starburst disks. Absent feedback, gas cools and collapses without limit. With feedback, the ISM reaches a multi-phase steady state in which GMCs continuously form, disperse, and re-form. Our primary results include: (1) Star forming galaxies generically self-regulate at Toomre Q~1. Most of the volume is in diffuse hot gas with most of the mass in dense GMC complexes. The phase structure and gas mass at high densities are much more sensitive probes of stellar feedback physics than integrated quantities (Toomre Q or gas velocity dispersion). (2) Different feedback mechanisms act on different scales: radiation & HII pressure are critical to prevent runaway collapse of dense gas in GMCs. SNe and stellar winds dominate the dynamics of volume-filling hot gas; however this primarily vents out of the disk. (3) The galaxy-averaged SFR is determined by feedback. For given feedback efficiency, restricting star formation to molecular gas or modifying the cooling function has little effect; but changing feedback mechanisms directly translates to shifts off the Kennicutt-Schmidt relation. (4) Self-gravity leads to marginally-bound GMCs with an ~M^-2 mass function with a cutoff at the Jeans mass; they live a few dynamical times before being disrupted by stellar feedback and turn ~1-10% of their mass into stars (increasing from dwarfs through starburst galaxies). Low-mass GMCs are preferentially unbound.

preprint2011arXiv

Broadband Transmission Spectroscopy of the super-Earth GJ 1214b suggests a Low Mean Molecular Weight Atmosphere

We used WIRCam on CFHT to observe four transits of the super-Earth GJ 1214b in the near-infrared. For each transit we observed in two bands nearly-simultaneously by rapidly switching the WIRCam filter wheel back and forth for the duration of the observations. By combining all our J-band (~1.25 microns) observations we find a transit depth in this band of 1.338\pm0.013% - a value consistent with the optical transit depth reported by Charbonneau and collaborators. However, our best-fit combined Ks-band (~2.15 microns) transit depth is deeper: 1.438\pm0.019%. Formally our Ks-band transits are deeper than the J-band transits observed simultaneously by a factor of 1.072\pm0.018 - a 4-sigma discrepancy. The most straightforward explanation for our deeper Ks-band depth is a spectral absorption feature from the limb of the atmosphere of the planet; for the spectral absorption feature to be this prominent the atmosphere of GJ 1214b must have a large scale height and a low mean molecular weight. That is, it would have to be hydrogen/helium dominated and this planet would be better described as a mini-Neptune. However, recently published observations from 0.78 - 1.0 microns, by Bean and collaborators, show a lack of spectral features and transit depths consistent with those obtained by Charbonneau and collaborators. The most likely atmospheric composition for GJ 1214b that arises from combining all these observations is less clear; if the atmosphere of GJ 1214b is hydrogen/helium dominated then it must have either a haze layer that is obscuring transit depth differences at shorter wavelengths, or significantly different spectral features than current models predict. Our observations disfavour a water-world composition, but such a composition will remain a possibility until observations reconfirm our deeper Ks-band transit depth or detect features at other wavelengths. [Abridged]

preprint2011arXiv

On the Dynamics and Tidal Dissipation Rate of the White Dwarf in 4U 1820-30

It has been suggested that the 170 day period in the light curve of the low mass X-ray binary 4U 1820-30 arises from the presence of a third body with a large inclination to the binary orbit. We show that this long period motion arises if the system is librating around the stable fixed point in a Kozai resonance. We demonstrate that mass transfer drives the system toward this fixed point, and calculate, both analytically and via numerical integrations, that the period of libration is of order 170 days when the mutual inclination is near the Kozai critical value. The non-zero eccentricity of the binary, combined with tidal dissipation, implies that the rate of change of the binary period would be slower than, or even of opposite sign to, that implied by standard mass transfer models. If the 170 day period results from libration, then, contrary to appearances, the orbital period of the inner binary is increasing with time; in that case, (e/0.009)^2Q/k_2 > 2.5 x 10^9, where k_2 = 0.01 is the tidal Love number and e = 0.009 is the fiducial eccentricity of the inner binary. It appears unlikely that the observed negative period derivative results from the smaller than expected (but positive) value of \dot P combined with the previously suggested acceleration of the system in the gravitational field of the host globular cluster NGC 6624. The discrepancy between the observed and expected period derivative requires further investigation.

preprint2011arXiv

Self-Regulated Star Formation in Galaxies via Momentum Input from Massive Stars

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 paper, 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.

preprint2011arXiv

Star Formation in Massive Clusters via Bondi Accretion

Essentially all stars form in giant molecular clouds (GMCs). However, inside GMCs, most of the gas does not participate in star formation; rather, denser gas accumulates in clumps in the GMC, with the bulk of the stars in a given GMC forming in a few of the most massive clumps. In the Milky Way, these clumps have masses $M_{\rm cl}\lesssim 5\times 10^{-2}$ of the GMC, radii $r_{\rm cl} \sim 1$pc, and free-fall times $τ_{\rm cl} \sim 2\times 10^5\yr$. We show that clumps inside giant molecular clouds should accrete at a modified Bondi accretion rate, which depends on clump mass as $\dot M_{\rm cl}\sim M_{\rm cl}^{5/4}$. This rate is initially rather slow, usually slower than the initial star formation rate inside the clump (we adopt the common assumption that inside the clump, $\dot M_*=ε_{\rm ff} M_{\rm cl}/τ_{\rm cl}$, with $ε_{\rm ff} \approx 0.017$). However, after $\sim 2$ GMC free-fall times $τ_{\rm GMC}$, the clump accretion rate accelerates rapidly; formally, the clump can accrete the entire GMC in $\sim 3τ_{\rm GMC}$. At the same time, the star formation rate accelerates, tracking the Bondi accretion rate. If the GMC is disrupted by feedback from the largest clump, half the stars in that clump form in the final $\taug$ before the GMC is disrupted. The theory predicts that the distribution of effective star formation rates, measured per GMC free-fall time, is broad, ranging from $\sim 0.001$ up to 0.1 or larger and that the mass spectrum of star clusters is flatter than that of clumps, consistent with observations.

preprint2010arXiv

Low-ionization Line Emission from Starburst Galaxies: A New Probe of Galactic-Scale Outflows

We study the kinematically narrow, low-ionization line emission from a bright, starburst galaxy at z = 0.69 using slit spectroscopy obtained with Keck/LRIS. The spectrum reveals strong absorption in MgII and FeII resonance transitions with Doppler shifts of -200 to -300 km/s, indicating a cool gas outflow. Emission in MgII near and redward of systemic velocity, in concert with the observed absorption, yields a P Cygni-like line profile similar to those observed in the Ly alpha transition in Lyman Break Galaxies. Further, the MgII emission is spatially resolved, and extends significantly beyond the emission from stars and HII regions within the galaxy. Assuming the emission has a simple, symmetric surface brightness profile, we find that the gas extends to distances > ~7 kpc. We also detect several narrow FeII* fine-structure lines in emission near the systemic velocity, arising from energy levels which are radiatively excited directly from the ground state. We suggest that the MgII and FeII* emission is generated by photon scattering in the observed outflow, and emphasize that this emission is a generic prediction of outflows. These observations provide the first direct constraints on the minimum spatial extent and morphology of the wind from a distant galaxy. Estimates of these parameters are crucial for understanding the impact of outflows in driving galaxy evolution.

preprint2010arXiv

Massive Star Forming Regions in the Galaxy using the Spitzer GLIMPSE Survey

We examine the thirteen most luminous sources in the WMAP free-free map using the Spitzer GLIMPSE and MSX surveys to identify massive star formation regions, emitting one-third of the Galactic free-free luminosity. We identify star forming regions by a combination of bubble morphology in 8 $\micronm$ (PAH) emission and radio recombination line radial velocities. We find 40 star forming regions associated with our WMAP sources, and determine unique distances to 31. We interpret the bubbles as evidence for radial expansion. The radial velocity distribution for each source allows us to measure the intrinsic speed of a region's expansion. This speed is consistent with the size and age of the bubbles. The high free-free luminosities, combined with negligible synchrotron emission, demonstrate that the bubbles are not driven by supernovae. The kinetic energy of the largest bubbles is a substantial fraction of that measured in the older superbubbles found by Heiles. We find that the energy injected into the ISM by our bubbles is similar to that required to maintain the turbulent motion in the gas disk inside 8 kpc. We report a number of new star forming regions powered by massive ($\textrm{M}_{*} > 10^4 \textrm{M}_\sun$) star clusters. We measure the scale height of the Galactic O stars to be $h_{\textrm{*}} = 35 \pm 5 \pc$. We determine an empirical relationship between the PAH and free-free emission of the form $F_{\textrm{PAH}} \propto F^2_{\textrm{ff}}$. Finally, we find that the bubble geometry is more consistent with a spherical shell rather than a flattened disk.

preprint2010arXiv

Near-infrared Thermal Emission from WASP-12b: detections of the secondary eclipse in Ks, H & J

We present Ks, H & J-band photometry of the very highly irradiated hot Jupiter WASP-12b using the Wide-field Infrared Camera on the Canada-France-Hawaii telescope. Our photometry brackets the secondary eclipse of WASP-12b in the Ks and H-bands, and in J-band starts in mid-eclipse and continues until well after the end of the eclipse. We detect its thermal emission in all three near-infrared bands. Our secondary eclipse depths are 0.309 +/- 0.013% in Ks-band (24-sigma), 0.176 +/- 0.020% in H-band (9-sigma) and 0.131 +/- 0.028% in J-band (4-sigma). All three secondary eclipses are best-fit with a consistent phase that is compatible with a circular orbit. By combining our secondary eclipse times with others published in the literature, as well as the radial velocity and transit timing data for this system, we show that there is no evidence that WASP-12b is precessing at a detectable rate, and show that its orbital eccentricity is likely zero. Our thermal emission measurements also allow us to constrain the characteristics of the planet's atmosphere; our Ks-band eclipse depth argues in favour of inefficient day to nightside redistribution of heat and a low Bond albedo for this very highly irradiated hot Jupiter. The J and H-band brightness temperatures are slightly cooler than the Ks-band brightness temperature, and thus hint at the possibility of a modest temperature inversion deep in the atmosphere of WASP-12b; the high pressure, deep atmospheric layers probed by our J and H-band observations are likely more homogenized than the higher altitude layer. Lastly, our best-fit Ks-band eclipse has a marginally longer duration than would otherwise be expected; this may be tentative evidence for material being tidally stripped from the planet - as was predicted for this system by Li & collaborators, and for which observational confirmation was recently arguably provided by Fossati & collaborators.

preprint2010arXiv

Radiation pressure from massive star clusters as a launching mechanism for super-galactic winds

Galactic outflows of low ionization, cool gas are ubiquitous in local starburst galaxies, and in the majority of galaxies at high redshift. How these cool outflows arise is still in question. Hot gas from supernovae has long been suspected as the primary driver, but this mechanism suffers from its tendency to destroy the cool gas as the latter is accelerated. We propose a modification of the supernova scenario that overcomes this difficulty. Star formation is observed to take place in clusters; in a given galaxy, the bulk of the star formation is found in the ~20 most massive clusters. We show that, for L* galaxies, the radiation pressure from clusters with M>10^6 M_sun is able to expel the surrounding gas at velocities in excess of the circular velocity of the disk galaxy. This cool gas can travel above the galactic disk in less than 2 Myr, well before any supernovae erupt in the driving cluster. Once above the disk, the cool outflowing gas is exposed to radiation, and supernovae induced hot gas outflows, from other clusters in the disk, which drive it to distances of several tens to hundreds of kpc. Because the radiatively driven clouds grow in size as they travel, and because the hot gas is more dilute at large distance, the clouds are less subject to destruction if they do eventually encounter hot gas. Therefore, unlike wind driven clouds, radiatively driven clouds can survive to distances ~50 kpc. We identify these cluster-driven winds with large-scale galactic outflows. Another implication of our model is that only starburst galaxies, where massive clusters reside, are able to drive winds cold outflows on galactic scales via this mechanism. We find that the critical star formation rates above which large scale cool outflows will be launched to be ~0.1 M_sun/yr/kpc^2, which is in good agreement with observations.

preprint2009arXiv

A Maximum Stellar Surface Density in Dense Stellar Systems

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

preprint2009arXiv

Compact High-Redshift Galaxies Are the Cores of the Most Massive Present-Day Spheroids

Observations suggest that effective radii of high-z massive spheroids are as much as a factor ~6 smaller than low-z galaxies of comparable mass. Given the apparent absence of low-z counterparts, this has often been interpreted as indicating that the high density, compact red galaxies must be 'puffed up' by some mechanism. We compare the ensemble of high-z observations with large samples of well-observed low-z ellipticals. At the same physical radii, the stellar surface mass densities of low and high-z systems are comparable. Moreover, the abundance of high surface density material at low redshift is comparable to or larger than that observed at z>1-2, consistent with the continuous buildup of spheroids over this time. The entire population of compact, high-z red galaxies may be the progenitors of the high-density cores of present-day ellipticals, with no need for a decrease in stellar density from z=2 to z=0. The primary difference between low and high-z systems is thus the observed low-density material at large radii in low-z spheroids (rather than the high-density material in high-z spheroids). Such low-density material may either (1) assemble at z<2 or (2) be present, but not yet detected, at z>2. Mock observations of low-z massive systems show that the high-z observations do not yet probe sufficiently low surface brightness material to detect the low surface density 'wings' (if present). Thus, if the high-z galaxies resemble the most massive systems today, their inferred effective radii could be under-estimated by factors ~2-4. This difference arises because massive systems at low redshift are not well-fit by single Sersic profiles. We discuss implications of our results for physical models of galaxy evolution.

preprint2009arXiv

Radio Emission from Supernova Remnants: Implications for Post-Shock Magnetic Field Amplification and the Magnetic Fields of Galaxies

Using observations from the literature, we show that the non-thermal radio luminosity (L) of supernova remnants (SNRs) is a strong function of the average gas surface density (Sigma) of the galaxy in which the remnants reside, from normal spirals to luminous starbursts. We combine a simple theory for electron cooling in SNRs with the observed radio luminosities to estimate the remnant magnetic field strength (B_SNR): the correlation between L and Sigma implies that B_SNR also increases with Sigma. We explore two interpretations of this correlation: (1) B_SNR is generated by post-shock magnetic field amplification, with B_SNR^2 proportional to Sigma and (2) B_SNR results from shock-compression of the ambient ISM magnetic field (B_ISM), with B_ISM being larger in denser galaxies. We find that shock compression is, on average, sufficient to produce the observed radio emission from SNRs in the densest starbursts; amplification of post-shock magnetic fields is not required. By contrast, in normal spirals post-shock field amplification (by a factor of a few - 10) is consistent with the data; we find tentative evidence that both the Alfven speed and the ratio of B_SNR^2 to the post-shock pressure ("epsilon_B") are constant in SNRs from galaxy to galaxy. We discuss observational tests that can be used to distinguish between these two interpretations of the radio luminosities of SNRs. Regardless of which is correct, the radio emission from SNRs provides an upper limit to B_ISM that is independent of the minimum energy assumption. For the densest starbursts, the ISM magnetic energy density is below that required for hydrostatic equilibrium; thus magnetic fields are not dynamically important on the largest scales in starbursts, in contrast with spiral galaxies like our own. This dichotomy may have implications for galactic dynamo theory.

preprint2009arXiv

The Disruption of Giant Molecular Clouds by Radiation Pressure and the Efficiency of Star Formation in Galaxies

Star formation is slow, in the sense that the gas consumption time is much longer than the dynamical time. It is also inefficient; essentially all star formation in local galaxies takes place in giant molecular clouds (GMCs), but the fraction of a GMC converted to stars is very small, ~5%. In the most luminous starbursts, the GMC lifetime is shorter than the main sequence lifetime of even the most massive stars, so that supernovae can play no role in GMC disruption. We investigate the disruption of GMCs across a wide range of galaxies, from normal spirals to the densest starbursts; we take into account the effects of HII gas pressure, shocked stellar winds, protostellar jets, and radiation pressure produced by the absorption and scattering of starlight on dust grains. In the Milky Way, we find that a combination of three mechanisms, jets, HII gas pressure, and radiation pressure, disrupts the clouds. In more rapidly star forming galaxies such as ``clump'' galaxies at high-redshift, ultra-luminous infrared galaxies (ULIRGs) and submillimeter galaxies, radiation pressure dominates natal cloud distribution. We predict the presence of 10-20 clusters with masses ~10^7 Msun in local ULIRGs such as Arp 220 and a similar number of clusters with M_* ~ 10^8 Msun in high redshift clump galaxies; submillimeter galaxies will have even more massive clusters. We find that the mass fraction of a GMC that ends up in stars is an increasing function of the gas surface density of a galaxy, reaching ~35% in the most luminous starbursts. Furthermore, the disruption of bubbles by radiation pressure stirs the interstellar medium to velocities of ~10 km/s in normal galaxies and to ~100 km/s in ULIRGs like Arp 220, consistent with observations. Thus, radiation pressure may play a dominant role in the ISM of star-forming galaxies.

preprint2009arXiv

The Small Scatter in BH-Host Correlations & The Case for Self-Regulated BH Growth

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.

preprint2006arXiv

The Ionization State of Sodium in Galactic Winds

Roughly 80% of Ultraluminous Infrared Galaxies (ULIRGs) show blue shifted absorption in the resonance lines of neutral sodium, indicating that cool winds are common in such objects, as shown by Rupke et al and by Martin. The neutral sodium (NaI) columns indicated by these absorption lines are ~ 10^{13}-3x10^{14}/cm^2, while the bolometric luminosity varies by a factor of only four. We show that the gas in ULIRG outflows is likely to be in photoionization equilibrium. The very small ULIRG sample of Goldader et al. demonstrates that the ratio of ultraviolet flux to far infrared flux varies by a factor $\sim100$ from object to object. While the Goldader sample does not overlap with those of Rupke et al. and Martin, we show that such a large variation in ultraviolet flux will produce a similar variation in the column of neutral sodium for a fixed mass flux and density. However, if the cold gas is in pressure equilibrium with a hot outflow with a mass loss rate similar to the star formation rate, the range of ionization state is significantly smaller. Measurements of the UV flux for objects in the Martin and Rupke et al. catalogs will definitively determine if photoionization effects are responsible for the wide variation seen in the sodium columns. If they are, a determination of the gas density and mass loss rate in the cool winds will follow, with attendant improvements in our understanding of wind driving mechanisms and of the effects of galaxies on their surroundings.

preprint2002arXiv

Hot Settling Accretion Flow onto a Spinning Black Hole

We study the structure and properties of hot MHD accretion onto a Kerr black hole. In such a system, the hole is magnetically coupled to the inflowing gas and exerts a torque onto the accretion flow. A hot settling flow can form around the hole and transport the angular momentum outward, to the outer edge of the flow. Unlike other hot flows, such as advection- and convection-dominated flows and inflow-outflow solutions (ADAFs, CDAFs, and ADIOS), the properties of the hot settling flow are determined by the spin of the central black hole, but are insensitive to the mass accretion rate. Therefore, it may be possible to identify rapidly spinning BHs simply from their broad-band spectra. Observationally, the hot settling flow around a Kerr hole is somewhat similar to other hot flows in that they all have hard, power-law spectra and relatively low luminosities. Thus, most black hole candidates in the low/hard and, perhaps, intermediate X-ray state may potentially accrete via the hot settling flow. However, a settling flow will be somewhat more luminous than ADAFs/CDAFs/ADIOS, will exhibit high variability in X-rays, and may have relativistic jets. This suggests that galactic microquasars and active galactic nuclei may be powered by hot settling flows. We identify several galactic X-ray sources as the best candidates.

Norman Murray

What is connected

Connect this record

See the researcher in context

Building this map preview

58 published item(s)

COMAP Early Science: IV. Power Spectrum Methodology and Results

COMAP Early Science: V. Constraints and Forecasts at $z \sim 3$

Gas infall and radial transport in cosmological simulations of Milky Way-mass disks

Kinematics and Feedback in H II regions in the Dwarf Starburst Galaxy IC 10

COMAP Early Science: I. Overview

COMAP Early Science: III. CO Data Processing

COMAP Early Science: VII. Prospects for CO Intensity Mapping at Reionization

The Mass of the Milky Way from the H3 Survey

Properties of the Circumgalactic Medium in Cosmic Ray-Dominated Galaxy Halos

Radiative Stellar Feedback in Galaxy Formation: Methods and Physics

Self-consistent proto-globular cluster formation in cosmological simulations of high-redshift galaxies

Swirls of FIRE: Spatially Resolved Gas Velocity Dispersions and Star Formation Rates in FIRE-2 Disk Environments

Binary Stars Can Provide the "Missing Photons" Needed for Reionization

Detection of lensing substructure using ALMA observations of the dusty galaxy SDP.81

Giant clumps in the FIRE simulations: a case study of a massive high-redshift galaxy

Observational Evidence of Dynamic Star Formation Rate in Milky Way Giant Molecular Clouds

Stellar and Quasar Feedback in Concert: Effects on AGN Accretion, Obscuration, and Outflows

Asynchronous rotation of Earth-mass planets in the habitable zone of lower-mass stars

Dynamics of Dusty Radiation Pressure Driven Shells and Clouds: Fast Outflows from Galaxies, Star Clusters, Massive Stars, and AGN

Forged in FIRE: cusps, cores, and baryons in low-mass dwarf galaxies

Gusty, gaseous flows of FIRE: galactic winds in cosmological simulations with explicit stellar feedback

Insights into Tidal Disruption of Stars from PS1-10jh

The Difficulty of Getting High Escape Fractions of Ionizing Photons from High-redshift Galaxies: a View from the FIRE Cosmological Simulations

Time Varying Dynamical Star Formation Rate

Detection of warm and diffuse baryons in large scale structure from the cross-correlation of gravitational lensing and the thermal Sunyaev-Zeldovich effect

Hot Galactic Winds Constrained by the X-Ray Luminosities of Galaxies

On the Dynamics of Ultra Compact X-ray Binaries: 4U 1850-087, 4U 0513-40 and M15 X-2

Radiation Feedback in ULIRGS: Are Photons Movers and Shakers?

Secular Effects of Tidal Damping in Compact Planetary Systems

What controls star formation in the central 500 pc of the Galaxy?

Accretion Does Not Drive the Turbulence in Galactic Disks

Black hole triple dynamics: breakdown of the orbit average approximation and implications for gravitational wave detections

Dense Molecular Gas: A Sensitive Probe of Stellar Feedback Models

On WD-WD Mergers in Triple Systems: The Role of Kozai Resonance with Tidal Friction

Resolving The Generation of Starburst Winds in Galaxy Mergers

Scattered Emission from z~1 Galactic Outflows

Star Formation in Galaxy Mergers with Realistic Models of Stellar Feedback & the Interstellar Medium

The Meaning and Consequences of Star Formation Criteria in Galaxy Models with Resolved Stellar Feedback

A Giant Metrewave Radio Telescope/Chandra view of IRAS 09104+4109: A type 2 QSO in a cooling flow

Dark Matter Substructure Detection Using Spatially Resolved Spectroscopy of Lensed Dusty Galaxies

Stellar Feedback & Bulge Formation in Clumpy Disks

Stellar Feedback in Galaxies and the Origin of Galaxy-scale Winds

The Structure of the Interstellar Medium of Star Forming Galaxies

Broadband Transmission Spectroscopy of the super-Earth GJ 1214b suggests a Low Mean Molecular Weight Atmosphere

On the Dynamics and Tidal Dissipation Rate of the White Dwarf in 4U 1820-30

Self-Regulated Star Formation in Galaxies via Momentum Input from Massive Stars

Star Formation in Massive Clusters via Bondi Accretion

Low-ionization Line Emission from Starburst Galaxies: A New Probe of Galactic-Scale Outflows

Massive Star Forming Regions in the Galaxy using the Spitzer GLIMPSE Survey

Near-infrared Thermal Emission from WASP-12b: detections of the secondary eclipse in Ks, H & J

Radiation pressure from massive star clusters as a launching mechanism for super-galactic winds

A Maximum Stellar Surface Density in Dense Stellar Systems

Compact High-Redshift Galaxies Are the Cores of the Most Massive Present-Day Spheroids

Radio Emission from Supernova Remnants: Implications for Post-Shock Magnetic Field Amplification and the Magnetic Fields of Galaxies

The Disruption of Giant Molecular Clouds by Radiation Pressure and the Efficiency of Star Formation in Galaxies

The Small Scatter in BH-Host Correlations & The Case for Self-Regulated BH Growth

The Ionization State of Sodium in Galactic Winds

Hot Settling Accretion Flow onto a Spinning Black Hole