Source author record

Zoltan Haiman

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

ResearcherUnclaimed source record

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

Catalog footprint

What is connected

76works

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.

preprint2026arXiv

Lomb-Scargle periodograms struggle with non-sinusoidal supermassive BH binary signatures in quasar lightcurves

Supermassive black hole binary (SMBHB) systems are expected to form as a consequence of galaxy mergers. At sub-parsec separations, SMBHBs can be identified as quasars with periodic variability with previous periodicity searches uncovering significant candidates. However, these searches focused primarily on sinusoidal signals, while theoretical models and hydrodynamical simulations predict that binaries produce more complex non-sinusoidal pulse shapes. Here we examine the efficacy of the Lomb-Scargle periodogram (LSP; one of the most popular tools for periodicity searches in unevenly sampled lightcurves) to detect periodicities with a sawtooth shape mimicking results of hydrodynamical simulations. We simulate idealised well-sampled lightcurves, lightcurves that mimic the data in the Palomar Transient Factory (PTF) analyzed in Charisi et al., 2016, and lightcurves that resemble our expectations for single-band data in the upcoming Legacy Survey of Space and Time (LSST) of the Rubin Observatory. We approximate quasar variability with a damped random walk (DRW) model, inject sinusoidal and sawtooth pulse shapes and assess their statistical significance. We find that in the presence of red noise the LSP detects a relatively low fraction of the sinusoidal signals (~45%, ~24% and ~23%, in the PTF-like, idealised, and LSST-like lightcurves, respectively). The fraction is significantly reduced for sawtooth periodicity (with only ~9% in PTF-like and ~1% in idealised and LSST-like lightcurves). These low recovery rates imply that previous searches have missed the large majority of binaries. They also have significant implications for the detection of SMBHBs in upcoming LSST necessitating the developement of advanced tools that go beyond the simple LSP.

preprint2022arXiv

Forming massive seed black holes in high-redshift quasar host progenitors

The presence of massive black holes (BHs) with masses of order $10^9\rm\, M_\odot$, powering bright quasars when the Universe was less than 1 Gyr old, poses strong constraints on their formation mechanism. Several scenarios have been proposed to date to explain massive BH formation, from the low-mass seed BH remnants of the first generation of stars to the massive seed BHs resulting from the rapid collapse of massive gas clouds. However, the plausibility of some of these scenarios to occur within the progenitors of high-z quasars has not yet been thoroughly explored. In this work, we investigate, by combining dark-matter only N-body simulations with a semi-analytic framework, whether the conditions for the formation of massive seed BHs from synchronised atomic-cooling halo pairs and/or dynamically-heated mini-haloes are fulfilled in the overdense regions where the progenitors of a typical high-redshift quasar host form and evolve. Our analysis shows that the peculiar conditions in such regions, i.e. strong halo clustering and high star formation rates, are crucial to produce a non-negligible number of massive seed BH host candidates: we find $\approx1400$ dynamically heated metal-free mini-haloes, including one of these which evolves to a synchronised pair and ends up in the massive quasar-host halo by $z=6$. This demonstrates that the progenitors of high-redshift quasar host haloes can harbour early massive seed BHs. Our results further suggest that multiple massive seed BHs may form in or near the quasar host's progenitors, potentially merging at lower redshifts and yielding gravitational wave events.

preprint2022arXiv

The Effect of Mission Duration on LISA Science Objectives

The science objectives of the LISA mission have been defined under the implicit assumption of a 4 yr continuous data stream. Based on the performance of LISA Pathfinder, it is now expected that LISA will have a duty cycle of $\approx 0.75$, which would reduce the effective span of usable data to 3 yr. This paper reports the results of a study by the LISA Science Group, which was charged with assessing the additional science return of increasing the mission lifetime. We explore various observational scenarios to assess the impact of mission duration on the main science objectives of the mission. We find that the science investigations most affected by mission duration concern the search for seed black holes at cosmic dawn, as well as the study of stellar-origin black holes and of their formation channels via multi-band and multi-messenger observations. We conclude that an extension to 6 yr of mission operations is recommended.

preprint2021arXiv

Gas-driven inspiral of binaries in thin accretion disks

Numerical studies of gas accretion onto supermassive black hole binaries (SMBHBs) have generally been limited to conditions where the circumbinary disk (CBD) is 10-100 times thicker than expected for disks in active galactic nuclei (AGN). This discrepancy arises from technical limitations, and also from publication bias toward replicating fiducial numerical models. Here we present the first systematic study of how the binary's orbital evolution varies with disk scale height. We report three key results: (1) Binary orbital evolution switches from outspiralling for warm disks (aspect ratio ~0.1), to inspiralling for more realistic cooler, thinner disks at a critical aspect ratio ~0.04, corresponding to orbital Mach number ~25. (2) The net torque on the binary arises from a competition between positive torque from gas orbiting close to the black holes, and negative torque from the inner edge of the CBD, which is denser for thinner disks. This leads to increasingly negative net torques on the binary for increasingly thin disks. (3) The accretion rate is modestly suppressed with increasing Mach number. We discuss how our results may influence modeling of the nano-Hz gravitational wave background, as well as estimates of the LISA merger event rate.

preprint2021arXiv

Multi-Messenger Constraints on Magnetic Fields in Merging Black Hole-Neutron Star Binaries

The LIGO-Virgo-KAGRA Collaboration recently detected gravitational waves (GWs) from the merger of black-hole-neutron-star (BHNS) binary systems GW200105 and GW200115. No coincident electromagnetic (EM) counterparts were detected. While the mass ratio and BH spin in both systems were not sufficient to tidally disrupt the NS outside of the BH event horizon, other, magnetospheric mechanisms for EM emission exist in this regime and depend sensitively on the NS magnetic field strength. Combining GW measurements with EM flux upper limits, we place upper limits on the NS surface magnetic field strength above which magnetospheric emission models would have generated an observable EM counterpart. We consider fireball models powered by the black-hole battery mechanism, where energy is output in gamma-rays over $\lesssim1$~second. Consistency with no detection by Fermi-GBM or INTEGRAL SPI-ACS constrains the NS surface magnetic field to $\lesssim10^{15}$~G. Hence, joint GW detection and EM upper limits rule out the theoretical possibility that the NSs in GW200105 and GW200115, and the putative NS in GW190814, retain $\gtrsim10^{15}$~G dipolar magnetic fields until merger. They also rule out formation scenarios where strongly magnetized magnetars quickly merge with BHs. We alternatively rule out operation of the BH-battery powered fireball mechanism in these systems. This is the first multi-messenger constraint on NS magnetic fields in BHNS systems and a novel approach to probe fields at this point in NS evolution. This demonstrates the constraining power that multi-messenger analyses of BHNS mergers have on BHNS formation scenarios, the magnetic-field evolution in NSs, and the physics of BHNS magnetospheric interactions.

preprint2020arXiv

Circumbinary Disks: Accretion and Torque as a Function of Mass Ratio and Disk Viscosity

Using numerical hydrodynamics calculations and a novel method for densely sampling parameter space, we measure the accretion and torque on a binary system from a circumbinary disk. In agreement with previous studies, we find that the net torque on the binary is positive for mass ratios close to unity, and that accretion always drives the binary towards equal mass. Accretion variability depends sensitively on the numerical sink prescription, but the torque and relative accretion onto each component do not depend on the sink timescale. Positive torque and highly variable accretion occurs only for mass ratios greater than around $0.05$. This means that for mass ratios below $0.05$, the binary would migrate inward until the secondary accreted sufficient mass, after which it would execute a U-turn and migrate outward. We explore a range of viscosities, from $α= 0.03$ to $α= 0.15$, and find that this outward torque is proportional to the viscous torque, simply proportional to viscosity in this range. Dependence of accretion and torque on mass ratio is explored in detail, densely sampling mass ratios between $0.01$ and unity. For mass ratio $q > 0.6$, accretion variability is found to exhibit a distinct sawtooth pattern, typically with a five-orbit cycle that provides a "smoking gun" prediction for variable quasars observed over long periods, as a potential means to confirm the presence of a binary.

preprint2020arXiv

Eccentric Black Hole Mergers in Active Galactic Nuclei

The astrophysical origin of gravitational wave (GW) transients is a timely open question in the wake of discoveries by LIGO/Virgo. In active galactic nuclei (AGNs), binaries form and evolve efficiently by interaction with a dense population of stars and the gaseous AGN disk. Previous studies have shown that stellar-mass black hole (BH) mergers in such environments can explain the merger rate and the number of suspected hierarchical mergers observed by LIGO/Virgo. The binary eccentricity distribution can provide further information to distinguish between astrophysical models. Here we derive the eccentricity distribution of BH mergers in AGN disks. We find that eccentricity is mainly due to binary-single (BS) interactions, which lead to most BH mergers in AGN disks having a significant eccentricity at $0.01\,\mathrm{Hz}$, detectable by LISA. If BS interactions occur in isotropic-3D directions, then $8$--$30\%$ of the mergers in AGN disks will have eccentricities at $10\,\mathrm{Hz}$ above $e_{10\,\rm Hz}\gtrsim 0.03$, detectable by LIGO/Virgo/KAGRA, while $5$--$17\%$ of mergers have $e_{10\,\rm Hz}\geq 0.3$. On the other hand, if BS interactions are confined to the AGN-disk plane due to torques from the disk, with 1-20 intermediate binary states during each interaction, or if BHs can migrate to $\lesssim10^{-3}\,\mathrm{pc}$ from the central supermassive black hole, then $10$--$70\%$ of the mergers will be highly eccentric ($e_{10\,\rm Hz} \geq 0.3$), consistent with the possible high eccentricity in GW190521.

preprint2020arXiv

Formation and Evolution of Compact Object Binaries in AGN Disks

The astrophysical origin of gravitational wave (GW) events discovered by LIGO/VIRGO remains an outstanding puzzle. In active galactic nuclei (AGN), compact-object binaries form, evolve, and interact with a dense star cluster and a gas disk. An important question is whether and how binaries merge in these environments. To address this question, we have performed one-dimensional $N$-body simulations combined with a semi-analytical model which includes the formation, disruption, and evolution of binaries self-consistently. We point out that binaries can form in single-single interactions by the dissipation of kinetic energy in a gaseous medium. This ``gas capture'' binary formation channel contributes up to $97\,\%$ of gas-driven mergers and leads to a high merger rate in AGN disks even without pre-existing binaries. We find the merger rate to be in the range $\sim 0.02-60\,\mathrm{Gpc^{-3}yr^{-1}}$. The results are insensitive to the assumptions on gaseous hardening processes: we find that once they are formed, binaries merge efficiently via binary-single interactions even if these gaseous processes are neglected. We find that the average number of mergers per BH is $0.4$, and the probability for repeated mergers in 30 Myr is $\sim 0.21-0.45$. High BH masses due to repeated mergers, high eccentricities, and a significant Doppler drift of GWs are promising signatures which distinguish this merger channel from others. Furthermore, we find that gas-capture binaries reproduce the distribution of LMXBs in the Galactic center, including an outer cutoff at $\sim1$ pc due to the competition between migration and hardening by gas torques.

preprint2020arXiv

Getting Ready for LISA: The Data, Support and Preparation Needed to Maximize US Participation in Space-Based Gravitational Wave Science

The NASA LISA Study Team was tasked to study how NASA might support US scientists to participate and maximize the science return from the Laser Interferometer Space Antenna (LISA) mission. LISA is gravitational wave observatory led by ESA with NASA as a junior partner, and is scheduled to launch in 2034. Among our findings: LISA science productivity is greatly enhanced by a full-featured US science center and an open access data model. As other major missions have demonstrated, a science center acts as both a locus and an amplifier of research innovation, data analysis, user support, user training and user interaction. In its most basic function, a US Science Center could facilitate entry into LISA science by hosting a Data Processing Center and a portal for the US community to access LISA data products. However, an enhanced LISA Science Center could: support one of the parallel independent processing pipelines required for data product validation; stimulate the high level of research on data analysis that LISA demands; support users unfamiliar with a novel observatory; facilitate astrophysics and fundamental research; provide an interface into the subtleties of the instrument to validate extraordinary discoveries; train new users; and expand the research community through guest investigator, postdoc and student programs. Establishing a US LISA Science Center well before launch can have a beneficial impact on the participation of the broader astronomical community by providing training, hosting topical workshops, disseminating mock catalogs, software pipelines, and documentation. Past experience indicates that successful science centers are established several years before launch; this early adoption model may be especially relevant for a pioneering mission like LISA.

preprint2020arXiv

Making a supermassive star by stellar bombardment

Approximately two hundred supermassive black holes (SMBHs) have been discovered within the first $\sim$Gyr after the Big Bang. One pathway for the formation of SMBHs is through the collapse of supermassive stars (SMSs). A possible obstacle to this scenario is that the collapsing gas fragments and forms a cluster of main-sequence stars. Here we raise the possibility that stellar collisions may be sufficiently frequent and energetic to inhibit the contraction of the massive protostar, avoiding strong UV radiation driven outflows, and allowing it to continue growing into an SMS. We investigate this scenario with semianalytic models incorporating star formation, gas accretion, dynamical friction from stars and gas, stellar collisions, and gas ejection. We find that when the collapsing gas fragments at a density of $\lesssim 3\times 10^{10}\,\mathrm{cm^{-3}}$, the central protostar contracts due to infrequent stellar mergers, and in turn photoevaporates the remaining collapsing gas, resulting in the formation of a $\lesssim 10^4~{\rm M_\odot}$ object. On the other hand, when the collapsing gas fragments at higher densities (expected for a metal-poor cloud with $Z\lesssim10^{-5}\,{\rm Z_\odot}$ with suppressed ${\rm H_2}$ abundance) the central protostar avoids contraction and keeps growing via frequent stellar mergers, reaching masses as high as $\sim 10^5-10^6\,{\rm M_\odot}$. We conclude that frequent stellar mergers represent a possible pathway to form massive BHs in the early universe.

preprint2020arXiv

Mass-gap Mergers in Active Galactic Nuclei

The recently discovered gravitational wave sources GW190521 and GW190814 have shown evidence of BH mergers with masses and spins that could be outside of the range expected from isolated stellar evolution. These merging objects could have undergone previous mergers. Such hierarchical mergers are predicted to be frequent in active galactic nuclei (AGN) disks, where binaries form and evolve efficiently by dynamical interactions and gaseous dissipation. Here we compare the properties of these observed events to the theoretical models of mergers in AGN disks, which are obtained by performing one-dimensional $N$-body simulations combined with semi-analytical prescriptions. The high BH masses in GW190521 are consistent with mergers of high-generation (high-g) BHs where the initial progenitor stars had high metallicity, 2g BHs if the original progenitors were metal-poor, or 1g BHs that had gained mass via super-Eddington accretion. Other measured properties related to spin parameters in GW190521 are also consistent with mergers in AGN disks. Furthermore, mergers in the lower mass gap or those with low mass ratio as found in GW190814 and GW190412 are also reproduced by mergers of 2g-1g or 1g-1g objects with significant accretion in AGN disks. Finally, due to gas accretion, the massive neutron star merger reported in GW190425 can be produced in an AGN disk.

preprint2020arXiv

Self-consistent Semi-analytic Modeling of Feedback During Primordial Star Formation and Reionization

We present a new semi-analytic model of the formation of the first stars. Our method takes dark matter halo merger trees (including 3-dimensional spatial information) from cosmological N-body simulations as input and applies analytic prescriptions to compute both the Population III (Pop III) and metal-enriched star formation histories. We have developed a novel method to accurately compute the major feedback processes affecting Pop III star formation: H$_2$ photodissociation from Lyman-Werner (LW) radiation, suppression of star formation due to inhomogeneous reionization, and metal enrichment via supernovae winds. Our method utilizes a grid-based approach relying on fast Fourier transforms (FFTs) to rapidly track the LW intensity, ionization fraction, and metallicity in 3-dimensions throughout the simulation box. We present simulations for a wide range of astrophysical model parameters from $z\approx 30-6$. Initially long-range LW feedback and local metal enrichment and reionization feedback dominate. However, for $z \lesssim 15$ we find that the star formation rate density (SFRD) of Pop III stars is impacted by the combination of external metal enrichment (metals from one halo polluting other pristine halos) and inhomogeneous reionization. We find that the interplay of these processes is particularly important for the Pop III SFRD at $z \lesssim 10$. Reionization feedback delays star formation long enough for metal bubbles to reach halos that would otherwise form Pop III stars. Including these effects can lead to more than an order of magnitude decrease in the Pop III SFRD at $z=6$ compared to LW feedback alone.

preprint2020arXiv

Spikey: Self-Lensing Flares from Eccentric SMBH Binaries

We examine the light curves of two quasars, motivated by recent suggestions that a supermassive black hole binary (SMBHB) can exhibit sharp lensing spikes. We model the variability of each light curve as due to a combination of two relativistic effects: the orbital relativistic Doppler boost and gravitational binary self-lensing. In order to model each system we extend previous Doppler plus self-lensing models to include eccentricity. The first quasar is identified in optical data as a binary candidate with a 20-yr period (Ark 120), and shows a prominent spike. For this source, we rule out the lensing hypothesis and disfavor the Doppler-boost hypothesis due to discrepancies in the measured vs. recovered values of the binary mass and optical spectral slope. The second source, which we nickname Spikey, is the rare case of an active galactic nucleus (AGN) identified in Kepler's high-quality, high-cadence photometric data. For this source, we find a model, consisting of a combination of a Doppler modulation and a narrow symmetric lensing spike, consistent with an eccentric SMBHB with a total mass of approximately 30 million solar masses, rest-frame orbital period T=418 days, eccentricity e=0.5, and seen at an inclination of 8 degrees from edge-on. This interpretation can be tested by monitoring Spikey for periodic behavior and recurring flares in the next few years. In preparation for such monitoring we present the first X-ray observations of this object taken by the Neil Gehrels Swift observatory.

preprint2020arXiv

Spin Evolution of Stellar-mass Black Hole Binaries in Active Galactic Nuclei

The astrophysical origin of gravitational wave (GW) events is one of the most timely problems in the wake of the LIGO/Virgo discoveries. In active galactic nuclei (AGN), binaries form and evolve efficiently by dynamical interactions and gaseous dissipation. Previous studies have suggested that binary black hole (BBH) mergers in AGN disks can contribute significantly to BBH mergers observed by GW interferometers. Here we examine the distribution of the effective spin parameter $χ_\mathrm{eff}$ of this GW source population. We extend our semi-analytical model of binary formation and evolution in AGN disks by following the evolution of the binary orbital angular momenta and black hole (BH) spins. BH spins change due to gas accretion and BH mergers, while the binary orbital angular momenta evolve due to gas accretion and binary-single interactions. We find that the distribution of $χ_\mathrm{eff}$ predicted by our AGN model is similar to the distribution observed during LIGO/Virgo O1 and O2. On the other hand, if radial migration of BHs is inefficient, $χ_\mathrm{eff}$ is skewed toward higher values compared with the observed distribution, because of the paucity of scattering events that would randomize spin directions relative to the orbital plane. We suggest that high binary masses and the positive correlation between binary mass and the standard deviation of $χ_\mathrm{eff}$ for chirp masses up to $\approx 20$ $\mathrm{M}_\odot$, can be possible signatures for mergers originating in AGN disks. Finally, hierarchical mergers in AGN disks naturally produce properties of the recent GW event GW190412, including a low mass ratio, a high primary BH spin, and a significant spin component in the orbital plane.

preprint2020arXiv

Testing the relativistic Doppler boost hypothesis for the binary candidate quasar PG1302-102 with multi-band Swift data

The bright quasar PG1302-102 has been identified as a candidate supermassive black hole binary from its near-sinusoidal optical variability. While the significance of its optical periodicity has been debated due to the stochastic variability of quasars, its multi-wavelength variability in the ultraviolet (UV) and optical bands is consistent with relativistic Doppler boost caused by the orbital motion in a binary. However, this conclusion was based previously on sparse UV data which were not taken simultaneously with the optical data. Here we report simultaneous follow-up observations of PG1302-102 with the Ultraviolet Optical Telescope on the Neil Gehrels Swift Observatory in six optical + UV bands. The additional nine Swift observations produce light curves roughly consistent with the trend under the Doppler boost hypothesis, which predicts that UV variability should track the optical, but with a ~2.2 times higher amplitude. We perform a statistical analysis to quantitatively test this hypothesis. We find that the data are consistent with the Doppler boost hypothesis when we compare the the amplitudes in optical B-band and UV light curves. However, the ratio of UV to V-band variability is larger than expected and is consistent with the Doppler model, only if either the UV/optical spectral slopes vary, the stochastic variability makes a large contribution in the UV, or the sparse new optical data underestimate the true optical variability. We have evidence for the latter from comparison with the optical light curve from ASAS-SN. Additionally, the simultaneous analysis of all four bands strongly disfavors the Doppler boost model whenever Swift V-band is involved. Additional, simultaneous optical + UV observations tracing out another cycle of the 5.2-year proposed periodicity should lead to a definitive conclusion.

preprint2020arXiv

The Quest for Dual and Binary Supermassive Black Holes: A Multi-Messenger View

The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field.

preprint2016arXiv

A transition in circumbinary accretion discs at a binary mass ratio of 1:25

We study circumbinary accretion discs in the framework of the restricted three-body problem (R3Bp) and via numerically solving the height-integrated equations of viscous hydrodynamics. Varying the mass ratio of the binary, we find a pronounced change in the behaviour of the disc near mass ratio $q \equiv M_s/M_p \sim 0.04$. For mass ratios above $q=0.04$, solutions for the hydrodynamic flow transition from steady, to strongly-fluctuating; a narrow annular gap in the surface density around the secondary's orbit changes to a hollow central cavity; and a spatial symmetry is lost, resulting in a lopsided disc. This phase transition is coincident with the mass ratio above which stable orbits do not exist around the L4 and L5 equilibrium points of the R3B problem. Using the DISCO code, we find that for thin discs, for which a gap or cavity can remain open, the mass ratio of the transition is relatively insensitive to disc viscosity and pressure. The $q=0.04$ transition has relevance for the evolution of massive black hole binary+disc systems at the centers of galactic nuclei, as well as for young stellar binaries and possibly planets around brown dwarfs.

preprint2016arXiv

Assisted Inspirals of Stellar Mass Black Holes Embedded in AGN Disks: Solving the "Final AU Problem"

We explore the evolution of stellar mass black hole binaries (BHBs) which are formed in the self-gravitating disks of active galactic nuclei (AGN). Hardening due to three-body scattering and gaseous drag are effective mechanisms that reduce the semi-major axis of a BHB to radii where gravitational waves take over, on timescales shorter than the typical lifetime of the AGN disk. Taking observationally-motivated assumptions for the rate of star formation in AGN disks, we find a rate of disk-induced BHB mergers ($\mathcal{R} \sim 3~{\rm yr}^{-1}~{\rm Gpc}^{-3}$, but with large uncertainties) that is comparable with existing estimates of the field rate of BHB mergers, and the approximate BHB merger rate implied by the recent Advanced LIGO detection of GW150914. BHBs formed thorough this channel will frequently be associated with luminous AGN, which are relatively rare within the sky error regions of future gravitational wave detector arrays. This channel could also possess a (potentially transient) electromagnetic counterpart due to super-Eddington accretion onto the stellar mass black hole following the merger.

preprint2016arXiv

Formation of Massive Population III Galaxies through Photoionization Feedback: A Possible Explanation for CR7

We explore the formation of massive high-redshift Population III (Pop III) galaxies through photoionization feedback. We consider dark matter halos formed from progenitors that have undergone no star formation as a result of early reionization and photoevaporation caused by a nearby galaxy. Once such a halo reaches $\approx 10^9~M_\odot$, corresponding to the Jeans mass of the photoheated intergalactic medium (IGM) at $z\approx 7$, pristine gas is able to collapse into the halo, potentially producing a massive Pop III starburst. We suggest that this scenario may explain the recent observation of strong He II 1640~Å~line emission in CR7, which is consistent with $\sim 10^7~M_\odot$ of young Pop III stars. Such a large mass of Pop III stars is unlikely without the photoionization feedback scenario, because star formation is expected to inject metals into halos above the atomic cooling threshold ($\sim 10^8~M_\odot$ at $z \approx 7$). We use merger trees to analytically estimate the abundance of observable Pop III galaxies formed through this channel, and find a number density of $\approx 10^{-7}~{\rm Mpc^{-3}}$ at $z=6.6$ (the redshift of CR7). This is approximately a factor of ten lower than the density of Ly$α$ emitters as bright as CR7.

preprint2016arXiv

Gravitational wave background from Population III binary black holes consistent with cosmic reionization

The recent discovery of the gravitational wave source GW150914 has revealed a coalescing binary black hole (BBH) with masses of $\sim 30~M_\odot$. Previous proposals for the origin of such a massive binary include Population III (PopIII) stars. PopIII stars are efficient producers of BBHs and of a gravitational wave background (GWB) in the $10-100$ Hz band, and also of ionizing radiation in the early Universe. We quantify the relation between the amplitude of the GWB ($Ω_{\rm gw}$) and the electron scattering optical depth ($τ_{\rm e}$), produced by PopIII stars, assuming that $f_{\rm esc}\approx 10\%$ of their ionizing radiation escapes into the intergalactic medium. We find that PopIII stars would produce a GWB that is detectable by the future O5 LIGO/Virgo if $τ_{\rm e} \gtrsim 0.07$, consistent with the recent Planck measurement of $τ_e=0.055 \pm 0.09$. Moreover, the spectral index of the background from PopIII BBHs becomes as small as ${\rm d}\ln Ω_{\rm gw}/{\rm d}\ln f\lesssim 0.3$ at $f \gtrsim 30$ Hz, which is significantly flatter than the value $\sim 2/3$ generically produced by lower-redshift and less-massive BBHs. A detection of the unique flattening at such low frequencies by the O5 LIGO/Virgo will indicate the existence of a high-chirp mass, high-redshift BBH population, which is consistent with the PopIII origin. A precise characterization of the spectral shape near $30-50$ Hz by the Einstein Telescope could also constrain the PopIII initial mass function and star formation rate.

preprint2016arXiv

Hyper-Eddington accretion flows onto massive black holes

We study very-high rate spherically symmetric accretion flows onto a massive black hole (BH; 10^2 < M_BH < 10^6 Msun) embedded in a dense gas cloud with a low abundance of metals, performing one-dimensional hydrodynamical simulations which include multi-frequency radiation transfer and non-equilibrium primordial chemistry. We find that rapid gas supply from the Bondi radius at a hyper-Eddington rate can occur without being impeded by radiation feedback when (n/10^5 cm^-3) > (M_BH/10^4Msun)^{-1}(T/10^4 K)^{3/2}, where n and T are the density and temperature of ambient gas outside of the Bondi radius. The resulting accretion rate in this regime is steady, and larger than 3000 times the Eddington rate. At lower Bondi rates, the accretion is episodic due to radiative feedback and the average rate is limited below the Eddington rate. For the hyper-Eddington case, the steady solution consists of two parts: a radiation-dominated central core, where photon trapping due to electron scattering is important, and an accreting envelope which follows a Bondi profile with T~8000 K. When the emergent luminosity is limited below the Eddington luminosity because of photon trapping, radiation from the central region does not affect the gas dynamics at larger scales. We apply our result to the rapid formation of massive BHs in protogalaxies with a virial temperature of T_vir> 10^4 K. Once a seed BH forms at the center of the galaxy, it can grow up to a maximum ~10^5 (T_vir/10^4 K) Msun via gas accretion independent of the initial BH mass. Finally, we discuss possible observational signatures of rapidly accreting BHs with/without allowance for dust. We suggest that these systems could explain Lya emitters without X-rays and luminous infrared sources with hot dust emission, respectively.

preprint2016arXiv

Intermediate-mass black holes from Population III remnants in the first galactic nuclei

We report the formation of intermediate-mass black holes (IMBHs) in suites of numerical $N$-body simulations of Population III remnant black holes (BHs) embedded in gas-rich protogalaxies at redshifts $z\gtrsim10$. We model the effects of gas drag on the BHs' orbits, and allow BHs to grow via gas accretion, including a mode of hyper-Eddington accretion in which photon trapping and rapid gas inflow suppress any negative radiative feedback. Most initial BH configurations lead to the formation of one (but never more than one) IMBH in the center of the protogalaxy, reaching a mass of $10^{3-5}\mathrm{M}_{\odot}$ through hyper-Eddington growth. Our results suggest a viable pathway to forming the earliest massive BHs in the centers of early galaxies. We also find that the nuclear IMBH typically captures a stellar-mass BH companion, making these systems observable in gravitational waves as extreme mass-ratio inspirals (EMRIs) with \textit{eLISA}.

preprint2016arXiv

Is there a maximum mass for black holes in galactic nuclei?

The largest observed supermassive black holes (SMBHs) have a mass of M_BH ~ 10^{10} M_sun, nearly independent of redshift, from the local (z~0) to the early (z>6) Universe. We suggest that the growth of SMBHs above a few 10^{10} M_sun is prevented by small-scale accretion physics, independent of the properties of their host galaxies or of cosmology. Growing more massive BHs requires a gas supply rate from galactic scales onto a nuclear region as high as >10^3 M_sun/yr. At such a high accretion rate, most of the gas converts to stars at large radii (~10-100 pc), well before reaching the BH. We adopt a simple model (Thompson et al. 2005) for a star-forming accretion disk, and find that the accretion rate in the sub-pc nuclear region is reduced to the smaller value of at most a few M_sun/yr. This prevents SMBHs from growing above ~10^{11} M_sun in the age of the Universe. Furthermore, once a SMBH reaches a sufficiently high mass, this rate falls below the critical value at which the accretion flow becomes advection dominated. Once this transition occurs, BH feeding can be suppressed by strong outflows and jets from hot gas near the BH. We find that the maximum SMBH mass, given by this transition, is between M_{BH,max} ~ (1-6) * 10^{10} M_sun, depending primarily on the efficiency of angular momentum transfer inside the galactic disk, and not on other properties of the host galaxy.

preprint2016arXiv

The Origin of Weak Lensing Convergence Peaks

Weak lensing convergence peaks are a promising tool to probe nonlinear structure evolution at late times, providing additional cosmological information beyond second-order statistics. Previous theoretical and observational studies have shown that the cosmological constraints on $Ω_m$ and $σ_8$ are improved by a factor of up to ~ 2 when peak counts and second-order statistics are combined, compared to using the latter alone. We study the origin of lensing peaks using observational data from the 154 deg$^2$ Canada-France-Hawaii Telescope Lensing Survey. We found that while high peaks (with height $κ$ >3.5 $σ_κ$, where $σ_κ$ is the r.m.s. of the convergence $κ$) are typically due to one single massive halo of ~$10^{15}M_\odot$, low peaks ($κ$ <~ $σ_κ$) are associated with constellations of 2-8 smaller halos (<~$10^{13}M_\odot$). In addition, halos responsible for forming low peaks are found to be significantly offset from the line-of-sight towards the peak center (impact parameter >~ their virial radii), compared with ~0.25 virial radii for halos linked with high peaks, hinting that low peaks are more immune to baryonic processes whose impact is confined to the inner regions of the dark matter halos. Our findings are in good agreement with results from the simulation work by Yang el al. (2011).

preprint2015arXiv

A Statistical Method to Search for Recoiling Supermassive Black Holes in Active Galactic Nuclei

We propose an observational test for gravitationally recoiling supermassive black holes (BHs) in active galactic nuclei, based on a correlation between the velocities of BHs relative to their host galaxies, |Δv|, and their obscuring dust column densities, Σ_{dust} (both measured along the line of sight). We use toy models for the distribution of recoil velocities, BH trajectories, and the geometry of obscuring dust tori in galactic centres, to simulate 2.5x10^5 random observations of recoiling quasars. BHs with recoil velocities comparable to the escape velocity from the galactic centre remain bound to the nucleus, and do not fully settle back to the centre of the torus due to dynamical friction in a typical quasar lifetime. We find that |Δv| and Σ_ {dust} for these BHs are positively correlated. For obscured (Σ_{dust}>0) and for partially obscured (0<Σ_{dust}<~2.3 g/m^2) quasars with |Δv|>=45 km/s, the sample correlation coefficient between log10(|Δv|) and Σ_{dust} is r_{45} = 0.28+/-0.02 and r_{45} = 0.13+/-0.02, respectively. Allowing for random +/-100 km/s errors in |Δv| unrelated to the recoil dilutes the correlation for the partially obscured quasars to r_{45} = 0.026+/-0.004 measured between |Δv| and Σ_{dust}. A random sample of >~3,500 obscured quasars with |Δv|>=45 km/s would allow rejection of the no-correlation hypothesis with 3 sigma significance 95% of the time. Finally, we find that the fraction of obscured quasars, F_{obs}(|Δv|), decreases with |Δv| from F_{obs}(<10 km/s)>~0.8 to F_{obs}(>10^3 km/s)<~0.4. This predicted trend can be compared to the observed fraction of type II quasars, and can further test combinations of recoil, trajectory, and dust torus models.

preprint2015arXiv

Cosmology Constraints from the Weak Lensing Peak Counts and the Power Spectrum in CFHTLenS

Lensing peaks have been proposed as a useful statistic, containing cosmological information from non-Gaussianities that is inaccessible from traditional two-point statistics such as the power spectrum or two-point correlation functions. Here we examine constraints on cosmological parameters from weak lensing peak counts, using the publicly available data from the 154 deg$^2$ CFHTLenS survey. We utilize a new suite of ray-tracing N-body simulations on a grid of 91 cosmological models, covering broad ranges of the three parameters $Ω_m$, $σ_8$, and $w$, and replicating the Galaxy sky positions, redshifts, and shape noise in the CFHTLenS observations. We then build an emulator that interpolates the power spectrum and the peak counts to an accuracy of $\leq 5\%$, and compute the likelihood in the three-dimensional parameter space ($Ω_m$, $σ_8$, $w$) from both observables. We find that constraints from peak counts are comparable to those from the power spectrum, and somewhat tighter when different smoothing scales are combined. Neither observable can constrain $w$ without external data. When the power spectrum and peak counts are combined, the area of the error "banana'' in the ($Ω_m$, $σ_8$) plane reduces by a factor of $\approx2$, compared to using the power spectrum alone. For a flat $Λ$ cold dark matter model, combining both statistics, we obtain the constraint $σ_8(Ω_m/0.27)^{0.63}=0.85\substack{+0.03 \\ -0.03}$.

preprint2015arXiv

Direct collapse black hole formation from synchronized pairs of atomic cooling halos

High-redshift quasar observations imply that supermassive black holes (SMBHs) larger than $\sim 10^9 ~ M_\odot$ formed before $z=6$. That such large SMBHs formed so early in the Universe remains an open theoretical problem. One possibility is that gas in atomic cooling halos exposed to strong Lyman-Werner (LW) radiation forms $10^4-10^6 ~ M_\odot$ supermassive stars which quickly collapse into black holes. We propose a scenario for direct collapse black hole (DCBH) formation based on synchronized pairs of pristine atomic cooling halos. We consider halos at very small separation with one halo being a subhalo of the other. The first halo to surpass the atomic cooling threshold forms stars. Soon after these stars are formed, the other halo reaches the cooling threshold and due to its small distance from the newly formed galaxy, is exposed to the critical LW intensity required to form a DCBH. The main advantage of this scenario is that synchronization can potentially prevent photoevaporation and metal pollution in DCBH-forming halos. Since the halos reach the atomic cooling threshold at nearly the same time, the DCBH-forming halo is only exposed to ionizing radiation for a brief period. Tight synchronization could allow the DCBH to form before stars in the nearby galaxy reach the end of their lives and generate supernovae winds. We use N-body simulations to estimate the abundance of DCBHs formed in this way. The largest source of uncertainty in our estimate is the initial mass function (IMF) of metal free stars formed in atomic cooling halos. We find that even for tight synchronization, the density of DCBHs formed in this scenario could explain the SMBHs implied by $z=6$ quasar observations. Metal pollution and photoevaporation could potentially reduce the abundance of DCBHs below that required to explain the observations in other models that rely on a high LW flux.

preprint2015arXiv

Emulating the CFHTLenS Weak Lensing data: Cosmological Constraints from moments and Minkowski functionals

Weak gravitational lensing is a powerful cosmological probe, with non--Gaussian features potentially containing the majority of the information. We examine constraints on the parameter triplet $(Ω_m,w,σ_8)$ from non-Gaussian features of the weak lensing convergence field, including a set of moments (up to $4^{\rm th}$ order) and Minkowski functionals, using publicly available data from the 154deg$^2$ CFHTLenS survey. We utilize a suite of ray--tracing N-body simulations spanning 91 points in $(Ω_m,w,σ_8)$ parameter space, replicating the galaxy sky positions, redshifts and shape noise in the CFHTLenS catalogs. We then build an emulator that interpolates the simulated descriptors as a function of $(Ω_m,w,σ_8)$, and use it to compute the likelihood function and parameter constraints. We employ a principal component analysis to reduce dimensionality and to help stabilize the constraints with respect to the number of bins used to construct each statistic. Using the full set of statistics, we find $Σ_8\equivσ_8(Ω_m/0.27)^{0.55}=0.75\pm0.04$ (68% C.L.), in agreement with previous values. We find that constraints on the $(Ω_m,σ_8)$ doublet from the Minkowski functionals suffer a strong bias. However, high-order moments break the $(Ω_m,σ_8)$ degeneracy and provide a tight constraint on these parameters with no apparent bias. The main contribution comes from quartic moments of derivatives.

preprint2015arXiv

High-Redshift Star Formation in a Time-Dependent Lyman-Werner Background

The first generation of stars produces a background of Lyman-Werner (LW) radiation which can photo-dissociate molecular hydrogen, increasing the mass of dark matter halos required to host star formation. Previous studies have determined the critical mass required for efficient molecular cooling with a constant LW background. However, the true background is expected to increase rapidly at early times. Neglecting this evolution could underestimate star formation in small halos that may have started to cool in the past when the LW intensity was much lower. Background evolution is a large source of uncertainty in pre-reionization predictions of the cosmological 21cm signal, which can be observed with future radio telescopes. To address this, we perform zero-dimentional one-zone calculations that follow the density, chemical abundances, and temperature of gas in the central regions of dark matter halos, including hierarchical growth and an evolving LW background. We begin by studying the physics of halos subjected to a background that increases exponentially with redshift. We find that when the intensity increases more slowly than $J_{\rm LW}(z) \propto 10^{-z/5}$, cooling in the past is a relatively small effect. We then self-consistently compute the cosmological LW background over $z=15-50$ and find that cooling in the past due to an evolving background has a modest impact. Finally, we compare these results to three-dimensional hydrodynamical cosmological simulations with varying LW histories. While only a small number of halos were simulated, the results are consistent with our one-zone calculations.

preprint2015arXiv

Limits on Population III star formation in minihaloes implied by Planck

Recently, Planck measured a value of the cosmic microwave background (CMB) optical depth due to electron scattering of $τ=0.066 \pm 0.016$. Here we show that this low value leaves essentially no room for an early partial reionisation of the intergalactic medium (IGM) by high-redshift Population III (Pop III) stars, expected to have formed in low-mass minihaloes. We perform semi-analytic calculations of reionisation which include the contribution from Pop II stars in atomic cooling haloes, calibrated with high-redshift galaxy observations, and Pop III stars in minihaloes with feedback due to Lyman-Werner (LW) radiation and metal enrichment. We find that without LW feedback or prompt metal enrichment (and assuming a minihalo escape fraction of 0.5) the Pop III star formation efficiency cannot exceed $\sim{\rm a~few}\times 10^{-4}$, without violating the constraints set by Planck data. This excludes massive Pop III star formation in typical $10^6 M_\odot$ minihaloes. Including LW feedback and metal enrichment alleviates this tension, allowing large Pop III stars to form early on before they are quenched by feedback. We find that the total density of Pop III stars formed across cosmic time is $\lesssim 10^{4-5}~M_\odot~{\rm Mpc^{-3}}$ and does not depend strongly on the feedback prescription adopted. Additionally, we perform a simple estimate of the possible impact on reionisation of X-rays produced by accretion onto black hole remnants of Pop III stars. We find that unless the accretion duty cycle is very low ($\lesssim 0.01$), this could lead to an optical depth inconsistent with Planck.

preprint2015arXiv

Looking for Population III stars with He II line intensity mapping

Constraining the properties of Population III (Pop III) stars will be very challenging because they reside in small galaxies at high redshift which will be difficult to directly detect. In this paper, we suggest that intensity mapping may be a promising method to study Pop III stars. Intensity mapping is a technique proposed to measure large-scale fluctuations of galaxy line emission in three dimensions without resolving individual sources. This technique is well suited for observing many faint galaxies because it can measure their cumulative emission even if they cannot be directly detected. We focus on intensity mapping of He II recombination lines. These lines are much stronger in Pop III stars than Pop II stars because the harder spectra of Pop III stars are expected to produce many He II ionizing photons. Measuring the He II 1640 Å intensity mapping signal, along with the signals from other lines such as Ly$α$, H$α$, and metal lines, could give constraints on the initial mass function (IMF) and star formation rate density (SFRD) of Pop III stars as a function of redshift. To demonstrate the feasibility of these observations, we estimate the strength of the Pop III He II 1640 Å intensity mapping signal from $z=10-20$. We show that at $z\approx10$, the signal could be measured accurately by two different hypothetical future instruments, one which cross-correlates He II 1640 Å with CO(1-0) line emission from galaxies and the other with 21 cm emission from the intergalactic medium (IGM).

preprint2015arXiv

The impact of spurious shear on cosmological parameter estimates from weak lensing observables

Residual errors in shear measurements, after corrections for instrument systematics and atmospheric effects, can impact cosmological parameters derived from weak lensing observations. Here we combine convergence maps from our suite of ray-tracing simulations with random realizations of spurious shear. This allows us to quantify the errors and biases of the triplet $(Ω_m,w,σ_8)$ derived from the power spectrum (PS), as well as from three different sets of non-Gaussian statistics of the lensing convergence field: Minkowski functionals (MF), low--order moments (LM), and peak counts (PK). Our main results are: (i) We find an order of magnitude smaller biases from the PS than in previous work. (ii) The PS and LM yield biases much smaller than the morphological statistics (MF, PK). (iii) For strictly Gaussian spurious shear with integrated amplitude as low as its current estimate of $σ^2_{sys}\approx 10^{-7}$, biases from the PS and LM would be unimportant even for a survey with the statistical power of LSST. However, we find that for surveys larger than $\approx 100$ deg$^2$, non-Gaussianity in the noise (not included in our analysis) will likely be important and must be quantified to assess the biases. (iv) The morphological statistics (MF,PK) introduce important biases even for Gaussian noise, which must be corrected in large surveys. The biases are in different directions in $(Ω_m,w,σ_8)$ parameter space, allowing self-calibration by combining multiple statistics. Our results warrant follow-up studies with more extensive lensing simulations and more accurate spurious shear estimates.

preprint2014arXiv

A No-Go Theorem for Direct Collapse Black Holes Without a Strong Ultraviolet Background

Explaining the existence of supermassive black holes (SMBHs) larger than $\sim 10^9 M_\odot$ at redshifts $z >\sim 6$ remains an open theoretical question. One possibility is that gas collapsing rapidly in pristine atomic cooling halos ($T_{\rm vir} >\sim 10^4 \rm{K}$) produces $10^4-10^6 M_\odot$ black holes. Previous studies have shown that the formation of such a black hole requires a strong UV background to prevent molecular hydrogen cooling and gas fragmentation. Recently it has been proposed that a high UV background may not be required for halos that accrete material extremely rapidly or for halos where gas cooling is delayed due to a high baryon-dark matter streaming velocity. In this work, we point out that building up a halo with $T_{\rm vir} >\sim 10^4 \rm{K}$ before molecular cooling becomes efficient is not sufficient for forming a direct collapse black hole (DCBH). Though molecular hydrogen formation may be delayed, it will eventually form at high densities leading to efficient cooling and fragmentation. The only obvious way that molecular cooling could be avoided in the absence of strong UV radiation, is for gas to reach high enough density to cause collisional dissociation of molecular hydrogen ($\sim 10^4 ~ {\rm cm}^{-3}$) before cooling occurs. However, we argue that the minimum core entropy, set by the entropy of the intergalactic medium (IGM) when it decouples from the CMB, prevents this from occurring for realistic halo masses. This is confirmed by hydrodynamical cosmological simulations without radiative cooling. We explain the maximum density versus halo mass in these simulations with simple entropy arguments. The low densities found suggest that DCBH formation indeed requires a strong UV background.

preprint2014arXiv

Binary Black Hole Accretion During Inspiral and Merger

We present the results of 2D, moving mesh, viscous hydrodynamical simulations of accretion onto merging supermassive black hole (SMBH) binaries. We include viscous heating, shock heating, and radiative cooling, and simulate the transition from the "pre-decoupling" epoch, where the inspiral timescale is longer than the viscous timescale, to the "post-decoupling" epoch, where the inspiral timescale is shorter than the viscous timescale. We find that there is no abrupt halt to the accretion at decoupling, but rather the accretion shows a slow decay, with significant accretion well after the expected decoupling. Moreover, we find that the luminosity in X-rays is significantly higher prior to the merger, as orbital energy from the SMBH binary is converted to heat via strong shocks inside the cavity, and radiated away. Following the merger, the cavity refills viscously and the accretion rate relaxes to the Shakura-Sunyaev value, while the X-ray luminosity drops as the shocks quickly dissipate.

preprint2014arXiv

Characteristic Signatures in the Thermal Emission from Accreting Binary Black Holes

We present the results of a calculation of the thermal spectrum from a 2D, moving mesh, high-accuracy, viscous hydrodynamical simulation of an accreting supermassive black hole binary. We include viscous heating, shock heating, and radiative cooling, evolving for longer than a viscous time so that we reach a quasi-steady accretion state. In agreement with previous work, we find that gas is efficiently stripped from the inner edge of the circumbinary disk and enters the cavity along accretion streams, which feed persistent "mini-disks" surrounding each black hole. We also find that emission from the shock-heated mini-disks and accretion streams prevents any deficit in high-energy emission that may be expected inside the circumbinary cavity, and instead leads to a characteristic brightening of the spectrum beginning in soft X-rays.

preprint2014arXiv

Close companions to two high-redshift quasars

We report the serendipitous discoveries of companion galaxies to two high-redshift quasars. SDSS J025617.7+001904 is a z=4.79 quasar included in our recent survey of faint quasars in the SDSS Stripe 82 region. The initial MMT slit spectroscopy shows excess Lyman alpha emission extending well beyond the quasar's light profile. Further imaging and spectroscopy with LBT/MODS1 confirms the presence of a bright galaxy (i_AB = 23.6) located 2arcsec (12 kpc projected) from the quasar with strong Lyman alpha emission (EW_0 ~ 100Ang) at the redshift of the quasar, as well as faint continuum. The second quasar, CFHQS J005006.6+344522 (z=6.25), is included in our recent HST SNAP survey of z~6 quasars searching for evidence of gravitational lensing. Deep imaging with ACS and WFC3 confirms an optical dropout ~4.5 mag fainter than the quasar (Y_AB=25) at a separation of 0.9 arcsec. The red i_775-Y_105 color of the galaxy and its proximity to the quasar (5 kpc projected if at the quasar redshift) strongly favor an association with the quasar. Although it is much fainter than the quasar it is remarkably bright when compared to field galaxies at this redshift, while showing no evidence for lensing. Both systems may represent late-stage mergers of two massive galaxies, with the observed light for one dominated by powerful ongoing star formation and for the other by rapid black hole growth. Observations of close companions are rare; if major mergers are primarily responsible for high-redshift quasar fueling then the phase when progenitor galaxies can be observed as bright companions is relatively short.

preprint2014arXiv

Does disk fragmentation prevent the formation of supermassive stars in protogalaxies?

Supermassive stars (SMSs; >10^5 Msun) formed in the first protogalaxies with virial temperature T_vir>10^4 K are expected to collapse into seeds of supermassive black hole (SMBHs) in the high-redshift universe (z>7). Fragmentation of the primordial gas is, however, a possible obstacle to SMS formation. We discuss the expected properties of a compact, metal-free, marginally unstable nuclear protogalactic disk, and the fate of the clumps formed in the disk by gravitational instability. Interior to a characteristic radius R_f=few*10^{-2} pc, the disk fragments into massive clumps with M_c~30 Msun. The clumps grow via accretion and migrate inward rapidly on a timescale of ~10^4 yr, which is comparable or shorter than the Kelvin-Helmholz time >10^4 yr. Some clumps may evolve to zero-age main sequence stars and halt gas accretion by radiative feedback, but most of the clumps can migrate inward and merge with the central protostar before forming massive stars. Moreover, we found that dust-induced-fragmentation in metal-enriched gas does not modify these conclusions unless Z> 3*10^{-4} Zsun, because clump migration below this metallicity remains as rapid as in the primordial case. Our results suggest that fragmentation of a compact, metal--poor disk can not prevent the formation of a SMS.

preprint2014arXiv

Evolution in the Escape Fraction of Ionizing Photons and the Decline in Strong Lya Emission from z>6 Galaxies

The rapid decline in the number of strong Lyman Alpha (Lya) emitting galaxies at z > 6 provides evidence for neutral hydrogen in the IGM, but is difficult to explain with plausible models for reionization. We demonstrate that the observed reduction in Lya flux from galaxies at z > 6 can be explained by evolution in the escape fraction of ionizing photons, f_esc. We find that the median observed drop in the fraction of galaxies showing strong Lya emission, as well as the observed evolution of the Lya luminosity function both follow from a small increase in f_esc of Delta f_esc ~ 0.1 from f_esc ~ 0.6 at z ~ 6. This high escape fraction may be at odds with current constraints on the ionising photon escape fraction, which favor smaller values of f_esc < 20%. However, models that invoke a redshift evolution of f_ esc that is consistent with these constraints can suppress the z~7 Lya flux to the observed level, if they also include a small evolution in global neutral fraction of Delta x_HI ~ 0.2. Thus, an evolving escape fraction of ionising photons can be a plausible part of the explanation for evolution in the Lya emission of high redshift galaxies. More generally, our analysis also shows that the drop in the Lya fraction is quantitatively consistent with the observed evolution in the Lya luminosity functions of Lya Emitters.

preprint2014arXiv

H_2 Suppression with Shocking Inflows: Testing a Pathway for Supermassive Black Hole Formation

The presence of quasars at redshifts z > 6 indicates the existence of supermassive black holes (SMBHs) as massive as a few times 10^9 Msun, challenging models for SMBH formation. One pathway is through the direct collapse of gas in T_{vir} > 10^4 K halos; however, this requires the suppression of H_2 cooling to prevent fragmentation. In this paper, we examine a proposed new mechanism for this suppression which relies on cold-mode accretion flows leading to shocks at high densities (n > 10^4 cm^{-3}) and temperatures (T > 10^4 K). In such gas, H_2 is efficiently collisionally dissociated. We use high-resolution numerical simulations to test this idea, demonstrating that such halos typically have lower temperature progenitors, in which cooling is efficient. Those halos do show filamentary flows; however, the gas shocks at or near the virial radius (at low densities), thus preventing the proposed collisional mechanism from operating. We do find that, if we artificially suppress H_2 formation with a high UV background, so as to allow gas in the halo center to enter the high-temperature, high-density "zone of no return", it will remain there even if the UV flux is turned off, collapsing to high density at high temperature. Due to computational limitations, we simulated only three halos. However, we demonstrate, using Monte Carlo calculations of 10^6 halo merger histories, that a few rare halos could assemble rapidly enough to avoid efficient H_2 cooling in all of their progenitor halos, provided that the UV background exceeds J_{21} ~ few at redshifts as high as z ~ 20.

preprint2014arXiv

Lyman edges in supermassive black hole binaries

We propose a new spectral signature for supermassive black hole binaries (SMBHBs) with circumbinary gas disks: a sharp drop in flux blueward of the Lyman limit. A prominent edge is produced if the gas dominating the emission in the Lyman continuum region of the spectrum is sufficiently cold (T < 20,000 K) to contain significant neutral hydrogen. Circumbinary disks may be in this regime if the binary torques open a central cavity in the disk and clear most of the hot gas from the inner region, and if any residual UV emission from the individual BHs is either dim or intermittent. We model the vertical structure and spectra of circumbinary disks using the radiative transfer code TLUSTY, and identify the range of BH masses and binary separations producing a Lyman edge. We find that compact supermassive binaries with orbital periods of ~0.1 - 10 yr, whose gravitational waves (GWs) are expected to be detectable by pulsar timing arrays (PTAs), could have prominent Lyman edges. Such strong spectral edge features are not typically present in AGN spectra and could serve as corroborating evidence for the presence of a SMBHB.

preprint2014arXiv

The Migration of Gap-Opening Planets is not Locked to Viscous Disk Evolution

Most standard descriptions of Type II migration state that massive, gap-opening planets must migrate at the viscous drift rate. This is based on the idea that the disk is separated into an inner and outer region and gas is considered unable to cross the gap. In fact, gas easily crosses the gap on horseshoe orbits, nullifying this necessary premise which would set the migration rate. In this work, it is demonstrated using highly accurate numerical calculations that the actual migration rate is dependent on disk and planet parameters, and can be significantly larger or smaller than the viscous drift rate. In the limiting case of a disk much more massive than the secondary, the migration rate saturates to a constant which is sensitive to disk parameters and is not necessarily of order viscous rate. In the opposite limit of a low-mass disk, the migration rate decreases linearly with disk mass. Steady-state solutions in the low disk mass limit show no pile-up outside the secondary's orbit, and no corresponding drainage of the inner disk.

preprint2013arXiv

A Lyman Break Galaxy in the Epoch of Reionization from HST Grism Spectroscopy

We present observations of a luminous galaxy at redshift z=6.573 --- the end of the reioinization epoch --- which has been spectroscopically confirmed twice. The first spectroscopic confirmation comes from slitless HST ACS grism spectra from the PEARS survey (Probing Evolution And Reionization Spectroscopically), which show a dramatic continuum break in the spectrum at restframe 1216 A wavelength. The second confirmation is done with Keck + DEIMOS. The continuum is not clearly detected with ground-based spectra, but high wavelength resolution enables the Lyman alpha emission line profile to be determined. We compare the line profile to composite line profiles at redshift z=4.5. The Lyman alpha line profile shows no signature of a damping wing attenuation, confirming that the intergalactic gas is ionized at redshift z=6.57. Spectra of Lyman breaks at yet higher redshifts will be possible using comparably deep observations with IR-sensitive grisms, even at redshifts where Lyman alpha is too attenuated by the neutral IGM to be detectable using traditional spectroscopy from the ground.

preprint2013arXiv

Binary Black Hole Accretion From a Circumbinary Disk: Gas Dynamics Inside the Central Cavity

We present the results of 2D hydrodynamical simulations of circumbinary disk accretion using the finite-volume code DISCO. This code solves the 2D viscous Navier-Stokes equations on a high-resolution moving mesh which shears with the fluid flow, greatly reducing advection errors in comparison with a fixed grid. We perform a series of simulations for binary mass ratios in the range 0.026 < q < 1.0, each lasting longer than a viscous time so that we reach a quasi-steady accretion state. In each case, we find that gas is efficiently stripped from the inner edge of the circumbinary disk and enters the cavity along accretion streams, which feed persistent "mini-disks" surrounding each black hole. We find that for q > 0.1, the binary excites eccentricity in the inner region of the circumbinary disk, creating an overdense lump which gives rise to enhanced periodicity in the accretion rate. The dependence of the periodicity on mass ratio may provide a method for observationally inferring mass ratios from measurements of the accretion rate. We also find that for all mass ratios studied, the magnitude of the accretion onto the secondary is sufficient to drive the binary toward larger mass ratio. This suggests a mechanism for biasing mass ratio distributions toward equal mass.

preprint2013arXiv

Constraints on Warm Dark Matter models from high-redshift long gamma-ray bursts

Structures in Warm Dark Matter (WDM) models are exponentially suppressed below a certain scale, characterized by the dark matter particle mass, $m_{\rm x}$. Since structures form hierarchically, the presence of collapsed objects at high-redshifts can set strong lower limits on $m_{\rm x}$. We place robust constraints on $m_{\rm x}$ using recent results from the {\it Swift} database of high-redshift gamma-ray bursts (GRBs). We parameterize the redshift evolution of the ratio between the cosmic GRB rate and star formation rate (SFR) as $\propto (1+z)^α$, thereby allowing astrophysical uncertainties to partially mimic the cosmological suppression of structures in WDM models. Using a maximum likelihood estimator on two different $z>4$ GRB subsamples (including two bursts at $z>8$), we constrain $m_{\rm x} \gtrsim 1.6$-1.8 keV at 95% C.L., when marginalized over a flat prior in $α$. We further estimate that 5 years of a SVOM-like mission would tighten these constraints to $m_{\rm x} \gtrsim 2.3 $ keV. Our results show that GRBs are a powerful probe of high-redshift structures, providing robust and competitive constraints on $m_{\rm x}$.

preprint2013arXiv

Focusing on Warm Dark Matter with Lensed High-redshift Galaxies

We propose a novel use of high-redshift galaxies, discovered in deep Hubble Space Telescope (HST) fields around strong lensing clusters. These fields probe small comoving volumes (about 1000 cubic Mpc) at high magnification (μ > 10), and can detect otherwise inaccessible ultra-faint galaxies. Even a few galaxies found in such small volumes require a very high number density of collapsed dark matter (DM) halos. This implies significant primordial power on small scales, allowing these observations to rule out popular alternatives to standard cold dark matter (CDM) models, such as warm dark matter (WDM). In this work, we analytically compute WDM halo mass functions at z = 10, including the effects of both particle free-streaming and residual velocity dispersion. The preliminary number density corresponding to the two galaxies at z about 10 already detected by the Cluster Lensing And Supernova survey with Hubble (CLASH) constrains the WDM particle mass to m_x > 1 (0.9) keV at 84% (99.9%) confidence limit (for a thermal relic relativistic at decoupling). This limit depends only on the WDM halo mass function and, unlike previous constraints on m_x, is independent of any astrophysical modeling. The forthcoming HST Frontier Fields can significantly tighten these constraints.

preprint2013arXiv

G2 can Illuminate the Black Hole Population near the Galactic Center

Galactic nuclei are expected to be densely populated with stellar and intermediate mass black holes. Exploring this population will have important consequences for the observation prospects of gravitational waves as well as understanding galactic evolution. The gas cloud G2 currently approaching Sgr A* provides an unprecedented opportunity to probe the black hole and neutron star population of the Galactic nucleus. We examine the possibility of a G2-black hole encounter and its detectability with current X-ray satellites, such as Chandra and NuSTAR. We find that multiple encounters are likely to occur close to the pericenter, which may be detectable upon favorable circumstances. This opportunity provides an additional, important science case for leading X-ray observatories to closely follow G2 on its way to the nucleus.

preprint2013arXiv

Gas infall into atomic cooling haloes: on the formation of protogalactic disks and supermassive black holes at z > 10

We have performed cosmo-hydro simulations using the RAMSES code to study atomic cooling (ACHs) haloes at z=10 with masses 5E7Msun<~M<~2E9Msun. We assume primordial gas and H2-cooling and prior star-formation have been suppressed. We analysed 19 haloes (gas and DM) at a resolution of ~10 (proper) pc, selected from a total volume of ~2E3 (comoving) Mpc3. This is the largest statistical hydro-sim. study of ACHs at z>10 to date. We examine the morphology, angular momentum (AM), thermodynamic, and turbulence of these haloes, in order to assess the prevalence of disks and supermassive black holes (SMBHs). We find no correlation between either the magnitude or the direction of the AM of the gas and its parent DM halo. Only 3 haloes form rotationally supported cores. Two of the most massive haloes form massive, compact overdense blobs. These blobs have an accretion rate ~0.5 Msun/yr (at a distance of 100 pc), and are possible sites of SMBH formation. Our results suggest that the degree of rotational support and the fate of the gas in a halo is determined by its large-scale environment and merger history. In particular, the two haloes forming blobs are located at knots of the cosmic web, cooled early on, and experienced many mergers. The gas in these haloes is lumpy and highly turbulent, with Mach N. >~ 5. In contrast, the haloes forming rotationally supported cores are relatively more isolated, located midway along filaments, cooled more recently, and underwent fewer mergers. Thus, the gas in these haloes is less lumpy and less turbulent (Mach <~ 4), and could retain most of its AM. The remaining 14 haloes have intermediate properties. If verified in a larger sample of haloes and with additional physics, our results will have implications for observations of the highest-redshift galaxies and quasars with JWST.

preprint2012arXiv

Evidence of Gunn-Peterson damping wings in high-z quasar spectra: strengthening the case for incomplete reionization

The spectra of several high-redshift (z>6) quasars have shown evidence for a Gunn-Peterson (GP) damping wing, indicating a substantial mean neutral hydrogen fraction (x_HI > 0.03) in the z ~ 6 intergalactic medium (IGM). However, previous analyses assumed that the IGM was uniformly ionized outside of the quasar's HII region. Here we relax this assumption and model patchy reionization scenarios for a range of IGM and quasar parameters. We quantify the impact of these differences on the inferred x_HI, by fitting the spectra of three quasars: SDSS J1148+5251 (z=6.419), J1030+0524 (z=6.308), and J1623+3112 (z=6.247). We find that the best-fit values of x_HI in the patchy models agree well with the uniform case. More importantly, we confirm that the observed spectra favor the presence of a GP damping wing, with peak likelihoods decreasing by factors of > few - 10 when the spectra are modeled without a damping wing. We also find that the Ly alpha absorption spectra, by themselves, cannot distinguish the damping wing in a relatively neutral IGM from a damping wing in a highly ionized IGM, caused either by an isolated neutral patch, or by a damped Ly alpha absorber (DLA). However, neutral patches in a highly ionized universe (x_HI < 0.01), and DLAs with the large required column densities (N_HI > few x 10^{20} cm^{-2}) are both rare. As a result, when we include reasonable prior probabilities for the line of sight (LOS) to intercept either a neutral patch or a DLA at the required distance of ~ 40-60 comoving Mpc away from the quasar, we find strong lower limits on the neutral fraction in the IGM, x_HI > 0.1 (at 95% confidence). This strengthens earlier claims that a substantial global fraction of hydrogen in the z~6 IGM is in neutral form.

preprint2012arXiv

Gas pile-up, gap overflow, and Type 1.5 migration in circumbinary disks: application to supermassive black hole binaries

We study the interaction of a supermassive black hole (SMBH) binary and a standard radiatively efficient thin accretion disk. We examine steady-state configurations of the disk and migrating SMBH system, self-consistently accounting for tidal and viscous torques and heating, radiative diffusion limited cooling, gas and radiation pressure, and the decay of the binary's orbit. We obtain a "phase diagram" of the system as a function of binary parameters, showing regimes in which both the disk structure and migration have a different character. Although massive binaries can create a central gap in the disk at large radii, the tidal barrier of the secondary causes a significant pile-up of gas outside of its orbit, which can lead to the closing of the gap. We find that this spillover occurs at an orbital separation as large as ~200 M_7^{-1/2} gravitational radii, where M = 10^7 M_7 Msun is the total binary mass. If the secondary is less massive than ~10^6 Msun, then the gap is closed before gravitational waves (GWs) start dominating the orbital decay. In this regime, the disk is still strongly perturbed, but the piled-up gas continuously overflows as in a porous dam, and crosses inside the secondary's orbit. The corresponding migration rate, which we label Type 1.5, is slower than the usual limiting cases known as Type I and II migration. Compared to an unperturbed disk, the steady-state disk in the overflowing regime is up to several hundred times brighter in the optical bands. Surveys such as PanSTARRS or LSST may discover the periodic variability of this population of binaries. Our results imply that the circumbinary disks around SMBHs can extend to small radii during the last stages of their merger, when they are detectable by LISA, and may produce coincident electromagnetic (EM) emission similar to active galactic nuclei (AGN).

preprint2012arXiv

Gas pile-up, gap overflow, and Type 1.5 migration in circumbinary disks: general theory

Many astrophysical binaries, from planets to black holes, exert strong torques on their circumbinary accretion disks, and are expected to significantly modify the disk structure. Despite the several decade long history of the subject, the joint evolution of the binary + disk system has not been modeled with self-consistent assumptions for arbitrary mass ratios and accretion rates. Here we solve the coupled binary-disk evolution equations analytically in the strongly perturbed limit, treating the azimuthally-averaged angular momentum exchange between the disk and the binary and the modifications to the density, scale-height, and viscosity self-consistently, including viscous and tidal heating, diffusion limited cooling, radiation pressure, and the orbital decay of the binary. We find a solution with a central cavity and a migration rate similar to those previously obtained for Type-II migration, applicable for large masses and binary separations, and near-equal mass ratios. However, we identify a distinct new regime, applicable at smaller separations and masses, and mass ratio in the range 0.001< q < 0.1. For these systems, gas piles up outside the binary's orbit, but rather than creating a cavity, it continuously overflows as in a porous dam. The disk profile is intermediate between a weakly perturbed disk (producing Type-I migration) and a disk with a gap (with Type-II migration). However, the migration rate of the secondary is typically slower than both Type-I and Type-II rates. We term this new regime "Type-1.5" migration.

preprint2011arXiv

Cosmological Information in Weak Lensing Peaks

Recent studies have shown that the number counts of convergence peaks N(kappa) in weak lensing (WL) maps, expected from large forthcoming surveys, can be a useful probe of cosmology. We follow up on this finding, and use a suite of WL convergence maps, obtained from ray-tracing N-body simulations, to study (i) the physical origin of WL peaks with different heights, and (ii) whether the peaks contain information beyond the convergence power spectrum P_ell. In agreement with earlier work, we find that high peaks (with amplitudes >~ 3.5 sigma, where sigma is the r.m.s. of the convergence kappa) are typically dominated by a single massive halo. In contrast, medium-height peaks (~0.5-1.5 sigma) cannot be attributed to a single collapsed dark matter halo, and are instead created by the projection of multiple (typically, 4-8) halos along the line of sight, and by random galaxy shape noise. Nevertheless, these peaks dominate the sensitivity to the cosmological parameters w, sigma_8, and Omega_m. We find that the peak height distribution and its dependence on cosmology differ significantly from predictions in a Gaussian random field. We directly compute the marginalized errors on w, sigma_8, and Omega_m from the N(kappa) + P_ell combination, including redshift tomography with source galaxies at z_s=1 and z_s=2. We find that the N(kappa) + P_ell combination has approximately twice the cosmological sensitivity compared to P_ell alone. These results demonstrate that N(kappa) contains non-Gaussian information complementary to the power spectrum.

preprint2011arXiv

Probing Cosmology with Weak Lensing Minkowski Functionals

In this paper, we show that Minkowski Functionals (MFs) of weak gravitational lensing (WL) convergence maps contain significant non-Gaussian, cosmology-dependent information. To do this, we use a large suite of cosmological ray-tracing N-body simulations to create mock WL convergence maps, and study the cosmological information content of MFs derived from these maps. Our suite consists of 80 independent 512^3 N-body runs, covering seven different cosmologies, varying three cosmological parameters Omega_m, w, and sigma_8 one at a time, around a fiducial LambdaCDM model. In each cosmology, we use ray-tracing to create a thousand pseudo-independent 12 deg^2 convergence maps, and use these in a Monte Carlo procedure to estimate the joint confidence contours on the above three parameters. We include redshift tomography at three different source redshifts z_s=1, 1.5, 2, explore five different smoothing scales theta_G=1, 2, 3, 5, 10 arcmin, and explicitly compare and combine the MFs with the WL power spectrum. We find that the MFs capture a substantial amount of information from non-Gaussian features of convergence maps, i.e. beyond the power spectrum. The MFs are particularly well suited to break degeneracies and to constrain the dark energy equation of state parameter w (by a factor of ~ three better than from the power spectrum alone). The non-Gaussian information derives partly from the one-point function of the convergence (through V_0, the "area" MF), and partly through non-linear spatial information (through combining different smoothing scales for V_0, and through V_1 and V_2, the boundary length and genus MFs, respectively). In contrast to the power spectrum, the best constraints from the MFs are obtained only when multiple smoothing scales are combined.

preprint2010arXiv

Constraining Intra-cluster Gas Models with AMiBA13

Clusters of galaxies have been used extensively to determine cosmological parameters. A major difficulty in making best use of Sunyaev-Zel'dovich (SZ) and X-ray observations of clusters for cosmology is that using X-ray observations it is difficult to measure the temperature distribution and therefore determine the density distribution in individual clusters of galaxies out to the virial radius. Observations with the new generation of SZ instruments are a promising alternative approach. We use clusters of galaxies drawn from high-resolution adaptive mesh refinement (AMR) cosmological simulations to study how well we should be able to constrain the large-scale distribution of the intra-cluster gas (ICG) in individual massive relaxed clusters using AMiBA in its configuration with 13 1.2-m diameter dishes (AMiBA13) along with X-ray observations. We show that non-isothermal beta models provide a good description of the ICG in our simulated relaxed clusters. We use simulated X-ray observations to estimate the quality of constraints on the distribution of gas density, and simulated SZ visibilities (AMiBA13 observations) for constraints on the large-scale temperature distribution of the ICG. We find that AMiBA13 visibilities should constrain the scale radius of the temperature distribution to about 50% accuracy. We conclude that the upgraded AMiBA, AMiBA13, should be a powerful instrument to constrain the large-scale distribution of the ICG.

preprint2010arXiv

Cosmology with Standard Sirens: the Importance of the Shape of the Lensing Magnification Distribution

The gravitational waves (GWs) emitted by inspiraling binary black holes, expected to be detected by the Laser Interferometer Space Antenna (LISA), could be used to determine the luminosity distance to these sources with the unprecedented precision of <~ 1%. We study cosmological parameter constraints from such standard sirens, in the presence of gravitational lensing by large-scale structure. Lensing introduces magnification with a probability distribution function (PDF) whose shape is highly skewed and depends on cosmological parameters. We use Monte-Carlo simulations to generate mock samples of standard sirens, including a small intrinsic scatter, as well as the additional, larger scatter from lensing, in their inferred distances. We derive constraints on cosmological parameters, by simultaneously fitting the mean and the distribution of the residuals on the distance vs redshift (d_L - z) Hubble diagram. We find that for standard sirens at redshift z ~ 1, the sensitivity to a single cosmological parameter, such as the matter density Omega_m, or the dark energy equation of state w, is ~ 50%-80% tighter when the skewed lensing PDF is used, compared to the sensitivity derived from a Gaussian PDF with the same variance. When these two parameters are constrained simultaneously, the skewness yields a further enhanced improvement (by ~ 120%), owing to the correlation between the parameters. The sensitivity to the amplitude of the matter power spectrum, sigma_8 from the cosmological dependence of the PDF alone, however, is ~ 20% worse than that from the Gaussian PDF. At higher redshifts, the PDF resembles a Gaussian more closely, and the effects of the skewness become less prominent. These results highlight the importance of obtaining an accurate and reliable PDF of the lensing convergence, in order to realize the full potential of standard sirens as cosmological probes.

preprint2009arXiv

Astro2010 Decadal Survey Whitepaper: Coordinated Science in the Gravitational and Electromagnetic Skies

It is widely expected that the coming decade will witness the first direct detection of gravitational waves (GWs). The ground-based LIGO and Virgo GW observatories are being upgraded to advanced sensitivity, and are expected to observe a significant binary merger rate. The launch of The Laser Interferometer Space Antenna (LISA) would extend the GW window to low frequencies, opening new vistas on dynamical processes involving massive (M >~ 10^5 M_Sun) black holes. GW events are likely to be accompanied by electromagnetic (EM) counterparts and, since information carried electromagnetically is complementary to that carried gravitationally, a great deal can be learned about an event and its environment if it becomes possible to measure both forms of radiation in concert. Measurements of this kind will mark the dawn of trans-spectral astrophysics, bridging two distinct spectral bands of information. The aim of this whitepaper is to articulate future directions in both theory and observation that are likely to impact broad astrophysical inquiries of general interest. What will EM observations reflect on the nature and diversity of GW sources? Can GW sources be exploited as complementary probes of cosmology? What cross-facility coordination will expand the science returns of gravitational and electromagnetic observations?

preprint2009arXiv

Probing Cosmology and Galaxy Cluster Structure with the Sunyaev-Zel'dovich Decrement vs. X-ray Temperature Scaling Relation

Scaling relations among galaxy cluster observables, which will become available in large future samples of galaxy clusters, could be used to constrain not only cluster structure, but also cosmology. We study the utility of this approach, employing a physically motivated parametric model to describe cluster structure, and applying it to the expected relation between the Sunyaev-Zel'dovich decrement (S_ν) and the emission-weighted X-ray temperature (T_ew). The slope and normalization of the entropy profile, the concentration of the dark matter potential, the pressure at the virial radius, and the level of non-thermal pressure support, as well as the mass and redshift-dependence of these quantities are described by free parameters. With a suitable choice of fiducial parameter values, the cluster model satisfies several existing observational constraints. We employ a Fisher matrix approach to estimate the joint errors on cosmological and cluster structure parameters from a measurement of S_νvs. T_ew in a future survey. We find that different cosmological parameters affect the scaling relation differently: predominantly through the baryon fraction (Ω_m and Ω_b), the virial overdensity (w_0 and w_a for low-z clusters) and the angular diameter distance (w_0, w_a for high-z clusters; Ω_DE and h). We find that the cosmology constraints from the scaling relation are comparable to those expected from the number counts (dN/dz) of the same clusters. The scaling relation approach is relatively insensitive to selection effects and it offers a valuable consistency check; combining the information from the scaling relation and dN/dz is also useful to break parameter degeneracies and help disentangle cluster physics from cosmology.

preprint2009arXiv

Probing Reionization with Quasar Spectra: the Impact of the Intrinsic Lyman-alpha Emission Line Shape Uncertainty

Arguably the best hope of understanding the tail end of the reionization of the intergalactic medium (IGM) at redshift z > 6 is through the detection and characterization of the Gunn-Peterson (GP) damping wing absorption of the IGM in bright quasar spectra. However, the use of quasar spectra to measure the IGM damping wing requires a model of the quasar's intrinsic Lyman-alpha emission line. Here we quantify the uncertainties in the intrinsic line shapes, and how those uncertainties affect the determination of the IGM neutral fraction. We have assembled a catalog of high-resolution HST spectra of the emission lines of unobscured low-redshift quasars, and have characterized the variance in the shapes of their lines. We then add simulated absorption from the high-redshift IGM to these quasar spectra in order to determine the corresponding uncertainties in reionization constraints using current and future samples of z > 6 quasar spectra. We find that, if the redshift of the Lyman-alpha emission line is presumed to coincide with the systemic redshift determined from metal lines, the inferred IGM neutral fraction is systematically biased to low values due to a systematic blueshift of the Lyman-alpha line relative to the metal lines. If a similar blueshift persists in quasars at z > 6, this bias strengthens previous claims of a significant neutral hydrogen fraction at z ~ 6. This technique is capable of making a robust distinction between a highly ionized (x_IGM ~ 10^-3) and a neutral (x_IGM = 1) IGM with even a few bright quasars.

preprint2009arXiv

Radiative Feedback in Relic HII Regions at High-Redshift

UV radiation from early astrophysical sources could have a large impact on subsequent star formation in nearby protogalaxies. Here we study the radiative feedback from the first, short-lived stars using hydrodynamical simulations with transient UV backgrounds (UVBs) and persistent Lyman-Werner backgrounds (LWBs) of varying intensity. We extend our prior work in Mesinger et al. (2006), by studying a more typical region whose proto-galaxies form at lower redshifts, z~13-20, in the epoch likely preceding the bulk of reionization. We confirm our previous results that feedback in the relic HII regions resulting from such transient radiation, is itself transient. Feedback effects dwindle away after ~30% of the Hubble time, and the same critical specific intensity of J_UV~0.1 x 10^{-21} ergs/s/cm^2/Hz/sr separates positive and negative feedback regimes. Additionally, we discover a second episode of eventual positive feedback in halos which have not yet collapsed when their progenitor regions were exposed to the transient UVB. This eventual positive feedback appears in all runs, regardless of the strength of the UVB. However, this feedback regime is very sensitive to the presence of Lyman-Werner radiation, and notable effects disappear under fairly modest background intensities of J_LW>10^{-3} x 10^{-21} ergs/s/cm^2/Hz/sr. We conclude that UV radiative feedback in relic HII regions, although a complicated process, seems unlikely to have a major impact on the progress of cosmological reionization, provided that present estimates of the lifetime and luminosity of a PopIII star are accurate. More likely is that the build-up of the LWB ultimately governs the feedback strength until a persistent UV background can be established. [abridged]

preprint2009arXiv

Supermassive Black Hole Formation by Direct Collapse: Keeping Protogalactic Gas H_2--Free in Dark Matter Halos with Virial Temperatures T_vir >~ 10^4 K

In the absence of H_2 molecules, the primordial gas in early dark matter halos with virial temperatures just above T_vir >~ 10^4 K cools by collisional excitation of atomic H. Although it cools efficiently, this gas remains relatively hot, at a temperature near T ~ 8000 K, and consequently might be able to avoid fragmentation and collapse directly into a supermassive black hole (SMBH). In order for H_2--formation and cooling to be strongly suppressed, the gas must be irradiated by a sufficiently intense ultraviolet (UV) flux. We performed a suite of three--dimensional hydrodynamical adaptive mesh refinement (AMR) simulations of gas collapse in three different protogalactic halos with T_vir >~ 10^4 K, irradiated by a UV flux with various intensities and spectra. We determined the critical specific intensity, Jcrit, required to suppress H_2 cooling in each of the three halos. For a hard spectrum representative of metal--free stars, we find (in units of 10^{-21} erg s^{-1} Hz^{-1} sr^{-1} cm^{-2}) 10^4<Jcrit<10^5, while for a softer spectrum, which is characteristic of a normal stellar population, and for which H^{-} --dissociation is important, we find 30<Jcrit<300. These values are a factor of 3--10 lower than previous estimates. We attribute the difference to the higher, more accurate H_2 collisional dissociation rate we adopted. The reduction in Jcrit exponentially increases the number of rare halos exposed to super--critical radiation. When H_2 cooling is suppressed, gas collapse starts with a delay, but it ultimately proceeds more rapidly. The infall velocity is near the increased sound speed, and an object as massive as M ~ 10^5 solar mass may form at the center of these halos, compared to the M ~ 10^2 solar mass stars forming when H_2--cooling is efficient.

preprint2008arXiv

Challenges to the DGP Model from Horizon-Scale Growth and Geometry

We conduct a Markov Chain Monte Carlo study of the Dvali-Gabadadze-Porrati (DGP) self-accelerating braneworld scenario given the cosmic microwave background (CMB) anisotropy, supernovae and Hubble constant data by implementing an effective dark energy prescription for modified gravity into a standard Einstein-Boltzmann code. We find no way to alleviate the tension between distance measures and horizon scale growth in this model. Growth alterations due to perturbations propagating into the bulk appear as excess CMB anisotropy at the lowest multipoles. In a flat cosmology, the maximum likelihood DGP model is nominally a 5.3 sigma poorer fit than Lambda CDM. Curvature can reduce the tension between distance measures but only at the expense of exacerbating the problem with growth leading to a 4.8 sigma result that is dominated by the low multipole CMB temperature spectrum. While changing the initial conditions to reduce large scale power can flatten the temperature spectrum, this also suppresses the large angle polarization spectrum in violation of recent results from WMAP5. The failure of this model highlights the power of combining growth and distance measures in cosmology as a test of gravity on the largest scales.

preprint2008arXiv

Observing the First Stars and Black Holes

The high sensitivity of JWST will open a new window on the end of the cosmological dark ages. Small stellar clusters, with a stellar mass of several 10^6 M_sun, and low-mass black holes (BHs), with a mass of several 10^5 M_sun should be directly detectable out to redshift z=10, and individual supernovae (SNe) and gamma ray burst (GRB) afterglows are bright enough to be visible beyond this redshift. Dense primordial gas, in the process of collapsing from large scales to form protogalaxies, may also be possible to image through diffuse recombination line emission, possibly even before stars or BHs are formed. In this article, I discuss the key physical processes that are expected to have determined the sizes of the first star-clusters and black holes, and the prospect of studying these objects by direct detections with JWST and with other instruments. The direct light emitted by the very first stellar clusters and intermediate-mass black holes at z>10 will likely fall below JWST's detection threshold. However, JWST could reveal a decline at the faint-end of the high-redshift luminosity function, and thereby shed light on radiative and other feedback effects that operate at these early epochs. JWST will also have the sensitivity to detect individual SNe from beyond z=10. In a dedicated survey lasting for several weeks, thousands of SNe could be detected at z>6, with a redshift distribution extending to the formation of the very first stars at z>15. Using these SNe as tracers may be the only method to map out the earliest stages of the cosmic star-formation history. Finally, we point out that studying the earliest objects at high redshift will also offer a new window on the primordial power spectrum, on 100 times smaller scales than probed by current large-scale structure data.

preprint2008arXiv

The Evolution of the M_BH-sigma relation Inferred from the Age Distribution of Local Early-Type Galaxies and AGN Evolution

We utilize the local velocity dispersion function (VDF) of spheroids, together with their inferred age--distributions, to predict the VDF at higher redshifts (0<z<6), under the assumption that (i) most of the stars in each nearby spheroid formed in a single episode, and (ii) the velocity dispersion sigma remained nearly constant afterward. We assume further that a supermassive black hole (BH) forms concurrently with the stars, and within ~1 Gyr of the formation of the potential well of the spheroid, and that the relation between the mass of the BH and host velocity dispersion maintains the form M_BH ~ sigma^{beta} with beta~4, but with the normalization allowed to evolve with redshift as ~(1+z)^{alpha}. We compute the BH mass function associated with the VDF at each redshift, and compare the accumulated total BH mass density with that inferred from the integrated quasar luminosity function (LF; the so--called Soltan argument). This comparison is insensitive to the assumed duty cycle or Eddington ratio of quasar activity, and we find that the match between the two BH mass densities favors a relatively mild redshift evolution, with alpha ~ 0.26, with a positive evolution as strong as alpha>1.3 excluded at the 99% confidence level. A direct match between the characteristic BH mass in the VDF--based and quasar LF--based BH mass functions also yields a mean Eddington ratio of lambda ~ 0.5-1 that is roughly constant within 0<z<3. A strong positive evolution in the M_BH-sigma relation is still allowed by the data if galaxies increase, on average, their velocity dispersions since the moment of formation, due to dissipative processes. If we assume that the mean velocity dispersion of the host galaxies evolves as sigma(z)=sigma(0)*(1+z)^{-gamma}, we find a lower limit of gamma>0.23 for alpha>1.5. abridged

preprint2007arXiv

Constraints on Reionization and Source Properties from the Absorption Spectra of z>6.2 Quasars

We make use of hydrodynamical simulations of the intergalactic medium (IGM) to create model quasar absorption spectra. We compare these model spectra with the observed Keck spectra of three z>6.2 quasars with full Gunn-Peterson troughs: SDSS J1148+5251 (z=6.42), SDSS J1030+0524 (z=6.28), and SDSS J1623+3112 (z=6.22). We fit the probability density distributions (PDFs) of the observed Ly alpha optical depths with those generated from the simulation, by exploring a range of values for the size of the quasar's surrounding HII region, R_S, the volume-weighted mean neutral hydrogen fraction in the ambient IGM, x_H, and the quasar's ionizing photon emissivity, N_Q. In order to avoid averaging over possibly large sightline-to-sightline fluctuations in IGM properties, we analyze each observed quasar independently. We find the following results for J1148+5251, J1030+0524, and J1623+3112: The best-fit sizes R_S are 40, 41, and 29 (comoving) Mpc, respectively. These constraints are tight, with only ~ 10% uncertainties, comparable to those caused by redshift--determination errors. The best-fit values of N_Q are 2.1, 1.3, and 0.9 x 10^{57} s^{-1}, respectively, with a factor of ~ 2 uncertainty in each case. Finally, the best--fit values of x_H are 0.16, 1.0, and 1.0, respectively. The uncertainty in the case of J1148+5251 is large, and x_H is not well constrained. However, for both J1030+0524 and J1623+3112, we find a significant lower limit of x_H >~ 0.033. Our method is different from previous analyses of the GP absorption spectra of these quasars, and our results strengthen the evidence that the rapid end--stage of reionization is occurring near z ~ 6.

preprint2006arXiv

Reconstructing the Cosmic Evolution of Quasars from the Age Distribution of Local Early-Type Galaxies

We use the spectra of 22,000 nearby early-type galaxies from the Sloan Digital Sky Survey (SDSS) to determine the age distribution of these galaxies as a function of their velocity dispersion sigma_v in the range 100 km/s < sigma_v < 280 km/s. We then combine the inferred age-distributions with the local abundance of spheroids, including early-type galaxies and late-type bulges, to predict the evolution of the quasar luminosity function (LF) in the redshift range 0<z<6. We make the following simple assumptions: (i) the formation of stars in each galaxy, at the epoch identified with the mean mass-weighted stellar age, is accompanied by the prompt assembly of the nuclear supermassive black hole (SMBH); (ii) the mass of the SMBH obeys the M_bh-sigma_v correlation observed in nearby galaxies; (iii) the SMBH radiates at a fraction f_Edd of the Eddington limit for a fixed duration t_Q, and is identified as a luminous quasar during this epoch, (iv) the intrinsic dispersions in the Eddington ratio and the M_bh-sigma_v relation produce a combined scatter of Delta(log L_Q) around the mean logarithmic quasar luminosity <log L_Q> at fixed sigma_v. These assumptions require that the SMBH remnants of quasars with bolometric luminosity below L_bol=10^{12.5} f_Edd L_sun reside predominantly in bulges of late type galaxies. We find that evolution of the observed quasar LF can be fit over the entire redshift range in this simple model, 0<z<6 with the choices of Delta(log L_Q)=0.6-0.9, t_Q= (6-8)x10^7 yr, and <f_Edd>=0.3-0.5. We find no evidence that any of the model parameters evolves with redshift, supporting the strong connection between the formation of stars and nuclear SMBHs in spheroids.

preprint2005arXiv

A Snapshot Survey for Gravitational Lenses Among z>=4.0 Quasars: II. Constraints on the 4.0<z<4.5 Quasar Population

We report on i-band snapshot observations of 157 Sloan Digital Sky Survey (SDSS) quasars at 4<z<5.4 using the Advanced Camera for Surveys on the Hubble Space Telescope (HST) to search for evidence of gravitational lensing of these sources. None of the quasars appear to be strongly lensed and multiply imaged at the angular resolution (~0.1") and sensitivity of HST. The non-detection of strong lensing in these systems constrains the z=4-5 luminosity function to an intrinsic slope of beta>-3.8 (3 sigma), assuming a break in the quasar luminosity function at M_1450^star=-24.5. This constraint is considerably stronger than the limit of beta>-4.63 obtained from the absence of lensing in four z>5.7 quasars. Such constraints are important for our understanding of the true space density of high-redshift quasars and the ionization state of the early universe.

preprint2005arXiv

Weighing Neutrinos with Galaxy Cluster Surveys

Large future galaxy cluster surveys, combined with cosmic microwave background observations, can achieve a high sensitivity to the masses of cosmologically important neutrinos. We show that a weak lensing selected sample of ~100,000 clusters could tighten the current upper bound on the sum of masses of neutrino species by an order of magnitude, to a level of 0.03 eV. Since this statistical sensitivity is below the best existing lower limit on the mass of at least one neutrino species, a future detection is likely, provided that systematic errors can be controlled to a similar level.

preprint2004arXiv

Rate Coefficient for H^+ + H_2(X ^1Sigma_g^+, nu=0, J=0) --> H(1s) + H_2^+ Charge Transfer and Some Cosmological Implications

Krstic has carried out the first quantum mechanical calculations near threshold for the charge transfer (CT) process H^+ + H_2(X ^1Sigma_g^+, nu=0, J=0) --> H(1s) + H_2^+. These results are relevant for models of primordial galaxy and first star formation that require reliable atomic and molecular data for obtaining the early universe hydrogen chemistry. Using the results of Krstic, we calculate the relevant CT rate coefficient for temperatures between 100 and 30,000 K. We also present a simple fit which can be readily implemented into early universe chemical models. Additionally, we explore how the range of previously published data for this reaction translates into uncertainties in the predicted gas temperature and H_2 relative abundance in a collapsing primordial gas cloud. Our new data significantly reduce these cosmological uncertainties that are due to the uncertainties in the previously published CT rate coefficients.

preprint2004arXiv

The Formation and Evolution of the First Massive Black Holes

The first massive astrophysical black holes likely formed at high redshifts (z>10) at the centers of low mass (~10^6 Msun) dark matter concentrations. These black holes grow by mergers and gas accretion, evolve into the population of bright quasars observed at lower redshifts, and eventually leave the supermassive black hole remnants that are ubiquitous at the centers of galaxies in the nearby universe. The astrophysical processes responsible for the formation of the earliest seed black holes are poorly understood. The purpose of this review is threefold: (1) to describe theoretical expectations for the formation and growth of the earliest black holes within the general paradigm of hierarchical cold dark matter cosmologies, (2) to summarize several relevant recent observations that have implications for the formation of the earliest black holes, and (3) to look into the future and assess the power of forthcoming observations to probe the physics of the first active galactic nuclei.

preprint2003arXiv

A Snapshot Survey for Gravitational Lenses Among z>=4.0 Quasars: I. The z>5.7 Sample

Over the last few years, the Sloan Digital Sky Survey (SDSS) has discovered several hundred quasars with redshift between 4.0 and 6.4. Including the effects of magnification bias, one expects a priori that an appreciable fraction of these objects are gravitationally lensed. We have used the Advanced Camera for Surveys on the Hubble Space Telescope to carry out a snapshot imaging survey of high-redshift SDSS quasars to search for gravitationally split lenses. This paper, the first in a series reporting the results of the survey, describes snapshot observations of four quasars at z = 5.74, 5.82, 5.99 and 6.30, respectively. We find that none of these objects has a lensed companion within 5 magnitudes with a separation larger than 0.3 arcseconds; within 2.5 magnitudes, we can rule out companions within 0.1 arcseconds. Based on the non-detection of strong lensing in these four systems, we constrain the z~6 luminosity function to a slope of beta>-4.63 (3 sigma), assuming a break in the quasar luminosity function at M_{1450}^*=-24.0. We discuss the implications of this constraint on the ionizing background due to quasars in the early universe. Given that these quasars are not highly magnified, estimates of the masses of their central engines by the Eddington argument must be taken seriously, possibly challenging models of black hole formation.

preprint2003arXiv

Do We Need Stars to Reionize the Universe at High Redshifts? Early Reionization by Decaying Heavy Sterile Neutrinos

A remarkable result of the Wilkinson Microwave Anisotropy Probe (WMAP) observations is that the universe was significantly reionized at large redshifts. The standard explanation is that massive stars formed early and reionized the universe around redshift z=17. Here we explore an alternative possibility, in which the universe was reionized in two steps. An early boost of reionization is provided by a decaying sterile neutrino, whose decay products, relativistic electrons, result in partial ionization of the smooth gas. We demonstrate that a neutrino with a mass of m_nu ~ 200 MeV and a decay time of t ~ 4 * 10^{15} s can account for the electron scattering optical depth tau=0.16 measured by WMAP without violating existing astrophysical limits on the cosmic microwave and gamma ray backgrounds. Reionization is then completed by subsequent star formation at lower redshifts. This scenario alleviates constraints on structure formation models with reduced small-scale power, such as those with a running or tilted scalar index, or warm dark matter models.

preprint2001arXiv

Constraints on Omega_m, Omega_L, and Sigma_8, from Galaxy Cluster Redshift Distributions

We show that the counts of galaxy clusters in future deep cluster surveys can place strong constraints on the matter density, Omega_m, the vacuum energy density, Omega_L, and the normalization of the matter power spectrum, sigma_8. Degeneracies between these parameters are different from those in studies of either high--redshift type Ia Supernovae (SNe), or cosmic microwave background (CMB) anisotropies. Using a mass threshold for cluster detection expected to be typical for upcoming SZE surveys, we find that constraints on Omega_m and sigma_8 at the level of roughly 5% or better can be expected, assuming redshift information is known at least to z=0.5 and in the absence of significant systematic errors. Without information past this redshift, Omega_L is constrained to 25%. With complete redshift information, deep (M_{lim}= 10^{14}h^{-1}{M_sun}), relatively small solid angle (roughly 12 {deg}^2) surveys can further constrain Omega_L to an accuracy of 15%, while large solid angle surveys with ground-based large-format bolometer arrays could measure Omega_L to a precision of 4% or better.

preprint2001arXiv

What does the local black hole mass distribution tell us about the evolution of the quasar luminosity function?

We present a robust method to derive the duty cycle of QSO activity based on the empirical QSO luminosity function and on the present-day linear relation between the masses of supermassive black holes and those of their spheroidal host stellar systems. It is found that the duty cycle is substantially less than unity, with characteristic values in the range $3-6\times 10^{-3}$. Finally, we tested the expectation that the QSO luminosity evolution and the star formation history should be roughly parallel, as a consequence of the above--mentioned relation between BH and galaxy masses.

preprint2000arXiv

Constraining the Lifetime of Quasars from their Spatial Clustering

The lifetime t_Q of the luminous phase of quasars is constrained by current observations to be between 10^6 and 10^8 years, but is otherwise unkown. We model the quasar luminosity function in detail in the optical and X-ray bands using the Press-Schechter formalism, and show that the expected clustering of quasars depends strongly on their assumed lifetime. We quantify this dependence, and find that existing measurements of the correlation length of quasars are consistent with the range 10^6 < t_Q < 10^8 years. We then show that future measurements of the power spectrum of quasars out to z=3, from the 2dF or Sloan Digital Sky Survey, can significantly improve this constraint, and in principle allow a precise determination of t_Q. We estimate the systematic errors introduced by uncertainties in the modeling of the quasar-halo relationship, as well as by the possible existence of obscured quasars.

preprint2000arXiv

Constraints on Cosmological Parameters from Future Galaxy Cluster Surveys

We study the expected redshift evolution of galaxy cluster abundance between 0 < z < 3 in different cosmologies, including the effects of the cosmic equation of state parameter w=p/rho. Using the halo mass function obtained in recent large scale numerical simulations, we model the expected cluster yields in a 12 deg^2 Sunyaev-Zeldovich Effect (SZE) survey and a deep 10^4 deg^2 X-ray survey over a wide range of cosmological parameters. We quantify the statistical differences among cosmologies using both the total number and redshift distribution of clusters. Provided that the local cluster abundance is known to a few percent accuracy, we find only mild degeneracies between w and either Omega_m or h. As a result, both surveys will provide improved constraints on Omega_m and w. The Omega_m-w degeneracy from both surveys is complementary to those found either in studies of CMB anisotropies or of high-redshift Supernovae (SNe). As a result, combining these surveys together with either CMB or SNe studies can reduce the statistical uncertainty on both w and Omega_m to levels below what could be obtained by combining only the latter two data sets. Our results indicate a formal statistical uncertainty of about 3% (68% confidence) on both Omega_m and w when the SZE survey is combined with either the CMB or SN data; the large number of clusters in the X-ray survey further suppresses the degeneracy between w and both Omega_m and h. Systematics and internal evolution of cluster structure at the present pose uncertainties above these levels. We briefly discuss and quantify the relevant systematic errors. By focusing on clusters with measured temperatures in the X-ray survey, we reduce our sensitivity to systematics such as non-standard evolution of internal cluster structure.

preprint2000arXiv

Galaxy Cluster Baryon Fractions, Cluster Surveys and Cosmology

The properties of nearby galaxy clusters limit the range of cosmological parameters consistent with our universe. We describe the limits which arise from studies of the intracluster medium (ICM) mass fraction fICM and consideration of the possible sources of systematic error: Omega_M<0.44h_{50}^{-1/2} at 95% confidence. We emphasize that independent of Type Ia supernovae (SNe Ia) observations, this cluster study, taken together with published cosmic microwave background (CMB) anisotropy studies, indicates a non-zero quintessence or dark energy component Omega_Q>0. We then discuss future galaxy cluster surveys which will probe the abundance of galaxy clusters to intermediate and high redshift. We investigate the sensitivity of these surveys to the cosmological density parameter Omega_M and the equation of state parameter w of any quintessence component. In particular, we show that cluster survey constraints from a proposed large solid angle X-ray survey are comparable in precision and complementary in nature to constraints expected from future CMB anisotropy and SNe Ia studies.

preprint1999arXiv

The Radiative Feedback of the First Cosmological Objects

In hierarchical models of structure formation, an early cosmic UV background (UVB) is produced by the small (T_vir < 10^4 K) halos that collapse before reionization. The UVB at energies below 13.6eV suppresses the formation of stars or black holes inside small halos, by photo-dissociating their only cooling agent, molecular H2. We self-consistently compute the buildup of the early UVB in Press-Schechter models, coupled with H2 photo-dissociation both in the intergalactic medium (IGM), and inside virialized halos. We find that the intergalactic H2 has a negligible effect on the UVB, both because its initial optical depth is small (tau<0.1), and because it is photo-dissociated at an early stage. If the UV sources in the first collapsed halos are stars, then their UV flux suppresses further star-formation inside small halos. This results in a pause in the buildup of the UVB, and reionization is delayed until larger halos (T_vir> 10^4 K) collapse. If the small halos host mini-quasars with hard spectra extending to approximately 1 keV, then their X-rays balance the effects of the UVB, the negative feedback does not occur, and reionization can be caused by the small halos.

Zoltan Haiman

What is connected

Connect this record

See the researcher in context

Building this map preview

76 published item(s)

Lomb-Scargle periodograms struggle with non-sinusoidal supermassive BH binary signatures in quasar lightcurves

Forming massive seed black holes in high-redshift quasar host progenitors

The Effect of Mission Duration on LISA Science Objectives

Gas-driven inspiral of binaries in thin accretion disks

Multi-Messenger Constraints on Magnetic Fields in Merging Black Hole-Neutron Star Binaries

Circumbinary Disks: Accretion and Torque as a Function of Mass Ratio and Disk Viscosity

Eccentric Black Hole Mergers in Active Galactic Nuclei

Formation and Evolution of Compact Object Binaries in AGN Disks

Getting Ready for LISA: The Data, Support and Preparation Needed to Maximize US Participation in Space-Based Gravitational Wave Science

Making a supermassive star by stellar bombardment

Mass-gap Mergers in Active Galactic Nuclei

Self-consistent Semi-analytic Modeling of Feedback During Primordial Star Formation and Reionization

Spikey: Self-Lensing Flares from Eccentric SMBH Binaries

Spin Evolution of Stellar-mass Black Hole Binaries in Active Galactic Nuclei

Testing the relativistic Doppler boost hypothesis for the binary candidate quasar PG1302-102 with multi-band Swift data

The Quest for Dual and Binary Supermassive Black Holes: A Multi-Messenger View

A transition in circumbinary accretion discs at a binary mass ratio of 1:25

Assisted Inspirals of Stellar Mass Black Holes Embedded in AGN Disks: Solving the "Final AU Problem"

Formation of Massive Population III Galaxies through Photoionization Feedback: A Possible Explanation for CR7

Gravitational wave background from Population III binary black holes consistent with cosmic reionization

Hyper-Eddington accretion flows onto massive black holes

Intermediate-mass black holes from Population III remnants in the first galactic nuclei

Is there a maximum mass for black holes in galactic nuclei?

The Origin of Weak Lensing Convergence Peaks

A Statistical Method to Search for Recoiling Supermassive Black Holes in Active Galactic Nuclei

Cosmology Constraints from the Weak Lensing Peak Counts and the Power Spectrum in CFHTLenS

Direct collapse black hole formation from synchronized pairs of atomic cooling halos

Emulating the CFHTLenS Weak Lensing data: Cosmological Constraints from moments and Minkowski functionals

High-Redshift Star Formation in a Time-Dependent Lyman-Werner Background

Limits on Population III star formation in minihaloes implied by Planck

Looking for Population III stars with He II line intensity mapping

The impact of spurious shear on cosmological parameter estimates from weak lensing observables

A No-Go Theorem for Direct Collapse Black Holes Without a Strong Ultraviolet Background

Binary Black Hole Accretion During Inspiral and Merger

Characteristic Signatures in the Thermal Emission from Accreting Binary Black Holes

Close companions to two high-redshift quasars

Does disk fragmentation prevent the formation of supermassive stars in protogalaxies?

Evolution in the Escape Fraction of Ionizing Photons and the Decline in Strong Lya Emission from z>6 Galaxies

H_2 Suppression with Shocking Inflows: Testing a Pathway for Supermassive Black Hole Formation

Lyman edges in supermassive black hole binaries

The Migration of Gap-Opening Planets is not Locked to Viscous Disk Evolution

A Lyman Break Galaxy in the Epoch of Reionization from HST Grism Spectroscopy

Binary Black Hole Accretion From a Circumbinary Disk: Gas Dynamics Inside the Central Cavity

Constraints on Warm Dark Matter models from high-redshift long gamma-ray bursts

Focusing on Warm Dark Matter with Lensed High-redshift Galaxies

G2 can Illuminate the Black Hole Population near the Galactic Center

Gas infall into atomic cooling haloes: on the formation of protogalactic disks and supermassive black holes at z > 10

Evidence of Gunn-Peterson damping wings in high-z quasar spectra: strengthening the case for incomplete reionization

Gas pile-up, gap overflow, and Type 1.5 migration in circumbinary disks: application to supermassive black hole binaries

Gas pile-up, gap overflow, and Type 1.5 migration in circumbinary disks: general theory

Cosmological Information in Weak Lensing Peaks

Probing Cosmology with Weak Lensing Minkowski Functionals

Constraining Intra-cluster Gas Models with AMiBA13

Cosmology with Standard Sirens: the Importance of the Shape of the Lensing Magnification Distribution

Astro2010 Decadal Survey Whitepaper: Coordinated Science in the Gravitational and Electromagnetic Skies

Probing Cosmology and Galaxy Cluster Structure with the Sunyaev-Zel'dovich Decrement vs. X-ray Temperature Scaling Relation

Probing Reionization with Quasar Spectra: the Impact of the Intrinsic Lyman-alpha Emission Line Shape Uncertainty

Radiative Feedback in Relic HII Regions at High-Redshift

Supermassive Black Hole Formation by Direct Collapse: Keeping Protogalactic Gas H_2--Free in Dark Matter Halos with Virial Temperatures T_vir >~ 10^4 K

Challenges to the DGP Model from Horizon-Scale Growth and Geometry

Observing the First Stars and Black Holes

The Evolution of the M_BH-sigma relation Inferred from the Age Distribution of Local Early-Type Galaxies and AGN Evolution

Constraints on Reionization and Source Properties from the Absorption Spectra of z>6.2 Quasars

Reconstructing the Cosmic Evolution of Quasars from the Age Distribution of Local Early-Type Galaxies

A Snapshot Survey for Gravitational Lenses Among z>=4.0 Quasars: II. Constraints on the 4.0<z<4.5 Quasar Population

Weighing Neutrinos with Galaxy Cluster Surveys

Rate Coefficient for H^+ + H_2(X ^1Sigma_g^+, nu=0, J=0) --> H(1s) + H_2^+ Charge Transfer and Some Cosmological Implications

The Formation and Evolution of the First Massive Black Holes

A Snapshot Survey for Gravitational Lenses Among z>=4.0 Quasars: I. The z>5.7 Sample

Do We Need Stars to Reionize the Universe at High Redshifts? Early Reionization by Decaying Heavy Sterile Neutrinos

Constraints on Omega_m, Omega_L, and Sigma_8, from Galaxy Cluster Redshift Distributions

What does the local black hole mass distribution tell us about the evolution of the quasar luminosity function?

Constraining the Lifetime of Quasars from their Spatial Clustering

Constraints on Cosmological Parameters from Future Galaxy Cluster Surveys