Source author record

Priyamvada Natarajan

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

60works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

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

60 published item(s)

preprint2026arXiv

Evolutionary Tracks and Spectral Properties of Quasi-stars and Their Correlation with Little Red Dots

JWST has revealed a population of red, compact, high-redshift (${z\sim3-10}$) objects referred to as ``Little Red Dots'' (LRDs). These objects exhibit unusual spectral features reminiscent of stellar spectra with blackbody-like SEDs, large hydrogen Balmer breaks, Balmer line absorption, and classical stellar absorption features such as calcium H&K and the calcium triplet. Following the recent suggestion that these may be actively accreting direct-collapse black holes in the process of assembly, i.e. quasi-stars, we present evolutionary models of quasi-stars using our recently released, publicly available MESA-QUEST modeling framework. We compute a grid of models spanning a range of black hole masses and predict the luminosities, temperatures, surface gravities, and lifetimes of these objects. We find that these models lie along a Hayashi track once they hit their ``late-stage'' which constitutes the majority of their lives ($\sim 20$~Myr). We present scaling relations for estimating the mass of a quasi-star as a function of the bolometric luminosity, as well as the bolometric luminosity as a function of the effective temperature for the Hayashi track. The short lifetimes in tandem with the observed number density of LRDs imply the possibility that every supermassive black hole was once a quasi-star. We compare synthetic spectra of our quasi-star models to observations of LRDs, and show that these models are broadly capable of reproducing the continuum spectra of observed LRDs. These results indicate that quasi-stars are promising candidates for the origin of supermassive black holes via direct collapse in the early universe.

preprint2026arXiv

MESA-QUEST: Tracing the formation of direct collapse black hole seeds via quasi-stars

The origin of the first supermassive black holes (SMBHs) observed at redshifts $z\geq 9$ remains one of the most challenging open questions in astrophysics. Their rapid emergence suggests that massive ``heavy seeds'' must have formed early, possibly through the direct collapse of pristine gas clouds in the first galaxies. We present MESA-QUEST, a new framework built upon the Modules for Experiments in Stellar Astrophysics (MESA) code, designed to model the structure and evolution of quasi-stars -- massive, radiation-supported envelopes hosting accreting black holes at their cores -- believed to be the progenitors of direct-collapse black hole (DCBH) seeds. Our implementation introduces flexible boundary conditions representing both Bondi accretion and saturated-convection regimes, and explores the impact of several stellar wind and mass-loss prescriptions, including Reimers, Dutch, and super-Eddington radiation-driven winds. We find that quasi-stars can grow central black holes to $\geq 10^3\,M_{\odot}$ under favorable conditions, with saturated-convection models yielding BH-to-total mass ratios up to 0.55$M_*$ -- five times higher than Bondi-limited cases. However, strong radiation-driven winds can dramatically curtail growth, potentially quenching heavy-seed formation unless balanced by sustained envelope accretion. Our results delineate the physical limits under which quasi-stars can remain stable and produce heavy seeds capable of evolving into the earliest SMBHs detected by JWST and Chandra. Future extensions will incorporate rotation, magnetic fields, and GR-radiation hydrodynamics to refine accretion physics and constrain the viability of the quasi-star pathway for early SMBH formation.

preprint2025arXiv

Constraining the Nature of Dark Matter from Tidal Radii of Cluster Galaxy Subhalos

Gravitational lensing by galaxy clusters provides a powerful probe of the spatial distribution of dark matter and its microphysical properties. Strong and weak lensing constraints on the density profiles of subhalos and their truncation radii offer key diagnostics for distinguishing between collisionless cold dark matter (CDM) and self-interacting dark matter (SIDM). Notably, in the strongly collisional SIDM regime, subhalo core collapse and enhanced mass loss from ram-pressure stripping predict steeper central density slopes and more compact truncation radii--features that are directly testable with current lensing data. We analyze subhalo truncation in eight lensing clusters (Abell 2218, 383, 963, 209, 2390, and MACS J0416.1, J1206.2, J1149.6) that span the redshift range <$z_\text{spec}$>$ \simeq 0.17$-$0.54$ with virial masses $M_{200} \simeq0.41$-$2.2\times 10^{15}$ M$_\odot$ to constrain SIDM versus CDM. Our results indicate that the outer spatial extents of subhalos are statistically consistent with CDM, corroborated by redshift- and mass-matched analogs from the Illustris-TNG simulations. We conclude that the tidal radii of cluster galaxy subhalos serve as an important and complementary diagnostic of the nature of dark matter in these violent, dense environments.

preprint2022arXiv

Astrophysical and Cosmological Probes of Dark Matter

While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for these strides is being laid now, yet much remains to be done. Chief among the upcoming challenges is establishing the theoretical foundation needed to harness the full potential of new observables in the astrophysical and cosmological domains, spanning the early Universe to the inner portions of galaxies and the stars therein. Identifying the nature of dark matter will also entail repurposing and implementing a wide range of theoretical techniques from outside the typical toolkit of astrophysics, ranging from effective field theory to the dramatically evolving world of machine learning and artificial-intelligence-based statistical inference. Through this work, the theory frontier will be at the heart of dark matter discoveries in the upcoming decade.

preprint2022arXiv

Galaxies in the central regions of simulated galaxy clusters

In this paper, we assess the impact of numerical resolution and of the implementation of energy input from AGN feedback models on the inner structure of cluster sub-haloes in hydrodynamic simulations. We compare several zoom-in re-simulations of a sub-sample of the cluster-sized haloes studied in Meneghetti et al. (2020), obtained by varying mass resolution, softening length and AGN energy feedback scheme. We study the impact of these different setups on the subhalo abundances, their radial distribution, their density and mass profiles and the relation between the maximum circular velocity, which is a proxy for subhalo compactness. Regardless of the adopted numerical resolution and feedback model, subhaloes with masses Msub < 1e11Msun/h, the most relevant mass-range for galaxy-galaxy strong lensing, have maximum circular velocities ~30% smaller than those measured from strong lensing observations of Bergamini et al. (2019). We also find that simulations with less effective AGN energy feedback produce massive subhaloes (Msub> 1e11 Msun/h ) with higher maximum circular velocity and that their Vmax - Msub relation approaches the observed one. However the stellar-mass number count of these objects exceeds the one found in observations and we find that the compactness of these simulated subhaloes is the result of an extremely over-efficient star formation in their cores, also leading to larger-than-observed subhalo stellar mass. We conclude that simulations are unable to simultaneously reproduce the observed stellar masses and compactness (or maximum circular velocities) of cluster galaxies. Thus, the discrepancy between theory and observations that emerged from the analysis of Meneghetti et al. (2020) persists. It remains an open question as to whether such a discrepancy reflects limitations of the current implementation of galaxy formation models or the LCDM paradigm.

preprint2022arXiv

Toward determining the number of observable supermassive black hole shadows

We present estimates for the number of shadow-resolved supermassive black hole (SMBH) systems that can be detected using radio interferometers, as a function of angular resolution, flux density sensitivity, and observing frequency. Accounting for the distribution of SMBHs across mass, redshift, and accretion rate, we use a new semi-analytic spectral energy distribution model to derive the number of SMBHs with detectable and optically thin horizon-scale emission. We demonstrate that (sub)millimeter interferometric observations with ${\sim}0.1$ $μ$as resolution and ${\sim}1$ $μ$Jy sensitivity could access ${>}10^6$ SMBH shadows. We then further decompose the shadow source counts into the number of black holes for which we could expect to observe the first- and second-order lensed photon rings. Accessing the bulk population of first-order photon rings requires ${\lesssim}2$ $μ$as resolution and ${\lesssim}0.5$ mJy sensitivity, while doing the same for second-order photon rings requires ${\lesssim}0.1$ $μ$as resolution and ${\lesssim}5$ $μ$Jy sensitivity. Our model predicts that with modest improvements to sensitivity, as many as $\sim$5 additional horizon-resolved sources should become accessible to the current Event Horizon Telescope (EHT), while a next-generation EHT observing at 345 GHz should have access to ${\sim}$3 times as many sources. More generally, our results can help guide enhancements of current arrays and specifications for future interferometric experiments that aim to spatially resolve a large population of SMBH shadows or higher-order photon rings.

preprint2021arXiv

Exploring the high-redshift PBH-$Λ$CDM Universe: early black hole seeding, the first stars and cosmic radiation backgrounds

We explore the observational implications of a model in which primordial black holes (PBHs) with a broad birth mass function ranging in mass from a fraction of a solar mass to $\sim$10$^6$ M$_{\odot}$, consistent with current observational limits, constitute the dark matter component in the Universe. The formation and evolution of dark matter and baryonic matter in this PBH-\LambdaCDM~Universe are presented. In this picture, PBH DM mini-halos collapse earlier than in standard \LambdaCDM, baryons cool to form stars at $z\sim15-20$, and growing PBHs at these early epochs start to accrete through Bondi capture. The volume emissivity of these sources peaks at $z\sim20$ and rapidly fades at lower redshifts. As a consequence, PBH DM could also provide a channel to make early black hole seeds and naturally account for the origin of an underlying dark matter halo - host galaxy and central black hole connection that manifests as the $M_{\rm bh}-σ$ correlation. To estimate the luminosity function and contribution to integrated emission power spectrum from these high-redshift PBH DM halos, we develop a Halo Occupation Distribution (HOD) model. In addition to tracing the star formation and reionizaton history, it permits us to evaluate the Cosmic Infrared and X-ray Backgrounds (CIB and CXB). We find that accretion onto PBHs/AGN successfully accounts for the detected backgrounds and their cross-correlation, with the inclusion of an additional IR stellar emission component. Detection of the deep IR source count distribution by the JWST could reveal the existence of this population of high-redshift star-forming and accreting PBH DM.

preprint2020arXiv

A Link Between Ram Pressure Stripping and Active Galactic Nuclei

The dense environment of a galaxy cluster can radically transform the content of in-falling galaxies. Recent observations have found a significant population of active galactic nuclei (AGN) within "jellyfish galaxies," galaxies with trailing tails of gas and stars that indicate significant ram pressure stripping. The relationship between AGN and ram pressure stripping is not well understood. In this letter, we investigate the connection between AGN activity and ram pressure in a fully cosmological setting for the first time using the RomulusC simulation, one of the highest resolution simulations of a galaxy cluster to date. We find unambiguous morphological evidence for ram pressure stripping. For lower mass galaxies (with stellar masses < 10^9.5 solar masses) both star formation and black hole accretion are suppressed by ram pressure before they reach pericenter, whereas for more massive galaxies accretion onto the black hole is enhanced during pericentric passage. Our analysis also indicates that as long as the galaxy retains gas, AGN with higher Eddington ratios are more likely to be the found in galaxies experiencing higher ram pressure. We conclude that prior to quenching star formation, ram pressure triggers enhanced accretion onto the black hole, which then produces heating and outflows due to AGN feedback. AGN feedback may in turn serve to aid in the quenching of star formation in tandem with ram pressure.

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

hybrid-Lenstool: A self-consistent algorithm to model galaxy clusters with strong- and weak-lensing simultaneously

We present a new galaxy cluster lens modeling approach, hybrid-Lenstool, that is implemented in the publicly available modeling software Lenstool. hybrid-Lenstool combines a parametric approach to model the core of the cluster, and a non-parametric (free-form) approach to model the outskirts. hybrid-Lenstool optimizes both strong- and weak-lensing constraints simultaneously (Joint-Fit), providing a self-consistent reconstruction of the cluster mass distribution on all scales. In order to demonstrate the capabilities of the new algorithm, we tested it on a simulated cluster. hybrid-Lenstool yields more accurate reconstructed mass distributions than the former Sequential-Fit approach where the parametric and the non-parametric models are optimized successively. Indeed, we show with the simulated cluster that the mass density profile reconstructed with a Sequential-Fit deviates form the input by $2-3σ$ at all scales while the Joint-Fit gives a profile that is within $1-1.5σ$ of the true value. This gain in accuracy is consequential for recovering mass distributions exploiting cluster lensing and therefore for all applications of clusters as cosmological probes. Finally we found that the Joint-Fit approach yields shallower slope of the inner density profile than the Sequential-Fit approach, thus revealing possible biases in previous lensing studies.

preprint2020arXiv

The BUFFALO HST Survey

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

preprint2020arXiv

The distribution of dark matter and gas spanning six megaparsecs around the post-merger galaxy cluster MS0451-03

Using the largest mosaic of Hubble Space Telescope images around a galaxy cluster, we map the distribution of dark matter throughout a $\sim$$6\times6$ Mpc$^2$ area centred on the cluster MS 0451-03 ($z=0.54$, $M_{200}=1.65\times10^{15} \rm{M}_\odot$). Our joint strong- and weak-lensing analysis shows three possible filaments extending from the cluster, encompassing six group-scale substructures. The dark-matter distribution in the cluster core is elongated, consists of two distinct components, and is characterized by a concentration parameter of $c_{200}=3.79\pm0.36$. By contrast, XMM-Newton observations show the gas distribution to be more spherical, with excess entropy near the core, and a lower concentration of $c_{200}=2.35^{+0.89}_{-0.70}$ (assuming hydrostatic equilibrium). Such a configuration is predicted in simulations of major mergers 2-7Gyr after the first core passage, when the two dark-matter halos approach second turnaround, and before their gas has relaxed. This post-merger scenario finds further support in optical spectroscopy of the cluster's member galaxies, which shows that star formation was abruptly quenched 5 Gyr ago. MS 0451-03 will be an ideal target for future studies of the growth of structure along filaments, star-formation processes after a major merger, and the late-stage evolution of cluster collisions.

preprint2019arXiv

Gravitational wave probes of dark matter: challenges and opportunities

In this white paper, we discuss the prospects for characterizing and identifying dark matter using gravitational waves, covering a wide range of dark matter candidate types and signals. We argue that present and upcoming gravitational wave probes offer unprecedented opportunities for unraveling the nature of dark matter and we identify the most urgent challenges and open problems with the aim of encouraging a strong community effort at the interface between these two exciting fields of research.

preprint2016arXiv

A new method to break the mass sheet degeneracy using aperture moments

Mass determinations from gravitational lensing shear and the higher order estimator flexion are both subject to the mass sheet degeneracy. Mass sheet degeneracy refers to a transformation that leaves the reduced shear and flexion invariant. In general, this transformation can be approximated by the addition of a constant surface mass density sheet. We propose a new technique to break the mass sheet degeneracy. The method uses mass moments of the shear or flexion fields in combination with convergence information derived from number counts which exploit the magnification bias. The difference between the measured mass moments provides an estimator for the magnitude of the additive constant that is the mass-sheet. For demonstrating this, we derive relations that hold true in general for n-th order moments and show how they can be employed effectively to break the degeneracy. We investigate the detectability of this degeneracy parameter from our method and find that the degeneracy parameter can be feasibly determined from stacked galaxy-galaxy lensing data and cluster lensing data. Furthermore, we compare the signal-to-noise ratios of convergence information from number counts with shear and flexion for SIS and NFW models. We find that the combination of shear and flexion performs best on galaxy and cluster scales and the convergence information can therefore be used to break the mass sheet degeneracy without quality loss in the mass reconstruction. In summary, there is power in the combination of shear, flexion, convergence and their higher order moments. With the anticipated wealth of lensing data from upcoming and future satellite missions - EUCLID and WFIRST - this technique will be feasible.

preprint2016arXiv

Detecting Direct Collapse Black Holes: making the case for CR7

We propose that one of the sources in the recently detected system CR7 by Sobral et al. (2015) through spectro-photometric measurements at $z = 6.6$ harbors a direct collapse blackhole (DCBH). We argue that the LW radiation field required for direct collapse in source A is provided by sources B and C. By tracing the LW production history and star formation rate over cosmic time for the halo hosting CR7 in a $Λ$CDM universe, we demonstrate that a DCBH could have formed at $z\sim 20$. The spectrum of source A is well fit by nebular emission from primordial gas around a BH with MBH $\sim 4.4 \times 10^6 \ \rm M_{\odot}$ accreting at a 40 % of the Eddington rate, which strongly supports our interpretation of the data. Combining these lines of evidence, we argue that CR7 might well be the first DCBH candidate.

preprint2016arXiv

New constraints on direct collapse black hole formation in the early Universe

Direct collapse black holes (DCBH) have been proposed as a solution to the challenge of assembling supermassive black holes by $z>6$ to explain the bright quasars observed at this epoch. The formation of a DCBH seed with $\rm M_{BH}\sim10^{4-5}\ \rm M_{\odot}$ requires a pristine atomic-cooling halo to be illuminated by an external radiation field that is sufficiently strong to entirely suppress H$_{2}$ cooling in the halo. Many previous studies have attempted to constrain the critical specific intensity that is likely required to suppress H$_{2}$ cooling, denoted as $J_{\rm crit}$. However, these studies have typically assumed that the incident external radiation field can be modeled with a black-body spectrum. Under this assumption, it is possible to derive a {unique} value for $J_{\rm crit}$ that depends only on the temperature of the black-body. In this study we consider a more realistic spectral energy distribution (SED) for the external source of radiation that depends entirely on its star formation history and age. The rate of destruction of the species responsible for suppressing molecular hydrogen cooling depends on the detailed shape of the SED. Therefore the value of $J_{\rm crit}$ is tied to the shape of the incident SED of the neighbouring galaxy. We fit a parametric form to the rates of destruction of H$_2$ and H$^-$ that permit direct collapse. Owing to this, we find that $J_{\rm crit}$ is not a fixed threshold but can lie anywhere in the range $J_{\rm crit} \sim 0.5$--$10^{3}$, depending on the details of the source stellar population.

preprint2015arXiv

Are Ultra-faint Galaxies at z=6-8 Responsible for Cosmic Reionization ? Combined Constraints from the Hubble Frontier Fields Clusters and Parallels

We use deep Hubble Space Telescope imaging of the Frontier Fields to accurately measure the galaxy rest-frame ultraviolet luminosity function (UV LF) in the redshift range $z \sim 6-8$. We combine observations in three lensing clusters A2744, MACS0416, MACS0717 and their associated parallels fields to select high-redshift dropout candidates. We use the latest lensing models to estimate the flux magnification and the effective survey volume in combination with completeness simulations performed in the source plane. We report the detection of 227 galaxy candidates at $z=6-7$ and 25 candidates at $z \sim 8$. While the total survey area is about 4 arcmin$^{2}$ in each parallel field, it drops to about 0.6 to 1 arcmin$^{2}$ in the cluster core fields because of the strong lensing. We compute the UV luminosity function at $z \sim 7$ using the combined galaxy sample and perform Monte Carlo simulations to determine the best fit Schechter parameters. We are able to reliably constrain the LF down to an absolute magnitude of $M_{UV}=-15.25$, which corresponds to 0.005$L^{\star}$. More importantly, we find that the faint-end slope remains steep down to this magnitude limit with $α=-2.04_{-0.17}^{+0.13}$. Our results confirm the most recent results in deep blank fields but extend the LF measurements more than two magnitudes deeper. The UV LF at $z \sim 8$ is not very well constrained below $M_{UV}=-18$ due to the small number statistics and incompleteness uncertainties. To assess the contribution of galaxies to cosmic reionization we derive the UV luminosity density at $z\sim7$ by integrating the UV LF down to an observationally constrained limit of $M_{UV} = -15$. We show that our determination of Log($ρ_{UV}$)=$26.2\pm0.13$ (erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$) can be sufficient to maintain the IGM ionized.

preprint2015arXiv

Bulge-driven Fueling of Seed Black Holes

We examine radiation-regulated accretion onto intermediate-mass and massive black holes (BHs) embedded in a bulge component. Using spherically symmetric one-dimensional radiation-hydrodynamics simulations, we track the growth of BHs accreting from a cold, neutral gas reservoir with temperature T=10^4 K. We find that the accretion rate of BHs embedded in bulges is proportional to r_{B,eff}/r_B, where r_{B,eff} is the increased effective Bondi radius that includes the gravitational potential of the bulge, and r_B is the Bondi radius of the BH. The radiative feedback from the BH suppresses the cold accretion rate to ~1 percent of the Bondi rate when a bulge is not considered. However, we find that the BH fueling rate increases rapidly when the bulge mass M_bulge is greater than the critical value of 10^6 M_sun and is proportional to M_bulge. Since the critical bulge mass is independent of the central BH mass M_{BH}, the growth rate of BHs with masses of 10^2, 10^4, and 10^6 M_sun exhibits distinct dependencies on the bulge-to-BH mass ratio. Our results imply that light seed BHs (<= 10^2 M_sun) which might be the remnants of the Pop III stars, cannot grow through accretion coevally with the early assembly of the bulge of the host galaxies until the bulge reaches the critical mass. However, massive BH seeds (>= 10^5 M_sun) that may form via direct collapse, are more likely to be embedded in a supercritical bulge and thus can grow efficiently coupling to the host galaxies and driving the early evolution of the M_{BH}-$σ$ relationship.

preprint2015arXiv

Hubble Frontier Fields : A High-Precision Strong-Lensing Mass Model of the Massive Galaxy Cluster Abell 2744 using 180 Multiple Images

We present a high-precision mass model of galaxy cluster Abell 2744, based on a strong-gravitational-lensing analysis of the \emph{Hubble Space Telescope Frontier Fields} (HFF) imaging data, which now include both \emph{Advanced Camera for Surveys} and \emph{Wide-Field Camera 3} observations to the final depth. Taking advantage of the unprecedented depth of the visible and near-infrared data, we identify 34 new multiply imaged galaxies, bringing the total to 61, comprising 181 individual lensed images. In the process, we correct previous erroneous identifications and positions of multiple systems in the northern part of the cluster core. With the \textsc{Lenstool} software and the new sets of multiple images, we model the cluster using two cluster-scale dark matter halos plus galaxy-scale halos for the cluster members. Our best-fit model predicts image positions with an \emph{RMS} error of 0.69$\arcsec$, which constitutes an improvement by almost a factor of two over previous parametric models of this cluster. We measure the total projected mass inside a 200~kpc aperture as ($2.162\pm 0.005$)$\times 10^{14}M_{\odot}$, thus reaching 1\% level precision for the second time, following the recent HFF measurement of MACSJ0416.1-2403. Importantly, the higher quality of the mass model translates into an overall improvement by a factor of 4 of the derived magnification factor. % for the high-redshift lensed background galaxies. Together with our previous HFF gravitational lensing analysis, this work demonstrates that the HFF data enables high-precision mass measurements for massive galaxy clusters and the derivation of robust magnification maps to probe the early Universe.

preprint2015arXiv

On R-W1 as A Diagnostic to Discover Obscured Active Galactic Nuclei in Wide-Area X-ray Surveys

Capitalizing on the all-sky coverage of {\it WISE}, and the 35\% and 50\% sky coverage from SDSS and Pan-STARRS, respectively, we explore the efficacy of $m_{R}$ (optical) - $m_{3.4 μm}$ (mid-infrared), hereafter $R-W1$, as a color diagnostic to identify obscured supermassive black hole accretion in wide-area X-ray surveys. We use the $\sim$16.5 deg$^2$ Stripe 82 X-ray survey data as a test-bed to compare $R-W1$ with $R-K$, an oft-used obscured AGN selection criterion, and examine where different classes of objects lie in this parameter space. Most stars follow a well-defined path in $R-K$ vs. $R-W1$ space. We demonstrate that optically normal galaxies hosting X-ray AGN at redshifts $0.5<z<1$ can be recovered with an $R-W1>4$ color-cut, while they typically are not selected as AGN based on their $W1-W2$ colors. Additionally, different observed X-ray luminosity bins favor different regions in $R-W1$ parameter space: moderate luminosity AGN ($10^{43}$ erg s$^{-1} < L_{\rm 0.5-10 keV} < 10^{44}$ erg s$^{-1}$) tend to have red colors while the highest luminosity AGN ($L_{\rm 0.5-10 keV} > 10^{45}$ erg s$^{-1}$) have bluer colors; higher spectroscopic completeness of the Stripe 82X sample is needed to determine whether this is a selection effect or an intrinsic property. Finally, we parameterize X-ray obscuration of Stripe 82X AGN by calculating their hardness ratios (HRs) and find no clear trends between HR and optical reddening. Our results will help inform best-effort practices in following-up obscured AGN candidates in current and future wide-area, shallow X-ray surveys, including the all-sky {\it eROSITA} mission.

preprint2015arXiv

Tidal Disruption Events by a Massive Black Hole Binary

Massive black hole binaries (MBHBs) are a natural byproduct of galaxy mergers. Previous studies have shown that flares from stellar tidal disruption events (TDEs) are modified by the presence of a secondary perturber, causing interruptions in the light curve. We study the dynamics of TDE debris in the presence of a milliparsec-separated MBHB by integrating ballistic particle orbits in the time-varying potential of the binary. We find that gaps in the light curve appear when material misses the accretion radius on its first return to pericentre. Subsequent recurrences can be decomposed into "continuous" and "delayed" components, which exhibit different behaviour. We find that this potential can substantially alter the locations of stream self-intersections. When debris is confined to the plane, we find that close encounters with the secondary BH leave noticeable signatures on the fallback rate and can result in significant accretion onto the secondary BH. Tight, equal-mass MBHBs accrete equally, periodically trading the infalling stream.

preprint2014arXiv

Hubble Frontier Fields : A High Precision Strong Lensing Analysis of Galaxy Cluster MACSJ0416.1-2403 using ~200 Multiple Images

We present a high-precision mass model of the galaxy cluster MACSJ0416.1-2403, based on a strong-gravitational-lensing analysis of the recently acquired Hubble Space Telescope Frontier Fields (HFF) imaging data. Taking advantage of the unprecedented depth provided by HST/ACS observations in three passbands, we identify 51 new multiply imaged galaxies, quadrupling the previous census and bringing the grand total to 68, comprising 194 individual lensed images. Having selected a subset of the 57 most securely identified multiply imaged galaxies, we use the Lenstool software package to constrain a lens model comprised of two cluster-scale dark-matter halos and 98 galaxy-scale halos. Our best-fit model predicts image positions with an $RMS$ error of 0.68'', which constitutes an improvement of almost a factor of two over previous, pre-HFF models of this cluster. We find the total projected mass inside a 200~kpc aperture to be $(1.60\pm0.01)\times 10^{14}\ M_\odot$, a measurement that offers a three-fold improvement in precision, reaching the percent level for the first time in any cluster. Finally, we quantify the increase in precision of the derived gravitational magnification of high-redshift galaxies and find an improvement by a factor of $\sim$2.5 in the statistical uncertainty. Our findings impressively confirm that HFF imaging has indeed opened the domain of high-precision mass measurements for massive clusters of galaxies.

preprint2014arXiv

Hubble Frontier Fields: The Geometry and Dynamics of the Massive Galaxy Cluster Merger MACSJ0416.1-2403

We use a joint optical/X-ray analysis to constrain the geometry and history of the ongoing merging event in the massive galaxy cluster MACSJ0416.1-2403 (z=0.397). Our investigation of cluster substructure rests primarily on a combined strong- and weak-lensing mass reconstruction based on the deep, high-resolution images obtained for the Hubble Frontier Fields initiative. To reveal the system's dynamics, we complement this lensing analysis with a study of the intra-cluster gas using shallow Chandra data, and a three-dimensional model of the distribution and motions of cluster galaxies derived from over 100 spectroscopic redshifts. The multi-scale grid model obtained from our combined lensing analysis extends the high-precision strong-lensing mass reconstruction recently performed to cluster-centric distances of almost 1 Mpc. Our analysis detects the two well known mass concentrations in the cluster core. A pronounced offset between collisional and collisionless matter is only observed for the SW cluster component, while excellent alignment is found for the NE cluster. Both the lensing analysis and the distribution of cluster light strongly suggest the presence of a third massive structure, almost 2 arcmin SW of the cluster centre. Since no X-ray emission is detected in this region, we conclude that this structure is non-virialised and speculate that it might be part of a large-scale filament almost aligned with our line of sight. Combining all evidence from the distribution of dark and luminous matter, we propose two alternative scenarios for the trajectories of the components of MACSJ0416.1-2403. Upcoming deep X-ray observations that allow the detection of shock fronts, cold cores, and sloshing gas (all key diagnostics for studies of cluster collisions) will allow us to test, and distinguish between these two scenarios.

preprint2014arXiv

New Constraints on the Faint-end of the UV Luminosity Function at z~7-8 using the Gravitational Lensing of the Hubble Frontier Fields Cluster A2744

Exploiting the power of gravitational lensing, the Hubble Frontier Fields (HFF) program aims at observing six massive galaxy clusters to explore the distant Universe far beyond the depth limits of blank field surveys. Using the complete Hubble Space Telescope observations of the first HFF cluster Abell 2744, we report the detection of 50 galaxy candidates at $z \sim 7$ and eight candidates at $z \sim 8$ in a total survey area of 0.96 arcmin$^{2}$ in the source plane. Three of these galaxies are multiply-imaged by the lensing cluster. Using an updated model of the mass distribution in the cluster we were able to calculate the magnification factor and the effective survey volume for each galaxy in order to compute the ultraviolet galaxy luminosity function at both redshifts 7 and 8. Our new measurements extend the $z \sim 7$ UV LF down to an absolute magnitude of $M_{UV} \sim -15.5$. We find a characteristic magnitude of $M^{\star}_{UV}=-20.63^{+0.69}_{-0.56}$ mag and a faint-end slope $α= -1.88^{+0.17}_{-0.20}$ close to previous determinations in blank fields. We show here for the first time that this slope remains steep down to very faint luminosities of 0.01$L^{\star}$. Although prone to large uncertainties, our results at $z \sim 8$ also seem to confirm a steep faint-end slope below 0.1$L^{\star}$. The HFF program is therefore providing an extremely efficient way to study the faintest galaxy populations at $z > 7$ that would otherwise be inaccessible with current instrumentation. The full sample of six galaxy clusters will provide yet better constraints on the build-up of galaxies at early epochs and their contribution to cosmic reionization.

preprint2014arXiv

Probing the z>6 Universe with the first Hubble Frontier Fields cluster Abell 2744

The Hubble Frontier Fields (HFF) program combines the capabilities of the Hubble Space Telescope (HST) with the gravitational lensing of massive galaxy clusters to probe the distant Universe to an unprecedented depth. Here, we present the results of the first combined HST and Spitzer observations of the cluster Abell 2744. We combine the full near-infrared data with ancillary optical images to search for gravitationally lensed high-redshift (z > 6) galaxies. We report the detection of 15 I814-dropout candidates at z ~ 6-7 and one Y105-dropout at z ~ 8 in a total survey area of 1.43 arcmin^2 in the source plane. The predictions of our lens model allow us to also identify five multiply-imaged systems lying at redshifts between z ~ 6 and z ~ 8. Thanks to constraints from the mass distribution in the cluster, we were able to estimate the effective survey volume corrected for completeness and magnification effects. This was in turn used to estimate the rest-frame ultraviolet luminosity function (LF) at z ~ 6-8. Our LF results are generally in agreement with the most recent blank field estimates, confirming the feasibility of surveys through lensing clusters. Although based on a shallower observations than what will be achieved in the final dataset including the full ACS observations, the LF presented here extends down to Muv ~ -18.5 at z ~ 7 with one identified object at Muv ~ -15 thanks to the highly-magnified survey areas. This early study forecasts the power of using massive galaxy clusters as cosmic telescopes and its complementarity to blank fields.

preprint2014arXiv

Rapid growth of seed black holes in the early universe by supra-exponential accretion

Mass accretion by black holes (BHs) is typically capped at the Eddington rate, when radiation's push balances gravity's pull. However, even exponential growth at the Eddington-limited e-folding time t_E ~ few x 0.01 billion years, is too slow to grow stellar-mass BH seeds into the supermassive luminous quasars that are observed when the universe is 1 billion years old. We propose a dynamical mechanism that can trigger supra-exponential accretion in the early universe, when a BH seed is trapped in a star cluster fed by the ubiquitous dense cold gas flows. The high gas opacity traps the accretion radiation, while the low-mass BH's random motions suppress the formation of a slowly-draining accretion disk. Supra-exponential growth can thus explain the puzzling emergence of supermassive BHs that power luminous quasars so soon after the Big Bang.

preprint2014arXiv

The Effect of Large-Scale Structure on the Magnification of High-Redshift Sources by Cluster-Lenses

Cluster gravitational lensing surveys like the Hubble Space Telescope Frontier Fields survey will detect distant galaxies 10-50 times fainter than any yet discovered. Using these surveys to measure the luminosity function of such faint, distant galaxies, however, requires that magnification maps built from the constraints of strongly-lensed images be accurate. For models that assume the cluster and nearby (correlated) structures are the only significant sources of lensing, a potential source of error in these maps comes from the fact that light rays also suffer weak deflections by uncorrelated large-scale structure along the line-of-sight, i.e. cosmic weak lensing (CWL). To demonstrate the magnitude of this effect, we calculate the magnification change which results when the same cluster-lens is placed along different lines of sight. Using a simple density profile for a cluster-lens at z~0.3-0.5 and the power spectrum of the matter density fluctuations responsible for CWL, we show that the typical magnifications of ~5(10) of sources at z=6-10 can differ by ~10-20(20-30)% from one line-of-sight to another. However, these fluctuations rise to greater than order unity near critical curves, indicating that CWL tends to make its greatest contribution to the most magnified images. We conclude that the neglect of CWL in determining the intrinsic luminosities of highly-magnified galaxies may introduce errors significant enough to warrant further effort to include this contribution in cluster-lens modeling. We suggest that methods of modeling CWL in galaxy-strong-lensing systems should be generalized to cluster-lensing systems.

preprint2014arXiv

Understanding black hole mass assembly via accretion and mergers at late times in cosmological simulations

Accretion is thought to primarily contribute to the mass accumulation history of supermassive black holes throughout cosmic time. While this may be true at high redshifts, at lower redshifts and for the most massive black holes mergers themselves might add significantly to the mass budget. We evolve SMBHs from $4 > z > 0$ using merger trees derived from hydrodynamical cosmological simulations of a cluster and void region, scaled to the observed value of the stellar mass fraction to account for overcooling. Mass gains from gas accretion proportional to bulge growth and BH-BH mergers are tracked, as are black holes that remain "orbiting" due to insufficient dynamical friction in a merger remnant, as well as those that are ejected due to gravitational recoil. We find that gas accretion remains the dominant source of mass accumulation in almost all SMBHs; mergers contribute $2.5\pm0.1\%$ for all SMBHs in the cluster and $1.0\pm0.1\%$ in the void since $z = 4$. However, mergers are significant for massive SMBHs. The fraction of mass accumulated from mergers for central BHs generally increases for larger values of the host bulge mass: in the void, the fraction is $2\%$ at $M_{*, bul} = 10^{10} M_{\odot}$, increasing to $4\%$ at $M_{*, bul} \gtrsim 10^{11} M_{\odot}$, and in the cluster it is $4\%$ at $M_{*, bul} = 10^{10} M_{\odot}$ and $23\%$ at $10^{12} M_{\odot}$. We find that $40\%$ of SMBHs and $\approx 8\%$ of the total SMBH mass is found orbiting in the cluster region at $z = 0$. The existence of orbiting and ejected SMBHs requires modification of the Soltan argument. We estimate this correction to the integrated accreted mass density of SMBHs to be in the range $6-21\%$, with a mean value of $11\pm3\%$. We also calculate the total energy output and strain from gravitational waves emitted by merging SMBHs, and obtain a signal potentially detectable by pulsar timing arrays.

preprint2013arXiv

Unravelling obese black holes in the first galaxies

We predict the existence and observational signatures of a new class of objects that assembled early, during the first billion years of cosmic time: Obese Black-hole Galaxies (OBGs). OBGs are objects in which the mass of the central black hole initially exceeds that of the stellar component of the host galaxy, and the luminosity from black-hole accretion dominates the starlight. From a cosmological simulation, we demonstrate that there are sites where star formation is initially inhibited and direct-collapse black holes (DCBHs) form due to the photo-dissociating effect of Lyman-Werner radiation on molecular hydrogen. We show that the formation of OBGs is inevitable, because the probability of finding the required extra-galactic environment and the right physical conditions in a halo conducive to DCBH formation is quite high in the early universe. We estimate an OBG number density of 0.009/Mpc^3 at z~8 and 0.03/Mpc^3 at z~6. Extrapolating from our simulation volume, we infer that the most luminous quasars detected at z~6 likely transited through an earlier OBG phase. We find that these primordial galaxies start off with an over-massive BH and acquire their stellar component from subsequent merging as well as in-situ star formation. In doing so, they inevitably go through an OBG phase dominated by the accretion luminosity at the Eddington rate or below, released from the growing BH. The OBG phase is characterised by an ultra-violet (UV) spectrum with slope of beta ~ -2.3 and the absence of a Balmer Break. OBGs should also be spatially unresolved, and are expected to be brighter than the majority of known high-redshift galaxies. OBGs could potentially be revealed via HST follow-up imaging of samples of brighter Lyman-break galaxies provided by wide-area ground-based surveys such as UltraVISTA, and should be easily uncovered and studied with instruments aboard JWST...(abridged)

preprint2012arXiv

Cluster Lenses

Clusters of galaxies are the most recently assembled, massive, bound structures in the Universe. As predicted by General Relativity, given their masses, clusters strongly deform space-time in their vicinity. Clusters act as some of the most powerful gravitational lenses in the Universe. Light rays traversing through clusters from distant sources are hence deflected, and the resulting images of these distant objects therefore appear distorted and magnified. Lensing by clusters occurs in two regimes, each with unique observational signatures. The strong lensing regime is characterized by effects readily seen by eye, namely, the production of giant arcs, multiple-images, and arclets. The weak lensing regime is characterized by small deformations in the shapes of background galaxies only detectable statistically. Cluster lenses have been exploited successfully to address several important current questions in cosmology: (i) the study of the lens(es) - understanding cluster mass distributions and issues pertaining to cluster formation and evolution, as well as constraining the nature of dark matter; (ii) the study of the lensed objects - probing the properties of the background lensed galaxy population - which is statistically at higher redshifts and of lower intrinsic luminosity thus enabling the probing of galaxy formation at the earliest times right up to the Dark Ages; and (iii) the study of the geometry of the Universe - as the strength of lensing depends on the ratios of angular diameter distances between the lens, source and observer, lens deflections are sensitive to the value of cosmological parameters and offer a powerful geometric tool to probe Dark Energy. In this review, we present the basics of cluster lensing and provide a current status report of the field.

preprint2012arXiv

The Effects of Primordial Non-Gaussianity on Giant-Arc Statistics: A Scale Dependent Example

In a recently published article, we quantified the impact of primordial non-Gaussianity on the probability of giant-arc formation. In that work, we focused on the local form of non-Gaussianity and found that it can have only a modest effect given the most recent constraints from Cosmic Microwave Background (CMB) measurements. Here, we present new calculations using a parameterization of scale-dependent non-Gaussianity in which the primordial bispectrum has the equilateral shape and the effective f_NL parameter depends on scale. We find that non-Gaussianity of this type can yield a larger effect on the giant-arc abundance compared to the local form due to both the scale dependence and the relatively weaker constraints on the equilateral shape from CMB measurements. In contrast to the maximum ~40% effect (within the latest CMB constraints) previously found for the local form, we find that the predicted giant-arc abundance for the scale-dependent equilateral form can differ by a factor of a few with respect to the Gaussian case.

preprint2012arXiv

The mass function of black holes at 1<z<4.5: comparison of models with observations

In this paper, we compare the observationally derived black hole mass function (BHMF) of luminous broad-line quasars (BLQSOs) at 1<z<4.5 drawn from the Sloan Digital Sky Survey (SDSS) presented in Kelly et al. (2010), with models of merger driven BH growth in the context of standard hierarchical structure formation models. In the models, we explore two distinct black hole seeding prescriptions at the highest redshifts: "light seeds" - remnants of Population III stars and "massive seeds" that form from the direct collapse of pre-galactic disks. The subsequent merger triggered mass build-up of the black hole population is tracked over cosmic time under the assumption of a fixed accretion rate as well as rates drawn from the distribution derived by Merloni & Heinz. Our model snapshots are compared to the SDSS derived BHMFs of BLQSOs. Our key findings are that the duty cycle of SMBHs powering BLQSOs increases with increasing redshift for all models and models with Pop III remnants as black hole seeds are unable to fit the observationally derived BHMFs for BLQSOs, lending strong support for the massive seeding model (abridged).

preprint2012arXiv

The polytropic approximation and X-ray scaling relations: constraints on gas and dark matter profiles for galaxy groups and clusters

We constrain gas and dark matter (DM) parameters of galaxy groups and clusters, by comparing X-ray scaling relations to theoretical expectations, obtained assuming that the gas is in hydrostatic equilibrium with the DM and follows a polytropic relation. We vary four parameters: the gas polytropic index Gamma, its temperature at large radii T_xi, the DM logarithmic slope at large radii zeta and its concentration c_vir. When comparing the model to the observed mass-temperature (M-T) relation of local clusters, our results are independent of both T_xi and c_vir. We thus obtain constraints on Gamma, by fixing the DM profile, and on zeta, by fixing the gas profile. For an NFW DM profile, we find that 6/5<Gamma<13/10, which is consistent with numerical simulations and observations of individual clusters. Taking 6/5<Gamma<13/10 allows the DM profile to be slightly steeper than the NFW profile at large radii. Upon including local groups, we constrain the mass-dependence of Gamma and the value of T_xi. Interestingly, with Gamma=6/5 and zeta=-3, we reproduce the observed steepening/breaking of the M-T relation at low M, if 10^6 K<T_xi<10^7 K, consistent with simulations and observations of the warm-hot intergalactic medium. When extrapolated to high redshift z, the model with a constant Gamma reproduces the expected self-similar behaviour. We also account for the observed, non-self-similar relations provided by some high-z clusters, as they provide constraints on the evolution of Gamma. Comparing our model to the observed luminosity-temperature relation, we discriminate between different M-c_vir relations: a weak dependence of c_vir on M is currently preferred by data. This simple theoretical model accounts for much of the complexity of recent, improved X-ray scaling relations, provided that we allow for a mild dependence of Gamma on M or for T_xi consistent with intercluster values. [abridged]

preprint2011arXiv

Black hole growth in the early Universe is self-regulated and largely hidden from view

The formation of the first massive objects in the infant Universe remains impossible to observe directly and yet it sets the stage for the subsequent evolution of galaxies. While some black holes with masses > billion solar masses? have been detected in luminous quasars less than one billion years after the Big Bang, these individual extreme objects have limited utility in constraining the channels of formation of the earliest black holes. The initial conditions of black hole seed properties are quickly erased during the growth process. From deep, optimally stacked, archival X-ray observations, we measure the amount of black hole growth in z=6-8 galaxies (0.7-1 billion years after the Big Bang). Our results imply that black holes grow in tandem with their hosts throughout cosmic history, starting from the earliest times. We find that most copiously accreting black holes at these epochs are buried in significant amounts of gas and dust that absorb most radiation except for the highest energy X-rays. This suggests that black holes grow significantly more than previously thought during these early bursts, and due to obscuration they do not contribute to the re-ionization of the Universe with their ultraviolet emission.

preprint2011arXiv

Combined strong and weak lensing analysis of 28 clusters from the Sloan Giant Arcs Survey

We study the mass distribution of a sample of 28 galaxy clusters using strong and weak lensing observations. The clusters are selected via their strong lensing properties as part of the Sloan Giant Arcs Survey (SGAS) from the Sloan Digital Sky Survey (SDSS). Mass modelling of the strong lensing information from the giant arcs is combined with weak lensing measurements from deep Subaru/Suprime-cam images to primarily obtain robust constraints on the concentration parameter and the shape of the mass distribution. We find that the concentration c_vir is a steep function of the mass, c_vir \propto M_vir^-0.59\pm0.12, with the value roughly consistent with the lensing-bias-corrected theoretical expectation for high mass (10^15 h^-1 M_sun) clusters. However, the observationally inferred concentration parameters appear to be much higher at lower masses (10^14 h^-1 M_sun), possibly a consequence of the modification to the inner density profiles provided by baryon cooling. The steep mass-concentration relation is also supported from direct stacking analysis of the tangential shear profiles. In addition, we explore the two-dimensional shape of the projected mass distribution by stacking weak lensing shear maps of individual clusters with prior information on the position angle from strong lens modelling, and find significant evidence for a large mean ellipticity with the best-fit value of e = 0.47 \pm 0.06 for the mass distribution of the stacked sample. We find that the luminous cluster member galaxy distribution traces the overall mass distribution very well, although the distribution of fainter cluster galaxies appears to be more extended than the total mass.

preprint2011arXiv

Cosmography with cluster strong lenses: the influence of substructure and line-of-sight halos

We explore the use of strong lensing by galaxy clusters to constrain the dark energy equation of state and its possible time variation. The cores of massive clusters often contain several multiply imaged systems of background galaxies at different redshifts. The locations of lensed images can be used to constrain cosmological parameters due to their dependence on the ratio of angular diameter distances. We employ Monte-Carlo simulations of cluster lenses, including the contribution from substructures, to assess the feasibility of this potentially powerful technique. At the present, parametric lens models use well motivated scaling relations between mass and light to incorporate cluster member galaxies, and do not explicitly model line-of-sight structure. Here, we quantify modeling errors due to scatter in the cluster galaxy scaling relations and un-modeled line-of-sight halos. These errors are of the order of a few arcseconds on average for clusters located at typical redshifts (z ~ 0.2 - 0.3). Using Bayesian Markov Chain Monte-Carlo techniques, we show that the inclusion of these modeling errors is critical to deriving unbiased constraints on dark energy. However, when the uncertainties are properly quantified, we show that constraints competitive with other methods may be obtained by combining results from a sample of just 10 simulated clusters with 20 families each. Cosmography with a set of well studied cluster lenses may provide a powerful complementary probe of the dark energy equation of state. Our simulations provide a convenient method of quantifying modeling errors and assessing future strong lensing survey strategies.

preprint2011arXiv

Evidence for Three Accreting Black Holes in a Galaxy at z~1.35: A Snapshot of Recently Formed Black Hole Seeds?

One of the key open questions in cosmology today pertains to understanding when, where and how super massive black holes form, while it is clear that mergers likely play a significant role in the growth cycles of black holes, how supermassive black holes form, and how galaxies grow around them. Here, we present Hubble Space Telescope WFC3/IR grism observations of a clumpy galaxy at z=1.35, with evidence for 10^6 - 10^7 Msun rapidly growing black holes in separate sub-components of the host galaxy. These black holes could have been brought into close proximity as a consequence of a rare multiple galaxy merger or they could have formed in situ. Such holes would eventually merge into a central black hole as the stellar clumps/components presumably coalesce to form a galaxy bulge. If we are witnessing the in-situ formation of multiple black holes, their properties can inform seed formation models and raise the possibility that massive black holes can continue to emerge in star-forming galaxies as late as z=1.35 (4.8 Gyr after the Big Bang).

preprint2011arXiv

How important is the dark matter halo for black hole growth?

In this paper, we examine if the properties of central black holes in galactic nuclei correlate with their host dark matter halo. We analyze the entire sample of galaxies where black hole mass, velocity dispersion, sigma, and asymptotic circular velocity, Vc, have all been measured. We fit M-sigma and M-Vc to a power law, and find that in both relationships the scatter and slope are similar. This model-independent analysis suggests that although the black hole masses are not uniquely determined by dark matter halo mass, when considered for the current sample as a whole, the M-Vc correlation may be as strong (or as weak) as M-sigma. Although the data are sparse, there appears to be more scatter in the correlation for both sigma and Vc at the low--mass end. This is not unexpected given our current understanding of galaxy and black hole assembly. In fact, there are several compelling reasons that account for this: (i) SMBH formation is likely less efficient in low-mass galaxies with large angular momentum content; (ii) SMBH growth is less efficient in low-mass disk galaxies that have not experienced major mergers; (iii) dynamical effects, such as gravitational recoil, increase scatter preferentially at the low-mass end. Therefore, the recent observational claim of the absence of central SMBHs in bulgeless, low mass galaxies, or deviations from the correlations defined by high-mass black holes in large galaxies today is, in fact, predicated by current models of black hole growth. We show how this arises as a direct consequence of the coupling between dark matter halos and central black holes at the earliest epochs.

preprint2011arXiv

The effects of primordial non-Gaussianity on giant-arc statistics

For over a decade, it has been debated whether the concordance LCDM model is consistent with the observed abundance of giant arcs in clusters. While previous theoretical studies have focused on properties of the lens and source populations, as well as cosmological effects such as dark energy, the impact of initial conditions on the giant-arc abundance is relatively unexplored. Here, we quantify the impact of non-Gaussian initial conditions with the local bispectrum shape on the predicted frequency of giant arcs. Using a path-integral formulation of the excursion set formalism, we extend a semi-analytic model for calculating halo concentrations to the case of primordial non-Gaussianity, which may be useful for applications outside of this work. We find that massive halos tend to collapse earlier in models with positive f_NL, relative to the Gaussian case, leading to enhanced concentration parameters. The converse is true for f_NL < 0. In addition to these effects, which change the lensing cross sections, non-Gaussianity also modifies the abundance of supercritical clusters available for lensing. These combined effects work together to either enhance (f_NL > 0) or suppress (f_NL < 0) the probability of giant-arc formation. Using the best value and 95% confidence levels currently available from the Wilkinson Microwave Anisotropy Probe, we find that the giant-arc optical depth for sources at z_s~2 is enhanced by ~20% and ~45% for f_NL = 32 and 74 respectively. In contrast, we calculate a suppression of ~5% for f_NL = -10. These differences translate to similar relative changes in the predicted all-sky number of giant arcs.

preprint2011arXiv

The formation and evolution of massive black hole seeds in the early Universe

Tracking the evolution of high redshift seed black hole masses to late times, we examine the observable signatures today. These massive initial black hole seeds form at extremely high redshifts from the direct collapse of pre-galactic gas discs. Populating dark matter halos with seeds formed in this fashion, we follow the mass assembly history of these black holes to the present time using a Monte-Carlo merger tree approach. Utilizing this formalism, we predict the black hole mass function at high redshifts and at the present time; the integrated mass density of black holes in the Universe; the luminosity function of accreting black holes as a function of redshift and the scatter in observed, local M-sigma relation. Comparing the predictions of the `light' seed model with these massive seeds we find that significant differences appear predominantly at the low mass end of the present day black hole mass function. However, all our models predict that low surface brightness, bulge-less galaxies with large discs are least likely to be sites for the formation of massive seed black holes at high redshifts. The efficiency of seed formation at high redshifts has a direct influence on the black hole occupation fraction in galaxies at z=0. This effect is more pronounced for low mass galaxies. This is the key discriminant between the models studied here and the Population III remnant `light' seed model. We find that there exists a population of low mass galaxies that do not host nuclear black holes. Our prediction of the shape of the M-sigma relation at the low mass end and increased scatter has recently been corroborated by observations.

preprint2011arXiv

The mass assembly history of black holes in the Universe

We track the growth and evolution of high redshift seed black holes over cosmic time. This population of massive, initial black hole seeds form at these early epochs from the direct collapse of pre-galactic gas discs. Populating dark matter halos with seeds formed in this fashion, we follow their mass assembly history to the present time using a Monte-Carlo merger tree approach. Using this formalism, we predict the black hole mass function at the present time; the integrated mass density of black holes in the Universe; the luminosity function of accreting black holes as a function of redshift and the scatter in observed, local Mbhσs relation. Signatures of these massive seed models appear predominantly at the low mass end of the present day black hole mass function. In fact, our prediction of the shape of the Mbhσs relation at the low mass end and increased scatter has recently been corroborated by observations. These models predict that low surface brightness, bulge-less galaxies with large discs are least likely to be sites for the formation of massive seed black holes at high redshifts. The efficiency of seed formation at high redshifts also has a direct influence on the black hole occupation fraction in galaxies at z = 0. This effect is more pronounced for low mass galaxies today as we predict the existence of a population of low mass galaxies that do not host nuclear black holes. This is the key discriminant between the models studied here and the Population-III remnant seed model.

preprint2011arXiv

Virialization of high redshift dark matter haloes

We present results of a study of the virial state of high redshift dark matter haloes in an N-body simulation. We find that the majority of collapsed, bound haloes are not virialized at any redshift slice in our study ($z=15-6$) and have excess kinetic energy. At these redshifts, merging is still rampant and the haloes cannot strictly be treated as isolated systems. To assess if this excess kinetic energy arises from the environment, we include the surface pressure term in the virial equation explicitly and relax the assumption that the density at the halo boundary is zero. Upon inclusion of the surface term, we find that the haloes are much closer to virialization, however, they still have some excess kinetic energy. We report trends of the virial ratio including the extra surface term with three key halo properties: spin, environment, and concentration. We find that haloes with closer neighbors depart more from virialization, and that haloes with larger spin parameters do as well. We conclude that except at the lowest masses ($M < 10^6 \Msun$), dark matter haloes at high redshift are not fully virialized. This finding has interesting implications for galaxy formation at these high redshifts, as the excess kinetic energy will impact the subsequent collapse of baryons and the formation of the first disks and/or baryonic structures.

preprint2010arXiv

Cosmological Constraints from Strong Gravitational Lensing in Galaxy Clusters

Current efforts in observational cosmology are focused on characterizing the mass-energy content of the Universe. We present results from a geometric test based on strong lensing in galaxy clusters. Based on Hubble Space Telescope images and extensive ground-based spectroscopic follow-up of the massive galaxy cluster Abell 1689, we used a parametric model to simultaneously constrain the cluster mass distribution and dark energy equation of state. Combining our cosmological constraints with those from X-ray clusters and the Wilkinson Microwave Anisotropy Probe 5-year data gives Ωm = 0.25 +/- 0.05 and wx = -0.97 +/- 0.07 which are consistent with results from other methods. Inclusion of our method with all other techniques available brings down the current 2σ contours on the dark energy equation of state parameter wx by about 30%.

preprint2010arXiv

Hydrostatic equilibrium profiles for gas in elliptical galaxies

We present an analytic formulation for the equilibrium gas density profile of early-type galaxies that explicitly includes the contribution of stars in the gravitational potential. We build a realistic model for an isolated elliptical galaxy and explore the equilibrium gas configurations as a function of multiple parameters. For an assumed central gas temperature k_B*T_0=0.6 keV, we find that neglecting the gravitational effects of stars, which can contribute substantially in the innermost regions, leads to an underestimate of the enclosed baryonic gas mass by up to ~65% at the effective radius, and by up to ~15% at the NFW scale radius, depending on the stellar baryon fraction. This formula is therefore important for estimating the baryon fraction in an unbiased fashion. These new hydrostatic equilibrium solutions, derived for the isothermal and polytropic cases, can also be used to generate more realistic initial conditions for simulations of elliptical galaxies. Moreover, the new formulation is relevant when interpreting X-ray data. We compare our composite isothermal model to the standard beta-model used to fit X-ray observations of early-type galaxies, to determine the value of the NFW scale radius r_s. Assuming a 10% stellar baryon fraction, we find that the exclusion of stars from the gravitational potential leads to (i) an underestimate of r_s by ~80%, and to (ii) an overestimate of the enclosed dark matter at r_s by a factor of ~2, compared to the equivalent beta-model fit results when stars are not taken into account. For higher stellar mass fractions, a beta-model is unable to accurately reproduce our solution, indicating that when the observed surface brightness profile of an isolated elliptical galaxy is found to be well fitted by a beta-model, the stellar mass fraction cannot be much greater than ~10%.

preprint2010arXiv

Major Galaxy Mergers and the Growth of Supermassive Black Holes in Quasars

Despite observed strong correlations between central supermassive black holes (SMBHs) and star-formation in galactic nuclei, uncertainties exist in our understanding of their coupling. We present observations of the ratio of heavily-obscured to unobscured quasars as a function of cosmic epoch up to z~3, and show that a simple physical model describing mergers of massive, gas-rich galaxies matches these observations. In the context of this model, every obscured and unobscured quasar represent two distinct phases that result from a massive galaxy merger event. Much of the mass growth of the SMBH occurs during the heavily-obscured phase. These observations provide additional evidence for a causal link between gas-rich galaxy mergers, accretion onto the nuclear SMBH and coeval star formation.

preprint2010arXiv

Observed Scaling Relations for Strong Lensing Clusters: Consequences for Cosmology and Cluster Assembly

Scaling relations of observed galaxy cluster properties are useful tools for constraining cosmological parameters as well as cluster formation histories. One of the key cosmological parameters, sigma8, is constrained using observed clusters of galaxies, although current estimates of sigma8 from the scaling relations of dynamically relaxed galaxy clusters are limited by the large scatter in the observed cluster mass-temperature (M-T) relation. With a sample of eight strong lensing clusters at 0.3 < z <0.8, we find that the observed cluster concentration-mass relation can be used to reduce the M-T scatter by a factor of 6. Typically only relaxed clusters are used to estimate sigma8, but combining the cluster concentration-mass relation with the M-T relation enables the inclusion of unrelaxed clusters as well. Thus, the resultant gains in the accuracy of sigma8 measurements from clusters are twofold: the errors on sigma8 are reduced and the cluster sample size is increased. Therefore, the statistics on sigma8 determination from clusters are greatly improved by the inclusion of unrelaxed clusters. Exploring cluster scaling relations further, we find that the correlation between brightest cluster galaxy (BCG) luminosity and cluster mass offers insight into the assembly histories of clusters. We find preliminary evidence for a steeper BCG luminosity - cluster mass relation for strong lensing clusters than the general cluster population, hinting that strong lensing clusters may have had more active merging histories.

preprint2010arXiv

Spin and structural halo properties at high redshift in a LCDM Universe

In this paper, we examine in detail the key structural properties of high redshift dark matter haloes as a function of their spin parameter. We perform and analyze high resolution cosmological simulations of the formation of structure in a LCDM Universe. We study the mass function, ellipticities, shapes, density profiles, rotation curves and virialization for a large sample of dark matter haloes from z = 15 - 6. We also present detailed convergence tests for individual haloes. We find that high spin haloes have stronger clustering strengths (up to 25%) at all mass and redshift ranges at these early epochs. High redshift spherical haloes are also up to 50% more clustered than aspherical haloes. High spin haloes at these redshifts are also preferentially found in high density environments, and have more neighbors than their low spin counterparts. We report a systematic offset in the peak of the circular velocity curves for high and low spin haloes of the same mass. Therefore, estimating halo masses without knowledge of the spin, using only the circular velocity can yield errors of up to 40%. The strong dependence of key structural properties on spin that we report here likely have important implications for studies of star formation and feedback from these galaxies.

preprint2010arXiv

The Sudden Death of the Nearest Quasar

Galaxy formation is significantly modulated by energy output from supermassive black holes at the centers of galaxies which grow in highly efficient luminous quasar phases. The timescale on which black holes transition into and out of such phases is, however, unknown. We present the first measurement of the shutdown timescale for an individual quasar using X-ray observations of the nearby galaxy IC 2497, which hosted a luminous quasar no more than 70,000 years ago that is still seen as a light echo in `Hanny's Voorwerp', but whose present-day radiative output is lower by at least 2 and more likely by over 4 orders of magnitude. This extremely rapid shutdown provides new insights into the physics of accretion in supermassive black holes, and may signal a transition of the accretion disk to a radiatively inefficient state.

preprint2009arXiv

Cosmography with cluster strong lensing

By stacking an ensemble of strong lensing clusters, we demonstrate the feasibility of placing constraints on the dark energy equation of state. This is achieved by using multiple images of sources at two or more distinct redshift planes. The sample of smooth clusters in our simulations is based on observations of massive clusters and the distribution of background galaxies is constructed using the Hubble Deep Field. Our source distribution reproduces the observed redshift distribution of multiply imaged sources in Abell 1689. The cosmology recovery depends on the number of image families with known spectroscopic redshifts and the number of stacked clusters. Our simulations suggest that constraints comparable to those derived from other competing established techniques on a constant dark energy equation of state can be obtained using 10 to 40 clusters with 5 or more families of multiple images. We have also studied the observational errors in the image redshifts and positions. We find that spectroscopic redshifts and high resolution {\it Hubble Space Telescope} images are required to eliminate confidence contour relaxation relative to the ideal case in our simulations. This suggests that the dark energy equation of state, and other cosmological parameters, can be constrained with existing {\it Hubble Space Telescope} images of lensing clusters coupled with dedicated ground-based arc spectroscopy.

preprint2009arXiv

Journey to the M_BH -sigma relation: the fate of low mass black holes in the Universe

In this paper, we explore the establishment and evolution of the empirical correlation between black hole mass and velocity dispersion with redshift. We track the growth and accretion history of massive black holes starting from high redshift using two seeding models:(i) Population III remnants, and (ii) massive seeds from direct gas collapse. Although the seeds do not initially satisfy the $M_{\rm BH} - σ$ relation, the correlation is established and maintained at all times if self-regulating accretion episodes are associated with major mergers. The massive end of the $M_{\rm BH} - σ$ relation is established early, and lower mass MBHs migrate over time. How MBHs migrate toward the relation, the slope and the scatter of the relation all depend critically on the seeding model as well as the adopted self-regulation prescription. We expect flux limited AGN surveys and LISA to select accreting and merging MBHs respectively that have already migrated onto the $\msigma$ relation. This is a consequence of major mergers being more common at high redshift for the most massive, biased, galaxies that anchor the $\msigma$ relation early. We also predict the existence of a large population of low mass `hidden' MBHs at high redshift which can easily escape detection. Additionally, we find that if MBH seeds are massive, $\sim 10^5 M_{\odot}$, the low-mass end of the $\msigma$ flattens towards this asymptotic value, creating a characteristic `plume'.

preprint2009arXiv

Subaru Weak Lensing Measurements of Four Strong Lensing Clusters: Are Lensing Clusters Over-Concentrated?

We derive radial mass profiles of four strong lensing selected clusters which show prominent giant arcs (Abell 1703, SDSS J1446+3032, SDSS J1531+3414, and SDSS J2111-0115), by combining detailed strong lens modeling with weak lensing shear measured from deep Subaru Suprime-cam images. Weak lensing signals are detected at high significance for all four clusters, whose redshifts range from z=0.28 to 0.64. We demonstrate that adding strong lensing information with known arc redshifts significantly improves constraints on the mass density profile, compared to those obtained from weak lensing alone. While the mass profiles are well fitted by the universal form predicted in N-body simulations of the LCDM model, all four clusters appear to be slightly more centrally concentrated (the concentration parameters c_vir \sim 8) than theoretical predictions, even after accounting for the bias toward higher concentrations inherent in lensing selected samples. Our results are consistent with previous studies which similarly detected a concentration excess, and increases the total number of clusters studied with the combined strong and weak lensing technique to ten. Combining our sample with previous work, we find that clusters with larger Einstein radii are more anomalously concentrated. We also present a detailed model of the lensing cluster A1703 with constraints from multiple image families, and find the dark matter inner density profile to be cuspy with the slope consistent with -1, in agreement with expectations.

preprint2004arXiv

Modeling the accretion history of supermassive black holes

There is overwhelming evidence for the presence of supermassive black holes (SMBHs) in the centers of most nearby galaxies. The mass estimates for these remnant black holes from the stellar kinematics of local galaxies and the quasar phenomenon at high redshifts point to the presence of assembled SMBHs. The accretion history of SMBHs can be reconstructed using observations at high and low redshifts as model constraints. Observations of galaxies and quasars in the submillimeter, infrared, optical, and X-ray wavebands are used as constraints, along with data from the demography of local black holes. Theoretical modeling of the growth of black hole mass with cosmic time has been pursued thus far in two distinct directions: a phenomenological approach that utilizes observations in various wavebands, and a semi-analytic approach that starts with a theoretical framework and a set of assumptions with a view to matching observations. Both techniques have been pursued in the context of the standard paradigm for structure formation in a Cold Dark Matter dominated universe. Here, we examine the key issues and uncertainties in the theoretical understanding of the growth of SMBHs.

preprint2002arXiv

Accretion during the merger of supermassive black holes

We study the evolution of disk accretion during the merger of supermassive black hole binaries in galactic nuclei. In hierarchical galaxy formation models, the most common binaries are likely to arise from minor galactic mergers, and have unequal mass black holes. Once such a binary becomes embedded in an accretion disk at a separation of the order of 0.1 pc, the merger proceeds in two distinct phases. During the first phase, the loss of orbital angular momentum to the gaseous disk shrinks the binary on a timescale of ~10 Myr. The accretion rate onto the primary black hole is not increased, and can be substantially reduced, during this disk-driven migration. At smaller separations, gravitational radiation becomes the dominant angular momentum loss process, and any gas trapped inside the orbit of the secondary is driven inwards by the inspiralling black hole. The implied accretion rate just prior to coalescence exceeds the Eddington limit, so the final merger is likely to occur within a common envelope formed from the disrupted inner disk, and be accompanied by high velocity (~10,000 km/s) outflows.

preprint2000arXiv

Discriminating weak lensing from intrinsic spin correlations using the curl-gradient decomposition

The distortion field defined by the ellipticities of galaxy shapes projected on the sky can be uniquely decomposed into a gradient and a curl component. If the observed ellipticities are induced by weak gravitational lensing, then the distortion field is curl free. Here we show that, in contrast, the distortion field resulting from intrinsic spin alignments is not curl free. This provides a powerful discriminant between lensing and intrinsic contributions to observed ellipticity correlations. We also show how these contributions can be disentangled statistically from the ellipticity correlations or computed locally from circular integrals of the ellipticity field. This allows for an unambiguous detection of intrinsic galaxy alignments in the data. When the distortions are dominated by lensing, as occurs at high redshifts, the decomposition provides a valuable tool for understanding properties of the noise and systematic errors. These techniques can be applied equally well to the polarization of the microwave background, where it can be used to separate curl-free scalar perturbations from those produced by gravity waves or defects.

preprint1997arXiv

Estimating the mass density of neutral gas at $z < 1$

We use the relationships between galactic HI mass and B-band luminosity determined by Rao & Briggs to recalculate the mass density of neutral gas at the present epoch based on more recent measures of the galaxy luminosity function than were available to those authors. We find $Ω_{gas}(z=0) \simeq 5 \times 10^{-4}$ in good agreement with the original Rao & Briggs value, suggesting that this quantity is now reasonably secure. We then show that, if the scaling between H I mass and B-band luminosity has remained approximately constant since $z = 1$, the evolution of the luminosity function found by the Canada-France redshift survey translates to an increase of $Ω_{gas}$ by a factor of $\approx 3$ at $z = 0.5 - 1$ . A similar value is obtained quite independently from consideration of the luminosity function of Mg II absorbers at $z = 0.65$. By combining these new estimates with data from damped \lya systems at higher redshift, it is possible to assemble a rough sketch of the evolution of $Ω_{gas}$ over the last 90% of the age of the universe. The consumption of H I gas with time is in broad agreement with models of chemical evolution which include the effects of dust, although more extensive samples of damped \lya systems at low and intermediate redshift are required for a quantitative assessment of the dust bias.

preprint1997arXiv

High-redshift galaxies, their active nuclei and central black holes

We demonstrate that the luminosity function of the recently detected population of star-forming galaxies and the QSO luminosity function at z=3 can be matched with the mass function of dark matter haloes predicted by hierarchical cosmogonies for lifetimes of optically bright QSOs in the range 10^6 to 10^8 yr. We suggest that the mass of supermassive black holes may be limited by the back-reaction of the emitted energy on the accretion flow in a self-gravitating disc. This would imply a relation of black hole to halo mass of the form M_bh ~(v_halo)^5 and a typical duration of the optically bright QSO phase of a few times 10^7 yr. The high integrated mass density of black holes inferred from recent black hole mass estimates in nearby galaxies may indicate that the overall efficiency of supermassive black holes for producing blue light is smaller than previously assumed. We discuss three possible accretion modes with low optical emission efficiency: (i) accretion at far above the Eddington rate, (ii) accretion obscured by dust, and (iii) accretion below the critical rate leading to an advection dominated accretion flow lasting for a Hubble time. We further argue that accretion with low optical efficiency might be closely related to the origin of the hard X-ray background and that the ionizing background might be progressively dominated by stars rather than QSOs at higher redshift.

preprint1997arXiv

Lensing by galaxy halos in clusters of galaxies

Weak shear maps of the outer regions of clusters have been successfully used to map the distribution of mass at large radii from the cluster center. The typical smoothing lengths employed thus far preclude the systematic study of the effects of galactic-scale substructure on the measured weak lensing signal. In this paper, we present two methods to infer the possible existence and extent of dark halos around bright cluster galaxies by quantifying the `local' weak lensing induced by them. The proposed methods are: direct radial averaging of the shear field in the vicinity of bright cluster members and a maximum-likelihood method to extract fiducial parameters characterizing galaxy halos. We demonstrate using simulations that these observed local weak-shear effects on galaxy scales within the cluster can be used to statistically constrain the mean mass-to-light ratio, and fiducial parameters like the halo size, velocity dispersion and hence mass of cluster galaxies. We compare the two methods and investigate their relative drawbacks and merits in the context of feasibility of application to HST cluster data, whereby we find that the prospects are promising for detection on stacking a minimum of 20 WFPC2 deep cluster fields.

preprint1996arXiv

Distribution functions for clusters of galaxies from N-body simulations

We present the results of an attempt to adapt the distribution function formalism to characterize large-scale structures like clusters of galaxies that form in a cosmological N-body simulation. While galaxy clusters are systems that are not strictly in equilibrium, we show that their evolution can nevertheless be studied using a physically motivated extension of the language of equilibrium stellar dynamics. Restricting our analysis to the virialized region, a prescription to limit the accessible phase-space is presented, which permits the construction of both the isotropic and the anisotropic distribution functions $f(E)$ and $f(E,L)$. The method is applied to models extracted from a catalogue of simulated clusters. Clusters evolved in open and flat background cosmologies are followed during the course of their evolution, and are found to transit through a sequence of what we define as `quasi-equilibrium' states. An interesting feature is that the computed $f(E)$ is well fit by an exponential form. We conclude that the dynamical evolution of a cluster, undergoing relaxation punctuated by interactions and violent mergers with consequent energy-exchange, can be studied both in a qualitative and quantitative fashion by following the time evolution of $f(E)$.

preprint1996arXiv

Probing the dynamics of cluster-lenses

We propose a new approach to study the dynamical implications of mass models of clusters for the velocity structure of galaxies in the core. Strong and weak lensing data are used to construct the total mass profile of the cluster, which is used in conjunction with the optical galaxy data to solve in detail for the nature of galaxy orbits and the velocity anisotropy in the central regions. We also examine other observationally and physically motivated mass models, specifically those obtained from X-ray observations and N-body simulations. The aim of this analysis is to understand qualitatively the structure of the core and test some of the key assumptions of the standard picture of cluster formation regarding relaxation, virialization and equilibrium. This technique is applied to the cluster Abell 2218, where we find evidence for an anisotropic core, which we interpret to indicate the existence of a dynamically disturbed central region. We find that the requirement of physically meaningful solutions for the velocity anisotropy places stringent bounds on the slope of cluster density profiles in the inner regions.