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

Brandon C. Kelly

Brandon C. Kelly 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.GAastro-ph.HEastro-ph.IMastro-phphysics.data-an

Catalog footprint

What is connected

32works
6topics
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

32 published item(s)

preprint2014arXiv

Data Mining for Gravitationally Lensed Quasars

Gravitationally lensed (GL) quasars are brighter than their unlensed counterparts and produce images with distinctive morphological signatures. Past searches and target selection algorithms, in particular the Sloan Quasar Lens Search (SQLS), have relied on basic morphological criteria, which were applied to samples of bright, spectroscopically confirmed quasars. The SQLS techniques are not sufficient for searching into new surveys (e.g. DES, PS1, LSST), because spectroscopic information is not readily available and the large data volume requires higher purity in target/candidate selection. We carry out a systematic exploration of machine learning techniques and demonstrate that a two step strategy can be highly effective. In the first step we use catalog-level information ($griz$+WISE magnitudes, second moments) to preselect targets, using artificial neural networks. The accepted targets are then inspected with pixel-by-pixel pattern recognition algorithms (Gradient-Boosted Trees), to form a final set of candidates. The results from this procedure can be used to further refine the simpler SQLS algorithms, with a twofold (or threefold) gain in purity and the same (or $80\%$) completeness at target-selection stage, or a purity of $70\%$ and a completeness of $60\%$ after the candidate-selection step. Simpler photometric searches in $griz$+WISE based on colour cuts would provide samples with $7\%$ purity or less. Our technique is extremely fast, as a list of candidates can be obtained from a stage III experiment (e.g. DES catalog/database) in {a few} CPU hours. The techniques are easily extendable to Stage IV experiments like LSST with the addition of time domain information.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Flexible and Scalable Methods for Quantifying Stochastic Variability in the Era of Massive Time-Domain Astronomical Data Sets

We present the use of continuous-time autoregressive moving average (CARMA) models as a method for estimating the variability features of a light curve, and in particular its power spectral density (PSD). CARMA models fully account for irregular sampling and measurement errors, making them valuable for quantifying variability, forecasting and interpolating light curves, and for variability-based classification. We show that the PSD of a CARMA model can be expressed as a sum of Lorentzian functions, which makes them extremely flexible and able to model a broad range of PSDs. We present the likelihood function for light curves sampled from CARMA processes, placing them on a statistically rigorous foundation, and we present a Bayesian method to infer the probability distribution of the PSD given the measured lightcurve. Because calculation of the likelihood function scales linearly with the number of data points, CARMA modeling scales to current and future massive time-domain data sets. We conclude by applying our CARMA modeling approach to light curves for an X-ray binary, two AGN, a long-period variable star, and an RR-Lyrae star, in order to illustrate their use, applicability, and interpretation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stochastic Modeling of the Fermi/LAT Gamma-ray Blazar Variability

We study the gamma-ray variability of 13 blazars observed with the Fermi Large Area Telescope (LAT). These blazars have the most complete light curves collected during the first 4 years of the Fermi sky survey. We model them with the Ornstein-Uhlenbeck (OU) process or a mixture of the OU processes. The OU process has power spectral density (PSD) proportional to 1/f^alpha with alpha changing at a characteristic time scale, tau_0, from 0 (tau>>tau_0) to 2 (tau<<tau_0). The PSD of the mixed OU process has two characteristic time scales and an additional intermediate region with 0<alpha<2. We show that the OU model provides a good description of the Fermi/LAT light curves of three blazars in our sample. For the first time we constrain a characteristic gamma-ray time scale of variability in two BL Lac sources, 3C 66A and PKS 2155-304 (tau_0=25 day and tau_0=43 day, respectively, in the observer's frame), which are longer than the soft X-ray time scales detected in blazars and Seyfert galaxies. We find that the mixed OU process approximates the light curves of the remaining 10 blazars better than the OU process. We derive limits on their long and short characteristic time scales, and infer that their Fermi/LAT PSDs resemble a power-law function. We constrain the PSD slopes for all but one source in the sample. We find hints for sub-hour Fermi/LAT variability in four flat spectrum radio quasars. We discuss the implications of our results for theoretical models of blazar variability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Luminosity Function at z~8 from 97 Y-band dropouts: Inferences About Reionization

[Abbreviated] We present the largest search to date for $z\sim8$ Lyman break galaxies (LBGs) based on 350 arcmin$^2$ of HST observations in the V-, Y-, J- and H-bands from the Brightest of Reionizing Galaxies (BoRG) survey. The BoRG dataset includes $\sim$50 arcmin$^2$ of new data and deeper observations of two previous BoRG pointings, from which we present 9 new $z\sim8$ LBG candidates, bringing the total number of BoRG LBGs to 38 with $25.5\leqslant m_{J} \leqslant 27.6$ (AB system). We introduce a new Bayesian formalism for estimating the galaxy luminosity function (LF), which does not require binning (and thus smearing) of the data and includes a likelihood based on the formally correct binomial distribution as opposed to the often used approximate Poisson distribution. We demonstrate the utility of the new method on a sample of $97$ LBGs that combines the bright BoRG galaxies with the fainter sources published in Bouwens et al. (2012) from the HUDF and ERS programs. We show that the $z\sim8$ LF is well described by a Schechter function with a characteristic magnitude $M^\star = -20.15^{+0.29}_{-0.38}$, a faint-end slope of $α= -1.87^{+0.26}_{-0.26}$, and a number density of $\log_{10} ϕ^\star [\textrm{Mpc}^{-3}] = -3.24^{+0.25}_{-0.24}$. Integrated down to $M=-17.7$ this LF yields a luminosity density, $\log_{10} ε[\textrm{erg}/\textrm{s/Hz/Mpc}^{3}] = 25.52^{+0.05}_{-0.05}$. Our LF analysis is consistent with previously published determinations within 1$σ$. We discuss the implication of our study for the physics of reionization. By assuming theoretically motivated priors on the clumping factor and the photon escape fraction we show that the UV LF from galaxy samples down to $M=-17.7$ can ionize only 10-50% of the neutral hydrogen at $z\sim8$. Full reionization would require extending the LF down to $M=-15$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Technical Overview

The Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM) is a dedicated multi-object RM experiment that has spectroscopically monitored a sample of 849 broad-line quasars in a single 7 deg$^2$ field with the SDSS-III BOSS spectrograph. The RM quasar sample is flux-limited to i_psf=21.7 mag, and covers a redshift range of 0.1<z<4.5. Optical spectroscopy was performed during 2014 Jan-Jul dark/grey time, with an average cadence of ~4 days, totaling more than 30 epochs. Supporting photometric monitoring in the g and i bands was conducted at multiple facilities including the CFHT and the Steward Observatory Bok telescopes in 2014, with a cadence of ~2 days and covering all lunar phases. The RM field (RA, DEC=14:14:49.00, +53:05:00.0) lies within the CFHT-LS W3 field, and coincides with the Pan-STARRS 1 (PS1) Medium Deep Field MD07, with three prior years of multi-band PS1 light curves. The SDSS-RM 6-month baseline program aims to detect time lags between the quasar continuum and broad line region (BLR) variability on timescales of up to several months (in the observed frame) for ~10% of the sample, and to anchor the time baseline for continued monitoring in the future to detect lags on longer timescales and at higher redshift. SDSS-RM is the first major program to systematically explore the potential of RM for broad-line quasars at z>0.3, and will investigate the prospects of RM with all major broad lines covered in optical spectroscopy. SDSS-RM will provide guidance on future multi-object RM campaigns on larger scales, and is aiming to deliver more than tens of BLR lag detections for a homogeneous sample of quasars. We describe the motivation, design and implementation of this program, and outline the science impact expected from the resulting data for RM and general quasar science.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

X-ray constraints on the local supermassive black hole occupation fraction

Distinct seed formation mechanisms are imprinted upon the fraction of dwarf galaxies currently containing a central supermassive black hole. Seeding by Pop III remnants is expected to produce a higher occupation fraction than is generated with direct gas collapse precursors. Chandra observations of nearby early-type galaxies can directly detect even low-level supermassive black hole activity, and the active fraction immediately provides a firm lower limit to the occupation fraction. Here, we use the volume-limited AMUSE surveys of ~200 optically-selected early-type galaxies to characterize simultaneously, for the first time, the occupation fraction and the scaling of nuclear X-ray luminosity with stellar mass, accounting for intrinsic scatter, measurement uncertainties, and X-ray limits. For early-type galaxies with log(M_star/M_sun)<10, we obtain a lower limit to the occupation fraction of >20% (at 95% confidence), but full occupation cannot be excluded. The preferred dependence of log(L_X) upon log(M_star) has a slope of about 0.7-0.8, consistent with the "downsizing" trend previously identified from the AMUSE dataset, and a uniform Eddington efficiency is disfavored at ~2 sigma. We provide guidelines for the future precision with which these parameters may be refined with larger or more sensitive samples.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Quasar-Galaxy Mixing Diagram: Quasar Spectral Energy Distribution Shapes in the Optical to Near-Infrared

We define a quasar-galaxy mixing diagram using the slopes of their spectral energy distributions (SEDs) from 1μm to 3000Å and from 1μm to 3μm in the rest frame. The mixing diagram can easily distinguish among quasar-dominated, galaxy-dominated and reddening-dominated SED shapes. By studying the position of the 413 XMM selected Type 1 AGN in the wide-field "Cosmic Evolution Survey" (COSMOS) in the mixing diagram, we find that a combination of the Elvis et al. (1994, hereafter E94) quasar SED with various contributions from galaxy emission and some dust reddening is remarkably effective in describing the SED shape from 0.3-3μm for large ranges of redshift, luminosity, black hole mass and Eddington ratio of type 1 AGN. In particular, the location in the mixing diagram of the highest luminosity AGN is very close (within 1σ) to that of the E94 SED. The mixing diagram can also be used to estimate the host galaxy fraction and reddening in quasar. We also show examples of some outliers which might be AGN in different evolutionary stages compared to the majority of AGN in the quasar-host galaxy co-evolution cycle.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Indications of a sub-linear and non-universal Kennicutt-Schmidt relationship

We estimate the parameters of the Kennicutt-Schmidt (KS) relationship, linking the star formation rate (Sigma_SFR) to the molecular gas surface density (Sigma_mol), in the STING sample of nearby disk galaxies using a hierarchical Bayesian method. This method rigorously treats measurement uncertainties, and provides accurate parameter estimates for both individual galaxies and the entire population. Assuming standard conversion factors to estimate Sigma_SFR and Sigma_mol from the observations, we find that the KS parameters vary between galaxies, indicating that no universal relationship holds for all galaxies. The KS slope of the whole population is 0.76, with the 2sigma range extending from 0.58 to 0.94. These results imply that the molecular gas depletion time is not constant, but varies from galaxy to galaxy, and increases with the molecular gas surface density. Therefore, other galactic properties besides just Sigma_mol affect Sigma_SFR, such as the gas fraction or stellar mass. The non-universality of the KS relationship indicates that a comprehensive theory of star formation must take into account additional physical processes that may vary from galaxy to galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Jet Emission in Young Radio Sources: a Fermi-LAT Gamma-ray View

We investigate the contribution of the beamed jet component to the high energy emission in young and compact extragalactic radio sources, focusing for the first time on the gamma-ray band. We derive predictions on the gamma-ray luminosities associated with the relativistic jet assuming a leptonic radiative model. The high energy emission is produced via Compton scattering by the relativistic electrons in a spherical region at the considered scales ($\lesssim$10 kpc). Simulations show a wide range of gamma-ray luminosities, with intensities up to $\sim10^{46}-10^{48}$ erg s$^{-1}$ depending on the assumed jet parameters. We find a highly linear relation between the simulated X-ray and gamma-ray luminosities that can be used to select candidates for a gamma-ray detection. We compare the simulated luminosity distributions in the radio, X-ray and gamma-ray regimes with observations for the largest sample of X-ray detected young radio quasars. Our analysis of $\sim$4-year Fermi Large Area Telescope (LAT) data does not give any statistically significant detection. However, the majority of the model-predicted gamma-ray fluxes for the sample are near or below the current Fermi-LAT flux threshold and compatible with the derived upper limits. Our study gives constraints on the minimum jet power ($L_{jet,kin}/L_{disk}>0.01$), on a potential jet contribution to the X-ray emission in the most compact sources ($\lesssim1$ kpc) and on the particles to magnetic field energy density ratio in broad agreement with equipartition assumption.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quasar Black Hole Mass Estimates in the Era of Time Domain Astronomy

We investigate the dependence of the normalization of the high-frequency part of the X-ray and optical power spectral densities (PSD) on black hole mass for a sample of 39 active galactic nuclei (AGN) with black hole masses estimated from reverberation mapping or dynamical modeling. We obtained new Swift observations of PG 1426+015, which has the largest estimated black hole mass of the AGN in our sample. We develop a novel statistical method to estimate the PSD from a lightcurve of photon counts with arbitrary sampling, eliminating the need to bin a lightcurve to achieve Gaussian statistics, and we use this technique to estimate the X-ray variability parameters for the faint AGN in our sample. We find that the normalization of the high-frequency X-ray PSD is inversely proportional to black hole mass. We discuss how to use this scaling relationship to obtain black hole mass estimates from the short time-scale X-ray variability amplitude with precision ~ 0.38 dex. The amplitude of optical variability on time scales of days is also anti-correlated with black hole mass, but with larger scatter. Instead, the optical variability amplitude exhibits the strongest anti-correlation with luminosity. We conclude with a discussion of the implications of our results for estimating black hole mass from the amplitude of AGN variability.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Description of Quasar Variability Measured Using Repeated SDSS and POSS Imaging

We provide a quantitative description and statistical interpretation of the optical continuum variability of quasars. The Sloan Digital Sky Survey (SDSS) has obtained repeated imaging in five UV-to-IR photometric bands for 33,881 spectroscopically confirmed quasars. About 10,000 quasars have an average of 60 observations in each band obtained over a decade along Stripe 82 (S82), whereas the remaining ~25,000 have 2-3 observations due to scan overlaps. The observed time lags span the range from a day to almost 10 years, and constrain quasar variability at rest-frame time lags of up to 4 years, and at rest-frame wavelengths from 1000A to 6000A. We publicly release a user-friendly catalog of quasars from the SDSS Data Release 7 that have been observed at least twice in SDSS or once in both SDSS and the Palomar Observatory Sky Survey, and we use it to analyze the ensemble properties of quasar variability. Based on a damped random walk (DRW) model defined by a characteristic time scale and an asymptotic variability amplitude that scale with the luminosity, black hole mass, and rest wavelength for individual quasars calibrated in S82, we can fully explain the ensemble variability statistics of the non-S82 quasars such as the exponential distribution of large magnitude changes. All available data are consistent with the DRW model as a viable description of the optical continuum variability of quasars on time scales of ~5-2000 days in the rest frame. We use these models to predict the incidence of quasar contamination in transient surveys such as those from PTF and LSST.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Disclosing the Radio Loudness Distribution Dichotomy in Quasars: An Unbiased Monte Carlo Approach Applied to the SDSS-FIRST Quasar Sample

We investigate the dichotomy in the radio loudness distribution of quasars by modelling their radio emission and various selection effects using a Monte Carlo approach. The existence of two physically distinct quasar populations, the radio-loud and radio-quiet quasars, is controversial and over the last decade a bimodal distribution of radio loudness of quasars has been both affirmed and disputed. We model the quasar radio luminosity distribution with simple unimodal and bimodal distribution functions. The resulting simulated samples are compared to a fiducial sample of 8,300 quasars drawn from the SDSS DR7 Quasar Catalog and combined with radio observations from the FIRST survey. Our results indicate that the SDSS-FIRST sample is best described by a radio loudness distribution which consists of two components, with 12+/-1 % of sources in the radio-loud component. On the other hand, the evidence for a local minimum in the loudness distribution (bimodality) is not strong and we find that previous claims for its existence were probably affected by the incompleteness of the FIRST survey close to its faint limit. We also investigate the redshift and luminosity dependence of the radio loudness distribution and find tentative evidence that at high redshift radio-loud quasars were rarer, on average "louder", and exhibited a smaller range in radio loudness. In agreement with other recent work, we conclude that the SDSS-FIRST sample strongly suggests that the radio loudness distribution of quasars is not a universal function, and that more complex models than presented here are needed to fully explain available observations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dust SEDs in the era of Herschel and Planck: a Hierarchical Bayesian fitting technique

We present a hierarchical Bayesian method for fitting infrared spectral energy distributions (SEDs) of dust emission to observed fluxes. Under the standard assumption of optically thin single temperature (T) sources the dust SED as represented by a power--law modified black body is subject to a strong degeneracy between T and the spectral index beta. The traditional non-hierarchical approaches, typically based on chi-square minimization, are severely limited by this degeneracy, as it produces an artificial anti-correlation between T and beta even with modest levels of observational noise. The hierarchical Bayesian method rigorously and self-consistently treats measurement uncertainties, including calibration and noise, resulting in more precise SED fits. As a result, the Bayesian fits do not produce any spurious anti-correlations between the SED parameters due to measurement uncertainty. We demonstrate that the Bayesian method is substantially more accurate than the chi-square fit in recovering the SED parameters, as well as the correlations between them. As an illustration, we apply our method to Herschel and sub millimeter ground-based observations of the star-forming Bok globule CB244. This source is a small, nearby molecular cloud containing a single low-mass protostar and a starless core. We find that T and beta are weakly positively correlated -- in contradiction with the chi-square fits, which indicate a T-beta anti-correlation from the same data-set. Additionally, in comparison to the chi-square fits the Bayesian SED parameter estimates exhibit a reduced range in values.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Evidence for a non-universal Kennicutt-Schmidt relationship using hierarchical Bayesian linear regression

For investigating the relationship between the star formation rate and gas surface density, we develop a Bayesian linear regression method that rigorously treats measurement uncertainties and accounts for hierarchical data structure. The hierarchical method simultaneously estimates the intercept, slope, and scatter about the regression line of each individual subject and the population. Using synthetic datasets, we demonstrate that the method recovers the underlying parameters of both the individuals and the population, especially when compared to commonly employed ordinary least squares techniques, such as the bisector fit. We apply the hierarchical method to estimate the Kennicutt-Schmidt (KS) parameters of a sample of spiral galaxies compiled by Bigiel et al. (2008). We find significant variation in the KS parameters, indicating that no single relationship holds for all galaxies. This suggests that the relationship between molecular gas and star formation differs between galaxies, possibly due to the influence of other properties within a given galaxy, such as metallicity, molecular gas fraction, or stellar mass. In four galaxies the slope estimates are sub-linear (at the 95% level). We estimate the mean index of the KS relationship for the population to be 0.84, with 95% range [0.63, 1.0]. For the galaxies with sub-linear KS relationships, a possible interpretation is that CO emission is tracing some molecular gas that is not directly associated with star formation. Equivalently, a sub-linear KS relationship may be indicative of an increasing gas depletion time at higher surface densities, as traced by CO emission. The Bayesian method can account for all sources of uncertainties, including variations in the conversion of observed intensities to star formation rates and gas surface densities, and is therefore well suited for a thorough statistical analysis of the KS law.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Recalibration of the virial factor and M-sigma relation for local active galaxies

Determining the virial factor of the broad-line region (BLR) gas is crucial for calibrating AGN black hole mass estimators, since the measured line-of-sight velocity needs to be converted into the intrinsic virial velocity. The average virial factor has been empirically calibrated based on the M-sigma relation of quiescent galaxies, but the claimed values differ by a factor of two in recent studies. We investigate the origin of the difference by measuring the M-sigma relation using an updated galaxy sample from the literature, and explore the dependence of the virial factor on various fitting methods. We find that the discrepancy is primarily caused by the sample selection, while the difference stemming from the various regression methods is marginal. However, we generally prefer the FITEXY and Bayesian estimators based on Monte Carlo simulations for the M-sigma relation. In addition, the choice of independent variable in the regression leads to ~0.2 dex variation in the virial factor inferred from the calibration process. Based on the determined virial factor, we present the updated M-sigma relation of local active galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Demographics of Broad Line Quasars in the Mass-Luminosity Plane II. Black Hole Mass and Eddington Ratio Functions

We employ a flexible Bayesian technique to estimate the black hole mass and Eddington ratio functions for Type 1 (i.e., broad line) quasars from a uniformly-selected data set of ~58,000 quasars from the SDSS DR7. We find that the SDSS becomes significantly incomplete at M_{BH} < 3 x 10^8 M_{Sun} or L / L_{Edd} < 0.07, and that the number densities of Type 1 quasars continue to increase down to these limits. Both the mass and Eddington ratio functions show evidence of downsizing, with the most massive and highest Eddington ratio black holes experiencing Type 1 quasar phases first, although the Eddington ratio number densities are flat at z < 2. We estimate the maximum Eddington ratio of Type 1 quasars in the observable Universe to be L / L_{Edd} ~ 3. Consistent with our results in Paper I, we do not find statistical evidence for a so-called "sub-Eddington boundary" in the mass-luminosity plane of broad line quasars, and demonstrate that such an apparent boundary in the observed distribution can be caused by selection effect and errors in virial BH mass estimates. Based on the typical Eddington ratio in a given mass bin, we estimate typical growth times for the black holes in Type 1 quasars and find that they are typically comparable to or longer than the age of the universe, implying an earlier phase of accelerated (i.e., with higher Eddington ratios) and possibly obscured growth. The large masses probed by our sample imply that most of our black holes reside in what are locally early type galaxies, and we interpret our results within the context of models of self-regulated black hole growth.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Demographics of Broad-Line Quasars in the Mass-Luminosity Plane. I. Testing FWHM-Based Virial Black Hole Masses

We jointly constrain the luminosity function (LF) and black hole mass function (BHMF) of broad-line quasars with forward Bayesian modeling in the quasar mass-luminosity plane, based on a homogeneous sample of ~ 58,000 SDSS DR7 quasars at z ~ 0.3-5. We take into account the selection effect of the sample flux limit; more importantly, we deal with the statistical scatter between true BH masses and FWHM-based single-epoch virial mass estimates, as well as potential luminosity-dependent biases of these mass estimates. The LF is tightly constrained in the regime sampled by SDSS, and makes reasonable predictions when extrapolated to ~ 3 magnitudes fainter. Downsizing is seen in the model LF. On the other hand, we find it difficult to constrain the BHMF to within a factor of a few at z>~0.7 (with MgII and CIV-based virial BH masses). This is mainly driven by the unknown luminosity-dependent bias of these mass estimators and its degeneracy with other model parameters, and secondly driven by the fact that SDSS quasars only sample the tip of the active BH population at high redshift. Nevertheless, the most likely models favor a positive luminosity-dependent bias for MgII and possibly for CIV, such that at fixed true BH mass, objects with higher-than-average luminosities have over-estimated FWHM-based virial masses. There is tentative evidence that downsizing also manifests itself in the active BHMF, and the BH mass density in broad-line quasars contributes an insignificant amount to the total BH mass density at all times. Within our model uncertainties, we do not find a strong BH mass dependence of the mean Eddington ratio; but there is evidence that the mean Eddington ratio (at fixed BH mass) increases with redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The linewidth-size relationship in the dense ISM of the Central Molecular Zone

The linewidth (sigma) - size (R) relationship has been extensively measured and analysed, in both the local ISM and in nearby normal galaxies. Generally, a power-law describes the relationship well with an index ranging from 0.2-0.6, now referred to as one of "Larson's Relationships." The nature of turbulence and star formation is considered to be intimately related to these relationships, so evaluating the sigma-R correlations in various environments is important for developing a comprehensive understanding of the ISM. We measure the sigma-R relationship in the Central Molecular Zone (CMZ) of the Galactic Centre using spectral line observations of the high density tracers N2H+, HCN, H13CN, and HCO+. We use dendrograms, which map the hierarchical nature of the position-position-velocity (PPV) data, to compute sigma and R of contiguous structures. The dispersions range from ~2-30 km/s in structures spanning sizes 2-40 pc, respectively. By performing Bayesian inference, we show that a power-law with exponent 0.3-1.1 can reasonably describe the sigma-R trend. We demonstrate that the derived sigma-R relationship is independent of the locations in the PPV dataset where sigma and R are measured. The uniformity in the sigma-R relationship suggests turbulence in the CMZ is driven on the large scales beyond >30 pc. We compare the CMZ sigma-R relationship to that measured in the Galactic molecular cloud Perseus. The exponents between the two systems are similar, suggestive of a connection between the turbulent properties within a cloud to its ambient medium. Yet, the velocity dispersion in the CMZ is systematically higher, resulting in a coefficient that is nearly five times larger. The systematic enhancement of turbulent velocities may be due to the combined effects of increased star formation activity, larger densities, and higher pressures relative to the local ISM.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Stochastic Model for the Luminosity Fluctuations of Accreting Black Holes

In this work we have developed a new stochastic model for the fluctuations in lightcurves of accreting black holes. The model is based on a linear combination of stochastic processes and is also the solution to the linear diffusion equation perturbed by a spatially correlated noise field. This allows flexible modeling of the power spectral density (PSD), and we derive the likelihood function for the process, enabling one to estimate the parameters of the process, including break frequencies in the PSD. Our statistical technique is computationally efficient, unbiased by aliasing and red noise leak, and fully accounts for irregular sampling and measurement errors. We show that our stochastic model provides a good approximation to the X-ray lightcurves of galactic black holes, and the optical and X-ray lightcurves of AGN. We use the estimated time scales of our stochastic model to recover the correlation between characteristic time scale of the high frequency X-ray fluctuations and black hole mass for AGN, including two new `detections' of the time scale for Fairall 9 and NGC 5548. We find a tight anti-correlation between the black hole mass and the amplitude of the driving noise field, which is proportional to the amplitude of the high frequency X-ray PSD, and we estimate that this parameter gives black hole mass estimates to within ~ 0.2 dex precision, potentially the most accurate method for AGN yet. We also find evidence that ~ 13% of AGN optical PSDs fall off flatter than 1 / f^2, and, similar to previous work, find that the optical fluctuations are more suppressed on short time scales compared to the X-rays, but are larger on long time scales, suggesting the optical fluctuations are not solely due to reprocessing of X-rays.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Accretion Rate and the Physical Nature of Unobscured Active Galaxies

We show how accretion rate governs the physical properties of a sample of unobscured broad-line, narrow-line, and lineless active galactic nuclei (AGNs). We avoid the systematic errors plaguing previous studies of AGN accretion rate by using accurate accretion luminosities (L_int) from well-sampled multiwavelength SEDs from the Cosmic Evolution Survey (COSMOS), and accurate black hole masses derived from virial scaling relations (for broad-line AGNs) or host-AGN relations (for narrow-line and lineless AGNs). In general, broad emission lines are present only at the highest accretion rates (L_int/L_Edd > 0.01), and these rapidly accreting AGNs are observed as broad-line AGNs or possibly as obscured narrow-line AGNs. Narrow-line and lineless AGNs at lower specific accretion rates (L_int/L_Edd < 0.01) are unobscured and yet lack a broad line region. The disappearance of the broad emission lines is caused by an expanding radiatively inefficient accretion flow (RIAF) at the inner radius of the accretion disk. The presence of the RIAF also drives L_int/L_Edd < 10^-2 narrow-line and lineless AGNs to 10 times higher ratios of radio to optical/UV emission than L_int/L_Edd > 0.01 broad-line AGNs, since the unbound nature of the RIAF means it is easier to form a radio outflow. The IR torus signature also tends to become weaker or disappear from L_int/L_Edd < 0.01 AGNs, although there may be additional mid-IR synchrotron emission associated with the RIAF. Together these results suggest that specific accretion rate is an important physical "axis" of AGN unification, described by a simple model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mass Functions of Supermassive Black Holes Across Cosmic Time

The black hole mass function of supermassive black holes describes the evolution of the distribution of black hole mass. It is one of the primary empirical tools available for mapping the growth of supermassive black holes and for constraining theoretical models of their evolution. In this review we discuss methods for estimating the black hole mass function, including their advantages and disadvantages. We also review the results of using these methods for estimating the mass function of both active and inactive black holes. In addition, we review current theoretical models for the growth of supermassive black holes that predict the black hole mass function. We conclude with a discussion of directions for future research which will lead to improvement in both empirical and theoretical determinations of the mass function of supermassive black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Measurement Error Models in Astronomy

I discuss the effects of measurement error on regression and density estimation. I review the statistical methods that have been developed to correct for measurement error that are most popular in astronomical data analysis, discussing their advantages and disadvantages. I describe functional models for accounting for measurement error in regression, with emphasis on the methods of moments approach and the modified loss function approach. I then describe structural models for accounting for measurement error in regression and density estimation, with emphasis on maximum-likelihood and Bayesian methods. As an example of a Bayesian application, I analyze an astronomical data set subject to large measurement errors and a non-linear dependence between the response and covariate. I conclude with some directions for future research.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Spectropolarimetric Evidence for Radiatively Inefficient Accretion in an Optically Dull Active Galaxy

We present Subaru/FOCAS spectropolarimetry of two active galaxies in the Cosmic Evolution Survey. These objects were selected to be optically dull, with the bright X-ray emission of an AGN but missing optical emission lines in our previous spectroscopy. Our new observations show that one target has very weak emission lines consistent with an optically dull AGN, while the other object has strong emission lines typical of a host-diluted Type 2 Seyfert galaxy. In neither source do we observe polarized emission lines, with 3-sigma upper limits of P_BLR < 2%. This means that the missing broad emission lines (and weaker narrow emission lines) are not due to simple anisotropic obscuration, e.g., by the canonical AGN torus. The weak-lined optically dull AGN exhibits a blue polarized continuum with P = 0.78 +/- 0.07% at 4400 A < lambda_rest < 7200 A (P = 1.37 +/- 0.16% at 4400 A < lambda_rest < 5050 A). The wavelength dependence of this polarized flux is similar to that of an unobscured AGN continuum and represents the intrinsic AGN emission, either as synchrotron emission or the outer part of an accretion disk reflected by a clumpy dust scatterer. Because this intrinsic AGN emission lacks emission lines, this source is likely to have a radiatively inefficient accretion flow.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Using the Fundamental Plane of Black Hole Activity to Distinguish X-ray Processes from Weakly Accreting Black Holes

The fundamental plane of black hole activity is a relation between X-ray luminosity, radio luminosity, and black hole mass for hard state Galactic black holes and their supermassive analogs. The fundamental plane suggests that, at low-accretion rates, the physical processes regulating the conversion of an accretion flow into radiative energy could be universal across the entire black hole mass scale. However, there is still a need to further refine the fundamental plane in order to better discern the radiative processes and their geometry very close to the black hole, in particular the source of hard X-rays. Further refinement is necessary because error bars on the best-fit slopes of the fundamental plane are generally large, and also the inferred coefficients can be sensitive to the adopted sample of black holes. In this work, we regress the fundamental plane with a Bayesian technique. Our approach shows that sub-Eddington black holes emit X-ray emission that is predominantly optically thin synchrotron radiation from the jet, provided that their radio spectra are flat or inverted. X-ray emission dominated by very radiatively inefficient accretion flows are excluded at the >3σ level. We also show that it is difficult to place FR I galaxies onto the fundamental plane because their X-ray jet emission is highly affected by synchrotron cooling. On the other hand, BL Lac objects fit onto the fundamental plane. Including a uniform subset of high-energy peaked BL Lac objects from the SDSS, we find sub-Eddington black holes with flat/inverted radio spectra follow log L_x=(1.45\pm0.04)log L_R-(0.88\pm0.06)\logM_{BH}-6.07\pm1.10, with σ_{int}=0.07\pm0.05 dex. Finally, we discuss how the effects of synchrotron cooling of jet emission from the highest black hole masses can bias fundamental plane regressions, perhaps leading to incorrect inferences on X-ray radiation mechanisms.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Constraints on Black Hole Growth, Quasar Lifetimes, and Eddington Ratio Distributions from the SDSS Broad Line Quasar Black Hole Mass Function

We present an estimate of the black hole mass function (BHMF) of broad line quasars (BLQSOs) that self-consistently corrects for incompleteness and the statistical uncertainty in the mass estimates, based on a sample of 9886 quasars at 1 < z < 4.5 drawn from the Sloan Digital Sky Survey. We find evidence for `cosmic downsizing' of black holes in BLQSOs, where the peak in their number density shifts to higher redshift with increasing black hole mass. The cosmic mass density for black holes seen as BLQSOs peaks at z ~ 2. We estimate the completeness of the SDSS as a function of black hole mass and Eddington ratio, and find that at z > 1 it is highly incomplete at M_BH < 10^9 M_Sun and L / L_Edd < 0.5. We also estimate a lower limit on the lifetime of a single BLQSO phase and we place constraints on the maximum mass of a black hole in a BLQSO. Our estimated distribution of BLQSO Eddington ratios peaks at L / L_Edd ~ 0.05 and has a dispersion of ~ 0.4 dex, implying that most BLQSOs are not radiating at or near the Eddington limit; however the location of the peak is subject to considerable uncertainty. The steep increase in number density of BLQSOs toward lower Eddington ratios is expected if the BLQSO accretion rate monotonically decays with time. Furthermore, our estimated lifetime and Eddington ratio distributions imply that the majority of the most massive black holes spend a significant amount of time growing in an earlier obscured phase, a conclusion which is independent of the unknown obscured fraction. These results are consistent with models for self-regulated black hole growth, at least for massive systems at z > 1, where the BLQSO phase occurs at the end of a fueling event when black hole feedback unbinds the accreting gas, halting the accretion flow.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Submillimeter Variability and the Gamma-ray Connection in Fermi Blazars

We present multi-epoch observations from the Submillimeter Array (SMA) for a sample of 171 bright blazars, 43 of which were detected by Fermi during the first three months of observations. We explore the correlation between their gamma-ray properties and submillimeter observations of their parsec-scale jets, with a special emphasis on spectral index in both bands and the variability of the synchrotron component. Subclass is determined using a combination of Fermi designation and the Candidate Gamma-Ray Blazar Survey (CGRaBS), resulting in 35 BL Lac objects and 136 flat-spectrum radio quasars (FSRQs) in our total sample. We calculate submillimeter energy spectral indices using contemporaneous observations in the 1 mm and 850 micron bands during the months August-October 2008. The submillimeter light curves are modeled as first-order continuous autoregressive processes, from which we derive characteristic timescales. Our blazar sample exhibits no differences in submillimeter variability amplitude or characteristic timescale as a function of subclass or luminosity. All of the the light curves are consistent with being produced by a single process that accounts for both low and high states, and there is additional evidence that objects may be transitioning between blazar class during flaring epochs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Impact of the Uncertainty in Single-Epoch Virial Black Hole Mass Estimates on the Observed Evolution of the Black Hole - Bulge Scaling Relations

Recent observations of the black hole (BH) - bulge scaling relations usually report positive redshift evolution, with higher redshift galaxies harboring more massive BHs than expected from the local relations. All of these studies focus on broad line quasars with BH mass estimated from virial estimators based on single-epoch spectra. Since the sample selection is largely based on quasar luminosity, the cosmic scatter in the BH-bulge relation introduces a statistical bias leading to on average more massive BHs given galaxy properties at high redshift (Lauer et al. 2007). We here emphasize a previously under-appreciated statistical bias resulting from the uncertainty of single-epoch virial BH mass estimators and the shape of the underlying (true) BH mass function, which leads to on average overestimation of the true BH masses at the high-mass end (Shen et al. 2008). We demonstrate that the latter virial mass bias can contribute a substantial amount to the observed excess in BH mass at fixed bulge properties, comparable to the Lauer et al. bias. The virial mass bias is independent of the Lauer et al. bias, hence if both biases are at work, they can largely (or even fully) account for the observed BH mass excess at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Are the Variations in Quasar Optical Flux Driven by Thermal Fluctuations?

We analyze a sample of optical light curves for 100 quasars, 70 of which have black hole mass estimates. Our sample is the largest and broadest used yet for modeling quasar variability. The sources in our sample have z < 2.8 and 10^6 < M_BH < 10^10. We model the light curves as a continuous time stochastic process, providing a natural means of estimating the characteristic time scale and amplitude of quasar variations. We employ a Bayesian approach to estimate the characteristic time scale and amplitude of flux variations; our approach is not affected by biases introduced from discrete sampling effects. We find that the characteristic time scales strongly correlate with black hole mass and luminosity, and are consistent with disk orbital or thermal time scales. In addition, the amplitude of short time scale variations is significantly anti-correlated with black hole mass and luminosity. We interpret the optical flux fluctuations as resulting from thermal fluctuations that are driven by an underlying stochastic process, such as a turbulent magnetic field. In addition, the intranight variations in optical flux implied by our empirical model are < 0.02 mag, consistent with current microvariability observations of radio-quiet quasars. Our stochastic model is therefore able to unify both long and short time scale optical variations in radio-quiet quasars as resulting from the same underlying process, while radio-loud quasars have an additional variability component that operates on time scales < 1 day.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Nature of Optically Dull Active Galactic Nuclei in COSMOS

We present infrared, optical, and X-ray data of 48 X-ray bright, optically dull AGNs in the COSMOS field. These objects exhibit the X-ray luminosity of an active galactic nucleus (AGN) but lack broad and narrow emission lines in their optical spectrum. We show that despite the lack of optical emission lines, most of these optically dull AGNs are not well-described by a typical passive red galaxy spectrum: instead they exhibit weak but significant blue emission like an unobscured AGN. Photometric observations over several years additionally show significant variability in the blue emission of four optically dull AGNs. The nature of the blue and infrared emission suggest that the optically inactive appearance of these AGNs cannot be caused by obscuration intrinsic to the AGNs. Instead, up to ~70% of optically dull AGNs are diluted by their hosts, with bright or simply edge-on hosts lying preferentially within the spectroscopic aperture. The remaining ~30% of optically dull AGNs have anomalously high f_x/f_o ratios and are intrinsically weak, not obscured, in the optical. These optically dull AGNs are best described as a weakly accreting AGN with a truncated accretion disk from a radiatively inefficient accretion flow.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Determining Quasar Black Hole Mass Functions from their Broad Emission Lines: Application to the Bright Quasar Survey

We describe a Bayesian approach to estimating quasar black hole mass functions (BHMF) when using the broad emission lines to estimate black hole mass. We show how using the broad line mass estimates in combination with statistical techniques developed for luminosity function estimation leads to statistically biased results. We derive the likelihood function for the BHMF based on the broad line mass estimates, and derive the posterior distribution for the BHMF, given the observed data. We develop our statistical approach for a flexible model where the BHMF is modelled as a mixture of Gaussian functions. Statistical inference is performed using markov chain monte carlo (MCMC) methods. Our method has the advantage that it is able to constrain the BHMF even beyond the survey detection limits at the adopted confidence level, accounts for measurement errors and the intrinsic uncertainty in broad line mass estimates, and provides a natural way of estimating the probability distribution of any quantities derived from the BHMF. We conclude by using our method to estimate the local active BHMF using the z < 0.5 Bright Quasar Survey sources. At z = 0.2, the quasar BHMF falls off approximately as a power law with slope ~ 2 for M_{BH} > 10^8. Our analysis implies that z < 0.5 broad line quasars have a typical Eddington ratio of ~ 0.4 and a dispersion in Eddington ratio of < 0.5 dex (abridged).

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The COSMOS AGN Spectroscopic Survey I: XMM Counterparts

We present optical spectroscopy for an X-ray and optical flux-limited sample of 677 XMM-Newton selected targets covering the 2 deg^2 COSMOS field, with a yield of 485 high-confidence redshifts. The majority of the spectra were obtained over three seasons (2005-2007) with the IMACS instrument on the Magellan (Baade) telescope. We also include in the sample previously published Sloan Digital Sky Survey spectra and supplemental observations with MMT/Hectospec. We detail the observations and classification analyses. The survey is 90% complete to flux limits of f_{0.5-10 keV}>8 x 10^-16 erg cm^-2 s^-1 and i_AB+<22, where over 90% of targets have high-confidence redshifts. Making simple corrections for incompleteness due to redshift and spectral type allows for a description of the complete population to $i_AB+<23. The corrected sample includes 57% broad emission line (Type 1, unobscured) AGN at 0.13<z<4.26, 25% narrow emission line (Type 2, obscured) AGN at 0.07<z<1.29, and 18% absorption line (host-dominated, obscured) AGN at 0<z<1.22 (excluding the stars that made up 4% of the X-ray targets). We show that the survey's limits in X-ray and optical flux include nearly all X-ray AGN (defined by L_{0.5-10 keV}>3 x 10^42 erg s^-1) to z<1, of both optically obscured and unobscured types. We find statistically significant evidence that the obscured to unobscured AGN ratio at z<1 increases with redshift and decreases with luminosity.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Atmospheric Parameters of Field L and T Dwarfs

We present an analysis of the 0.95-14.5 micron spectral energy distributions of nine field ultracool dwarfs with spectral types ranging from L1 to T4.5. Effective temperatures, gravities, and condensate cloud sedimentation efficiencies are derived by comparing the data to synthetic spectra computed from atmospheric models that self-consistently include the formation of condensate clouds. Derived effective temperatures decrease steadily through the L1 to T4.5 spectral types and we confirm that the effective temperatures of ultracool dwarfs at the L/T transition are nearly constant, decreasing by only ~200 K from spectral types L7.5 to T4.5. The two objects in our sample with very red J-Ks colors are best fitted with synthetic spectra that have thick clouds which hints at a possible correlation between the near-infrared colors of L dwarfs and the condensate cloud properties. The fits to the two T dwarfs in our sample (T2 and T4.5) also suggest that the clouds become thinner in this spectral class, in agreement with previous studies. Restricting the fits to narrower wavelength ranges (i.e., individual photometric bands) almost always yields excellent agreement between the data and models. Limitations in our knowledge of the opacities of key absorbers such as FeH, VO, and CH4 at certain wavelengths remain obvious, however. The effective temperatures obtained by fitting the narrower wavelength ranges can show a large scatter compared to the values derived by fitting the full spectral energy distributions; deviations are typically ~200 K and in the worst cases, up to 700 K.

0 citations0 reviewsNo code linkedOpen access