We explore in this paper the use of neural networks designed for point-clouds and sets on a new meta-learning task. We present experiments on the astronomical challenge of characterizing the stellar population of stellar streams. Stellar streams are elongated structures of stars in the outskirts of the Milky Way that form when a (small) galaxy breaks up under the Milky Way's gravitational force. We consider that we obtain, for each stream, a small 'support set' of stars that belongs to this stream. We aim to predict if the other stars in that region of the sky are from that stream or not, similar to one-class classification. Each "stream task" could also be transformed into a binary classification problem in a highly imbalanced regime (or supervised anomaly detection) by using the much bigger set of "other" stars and considering them as noisy negative examples. We propose to study the problem in the meta-learning regime: we expect that we can learn general information on characterizing a stream's stellar population by meta-learning across several streams in a fully supervised regime, and transfer it to new streams using only positive supervision. We present a novel use of Deep Sets, a model developed for point-cloud and sets, trained in a meta-learning fully supervised regime, and evaluated in a one-class classification setting. We compare it against Random Forests (with and without self-labeling) in the classic setting of binary classification, retrained for each task. We show that our method outperforms the Random-Forests even though the Deep Sets is not retrained on the new tasks, and accesses only a small part of the data compared to the Random Forest. We also show that the model performs well on a real-life stream when including additional fine-tuning.

Harnessing the sheer discovery potential of gravitational wave astronomy will require bold, deliberate, and sustained efforts to train and develop the requisite workforce. The next decade requires a strategic plan to build -- from the ground up -- a robust, open, and well-connected gravitational wave astronomy community with deep participation from traditional astronomers, physicists, data scientists, and instrumentalists. This basic infrastructure is sorely needed as an enabling foundation for research. We outline a set of recommendations for funding agencies, universities, and professional societies to help build a thriving, diverse, and inclusive new field.

Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems.

The first terrestrial gravitational wave interferometers have dramatically underscored the scientific value of observing the Universe through an entirely different window, and of folding this new channel of information with traditional astronomical data for a multimessenger view. The Laser Interferometer Space Antenna (LISA) will broaden the reach of gravitational wave astronomy by conducting the first survey of the millihertz gravitational wave sky, detecting tens of thousands of individual astrophysical sources ranging from white-dwarf binaries in our own galaxy to mergers of massive black holes at redshifts extending beyond the epoch of reionization. These observations will inform - and transform - our understanding of the end state of stellar evolution, massive black hole birth, and the co-evolution of galaxies and black holes through cosmic time. LISA also has the potential to detect gravitational wave emission from elusive astrophysical sources such as intermediate-mass black holes as well as exotic cosmological sources such as inflationary fields and cosmic string cusps.

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state.

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts.

Low-frequency gravitational-wave astronomy can perform precision tests of general relativity and probe fundamental physics in a regime previously inaccessible. A space-based detector will be a formidable tool to explore gravity's role in the cosmos, potentially telling us if and where Einstein's theory fails and providing clues about some of the greatest mysteries in physics and astronomy, such as dark matter and the origin of the Universe.

Binary black hole (BBH) mergers provide a prime source for current and future interferometric GW observatories. Massive BBH mergers may often take place in plasma-rich environments, leading to the exciting possibility of a concurrent electromagnetic (EM) signal observable by traditional astronomical facilities. However, many critical questions about the generation of such counterparts remain unanswered. We explore mechanisms that may drive EM counterparts with magnetohydrodynamic simulations treating a range of scenarios involving equal-mass black-hole binaries immersed in an initially homogeneous fluid with uniform, orbitally aligned magnetic fields. We find that the time development of Poynting luminosity, which may drive jet-like emissions, is relatively insensitive to aspects of the initial configuration. In particular, over a significant range of initial values, the central magnetic field strength is effectively regulated by the gas flow to yield a Poynting luminosity of $10^{45}-10^{46} \rho_{-13} M_8^2 \, {\rm erg}\,{\rm s}^{-1}$, with BBH mass scaled to $M_8 \equiv M/(10^8 M_{\odot})$ and ambient density $\rho_{-13} \equiv \rho/(10^{-13} \, {\rm g} \, {\rm cm}^{-3})$. We also calculate the direct plasma synchrotron emissions processed through geodesic ray-tracing. Despite lensing effects and dynamics, we find the observed synchrotron flux varies little leading up to merger.

During the Space Telescope and Optical Reverberation Mapping Project (STORM) observations of NGC 5548, the continuum and emission-line variability became de-correlated during the second half of the 6-month long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as a part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuum) in the soft X-ray band appears before the start of the anomalous emission-line behavior, peaks during the period of the anomaly, and then declines. This is a model-independent result suggesting that the soft excess is related to the anomaly. We divide the Swift data into on- and off-anomaly spectra to characterize the soft excess via spectral fitting. The cause of the spectral differences is likely due to a change in the intrinsic spectrum rather than being due to variable obscuration or partial covering. The Chandra spectra have lower signal-to-noise ratios, but are consistent with Swift data. Our preferred model of the soft excess is emission from an optically thick, warm Comptonizing corona, the effective optical depth of which increases during the anomaly. This model simultaneously explains all the three observations: the UV emission line flux decrease, the soft-excess increase, and the emission line anomaly.

We calculate the electromagnetic signal of a gamma-ray flare coming from the surface of a neutron star shortly before merger with a black hole companion. Using a new version of the Monte Carlo radiation transport code Pandurata that incorporates dynamic spacetimes, we integrate photon geodesics from the neutron star surface until they reach a distant observer or are captured by the black hole. The gamma-ray light curve is modulated by a number of relativistic effects, including Doppler beaming and gravitational lensing. Because the photons originate from the inspiraling neutron star, the light curve closely resembles the corresponding gravitational waveform: a chirp signal characterized by a steadily increasing frequency and amplitude. We propose to search for these electromagnetic chirps using matched filtering algorithms similar to those used in LIGO data analysis.

We present the results of an optical spectroscopic monitoring program targeting NGC 5548 as part of a larger multi-wavelength reverberation mapping campaign. The campaign spanned six months and achieved an almost daily cadence with observations from five ground-based telescopes. The H$\beta$ and He II $\lambda$4686 broad emission-line light curves lag that of the 5100 $\AA$ optical continuum by $4.17^{+0.36}_{-0.36}$ days and $0.79^{+0.35}_{-0.34}$ days, respectively. The H$\beta$ lag relative to the 1158 $\AA$ ultraviolet continuum light curve measured by the Hubble Space Telescope is roughly $\sim$50% longer than that measured against the optical continuum, and the lag difference is consistent with the observed lag between the optical and ultraviolet continua. This suggests that the characteristic radius of the broad-line region is $\sim$50% larger than the value inferred from optical data alone. We also measured velocity-resolved emission-line lags for H$\beta$ and found a complex velocity-lag structure with shorter lags in the line wings, indicative of a broad-line region dominated by Keplerian motion. The responses of both the H$\beta$ and He II $\lambda$4686 emission lines to the driving continuum changed significantly halfway through the campaign, a phenomenon also observed for C IV, Ly $\alpha$, He II(+O III]), and Si IV(+O IV]) during the same monitoring period. Finally, given the optical luminosity of NGC 5548 during our campaign, the measured H$\beta$ lag is a factor of five shorter than the expected value implied by the $R_\mathrm{BLR} - L_\mathrm{AGN}$ relation based on the past behavior of NGC 5548.

We conduct a multiwavelength continuum variability study of the Seyfert 1 galaxy NGC 5548 to investigate the temperature structure of its accretion disk. The 19 overlapping continuum light curves (1158 to 9157 angstroms) combine simultaneous HST , Swift , and ground-based observations over a 180 day period from 2014 January to July. Light-curve variability is interpreted as the reverberation response of the accretion disk to irradiation by a central time-varying point source. Our model yields the disk inclination, i, temperature T1 at 1 light day from the black hole, and a temperature-radius slope, alpha. We also infer the driving light curve and find that it correlates poorly with both the hard and soft X-ray light curves, suggesting that the X-rays alone may not drive the ultraviolet and optical variability over the observing period. We also decompose the light curves into bright, faint, and mean accretion-disk spectra. These spectra lie below that expected for a standard blackbody accretion disk accreting at L/LEdd = 0.1

During an intensive Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS) UV monitoring campaign of the Seyfert~1 galaxy NGC 5548 performed from 2014 February to July, the normally highly correlated far-UV continuum and broad emission-line variations decorrelated for ~60 to 70 days, starting ~75 days after the first HST/COS observation. Following this anomalous state, the flux and variability of the broad emission lines returned to a more normal state. This transient behavior, characterised by significant deficits in flux and equivalent width of the strong broad UV emission lines, is the first of its kind to be unambiguously identified in an active galactic nucleus reverberation mapping campaign. The largest corresponding emission-line flux deficits occurred for the high-ionization collisionally excited lines, C IV and Si IV(+O IV]), and also He II(+O III]), while the anomaly in Ly-alpha was substantially smaller. This pattern of behavior indicates a depletion in the flux of photons with E_{\rm ph} > 54 eV, relative to those near 13.6 eV. We suggest two plausible mechanisms for the observed behavior: (i) temporary obscuration of the ionizing continuum incident upon BLR clouds by a moving veil of material lying between the inner accretion disk and inner BLR, perhaps resulting from an episodic ejection of material from the disk, or (ii) a temporary change in the intrinsic ionizing continuum spectral energy distribution resulting in a deficit of ionizing photons with energies > 54 eV, possibly due to a transient restructuring of the Comptonizing atmosphere above the disk. Current evidence appears to favor the latter explanation.

We present ground-based optical photometric monitoring data for NGC 5548, part of an extended multi-wavelength reverberation mapping campaign. The light curves have nearly daily cadence from 2014 January to July in nine filters (\emph{BVRI} and \emph{ugriz}). Combined with ultraviolet data from the \emph{Hubble Space Telescope} and \emph{Swift}, we confirm significant time delays between the continuum bands as a function of wavelength, extending the wavelength coverage from 1158\,\AA\ to the $z$ band ($\sim\!9160$\,\AA). We find that the lags at wavelengths longer than the {\it V} band are equal to or greater than the lags of high-ionization-state emission lines (such as He\,{\sc ii}\,$\lambda 1640$ and $\lambda 4686$), suggesting that the continuum-emitting source is of a physical size comparable to the inner broad-line region (BLR). The trend of lag with wavelength is broadly consistent with the prediction for continuum reprocessing by an accretion disk with $\tau \propto \lambda^{4/3}$. However, the lags also imply a disk radius that is 3 times larger than the prediction from standard thin-disk theory, assuming that the bolometric luminosity is 10\% of the Eddington luminosity ($L = 0.1L_{\rm Edd}$). Using optical spectra from the Large Binocular Telescope, we estimate the bias of the interband continuum lags due to BLR emission observed in the filters. We find that the bias for filters with high levels of BLR contamination ($\sim\! 20\%$) can be important for the shortest continuum lags, and likely has a significant impact on the {\it u} and {\it U} bands owing to Balmer continuum emission.

We present the selection algorithm and anticipated results for the Time Domain Spectroscopic Survey (TDSS). TDSS is an SDSS-IV eBOSS subproject that will provide initial identification spectra of approximately 220,000 luminosity-variable objects (variable stars and AGN) across 7,500 square degrees selected from a combination of SDSS and multi-epoch Pan-STARRS1 photometry. TDSS will be the largest spectroscopic survey to explicitly target variable objects, avoiding pre-selection on the basis of colors or detailed modeling of specific variability characteristics. Kernel Density Estimate (KDE) analysis of our target population performed on SDSS Stripe 82 data suggests our target sample will be 95% pure (meaning 95% of objects we select have genuine luminosity variability of a few magnitudes or more). Our final spectroscopic sample will contain roughly 135,000 quasars and 85,000 stellar variables, approximately 4,000 of which will be RR Lyrae stars which may be used as outer Milky Way probes. The variability-selected quasar population has a smoother redshift distribution than a color-selected sample, and variability measurements similar to those we develop here may be used to make more uniform quasar samples in large surveys. The stellar variable targets are distributed fairly uniformly across color space, indicating that TDSS will obtain spectra for a wide variety of stellar variables including pulsating variables, stars with significant chromospheric activity, cataclysmic variables and eclipsing binaries. TDSS will serve as a pathfinder mission to identify and characterize the multitude of variable objects that will be detected photometrically in even larger variability surveys such as LSST.

We describe the first results from a six-month long reverberation-mapping experiment in the ultraviolet based on 170 observations of the Seyfert 1 galaxy NGC 5548 with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Significant correlated variability is found in the continuum and broad emission lines, with amplitudes ranging from ~30% to a factor of two in the emission lines and a factor of three in the continuum. The variations of all the strong emission lines lag behind those of the continuum, with He II 1640 lagging behind the continuum by ~2.5 days and Lyman alpha 1215, C IV 1550, and Si IV 1400 lagging by ~5-6 days. The relationship between the continuum and emission lines is complex. In particular, during the second half of the campaign, all emission-line lags increased by a factor of 1.3-2 and differences appear in the detailed structure of the continuum and emission-line light curves. Velocity-resolved cross-correlation analysis shows coherent structure in lag versus line-of-sight velocity for the emission lines; the high-velocity wings of C IV respond to continuum variations more rapidly than the line core, probably indicating higher velocity BLR clouds at smaller distances from the central engine. The velocity-dependent response of Lyman alpha, however, is more complex and will require further analysis.

Recent intensive Swift monitoring of the Seyfert 1 galaxy NGC 5548 yielded 282 usable epochs over 125 days across six UV/optical bands and the X-rays. This is the densest extended AGN UV/optical continuum sampling ever obtained, with a mean sampling rate <0.5 day. Approximately daily HST UV sampling was also obtained. The UV/optical light curves show strong correlations (r_max = 0.57 - 0.90) and the clearest measurement to date of interband lags. These lags are well-fit by a \tau propto \lambda^4/3 wavelength dependence, with a normalization that indicates an unexpectedly large disk radius of 0.35 +/- 0.05 lt-day at 1367 A, assuming a simple face-on model. The U-band shows a marginally larger lag than expected from the fit and surrounding bands, which could be due to Balmer continuum emission from the broad-line region as suggested by Korista and Goad. The UV/X-ray correlation is weaker (r_max < 0.45) and less consistent over time. This indicates that while Swift is beginning to measure UV/optical lags in general agreement with accretion disk theory (although the derived size is larger than predicted), the relationship with X-ray variability is less well understood. Combining this accretion disk size estimate with those from quasar microlensing studies suggests that AGN disk sizes scale approximately linearly with central black hole mass over a wide range of masses.

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.

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.

Robust gauge conditions are critically important to the stability and accuracy of numerical relativity (NR) simulations involving compact objects. Most of the NR community use the highly robust---though decade-old---moving-puncture (MP) gauge conditions for such simulations. It has been argued that in binary black hole (BBH) evolutions adopting this gauge, noise generated near adaptive-mesh-refinement (AMR) boundaries does not converge away cleanly with increasing resolution, severely limiting gravitational waveform accuracy at computationally feasible resolutions. We link this noise to a sharp (short-wavelength), initial outgoing gauge wave crossing into progressively lower resolution AMR grids, and present improvements to the standard MP gauge conditions that focus on stretching, smoothing, and more rapidly settling this outgoing wave. Our best gauge choice greatly reduces gravitational waveform noise during inspiral, yielding less fluctuation in convergence order and $\sim 40%$ lower waveform phase and amplitude errors at typical resolutions. Noise in other physical quantities of interest is also reduced, and constraint violations drop by more than an order of magnitude. We expect these improvements will carry over to simulations of all types of compact binary systems, as well as other $N$+1 formulations of gravity for which MP-like gauge conditions can be chosen.

We present time series analyses of the full Kepler dataset of Zw 229-15. This Kepler light curve --- with a baseline greater than three years, composed of virtually continuous, evenly sampled 30-minute measurements --- is unprecedented in its quality and precision. We utilize two methods of power spectral analysis to investigate the optical variability and search for evidence of a bend frequency associated with a characteristic optical variability timescale. Each method yields similar results. The first interpolates across data gaps to use the standard Fourier periodogram. The second, using the CARMA-based time-domain modeling technique of Kelly et al. (2014), does not need evenly-sampled data. Both methods find excess power at high frequencies that may be due to Kepler instrumental effects. More importantly both also show strong bends ({\Delta}{\alpha} ~ 2) at timescales of ~5 days, a feature similar to those seen in the X-ray PSDs of AGN but never before in the optical. This observed ~5 day timescale may be associated with one of several physical processes potentially responsible for the variability. A plausible association could be made with light-crossing, dynamical or thermal timescales, depending on the assumed value of the accretion disk size and on unobserved disk parameters such as {\alpha} and H/R. This timescale is not consistent with the viscous timescale, which would be years in a ~10^7 Solar mass AGN such as Zw 229-15. However there must be a second bend on long (>~1 year) timescales, and that feature could be associated with the viscous timescale.

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.

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.

[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 $\alpha = -1.87^{+0.26}_{-0.26}$, and a number density of $\log_{10} \phi^\star [\textrm{Mpc}^{-3}] = -3.24^{+0.25}_{-0.24}$. Integrated down to $M=-17.7$ this LF yields a luminosity density, $\log_{10} \epsilon [\textrm{erg}/\textrm{s/Hz/Mpc}^{3}] = 25.52^{+0.05}_{-0.05}$. Our LF analysis is consistent with previously published determinations within 1$\sigma$. 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$.

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.

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.

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.

We compile and analyse long term (~10 year) submillimetre (1.3, 0.87, 0.43 mm, submm) wavelength light curves of the Galactic centre black hole, Sagittarius A*. The 0.87 and 0.43 mm data are taken from the literature, while the majority of the 1.3 mm light curve is from previously unpublished SMA and CARMA data. We use Monte Carlo simulations to show that on minute to few hour time-scales the variability is consistent with a red noise process with a 230 GHz power spectrum slope of 2.3+0.8-0.6 at 95% confidence. The light curve is de-correlated (white noise) on very long (month to year) times. In order to identify the transition time between red and white noise, we model the light curves as a stochastic damped random walk process. The models allow a quantitative estimate of this physical characteristic time-scale of 8-4+3 hours at 230 GHz at 95% confidence, with consistent results at 345 and 690 GHz. This corresponds to ~10 orbital times or ~1 inflow (viscous) time at R = 3 Rs, a typical radius producing the 230 GHz emission as measured by very long baseline interferometry and found in theoretical accretion flow and jet models. This time-scale is significantly shorter (longer) than those measured at radio (near-infrared, NIR) wavelengths, and is marginally inconsistent with the same variability mechanism operating in the submm and NIR for the expected t ~ R^3/2 scaling. It is crudely consistent with the analogous time-scale inferred in studies of quasar optical light curves after accounting for the difference in emission radius. We find evidence that the submm variability persists at least down to the ISCO, if not the event horizon. These results can be compared quantitatively with similar analyses at different wavebands to test for connections between the variability mechanisms, and with light curves from theoretical models of accreting black holes.

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.

We define a quasar-galaxy mixing diagram using the slopes of their spectral energy distributions (SEDs) from 1\mu m to 3000\AA\ and from 1\mu m to 3\mu 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\mu 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\sigma) 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.

Striking similarities have been seen between accretion signatures of Galactic X-ray binary (XRB) systems and active galactic nuclei (AGN). XRB spectral states show a V-shaped correlation between X-ray spectral hardness and Eddington ratio as they vary, and some AGN samples reveal a similar trend, implying analogous processes at vastly larger masses and timescales. To further investigate the analogies, we have matched 617 sources from the Chandra Source Catalog to SDSS spectroscopy, and uniformly measured both X-ray and optical spectral characteristics across a broad range of AGN and galaxy types. We provide useful tabulations of X-ray spectral slope for broad and narrow line AGN, star-forming and passive galaxies and composite systems, also updating relationships between optical (Ha and [OIII]) line emission and X-ray luminosity. We further fit broadband spectral energy distributions with a variety of templates to estimate bolometric luminosity. Our results confirm a significant trend in AGN between X-ray spectral hardness and Eddington ratio expressed in X-ray luminosity, albeit with significant dispersion. The trend is not significant when expressed in the full bolometric or template-estimated AGN luminosity. We also confirm a relationship between the X-ray/optical spectral slope aox, and Eddington ratio, but it may not follow the trend predicted by analogy with XRB accretion states.

Over the past 20years reverberation mapping has proved one of the most successful techniques for studying the local (<1pc) environment of super-massive black holes that drive active galactic nuclei. Key successes of reverberation mapping have been direct black hole mass estimates, the radius-luminosity relation for the Hbeta line and the calibration of single-epoch mass estimators commonly employed up to z~7. However, observing constraints mean that few studies have been successful at z>0.1, or for the more-luminous quasars that make up the majority of current spectroscopic samples, or for the rest-frame ultra-violet emission lines available in optical spectra of z>0.5 objects. Previously we described a technique for stacking cross correlations to obtain reverberation mapping results at high z. Here we present the first results from a campaign designed for this purpose. We construct stacked cross-correlation functions for the CIV and MgII lines and find a clear peak in both. We find the peak in the MgII correlation is at longer lags than CIV consistent with previous results at low redshift. For the CIV sample we are able to bin by luminosity and find evidence for increasing lags for more-luminous objects. This CIV radius-luminosity relation is consistent with previous studies but with a fraction of the observational cost.

SCUBA-2 is an innovative 10000 pixel bolometer camera operating at submillimetre wavelengths on the James Clerk Maxwell Telescope (JCMT). The camera has the capability to carry out wide-field surveys to unprecedented depths, addressing key questions relating to the origins of galaxies, stars and planets. With two imaging arrays working simultaneously in the atmospheric windows at 450 and 850 microns, the vast increase in pixel count means that SCUBA-2 maps the sky 100-150 times faster than the previous SCUBA instrument. In this paper we present an overview of the instrument, discuss the physical characteristics of the superconducting detector arrays, outline the observing modes and data acquisition, and present the early performance figures on the telescope. We also showcase the capabilities of the instrument via some early examples of the science SCUBA-2 has already undertaken. In February 2012, SCUBA-2 began a series of unique legacy surveys for the JCMT community. These surveys will take 2.5 years and the results are already providing complementary data to the shorter wavelength, shallower, larger-area surveys from Herschel. The SCUBA-2 surveys will also provide a wealth of information for further study with new facilities such as ALMA, and future telescopes such as CCAT and SPICA.

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.

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.

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.

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.

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.

Reverberation mapping offers one of the best techniques for studying the inner regions of QSOs. It is based on cross-correlating continuum and emission-line light curves. New time-resolved optical surveys will produce well sampled light curves for many thousands of QSOs. We explore the potential of stacking samples to produce composite cross-correlations for groups of objects that have well sampled continuum light curves, but only a few (~2) emission-line measurements. This technique exploits current and future wide-field optical monitoring surveys (e.g. Pan-STARRS, LSST) and the multiplexing capability of multi-object spectrographs (e.g. 2dF, Hectospec) to significantly reduce the observational expense of reverberation mapping, in particular at high redshift (0.5 to 2.5). We demonstrate the technique using simulated QSO light curves and explore the biases involved when stacking cross-correlations in some simplified situations. We show that stacked cross correlations have smaller amplitude peaks compared to well sampled correlation functions as the mean flux of the emission light curve is poorly constrained. However, the position of the peak remains intact. We find there can be `kinks' in stacked correlation functions due to different measurements contributing to different parts of the correlation function. Using the Pan-STARRS Medium-Deep Survey (MDS) as a template we show that cross-correlation lags should be measurable in a sample size of 500 QSOs that have weekly photometric monitoring and two spectroscopic observations. Finally we apply the technique to a small sample (42) of QSOs that have light curves from the MDS. We find no indication of a peak in the stacked cross-correlation. A larger spectroscopic sample is required to produce robust reverberation lags.

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.

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.

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.

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.

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.

We conduct a descriptive analysis of the multipolar structure of gravitational-radiation waveforms from equal-mass aligned-spin mergers, following an approach first presented in the complementary context of nonspinning black holes of varying mass ratio [J.G. Baker et al. Phys. Rev. D 78 044046 (2008)]. We find that, as with the nonspinning mergers, the dominant waveform mode phases evolve together in lock-step through inspiral and merger, supporting the previous waveform description in terms of an adiabatically rigid rotator driving gravitational-wave emission-an implicit rotating source. We further apply the late-time merger-ringdown model for the rotational frequency introduced in [J.G. Baker et al. Phys. Rev. D 78 044046 (2008)], along with an improved amplitude model appropriate for the dominant (2, \pm2) modes. This provides a quantitative description of the merger-ringdown waveforms, and suggests that the major features of these waveforms can be described with reference only to the intrinsic parameters associated with the state of the final black hole formed in the merger. We provide an explicit model for the merger-ringdown radiation, and demonstrate that this model agrees to fitting factors better than 95% with the original numerical waveforms for system masses above \sim150M\odot. This model may be directly applicable to gravitational-wave detection of intermediate-mass black-hole mergers.

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\sigma\ 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 \sigma_{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.

We have produced the next generation of quasar spectral energy distributions (SEDs), essentially updating the work of Elvis et al. (1994) by using high-quality data obtained with several space and ground-based telescopes, including NASA's Great Observatories. We present an atlas of SEDs of 85 optically bright, non-blazar quasars over the electromagnetic spectrum from radio to X-rays. The heterogeneous sample includes 27 radio-quiet and 58 radio-loud quasars. Most objects have quasi-simultaneous ultraviolet-optical spectroscopic data, supplemented with some far-ultraviolet spectra, and more than half also have Spitzer mid-infrared IRS spectra. The X-ray spectral parameters are collected from the literature where available. The radio, far-infrared, and near-infrared photometric data are also obtained from either the literature or new observations. We construct composite spectral energy distributions for radio-loud and radio-quiet objects and compare these to those of Elvis et al., finding that ours have similar overall shapes, but our improved spectral resolution reveals more detailed features, especially in the mid and near-infrared.

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.

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.

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.

Recent breakthroughs in the field of numerical relativity have led to dramatic progress in understanding the predictions of General Relativity for the dynamical interactions of two black holes in the regime of very strong gravitational fields. Such black-hole binaries are important astrophysical systems and are a key target of current and developing gravitational-wave detectors. The waveform signature of strong gravitational radiation emitted as the black holes fall together and merge provides a clear observable record of the process. After decades of slow progress, these mergers and the gravitational-wave signals they generate can now be routinely calculated using the methods of numerical relativity. We review recent advances in understanding the predicted physics of events and the consequent radiation, and discuss some of the impacts this new knowledge is having in various areas of astrophysics.

Understanding the predictions of general relativity for the dynamical interactions of two black holes has been a long-standing unsolved problem in theoretical physics. Black-hole mergers are monumental astrophysical events, releasing tremendous amounts of energy in the form of gravitational radiation, and are key sources for both ground- and space-based gravitational-wave detectors. The black-hole merger dynamics and the resulting gravitational waveforms can only be calculated through numerical simulations of Einstein's equations of general relativity. For many years, numerical relativists attempting to model these mergers encountered a host of problems, causing their codes to crash after just a fraction of a binary orbit could be simulated. Recently, however, a series of dramatic advances in numerical relativity has allowed stable, robust black-hole merger simulations. This remarkable progress in the rapidly maturing field of numerical relativity, and the new understanding of black-hole binary dynamics that is emerging is chronicled. Important applications of these fundamental physics results to astrophysics, to gravitational-wave astronomy, and in other areas are also discussed.

We model the time variability of ~9,000 spectroscopically confirmed quasars in SDSS Stripe 82 as a damped random walk. Using 2.7 million photometric measurements collected over 10 years, we confirm the results of Kelly et al. (2009) and Koz{\l}owski et al. (2010) that this model can explain quasar light curves at an impressive fidelity level (0.01-0.02 mag). The damped random walk model provides a simple, fast [O(N) for N data points], and powerful statistical description of quasar light curves by a characteristic time scale (tau) and an asymptotic rms variability on long time scales (SF_inf). We searched for correlations between these two variability parameters and physical parameters such as luminosity and black hole mass, and rest-frame wavelength. We find that tau increases with increasing wavelength with a power law index of 0.17, remains nearly constant with redshift and luminosity, and increases with increasing black hole mass with power law index of 0.21+/-0.07. The amplitude of variability is anti-correlated with the Eddington ratio, which suggests a scenario where optical fluctuations are tied to variations in the accretion rate. The radio-loudest quasars have systematically larger variability amplitudes by about 30%, when corrected for the other observed trends, while the distribution of their characteristic time scale is indistinguishable from that of the full sample. We do not detect any statistically robust differences in the characteristic time scale and variability amplitude between the full sample and the small subsample of quasars detected by ROSAT. Our results provide a simple quantitative framework for generating mock quasar light curves, such as currently used in LSST image simulations. (abridged)

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.

We investigate the precision with which the parameters describing the characteristics and location of nonspinning black hole binaries can be measured with the Laser Interferometer Space Antenna (LISA). By using complete waveforms including the inspiral, merger and ringdown portions of the signals, we find that LISA will have far greater precision than previous estimates for nonspinning mergers that ignored the merger and ringdown. Our analysis covers nonspinning waveforms with moderate mass ratios, q >= 1/10, and total masses 10^5 < M/M_{Sun} < 10^7. We compare the parameter uncertainties using the Fisher matrix formalism, and establish the significance of mass asymmetry and higher-order content to the predicted parameter uncertainties resulting from inclusion of the merger. In real-time observations, the later parts of the signal lead to significant improvements in sky-position precision in the last hours and even the final minutes of observation. For comparable mass systems with total mass M/M_{Sun} = ~10^6, we find that the increased precision resulting from including the merger is comparable to the increase in signal-to-noise ratio. For the most precise systems under investigation, half can be localized to within O(10 arcmin), and 10% can be localized to within O(1 arcmin).

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.

Although the gravitational wave kick velocity in the orbital plane of coalescing black holes has been understood for some time, apparently conflicting formulae have been proposed for the dominant out-of-plane kick, each a good fit to different data sets. This is important to resolve because it is only the out-of-plane kicks that can reach more than 500 km/s and can thus eject merged remnants from galaxies. Using a different ansatz for the out-of-plane kick, we show that we can fit almost all existing data to better than 5 %. This is good enough for any astrophysical calculation, and shows that the previous apparent conflict was only because the two data sets explored different aspects of the kick parameter space.

We extend a recently developed galaxy morphology classification method, Quantitative Multiwavelength Morphology (QMM), to connect galaxy morphologies to their underlying physical properties. The traditional classification of galaxies approaches the problem separately through either morphological classification or, in more recent times, through analysis of physical properties. A combined approach has significant potential in producing a consistent and accurate classification scheme as well as shedding light on the origin and evolution of galaxy morphology. Here we present an analysis of a volume limited sample of 31703 galaxies from the fourth data release of the Sloan Digital Sky Survey. We use an image analysis method called Pixel-z to extract the underlying physical properties of the galaxies, which is then quantified using the concentration, asymmetry and clumpiness (CAS) parameters. The galaxies also have their multiwavelength morphologies quantified using QMM, and these results are then related to the distributed physical properties through a regression analysis. We show that this method can be used to relate the spatial distribution of physical properties with the morphological properties of galaxies.

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.

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.

Coalescing massive black hole binaries are produced by the mergers of galaxies. The final stages of the black hole coalescence produce strong gravitational radiation that can be detected by the space-borne LISA. In cases where the black hole merger takes place in the presence of gas and magnetic fields, various types of electromagnetic signals may also be produced. Modeling such electromagnetic counterparts of the final merger requires evolving the behavior of both gas and fields in the strong-field regions around the black holes. We have taken a step towards solving this problem by mapping the flow of pressureless matter in the dynamic, 3-D general relativistic spacetime around the merging black holes. We find qualitative differences in collision and outflow speeds, including a signature of the merger when the net angular momentum of the matter is low, between the results from single and binary black holes, and between nonrotating and rotating holes in binaries. If future magnetohydrodynamic results confirm these differences, it may allow assessment of the properties of the binaries as well as yielding an identifiable electromagnetic counterpart to the attendant gravitational wave signal.

We present black hole masses and accretion rates for 182 Type 1 AGN in COSMOS. We estimate masses using the scaling relations for the broad Hb, MgII, and CIV emission lines in the redshift ranges 0.16<z<0.88, 1<z<2.4, and 2.7<z<4.9. We estimate the accretion rate using an Eddington ratio L_I/L_Edd estimated from optical and X-ray data. We find that very few Type 1 AGN accrete below L_I/L_Edd ~ 0.01, despite simulations of synthetic spectra which show that the survey is sensitive to such Type 1 AGN. At lower accretion rates the BLR may become obscured, diluted or nonexistent. We find evidence that Type 1 AGN at higher accretion rates have higher optical luminosities, as more of their emission comes from the cool (optical) accretion disk with respect to shorter wavelengths. We measure a larger range in accretion rate than previous works, suggesting that COSMOS is more efficient at finding low accretion rate Type 1 AGN. However the measured range in accretion rate is still comparable to the intrinsic scatter from the scaling relations, suggesting that Type 1 AGN accrete at a narrow range of Eddington ratio, with L_I/L_Edd ~ 0.1.

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.

We present 24 micron photometry of the intermediate-age open cluster Praesepe. We assemble a catalog of 193 probable cluster members that are detected in optical databases, the Two Micron All Sky Survey (2MASS), and at 24 micron, within an area of ~ 2.47 square degrees. Mid-IR excesses indicating debris disks are found for one early-type and for three solar-type stars. Corrections for sampling statistics yield a 24 micron excess fraction (debris disk fraction) of 6.5 +- 4.1% for luminous and 1.9 +- 1.2% for solar-type stars. The incidence of excesses is in agreement with the decay trend of debris disks as a function of age observed for other cluster and field stars. The values also agree with those for older stars, indicating that debris generation in the zones that emit at 24 micron falls to the older 1-10 Gyr field star sample value by roughly 750 Myr. We discuss our results in the context of previous observations of excess fractions for early- and solar-type stars. We show that solar-type stars lose their debris disk 24 micron excesses on a shorter timescale than early-type stars. Simplistic Monte Carlo models suggest that, during the first Gyr of their evolution, up to 15-30% of solar-type stars might undergo an orbital realignment of giant planets such as the one thought to have led to the Late Heavy Bombardment, if the length of the bombardment episode is similar to the one thought to have happened in our Solar System. In the Appendix, we determine the cluster's parameters via boostrap Monte Carlo isochrone fitting, yielding an age of 757 Myr (+- 36 Myr at 1 sigma confidence) and a distance of 179 pc (+- 2 pc at 1 sigma confidence), not allowing for systematic errors.

Recent advances in numerical relativity provide a detailed description of the waveforms of coalescing massive black hole binaries (MBHBs), expected to be the strongest detectable LISA sources. We present a preliminary study of LISA's sensitivity to MBHB parameters using a hybrid numerical/analytic waveform for equal-mass, non-spinning holes. The Synthetic LISA software package is used to simulate the instrument response and the Fisher information matrix method is used to estimate errors in the parameters. Initial results indicate that inclusion of the merger signal can significantly improve the precision of some parameter estimates. For example, the median parameter errors for an ensemble of systems with total redshifted mass of one million Solar masses at a redshift of one were found to decrease by a factor of slightly more than two for signals with merger as compared to signals truncated at the Schwarzchild ISCO.

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.

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

Recent numerical relativistic results demonstrate that the merger of comparable-mass spinning black holes has a maximum ``recoil kick'' of up to $\sim 4000 \kms$. However the scaling of these recoil velocities with mass ratio is poorly understood. We present new runs showing that the maximum possible kick perpendicular to the orbital plane does not scale as $\sim\eta^2$ (where $\eta$ is the symmetric mass ratio), as previously proposed, but is more consistent with $\sim\eta^3$, at least for systems with low orbital precession. We discuss the effect of this dependence on galactic ejection scenarios and retention of intermediate-mass black holes in globular clusters.

We describe a Bayesian approach to estimating luminosity functions. We derive the likelihood function and posterior probability distribution for the luminosity function, given the observed data, and we compare the Bayesian approach with maximum-likelihood by simulating sources from a Schechter function. For our simulations confidence intervals derived from bootstrapping the maximum-likelihood estimate can be too narrow, while confidence intervals derived from the Bayesian approach are valid. We develop our statistical approach for a flexible model where the luminosity function is modeled as a mixture of Gaussian functions. Statistical inference is performed using Markov chain Monte Carlo (MCMC) methods, and we describe a Metropolis-Hastings algorithm to perform the MCMC. The MCMC simulates random draws from the probability distribution of the luminosity function parameters, given the data, and we use a simulated data set to show how these random draws may be used to estimate the probability distribution for the luminosity function. In addition, we show how the MCMC output may be used to estimate the probability distribution of any quantities derived from the luminosity function, such as the peak in the space density of quasars. The Bayesian method we develop has the advantage that it is able to place accurate constraints on the luminosity function even beyond the survey detection limits, and that it provides a natural way of estimating the probability distribution of any quantities derived from the luminosity function, including those that rely on information beyond the survey detection limits.

We assess the impact of starburst and AGN feedback-driven winds on the CO emission from galaxy mergers, and, in particular, search for signatures of these winds in the simulated CO morphologies and emission line profiles. We do so by combining a 3D non-LTE molecular line radiative transfer code with smoothed particle hydrodynamics (SPH) simulations of galaxy mergers that include prescriptions for star formation, black hole growth, a multiphase interstellar medium (ISM), and the winds associated with star formation and black hole growth. Our main results are: (1) Galactic winds can drive outflows of masses ~10^8-10^9 Msun which may be imaged via CO emission line mapping. (2) AGN feedback-driven winds are able to drive imageable CO outflows for longer periods of time than starburst-driven winds owing to the greater amount of energy imparted to the ISM by AGN feedback compared to star formation. (3) Galactic winds can control the spatial extent of the CO emission in post-merger galaxies, and may serve as a physical motivation for the sub-kiloparsec scale CO emission radii observed in local advanced mergers. (4) Secondary emission peaks at velocities greater than the circular velocity are seen in the CO emission lines in all models. In models with winds, these high velocity peaks are seen to preferentially correspond to outflowing gas entrained in winds, which is not the case in the model without winds. The high velocity peaks seen in models without winds are typically confined to velocity offsets (from the systemic) < 1.7 times the circular velocity, whereas the models with AGN feedback-driven winds can drive high velocity peaks to ~2.5 times the circular velocity.

In this work we use a sample of 318 radio-quiet quasars (RQQ) to investigate the dependence of the ratio of optical/UV flux to X-ray flux, alpha_ox, and the X-ray photon index, Gamma_X, on black hole mass, UV luminosity relative to Eddington, and X-ray luminosity relative to Eddington. Our sample is drawn from the SDSS, with X-ray data from ROSAT and Chandra, and optical data mostly from the SDSS; 153 of these sources have estimates of Gamma_X from Chandra. We estimate M_BH using standard estimates derived from the Hbeta, Mg II, and C IV broad emission lines. Our sample spans a broad range in black hole mass (10^6 < M_BH / M_Sun < 10^10) and redshift (z < 4.8). We find that alpha_ox increases with increasing M_BH and L_UV / L_Edd, and decreases with increasing L_X / L_Edd. In addition, we confirm the correlation seen in previous studies between Gamma_X and M_BH and both L_UV / L_Edd and L_X / L_Edd; however, we also find evidence that the dependence of Gamma_X of these quantities is not monotonic, changing sign at M_BH ~ 3 x 10^8 M_Sun. We argue that the alpha_ox correlations imply that the fraction of bolometric luminosity emitted by the accretion disk, as compared to the corona, increases with increasing accretion rate relative to Eddington. In addition, we argue that the Gamma_X trends are caused by a dependence of X-ray spectral index on accretion rate. We discuss our results within the context of accretion models with comptonizing corona, and discuss the implications of the alpha_ox correlations for quasar feedback. To date, this is the largest study of the dependence of RQQ X-ray parameters on black hole mass and related quantities, and the first to attempt to correct for the large statistical uncertainty in the broad line mass estimates.

We present a multipolar analysis of the gravitational recoil computed in recent numerical simulations of binary black hole (BH) coalescence, for both unequal masses and non-zero, non-precessing spins. We show that multipole moments up to and including l=4 are sufficient to accurately reproduce the final recoil velocity (within ~2%) and that only a few dominant modes contribute significantly to it (within ~5%). We describe how the relative amplitudes, and more importantly, the relative phases, of these few modes control the way in which the recoil builds up throughout the inspiral, merger, and ringdown phases. We also find that the numerical results can be reproduced by an ``effective Newtonian'' formula for the multipole moments obtained by replacing the radial separation in the Newtonian formulae with an effective radius computed from the numerical data. Beyond the merger, the numerical results are reproduced by a superposition of three Kerr quasi-normal modes (QNMs). Analytic formulae, obtained by expressing the multipole moments in terms of the fundamental QNMs of a Kerr BH, are able to explain the onset and amount of ``anti-kick'' for each of the simulations. Lastly, we apply this multipolar analysis to help explain the remarkable difference between the amplitudes of planar and non-planar kicks for equal-mass spinning black holes.

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.

Several groups have recently computed the gravitational radiation recoil produced by the merger of two spinning black holes. The results suggest that spin can be the dominant contributor to the kick, with reported recoil speeds of hundreds to even thousands of kilometers per second. The parameter space of spin kicks is large, however, and it is ultimately desirable to have a simple formula that gives the approximate magnitude of the kick given a mass ratio, spin magnitudes, and spin orientations. As a step toward this goal, we perform a systematic study of the recoil speeds from mergers of black holes with mass ratio $q\equiv m_1/m_2=2/3$ and dimensionless spin parameters of $a_1/m_1$ and $a_2/m_2$ equal to 0 or 0.2, with directions aligned or anti-aligned with the orbital angular momentum. We also run an equal-mass $a_1/m_1=-a_2/m_2=0.2$ case, and find good agreement with previous results. We find that, for currently reported kicks from aligned or anti-aligned spins, a simple kick formula inspired by post-Newtonian analyses can reproduce the numerical results to better than $\sim$10%.

Coalescing black-hole binaries are expected to be the strongest sources of gravitational waves for ground-based interferometers as well as the space-based interferometer LISA. Recent progress in numerical relativity now makes it possible to calculate the waveforms from the strong-field dynamical merger and is revolutionizing our understanding of these systems. We review these dramatic developments, emphasizing applications to issues in gravitational wave observations. These new capabilities also make possible accurate calculations of the recoil or kick imparted to the final remnant black hole when the merging components have unequal masses, or unequal or unaligned spins. We highlight recent work in this area, focusing on results of interest to astrophysics.

We present Gemini near-infrared spectroscopic observations of six luminous quasars at z=5.8$\sim$6.3. Five of them were observed using Gemini-South/GNIRS, which provides a simultaneous wavelength coverage of 0.9--2.5 $\mu$m in cross dispersion mode. The other source was observed in K band with Gemini-North/NIRI. We calculate line strengths for all detected emission lines and use their ratios to estimate gas metallicity in the broad-line regions of the quasars. The metallicity is found to be supersolar with a typical value of $\sim$4 Z_{\sun}, and a comparison with low-redshift observations shows no strong evolution in metallicity up to z$\sim$6. The FeII/MgII ratio of the quasars is 4.9+/-1.4, consistent with low-redshift measurements. We estimate central BH masses of 10^9 to 10^{10} M_{\sun} and Eddington luminosity ratios of order unity. We identify two MgII $\lambda\lambda$2796,2803 absorbers with rest equivalent width W_0^{\lambda2796}>1 \AA at 2.2<z<3 and three MgII absorbers with W_0^{\lambda2796}>1.5 \AA at z>3 in the spectra, with the two most distant absorbers at z=4.8668 and 4.8823, respectively. The redshift number densities (dN/dz) of MgII absorbers with W_0^{\lambda2796}>1.5 \AA are consistent with no cosmic evolution up to z>4.

I describe a Bayesian method to account for measurement errors in linear regression of astronomical data. The method allows for heteroscedastic and possibly correlated measurement errors, and intrinsic scatter in the regression relationship. The method is based on deriving a likelihood function for the measured data, and I focus on the case when the intrinsic distribution of the independent variables can be approximated using a mixture of Gaussians. I generalize the method to incorporate multiple independent variables, non-detections, and selection effects (e.g., Malmquist bias). A Gibbs sampler is described for simulating random draws from the probability distribution of the parameters, given the observed data. I use simulation to compare the method with other common estimators. The simulations illustrate that the Gaussian mixture model outperforms other common estimators and can effectively give constraints on the regression parameters, even when the measurement errors dominate the observed scatter, source detection fraction is low, or the intrinsic distribution of the independent variables is not a mixture of Gaussians. I conclude by using this method to fit the X-ray spectral slope as a function of Eddington ratio using a sample of 39 z < 0.8 radio-quiet quasars. I confirm the correlation seen by other authors between the radio-quiet quasar X-ray spectral slope and the Eddington ratio, where the X-ray spectral slope softens as the Eddington ratio increases.

We present Spitzer observations of the dark globule CB190 (L771). We observe a roughly circular 24 micron shadow with a 70 arcsec radius. The extinction profile of this shadow matches the profile derived from 2MASS photometry at the outer edges of the globule and reaches a maximum of ~32 visual magnitudes at the center. The corresponding mass of CB190 is ~10 Msun. Our 12CO and 13CO J = 2-1 data over a 10 arcmin X 10 arcmin region centered on the shadow show a temperature ~10 K. The thermal continuum indicates a similar temperature for the dust. The molecular data also show evidence of freezeout onto dust grains. We estimate a distance to CB190 of 400 pc using the spectroscopic parallax of a star associated with the globule. Bonnor-Ebert fits to the density profile, in conjunction with this distance, yield xi_max = 7.2, indicating that CB190 may be unstable. The high temperature (56 K) of the best fit Bonnor-Ebert model is in contradiction with the CO and thermal continuum data, leading to the conclusion that the thermal pressure is not enough to prevent free-fall collapse. We also find that the turbulence in the cloud is inadequate to support it. However, the cloud may be supported by the magnetic field, if this field is at the average level for dark globules. Since the magnetic field will eventually leak out through ambipolar diffusion, it is likely that CB190 is collapsing or in a late pre-collapse stage.

Using a sample of optically-selected quasars from the Sloan Digital Sky Survey, we have determined the radio-loud fraction (RLF) of quasars as a function of redshift and optical luminosity. The sample contains more than 30,000 objects and spans a redshift range of 0<z<5 and a luminosity range of -30<M_i<-22. We use both the radio-to-optical flux ratio (R parameter) and radio luminosity to define radio-loud quasars. After breaking the correlation between redshift and luminosity due to the flux-limited nature of the sample, we find that the RLF of quasars decreases with increasing redshift and decreasing luminosity. The relation can be described in the form of log(RLF/(1-RLF)) = b_0 + b_z log(1+z) + b_M (M_{2500}+26), where M_{2500} is the absolute magnitude at rest-frame 2500A, and b_z, b_M<0. When using R>10 to define radio-loud quasars, we find that b_0=-0.132+/-0.116, b_z=-2.052+/-0.261, and b_M=-0.183+/-0.025. The RLF at z=0.5 declines from 24.3% to 5.6% as luminosity decreases from M_{2500}=-26 to M_{2500}=-22, and the RLF at M_{2500}=-26 declines from 24.3% to 4.1% as redshift increases from 0.5 to 3, suggesting that the RLF is a strong function of both redshift and luminosity. We also examine the impact of flux-related selection effects on the RLF determination using a series of tests, and find that the dependence of the RLF on redshift and luminosity is highly likely to be physical, and the selection effects we considered are not responsible for the dependence.

We report new Chandra observations of seven optically faint, z \sim 4 radio-quiet quasars. We have combined these new observations with previous Chandra observations of radio-quiet quasars to create a sample of 174 sources. These sources have 0.1 < z < 4.7, and 10^{44} ergs s^{-1} < nu L_{nu} (2500 \AA) < 10^{48} ergs s^{-1}. The X-ray detection fraction is 90%. We find that the X-ray loudness of radio-quiet quasars decreases with UV luminosity and increases with redshift. The model that is best supported by the data has a linear dependence of optical-to-X-ray ratio, alpha_{ox}, on cosmic time, and a quadratic dependence of alpha_{ox} on log L_{UV}, where alpha_{ox} becomes X-ray quiet more rapidly at higher log L_{UV}. We find no significant evidence for a relationship between the X-ray photon index, Gamma_X, and the UV luminosity, and we find marginally significant evidence that the X-ray continuum flattens with increasing z (2 sigma). The Gamma_X-z anti-correlation may be the result of X-ray spectral curvature, redshifting of a Compton reflection component into the observed Chandra band, and/or redshifting of a soft excess out of the observed Chandra band. Using the results for Gamma_X, we show that the alpha_{ox}-z relationship is unlikely to be a spurious result caused by redshifting of the observable X-ray spectral region. A correlation between alpha_{ox} and z implies evolution of the accretion process. We present a qualitative comparison of these new results with models for accretion disk emission.

We describe the general problem of estimating black hole masses of AGN by calculating the conditional probability distribution of M_BH given some set of observables. Special attention is given to the case where one uses the AGN continuum luminosity and emission line widths to estimate M_BH, and we outline how to set up the conditional probability distribution of M_BH given the observed luminosity, line width, and redshift. We show how to combine the broad line estimates of M_BH with information from an intrinsic correlation between M_BH and L, and from the intrinsic distribution of M_BH, in a manner that improves the estimates of M_BH. Simulation was used to assess how the distribution of M_BH inferred from the broad line mass estimates differs from the intrinsic distribution, and we find that this can lead to an inferred distribution that is too broad. We use these results and a sample of 25 sources that have recent reverberation mapping estimates of AGN black hole masses to investigate the effectiveness of using the C IV emission line to estimate M_BH and to indirectly probe the C IV region size--luminosity (R--L) relationship. We estimated M_BH from both C IV and H-Beta for a sample of 100 sources, including new spectra of 29 quasars. We find that the two emission lines give consistent estimates if one assumes R \propto L^{1/2}_{UV} for both lines.

We present long-term-stable and convergent evolutions of head-on black hole collisions and extraction of gravitational waves generated during the merger and subsequent ring-down. The new ingredients in this work are the use of fixed mesh-refinement and dynamical singularity excision techniques. We are able to carry out head-on collisions with large initial separations and demonstrate that our excision infrastructure is capable of accommodating the motion of the individual black holes across the computational domain as well as their their merger. We extract gravitational waves from these simulations using the Zerilli-Moncrief formalism and find the ring-down radiation to be, as expected, dominated by the l=2, m=0 quasi-normal mode. The total radiated energy is about 0.1 % of the total ADM mass of the system.

We describe the application of the `shapelet' linear decomposition of galaxy images to multi-wavelength morphological classification using the $u,g,r,i,$ and $z$-band images of 1519 galaxies from the Sloan Digital Sky Survey. We utilize elliptical shapelets to remove to first-order the effect of inclination on morphology. After decomposing the galaxies we perform a principal component analysis on the shapelet coefficients to reduce the dimensionality of the spectral morphological parameter space. We give a description of each of the first ten principal component's contribution to a galaxy's spectral morphology. We find that galaxies of different broad Hubble type separate cleanly in the principal component space. We apply a mixture of Gaussians model to the 2-dimensional space spanned by the first two principal components and use the results as a basis for classification. Using the mixture model, we separate galaxies into three classes and give a description of each class's physical and morphological properties. We find that the two dominant mixture model classes correspond to early and late type galaxies, respectively. The third class has, on average, a blue, extended core surrounded by a faint red halo, and typically exhibits some asymmetry. We compare our method to a simple cut on $u-r$ color and find the shapelet method to be superior in separating galaxies. Furthermore, we find evidence that the $u-r=2.22$ decision boundary may not be optimal for separation between early and late type galaxies, and suggest that the optimal cut may be $u-r \sim 2.4$.

We describe application of the `shapelet' linear decomposition of galaxy images to morphological classification using images of $\sim$ 3000 galaxies from the Sloan Digital Sky Survey. After decomposing the galaxies we perform a principal component analysis to reduce the number of dimensions of the shapelet space to nine. We find that each of these nine principal components contains unique morphological information, and give a description of each principal component's contribution to a galaxy's morphology. We find that galaxies of differing Hubble type separate cleanly in the shapelet space. We apply a Gaussian mixture model to the 9-dimensional space spanned by the principal components and use the results as a basis for classification. Using the mixture model, we separate galaxies into seven classes and give a description of each class's physical and morphological properties. We find that several of the mixture model classes correlate well with the traditional Hubble types both in their morphology and their physical parameters (e.g., color, velocity dispersions, etc.). In addition, we find an additional class of late-type morphology but with high velocity dispersions and very blue color; most of these galaxies exhibit post-starburst activity. This method provides an objective and quantitative alternative to traditional and subjective visual classification.

We introduce an algorithm for applying a cross-wavelet transform to analysis of quasiperiodic variations in a time-series, and introduce significance tests for the technique. We apply a continuous wavelet transform and the cross-wavelet algorithm to the Pearson-Readhead VLBI survey sources using data obtained from the University of Michigan 26-m parabloid at observing frequencies of 14.5, 8.0, and 4.8 GHz. Thirty of the sixty-two sources were chosen to have sufficient data for analysis, having at least 100 data points for a given time-series. Of these thirty sources, a little more than half exhibited evidence for quasiperiodic behavior in at least one observing frequency, with a mean characteristic period of 2.4 yr and standard deviation of 1.3 yr. We find that out of the thirty sources, there were about four time scales for every ten time series, and about half of those sources showing quasiperiodic behavior repeated the behavior in at least one other observing frequency.

A robust measurement of the clustering amplitude of the sub-mm population of starburst galaxies requires large-area surveys (>> 1 deg^2). The largest-format arrays subtend only 10 arcmin^2 on the sky and hence scan-mapping is a necessary observing mode. Providing realistic representations of the extragalactic sky and atmosphere, as the input to a detailed simulator of the telescope and instrument performance, allows important decisions to be made about the design of large-area fully-sampled surveys and observing strategies. In this paper we present preliminary simulations that include detector noise, time-constants and array geometry, telescope pointing errors, scan speeds and scanning angles, sky noise and sky rotation.