BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Zoltán Haiman

Zoltán 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.COastro-ph.HEastro-ph.GAgr-qcastro-phastro-ph.IMhep-phhep-th

Catalog footprint

What is connected

50works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

50 published item(s)

preprint2022arXiv

Cosmological Constraints from Weak Lensing Peaks: Can Halo Models Accurately Predict Peak Counts?

In order to extract full cosmological information from next-generation large and high-precision weak lensing (WL) surveys (e.g. Euclid, Roman, LSST), higher-order statistics that probe the small-scale, non-linear regime of large scale structure (LSS) need to be utilized. WL peak counts, which trace overdensities in the cosmic web, are one promising and simple statistic for constraining cosmological parameters. The physical origin of WL peaks have previously been linked to dark matter halos along the line of sight and this peak-halo connection has been used to develop various semi-analytic halo-based models for predicting peak counts. Here, we study the origin of WL peaks and the effectiveness of halo-based models for WL peak counts using a suite of ray-tracing N-body simulations. We compare WL peaks in convergence maps from the full simulations to those in maps created from only particles associated with halos -- the latter playing the role of a "perfect" halo model. We find that while halo-only contributions are able to replicate peak counts qualitatively well, halos do not explain all WL peaks. Halos particularly underpredict negative peaks, which are associated with local overdensities in large-scale underdense regions along the line of sight. In addition, neglecting non-halo contributions to peaks counts leads to a significant bias on the parameters ($Ω_{\rm m}$, $σ_{8}$) for surveys larger than $\geq$ 100 deg$^{2}$. We conclude that other elements of the cosmic web, outside and far away from dark matter halos, need to be incorporated into models of WL peaks in order to infer unbiased cosmological constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Long-term evolution of supercritical black hole accretion with outflows: a subgrid feedback model for cosmological simulations

We study the long-term evolution of the global structure of axisymmetric accretion flows onto a black hole (BH) at rates substantially higher than the Eddington value ($\dot{M}_{\rm Edd}$), performing two-dimensional hydrodynamical simulations with and without radiative diffusion. In the high-accretion optically-thick limit, where the radiation energy is efficiently trapped within the inflow, the accretion flow becomes adiabatic and comprises of turbulent gas in the equatorial region and strong bipolar outflows. As a result, the mass inflow rate decreases toward the center as $\dot{M}_{\rm in}\propto r^{p}$ with $p\sim 0.5-0.7$ and a small fraction of the inflowing gas feeds the nuclear BH. Thus, super-Eddington accretion is sustained only when a larger amount of gas is supplied from larger radii at $> 100-1000~\dot{M}_{\rm Edd}$. The global structure of the flow settles down to a quasi-steady state in millions of the orbital timescale at the BH event horizon, which is $> 10-100$ times longer than that addressed in previous (magneto-)RHD simulation studies. Energy transport via radiative diffusion accelerates the outflow near the poles in the inner region but does not change the overall properties of the accretion flow compared to the cases without diffusion. Based on our simulation results, we provide a mechanical feedback model for super-Eddington accreting BHs. This can be applied as a sub-grid model in large-scale cosmological simulations that do not sufficiently resolve galactic nuclei, and to the formation of the heaviest gravitational-wave sources via accretion in dense environments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Simultaneously constraining cosmology and baryonic physics via deep learning from weak lensing

Ongoing and planned weak lensing (WL) surveys are becoming deep enough to contain information on angular scales down to a few arcmin. To fully extract information from these small scales, we must capture non-Gaussian features in the cosmological WL signal while accurately accounting for baryonic effects. In this work, we account for baryonic physics via a baryonic correction model that modifies the matter distribution in dark matter-only $N$-body simulations, mimicking the effects of galaxy formation and feedback. We implement this model in a large suite of ray-tracing simulations, spanning a grid of cosmological models in $Ω_\mathrm{m}-σ_8$ space. We then develop a convolutional neural network (CNN) architecture to learn and constrain cosmological and baryonic parameters simultaneously from the simulated WL convergence maps. We find that in a Hyper-Suprime Cam (HSC)-like survey, our CNN achieves a 1.7$\times$ tighter constraint in $Ω_\mathrm{m}-σ_8$ space ($1σ$ area) than the power spectrum and 2.1$\times$ tighter than the peak counts, showing that the CNN can efficiently extract non-Gaussian cosmological information even while marginalizing over baryonic effects. When we combine our CNN with the power spectrum, the baryonic effects degrade the constraint in $Ω_\mathrm{m}-σ_8$ space by a factor of 2.4, compared to the much worse degradation by a factor of 4.7 or 3.7 from either method alone.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Fundamental Physics and Beyond the Standard Model

Gravitational wave detectors are formidable tools to explore strong-field gravity, especially black holes and neutron stars. These compact objects are extraordinarily efficient at producing electromagnetic and gravitational radiation. As such, they are ideal laboratories for fundamental physics and have an immense discovery potential. The detection of black hole binaries by third-generation Earth-based detectors, space-based detectors and pulsar timing arrays will provide exquisite tests of general relativity. Loud "golden" events and extreme mass-ratio inspirals can strengthen the observational evidence for horizons by mapping the exterior spacetime geometry, inform us on possible near-horizon modifications, and perhaps reveal a breakdown of Einstein's gravity. Measurements of the black-hole spin distribution and continuous gravitational-wave searches can turn black holes into efficient detectors of ultralight bosons across ten or more orders of magnitude in mass. A precise monitoring of the phase of inspiralling binaries can constrain the existence of additional propagating fields and characterize the environment in which the binaries live, bounding the local dark matter density and properties. Gravitational waves from compact binaries will probe general relativity and fundamental physics in previously inaccessible regimes, and allow us to address fundamental issues in our current understanding of the cosmos.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Supercritical growth pathway to overmassive black holes at cosmic dawn: coevolution with massive quasar hosts

Observations of the most luminous quasars at high redshifts ($z > 6$) have revealed that the largest supermassive black holes (SMBHs) at those epochs tend to be substantially overmassive relative to their host galaxies compared to the local relations, suggesting they experienced rapid early growth phases. We propose an assembly model for the SMBHs that end up in rare massive $\sim10^{12}~M_{\odot}$ host halos at $z \sim 6-7$, applying a kinetic feedback prescription for BHs accreting above the Eddington rate, provided by radiation hydrodynamic simulations for the long-term evolution of the accretion-flow structure. The large inflow rates into these halos during their assembly enable the formation of $>10^9~M_{\odot}$ SMBHs by $z \sim 6$, even starting from stellar-mass seeds at $z \sim 30$, and even in the presence of outflows that reduce the BH feeding rate, especially at early times. This mechanism also naturally yields a high BH-to-galaxy mass ratio of $> 0.01$ before the SMBH mass reaches $M_{\rm BH} > 10^9~M_{\odot}$ by $z \sim 6$. These fast-growing SMBH progenitors are bright enough to be detected by upcoming observations with the James Webb Space Telescope over a wide range of redshift ($7 < z < 15$), regardless of how they were seeded.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The science case and challenges of space-borne sub-millimeter interferometry

Ultra-high angular resolution in astronomy has always been an important vehicle for making fundamental discoveries. Recent results in direct imaging of the vicinity of the supermassive black hole in the nucleus of the radio galaxy M87 by the millimeter VLBI system Event Horizon Telescope and various pioneering results of the Space VLBI mission RadioAstron provided new momentum in high angular resolution astrophysics. In both mentioned cases, the angular resolution reached the values of about 10-20 microrcseconds. Further developments toward at least an order of magnitude "sharper" values are dictated by the needs of astrophysical studies and can only be achieved by placing millimeter and submillimeter wavelength interferometric systems in space. A concept of such the system, called Terahertz Exploration and Zooming-in for Astrophysics (THEZA), has been proposed in the framework of the ESA Call for White Papers for the Voayage 2050 long term plan in 2019. In the current paper we discuss several approaches for addressing technological challenges of the THEZA concept. In particular, we consider a novel configuration of a space-borne millimeter/sub-millimeter antenna which might resolve several bottlenecks in creating large precise mechanical structures. The paper also presents an overview of prospective space-qualified technologies of low-noise analogue front-end instrumentation for millimeter/sub-millimeter telescopes, data handling and processing. The paper briefly discusses approaches to the interferometric baseline state vector determination and synchronisation and heterodyning system. In combination with the original ESA Voyage 2050 White Paper, the current work sharpens the case for the next generation microarcsceond-level imaging instruments and provides starting points for further in-depth technology trade-off studies.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

$κ$TNG: Effect of Baryonic Processes on Weak Lensing with IllustrisTNG Simulations

We study the effect of baryonic processes on weak lensing (WL) observables with a suite of mock WL maps, the $κ$TNG, based on the cosmological hydrodynamic simulations IllustrisTNG. We quantify the baryonic effects on the WL angular power spectrum, one-point probability distribution function (PDF), and number counts of peaks and minima. We also show the redshift evolution of the effects, which is a key to distinguish the effect of baryons from fundamental physics such as dark energy, dark matter, and massive neutrinos. We find that baryonic processes reduce the small-scale power, suppress the tails of the PDF, peak and minimum counts, and change the total number of peaks and minima. We compare our results to existing semi-analytic models and hydrodynamic simulations, and discuss the source of discrepancies. The $κ$TNG suite includes 10,000 realisations of $5 \times 5 \, \mathrm{deg}^2$ maps for 40 source redshifts up to $z_s = 2.6$, well covering the range of interest for existing and upcoming weak lensing surveys. We also produce the $κ$TNG-Dark suite of maps, generated based on the corresponding dark matter only IllustrisTNG simulations. Our mock maps are suitable for developing analytic models that incorporate the effect of baryons, but also particularly useful for studies that rely on mass maps, such as non-Gaussian statistics and machine learning with convolutional neural networks. The suite of mock maps is publicly available at Columbia Lensing (http://columbialensing.org).

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Can stellar-mass black hole growth disrupt disks of active galactic nuclei? The role of mechanical feedback

Stellar-mass BHs (sBHs) are predicted to be embedded in active galactic nuclei (AGN) disks due to gravitational drag and in-situ star formation. However, we find that due to a high gas density in an AGN disk environment, compact objects may rapidly grow to intermediate-mass BHs and deplete matter from the AGN disk unless accretion is suppressed by some feedback process(es). These consequences are inconsistent with AGN observations and the dynamics of the Galactic center. Here we consider mechanical feedback mechanisms for the reduction of gas accretion. Rapidly accreting sBHs launch winds and/or jets via the Blandford-Znajek mechanism, which produce high-pressure shocks and cocoons. Such a shock and cocoon can spread laterally in the plane of the disk, eject the outer regions of a circum-sBH disk (CsBD) and puncture a hole in the AGN disk with horizontal size comparable to the disk scale-height. Since the depletion timescale of the bound CsBD is much shorter than the resupply timescale of gas to the sBH, the time-averaged accretion rate onto sBHs is reduced by this process by a factor of $\sim 10$--$100$. This feedback mechanism can therefore help alleviate the sBH over-growth and AGN-disk depletion problems. On the other hand, we find that cocoons of jets can unbind a large fraction of the gas accreting in the disks of less massive SMBHs, which may help explain the dearth of high-Eddington ratio AGNs with SMBH mass $\lesssim10^5{\rm M_\odot}$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Equilibrium eccentricity of accreting binaries

Using high-resolution hydrodynamics simulations, we show that equal-mass binaries accreting from a circumbinary disk evolve toward an orbital eccentricity of $e \simeq 0.45$, unless they are initialized on a nearly circular orbit with $e \lesssim 0.08$, in which case they further circularize. The implied bi-modal eccentricity distribution resembles that seen in post-AGB stellar binaries. Large accretion spikes around periapse impart a tell-tale, quasi-periodic, bursty signature on the light curves of eccentric binaries. We predict that intermediate-mass and massive black hole binaries at $z \lesssim 10$ entering the \emph{LISA} band will have measurable eccentricities in the range $e \simeq 10^{-3} - 10^{-2}$, if they have experienced a gas-driven phase. On the other hand, GW190521 would have entered the LIGO/Virgo band with undetectable eccentricity $\sim 10^{-6}$ if it had been driven into the gravitational wave regime by a gas disk.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Self-lensing flares from black hole binaries I: general-relativistic ray tracing of black hole binaries

The self-lensing of a massive black hole binary (MBHB), which occurs when the two BHs are aligned close to the line of sight, is expected to produce periodic, short-duration flares. Here we study the shapes of self-lensing flares (SLFs) via general-relativistic ray tracing in a superimposed binary BH metric, in which the emission is generated by geometrically thin accretion flows around each component. The suite of models covers eccentric binary orbits, black hole spins, unequal mass binaries, and different emission model geometries. We explore the above parameter space, and report how the light curves change as a function of, e.g., binary separation, inclination, and eccentricity. We also compare our light curves to those in the microlensing approximation, and show how strong deflections, as well as time-delay effects, change the size and shape of the SLF. If gravitational waves (GWs) from the inspiraling MBHB are observed by LISA, SLFs can help securely identify the source and localizing it on the sky, and to constrain the graviton mass by comparing the phasing of the SLFs and the GWs. Additionally, when these systems are viewed edge-on the SLF shows a distinct dip that can be directly correlated with the BH shadow size. This opens a new way to measure BH shadow sizes in systems that are unresolvable by current VLBI facilities.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Self-lensing flares from black hole binaries II: observing black hole shadows via light-curve tomography

Supermassive black hole (BH) binaries are thought to produce self-lensing flares (SLF) when the two BHs are aligned with the line-of-sight. If the binary orbit is observed nearly edge-on, we find a distinct feature in the light curve imprinted by the BH shadow from the lensed BH. We study this feature by ray-tracing in a binary model and predict that 1\% of the current binary candidates could show this feature. Our BH tomography method proposed here could make it possible to extract BH shadows that are spatially unresolvable by high-resolution VLBI.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Chandra Observations of Candidate Sub-Parsec Binary Supermassive Black Holes

We present analysis of Chandra X-ray observations of seven quasars that were identified as candidate sub-parsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey (CRTS) based on apparent periodicity in their optical light curves. Simulations predict close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs, including harder or softer X-ray spectra, ripple-like profiles in the Fe K-$α$ line, and distinct peaks in the spectrum due to the separation of the accretion disk into a circumbinary disk and mini-disks around each SMBH. We obtained Chandra observations to test these models and assess whether these quasars could contain binary SMBHs. We instead find that the quasar spectra are all well fit by simple absorbed power law models, with the rest frame 2-10 keV photon indices, $Γ$, and the X-ray-to-optical power slopes, $α_{\rm OX}$, indistinguishable from the larger quasar population. This may indicate that these seven quasars are not truly sub-parsec binary SMBH systems, or it may simply reflect that our sample size was too small to robustly detect any differences. Alternatively, the X-ray spectral changes might only be evident at higher energies than probed by Chandra. Given the available models and current data, no firm conclusions are drawn. These observations will help motivate and direct further work on theoretical models of binary SMBH systems, such as modeling systems with thinner accretion disks and larger binary separations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Correlation between Optical and UV Variability of Quasars

The variability of quasars across multiple wavelengths is a useful probe of physical conditions in active galactic nuclei. In particular, variable accretion rates, instabilities, and reverberation effects in the accretion disk of a supermassive black hole (SMBH) are expected to produce correlated flux variations in UV and optical bands. Recent work has further argued that binary quasars should exhibit strongly correlated UV and optical periodicities. Strong UV-optical correlations have indeed been established in small samples of up to approximately 30 quasars with well-sampled light curves, and have extended the "bluer-when-brighter" trend previously found within the optical bands. Here we further test the nature of quasar variability by examining the observed-frame UV-optical correlations in a large sample of 1,315 bright quasars with overlapping UV and optical light curves for the Galaxy Evolution Explorer (GALEX) and the Catalina Real-time Transient Survey (CRTS), respectively. We find that strong correlations exist in this much larger sample, but we rule out, at approximately 95% confidence, the simple hypothesis that the intrinsic UV and optical variations of all quasars are fully correlated. Our results therefore imply the existence of physical mechanism(s) that can generate uncorrelated optical and UV flux variations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Massive Star Formation in Metal-Enriched Haloes at High Redshift

The formation of supermassive stars has generally been studied under the assumption of rapid accretion of pristine metal-free gas. Recently it was found, however, that gas enriched to metallicities up to $Z \sim 10^{-3}$ Z$_{\odot}$ can also facilitate supermassive star formation, as long as the total mass infall rate onto the protostar remains sufficiently high. We extend the analysis further by examining how the abundance of supermassive star candidate haloes would be affected if all haloes with super-critical infall rates, regardless of metallicity were included. We investigate this scenario by identifying all atomic cooling haloes in the Renaissance simulations with central mass infall rates exceeding a fixed threshold. We find that among these haloes with central mass infall rates above 0.1 M$_{\odot}$ yr$^{-1}$ approximately two-thirds of these haloes have metallicities of $Z > 10^{-3}$ Z$_{\odot}$. If metal mixing within these haloes is inefficient early in their assembly and pockets of metal-poor gas can remain then the number of haloes hosting supermassive stars can be increased by at least a factor of four. Additionally the centres of these high infall-rate haloes provide ideal environments in which to grow pre-existing black holes. Further research into the (supermassive) star formation dynamics of rapidly collapsing haloes, with inhomogeneous metal distributions, is required to gain more insight into both supermassive star formation in early galaxies as well as early black hole growth.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Probing gaseous galactic halos through the rotational kSZ effect

The rotational kinematic Sunyaev-Zeldovich (rkSZ) signal, imprinted on the cosmic microwave background (CMB) by the gaseous halos (spinning "atmospheres") of foreground galaxies, would be a novel probe of galaxy formation. Although the signal is too weak to detect in individual galaxies, we analyze the feasibility of its statistical detection via stacking CMB data on many galaxies for which the spin orientation can be estimated spectroscopically. We use an "optimistic" model, in which fully ionized atmospheres contain the cosmic baryon fraction and spin at the halo's circular velocity $v_{\rm circ}$, and a more realistic model, based on hydrodynamical simulations, with multi-phase atmospheres spinning at a fraction of $v_{\rm circ}$. We incorporate realistic noise estimates into our analysis. Using low-redshift galaxy properties from the MaNGA spectroscopic survey (with median halo mass of $6.6\times10^{11}\,M_{\odot}$), and CMB data quality from Planck, we find that a $3σ$ detection would require a few$\times 10^4$ galaxies, even in the optimistic model. This is too high for current surveys, but upcoming higher-angular resolution CMB experiments will significantly reduce the requirements: stacking CMB data on galaxy spins in a $\sim10$ deg$^2$ can rule out the optimistic models, and $\approx$350 deg$^2$ will suffice for a $3σ$ detection with ACT. As a proof-of-concept, we stacked Planck data on the position of $\approx2,000$ MaNGA galaxies, aligned with the galaxies' projected spin, and scaled to their halos' angular size. We rule out average temperature dipoles larger than $\approx1.9\,μ$K around field spiral galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

CMB Lensing Beyond the Power Spectrum: Cosmological Constraints from the One-Point PDF and Peak Counts

Unprecedentedly precise cosmic microwave background (CMB) data are expected from ongoing and near-future CMB Stage-III and IV surveys, which will yield reconstructed CMB lensing maps with effective resolution approaching several arcminutes. The small-scale CMB lensing fluctuations receive non-negligible contributions from nonlinear structure in the late-time density field. These fluctuations are not fully characterized by traditional two-point statistics, such as the power spectrum. Here, we use $N$-body ray-tracing simulations of CMB lensing maps to examine two higher-order statistics: the lensing convergence one-point probability distribution function (PDF) and peak counts. We show that these statistics contain significant information not captured by the two-point function, and provide specific forecasts for the ongoing Stage-III Advanced Atacama Cosmology Telescope (AdvACT) experiment. Considering only the temperature-based reconstruction estimator, we forecast 9$σ$ (PDF) and 6$σ$ (peaks) detections of these statistics with AdvACT. Our simulation pipeline fully accounts for the non-Gaussianity of the lensing reconstruction noise, which is significant and cannot be neglected. Combining the power spectrum, PDF, and peak counts for AdvACT will tighten cosmological constraints in the $Ω_m$-$σ_8$ plane by $\approx 30\%$, compared to using the power spectrum alone.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cosmology with photometric weak lensing surveys: constraints with redshift tomography of convergence peaks and moments

Weak gravitational lensing is becoming a mature technique for constraining cosmological parameters, and future surveys will be able to constrain the dark energy equation of state $w$. When analyzing galaxy surveys, redshift information has proven to be a valuable addition to angular shear correlations. We forecast parameter constraints on the triplet $(Ω_m,w,σ_8)$ for an LSST-like photometric galaxy survey, using tomography of the shear-shear power spectrum, convergence peak counts and higher convergence moments. We find that redshift tomography with the power spectrum reduces the area of the $1σ$ confidence interval in $(Ω_m,w)$ space by a factor of 8 with respect to the case of the single highest redshift bin. We also find that adding non-Gaussian information from the peak counts and higher-order moments of the convergence field and its spatial derivatives further reduces the constrained area in $(Ω_m,w)$ by a factor of 3 and 4, respectively. When we add cosmic microwave background parameter priors from Planck to our analysis, tomography improves power spectrum constraints by a factor of 3. Adding moments yields an improvement by an additional factor of 2, and adding both moments and peaks improves by almost a factor of 3, over power spectrum tomography alone. We evaluate the effect of uncorrected systematic photometric redshift errors on the parameter constraints. We find that different statistics lead to different bias directions in parameter space, suggesting the possibility of eliminating this bias via self-calibration.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Do dark matter halos explain lensing peaks?

We have investigated a recently proposed halo-based model, Camelus, for predicting weak-lensing peak counts, and compared its results over a collection of 162 cosmologies with those from N-body simulations. While counts from both models agree for peaks with $\mathcal{S/N}>1$ (where $\mathcal{S/N}$ is the ratio of the peak height to the r.m.s. shape noise), we find $\approx 50\%$ fewer counts for peaks near $\mathcal{S/N}=0$ and significantly higher counts in the negative $\mathcal{S/N}$ tail. Adding shape noise reduces the differences to within $20\%$ for all cosmologies. We also found larger covariances that are more sensitive to cosmological parameters. As a result, credibility regions in the $\{Ω_m, σ_8\}$ are $\approx 30\%$ larger. Even though the credible contours are commensurate, each model draws its predictive power from different types of peaks. Low peaks, especially those with $2<\mathcal{S/N}<3$, convey important cosmological information in N-body data, as shown in \cite{DietrichHartlap, Kratochvil2010}, but \textsc{Camelus} constrains cosmology almost exclusively from high significance peaks $(\mathcal{S/N}>3)$. Our results confirm the importance of using a cosmology-dependent covariance with at least a 14\% improvement in parameter constraints. We identified the covariance estimation as the main driver behind differences in inference, and suggest possible ways to make Camelus even more useful as a highly accurate peak count emulator.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hyper-Eddington mass accretion onto a black hole with super-Eddington luminosity

We perform one-dimensional radiation hydrodynamical simulations to solve accretion flows onto massive black holes (BHs) with a very high rate. Assuming that photon trapping limits the luminosity emerging from the central region to $L\lesssim L_{\rm Edd}$, Inayoshi, Haiman & Ostriker (2016) have shown that an accretion flow settles to a "hyper-Eddington" solution, with a steady and isothermal ($T\simeq 8000$ K) Bondi profile reaching $\gtrsim 5000$ times the Eddington accretion rate $\dot{M}_{\rm Edd}\equiv L_{\rm Edd}/c^2$. Here we address the possibility that gas accreting with finite angular momentum forms a bright nuclear accretion disc, with a luminosity exceeding the Eddington limit ($1\lesssim L/L_{\rm Edd}\lesssim 100$). Combining our simulations with an analytic model, we find that a transition to steady hyper-Eddington accretion still occurs, as long as the luminosity remains below $L/L_{\rm Edd}\lesssim 35~(M_{\rm BH}/10^4~M_\odot)^{3/2}(n_\infty/10^5~{\rm cm^{-3}})(T_\infty/10^4~{\rm K})^{-3/2}(r_{\star}/10^{14}~{\rm cm})^{-1/2}$, where $n_\infty$ and $T_\infty$ are the density and temperature of the ambient gas, and $r_\star$ is the radius of the photosphere, at which radiation emerges. If the luminosity exceeds this value, accretion becomes episodic. Our results can be accurately recovered in a toy model of an optically thick spherical shell, driven by radiation force into a collapsing medium. When the central source is dimmer than the above critical value, the expansion of the shell is halted and reversed by ram pressure of the collapsing medium, and by shell's weight. Our results imply that rapid, unimpeded hyper-Eddington accretion is possible even if the luminosity of the central source far exceeds the Eddington limit, and can be either steady or strongly episodic.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Sample variance in weak lensing: how many simulations are required?

Constraining cosmology using weak gravitational lensing consists of comparing a measured feature vector of dimension $N_b$ with its simulated counterpart. An accurate estimate of the $N_b\times N_b$ feature covariance matrix $\mathbf{C}$ is essential to obtain accurate parameter confidence intervals. When $\mathbf{C}$ is measured from a set of simulations, an important question is how large this set should be. To answer this question, we construct different ensembles of $N_r$ realizations of the shear field, using a common randomization procedure that recycles the outputs from a smaller number $N_s\leq N_r$ of independent ray-tracing $N$--body simulations. We study parameter confidence intervals as a function of ($N_s,N_r$) in the range $1\leq N_s\leq 200$ and $1\leq N_r\lesssim 10^5$. Previous work has shown that Gaussian noise in the feature vectors (from which the covariance is estimated) lead, at quadratic order, to an $O(1/N_r)$ degradation of the parameter confidence intervals. Using a variety of lensing features measured in our simulations, including shear-shear power spectra and peak counts, we show that cubic and quartic covariance fluctuations lead to additional $O(1/N_r^2)$ error degradation that is not negligible when $N_r$ is only a factor of few larger than $N_b$. We study the large $N_r$ limit, and find that a single, 240Mpc$/h$ sized $512^3$-particle $N$--body simulation ($N_s=1$) can be repeatedly recycled to produce as many as $N_r={\rm few}\times10^4$ shear maps whose power spectra and high-significance peak counts can be treated as statistically independent. As a result, a small number of simulations ($N_s=1$ or $2$) is sufficient to forecast parameter confidence intervals at percent accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A reduced orbital period for the supermassive black hole binary candidate in the quasar PG 1302-102?

Graham et al. (2015) have detected a 5.2 year periodic optical variability of the quasar PG 1302-102 at redshift $z=0.3$, which they interpret as the redshifted orbital period $(1+z)t_{\rm bin}$ of a putative supermassive black hole binary (SMBHB). Here we consider the implications of a $3-8$ times shorter orbital period, suggested by hydrodynamical simulations of circumbinary discs (CBDs) with nearly equal--mass SMBHBs ($q\equiv M_2/M_1\gtrsim 0.3$). With the corresponding $2-4$ times tighter binary separation, PG 1302 would be undergoing gravitational wave dominated inspiral, and serve as a proof that the BHs can be fueled and produce bright emission even in this late stage of the merger. The expected fraction of binaries with the shorter $t_{\rm bin}$, among bright quasars, would be reduced by 1-2 orders of magnitude, compared to the 5.2 year period, in better agreement with the rarity of candidates reported by Graham et al. (2015). Finally, shorter periods would imply higher binary speeds, possibly imprinting periodicity on the light curves from relativistic beaming, as well as measurable relativistic effects on the Fe K$α$ line. The CBD model predicts additional periodic variability on timescales of $t_{\rm bin}$ and $\approx 0.5 t_{\rm bin}$, as well as periodic variation of broad line widths and offsets relative to the narrow lines, which are consistent with the observations. Future observations will be able to test these predictions and hence the binary+CBD hypothesis for PG 1302.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cosmic Reionization and the First Nonlinear Structures in the Universe

In this introductory chapter, we outline expectations for when and how the hydrogen and helium atoms in the universe turned from neutral to ionized, focusing on the earliest, least well understood stages, and emphasizing the most important open questions. We include a historical summary, and highlight the role of reionization as one of the few milestones in the evolution of the universe since the Big Bang, and its status as a unique probe of the beginning stages of structure formation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multiple periods in the variability of the supermassive black hole binary candidate quasar PG1302-102?

Graham et al. (2015) discovered a supermassive black hole binary (SMBHB) candidate and identified the detected 5.2 yr period of the optical variability as the orbital period of the binary. Hydrodynamical simulations predict multiple periodic components for the variability of SMBHBs, thus raising the possibility that the true period of the binary is different from 5.2 yr. We analyse the periodogram of PG1302 and find no compelling evidence for additional peaks. We derive upper limits on any additional periodic modulations in the available data, by modeling the light-curve as the sum of red noise and the known 5.2 yr periodic component, and injecting additional sinusoidal signals. We find that, with the current data, we would be able to detect with high significance (false alarm probability <1%) secondary periodic terms with periods in the range predicted by the simulations, if the amplitude of the variability was at least ~0.07 mag (compared to 0.14 mag for the main peak). A three-year follow-up monitoring campaign with weekly observations can increase the sensitivity for detecting secondary peaks ~50%, and would allow a more robust test of predictions from hydrodynamical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Relativistic boost as the cause of periodicity in a massive black-hole binary candidate

Because most large galaxies contain a central black hole, and galaxies often merge, black-hole binaries are expected to be common in galactic nuclei. Although they cannot be imaged, periodicities in the light curves of quasars have been interpreted as evidence for binaries, most recently in PG~1302-102, with a short rest-frame optical period of 4 years. If the orbital period matches this value, then for the range of estimated black hole masses the components would be separated by 0.007-0.017 pc, implying relativistic orbital speeds. There has been much debate over whether black hole orbits could be smaller than 1 pc. Here we show that the amplitude and the sinusoid-like shape of the variability of PG~1302-102 can be fit by relativistic Doppler boosting of emission from a compact, steadily accreting, unequal-mass binary. We predict that brightness variations in the ultraviolet light curve track those in the optical, but with a 2-3 times larger amplitude. This prediction is relatively insensitive to the details of the emission process, and is consistent with archival UV data. Follow-up UV and optical observations in the next few years can test this prediction and confirm the existence of a binary black hole in the relativistic regime.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cosmology with Minkowski functionals and moments of the weak lensing convergence field

We compare the efficiency of moments and Minkowski functionals (MFs) in constraining the subset of cosmological parameters (Omega_m,w,sigma_8) using simulated weak lensing convergence maps. We study an analytic perturbative expansion of the MFs in terms of the moments of the convergence field and of its spatial derivatives. We show that this perturbation series breaks down on smoothing scales below 5', while it shows a good degree of convergence on larger scales (15'). Most of the cosmological distinguishing power is lost when the maps are smoothed on these larger scales. We also show that, on scales comparable to 1', where the perturbation series does not converge, cosmological constraints obtained from the MFs are approximately 1.5-2 times better than the ones obtained from the first few moments of the convergence distribution --- provided that the latter include spatial information, either from moments of gradients, or by combining multiple smoothing scales. Including either a set of these moments or the MFs can significantly tighten constraints on cosmological parameters, compared to the conventional method of using the power spectrum alone.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The origin of spin in galaxies: clues from simulations of atomic cooling halos

In order to elucidate the origin of spin in both dark matter and baryons in galaxies, we have performed hydrodynamical simulations from cosmological initial conditions. We study atomic cooling haloes in the redshift range $100 > z > 9$ with masses of order $10^9{\rm M_{\odot}}$ at redshift $z=10$. We assume that the gas has primordial composition and that ${\rm H_2}$-cooling and prior star-formation in the haloes have been suppressed. We present a comprehensive analysis of the gas and dark matter properties of four halos with very low ($λ\approx 0.01$), low ($λ\approx 0.04$), high ($λ\approx 0.06$) and very high ($λ\approx 0.1$) spin parameter. Our main conclusion is that the spin orientation and magnitude is initially well described by tidal torque linear theory, but later on is determined by the merging and accretion history of each halo. We provide evidence that the topology of the merging region, i.e. the number of colliding filaments, gives an accurate prediction for the spin of dark matter and gas: halos at the center of knots will have low spin while those in the center of filaments will have high spin. The spin of a halo is given by $λ\approx 0.05 \times \left(\frac{7.6}{\rm number\,\,\, of \,\,\, filaments}\right)^{5.1}$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Accretion into the Central Cavity of a Circumbinary Disk

A near-equal mass binary black hole can clear a central cavity in a circumbinary accretion disk; however, previous works have revealed accretion streams entering this cavity. Here we use 2D hydrodynamical simulations to study the accretion streams and their periodic behavior. In particular, we perform a suite of simulations, covering different binary mass ratios $q=M_2/M_1$ in the range $0.01 \leq q \leq 1$. In each case, we follow the system for several thousand binary orbits, until it relaxes to a stable accretion pattern. We find the following results: (i) while the binary is efficient in maintaining a low-density cavity, the time-averaged mass accretion rate into the cavity, through narrow coherent accretion streams, is suppressed by at most a factor of $\sim 5$, compared to a disk with a single BH with the same mass; (ii) the largest suppression occurs for $q\approx 0.05$; binaries whose mass ratios are either lower or higher both suppress accretion less significantly; (iii) for $q \gtrsim 0.05$, the accretion rate is strongly modulated by the binary, and depending on the precise value of $q$, the power spectrum of the accretion rate shows either one, two, or three distinct periods; and (v) for $q \lesssim 0.05$, the accretion rate becomes steady, with no time-variations. Most binaries produced in galactic mergers are expected to have $q\gtrsim 0.05$. If the luminosity of these binaries tracks their accretion rate, then a periodogram of their light-curve could help in their identification, and to constrain their mass ratio and disk properties.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Electromagnetic signatures of supermassive black hole binaries resolved by PTAs

Pulsar timing arrays (PTAs) may eventually be able to detect not only the stochastic gravitational-wave (GW) background of SMBH binaries, but also individual, particularly massive binaries whose signals stick out above the background. In this contribution, we discuss the possibility of identifying and studying such `resolved' binaries through their electromagnetic emission. The host galaxies of such binaries are themselves expected to be also very massive and rare, so that out to redshifts z~2 a unique massive galaxy may be identified as the host. At higher redshifts, the PTA error boxes are larger and may contain as many as several hundred massive-galaxy interlopers. In this case, the true counterpart may be identified, if it is accreting gas efficiently, as an active galactic nucleus (AGN) with a peculiar spectrum and variable emission features. Specifically, the binary's tidal torques expel the gas from the inner part of the accretion disk, making it unusually dim in X-ray and UV bands and in broad optical emission lines. The tails of the broad wings of any FeK-alpha emission line may also be `clipped' and missing. The binary's orbital motion, as well as the gas motions it induces, may trigger quasiperiodic variations. These include coherent flux variability, such as luminous, multi-wavelength flares, as well as Doppler shifts of broad emission lines and `see-saw' oscillations in the FeK-alpha line. Additional features, such as evidence for a recent major merger or dual collimated jets, could also corroborate the counterpart. These properties would make resolved PTA sources stand out among AGN with similar overall luminosities and allow their identification.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The effect of baryonic streaming motions on the formation of the first supermassive black holes

Observations of quasars at redshifts z > 6 reveal that 10^9 Msol supermassive black holes (SMBHs) had already formed when the Universe was < 0.9 Gyr old. One hypothesis for the origins of these SMBHs is that they grew from the remnants of the first generation of massive stars, which formed in low-mass (~ 10^5 to 10^6 Msol) dark matter minihaloes at z > 20. This is the regime where baryonic streaming motions--the relative velocities of baryons with respect to dark matter in the early Universe--most strongly inhibit star formation by suppressing gas infall and cooling. We investigate the impact of this effect on the growth of the first SMBHs using a suite of high-fidelity, ellipsoidal-collapse Monte Carlo merger-tree simulations. We find that the suppression of seed BH formation by the streaming motions significantly reduces the number density of the most massive BHs at z > 15, but the residual effect at lower redshifts is essentially negligible. The streaming motions can reduce by a factor of few the number density of the most luminous quasars at z ~ 10-11, where such objects could be detected by the James Webb Space Telescope. We conclude, with minor theoretical caveats, that baryonic streaming motions are unlikely to pose a significant additional obstacle to the formation of the observed high-redshift quasar SMBHs. Nor do they appreciably affect the heating and reionization histories of the Universe or the merger rates of nuclear BHs in the mass and redshift ranges of interest for proposed gravitational-wave detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Impact of Magnification and Size Bias on Weak Lensing Power Spectrum and Peak Statistics

The weak lensing power spectrum is a powerful tool to probe cosmological parameters. Additionally, lensing peak counts contain cosmological information beyond the power spectrum. Both of these statistics can be affected by the preferential selection of source galaxies in patches of the sky with high magnification, as well as by the dilution in the source galaxy surface density in such regions. If not accounted for, these biases introduce systematic errors for cosmological measurements. Here we quantify these systematic errors, using convergence maps from a suite of ray-tracing N-body simulations. At the cut-off magnitude m of on-going and planned major weak lensing surveys, the logarithmic slope of the cumulative number counts s = dlog[n(>m)]/dlog(m) is in the range 0.1 < s < 0.5. At s = 0.2, expected in the I band for LSST, the inferred values of Omega_m, w and sigma_8 are biased by many sigma (where sigma denotes the marginalized error) and therefore the biases will need to be carefully modeled. We also find that the parameters are biased differently in the (Omega_m, w, sigma_8) parameter space when the power spectrum and when the peak counts are used. In particular, w derived from the power spectrum is less affected than w derived from peak counts, while the opposite is true for the best-constrained combination of [sigma_8 Omega_m^gamma] (with gamma=0.62 from the power spectrum and gamma = 0.48 from peak counts). This suggests that the combination of the power spectrum and peak counts can help mitigate the impact of magnification and size biases.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Baryon impact on weak lensing peaks and power spectrum: low-bias statistics and self-calibration in future surveys

Peaks in two-dimensional weak lensing (WL) maps contain significant cosmological information, complementary to the WL power spectrum. This has recently been demonstrated using N-body simulations which neglect baryonic effects. Here we employ ray-tracing N-body simulations in which we manually steepen the density profile of each dark matter halo, mimicking the cooling and concentration of baryons into dark matter potential wells. We find, in agreement with previous works, that this causes a significant increase in the amplitude of the WL power spectrum on small scales (spherical harmonic index l>1,000). We then study the impact of the halo concentration increase on the peak counts, and find the following. (i) Low peaks (with convergence 0.02 < kappa_peak < 0.08), remain nearly unaffected. These peaks are created by a constellation of several halos with low masses (10^12-10^13 M_sun) and large angular offsets from the peak center (> 0.5 R_vir); as a result, they are insensitive to the central halo density profiles. These peaks contain most of the cosmological information, and thus provide an unusually sensitive and unbiased probe. (ii) The number of high peaks (with convergence kappa_peak > 0.08) is increased. However, when the baryon effects are neglected in cosmological parameter estimation, then the high peaks lead to a modest bias, comparable to that from the power spectrum on relatively large-scales (l<2000), and much smaller than the bias from the power spectrum on smaller scales (l>2,000). (iii) In the 3D parameter space (sigma_8, Omega_m, w), the biases from the high peaks and the power spectra are in different directions. This suggests the possibility of "self-calibration": the combination of peak counts and power spectrum can simultaneously constrain baryonic physics and cosmological parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Feedback from the IR Background in the Early Universe

It is commonly believed that the earliest stages of star-formation in the Universe were self-regulated by global radiation backgrounds - either by the ultraviolet Lyman-Werner (LW) photons emitted by the first stars (directly photodissociating H_2), or by the X-rays produced by accretion onto the black hole (BH) remnants of these stars (heating the gas but catalyzing H_2 formation). Recent studies have suggested that a significant fraction of the first stars may have had low masses (a few M_sun). Such stars do not leave BH remnants and they have softer spectra, with copious infrared (IR) radiation at photon energies around 1eV. Similar to LW and X-ray photons, these photons have a mean-free path comparable to the Hubble distance, building up an early IR background. Here we show that if soft-spectrum stars, with masses of a few M_sun, contributed more than 1% of the UV background (or their mass fraction exceeded 90%), then their IR radiation dominated radiative feedback in the early Universe. The feedback is different from the UV feedback from high-mass stars, and occurs through the photo-detachment of H^- ions, necessary for efficient H_2 formation. Nevertheless, we find that the baryon fraction which must be incorporated into low-mass stars in order to suppress H_2-cooling is only a factor of few higher than for high-mass stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Formation of the First Massive Black Holes

Supermassive black holes (SMBHs) are common in local galactic nuclei, and SMBHs as massive as several billion solar masses already exist at redshift z=6. These earliest SMBHs may grow by the combination of radiation-pressure-limited accretion and mergers of stellar-mass seed BHs, left behind by the first generation of metal-free stars, or may be formed by more rapid direct collapse of gas in rare special environments where dense gas can accumulate without first fragmenting into stars. This chapter offers a review of these two competing scenarios, as well as some more exotic alternative ideas. It also briefly discusses how the different models may be distinguished in the future by observations with JWST, (e)LISA and other instruments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

X-ray emission from high-redshift miniquasars: self-regulating the population of massive black holes through global warming

Observations of high-redshift quasars at z>6 imply that supermassive black holes (SMBHs) with masses over a billion solar masses were in place less than 1 Gyr after the Big Bang. If these SMBHs assembled from "seed" BHs left behind by the first stars, then they must have accreted gas at close to the Eddington limit during a large fraction (>50%) of the time. A generic problem with this scenario, however, is that the mass density in million-solar-mass SMBHs at z=6 already exceeds the locally observed SMBH mass density by several orders of magnitude; in order to avoid this overproduction, BH seed formation and growth must become significantly less efficient in less massive protogalaxies, while proceeding uninterrupted in the most massive galaxies that formed first. Using Monte-Carlo realizations of the merger and growth history of BHs, we show that X-rays from the earliest accreting BHs can provide such a feedback mechanism. Our calculations paint a self-consistent picture of black-hole-made climate change, in which the first miniquasars---among them the ancestors of the z>6 quasar SMBHs---globally warm the IGM and suppress the formation and growth of subsequent generations of BHs. We present two specific models with global miniquasar feedback that provide excellent agreement with recent estimates of the z=6 SMBH mass function. For each of these models, we estimate the rate of BH mergers at z>6 that could be detected by the proposed gravitational-wave observatory eLISA/NGO.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Delayed Enrichment by Unseen Galaxies: Explaining the Rapid Rise in IGM CIV Absorption from z = 6-5

In the near future, measurements of metal absorption features in the intergalactic medium (IGM) will become an important constraint on models of the formation and evolution of the earliest galaxies, the properties of the first stars, and the reionization and enrichment of the IGM. The first measurement of a metal abundance in the IGM at a redshift approaching the epoch of reionization already offers intriguing hints. Between z=5.8 and 4.7 (a 0.3 Gyr interval only 1 Gyr after the big bang), the measured density of CIV absorbers in the IGM increased by a factor of ~ 3.5 (Ryan-Weber et al. 2009; Becker, Rauch & Sargent 2009). If these values prove to be accurate, they pose two puzzles: (1) The total amount of CIV at z=5.8 implies too little star formation to reionize the IGM by z=6 or to match the WMAP electron scattering optical depth (tau). (2) The rapid growth from z = 6-5 is faster than the buildup of stellar mass or the increase in the star formation rate density over the same interval. We show that a delay of ~ 0.4-0.7 Gyr between the instantaneous production of ionizing photons and the later production of metal absorption features (added to the delay due to stellar lifetimes) can provide the full explanation for both puzzles. We calculate the delay in metal production due to finite stellar lifetimes alone and find that it is too short to explain the rapid CIV density increase. The additional delay could naturally be explained as the result of ~ 200 km/s outflows carrying carbon to distances of ~ 100 kpc, the typical distance between galaxies and CIV absorbers in enrichment simulations, and the typical outflow or absorption region scale observed at z ~ 2-3.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Electromagnetic counterparts of supermassive black hole binaries resolved by pulsar timing arrays

Pulsar timing arrays (PTAs) are expected to detect gravitational waves (GWs) from individual low-redshift (z<1.5) compact supermassive (M>10^9 Msun) black hole (SMBH) binaries with orbital periods of approx. 0.1 - 10 yrs. Identifying the electromagnetic (EM) counterparts of these sources would provide confirmation of putative direct detections of GWs, present a rare opportunity to study the environments of compact SMBH binaries, and could enable the use of these sources as standard sirens for cosmology. Here we consider the feasibility of such an EM identification. We show that because the host galaxies of resolved PTA sources are expected to be exceptionally massive and rare, it should be possible to find unique hosts of resolved sources out to redshift z=0.2. At higher redshifts, the PTA error boxes are larger, and may contain as many as 100 massive-galaxy interlopers. The number of candidates, however, remains tractable for follow-up searches in upcoming wide-field EM surveys. We develop a toy model to characterize the dynamics and the thermal emission from a geometrically thin, gaseous disc accreting onto a PTA-source SMBH binary. Our model predicts that at optical and infrared frequencies, the source should appear similar to a typical luminous active galactic nucleus (AGN). However, owing to the evacuation of the accretion flow by the binary's tidal torques, the source might have an unusually low soft X-ray luminosity and weak UV and broad optical emission lines, as compared to an AGN powered by a single SMBH with the same total mass. For sources near z=1, the decrement in the rest-frame UV should be observable as an extremely red optical color. These properties would make the PTA sources stand out among optically luminous AGN, and could allow their unique identification.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Imprint of Accretion Disk-Induced Migration on Gravitational Waves from Extreme Mass Ratio Inspirals

We study the effects of a thin gaseous accretion disk on the inspiral of a stellar--mass black hole into a supermassive black hole. We construct a phenomenological angular momentum transport equation that reproduces known disk effects. Disk torques modify the gravitational wave phase evolution to detectable levels with LISA for reasonable disk parameters. The Fourier transform of disk-modified waveforms acquires a correction with a different frequency trend than post-Newtonian vacuum terms. Such inspirals could be used to detect accretion disks with LISA and to probe their physical parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Improved Models for Cosmic Infrared Background Anisotropies: New Constraints on the IR Galaxy Population

The power spectrum of cosmic infrared background (CIB) anisotropies is sensitive to the connection between star formation and dark matter halos over the entire cosmic star formation history. Here we develop a model that associates star-forming galaxies with dark matter halos and their subhalos. The model is based on a parameterized relation between the dust-processed infrared luminosity and (sub)halo mass. By adjusting 3 free parameters, we attempt to simultaneously fit the 4 frequency bands of the Planck measurement of the CIB anisotropy power spectrum. To fit the data, we find that the star-formation efficiency must peak on a halo mass scale of ~ 5x10^12 solar mass and the infrared luminosity per unit mass must increase rapidly with redshift. By comparing our predictions with a well-calibrated phenomenological model for shot noise, and with a direct observation of source counts, we show that the mean duty cycle of the underlying infrared sources must be near unity, indicating that the CIB is dominated by long-lived quiescent star formation, rather than intermittent short "star bursts". Despite the improved flexibility of our model, the best simultaneous fit to all four Planck channels remains relatively poor. We discuss possible further extensions to alleviate the remaining tension with the data. Our model presents a theoretical framework for a future joint analysis of both background anisotropy and source count measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Photodissociation of H2 in Protogalaxies: Modeling Self-Shielding in 3D Simulations

The ability of primordial gas to cool in proto-galactic haloes exposed to Lyman-Werner (LW) radiation is critically dependent on the self-shielding of H_2. We perform radiative transfer calculations of LW line photons, post-processing outputs from three-dimensional adaptive mesh refinement (AMR) simulations of haloes with T_vir > 10^4 K at redshifts around z=10. We calculate the optically thick photodissociation rate numerically, including the effects of density, temperature, and velocity gradients in the gas, as well as line overlap and shielding of H_2 by HI, over a large number of sight-lines. In low-density regions (n<10^4 cm^-3) the dissociation rates exceed those obtained using most previous approximations by more than an order of magnitude; the correction is smaller at higher densities. We trace the origin of the deviations primarily to inaccuracies of (i) the most common fitting formula (Draine & Bertoldi 1996) for the suppression of the dissociation rate and (ii) estimates for the effective shielding column density from local properties of the gas. The combined effects of gas temperature and velocity gradients are comparatively less important, typically altering the spherically averaged rate only by a factor of less than two. We present a simple modification to the DB96 fitting formula for the optically thick rate which improves agreement with our numerical results to within approx. 15 per cent, and can be adopted in future simulations. We find that estimates for the effective shielding column can be improved by using the local Sobolev length. Our correction to the H_2 self-shielding reduces the critical LW flux to suppress H_2-cooling in T_vir>10^4 K haloes by an order of magnitude; this increases the number of such haloes in which supermassive (approx. M=10^5 M_sun) black holes may have formed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Supermassive Black Hole Formation at High Redshifts Through a Primordial Magnetic Field

It has been proposed that primordial gas in early dark matter halos, with virial temperatures above 10^4 K, can avoid fragmentation and undergo rapid collapse, possibly resulting in a supermassive black hole (SMBH). This requires the gas to avoid cooling and to remain at temperatures near T=10^4 K. We show that this condition can be satisfied in the presence of a sufficiently strong primordial magnetic field, which heats the collapsing gas via ambipolar diffusion. If the field has a strength above B = 3.6 (comoving) nG, the collapsing gas is kept warm (T=10^4K) until it reaches the critical density n_crit=10^3 cm^{-3} at which the roto-vibrational states of H_2 approach local thermodynamic equilibrium. H_2-cooling then remains inefficient, and the gas temperature stays near 10^4K, even as it continues to collapse to higher densities. The critical magnetic field strength required to permanently suppress H_2-cooling is somewhat higher than upper limit of approx. 2 nG from the cosmic microwave background (CMB). However, it can be realized in the rare (2-3)-sigma regions of the spatially fluctuating B-field; these regions contain a sufficient number of halos to account for the z=6 quasar BHs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Suppression of HD-cooling in protogalactic gas clouds by Lyman-Werner radiation

It has been shown that HD molecules can form efficiently in metal-free gas collapsing into massive protogalactic halos at high redshift. The resulting radiative cooling by HD can lower the gas temperature to that of the cosmic microwave background, T_CMB=2.7(1+z)K, significantly below the temperature of a few 100 K achievable via H_2-cooling alone, and thus reduce the masses of the first generation of stars. Here we consider the suppression of HD-cooling by UV irradiation in the Lyman-Werner (LW) bands. We include photo-dissociation of both H_2 and HD, and explicitly compute the self-shielding and shielding of both molecules by neutral hydrogen as well as the shielding of HD by H_2. We use a simplified dynamical collapse model, and follow the chemical and thermal evolution of the gas, in the presence of a UV background. We find that a LW flux of J_crit = 1e-22 erg/cm^2/sr/s/Hz is able to suppress HD cooling and thus prevent collapsing primordial gas from reaching temperatures below 100 K. The main reason for the lack of HD cooling for J>J_crit is the partial photo-dissociation of H_2, which prevents the gas from reaching sufficiently low temperatures (T<150K) for HD to become the dominant coolant; direct HD photo-dissociation is unimportant except for a narrow range of fluxes and column densities. Since the prevention of HD-cooling requires only partial H_2 photo-dissociation, the critical flux J_crit is modest, and is below the UV background required to reionize the universe at redshift z=10-20. We conclude that HD-cooling can reduce the masses of typical stars only in rare halos forming well before the epoch of reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Origin and Detection of High-Redshift Supermassive Black Holes

Supermassive black holes (SMBHs) are common in local galactic nuclei, and SMBHs as massive as several billion solar masses already exist at redshift z=6. These earliest SMBHs may arise by the combination of Eddington-limited growth and mergers of stellar-mass seed BHs left behind by the first generation of metal-free stars, or by the rapid direct collapse of gas in rare special environments where the gas can avoid fragmenting into stars. In this contribution, I review these two competing scenarios. I also briefly mention some more exotic ideas and how the different models may be distinguished in the future by LISA and other instruments.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Witnessing the Birth of a Quasar

The coalescence of a supermassive black hole binary (SMBHB) is thought to be accompanied by an electromagnetic (EM) afterglow, produced by the viscous infall of the surrounding circumbinary gas disk after the merger. It has been proposed that once the merger has been detected in gravitational waves (GWs) by LISA, follow-up EM searches for this afterglow can help identify the EM counterpart of the LISA source. Here we study whether the afterglows may be sufficiently bright and numerous to be detectable in EM surveys alone. The viscous afterglow, which lasts for years to decades for SMBHBs in LISA's sensitivity window, is characterized by rapid increases in both the bolometric luminosity and in the spectral hardness of the source. If quasar activity is triggered by the same major galaxy mergers that produce SMBHBs, then the afterglow could be interpreted as a signature of the birth of a quasar. Using an idealized model for the post-merger viscous spreading of the circumbinary disk and the resulting light curve, and using the observed luminosity function of quasars as a proxy for the SMBHB merger rate, we delineate the survey requirements for identifying such birthing quasars. If circumbinary disks have a high disk surface density and viscosity, an all-sky soft X-ray survey with a sensitivity of ~<3x10^-14 erg s^-1 cm^-2 and a time resolution of ~months could identify dozens of birthing quasars with sustained brightening rates of >10%/yr. If >1% of the X-ray emission is reprocessed into optical frequencies, birthing quasars could also be identified in optical transient surveys such as the LSST. Distinguishing a birthing quasar from other variable sources may be facilitated by the monotonic hardening of its spectrum, but will likely remain challenging. This reinforces the notion that joint EM-plus-GW observations offer the best prospects for identifying the EM signatures of SMBHB mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Accretion Shocks in Clusters of Galaxies and their SZ Signature from Cosmological Simulations

Cold dark matter (CDM) hierarchical structure formation models predict the existence of large-scale accretion shocks between the virial and turnaround radii of clusters of galaxies. Kocsis et al. (2005) suggest that the Sunyaev-Zel'dovich (SZ) signal associated with such shocks might be observable with the next generation radio interferometer, ALMA. We study the three--dimensional distribution of accretion shocks around individual clusters of galaxies drawn from adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) simulations of LCDM (dark energy dominated CDM) models. In relaxed clusters, we find two distinct sets of shocks. One set ("virial shocks"), with Mach numbers of 2.5-4, is located at radii 0.9-1.3 Rvir, where Rvir is the spherical infall estimate of the virial radius, covering about 40-50% of the total surface area around clusters at these radii. Another set of stronger shocks ("external shocks") is located farther out, at about 3 Rvir, with large Mach numbers (~100), covering about 40-60% of the surface area. We simulate SZ surface brightness maps of relaxed massive galaxy clusters drawn from high resolution AMR runs, and conclude that ALMA should be capable of detecting the virial shocks in massive clusters of galaxies. More simulations are needed to improve estimates of astrophysical noise and to determine optimal observational strategies.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Hydrodynamical Response of a Circumbinary Gas Disk to Black Hole Recoil and Mass Loss

Finding electromagnetic (EM) counterparts of future gravitational wave (GW) sources would bring rich scientific benefits. A promising possibility, in the case of the coalescence of a super-massive black hole binary (SMBHB), is that prompt emission from merger-induced disturbances in a supersonic circumbinary disk may be detectable. We follow the post-merger evolution of a thin, zero-viscosity circumbinary gas disk with two-dimensional simulations, using the hydrodynamic code FLASH. We analyze perturbations arising from the 530 km/s recoil of a 10^6 M_sun binary, oriented in the plane of the disk, assuming either an adiabatic or a pseudo-isothermal equation of state for the gas. We find that a single-armed spiral shock wave forms and propagates outward, sweeping up about 20% of the mass of the disk. The morphology and evolution of the perturbations agrees well with those of caustics predicted to occur in a collisionless disk. Assuming that the disk radiates nearly instantaneously to maintain a constant temperature, we estimate the amount of dissipation and corresponding post-merger light-curve. The luminosity rises steadily on the time-scale of months, and reaches few times 10^{43} erg/s, corresponding to about 10% of the Eddington luminosity of the central SMBHB. We also analyze the case in which gravitational wave emission results in a 5% mass loss in the merger remnant. The mass-loss reduces the shock overdensities and the overall luminosity of the disk by 15-20%, without any other major effects on the spiral shock pattern.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Probing Cosmology with Weak Lensing Peak Counts

We propose counting peaks in weak lensing (WL) maps, as a function of their height, to probe models of dark energy and to constrain cosmological parameters. Because peaks can be identified in two-dimensional WL maps directly, they can provide constraints that are free from potential selection effects and biases involved in identifying and determining the masses of galaxy clusters. We have run cosmological N-body simulations to produce WL convergence maps in three models with different constant values of the dark energy equation of state parameter, w=-0.8, -1, and -1.2, with a fixed normalization of the primordial power spectrum (corresponding to present-day normalizations of sigma8=0.742, 0.798, and 0.839, respectively). By comparing the number of WL peaks in 8 convergence bins in the range of -0.1 < kappa < 0.2, in multiple realizations of a single simulated 3x3 degree field, we show that the first (last) pair of models can be distinguished at the 95% (85%) confidence level. A survey with depth and area (20,000 sq. degrees), comparable to those expected from LSST, should have a factor of approx. 50 better parameter sensitivity. We find that relatively low-amplitude peaks (kappa = 0.03), which typically do not correspond to a single collapsed halo along the line of sight, account for most of this sensitivity. We study a range of smoothing scales and source galaxy redshifts (z_s). With a fixed source galaxy density of 15/arcmin^2, the best results are provided by the smallest scale we can reliably simulate, 1 arcminute, and z_s=2 provides substantially better sensitivity than z_s< 1.5.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Assembly of Supermassive Black Holes at High Redshifts

The supermassive black holes (SMBHs) massive enough to power the bright redshift ~6 quasars observed in the Sloan Digital Sky Survey (SDSS) are thought to have assembled by mergers and/or gas accretion from less massive "seed" BHs. If the seeds are the ~100 Msol remnant BHs from the first generation of stars, they must be in place well before redshift z=6, and must avoid being ejected from their parent proto-galaxies by the large (few hundred km/s) kicks they suffer from gravitational-radiation induced recoil during mergers with other BHs. We simulate the SMBH mass function at redshift z=6 using dark matter (DM) halo merger trees, coupled with a prescription for the halo occupation fraction, accretion histories, and radial recoil trajectories of the growing BHs. Our purpose is (i) to map out plausible scenarios for successful assembly of the z~6 quasar BHs by exploring a wide region of parameter space, and (ii) to predict the rate of low-frequency gravitational wave events detectable by the Laser Interferometer Space Antenna (LISA) for each such scenario. Our main findings are as follows: (1) ~100 Msol seed BHs can grow into the SDSS quasar BHs without super--Eddington accretion, but only if they form in minihalos by z~30 and subsequently accrete ~60% of the time; (2) the scenarios with optimistic assumptions required to explain the SDSS quasar BHs overproduce the mass density in lower--mass (10^5 to 10^7 Msol) SMBHs by a factor of 100-1000, unless seeds stop forming, or accrete at a severely diminished rates or duty cycles (e.g. due to feedback), at z < 20-30. We also present several successful assembly models and their LISA detection rates, including a "maximal" model that gives the highest rate (~30/yr at z~6) without overproducing the total SMBH density.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Constraining Cosmology with High Convergence Regions in Weak Lensing Surveys

We propose to use a simple observable, the fractional area of "hot spots" in weak gravitational lensing mass maps which are detected with high significance, to determine background cosmological parameters. Because these high-convergence regions are directly related to the physical nonlinear structures of the universe, they derive cosmological information mainly from the nonlinear regime of density fluctuations. We show that in combination with future cosmic microwave background anisotropy measurements, this method can place constraints on cosmological parameters that are comparable to those from the redshift distribution of galaxy cluster abundances. The main advantage of the statistic proposed in this paper is that projection effects, normally the main source of uncertainty when determining the presence and the mass of a galaxy cluster, here serve as a source of information.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Correlations Between Central Massive Objects And Their Host Galaxies: From Bulgeless Spirals to Ellipticals

Recent observations by Ferrarese et al. (2006) and Wehner et al. (2006) reveal that a majority of galaxies contain a central massive object (CMO), either a supermassive black hole (SMBH) or a compact stellar nucleus, regardless of the galaxy mass or morphological type, and that there is a tight relation between the masses of CMOs and those of the host galaxies. Several recent studies show that feedback from black holes can successfully explain the $\msigma$ correlation in massive elliptical galaxies that contain SMBHs. However, puzzles remain in spirals or dwarf spheroids that do not appear to have black holes but instead harbor a compact central stellar cluster. Here we use three-dimensional, smoothed particle hydrodynamics simulations of isolated galaxies to study the formation and evolution of CMOs in bulgeless disk galaxies, and simulations of merging galaxies to study the transition of the CMO--host mass relation from late-type bulgeless spirals to early-type ellipticals. Our results suggest that the observed correlations may be established primarily by the depletion of gas in the central region by accretion and star-formation, and may hold for all galaxy types. A systematic search for CMOs in the nuclei of bulgeless disk galaxies would offer a test of this conclusion. (Abridged)

0 citations0 reviewsNo code linkedOpen access