The Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) program monitors 849 active galactic nuclei (AGN) both spectroscopically and photometrically. The photometric observations used in this work span over four years and provide an excellent baseline for variability studies of these objects. We present the photometric light curves from 2014 to 2017 obtained by the Steward Observatory's Bok telescope and the CFHT telescope with MegaCam. We provide details on the data acquisition and processing of the data from each telescope, the difference imaging photometry used to produce the light curves, and the calculation of a variability index to quantify each AGN's variability. We find that the Welch-Stetson J-index provides a useful characterization of AGN variability and can be used to select AGNs for further study.

It is well known that reverberation mapping of active galactic nuclei (AGN) reveals a relationship between AGN luminosity and the size of the broad-line region, and that use of this relationship, combined with the Doppler width of the broad emission line, enables an estimate of the mass of the black hole at the center of the active nucleus based on a single spectrum. This has been discussed in numerous papers over the last two decades. An unresolved key issue is the choice of parameter used to characterize the line width; generally, most researchers use FWHM in favor of line dispersion (the square root of the second moment of the line profile) because the former is easier to measure, less sensitive to blending with other features, and usually can be measured with greater precision. However, use of FWHM introduces a bias, stretching the mass scale such that high masses are overestimated and low masses are underestimated. Here we investigate estimation of black hole masses in AGN based on individual or "single epoch" observations, with a particular emphasis in comparing mass estimates based on line dispersion and FWHM. We confirm the recent findings that, in addition to luminosity and line width, a third parameter is required to obtain accurate masses and that parameter seems to be Eddington ratio. We present simplified empirical formulae for estimating black hole masses from the Hbeta (4861 A) and C IV (1549 A) emission lines.

The size of the dust torus in Active Galactic Nuclei (AGN) and their high-luminosity counterparts, quasars, can be inferred from the time delay between UV/optical accretion disk continuum variability and the response in the mid-infrared (MIR) torus emission. This dust reverberation mapping (RM) technique has been successfully applied to $\sim 70$ $z\lesssim 0.3$ AGN and quasars. Here we present first results of our dust RM program for distant quasars covered in the SDSS Stripe 82 region combining $\sim 20$-yr ground-based optical light curves with 10-yr MIR light curves from the WISE satellite. We measure a high-fidelity lag between W1-band (3.4 $\mu$m) and $g$ band for 587 quasars over $0.3\lesssim z\lesssim 2$ ($\left<z\right>\sim 0.8$) and two orders of magnitude in quasar luminosity. They tightly follow (intrinsic scatter $\sim 0.17$ dex in lag) the IR lag-luminosity relation observed for $z<0.3$ AGN, revealing a remarkable size-luminosity relation for the dust torus over more than four decades in AGN luminosity, with little dependence on additional quasar properties such as Eddington ratio and variability amplitude. This study motivates further investigations in the utility of dust RM for cosmology, and strongly endorses a compelling science case for the combined 10-yr Vera C. Rubin Observatory Legacy Survey of Space and Time (optical) and 5-yr Nancy Grace Roman Space Telescope 2$\mu$m light curves in a deep survey for low-redshift AGN dust RM with much lower luminosities and shorter, measurable IR lags. The compiled optical and MIR light curves for 7,384 quasars in our parent sample are made public with this work.

We report the identification of a low-mass AGN, DES J0218$-$0430, in a redshift $z = 0.823$ galaxy in the Dark Energy Survey (DES) Supernova field. We select DES J0218$-$0430 as an AGN candidate by characterizing its long-term optical variability alone based on DES optical broad-band light curves spanning over 6 years. An archival optical spectrum from the fourth phase of the Sloan Digital Sky Survey shows both broad Mg II and broad H$\beta$ lines, confirming its nature as a broad-line AGN. Archival XMM-Newton X-ray observations suggest an intrinsic hard X-ray luminosity of $L_{{\rm 2-12\,keV}}\sim7.6\pm0.4\times10^{43}$ erg s$^{-1}$, which exceeds those of the most X-ray luminous starburst galaxies, in support of an AGN driving the optical variability. Based on the broad H$\beta$ from SDSS spectrum, we estimate a virial BH mass of $M_{\bullet}\approx10^{6.43}$-$10^{6.72}M_{\odot}$ (with the error denoting 1$\sigma$ statistical uncertainties only), consistent with the estimation from OzDES, making it the lowest mass AGN with redshift $>$ 0.4 detected in optical. We estimate the host galaxy stellar mass to be $M_{\ast}\sim10^{10.5\pm0.3}M_{\odot}$ based on modeling the multi-wavelength spectral energy distribution. DES J0218$-$0430 extends the $M_{\bullet}$-$M_{\ast}$ relation observed in luminous AGNs at $z\sim1$ to masses lower than being probed by previous work. Our work demonstrates the feasibility of using optical variability to identify low-mass AGNs at higher redshift in deeper synoptic surveys with direct implications for the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory.

Quasar broad emission lines are largely powered by photoionization from the accretion continuum. Increased central luminosity will enhance line emissivity in more distant clouds, leading to increased average distance of the broad-line-emitting clouds and decreased averaged line width, known as the broad-line region (BLR) "breathing". However, different lines breathe differently, and some high-ionization lines, such as C IV, can even show "anti-breathing" where the line broadens when luminosity increases. Using multi-year photometric and spectroscopic monitoring data from the Sloan Digital Sky Survey Reverberation Mapping project, we quantify the breathing effect ($\Delta$log W=$\alpha\Delta$log L) of broad H$\alpha$, H$\beta$, Mg II, C IV, and C III] for statistical quasar samples over $z\approx 0.1-2.5$. We found that H$\beta$ displays the most consistent normal breathing expected from the virial relation ($\alpha\sim-0.25$), Mg II and H$\alpha$ on average show no breathing ($\alpha\sim 0$), and C IV (and similarly C III] and Si IV mostly shows anti-breathing ($\alpha>0$). The anti-breathing of C IV can be well understood by the presence of a non-varying core component in addition to a reverberating broad-base component, consistent with earlier findings. The deviation from canonical breathing introduces extra scatter (a luminosity-dependent bias) in single-epoch virial BH mass estimates due to intrinsic quasar variability, which underlies the long argued caveats of C IV single-epoch masses. Using the line dispersion instead of FWHM leads to less, albeit still substantial, deviations from canonical breathing in most cases. Our results strengthen the need for reverberation mapping to provide reliable quasar BH masses, and quantify the level of variability-induced bias in single-epoch BH masses based on various lines.

We present the results of a pilot Hubble Space Telescope (HST) imaging study of the host galaxies of ten quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Probing more than an order of magnitude in BH and stellar masses, our sample is the first statistical sample to study the BH-host correlations beyond z>0.3 with reliable BH masses from reverberation mapping rather than from single-epoch spectroscopy. We perform image decomposition in two HST bands (UVIS-F606W and IR-F110W) to measure host colors and estimate stellar masses using empirical relations between broad-band colors and the mass-to-light ratio. The stellar masses of our targets are mostly dominated by a bulge component. The BH masses and stellar masses of our sample broadly follow the same correlations found for local RM AGN and quiescent bulge-dominant galaxies, with no strong evidence of evolution in the M_BH-M_*,bulge relation to z~0.6. We further compare the host light fraction from HST imaging decomposition to that estimated from spectral decomposition. We found a good correlation between the host fractions derived with both methods. However, the host fraction derived from spectral decomposition is systematically smaller than that from imaging decomposition by ~30%, indicating different systematics in both approaches. This study paves the way for upcoming more ambitious host galaxy studies of quasars with direct RM-based BH masses at high redshift.

We present optical light curves from the Transiting Exoplanet Survey Satellite (TESS) for the archetypical dwarf AGN in the nearby galaxy NGC 4395 hosting a $\sim 10^5\,M_\odot$ supermassive black hole. Significant variability is detected on timescales from weeks to hours before reaching the background noise level. The $\sim$month-long, 30-minute-cadence, high-precision TESS light curve can be well fit by a simple damped random walk model, with the damping timescale $\tau_{\rm DRW}$ constrained to be $2.3_{-0.7}^{+1.8}$~days ($1\sigma$). \obj\ lies almost exactly on the extrapolation of the $\tau_{\rm DRW}-M_{\rm BH}$ relation measured for AGN with BH masses that are more than three orders of magnitude larger. The optical variability periodogram can be well fit by a broken power-law with the high-frequency slope ($-1.88\pm0.15$) and the characteristic timescale ($\tau_{\rm br}\equiv 1/(2\pi f_{\rm br})=1.4_{-0.5}^{+1.9}\,$days) consistent with the DRW model within 1$\sigma$. This work demonstrates the power of TESS light curves in identifying low-mass accreting SMBHs with optical variability, and a potential global $\tau_{\rm DRW}-M_{\rm BH}$ relation that can be used to estimate SMBH masses with optical variability measurements.

We present reverberation mapping results for the MgII 2800 A broad emission line in a sample of 193 quasars at 0.35<z<1.7 with photometric and spectroscopic monitoring observations from the Sloan Digital Sky Survey Reverberation Mapping project during 2014 - 2017. We find significant time lags between the MgII and continuum lightcurves for 57 quasars and define a "gold sample" of 24 quasars with the most reliable lag measurements. We estimate false-positive rates for each lag that range from 1-24%, with an average false-positive rate of 11% for the full sample and 8% for the gold sample. There are an additional ~40 quasars with marginal MgII lag detections which may yield reliable lags after additional years of monitoring. The MgII lags follow a radius -- luminosity relation with a best-fit slope that is consistent with alpha=0.5 but with an intrinsic scatter of 0.36dex that is significantly larger than found for the Hb radius -- luminosity relation. For targets with SDSS-RM lag measurements of other emission lines, we find that our MgII lags are similar to the Hb lags and ~2-3 times larger than the CIV lags. This work significantly increases the number of MgII broad-line lags and provides additional reverberation-mapped black hole masses, filling the redshift gap at the peak of supermassive black hole growth between the Hb and CIV emission lines in optical spectroscopy.

We report on small-amplitude optical variability and recent dissipation of the unusually persistent broad emission lines in the blue compact dwarf galaxy PHL 293B. The galaxy's unusual spectral features (P Cygni-like profiles with $\sim$800 km s$^{-1}$ blueshifted absorption lines) have resulted in conflicting interpretations of the nature of this source in the literature. However, analysis of new Gemini spectroscopy reveals the broad emission has begun to fade after being persistent for over a decade prior. Precise difference imaging light curves constructed with the Sloan Digital Sky Survey and the Dark Energy Survey reveal small-amplitude optical variability of $\sim$0.1 mag in the g band offset by $100\pm21$ pc from the brightest pixel of the host. The light curve is well-described by an active galactic nuclei (AGN)-like damped random walk process. However, we conclude that the origin of the optical variability and spectral features of PHL 293B is due to a long-lived stellar transient, likely a Type IIn supernova or non-terminal outburst, mimicking long-term AGN-like variability. This work highlights the challenges of discriminating between scenarios in such extreme environments, relevant to searches for AGNs in dwarf galaxies. This is the second long-lived transient discovered in a blue compact dwarf, after SDSS1133. Our result implies such long-lived stellar transients may be more common in metal-deficient galaxies. Systematic searches for low-level variability in dwarf galaxies will be possible with the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory.

Results from a few decades of reverberation mapping (RM) studies have revealed a correlation between the radius of the broad-line emitting region (BLR) and the continuum luminosity of active galactic nuclei. This "radius-luminosity" relation enables survey-scale black-hole mass estimates across cosmic time, using relatively inexpensive single-epoch spectroscopy, rather than intensive RM time monitoring. However, recent results from newer reverberation mapping campaigns challenge this widely used paradigm, reporting quasar BLR sizes that differ significantly from the previously established radius-luminosity relation. Using simulations of the radius--luminosity relation with the observational parameters of SDSS-RM, we find that this difference is not likely due to observational biases. Instead, it appears that previous RM samples were biased to a subset of quasar properties, and the broader parameter space occupied by the SDSS-RM quasar sample has a genuinely wider range of BLR sizes. We examine the correlation between the deviations from the radius-luminosity relation and several quasar parameters; the most significant correlations indicate that the deviations depend on UV/optical SED and the relative amount of ionizing radiation. Our results indicate that single-epoch black-hole mass estimates that do not account for the diversity of quasars in the radius-luminosity relation could be overestimated by an average of ~0.3 dex.

We study the optical light curves - primarily probing the variable emission from the accretion disk - of ~ 900 extreme variability quasars (EVQs, with maximum flux variations more than 1 mag) over an observed-frame baseline of ~ 16 years using public data from the SDSS Stripe 82, PanSTARRS-1 and the Dark Energy Survey. We classify the multi-year long-term light curves of EVQs into three categories roughly in the order of decreasing smoothness: monotonic decreasing or increasing (3.7%), single broad peak and dip (56.8%), and more complex patterns (39.5%). The rareness of monotonic cases suggests that the major mechanisms driving the extreme optical variability do not operate over timescales much longer than a few years. Simulated light curves with a damped random walk model generally under-predict the first two categories with smoother long-term trends. Despite the different long-term behaviors of these EVQs, there is little dependence of the long-term trend on the physical properties of quasars, such as their luminosity, BH mass, and Eddington ratio. The large dynamic range of optical flux variability over multi-year timescales of these EVQs allows us to explore the ensemble correlation between the short-term (< 6 months) variability and the seasonal-average flux across the decade-long baseline (the rms-mean flux relation). We find that unlike the results for X-ray variability studies, the linear short-term flux variations do not scale with the seasonal-average flux, indicating different mechanisms that drive the short-term flickering and long-term extreme variability of accretion disk emission. Finally, we present a sample of 16 EVQs, where the approximately bell-shaped large amplitude variation in the light curve can be reasonably well fit by a simple microlensing model.

A comprehensive analysis of plasma and composition characteristics inside magnetic clouds (MCs) observed by the Advanced Composition Explorer (ACE) spacecraft from 1998 February to 2011 August is presented. The results show that MCs have specific interior structures, and MCs of different speeds show differences in composition and structure. Compared with the slow MCs, fast MCs have enhanced mean charge states of iron, oxygen, silicon, magnesium, $\mathrm{O^{7+}/O^{6+}}$, $\mathrm{C^{6+}/C^{5+}}$, $\mathrm{C^{6+}/C^{4+}}$ and $\mathrm{Fe^{\geq16+}/Fe_{total}}$ values. For ionic species in fast MCs, a higher atomic number represents a greater enhancement of mean charge state than slow MCs. We also find that both the fast and slow MCs display bimodal structure distribution in the mean iron charge state ($\mathrm{\langle Q\rangle Fe}$), which suggests that the existence of flux rope prior to the eruption is common. Furthermore, the $\mathrm{\langle Q\rangle Fe} $, $\mathrm{Fe^{\geq16+}/Fe_{total}}$, and $\mathrm{O^{7+}/O^{6+}}$ ratio distribution inside fast MCs have the feature that the posterior peak is higher than the anterior one. This result agrees with the "standard model" for CME/flares, by which magnetic reconnection occurs beneath the flux rope, thereby ionizing the ions of the posterior part of flux rope sufficiently by high-energy electron collisions or by direct heating in the reconnection region.

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$\sigma$ goal of 70 $\mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.

We report a detailed observational study of two quasi-periodic fast-propagating (QFP) magnetosonic wave events occurred on 2011 March 09 and 10, respectively. Interestingly, both the two events have two wave trains (WTs): one main and strong (WT-1) whereas the second appears small and weak (WT-2). Peculiar and common characteristics of the two events are observed, namely: 1) the two QFP waves are accompanied with brightenings during the whole stage of the eruptions; 2) both the two main wave trains are nearly propagating along the same direction; 3) EUV waves are found to be associated with the two events. Investigating various aspects of the target events, we argue that: 1) the second event is accompanied with a flux rope eruption during the whole stage; 2) the second event eruption produces a new filament-like (FL) dark feature; 3) the ripples of the two WT-2 QFP waves seem to result from different triggering mechanisms. Based on the obtained observational results, we propose that the funnel-like coronal loop system is indeed playing an important role in the two WT-1 QFP waves. The development of the second WT-2 QFP wave can be explained as due to the dispersion of the main EUV front. The co-existence of the two events offer thereby a significant opportunity to reveal what driving mechanisms and structures are tightly related to the waves.

We report the discovery of SDSS J0849+1114 as the first known triple Type 2 Seyfert nucleus. It represents three active black holes that are identified from new spatially resolved optical slit spectroscopy using the Dual Imaging Spectrograph on the 3.5 m telescope at the Apache Point Observatory. We also present new complementary observations including the Hubble Space Telescope Wide Field Camera 3 U- and Y-band imaging, Chandra Advanced CCD Imaging Spectrometer S-array X-ray 0.5--8 keV imaging spectroscopy, and NSF Karl G. Jansky Very Large Array radio 9.0 GHz imaging in its most extended A configuration. These comprehensive multiwavelength observations, when combined together, strongly suggest that all three nuclei are active galactic nuclei. While they are now still at kiloparsec-scale separations, where the host-galaxy gravitational potential dominates, the black holes may evolve into a bound triple system in $\lesssim$2 Gyr. These triple merger systems may explain the overly massive stellar cores that have been observed in some elliptical galaxies such as M87, which are expected to be unique gravitational wave sources. Similar systems may be more common in the early universe, when galaxy mergers are thought to have been more frequent.

The broad Mg II line in quasars has distinct variability properties compared with broad Balmer lines: it is less variable, and usually does not display a "breathing" mode, the increase in the average cloud distance when luminosity increases. We demonstrate that these variability properties of Mg II can be reasonably well explained by simple Locally Optimally Emitting Cloud (LOC) photoionization models, confirming earlier photoionization results. In the fiducial LOC model, the Mg II-emitting gas is on average more distant from the ionizing source than the H$\alpha$/H$\beta$ gas, and responds with a lower amplitude to continuum variations. If the broad-line region (BLR) is truncated at a physical radius of $\sim 0.3$ pc (for a $10^{8.5}M_{\odot}$ BH accreting at Eddington ratio of 0.1), most of the Mg II flux will always be emitted near this outer boundary and hence will not display breathing. These results indicate that reverberation mapping results on broad Mg II, while generally more difficult to obtain due to the lower line responsivity, can still be used to infer the Mg II BLR size and hence black hole mass. But it is possible that Mg II does not have a well defined intrinsic BLR size-luminosity relation for individual quasars, even though a global one for the general population may still exist. The dramatic changes in broad H$\alpha$/H$\beta$ emission in the observationally-rare changing-look quasars are fully consistent with photoionization responses to extreme continuum variability, and the LOC model provides natural explanations for the persistence of broad Mg II in changing-look quasars defined on H$\alpha$/H$\beta$, and the rare population of broad Mg II emitters in the spectra of massive inactive galaxies.

Multi-messenger astrophysics is a fast-growing, interdisciplinary field that combines data, which vary in volume and speed of data processing, from many different instruments that probe the Universe using different cosmic messengers: electromagnetic waves, cosmic rays, gravitational waves and neutrinos. In this Expert Recommendation, we review the key challenges of real-time observations of gravitational wave sources and their electromagnetic and astroparticle counterparts, and make a number of recommendations to maximize their potential for scientific discovery. These recommendations refer to the design of scalable and computationally efficient machine learning algorithms; the cyber-infrastructure to numerically simulate astrophysical sources, and to process and interpret multi-messenger astrophysics data; the management of gravitational wave detections to trigger real-time alerts for electromagnetic and astroparticle follow-ups; a vision to harness future developments of machine learning and cyber-infrastructure resources to cope with the big-data requirements; and the need to build a community of experts to realize the goals of multi-messenger astrophysics.

A solar jet on 2014 July 31, which was accompanied by a GOES C1.3 flare and a mini-filament eruption at the jet base, was studied by using observations taken by the New Vacuum Solar Telescope and the Solar Dynamic Observatory. Magnetic field extrapolation revealed that the jet was confined in a fan-spine magnetic system that hosts a null point at the height of about 9 Mm from the solar surface. An inner flare ribbon surrounded by an outer circular ribbon and a remote ribbon were observed to be associated with the eruption, in which the inner and remote ribbons respectively located at the footprints of the inner and outer spines, while the circular one manifested the footprint of the fan structure. It is interesting that the circular ribbon's west part showed an interesting round-trip slipping motion, while the inner ribbon and the circular ribbon's east part displayed a northward slipping motion. Our analysis results indicate that the slipping motions of the inner and the circular flare ribbons reflected the slipping magnetic reconnection process in the fan quasi-separatrix layer, while the remote ribbon was associated with the magnetic reconnection at the null point. In addition, the filament eruption was probably triggered by the magnetic cancellation around its south end, which further drove the slipping reconnection in the fan quasi-separatrix layer and the reconnection at the null point.

The Kelvin-Helmholtz (KH) instability is commonly found in many astrophysical, laboratory, and space plasmas. It could mix plasma components of different properties and convert dynamic fluid energy from large scale structure to smaller ones. In this study, we combined the ground-based New Vacuum Solar Telescope (NVST) and the Solar Dynamic Observatories (SDO) / Atmospheric Imaging Assembly (AIA) to observe the plasma dynamics associated with active region 12673 on 09 September 2017. In this multi-temperature view, we identified three adjacent layers of plasma flowing at different speeds, and detected KH instabilities at their interfaces. We could unambiguously track a typical KH vortex and measure its motion. We found that the speed of this vortex suddenly tripled at a certain stage. This acceleration was synchronized with the enhancements in emission measure and average intensity of the 193 \AA{} data. We interpret this as evidence that KH instability triggers plasma heating. The intriguing feature in this event is that the KH instability observed in the NVST channel was nearly complementary to that in the AIA 193 \AA{}. Such a multi-thermal energy exchange process is easily overlooked in previous studies, as the cold plasma component is usually not visible in the extreme ultraviolet channels that are only sensitive to high temperature plasma emissions. Our finding indicates that embedded cold layers could interact with hot plasma as invisible matters. We speculate that this process could occur at a variety of length scales and could contribute to plasma heating.

Off-nucleus active galactic nuclei (AGN) can be signposts of inspiraling supermassive black holes (SMBHs) on galactic scales, or accreting SMBHs recoiling after the coalescence of a SMBH binary or slingshot from three-body interactions. Because of the stochastic variability of AGN, the measured photocenter of an unresolved AGN-host system will display astrometric jitter that depends on the off-nucleus distance of the AGN, the total photometric variability of the system, and the AGN-host contrast. Here we use the precision astrometry from Gaia DR2 to constrain the off-nucleus population of a low-redshift (0.3<z<0.8) sample of unobscured broad-line AGN drawn from the SDSS with significant host contribution and photometric variability. We find that Gaia DR2 already provides strong constraints on the projected off-nucleus distance in the sub-kpc regime at these redshifts: 99%, 90% and 40% of AGN must be well-centered to <1 kpc, <500 pc and <100 pc, respectively. Limiting the sample to the most variable subset constrains >99% of AGN to be well-centered below 500 parsec. These results suggest that genuine off-nucleus AGN (offset by > a few hundred pc) must be rare at low redshift. Future Gaia releases of time series of photocenter and flux measurements, improved treatments for extended sources and longer baselines will further tighten these constraints, and enable a systematic full-sky search for rare off-nucleus AGN on ~ 10-1000 pc scales.

We present Chandra ACIS-S X-ray imaging spectroscopy for five dual active galactic nucleus (AGN) candidates. Our targets were drawn from a sample of 1286 [O III]-selected AGN pairs systematically selected from the Sloan Digital Sky Survey Seventh Data Release. Each of the targets contains two nuclei separated by $\sim$3-9 kpc in projection, both of which are optically classified as Type 2 (obscured) AGNs based on diagnostic ratios of the narrow emission lines. Combined with independent, empirical star formation rate estimates based on the host-galaxy stellar continua, the new Chandra X-ray observations allow us to evaluate the dual-AGN hypothesis for each merging system. We confirm two (SDSS J0907+5203 and SDSS J1544+0446) of the five targets as bona-fide dual AGNs. For the other three targets, the existing data are consistent with the dual-AGN scenario, but we cannot rule out the possibility of stellar/shock heating and/or one AGN ionizing both gaseous components in the merger. The average X-ray-to-[O III] luminosity ratio in our targets seems to be systematically smaller than that observed in single AGNs but is higher than that seen in dual AGNs selected from AGNs with double-peaked narrow emission lines. We suggest that the systematically smaller X-ray-to-[O III] luminosity ratio observed in dual AGNs than in single AGNs is due to a high nuclear gas column likely from strong merger-induced inflows. Unlike double-peaked-[O III]-selected dual AGNs, the new sample selected from resolved galaxy pairs are not subject to the orientation bias caused by the double-peak line-of-sight velocity splitting selection, which also contributes to lowering the X-ray-to-[O III] luminosity ratio.

We investigate the performance of different methodologies that measure the time lag between broad-line and continuum variations in reverberation mapping data using simulated light curves that probe a range of cadence, time baseline, and signal-to-noise ratio in the flux measurements. We compare three widely-adopted lag measuring methods: the Interpolated Cross-Correlation Function (ICCF), the z-transformed Discrete Correlation Function (ZDCF) and the MCMC code JAVELIN, for mock data with qualities typical of multi-object spectroscopic reverberation mapping (MOS-RM) surveys that simultaneously monitor hundreds of quasars. We quantify the overall lag detection efficiency, the rate of false detections, and the quality of lag measurements for each of these methods and under different survey designs (e.g., observing cadence and depth) using mock quasar light curves. Overall JAVELIN and ICCF outperform ZDCF in essentially all tests performed. Compared with ICCF, JAVELIN produces higher quality lag measurements, is capable of measuring more lags with timescales shorter than the observing cadence, is less susceptible to seasonal gaps and S/N degradation in the light curves, and produces more accurate lag uncertainties. We measure the Hbeta broad-line region size-luminosity (R-L) relation with each method using the simulated light curves to assess the impact of selection effects of the design of MOS-RM surveys. The slope of the R-L relation measured by JAVELIN is the least biased among the three methods, and is consistent across different survey designs. These results demonstrate a clear preference for JAVELIN over the other two non-parametric methods for MOS-RM programs, particularly in the regime of limited light curve quality as expected from most MOS-RM programs.

The ubiquitous solar jets or jet-like activities are generally regarded as an important source of energy and mass input to the upper solar atmosphere and the solar wind. However, questions about their triggering and driving mechanisms are not completely understood. By taking advantage of high temporal and high spatial resolution stereoscopic observations taken by the Solar Dynamic Observatory (SDO) and the Solar Terrestrial Relations Observatory (STEREO), we report an intriguing two-sided-loop jet occurred on 2013 June 02, which was dynamically associated with the eruption of a mini-filament below an overlying large filament, and two distinct reconnection processes are identified during the formation stage. The SDO observations reveals that the two-sided-loop jet showed a concave shape with a projection speed of about 80 - 136. From the other view angle, the STEREO ahead observations clearly showed that the trajectory of the two arms of the two-sided-loop were along the cavity magnetic field lines hosting the large filament. Contrary to the well-accepted theoretical model, the present observation sheds new light on our understanding of the formation mechanism of two-sided-loop jets. Moreover, the eruption of the two-sided-loop jet not only supplied mass to the overlying large filament, but also provided a rare opportunity to diagnose the magnetic structure of the overlying large filament via the method of three-dimensional reconstruction.

Gaia's milli-arcsec (mas) astrometric precision allows systematic identification of optically-selected sub-kpc dual active galactic nuclei (AGN), off-nucleus AGN, and small-scale lensed quasars by `varstrometry' -- where variability-induced astrometric jitter, i.e., temporal displacements of photocenter in unresolved sources, can be reasonably well detected or constrained. This approach extends systematic searches for small-scale ($\gtrsim$ mas) dual and off-nucleus AGN to poorly explored regime between $\sim 10$ pc and $\sim 1$ kpc, with Gaia's full sky coverage and depth to $G\sim 21$. We outline the general principles of this method and calculate the expected astrometric signals from the full time series of photocenter measurements and light curves. We demonstrate the feasibility of varstrometry by using Gaia DR2 data on a sample of variable pre-main sequence stars with known close companions. We find that extended host galaxies have a significant impact on the accuracy of astrometric and photometric variability in Gaia DR2, a situation to be improved in future Gaia releases. Using spectroscopically confirmed SDSS quasars, we present several examples of candidate sub-kpc off-nucleus or dual AGN selected from Gaia DR2. We discuss the merits and limitations of this method and follow-up strategy for promising candidates. We highlight Gaia's potential of systematically discovering and characterizing the sub-kpc off-nucleus and dual AGN population in the entire optical sky.

We investigate the effects of extended multi-year light curves (9-year photometry and 5-year spectroscopy) on the detection of time lags between the continuum variability and broad-line response of quasars at z>~1.5, and compare with the results using 4-year photometry+spectroscopy presented in a companion paper. We demonstrate the benefits of the extended light curves in three cases: (1) lags that are too long to be detected by the shorter-duration data but can be detected with the extended data; (2) lags that are recovered by the extended light curves but missed in the shorter-duration data due to insufficient light curve quality; and (3) lags for different broad line species in the same object. These examples demonstrate the importance of long-term monitoring for reverberation mapping to detect lags for luminous quasars at high-redshift, and the expected performance of the final dataset from the Sloan Digital Sky Survey Reverberation Mapping project that will have 11-year photometric and 7-year spectroscopic baselines.

We report the discovery of an exceptional MIR flare in a Type 2 AGN, SDSS J165726.81+234528.1, at $z=0.059$. This object brightened by 3 mag in the Wide-field Infrared Survey Explorer (WISE) $W1$ and $W2$ bands between 2015 and 2017 (and is fading since 2018), without significant changes ($\lesssim$ 0.2 mag) in the optical over the same period of time. Based on the WISE light curves and near-IR imaging, the flare is more significant at longer wavelengths, suggesting an origin of hot dust emission. The estimated black hole mass ($\sim 10^{6.5}\,M_\odot$) from different methods places its peak bolometric luminosity around the Eddington limit. The high luminosity of the MIR flare and its multi-year timescale suggest that it most likely originated from reprocessed dust radiation in an extended torus surrounding the AGN, instead of from stellar explosions. The MIR color variability is consistent with known changing-look AGN and tidal disruption events (TDEs), but inconsistent with normal supernovae. We suggest that it is a turning-on Type 2 AGN or TDE, where the optical variability is obscured by the dust torus during the transition. This MIR flare event reveals a population of dramatic nuclear transients that are missed in the optical.

The Maunakea Spectroscopic Explorer is a next-generation massively multiplexed spectroscopic facility currently under development in Hawaii. It is completely dedicated to large-scale spectroscopic surveys and will enable transformative science. In this white paper we summarize the science case and describe the current state of the project.

We report the formation of a complicated coronal mass ejection (CME) on 2015 August 23 by using the high temporal and high spatial resolution multi-wavelength observations taken by the Solar Dynamic Observatory and the Solar and Heliospheric Observatory. The CME exhibited both jet-like and bubble-like components simultaneously, and therefore we call it a twin CME. Detailed imaging and kinematic analysis results indicate that the twin CME were evolved form the eruption of a mini-filament driven blowout jet at the east edge of an equatorial coronal hole, in which the activation of the mini-filament was tightly associated with the continuous flux cancellation and quasi-periodic jet-like activities in the filament channel. Due to the magnetic reconnection between the filament and the ambient open field lines, the filament broke partially at the northern part and resulted in an intriguing blowout jet in the south direction. It is interesting that the ejecting jet was deflected by a group of remote open field lines, which resulted in the significant direction change of the jet from southward to eastward. Based on the close temporal and spatial relationships among the jet, filament eruption, and the twin CME, we conclude that the jet-like CME should be the coronal extension of the jet plasma, while the bubble-like one should be originated from the eruption of the mini-filament confined by the closed magnetic fields at the jet-base.

We present accretion-disk structure measurements from continuum lags in the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Lags are measured using the \texttt{JAVELIN} software from the first-year SDSS-RM $g$ and $i$ photometry, resulting in well-defined lags for 95 quasars, 33 of which have lag SNR $>$ 2$\sigma$. We also estimate lags using the \texttt{ICCF} software and find consistent results, though with larger uncertainties. Accretion-disk structure is fit using a Markov Chain Monte Carlo approach, parameterizing the measured continuum lags as a function of disk size normalization, wavelength, black hole mass, and luminosity. In contrast with previous observations, our best-fit disk sizes and color profiles are consistent (within 1.5~$\sigma$) with the \citet{SS73} analytic solution. We also find that more massive quasars have larger accretion disks, similarly consistent with the analytic accretion-disk model. The data are inconclusive on a correlation between disk size and continuum luminosity, with results that are consistent with both no correlation and with the \citet{SS73} expectation. The continuum lag fits have a large excess dispersion, indicating that our measured lag errors are underestimated and/or our best-fit model may be missing the effects of orientation, spin, and/or radiative efficiency. We demonstrate that fitting disk parameters using only the highest-SNR lag measurements biases best-fit disk sizes to be larger than the disk sizes recovered using a Bayesian approach on the full sample of well-defined lags.

The width of the broad emission lines in quasars is commonly characterized either by the full-width-at-half-maximum (FWHM) or the square root of the second moment of the line profile ($\sigma_{\rm line}$), and used as an indicator of the virial velocity of the broad-line region (BLR) in the estimation of black hole (BH) mass. We measure FWHM and $\sigma_{\rm line}$ for H$\alpha$, H$\beta$ and Mg II broad lines in both the mean and root-mean-square (rms) spectra of a large sample of quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. We introduce a new quantitative recipe to measure $\sigma_{\rm line}$ that is reproducible, less susceptible to noise and blending in the wings, and scales with the intrinsic width of the line. We compare the four definitions of line width (FWHM and $\sigma_{\rm line}$ in mean and rms spectra, respectively) for each of the three broad lines and among different lines. There are strong correlations among different width definitions for each line, providing justification for using the line width measured in single-epoch spectroscopy as a virial velocity indicator. There are also strong correlations among different lines, suggesting alternative lines to H$\beta$ can be used to estimate virial BH masses. We further investigate the correlations between virial BH masses using different line width definitions and the stellar velocity dispersion of the host galaxies, and the dependence of line shape (characterized by the ratio FWHM/$\sigma_{\rm line}$) on physical properties of the quasar. Our results provide further evidence that FWHM is more sensitive to the orientation of a flattened BLR geometry than $\sigma_{\rm line}$, but the overall comparison between the virial BH mass and host stellar velocity dispersion does not provide conclusive evidence that one particular width definition is significantly better than the others.

We analyze extensive spectroscopic and photometric data of the hypervariable quasar SDSS J131424+530527 (RMID 017) at z=0.456, an optical "changing look" quasar from the Sloan Digital Sky Survey Reverberation Mapping project that increased in optical luminosity by a factor of 10 between 2014 and 2017. The observed broad emission lines all respond in luminosity and width to the changing optical continuum, as expected for photoionization in a stratified, virialized broad emission line region. The luminosity changes therefore result from intrinsic changes in accretion power rather than variable obscuration. The variability is continuous and apparently stochastic, disfavoring an origin as a discrete event such as a tidal disruption flare or microlensing event. It is coordinated on day timescales with blue leading red, consistent with reprocessing powering the entire optical SED. We show that this process cannot work in a standard thin disk geometry on energetic grounds, and would instead require a large covering factor reprocessor. Disk instability models could potentially also explain the data, provided that the instability sets in near the inner radius of a geometrically thick accretion disk.

We present reverberation-mapping lags and black-hole mass measurements using the CIV 1549 broad emission line from a sample of 349 quasars monitored as a part of the Sloan Digital Sky Survey Reverberation Mapping Project. Our data span four years of spectroscopic and photometric monitoring for a total baseline of 1300 days. We report significant time delays between the continuum and the CIV 1549 emission line in 52 quasars, with an estimated false-positive detection rate of 10%. Our analysis of marginal lag measurements indicates that there are on the order of 100 additional lags that should be recoverable by adding more years of data from the program. We use our measurements to calculate black-hole masses and fit an updated CIV radius-luminosity relationship. Our results significantly increase the sample of quasars with CIV RM results, with the quasars spanning two orders of magnitude in luminosity toward the high-luminosity end of the CIV radius-luminosity relation. In addition, these quasars are located at among the highest redshifts (z~1.4-2.8) of quasars with black hole masses measured with reverberation mapping. This work constitutes the first large sample of CIV reverberation-mapping measurements in more than a dozen quasars, demonstrating the utility of multi-object reverberation mapping campaigns.

Most bulge-dominated galaxies host black holes with masses that tightly correlate with the masses of their bulges. This may indicate that the black holes may regulate galaxy growth or vice versa, or that they may grow in lock-step. The quest to understand how, when, and where those black-holes formed motivates much of extragalactic astronomy. Here we focus on a population of galaxies with active black holes in their nuclei (active galactic nuclei or AGN), that are fully or partially hidden by dust and gas: the emission from the broad line region is either completely or partially obscured with a visual extinction of 1 or above. This limit, though not yet precise, appears to be the point at which the populations of AGN may evolve differently. We highlight the importance of finding and studying those dusty AGN at redshifts between 1 and 3, the epoch when the universe may have gone through its most dramatic changes. We emphasize the need for future large multiplexed spectroscopic instruments that can perform dedicated surveys in the optical and NIR to pin down the demographics of such objects and study their reddening properties, star-formation histories, and excitation conditions. These key studies will shed light on the role of black holes in galaxy evolution during the epoch of peak growth activity.

In this paper,we report a detailed analysis of recurrent jets originated from a location with emerging, canceling and converging negative magnetic field at the east edge of NOAA active region AR11166 from 2011 March 09 to 10. The event presented several interesting features. First, a satellite sunspot appeared and collided with a pre-existing opposite polarity magnetic field and caused a recurrent solar jet event. Second, the evolution of the jets showed blowout-like nature and standard characteristics. Third, the satellite sunspot exhibited a motion toward southeast of AR11166 and merged with the emerging flux near the opposite polarity sunspot penumbra, which afterward, due to flux convergence and cancellation episodes, caused recurrent jets. Fourth, three of the blowout jets associated with coronal mass ejections (CMEs), were observed from field of view of the Solar Terrestrial Relations Observatory. Fifth, almost all the blowout jet eruptions were accompanied with flares or with more intense brightening in the jet base region, while almost standard jets did not manifest such obvious feature during eruptions. The most important, the blowout jets were inclined to faster and larger scale than the standard jets. The standard jets instead were inclined to relative longer-lasting. The obvious shearing and twisting motions of the magnetic field may be interpreted as due to the shearing and twisting motions for a blowout jet eruption. From the statistical results, about 30% blowout jets directly developed into CMEs. It suggests that the blowout jets and CMEs should have a tight relationship.

We present new Gemini/GMOS optical spectroscopy of 16 extreme variability quasars (EVQs) that dimmed by more than 1.5 mag in the $g$ band between the Sloan Digital Sky Survey (SDSS) and the Dark Energy Survey (DES) epochs (separated by a few years in the quasar rest frame). The quasar sample covers a redshift range of $0.5 < z < 2.1$. Nearly half of these EVQs brightened significantly (by more than 0.5 mag in the $g$ band) in a few years after reaching their previous faintest state, and some EVQs showed rapid (non-blazar) variations of greater than 1-2 mag on timescales of only months. Leveraging on the large dynamic range in continuum variability between the earlier SDSS and the new GMOS spectra, we explore the associated variations in the broad Mg II,$\lambda2798$ line, whose variability properties have not been well studied before. The broad Mg II flux varies in the same direction as the continuum flux, albeit with a smaller amplitude, which indicates at least some portion of Mg II is reverberating to continuum changes. However, the width (FWHM) of Mg II does not vary accordingly as continuum changes for most objects in the sample, in contrast to the case of the broad Balmer lines. Using the width of broad Mg II to estimate the black hole mass therefore introduces a luminosity-dependent bias.

We present a detailed characterization of the 849 broad-line quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Our quasar sample covers a redshift range of 0.1<z<4.5 and is flux-limited to i_PSF<21.7 without any other cuts on quasar properties. The main sample characterization includes: 1) spectral measurements of the continuum and broad emission lines for individual objects from the coadded first-season spectroscopy in 2014; 2) identification of broad and narrow absorption lines in the spectra; 3) optical variability properties for continuum and broad lines from multi-epoch spectroscopy. We provide improved systemic redshift estimates for all quasars, and demonstrate the effects of signal-to-noise ratio on the spectral measurements. We compile measured properties for all 849 quasars along with supplemental multi-wavelength data for subsets of our sample from other surveys. The SDSS-RM sample probes a diverse range in quasar properties, and shows well detected continuum and broad-line variability for many objects from first-season monitoring data. The compiled properties serve as the benchmark for follow-up work based on SDSS-RM data. The spectral fitting tools are made public along with this work.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

Time-domain Astrophysics (TDA), a foundation of Astronomy, has become a major part of current and projected (2020s) astrophysics. While much has been derived from temporal measures of flux and color, the real physics comes from spectroscopy. With LSST coming on line in 2022, with TDA as one of its original drivers, the deluge of Transients and new types of variables will be truly astronomical. With multi-wavelength targeted EM surveys and multi-messenger (e.g. LIGO-international and advanced ICECUBE), and the possibility of full-sky and full-time X-ray imagers, the discovery of new Transients and Variables will flood telescopes on the ground and in space, and this just for multi-band imaging without spectroscopy. In this White Paper we briefly summarize several long-standing major science objectives that can be realized with TDA imaging and spectroscopy (near UV to mid IR) from space. We provide a brief description of how these can be achieved with the Time-domain Spectroscopic Observatory (TSO), a Probe-Class mission concept that ELTs on the ground and Flagship missions in space can not achieve on their own.

Long gamma-ray bursts (GRBs) are the most luminous known electromagnetic radiation sources in the Universe for the 3 to 300 sec of their prompt flashes (isotropic X/ gamma-ray luminosities up to 10^53 ergs/sec). Their afterglows have first day rest-frame UV/optical absolute magnitudes AB = -30 to -23. This luminous continuum nUV-nIR back-light provides the ultimate probe of the SFR(z) back to the first Pop III to II.5 stars, expected to be massive and GRB progenitors. GRB afterglow spectra in the first 1 to 3 hours will directly measure their host galaxy ionization fraction x_i vs. z in the Epoch of Reionization (EOR), tracing the growth of structure. Only 28% of Swift GRBs have measured redshifts due to limited followup at R, J >21. Some ~25% of GRBs are optically dark due to dust absorption in their host galaxies, but those with low NH in their X-ray spectra are likely at z >7. Current 8-10m telescopes and coming ELTs cannot pursue optically dark GRBs promptly, nor can JWST or WFIRST slew within 0.5 to 1 days of a GRB. The Time-domain Spectroscopic Observatory (TSO) is a proposed Probe-class 1.3m telescope at L2, with imaging and spectroscopy (R = 200, 1800) in 4 bands (0.3 to 5um) and rapid slew capability to 90% of sky. TSO would finally utilize z > 6 to 12 GRBs as the most direct probe of the SFR(z), EOR(z), and possibly the first direct detection of the core collapse of the very first (Pop III) stars.

With the epochal first detection of gravitational waves from a binary neutron star (NS) merger with the GW170817 event, and its direct confirmation that NS-NS mergers are significant sources of the of the r-process nucleosynthesis of heavy elements, an immense new arena for prompt EM (X-rays through IR and radio) studies of fundamental physics has been opened. Over the next decade, GW observatories will expand in scale and sensitivity so the need for facilities that can provide prompt, high sensitivity, broad-band EM followup becomes more urgent. NS-NS or NS-black hole (BH) mergers will be instantly recognized (and announced) by the LIGO-international collaboration. LSST will be a prime resource for rapid tiling of what will usually be large (~10-100 degree squared) error boxes. X-ray through IR Telescopes in space with (nearly) full-sky access that can rapidly image and tile are crucial for providing the earliest imaging and spectroscopic studies of the kilonova emission immediately following NS-NS mergers. The Time-domain Spectroscopic Observatory (TSO) is a proposed Probe-class 1.3 m telescope at L2, with imaging and spectroscopy (R = 200, 1800) in 4 bands (0.3 - 5 micron) and rapid slew capability to 90% of sky. TSO nUV-mid-IR spectra will enable new constraints on NS structure and nucleosynthesis.

The ubiquitous variability of quasars across a wide range of wavelengths and timescales encodes critical information about the structure and dynamics of the circumnuclear emitting regions that are too small to be directly resolved, as well as detailed underlying physics of accretion and feedback processes in these active supermassive black holes. We emphasize the importance of studying quasar variability with time-domain spectroscopy, focusing on two science cases: (1) reverberation mapping (RM) to measure the broad-line region sizes and black hole masses in distant quasars; (2) spectroscopic follow-up of extreme variability quasars that dramatically change their continuum and broad-line flux within several years. We highlight the need for dedicated optical-infrared spectroscopic survey facilities in the coming decades to accompany wide-area time-domain imaging surveys, including: (1) the next phase of the Sloan Digital Sky Survey (SDSS-V; ~2020-2025), an all-sky, time-domain multi-object spectroscopic survey with 2.5m-class telescopes; (2) the planned Maunakea Spectroscopic Explorer, a dedicated 10m-class spectroscopic survey telescope with a 1.5 sq. deg field-of-view and multiplex of thousands of fibers in both optical and near-IR (J+H) to begin operations in 2029; (3) the Time-domain Spectroscopic Observatory (TSO), a proposed Probe-class ~1.3m telescope at L2, with imaging and spectroscopy (R=200, 1800) in 4 bands (0.3 - 5 micron) and rapid slew capability to 90% of sky, which will extend the coverage of Hbeta to z=8.

The discovery of luminous quasars at redshifts up to 7.5 demonstrates the existence of several billion M_sun supermassive black holes (SMBHs) less than a billion years after the Big Bang. They are accompanied by intense star formation in their host galaxies, pinpointing sites of massive galaxy assembly in the early universe, while their absorption spectra reveal an increasing neutral intergalactic medium (IGM) at the epoch of reionization. Extrapolating from the rapid evolution of the quasar density at z=5-7, we expect that there is only one luminous quasar powered by a billion M_sun SMBH in the entire observable universe at z~9. In the next decade, new wide-field, deep near-infrared (NIR) sky surveys will push the redshift frontier to the first luminous quasars at z~9-10; the combination with new deep X-ray surveys will probe fainter quasar populations that trace earlier phases of SMBH growth. The identification of these record-breaking quasars, and the measurements of their BH masses and accretion properties require sensitive spectroscopic observations with next generation of ground-based and space telescopes at NIR wavelengths. High-resolution integral-field spectroscopy at NIR, and observations at millimeter and radio wavelengths, will together provide a panchromatic view of the quasar host galaxies and their galactic environment at cosmic dawn, connecting SMBH growth with the rise of the earliest massive galaxies. Systematic surveys and multiwavelength follow-up observations of the earliest luminous quasars will strongly constrain the seeding and growth of the first SMBHs in the universe, and provide the best lines of sight to study the history of reionization.

This report provides an overview of recent work that harnesses the Big Data Revolution and Large Scale Computing to address grand computational challenges in Multi-Messenger Astrophysics, with a particular emphasis on real-time discovery campaigns. Acknowledging the transdisciplinary nature of Multi-Messenger Astrophysics, this document has been prepared by members of the physics, astronomy, computer science, data science, software and cyberinfrastructure communities who attended the NSF-, DOE- and NVIDIA-funded "Deep Learning for Multi-Messenger Astrophysics: Real-time Discovery at Scale" workshop, hosted at the National Center for Supercomputing Applications, October 17-19, 2018. Highlights of this report include unanimous agreement that it is critical to accelerate the development and deployment of novel, signal-processing algorithms that use the synergy between artificial intelligence (AI) and high performance computing to maximize the potential for scientific discovery with Multi-Messenger Astrophysics. We discuss key aspects to realize this endeavor, namely (i) the design and exploitation of scalable and computationally efficient AI algorithms for Multi-Messenger Astrophysics; (ii) cyberinfrastructure requirements to numerically simulate astrophysical sources, and to process and interpret Multi-Messenger Astrophysics data; (iii) management of gravitational wave detections and triggers to enable electromagnetic and astro-particle follow-ups; (iv) a vision to harness future developments of machine and deep learning and cyberinfrastructure resources to cope with the scale of discovery in the Big Data Era; (v) and the need to build a community that brings domain experts together with data scientists on equal footing to maximize and accelerate discovery in the nascent field of Multi-Messenger Astrophysics.

We report initial results from a large Gemini program to observe z>~5.7 quasars with GNIRS near-IR spectroscopy. Our sample includes 50 quasars with simultaneous ~0.85-2.5 micron spectra covering the rest-frame ultraviolet and major broad emission lines from Ly-alpha to MgII. We present spectral measurements for these quasars and compare to their lower-redshift counterparts at z=1.5-2.3. We find that when quasar luminosity is matched, there are no significant differences between the rest-UV spectra of z>~5.7 quasars and the low-z comparison sample. High-z quasars have similar continuum and emission line properties and occupy the same region in the black hole mass and luminosity space as the comparison sample, accreting at an average Eddington ratio of ~0.3. There is no evidence for super-Eddington accretion or hypermassive (>10^10 Msun) black holes within our sample. We find a mild excess of quasars with weak CIV lines relative to the control sample. Our results, corroborating earlier studies but with better statistics, demonstrate that these high-z quasars are already mature systems of accreting supermassive black holes operating with the same physical mechanisms as those at lower redshifts.

We report the first unambiguous quasi-periodic large-scale extreme-ultraviolet (EUV) wave or shock that was detected by the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory. During the whip-like unwinding eruption of a small filament on 2012 April 24, multiple consecutive large-scale wavefronts emanating from AR11467 were observed simultaneously along the solar surface and a closed transequatorial loop system. In the meantime, an upward propagating dome-shaped wavefront was also observed, whose initial speed and deceleration are about 1392 km/s and 1.78 km/s^2, respectively. Along the solar surface, the quasi-peridoic wavefronts had a period of about 163 +/- 21 seconds and propagated at a nearly constant speed of 747 +/- 26 km/s; they interacted with active region AR11469 and launched a sympathetic upward propagating secondary EUV wave. The wavefronts along the loop system propagated at a speed of 897 km/s, and they were reflected back at the southern end of the loop system at a similar speed. In addition to the propagating waves, a standing kink wave was also present in the loop system simultaneously. Periodicity analysis reveals that the period of the wavefronts was consistent with that of the unwinding helical structures of the erupting filament. Based on these observational facts, we propose that the observed quasi-periodic EUV wavefronts were most likely excited by the periodic unwinding motion of the filament helical structures. In addition, two different seismological methods are applied to derive the magnetic field strength of the loop system, and for the first time the reliability of these inversion techniques are tested with the same magnetic structure.

Quasi-periodic fast-propagating (QFP) magnetosonic waves and extreme ultraviolet (EUV) waves were proposed to be driven by solar flares and coronal mass ejections (CMEs), respectively. In this Letter, we present a detailed analysis of an interesting event in which we find that both QFP magnetosonic waves and EUV waves are excited simultaneously in one solar eruption event. The co-existence of the two wave phenomena offers an excellent opportunity to explore their driving mechanisms. The QFP waves propagate in a funnel-like loop system with a speed of 682--837 \speed{} and a lifetime of 2 minutes. On the contrary, the EUV waves, which present a faster component and a slower component, propagate in a wide angular extent, experiencing reflection and refraction across a magnetic quasi-separatrix layer. The faster component of the EUV waves travels with a speed of 412--1287 \speed{}, whereas the slower component travels with a speed of 246--390 \speed{}. The lifetime of the EUV waves is $\sim$15 minutes. It is revealed that the faster component of the EUV waves is cospatial with the first wavefront of the QFP wave train. Besides, The QFP waves have a period of about 45$\pm$5 seconds, which is absent in the associated flares. All these results imply that QFP waves can also be excited by mass ejections, including CMEs or jets.

In this paper, we present a detailed analysis of a coronal blowout jet eruption which was associated with an obvious extreme-ultraviolet (EUV) wave and one complicated coronal mass ejection (CME) event based on the multi-wavelength and multi-view-angle observations from {\sl Solar Dynamics Observatory} and {\sl Solar Terrestrial Relations Observatory}. It is found that the triggering of the blowout jet was due to the emergence and cancellation of magnetic fluxes on the photosphere. During the rising stage of the jet, the EUV wave appeared just ahead of the jet top, lasting about 4 minutes and at a speed of 458 - \speed{762}. In addition, obvious dark material is observed along the EUV jet body, which confirms the observation of a mini-filament eruption at the jet base in the chromosphere. Interestingly, two distinct but overlapped CME structures can be observed in corona together with the eruption of the blowout jet. One is in narrow jet-shape, while the other one is in bubble-shape. The jet-shaped component was unambiguously related with the outwardly running jet itself, while the bubble-like one might either be produced due to the reconstruction of the high coronal fields or by the internal reconnection during the mini-filament ejection according to the double-CME blowout jet model firstly proposed by Shen et al. (2012b), suggesting more observational evidence should be supplied to clear the current ambiguity based on large samples of blowout jets in future studies.

This white paper specifies the footprints, cadence requirements, and total-depth requirements needed to allow the most-successful AGN studies in the four currently selected LSST Deep-Drilling Fields (DDFs): ELAIS-S1, XMM-LSS, CDF-S, and COSMOS. The information provided on cadence and total-depth requirements will also likely be applicable to enabling effective AGN science in any additional DDFs that are chosen.

We systematically investigate short-timescale ($<$10-day rest-frame) CIV broad absorption-line (BAL) variability to constrain quasar-wind properties and provide insights into BAL-variability mechanisms in quasars. We employ data taken by the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, as the rapid cadence of these observations provides a novel opportunity to probe BAL variability on shorter rest-frame timescales than have previously been explored. In a sample of 27 quasars with a median of 58 spectral epochs per quasar, we have identified 15 quasars ($55^{+18}_{-14}$%), 19 of 37 CIV BAL troughs ($51^{+15}_{-12}$%), and 54 of 1460 epoch pairs ($3.7 \pm 0.5$%) that exhibit significant CIV BAL equivalent-width variability on timescales of less than 10 days in the quasar rest frame. These frequencies indicate that such variability is common among quasars and BALs, though somewhat rare among epoch pairs. Thus, models describing BALs and their behavior must account for variability on timescales down to less than a day in the quasar rest frame. We also examine a variety of spectral characteristics and find that in some cases, BAL variability is best described by ionization-state changes, while other cases are more consistent with changes in covering fraction or column density. We adopt a simple model to constrain the density and radial distance of two outflows appearing to vary by ionization-state changes, yielding outflow density lower limits consistent with previous work.

Quasars whose broad emission lines show temporal, bulk radial velocity (RV) shifts have been proposed as candidate sub-parsec (sub-pc), binary supermassive black holes (BSBHs). We identified a sample of 16 BSBH candidates based on two-epoch spectroscopy among 52 quasars with significant RV shifts over a few rest-frame years. The candidates showed consistent velocity shifts independently measured from two broad lines (H$\beta$ and H$\alpha$ or Mg${\rm \,II}$) without significant changes in the broad-line profiles. Here in the third paper of the series, we present further third- and fourth-epoch spectroscopy for 12 of the 16 candidates for continued RV tests, spanning $\sim$5--15 yr in the quasars' rest frames. Cross-correlation analysis of the broad H$\beta$ calibrated against [O${\rm\,III}]\,\lambda 5007$ suggests that 5 of the 12 quasars remain valid as BSBH candidates. They show broad H$\beta$ RV curves that are consistent with binary orbital motion without significant changes in the broad line profiles. Their broad H$\alpha$ (or Mg${\rm \,II}$) lines display RV shifts that are either consistent with or smaller than those seen in broad H$\beta$. The RV shifts can be explained by a $\sim$0.05--0.1 pc BSBH with an orbital period of $\sim$40--130 yr, assuming a mass ratio of 0.5--2 and a circular orbit. However, the parameters are not well constrained given the few epochs that sample only a small portion of the hypothesized binary orbital cycle. The apparent occurrence rate of sub-pc BSBHs is $\lesssim$13$\pm$5% among all SDSS quasars, with no significant difference in the subsets with and without single-epoch broad line velocity offsets. Dedicated long-term spectroscopic monitoring is still needed to further confirm or reject these BSBH candidates.

We present reverberation mapping results from the first year of combined spectroscopic and photometric observations of the Sloan Digital Sky Survey Reverberation Mapping Project. We successfully recover reverberation time delays between the $g+i$-band emission and the broad H$\beta$ emission line for a total of 44 quasars, and for the broad H$\alpha$ emission line in 18 quasars. Time delays are computed using the JAVELIN and CREAM software and the traditional interpolated cross-correlation function (ICCF): Using well defined criteria, we report measurements of 32 H$\beta$ and 13 H$\alpha$ lags with JAVELIN, 42 H$\beta$ and 17 H$\alpha$ lags with CREAM, and 16 H$\beta$ and 8 H$\alpha$ lags with the ICCF. Lag values are generally consistent among the three methods, though we typically measure smaller uncertainties with JAVELIN and CREAM than with the ICCF, given the more physically motivated light curve interpolation and more robust statistical modeling of the former two methods. The median redshift of our H$\beta$-detected sample of quasars is 0.53, significantly higher than that of the previous reverberation-mapping sample. We find that in most objects, the time delay of the H$\alpha$ emission is consistent with or slightly longer than that of H$\beta$. We measure black hole masses using our measured time delays and line widths for these quasars. These black hole mass measurements are mostly consistent with expectations based on the local M-sigma relationship, and are also consistent with single-epoch black hole mass measurements. This work increases the current sample size of reverberation-mapped active galaxies by about two-thirds and represents the first large sample of reverberation mapping observations beyond the local universe (z < 0.3).

Galaxy clusters trace the largest structures of the Universe and provide ideal laboratories for studying galaxy evolution and cosmology. Clusters with extended X-ray emission have been discovered at redshifts up to z ~ 2.5. Meanwhile, there has been growing interest in hunting for protoclusters, the progenitors of clusters, at higher redshifts. It is, however, very challenging to find the largest protoclusters at early times when they start to assemble. Here we report a giant protocluster of galaxies at redshift z = 5.7, when the Universe was only one billion years old. This protocluster occupies a volume of about 35x35x35 cubic co-moving megaparsecs. It is embedded in an even larger overdense region with at least 41 spectroscopically confirmed, luminous Lyman-alpha emitting galaxies (Lyman-alpha Emitters, or LAEs), including several previously reported LAEs. Its LAE density is 6.6 times the average density at z ~ 5.7. It is the only one of its kind in an LAE survey in four square degrees on the sky. Such a large structure is also rarely seen in current cosmological simulations. This protocluster will collapse into a galaxy cluster with a mass of (3.6+/-0.9) x 10^{15} solar masses, comparable to those of the most massive clusters or protoclusters known to date.

We report on the discovery of a dramatic X-ray spectral variability event observed in a $z\sim 1$ broad line type-1 QSO. The XMM-Newton spectrum from the year 2000 is characterized by an unobscured power-law spectrum with photon index of $\Gamma\sim 2$, a column density of $N_{\mathrm{H}}\sim 5\times 10^{20}\,\mathrm{cm^{-2}}$, and no prominent reflection component. Five years later, Chandra captured the source in a heavily-obscured, reflection-dominated state. The observed X-ray spectral variability could be caused by a Compton-thick cloud with $N_{\mathrm{H}}\sim 2\times 10^{24}\,\mathrm{cm^{-2}}$ eclipsing the direct emission of the hot corona, implying an extreme $N_{\mathrm{H}}$ variation never before observed in a type-1 QSO. An alternative scenario is a corona that switched off in between the observations. In addition, both explanations require a significant change of the X-ray luminosity prior to the obscuration or fading of the corona and/or a change of the relative geometry of the source/reflector system. Dramatic X-ray spectral variability of this kind could be quite common in type-1 QSOs, considering the relatively few datasets in which such an event could have been identified. Our analysis implies that there may be a population of type-1 QSOs which are Compton-thick in the X-rays when observed at any given time.

The Next-Generation Very Large Array (ngVLA) has the potential to be a workhorse for the discovery and study of paired supermassive black holes either at large separations (dual) or in tightly bound systems (binary). In this chapter, we outline the science case for the study of these supermassive pairs, and summarize discovery methods that can be used at radio wavelengths to discover them: including morphological, spectral, and time-domain identifications. One critical aspect of this work is that multi-messenger binary black hole studies may be possible with the ngVLA when combined with gravitational-wave searches using pulsar timing array techniques. However, long-baseline interferometery (>>1000 km) will make this possibility more likely by expanding the redshift range at which radio emission arising from two separate black holes may be resolved and studied.

We present the observations of an extreme ultraviolet (EUV) wave, a quasi-periodic fast-propagating (QFP) magnetosonic wave, and a kink wave that were simultaneously associated with the impingement of a coronal jet upon a group of coronal loops. After the interaction, the coronal loop showed obvious kink oscillation that had a period of about 428 seconds. In the meantime, a large-scale EUV wave and a QFP wave are observed on the west of the interaction position. It is interesting that the QFP wave showed refraction effect during its passing through two strong magnetic regions. The angular extent, speed, and lifetime of the EUV (QFP) wave were about 140 (40) degree, 423 (322) km/s, and 6 (26) minutes, respectively. It is measured that the period of the QFP wave was about 390 +/- 100 second. Based on the observational analysis results, we propose that the kink wave was probably excited by the interaction of the jet; the EUV was probably launched by the sudden expansion of the loop system due to the impingement of the coronal jet; and the QFP wave was possibly formed through the dispersive evolution of the disturbance caused by the jet-loop interaction.

We present X-ray and multiwavelength analyses of 32 weak emission-line quasars (WLQs) selected in a consistent and unbiased manner. New $Chandra$ 3.1-4.8 ks observations were obtained for 14 of these WLQs with C IV rest-frame equivalent widths (REWs) of 5-15 \r{A}, and these serve as an X-ray observational "bridge" between previously studied WLQs with C IV REW $\lesssim$ 5 \r{A} and more-typical quasars with C IV REW $\approx$ 15-100 \r{A}. We have identified and quantified a strong dependence of the fraction of X-ray weak quasars upon C IV REW; this fraction declines by a factor of $\approx 13$ (from $\approx 44$% to $\approx 3$%) for C IV REW ranging from 4-50 \r{A}, and the rate of decline appears particularly strong in the 10-20 \r{A} range. The dependence broadly supports the proposed "shielding" model for WLQs, in which a geometrically and optically thick inner accretion disk, expected for a quasar accreting at a high Eddington ratio, both prevents ionizing EUV/X-ray photons from reaching the high-ionization broad emission-line region and also sometimes blocks the line-of-sight to the central X-ray emitting region. This model is also supported by the hard average spectral shape of X-ray weak WLQs (with a power-law effective photon index of $\Gamma_{\rm eff}=1.19^{+0.56}_{-0.45}$). Additionally, we have examined UV continuum/emission-line properties that might trace X-ray weakness among WLQs, confirming that red UV continuum color is the most-effective tracer.

We present high-resolution observations of two kinds of dynamic behavior in a quiescent prominence using the New Vacuum Solar Telescope, i.e., Kelvin-Helmholtz instabilities (KHIs) and small-scale oscillations. The KHIs were identified as rapidly developed vortex-like structures with counter-clockwise/clockwise rotations in the Ha red-wing images at +0.3 A, which were produced by the strong shear-flows motions on the surface/interface of prominence plumes. The KHI growth rates are estimated to be about 0.0135 +(-)0.0004 and 0.0138 +(-) 0.0004. Our observational results further suggest that the shear velocities (i.e, supersonic) of the mass flows are fast enough to produce the strong deformation of the boundary and overcome the restraining surface tension force. This flow-driven instability might play a significant role in the process of plasma transfer in solar prominences. The small-scale oscillations perpendicular to the prominence threads are observed in the Ha line-center images. The oscillatory periods changed non-monotonically and showed two changing patterns, in which one firstly decreased slowly and then it changed to increase, while the other grew fast at the beginning and then it changed to decrease. Both of these two thread oscillations with changing periods were observed to be unstable for an entire cycle, and they were local in nature. All our findings indicate that the small-scale thread oscillations could be magnetohydrodynamic waves in the solar corona.

Kiloparsec-scale dual active galactic nuclei (AGNs) are active supermassive black hole pairs co-rotating in galaxies with separations of less than a few kpc. Expected to be a generic outcome of hierarchical galaxy formation, their frequency and demographics remain uncertain. We have carried out an imaging survey with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) of AGNs with double-peaked narrow [O III] emission lines. HST/WFC3 offers high image quality in the near-infrared (NIR) to resolve the two stellar nuclei, and in the optical to resolve [O III] from ionized gas in the narrow-line regions. This combination has proven to be key in sorting out alternative scenarios. With HST/WFC3 we are able to explore a new population of close dual AGNs at more advanced merger stages than can be probed from the ground. Here we show that the AGN SDSS J0924+0510, which had previously shown two stellar bulges, contains two spatially distinct [O III] regions consistent with a dual AGN. While we cannot completely exclude cross-ionization from a single central engine, the nearly equal ratios of [O III] strongly suggest a dual AGN with a projected angular separation of 0."4, corresponding to a projected physical separation of $r_p$ = 1 kpc at redshift z = 0.1495. This serves as a proof of principle for combining high-resolution NIR and optical imaging to identify close dual AGNs. Our result suggests that studies based on low-resolution and/or low-sensitivity observations may miss close dual AGNs and thereby may underestimate their occurrence rate on $\lesssim$ kpc scales.

Using the X-ray-selected active galactic nuclei (AGN) from the XMM-XXL north survey and the SDSS Baryon Oscillation Spectroscopic Survey (BOSS) spectroscopic follow-up of them, we compare the properties of X-ray unobscured and obscured broad-line AGN (BLAGN1 and BLAGN2; $N_\textrm{H}$below and above $10^{21.5}$ cm$^{-2}$), including their X-ray luminosity $L_X$, black hole mass, Eddington ratio $\lambda_{\textrm{Edd}}$, optical continuum and line features. We find that BLAGN2 have systematically larger broad line widths and hence apparently higher (lower) $M_{\textrm{BH}}$ ($\lambda_{\textrm{Edd}}$) than BLAGN1. We also find that the X-ray obscuration in BLAGN tends to coincide with optical dust extinction, which is optically thinner than that in narrow-line AGN (NLAGN) and likely partial-covering to the broad line region. All the results can be explained in the framework of a multi-component, clumpy torus model by interpreting BLAGN2 as an intermediate type between BLAGN1 and NLAGN in terms of an intermediate inclination angle.

We study the relation between the metallicity of quasar broad line region (BLR) and black hole (BH) mass ($10^{7.5}M_{\odot} \sim 10^{10}M_{\odot}$) and quasar bolometric luminosity ($10^{44.6}erg/s \sim 10^{48} erg/s$) using a sample of $\sim$130,000 quasars at $2.25\leq z\leq 5.25$ from Sloan Digital Sky Survey (SDSS) Data Release 12 (DR12). We generate composite spectra by stacking individual spectra in the same BH mass (bolometric luminosity) and redshift bins and then estimate the metallicity of quasar BLR using metallicity-sensitive broad emission-line flux ratios based on the photoionization models. We find a significant correlation between quasar BLR metallicity and BH mass (bolometric luminosity) but no correlation between quasar BLR metallicity and redshift. We also compare the metallicity of quasar BLR and that of host galaxies inferred from the mass-metallicity relation of star-forming galaxy at $z\sim2.3$ and $z\sim3.5$. We find quasar BLR metallicity is 0.3 $\sim$ 1.0 dex higher than their host galaxies. This discrepancy cannot be interpreted by the uncertainty due to different metallicity diagnostic methods, mass-metallicity relation of galaxy, metallicity gradient in quasar host galaxies, BH mass estimation, the effect of different spectral energy distribution (SED) models, and a few other potential sources of uncertainties. We propose a possibility that the high metallicity in quasar BLR might be caused by metal enrichment from massive star formation in the nucleus region of quasars or even the accretion disk.

We present the rest-frame UV and optical photometry and morphology of low-redshift broad-line quasar host galaxies from the Sloan Digital Sky Survey Reverberation Mapping project. Our sample consists of 103 quasars at $z<0.8$, spanning a luminosity range of $-25\le M_g\le -17$ mag. We stack the multi-epoch images in the $g$ and $i$ bands taken by the Canada-France-Hawaii Telescope. The combined $g$-band ($i$-band) images reach a $5\sigma$ depth of 26.2 (25.2) mag, with a typical PSF size of $0.7"$ ($0.6"$). Each quasar is decomposed into a PSF and a S\'ersic profile, representing the central AGN and the host galaxy components, respectively. The systematic errors of the measured host galaxy flux in the two bands are 0.23 and 0.18 mag. The relative errors of the measured galaxy half-light radii ($R_e$) are about 13%. We estimate the rest-frame $u$- and $g$-band flux of the host galaxies, and find that the AGN-to-galaxy flux ratios in the $g$ band are between 0.9 to 4.4 (68.3% confidence). These galaxies have high stellar masses $M_\ast = 10^{10}\sim10^{11}\, M_\odot$. They have similar color with star-forming galaxies at similar redshifts, in consistent with AGN positive feedback in these quasars. We find that the $M_*-M_\text{BH}$ relation in our sample is shallower than the local $M_\text{Bulge}-M_\text{BH}$ relation. The S\'ersic indices and the $M_*-R_e$ relation indicate that the majority of the host galaxies are disk-like.

We present two events to study the driving mechanism of extreme-ultraviolet (EUV) waves that are not associated with coronal mass ejections (CMEs), by using high resolution observations taken by the Atmospheric Imaging Assembly (AIA) on board Solar Dynamics Observatory. Observational results indicate that the observed EUV waves were accompanied by ares and coronal jets, but without CMEs that were regarded as the driver of most EUV waves in previous studies. In the first case, it is observed that a coronal jet ejected along a transequatorial loop system at a plane-of-the-sky (POS) speed of $335 \pm 22$ km/s, in the meantime, an arc-shaped EUV wave appeared on the eastern side of the loop system. In addition, the EUV wave further interacted with another interconnecting loop system and launched a fast propagating (QFP) magnetosonic wave along the loop system, which had a period of 200 s and a speed of $388 \pm 65$ km/s, respectively. In the second case, we also observed a coronal jet ejected at a POS speed of $282 \pm 44$ km/s along a transequatorial loop system and the generation of bright EUV wave on the eastern side of the loop system. Based on the observational results, we propose that the observed EUV waves on the eastern side of the transequatorial loop systems are fast-mode magnetosonic waves, and they were driven by the sudden lateral expansion of the transequatorial loop systems due to the direct impingement of the associated coronal jets, while the QFP wave in the first case formed due to the dispersive evolution of the disturbance caused by the interaction between the EUV wave and the interconnecting coronal loops. It is noted that EUV waves driven by sudden loop expansions have shorter lifetimes than those driven by CMEs.

The detailed observational analysis of a homologous Extreme-ultraviolet (EUV) wave event is presented to study the driving mechanism and the physical property of the EUV waves, combining high resolution data taken by the Solar Dynamics Observatory and the Solar TErrestrial RElations Observatory. It is observed that four homologous EUV waves originated from the same active region AR11476 within about one hour, and the time separations between consecutive waves were of 8 - 20 minutes. The waves showed narrow arc-shaped wavefronts and propagated in the same direction along a large-scale transequatorial loop system at a speed of 648 - 712 km/s and a deceleration of 0.985 - 1.219 km/s2. The EUV waves were accompanied by weak flares, coronal jets, and radio type III bursts, in which the EUV waves were delayed with respect to the start times of the radio type III bursts and coronal jets about 2 - 13 and 4 - 9 minutes, respectively. Different to previous studies of homologous EUV waves, no coronal mass ejections were found in the present event. Based on the observational results and the close temporal the spatial relationship between the EUV waves and the coronal jets, for the first time, we propose that the observed homologous EUV waves were large-amplitude nonlinear fast-mode magnetosonic waves or shocks driven by the associated recurrent coronal jets, resemble the generation mechanism of a piston shock in a tube. In addition, it is found that the recurrent jets were tightly associated with the alternating flux cancellation and emergence in the eruption source region and radio type III bursts.

We present the observations of a double-decker filament to study its formation, triggering, and eruption physics. It is observed that the double-decker filament was formed by splitting of an original single filament. During the splitting process, intermittent bright point bursts are observed in the filament channel, which resulted in the generation of the upper filament branch. The eruption of the newly formed double-decker filament was possibly triggered by two recurrent two-sided loop jets in the filament channel and the continuous mass unloading from the upper filament body. The interaction between the first jet and the filament directly resulted in the unstable of the lower branch and the fast rising phase of the upper branch. The second jet occurred at the same site about three hours after the first one, which further disturbed and accelerated the rising of the lower filament branch. It is interesting that the rising lower branch overtook the upper one, and then the two branches probably merged into one filament. Finally, the whole filament erupted violently and caused a large-scale coronal mass ejection, leaving behind a pair of flare ribbons and two dimming regions on the both sides of the filament channel. We think that the intermittent bursts may directly result in the rearrangement of the filament magnetic field and therefore the formation of the double-decker filament, then the recurrent jets further caused the fully eruption of the entire filament system. The study provides convincing evidence for supporting the scenario that a double-decker filament can be formed by splitting a single filament into two branches.

A small fraction of quasars show an unusually high nitrogen-to-carbon ratio (N/C) in their spectra. These "nitrogen-rich" (N-rich) quasars are a long-standing puzzle because their interstellar medium implies stellar populations with abnormally high metallicities. It has recently been proposed that N-rich quasars may result from tidal disruption events (TDEs) of stars by supermassive black holes. The rapid enhancement of nitrogen and the depletion of carbon due to the carbon--nitrogen--oxygen cycle in supersolar mass stars could naturally produce high N/C. However, the TDE hypothesis predicts that the N/C should change with time, which has never hitherto been observed. Here we report the discovery of the first N-rich quasar with rapid N/C variability that could be caused by a TDE. Two spectra separated by 1.7 years (rest-frame) show that the N III]\lambda 1750/C III]\lambda 1909 intensity ratio decayed by ~86%+/-14% (1 \sigma). Optical (rest-frame UV) light-curve and X-ray observations are qualitatively consistent with the TDE hypothesis; though, the time baseline falls short of a definitive proof. Putting the single-object discovery into context, statistical analyses of the ~80 known N-rich quasars with high-quality archival spectra show evidence (at a 5\sigma\ significance level) of a decrease in N/C on timescales of >1 year (rest-frame) and a constant level of ionization (indicated by the C III]\lambda 1909/C IV\lambda 1549 intensity ratio). If confirmed, our results demonstrate the method of identifying TDE candidates in quasars via abundance ratio variability, opening a new window of TDE observations at high redshift (z>2) with upcoming large-scale time-domain spectroscopic surveys.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

A stationary and spherically symmetric black hole (For example, Reissner-Nordstrom black hole or Kerr-Newman black hole) has at most one singularity and two horizons. One horizon is the outer event horizon and the other is the inner Cauchy horizon. Can we construct static and spherically symmetric black hole solutions with N horizons and M singularities? De Sitter cosmos has only one apparent horizon. Can we construct cosmos solutions with N horizons? In this article, we present the static and spherically symmetric black hole and cosmos solutions with N horizons and M singularities in the vector-tensor theories. Following these motivations, we also construct the black hole solutions with a firewall. The deviation of these black hole solutions from the usual ones can be potentially tested by future measurements of gravitational waves.

The mean white dwarf (WD) mass in the Galactic bulge cataclysmic variables (CVs) was measured by applying the the shock temperature-WD mass correlation of magnetic cataclysmic variables (mCVs) to the Galactic Bulge X-ray Emission (GBXE) spectra. However, the resulting mean WD mass is lower than that of the local CVs. This discrepancy could be explained by the dominating sources in the GBXE are non-mCVs instead of mCVs. In this work, we conduct an thorough investigation on the X-ray spectra of local DNe from \suzaku\ archives, and derive semi-empirical correlations between the shock temperature $T_{\rm max}$, the flux ratio of \fec\ to \feb\ lines, and WD mass for quiescent, non-magnetic CVs. By applying these correlations to the GBXE, we derive the average WD mass of CVs in the Galactic bulge to be $0.81\pm 0.07M_\odot$. This value is consistent with previous optical measurements of WD mass in local CVs.

The observational analysis is performed to study the excitation mechanism and the propagation properties of a quasi-periodic fast-propagating (QFP) magnetosonic wave. The QFP wave was associated with the eruption of a nearby mini-filament and a small B4 GOES flare, which may indicate that the generation of a QFP wave do not need too much flare energy. The propagation of the QFP wave was along a bundle of funnel-shaped open loops with a speed of about 1100+/-78, and an acceleration of -2.2+/-1.1. Periodicity analysis indicates that the periods of the QFP wave are 43+/-6, 79+/-18 second. For the first time, we find that the periods of the QFP wave and the accompanying flare are inconsistent, which is different from the findings as reported in previous studies. We propose that the present QFP wave was possibly caused by the mechanism of dispersive evolution of an initially broadband disturbance resulted from the nearby mini-filament eruption.

This paper presents Very Long Baseline Array (VLBA) observations of 13 double-peaked [O III] emission-line type-2 Active Galactic Nuclei (AGNs) at redshifts 0.06 < z < 0.41 (with a median redshift of z~0.15) identified in the Sloan Digital Sky Survey. Such double-peaked emission-line objects may result from jets or outflows from the central engine or from a dual AGN. The VLBA provides an angular resolution of <~10 pc at the distance of many of these galaxies, sufficient to resolve the radio emission from extremely close dual AGNs and to contribute to understanding the origin of double-peaked [O III] emission lines. Of the 13 galaxies observed at 3.6 cm (8.4 GHz), we detect six at a 1\sigma\ sensitivity level of ~0.15 mJy/beam, two of which show clear jet structures on scales ranging from a few milliarcseconds to tens of milliarcseconds (corresponding to a few pc to tens of pc at a median redshift of 0.15). We suggest that radio-loud double-peaked emission-line type-2 AGNs may be indicative of jet produced structures, but a larger sample of double-peaked [O III] AGNs with high angular resolution radio observations will be required to confirm this suggestion.

We use the multi-epoch spectra of 362 quasars from the Sloan Digital Sky Survey Reverberation Mapping project to investigate the dependence of the blueshift of CIV relative to MgII on quasar properties. We confirm that high-blueshift sources tend to have low CIV equivalent widths (EWs), and that the low-EW sources span a range of blueshift. Other high-ionization lines, such as HeII, also show similar blueshift properties. The ratio of the line width (measured as both the full-width at half maximum and the velocity dispersion) of CIV to that of MgII increases with blueshift. Quasar variability might enhance the connection between the CIV blueshift and quasar properties (e.g., EW). The variability of the MgII line center (i.e., the wavelength that bisects the cumulative line flux) increases with blueshift. In contrast, the CIV line center shows weaker variability at the extreme blueshifts. Quasars with the high-blueshift CIV lines tend to have less variable continuum emission, when controlling for EW, luminosity, and redshift. Our results support the scenario that high-blueshift sources tend to have large Eddington ratios.

We present the first rest-frame UV population study of 17 heavily reddened, high-luminosity (E(B-V)$_{\rm{QSO}}\gtrsim$ 0.5; L$_{\rm{bol}}>$ 10$^{46}$ergs$^{-1}$) broad-line quasars at $1.5 < z < 2.7$. We combine the first year of deep, optical, ground-based observations from the Dark Energy Survey (DES) with the near infrared VISTA Hemisphere Survey (VHS) and UKIDSS Large Area Survey (ULAS) data, from which the reddened quasars were initially identified. We demonstrate that the significant dust reddening towards the quasar in our sample allows host galaxy emission to be detected at the rest-frame UV wavelengths probed by the DES photometry. By exploiting this reddening effect, we disentangle the quasar emission from that of the host galaxy via spectral energy distribution (SED) fitting. We find evidence for a relatively unobscured, star-forming host galaxy in at least ten quasars, with a further three quasars exhibiting emission consistent with either star formation or scattered light. From the rest-frame UV emission, we derive instantaneous, dust-corrected star formation rates (SFRs) in the range 25 < SFR$_{\rm{UV}}$ < 365 M$_{\odot}$yr$^{-1}$, with an average SFR$_{\rm{UV}}$ = 130 $\pm$ 95 M$_{\odot}$yr$^{-1}$. We find a broad correlation between SFR$_{\rm{UV}}$ and the bolometric quasar luminosity. Overall, our results show evidence for coeval star formation and black hole accretion occurring in luminous, reddened quasars at the peak epoch of galaxy formation.

Using high temporal and high spatial resolution observations taken by the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory, we present the detailed observational analysis of a high quality quasi-periodic fast- propagating (QFP) magnetosonic wave that was associated with the eruption of a magnetic flux rope and a GOES C5.0 flare. For the first time, we find that the QFP wave lasted during the entire flare lifetime rather than only the rising phase of the accompanying flare as reported in previous studies. In addition, the propagation of the different parts of the wave train showed different kinematics and morphologies. For the southern (northern) part, the speed, duration, intensity variation are about 875 +/- 29 (1485 +/- 233) km/s, 45 (60) minutes, and 4% (2%), and the pronounced periods of them are 106 +/- 12 and 160 +/- 18 (75 +/- 10 and 120 +/- 16) seconds, respectively. It is interesting that the northern part of the wave train showed obvious refraction effect when they pass through a region of strong magnetic field. Periodicity analysis result indicates that all the periods of the QFP wave can be found in the period spectrum of the accompanying flare, suggesting their common physical origin. We propose that the quasi-periodic nonlinear magnetohydrodynamics process in the magnetic reconnection that produces the accompanying flare should be important for exciting of QFP wave, and the different magnetic distribution along different paths can account for the different speeds and morphology evolution of the wave fronts.

The Hawking-Penrose singularity theorem states that a singularity forms inside a black hole in general relativity. To remove this singularity one must resort to a more fundamental theory. Using a corrected dynamical equation arising in loop quantum cosmology and braneworld models, we study the gravitational collapse of a perfect fluid sphere with a rather general equation of state. In the frame of an observer comoving with this fluid, the sphere pulsates between a maximum and a minimum size, avoiding the singularity. The exterior geometry is also constructed. There are usually an outer and an inner apparent horizon, resembling the Reissner-Nordstr\"om situation. For a distant observer the {horizon} crossing occurs in an infinite time and the pulsations of the black hole quantum "beating heart" are completely unobservable. However, it may be observable if the black hole is not spherical symmetric and radiates gravitational wave due to the quadrupole moment, if any.

We present composite broad-line region (BLR) reverberation-mapping lag measurements for \halpha, \hbeta, \HeII\,$\lambda4686$ and \MgII\ for a sample of 144, $z\lesssim 1$ quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Using only the 32-epoch spectroscopic light curves in the first 6-month season of SDSS-RM observations, we compile correlation-function measurements for individual objects and then coadd them to allow the measurement of the average lags for our sample at mean redshifts of $0.4$ (for \halpha) and $\sim 0.65$ (for the other lines). At similar quasar luminosities and redshifts, the sample-averaged lag decreases in the order of \MgII, \halpha, \hbeta\ and \HeII. This decrease in lags is accompanied by an increase in the mean line width of the four lines, and is roughly consistent with the virialized motion for BLR gas in photoionization equilibrium. These are among the first RM measurements of stratified BLR structure at $z>0.3$. Dividing our sample by luminosity, \halpha\ shows clear evidence of increasing lags with luminosity, consistent with the expectation from the measured BLR size-luminosity relation based on \hbeta. The other three lines do not show a clear luminosity trend in their average lags due to the limited dynamic range of luminosity probed and the poor average correlation signals in the divided samples, a situation that will be improved with the incorporation of additional photometric and spectroscopic data from SDSS-RM. We discuss the utility and caveats of composite-lag measurements for large statistical quasar samples with reverberation-mapping data.

We present observations of a small-scale Extreme-ultraviolet (EUV) wave that was associated with a mini-filament eruption and a GOES B1.9 micro-fare in the quiet Sun region. The initiation of the event was due to the photospheric magnetic emergence and cancellation in the eruption source region, which first caused the ejection of a small plasma ejecta, then the ejecta impacted on a nearby mini-filament and thereby led to the filament's eruption and the associated fare. During the filament eruption, an EUV wave at a speed of 182 { 317 km/s was formed ahead of an expanding coronal loop, which propagated faster than the expanding loop and showed obvious deceleration and refection during the propagation. In addition, the EUV wave further resulted in the transverse oscillation of a remote filament whose period and damping time are 15 and 60 minutes, respectively. Based on the observational results, we propose that the small-scale EUV wave should be a fast-mode magnetosonic wave that was driven by the the expanding coronal loop. Moreover, with the application of filament seismology, it is estimated that the radial magnetic field strength is about 7 Gauss. The observations also suggest that small-scale EUV waves associated with miniature solar eruptions share similar driving mechanism and observational characteristics with their large-scale counterparts.

As astronomers increasingly exploit the information available in the time domain, spectroscopic variability in particular opens broad new channels of investigation. Here we describe the selection algorithms for all targets intended for repeat spectroscopy in the Time Domain Spectroscopic Survey (TDSS), part of the extended Baryon Oscillation Spectroscopic Survey within the Sloan Digital Sky Survey-IV. Also discussed are the scientific rationale and technical constraints leading to these target selections. The TDSS includes a large "Repeat Quasar Spectroscopy" (RQS) program delivering ~13,000 repeat spectra of confirmed SDSS quasars, and several smaller "Few-Epoch Spectroscopy" (FES) programs targeting specific classes of quasars as well as stars. The RQS program aims to provide a large and diverse quasar data set for studying variations in quasar spectra on timescales of years, a comparison sample for the FES quasar programs, and opportunity for discovering rare, serendipitous events. The FES programs cover a wide variety of phenomena in both quasars and stars. Quasar FES programs target broad absorption line quasars, high signal-to-noise ratio normal broad line quasars, quasars with double-peaked or very asymmetric broad emission line profiles, binary supermassive black hole candidates, and the most photometrically variable quasars. Strongly variable stars are also targeted for repeat spectroscopy, encompassing many types of eclipsing binary systems, and classical pulsators like RR Lyrae. Other stellar FES programs allow spectroscopic variability studies of active ultracool dwarf stars, dwarf carbon stars, and white dwarf/M dwarf spectroscopic binaries. We present example TDSS spectra and describe anticipated sample sizes and results.

SDSS-V will be an all-sky, multi-epoch spectroscopic survey of over six million objects. It is designed to decode the history of the Milky Way, trace the emergence of the chemical elements, reveal the inner workings of stars, and investigate the origin of planets. It will also create an integral-field spectroscopic map of the gas in the Galaxy and the Local Group that is 1,000x larger than the current state of the art and at high enough spatial resolution to reveal the self-regulation mechanisms of galactic ecosystems. SDSS-V will pioneer systematic, spectroscopic monitoring across the whole sky, revealing changes on timescales from 20 minutes to 20 years. The survey will thus track the flickers, flares, and radical transformations of the most luminous persistent objects in the universe: massive black holes growing at the centers of galaxies. The scope and flexibility of SDSS-V will be unique among extant and future spectroscopic surveys: it is all-sky, with matched survey infrastructures in both hemispheres; it provides near-IR and optical multi-object fiber spectroscopy that is rapidly reconfigurable to serve high target densities, targets of opportunity, and time-domain monitoring; and it provides optical, ultra-wide-field integral field spectroscopy. SDSS-V, with its programs anticipated to start in 2020, will be well-timed to multiply the scientific output from major space missions (e.g., TESS, Gaia, eROSITA) and ground-based projects. SDSS-V builds on the 25-year heritage of SDSS's advances in data analysis, collaboration infrastructure, and product deliverables. The project is now refining its science scope, optimizing the survey strategies, and developing new hardware that builds on the SDSS-IV infrastructure. We present here an overview of the current state of these developments as we seek to build our worldwide consortium of institutional and individual members.

We present the observations of a blowout jet that experienced two distinct ejection stages. The first stage started from the emergence of a small positive magnetic polarity, which cancelled with the nearby negative magnetic field and caused the rising of a mini-filament and its confining loops. This further resulted in a small jet due to the magnetic reconnection between the rising confining loops and the overlying open field. The second ejection stage was mainly due to the successive removal of the confining field by the reconnection. Thus that the filament erupted and the erupting cool filament material directly combined with the hot jet originated form the reconnection region and therefore formed the cool and hot components of the blowout jet. During the two ejection stages, cool H{\alpha} jets are also observed cospatial with their coronal counterparts, but their appearance times are earlier than the hot coronal jets a few minutes. Therefore, the hot coronal jets are possibly caused by the heating of the cool H{\alpha} jets, or the rising of the reconnection height from chromosphere to the corona. The scenario that magnetic reconnection occurred between the confining loops and the overlying open loops are supported by many observational facts, including the bright patches on the both sides of the mini-filament, hot plasma blobs along the jet body, and periodic metric radio type III bursts at the very beginnings of the two stages. The evolution and characteristics of these features manifest the detailed non-linear process in the magnetic reconnection.

We present a spectroscopic survey of high-redshift, luminous galaxies over four square degrees on the sky, aiming to build a large and homogeneous sample of Ly$\alpha$ emitters (LAEs) at $z\approx5.7$ and 6.5, and Lyman-break galaxies (LBGs) at $5.5<z<6.8$. The fields that we choose to observe are well-studied, such as SXDS and COSMOS. They have deep optical imaging data in a series of broad and narrow bands, allowing efficient selection of galaxy candidates. Spectroscopic observations are being carried out using the multi-object spectrograph M2FS on the Magellan Clay telescope. M2FS is efficient to identify high-redshift galaxies, owing to its 256 optical fibers deployed over a circular field-of-view 30 arcmin in diameter. We have observed $\sim2.5$ square degrees. When the program is completed, we expect to identify more than 400 bright LAEs at $z\approx5.7$ and 6.5, and a substantial number of LBGs at $z\ge6$. This unique sample will be used to study a variety of galaxy properties and to search for large protoclusters. Furthermore, the statistical properties of these galaxies will be used to probe cosmic reionization. We describe the motivation, program design, target selection, and M2FS observations. We also outline our science goals, and present a sample of the brightest LAEs at $z\approx5.7$ and 6.5. This sample contains 32 LAEs with Ly$\alpha$ luminosities higher than 10$^{43}$ erg s$^{-1}$. A few of them reach $\ge3\times10^{43}$ erg s$^{-1}$, comparable to the two most luminous LAEs known at $z\ge6$, `CR7' and `COLA1'. These LAEs provide ideal targets to study extreme galaxies in the distant universe.

We present observational analysis of two successive two-sided loop jets observed by the ground-based New Vacuum Solar Telescope (NVST) and the space-borne Solar Dynamics Observatory ( SDO). The two successive two-sided loop jets manifested similar evolution process and both were associated with the interaction of two small-scale adjacent filamentary threads, magnetic emerging and cancellation processes at the jet's source region. High temporal and high spatial resolution observations reveal that the two adjacent ends of the two filamentary threads are rooted in opposite magnetic polarities within the source region. The two threads approached to each other, and then an obvious brightening patch is observed at the interaction position. Subsequently, a pair of hot plasma ejections are observed heading to opposite directions along the paths of the two filamentary threads, and with a typical speed of two-sided loop jets of the order 150 km/s. Close to the end of the second jet, we report the formation of a bright hot loop structure at the source region, which suggests the formation of new loops during the interaction. Based on the observational results, we propose that the observed two-sided loop jets are caused by the magnetic reconnection between the two adjacent filamentary threads, largely different from the previous scenario that a two-sided loop jet is generated by magnetic reconnection between an emerging bipole and the overlying horizontal magnetic fields.

PHL 6625 is a luminous quasi-stellar object (QSO) at z = 0.3954 located behind the nearby galaxy NGC 247 (z = 0.0005). Hubble Space Telescope (HST) observations revealed an arc structure associated with it. We report on spectroscopic observations with the Very Large Telescope (VLT) and multiwavelength observations from the radio to the X-ray band for the system, suggesting that PHL 6625 and the arc are a close pair of merging galaxies, instead of a strong gravitational lens system. The QSO host galaxy is estimated to be (4-28) x 10^10 M_sun, and the mass of the companion galaxy of is estimated to be M_* = (6.8 +/- 2.4) x 10^9 M_sun, suggesting that this is a minor merger system. The QSO displays typical broad emission lines, from which a black hole mass of about (2-5) x 10^8 M_sun and an Eddington ratio of about 0.01-0.05 can be inferred. The system represents an interesting and rare case where a QSO is associated with an ongoing minor merger, analogous to Arp 142.

We perform a systematic search for long-term extreme variability quasars (EVQs) in the overlapping Sloan Digital Sky Survey (SDSS) and 3-Year Dark Energy Survey (DES) imaging, which provide light curves spanning more than 15 years. We identified ~1000 EVQs with a maximum g band magnitude change of more than 1 mag over this period, about 10% of all quasars searched. The EVQs have L_bol~10^45-10^47 erg/s and L/L_Edd~0.01-1. Accounting for selection effects, we estimate an intrinsic EVQ fraction of ~30-50% among all g<~22 quasars over a baseline of ~15 years. These EVQs are good candidates for so-called "changing-look quasars", where a spectral transition between the two types of quasars (broad-line and narrow-line) is observed between the dim and bright states. We performed detailed multi-wavelength, spectral and variability analyses for the EVQs and compared to their parent quasar sample. We found that EVQs are distinct from a control sample of quasars matched in redshift and optical luminosity: (1) their UV broad emission lines have larger equivalent widths; (2) their Eddington ratios are systematically lower; and (3) they are more variable on all timescales. The intrinsic difference in quasar properties for EVQs suggest that internal processes associated with accretion are the main driver for the observed extreme long-term variability. However, despite their different properties, EVQs seem to be in the tail of a continuous distribution of quasar properties, rather than standing out as a distinct population. We speculate that EVQs are normal quasars accreting at relatively low accretion rates, where the accretion flow is more likely to experience instabilities that drive the factor of few changes in flux on multi-year timescales.

The wide-area XMM-XXL X-ray survey is used to explore the fraction of obscured AGN at high accretion luminosities, $L_X (\rm 2-10 \, keV) > 10^{44} \, erg \,s ^{-1}$, and out to redshift $z\approx1.5$. The sample covers an area of about $\rm14\,deg^2$ and provides constraints on the space density of powerful AGN over a wide range of neutral hydrogen column densities extending beyond the Compton-thick limit, $\rm N_H\approx10^{24}\,cm^{-2}$. The fraction of obscured Compton-thin ($\rm N_H=10^{22}-10^{24}\,cm^{-2}$) AGN is estimated to be $\approx0.35$ for luminosities $L_X(\rm 2-10\,keV)>10^{44}\,erg\,s^{-1}$ independent of redshift. For less luminous sources the fraction of obscured Compton-thin AGN increases from $0.45\pm0.10$ at $z=0.25$ to $0.75\pm0.05$ at $z=1.25$. Studies that select AGN in the infrared via template fits to the observed Spectral Energy Distribution of extragalactic sources estimate space densities at high accretion luminosities consistent with the XMM-XXL constraints. There is no evidence for a large population of AGN (e.g. heavily obscured) identified in the infrared and missed at X-ray wavelengths. We further explore the mid-infrared colours of XMM-XXL AGN as a function of accretion luminosity, column density and redshift. The fraction of XMM-XXL sources that lie within the mid-infrared colour wedges defined in the literature to select AGN is primarily a function of redshift. This fraction increases from about 20-30% at z=0.25 to about 50-70% at $z=1.5$.

We present $Chandra$ observations of seven broad absorption line (BAL) quasars at $z=0.863$-2.516 with redshifted BAL troughs (RSBALs). Five of our seven targets were detected by $Chandra$ in 4-13 ks exposures with ACIS-S. The $\alpha_{\rm ox}$ values, $\Delta\alpha_{\rm ox}$ values, and spectral energy distributions of our targets demonstrate they are all X-ray weak relative to expectations for non-BAL quasars, and the degree of X-ray weakness is consistent with that of appropriately-matched BAL quasars generally. Furthermore, our five detected targets show evidence for hard X-ray spectral shapes with a stacked effective power-law photon index of $\Gamma_{\rm eff}=0.5^{+0.5}_{-0.4}$. These findings support the presence of heavy X-ray absorption ($N_{\rm H}\approx 2 \times 10^{23}$ cm$^{-2}$) in RSBAL quasars, likely by the shielding gas found to be common in BAL quasars more generally. We use these X-ray measurements to assess models for the nature of RSBAL quasars, finding that a rotationally-dominated outflow model is favored while an infall model also remains plausible with some stipulations. The X-ray data disfavor a binary quasar model for RSBAL quasars in general.

To investigate AGN outflows as a tracer of AGN feedback on star-formation, we perform integral-field spectroscopy of 20 type 2 AGNs at z<0.1, which are luminous AGNs with the [O III] luminosity >10$^{41.5}$ erg/s, and exhibit strong outflow signatures in the [O III] kinematics. By decomposing the emission-line profile, we obtain the maps of the narrow and broad components of [O III] and H$\alpha$ lines, respectively. The broad components in both [O III] and H$\alpha$ represent the non-gravitational kinematics, i.e., gas outflows, while the narrow components, especially in H$\alpha$, represent the gravitational kinematics, i.e., rotational disk. By using the integrated spectra within the flux-weighted size of the narrow-line region, we estimate the energetics of the gas outflows. The ionized gas mass is 1.0-38.5$\times 10^5$ $M_{\odot}$, and the mean mass outflow rate is 4.6$\pm$4.3 $M_{\odot}$/yr, which is a factor of ~260 higher than the mean mass accretion rate 0.02$\pm$0.01 $M_{\odot}$/yr. The mean energy injection rate of the sample is 0.8$\pm$0.6% of the AGN bolometric luminosity, while the momentum flux is (5.4$\pm$3.6)$\times$ $L_{bol}$/c on average, except for two most kinematically energetic AGNs with low $L_{bol}$, which are possibly due to the dynamical timescale of the outflows. The estimated outflow energetics are consistent with the theoretical expectations for energy-conserving outflows from AGNs, yet we find no supporting evidence of instantaneous quenching of star formation due to the outflows.

We study the time lags between the continuum emission of quasars at different wavelengths, based on more than four years of multi-band ($g$, $r$, $i$, $z$) light-curves in the Pan-STARRS Medium Deep Fields. As photons from different bands emerge from different radial ranges in the accretion disk, the lags constrain the sizes of the accretion disks. We select 240 quasars with redshifts $z \approx 1$ or $z \approx 0.3$ that are relatively emission line free. The light curves are sampled from day to month timescales, which makes it possible to detect lags on the scale of the light crossing time of the accretion disks. With the code JAVELIN, we detect typical lags of several days in the rest frame between the $g$ band and the $riz$ bands. The detected lags are $\sim 2-3$ times larger than the light crossing time estimated from the standard thin disk model, consistent with the recently measured lag in NGC5548 and micro-lensing measurements of quasars. The lags in our sample are found to increase with increasing luminosity. Furthermore, the increase in lags going from $g-r$ to $g-i$ and then to $g-z$ is slower than predicted in the thin disk model, particularly for high luminosity quasars. The radial temperature profile in the disk must be different from what is assumed. We also find evidence that the lags decrease with increasing line ratios between ultraviolet FeII lines and MgII, which may point to changes in the accretion disk structure at higher metallicity.

Quasar emission lines are often shifted from the systemic velocity due to various dynamical and radiative processes in the line-emitting region. The level of these velocity shifts depends both on the line species and on quasar properties. We study velocity shifts for the line peaks of various narrow and broad quasar emission lines relative to systemic using a sample of 849 quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. The coadded (from 32 epochs) spectra of individual quasars have sufficient signal-to-noise ratio (SNR) to measure stellar absorption lines to provide reliable systemic velocity estimates, as well as weak narrow emission lines. The sample also covers a large dynamic range in quasar luminosity (~2 dex), allowing us to explore potential luminosity dependence of the velocity shifts. We derive average line peak velocity shifts as a function of quasar luminosity for different lines, and quantify their intrinsic scatter. We further quantify how well the peak velocity can be measured for various lines as a function of continuum SNR, and demonstrate there is no systematic bias in the line peak measurements when the spectral quality is degraded to as low as SNR~3 per SDSS pixel. Based on the observed line shifts, we provide empirical guidelines on redshift estimation from [OII]3728, [OIII]5008, [NeV]3426, MgII, CIII], HeII1640, broad Hbeta, CIV, and SiIV, which are calibrated to provide unbiased systemic redshifts in the mean, but with increasing intrinsic uncertainties of 46, 56, 119, 205, 233, 242, 400, 415, and 477 km/s, in addition to the measurement uncertainties. These more realistic redshift uncertainties are generally much larger than the formal uncertainties reported by the redshift pipelines for spectroscopic quasar surveys, and demonstrate the infeasibility of measuring quasar redshifts to better than ~200 km/s with only broad lines.

In this paper, we explore the quantum spacetimes that are potentially connected with the generalized uncertainty principles. By analyzing the gravity-induced quantum interference pattern and the Gedanken for weighting photon, we find that the generalized uncertainty principles inspire the effective Newton constant as same as our previous proposal. A characteristic momentum associated with the tidal effect is suggested, which incorporates the quantum effect with the geometric nature of gravity. When the simplest generalized uncertainty principle is considered, the minimal model of the regular black holes is reproduced by the effective Newton constant. The black hole's tunneling probability, accurate to the second order correction, is carefully analyzed. We find that the tunneling probability is regularized by the size of the black hole remnant. Moreover, the black hole remnant is the final state of a tunneling process that the probability is minimized. A theory of modified gravity is suggested, by substituting the effective Newton constant into the Hilbert-Einstein action.

We report on the co-existence of longitudinal and transverse oscillations in an active filament. On March 15$^{th}$ 2013, a M1.1 class flare was observed in the active region AR 11692. A CME was found to be associated with the flare. {The CME generated a shock wave that triggered the oscillations in a nearby filament}, situated at the south-west of the active region as observed from National Solar Observatory (NSO)\textit{Global Oscillation Network Group}(GONG) H$\alpha$ images. In this work we report the longitudinal oscillations in the two ends of the filament, co-existing with the transverse oscillations. We propose a scenario in which {an} incoming shock wave hits the filament obliquely and triggers both longitudinal and transverse oscillations. Using the observed parameters, we estimate the lower limit of the magnetic field strength. We use simple pendulum model with gravity as the restoring force to estimate the radius of curvature. We also calculate the mass accretion rate which causes the filament motions to damp quite fast.

We present the discovery of nine quasars at $z\sim6$ identified in the Sloan Digital Sky Survey (SDSS) imaging data. This completes our survey of $z\sim6$ quasars in the SDSS footprint. Our final sample consists of 52 quasars at $5.7<z\le6.4$, including 29 quasars with $z_{\rm AB}\le20$ mag selected from 11,240 deg$^2$ of the SDSS single-epoch imaging survey (the main survey), 10 quasars with $20\le z_{\rm AB}\le20.5$ selected from 4223 deg$^2$ of the SDSS overlap regions (regions with two or more imaging scans), and 13 quasars down to $z_{\rm AB}\approx22$ mag from the 277 deg$^2$ in Stripe 82. They span a wide luminosity range of $-29.0\le M_{1450}\le-24.5$. This well-defined sample is used to derive the quasar luminosity function (QLF) at $z\sim6$. After combining our SDSS sample with two faint ($M_{1450}\ge-23$ mag) quasars from the literature, we obtain the parameters for a double power-law fit to the QLF. The bright-end slope $\beta$ of the QLF is well constrained to be $\beta=-2.8\pm0.2$. Due to the small number of low-luminosity quasars, the faint-end slope $\alpha$ and the characteristic magnitude $M_{1450}^{\ast}$ are less well constrained, with $\alpha=-1.90_{-0.44}^{+0.58}$ and $M^{\ast}=-25.2_{-3.8}^{+1.2}$ mag. The spatial density of luminous quasars, parametrized as $\rho(M_{1450}<-26,z)=\rho(z=6)\,10^{k(z-6)}$, drops rapidly from $z\sim5$ to 6, with $k=-0.72\pm0.11$. Based on our fitted QLF and assuming an IGM clumping factor of $C=3$, we find that the observed quasar population cannot provide enough photons to ionize the $z\sim6$ IGM at $\sim90$\% confidence. Quasars may still provide a significant fraction of the required photons, although much larger samples of faint quasars are needed for more stringent constraints on the quasar contribution to reionization.

We find a method to rewrite the equations of motion of scalar fields, generalized DBI field and quintessence, in the autonomous form for\emph{arbitrary} scalar potentials. With the aid of this method, we explore the cosmic evolution of generalized DBI field and quintessence with the potential of multiple vacua. Then we find that the scalars are always frozen in the false or true vacuum in the end. Compared to the evolution of quintessence, the generalized DBI field has more times of oscillations around the vacuum of the potential. The reason for this point is that, with the increasing of speed $\dot{\phi}$, the friction term of generalized DBI field is greatly decreased. Thus the generalized DBI field acquires more times of oscillations.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

The uncertain origin of the recently-discovered `changing-looking' quasar phenomenon -- in which a luminous quasar dims significantly to a quiescent state in repeat spectroscopy over ~10 year timescales -- may present unexpected challenges to our understanding of quasar accretion. To better understand this phenomenon, we take a first step to building a sample of changing-look quasars with a systematic but simple archival search for these objects in the Sloan Digital Sky Survey Data Release 12. By leveraging the >10 year baselines for objects with repeat spectroscopy, we uncover two new changing-look quasars, and a third discovered previously. Decomposition of the multi-epoch spectra and analysis of the broad emission lines suggest that the quasar accretion disk emission dims due to rapidly decreasing accretion rates (by factors of >2.5), while disfavoring changes in intrinsic dust extinction for the two objects where these analyses are possible. Broad emission line energetics also support intrinsic dimming of quasar emission as the origin for this phenomenon rather than transient tidal disruption events or supernovae. Although our search criteria included quasars at all redshifts and transitions from either quasar-like to galaxy-like states or the reverse, all of the clear cases of changing-look quasars discovered were at relatively low-redshift (z ~ 0.2 - 0.3) and only exhibit quasar-like to galaxy-like transitions.

We report a candidate centi-parsec supermassive black hole binary (SMBHB) in the radio-quiet quasar SDSS J0159+0105 at z=0.217. With a modified lomb-scargle code GLSdeDRW and the auto-correlation analysis ACF, we detect two significant (at P>99%) periodic signals at ~741 day and ~1500 day from the 8.1-year Catalina V-band light curve of this quasar. The period ratio, which is close to 1:2, is typical of a black-hole binary system with a mass ratio of 0.05<q<0.8 according to recent numerical simulations. SDSS J0159+0105 has two SDSS spectroscopic observations separated by ~10 years. There is a significant change in the broad H-beta profile between the two epochs, which can be explained by a single broad-line region (BLR) around the binary system illuminated by the aforementioned mini-disks, or a stream of gas flowing from the circumbinary disk to one of the SMBHs. From the single BLR assumption and the orbital period t_orb ~1500 day, we estimate the total virial masses of M_SMBHB ~ 1.3x10^8 M_sun, the average distances of BLR of ~0.04pc (~50 light-day, with +/-0.3 dex uncertainty), and a SMBHB separation of d= (0.01pc)M_{8,tot}^{1/3} (T_rest/3.3yr)^{2/3} ~ 0.013 pc (15 light-day). Based on analytical work, the postulated circumbinary disk has an inner radius of 2d = 0.026 pc (30 light-day). SDSS J0159+0105 also displays unusual spectral energy distribution. The unique properties of SDSS J0159+0105 are consistent with it being a centi-parsec SMBHB. The GLSdeDRW code link: http://butler.lab.asu.edu/qso_period/

In this paper we describe and publicly release a catalogue consisting of 8445 point-like X-ray sources detected in the XMM-XXL north survey. For the 2512 AGN which have reliable spectroscopy, we present the X-ray spectral fitting which has been computed with a Bayesian approach. We have also applied an X-ray spectral stacking method to different sub-samples. We confirm the well-known Iwasawa-Taniguchi effect in our luminosity-redshift sub-samples, and argue that such an effect is due to a decrease in the covering factor of a distant obscuring torus with increasing X-ray luminosity. By comparing the distribution of the reflection fraction, the ratio of the normalization of the reflected component to the direct radiation, we find that the low-luminosity, low-redshift sub-sample had systematically higher reflection fraction values than the high-redshift, high-luminosity one. On the other hand, no significant difference is found between samples having similar luminosity but different redshift, suggesting that the structure of the torus does not evolve strongly with redshift. Contrary to previous works, we do not find evidence for an increasing photon index at high Eddington ratio. This may be an indication that the structure of the accretion disc changes as the Eddington ratio approaches unity. Comparing our X-ray spectral analysis results with the optical spectral classification, we find that $\sim20$ per cent of optical type-1 AGN show an X-ray absorbing column density higher than $10^{21.5}\,\rm{cm^{-2}}$, and about 50 per cent of type-2 AGN have an X-ray absorbing column density less than $10^{21.5}\,\rm{cm^{-2}}$. We suggest that the excess X-ray absorption shown in the high-luminosity optical type-1 AGN can be due to small-scale dust-free gas within (or close to) the broad line region, while in the low-luminosity ones it can be due to a clumpy torus with a large covering factor.

Binary active galactic nuclei (AGNs) provide clues to how gas-rich mergers trigger and fuel AGNs and how supermassive black hole (SMBH) pairs evolve in a gas-rich environment. While significant effort has been invested in their identification, the detailed properties of binary AGNs and their host galaxies are still poorly constrained. In a companion paper, we examined the nature of ionizing sources in the double nuclei of four kpc-scale binary AGNs with redshifts between 0.1~0.2. Here, we present their host galaxy morphology based on F336W (U-band) and F105W (Y-band) images taken by the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope. Our targets have double-peaked narrow emission lines and were confirmed to host binary AGNs with follow up observations. We find that kpc-scale binary AGNs occur in galaxy mergers with diverse morphological types. There are three major mergers with intermediate morphologies and a minor merger with a dominant disk component. We estimate the masses of the SMBHs from their host bulge stellar masses and obtain Eddington ratios for each AGN. Compared with a representative control sample drawn at the same redshift and stellar mass, the AGN luminosities and Eddington ratios of our binary AGNs are similar to those of single AGNs. The U-Y color maps indicate that clumpy star forming regions could significantly affect the X-ray detection of binary AGNs, e.g., the hardness ratio. Considering the weak X-ray emission in AGNs triggered in merger systems, we suggest that samples of X-ray selected AGNs may be biased against gas-rich mergers.

We present the first data release (DR1) from our UV-bright Quasar Survey (UVQS) for new $z \sim 1$ active galactic nuclei (AGN) across the sky. Using simple GALEX UV and WISE near-IR color selection criteria, we generated a list of 1450 primary candidates with $FUV < 18.5$ mag. We obtained discovery spectra, primarily on 3m-class telescopes, for 1040 of these candidates and confirmed 86% as AGN with redshifts generally at $z>0.5$. Including a small set of observed secondary candidates, we report the discovery of 217 AGN with $FUV < 18$ mag that had no previously reported spectroscopic redshift. These are excellent potential targets for UV spectroscopy before the end of the {\it Hubble Space Telescope} mission. The main data products are publicly released through the Mikulski Archive for Space Telescopes.

The Time-Domain Spectroscopic Survey (TDSS) is an SDSS-IV eBOSS subproject primarily aimed at obtaining identification spectra of ~220,000 optically-variable objects systematically selected from SDSS/Pan-STARRS1 multi-epoch imaging. We present a preview of the science enabled by TDSS, based on TDSS spectra taken over ~320 deg^2 of sky as part of the SEQUELS survey in SDSS-III, which is in part a pilot survey for eBOSS in SDSS-IV. Using the 15,746 TDSS-selected single-epoch spectra of photometrically variable objects in SEQUELS, we determine the demographics of our variability-selected sample, and investigate the unique spectral characteristics inherent in samples selected by variability. We show that variability-based selection of quasars complements color-based selection by selecting additional redder quasars, and mitigates redshift biases to produce a smooth quasar redshift distribution over a wide range of redshifts. The resulting quasar sample contains systematically higher fractions of blazars and broad absorption line quasars than from color-selected samples. Similarly, we show that M-dwarfs in the TDSS-selected stellar sample have systematically higher chromospheric active fractions than the underlying M-dwarf population, based on their H-alpha emission. TDSS also contains a large number of RR Lyrae and eclipsing binary stars with main-sequence colors, including a few composite-spectrum binaries. Finally, our visual inspection of TDSS spectra uncovers a significant number of peculiar spectra, and we highlight a few cases of these interesting objects. With a factor of ~15 more spectra, the main TDSS survey in SDSS-IV will leverage the lessons learned from these early results for a variety of time-domain science applications.

We investigate the dependence on data quality of quasar properties measured from the CIV emission line region at high redshifts. Our measurements come from 32 epochs of Sloan Digital Sky Survey (SDSS) Reverberation Mapping Project spectroscopic observations of 482 z>1.46 quasars. We compare the differences between measurements made from the single-epoch and coadded spectra, focusing on the CIV emission line because of its importance for studies of high-redshift quasar demographics and physical properties, including black hole masses. In addition to statistical errors increasing (by factors of ~2-4), we find increasing systematic offsets with decreasing S/N. The systematic difference (measurement uncertainty) in our lowest S/N (<5) subsample between the single-epoch and coadded spectrum (i) CIV equivalent width is 17A (31A), (ii) centroid wavelength is <1A (2A), and fractional velocity widths, Delta V/V, characterized by (iii) the line dispersion, sigma_l, is 0.104 (0.12), and (iv) the mean absolute deviation (MAD) is 0.072 (0.11). These remain smaller than the 1-sigma measurement uncertainties for all subsamples considered. The MAD is found to be the most robust line-width characterization. Offsets in the CIV full-width at half maximum (FWHM) velocity width and the CIV profile characterized by FWHM/sigma_l are only smaller than the statistical uncertainties when S/N>10, although offsets in lower S/N spectra exceed the statistical uncertainties by only a factor of ~1.5. Characterizing the CIV line profile by the kurtosis is the least robust property investigated, as the median systematic coadded--single-epoch measurement differences are larger than the statistical uncertainties for all S/N subsamples.

Motivated by the mathematic theory of split-complex numbers (or hyperbolic numbers, also perplex numbers) and the split-quaternion numbers (or coquaternion numbers), we define the notion of split-complex scalar field and the split-quaternion scalar field. Then we explore the cosmic evolution of these scalar fields in the background of spatially flat Friedmann-Robertson-Walker Universe. We find that both the quintessence field and the phantom field could naturally emerge in these scalar fields. Introducing the metric of field space, these theories fall into a subclass of the multi-field theories which have been extensively studied in inflationary cosmology.

We present a reverberation mapping (RM) experiment that combines broad- and intermediate-band photometry; it is the first such attempt targeting a sample of 13 quasars at $0.2<z<0.9$. The quasars were selected to have strong H$\alpha$ or H$\beta$ emission lines that are located in one of three intermediate bands (with FWHM around 200 \AA) centered at 8045, 8505, and 9171 \AA. The imaging observations were carried out in the intermediate bands and the broad $i$ and $z$ bands using the prime-focus imager 90Prime on the 2.3m Bok telescope. Because of the large ($\sim$1 deg$^2$) field-of-view (FoV) of 90Prime, we were able to include the 13 quasars within only five telescope pointings or fields. The five fields were repeatedly observed over 20--30 epochs that were unevenly distributed over a duration of 5--6 months. The combination of the broad- and intermediate-band photometry allows us to derive accurate light curves for both optical continuum (from the accretion disk) and line (from the broad-line region, or BLR) emission. We detect H$\alpha$ time lags between the continuum and line emission in 6 quasars. These quasars are at a relatively low redshift range $0.2<z<0.4$. The measured lags are consistent with the current BLR size-luminosity relation for H$\beta$ at $z<0.3$. While this experiment appears successful in detecting lags of the bright H$\alpha$ line, further investigation is required to see if it can also be applied to the fainter H$\beta$ line for quasars at higher redshifts. Finally we demonstrate that by using a small telescope with a large FoV, intermediate-band photometric RM can be efficiently executed for a large sample of quasars at $z>0.2$.

Reverberation mapping (RM) measurements of broad-line region (BLR) lags in z>0.3 quasars are important for directly measuring black hole masses in these distant objects, but so far there have been limited attempts and success given the practical difficulties of RM in this regime. Here we report preliminary results of 15 BLR lag measurements from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, a dedicated RM program with multi-object spectroscopy designed for RM over a wide redshift range. The lags are based on the 2014 spectroscopic light curves alone (32 epochs over 6 months) and focus on the Hbeta and MgII broad lines in the 100 lowest-redshift (z<0.8) quasars included in SDSS-RM; they represent a small subset of the lags that SDSS-RM (including 849 quasars to z~4.5) is expected to deliver. The reported preliminary lag measurements are for intermediate-luminosity quasars at 0.3<~z<0.8, including 9 Hbeta lags and 6 MgII lags, for the first time extending RM results to this redshift-luminosity regime and providing direct quasar black hole mass estimates over ~ half of cosmic time. The MgII lags also increase the number of known MgII lags by several-fold, and start to explore the utility of MgII for RM at high redshift. The location of these new lags at higher redshifts on the observed BLR size-luminosity relationship is statistically consistent with previous Hbeta results at z<0.3. However, an independent constraint on the relationship slope at z>0.3 is not yet possible due to the limitations in our current sample. Our results demonstrate the general feasibility and potential of multi-object RM for z>0.3 quasars.

We study the rest-frame optical properties of 74 luminous (L_bol=10^46.2-48.2 erg/s), 1.5<z<3.5 broad-line quasars with near-IR (JHK) slit spectroscopy. Systemic redshifts based on the peak of the [OIII]5007 line reveal that redshift estimates from the rest-frame UV broad emission lines (mostly MgII) are intrinsically uncertain by ~ 200 km/s (measurement errors accounted for). The overall full-width-at-half-maximum of the narrow [OIII] line is ~ 1000 km/s on average. A significant fraction of the total [OIII] flux (~ 40%) is in a blueshifted wing component with a median velocity offset of ~ 700 km/s, indicative of ionized outflows within a few kpc from the nucleus; we do not find evidence of significant [OIII] flux beyond ~ 10 kpc in our slit spectroscopy. The [OIII] line is noticeably more asymmetric and weaker than that in typical less luminous low-z quasars. However, when matched in quasar continuum luminosity, low-z quasars have similar [OIII] profiles and strengths as these high-z systems. Therefore the exceptionally large width and blueshifted wing, and the relatively weak strength of [OIII] in high-z luminous quasars are mostly a luminosity effect rather than redshift evolution. The Hbeta-[OIII] region of these high-z quasars displays a similar spectral diversity and Eigenvector 1 correlations with anti-correlated [OIII] and optical FeII strengths, as seen in low-z quasars; but the average broad Hbeta width is larger by 25% than typical low-z quasars, indicating more massive black holes in these high-z systems. These results highlight the importance of understanding [OIII] emission in the general context of quasar parameter space in order to understand quasar feedback in the form of [OIII] outflows. The calibrated one-dimensional near-IR spectra are made publicly available, along with a composite spectrum.

We present high-resolution observations of a quiescent solar prominence which was consisted of a vertical and a horizontal foot encircled by an overlying spine, and counter-streaming mass flows were ubiquitous in the prominence. While the horizontal foot and the spine were connecting to the solar surface, the vertical foot was suspended above the solar surface and supported by a semicircular bubble structure. The bubble first collapsed and then reformed at a similar height, finally, it started to oscillate for a long time. We find that the collapsing and oscillation of the bubble boundary were tightly associated with a flare-like feature located at the bottom of the bubble. Based on the observational results, we propose that the prominence should be composed of an overlying horizontal spine encircling a low-lying horizontal and a vertical foot, in which the horizontal foot was consisted of shorter field lines running partially along the spine and with the both ends connecting to the solar surface, while the vertical foot was consisted of piling-up dips due to the sagging of the spine fields and supported by a bipolar magnetic system formed by parasitic polarities (i.e., the bubble). The upflows in the vertical foot were possibly caused by the magnetic reconnection at the separator between the bubble and the overlying dips, which intruded into the persistent downflow field and formed the picture of counter-streaming mass flows. In addition, the counter-streaming flows in the horizontal foot were possibly caused by the imbalanced pressure at the both ends.

We report the discovery of a new "changing-look" quasar, SDSS J101152.98+544206.4, through repeat spectroscopy from the Time Domain Spectroscopic Survey. This is an addition to a small but growing set of quasars whose blue continua and broad optical emission lines have been observed to decline by a large factor on a time scale of approximately a decade. The 5100 Angstrom monochromatic continuum luminosity of this quasar drops by a factor of > 9.8 in a rest-frame time interval of < 9.7 years, while the broad H-alpha luminosity drops by a factor of 55 in the same amount of time. The width of the broad H-alpha line increases in the dim state such that the black hole mass derived from the appropriate single-epoch scaling relation agrees between the two epochs within a factor of 3. The fluxes of the narrow emission lines do not appear to change between epochs. The light curve obtained by the Catalina Sky Survey suggests that the transition occurs within a rest-frame time interval of approximately 500 days. We examine three possible mechanisms for this transition suggested in the recent literature. An abrupt change in the reddening towards the central engine is disfavored by the substantial difference between the timescale to obscure the central engine and the observed timescale of the transition. A decaying tidal disruption flare is consistent with the decay rate of the light curve but not with the prolonged bright state preceding the decay, nor can this scenario provide the power required by the luminosities of the emission lines. An abrupt drop in the accretion rate onto the supermassive black hole appears to be the most plausible explanation for the rapid dimming.

Quasar host galaxies are key for understanding the relation between galaxies and the supermassive black holes (SMBHs) at their centers. We present a study of 191 broad-line quasars and their host galaxies at z < 1, using high signal-to-noise ratio (SNR) spectra produced by the Sloan Digital Sky Survey Reverberation Mapping project. Clear detection of stellar absorption lines allows a reliable decomposition of the observed spectra into nuclear and host components, using spectral models of quasar and stellar radiations as well as emission lines from the interstellar medium. We estimate age, mass (M*), and velocity dispersion (sigma*) of the host stars, the star formation rate (SFR), quasar luminosity, and SMBH mass (Mbh), for each object. The quasars are preferentially hosted by massive galaxies with M* ~ 10^{11} Msun characterized by stellar ages around a billion years, which coincides with the transition phase of normal galaxies from the blue cloud to the red sequence. The host galaxies have relatively low SFRs and fall below the main sequence of star-forming galaxies at similar redshifts. These facts suggest that the hosts have experienced an episode of major star formation sometime in the past billion years, which was subsequently quenched or suppressed. The derived Mbh - sigma* and Mbh - M* relations agree with our past measurements and are consistent with no evolution from the local Universe. The present analysis demonstrates that reliable measurements of stellar properties of quasar host galaxies are possible with high-SNR fiber spectra, which will be acquired in large numbers with future powerful instruments such as the Subaru Prime Focus Spectrograph.

We explore the variability of quasars in the MgII and Hbeta broad emission lines and UV/optical continuum emission using the Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM). This is the largest spectroscopic study of quasar variability to date: our study includes 29 spectroscopic epochs from SDSS-RM over $6$ months, containing 357 quasars with MgII and 41 quasars with Hbeta . On longer timescales, the study is also supplemented with two-epoch data from SDSS-I/II. The SDSS-I/II data include an additional $2854$ quasars with MgII and 572 quasars with Hbeta. The MgII emission line is significantly variable ($\Delta f/f$ 10% on 100-day timescales), a necessary prerequisite for its use for reverberation mapping studies. The data also confirm that continuum variability increases with timescale and decreases with luminosity, and the continuum light curves are consistent with a damped random-walk model on rest-frame timescales of $\gtrsim 5$ days. We compare the emission-line and continuum variability to investigate the structure of the broad-line region. Broad-line variability shows a shallower increase with timescale compared to the continuum emission, demonstrating that the broad-line transfer function is not a $\delta$-function. Hbeta is more variable than MgII (roughly by a factor of $1.5$), suggesting different excitation mechanisms, optical depths and/or geometrical configuration for each emission line. The ensemble spectroscopic variability measurements enabled by the SDSS-RM project have important consequences for future studies of reverberation mapping and black hole mass estimation of $1<z<2$ quasars.

We measure the two-point clustering of spectroscopically confirmed quasars from the final sample of the Baryon Oscillation Spectroscopic Survey (BOSS) on comoving scales of 4 < s < 22 Mpc/h. The sample covers 6950 deg^2 (~ 19 (Gpc/h)^3) and, over the redshift range 2.2 < z < 2.8, contains 55,826 homogeneously selected quasars, which is twice as many as in any similar work. We deduce b_Q = 3.54 +/- 0.10 ; the most precise measurement of quasar bias to date at these redshifts. This corresponds to a host halo mass of ~ 2 x 10^12 ~ M_sun/h with an implied quasar duty cycle of ~1 percent. The real-space projected correlation function is well-fit by a power law of index -2 and correlation length r0 = (8.12 +/- 0.22), Mpc/h over scales of 4 < rp < 25 ~ Mpc/h. To better study the evolution of quasar clustering at moderate redshift, we extend the redshift range of our study to z ~ 3.4 and measure the bias and correlation length of three subsamples over 2.2 < z < 3.4. We find no significant evolution of r0 or bias over this range, implying that the host halo mass of quasars decreases somewhat with increasing redshift. We find quasar clustering remains similar over a decade in luminosity, contradicting a scenario in which quasar luminosity is monotonically related to halo mass at z ~ 2.5. Our results are broadly consistent with previous BOSS measurements, but they yield more precise constraints based upon a larger and more uniform data set.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

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 report the discovery of rapid variations of a high-velocity CIV broad absorption line trough in the quasar SDSS J141007.74+541203.3. This object was intensively observed in 2014 as a part of the Sloan Digital Sky Survey Reverberation Mapping Project, during which 32 epochs of spectroscopy were obtained with the Baryon Oscillation Spectroscopic Survey spectrograph. We observe significant (>4sigma) variability in the equivalent width of the broad (~4000 km/s wide) CIV trough on rest-frame timescales as short as 1.20 days (~29 hours), the shortest broad absorption line variability timescale yet reported. The equivalent width varied by ~10% on these short timescales, and by about a factor of two over the duration of the campaign. We evaluate several potential causes of the variability, concluding that the most likely cause is a rapid response to changes in the incident ionizing continuum. If the outflow is at a radius where the recombination rate is higher than the ionization rate, the timescale of variability places a lower limit on the density of the absorbing gas of n_e > 3.9 x 10^5 cm^-3. The broad absorption line variability characteristics of this quasar are consistent with those observed in previous studies of quasars, indicating that such short-term variability may in fact be common and thus can be used to learn about outflow characteristics and contributions to quasar/host-galaxy feedback scenarios.

We present host stellar velocity dispersion measurements for a sample of 88 broad-line quasars at 0.1<z<1 (46 at z>0.6) from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. High signal-to-noise ratio coadded spectra (average S/N~30 per 69 km/s pixel) from SDSS-RM allowed decomposition of the host and quasar spectra, and measurement of the host stellar velocity dispersions and black hole (BH) masses using the single-epoch (SE) virial method. The large sample size and dynamic range in luminosity (L5100=10^(43.2-44.7) erg/s) lead to the first clear detection of a correlation between SE virial BH mass and host stellar velocity dispersion far beyond the local universe. However, the observed correlation is significantly flatter than the local relation, suggesting that there are selection biases in high-z luminosity-threshold quasar samples for such studies. Our uniform sample and analysis enable an investigation of the redshift evolution of the M-sigma relation free of caveats by comparing different samples/analyses at disjoint redshifts. We do not observe evolution of the M-sigma relation in our sample, up to z~1, but there is an indication that the relation flattens towards higher redshifts. Coupled with the increasing threshold luminosity with redshift in our sample, this again suggests certain selection biases are at work, and simple simulations demonstrate that a constant M-sigma relation is favored to z~1. Our results highlight the scientific potential of deep coadded spectroscopy from quasar monitoring programs, and offer a new path to probe the co-evolution of BHs and galaxies at earlier times.

The diverse properties of broad-line quasars appear to follow a well-defined main sequence along which the optical FeII strength increases. It has been suggested that this sequence is mainly driven by the Eddington ratio (L/L_Edd) of the black hole (BH) accretion. Shen & Ho demonstrated with quasar clustering analysis that the average BH mass decreases with increasing FeII strength when quasar luminosity is fixed, consistent with this suggestion. Here we perform an independent test by measuring the stellar velocity dispersion sigma* (hence the BH mass via the M-sigma* relation) from decomposed host spectra in low-redshift Sloan Digital Sky Survey quasars. We found that at fixed quasar luminosity, sigma* systematically decreases with increasing FeII strength, confirming that Eddington ratio increases with FeII strength. We also found that at fixed luminosity and FeII strength, there is little dependence of sigma* on the broad Hbeta FWHM. These new results reinforce the framework put forward by Shen & Ho that Eddington ratio and orientation govern most of the diversity seen in broad-line quasar properties.

We present an X-ray and multiwavelength study of 33 weak emission-line quasars (WLQs) and 18 quasars that are analogs of the extreme WLQ, PHL 1811, at z ~ 0.5-2.9. New Chandra 1.5-9.5 ks exploratory observations were obtained for 32 objects while the others have archival X-ray observations. Significant fractions of these luminous type 1 quasars are distinctly X-ray weak compared to typical quasars, including 16 (48%) of the WLQs and 17 (94%) of the PHL 1811 analogs with average X-ray weakness factors of 17 and 39, respectively. We measure a relatively hard ($\Gamma=1.16_{-0.32}^{+0.37}$) effective power-law photon index for a stack of the X-ray weak subsample, suggesting X-ray absorption, and spectral analysis of one PHL 1811 analog, J1521+5202, also indicates significant intrinsic X-ray absorption. We compare composite SDSS spectra for the X-ray weak and X-ray normal populations and find several optical-UV tracers of X-ray weakness; e.g., Fe II rest-frame equivalent width and relative color. We describe how orientation effects under our previously proposed "shielding-gas" scenario can likely unify the X-ray weak and X-ray normal populations. We suggest that the shielding gas may naturally be understood as a geometrically thick inner accretion disk that shields the broad line region from the ionizing continuum. If WLQs and PHL 1811 analogs have very high Eddington ratios, the inner disk could be significantly puffed up (e.g., a slim disk). Shielding of the broad emission-line region by a geometrically thick disk may have a significant role in setting the broad distributions of C IV rest-frame equivalent width and blueshift for quasars more generally.

Fast radio bursts (FRBs) are one of the most tantalizing mysteries of the radio sky; their progenitors and origins remain unknown and until now no rapid multiwavelength follow-up of an FRB has been possible. New instrumentation has decreased the time between observation and discovery from years to seconds, and enables polarimetry to be performed on FRBs for the first time. We have discovered an FRB (FRB 140514) in real-time on 14 May, 2014 at 17:14:11.06 UTC at the Parkes radio telescope and triggered follow-up at other wavelengths within hours of the event. FRB 140514 was found with a dispersion measure (DM) of 562.7(6) cm$^{-3}$ pc, giving an upper limit on source redshift of $z \lesssim 0.5$. FRB 140514 was found to be 21$\pm$7% (3-$\sigma$) circularly polarized on the leading edge with a 1-$\sigma$ upper limit on linear polarization $<10%$. We conclude that this polarization is intrinsic to the FRB. If there was any intrinsic linear polarization, as might be expected from coherent emission, then it may have been depolarized by Faraday rotation caused by passing through strong magnetic fields and/or high density environments. FRB 140514 was discovered during a campaign to re-observe known FRB fields, and lies close to a previous discovery, FRB 110220; based on the difference in DMs of these bursts and time-on-sky arguments, we attribute the proximity to sampling bias and conclude that they are distinct objects. Follow-up conducted by 12 telescopes observing from X-ray to radio wavelengths was unable to identify a variable multiwavelength counterpart, allowing us to rule out models in which FRBs originate from nearby ($z < 0.3$) supernovae and long duration gamma-ray bursts.

Quasars are rapidly accreting supermassive black holes at the center of massive galaxies. They display a broad range of properties across all wavelengths, reflecting the diversity in the physical conditions of the regions close to the central engine. These properties, however, are not random, but form well-defined trends. The dominant trend is known as Eigenvector 1, where many properties correlate with the strength of optical iron and [OIII] emission. The main physical driver of Eigenvector 1 has long been suspected to be the quasar luminosity normalized by the mass of the hole (the Eddington ratio), an important quantity of the black hole accretion process. But a definitive proof has been missing. Here we report an analysis of archival data that reveals that Eddington ratio indeed drives Eigenvector 1. We also find that orientation plays a significant role in determining the observed kinematics of the gas, implying a flattened, disklike geometry for the fast-moving clouds close to the hole. Our results show that most of the diversity of quasar phenomenology can be unified with two simple quantities, Eddington ratio and orientation.

We present the first stereoscopic and Doppler observations of simultaneous transverse oscillations of a prominence and a filament and longitudinal oscillation of another filament launched by a single shock wave. Using H-alpha Doppler observations, we derive the three-dimensional oscillation velocities at different heights along the prominence axis. The results indicate that the prominence has a larger oscillation amplitude and damping time at higher altitude, but the periods at different heights are the same (i.e., 13.5 minutes). This suggests that the prominence oscillates like a linear vertical rigid body with one end anchored on the Sun. One of the filaments shows weak transverse oscillation after the passing of the shock, which is possibly due to the low altitude of the filament and the weakening (due to reflection) of the shock wave before the interaction. Large amplitude longitudinal oscillation is observed in the other filament after the passing of the shock wave. The velocity amplitude and period are about 26.8 km/s and 80.3 minutes, respectively. We propose that the orientation of a filament or prominence relative to the normal vector of the incoming shock should be an important factor for launching transverse or longitudinal filament oscillations. In addition, the restoring forces of the transverse prominence are most likely due to the coupling of gravity and magnetic tension of the supporting magnetic field, while that for the longitudinal filament oscillation is probably the resultant force of gravity and magnetic pressure.

Based on an updated Hbeta reverberation mapping (RM) sample of 44 nearby active galactic nuclei (AGN), we propose a novel approach for black hole (BH) mass estimation using two filtered luminosities computed from single-epoch (SE) AGN spectra around the Hbeta region. We found that the two optimal-filter luminosities extract virial information (size and virial velocity of the broad line region, BLR) from the spectra, justifying their usage in this empirical BH mass estimator. The major advantages of this new recipe over traditional SE BH mass estimators utilizing continuum luminosity and broad line width are: 1) it has a smaller intrinsic scatter of 0.28 dex calibrated against RM masses; 2) it is extremely simple to use in practice, without any need to decompose the spectrum; 3) it produces unambiguous and highly repeatable results even with low signal-to-noise spectra. The combination of the two luminosities can also cancel out, to some extent, systematic luminosity errors potentially introduced by uncertainties in distance or flux calibration. In addition, we recalibrated the traditional SE mass estimators using broad Hbeta full width at half maximum (FWHM) and monochromatic continuum luminosity at 5100A (L_5100). We found that using the best-fit slopes on FWHM and L_5100 (derived from fitting the BLR radius-luminosity relation and the correlation between rms line dispersion and SE FWHM, respectively) rather than simple assumptions (e.g., 0.5 for L_5100 and 2 for FWHM) leads to more precise SE mass estimates, improving the intrinsic scatter from 0.41 dex to 0.36 dex with respect to the RM masses. We compared different estimators and discussed their applications to the Sloan Digital Sky Survey (SDSS) quasar sample. Due to the limitations of the current RM sample, application of any SE recipe calibrated against RM masses to distant quasars should be treated with caution.

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 explore a cyclic universe due to phantom and quintessence fields. We find that, in every cycle of the evolution of the universe, the phantom dominates the cosmic early history and quintessence dominates the cosmic far future. In this model of universe, there are infinite cycles of expansion and contraction. Different from the inflationary universe, the corresponding cosmic space-time is geodesically complete and quantum stable. But similar to the Cyclic Model, the flatness problem, the horizon problem and the large scale structure of the universe can be explained in this cyclic universe.

We consider how the profile and multi-year variability properties of a large sample of C IV Broad Absorption Line (BAL) troughs change when BALs from Si IV and/or Al III are present at corresponding velocities, indicating that the line-of-sight intercepts at least some lower ionization gas. We derive a number of observational results for C IV BALs separated according to the presence or absence of accompanying lower ionization transitions, including measurements of composite profile shapes, equivalent width (EW), characteristic velocities, composite variation profiles, and EW variability. We also measure the correlations between EW and fractional-EW variability for C IV, Si IV, and Al III. Our measurements reveal the basic correlated changes between ionization level, kinematics, and column density expected in accretion-disk wind models; e.g., lines-of-sight including lower ionization material generally show deeper and broader C IV troughs that have smaller minimum velocities and that are less variable. Many C IV BALs with no accompanying Si IV or Al III BALs may have only mild or no saturation.

The Sky Brightness Monitor (SBM) is an important instrument to measure the brightness level for the sky condition, which is a critical parameter for judging a site for solar coronal observations. In this paper we present an automatic method for the processing of SBM data in large quantity, which can separate the regions of the Sun and the nearby sky as well as recognize the regions of the supporting arms in the field of view. These processes are implemented on the data acquired by more than one SBM instruments during our site survey project in western China. An analysis applying the result from our processes has been done for the assessment of the scattered-light levels by the instrument. Those results are considerably significant for further investigations and studies, notably to derive a series of the other important atmospheric parameters such as extinctions, aerosol content and precipitable water vapor content for candidate sites. Our processes also provide a possible way for full-disk solar telescopes to track the Sun without an extra guiding system.

A small fraction of quasars have long been known to show bulk velocity offsets in the broad Balmer lines with respect to the systemic redshift of the host galaxy. Models to explain these offsets usually invoke broad-line region gas kinematics/asymmetry around single black holes (BHs), orbital motion of massive (~sub-pc) binary black holes (BBHs), or recoil BHs, but single-epoch spectra are unable to distinguish between these scenarios. The line-of-sight (LOS) radial velocity (RV) shifts from long-term spectroscopic monitoring can be used to test the BBH hypothesis. We have selected a sample of 399 quasars with offset broad H-beta lines from the SDSS DR7 quasar catalog, and have conducted second-epoch optical spectroscopy for 50 of them. Combined with the existing SDSS spectra, the new observations enable us to constrain the LOS RV shifts of broad H-beta lines with a rest-frame baseline of a few years to nearly a decade. Using cross-correlation analysis, we detect significant (99% confidence) radial accelerations in the broad H-beta lines in 24 of the 50 objects. We suggest that 9 of the 24 detections are sub-pc BBH candidates, which show consistent velocity shifts independently measured from a second broad line (either H-alpha or Mg II) without significant changes in the broad-line profiles. Combining the results on the general quasar population studied in Paper I, we find a tentative anti-correlation between the velocity offset in the first-epoch spectrum and the average acceleration between two epochs, which could be explained by orbital phase modulation when the time separation between two epochs is a non-negligible fraction of the orbital period of the motion causing the line displacement. We discuss the implications of our results for the identification of sub-pc BBH candidates in offset-line quasars and for the constraints on their frequency and orbital parameters. [Abridged]

The broad emission lines commonly seen in quasar spectra have velocity widths of a few per cent of the speed of light, so special- and general-relativistic effects have a significant influence on the line profile. We have determined the redshift of the broad H-beta line in the quasar rest frame (determined from the core component of the [OIII] line) for over 20,000 quasars from the Sloan Digital Sky Survey DR7 quasar catalog. The mean redshift as a function of line width is approximately consistent with the relativistic redshift that is expected if the line originates in a randomly oriented Keplerian disk that is obscured when the inclination of the disk to the line of sight exceeds ~30-45 degrees, consistent with simple AGN unification schemes. This result also implies that the net line-of-sight inflow/outflow velocities in the broad-line region are much less than the Keplerian velocity when averaged over a large sample of quasars with a given line width.

Winking (oscillating) filaments have been observed for many years. However, observations of successive winking filaments in one event have not been reported yet. In this paper, we present the observations of a chain of winking filaments and a subsequent jet that are observed right after the X2.1 are in AR11283. The event also produced an Extreme-ultraviolet (EUV) wave that has two components: upward dome-like wave (850 km/s) and lateral surface wave (554 km/s) which was very weak (or invisible) in imaging observations. By analyzing the temporal and spatial relationships between the oscillating filaments and the EUV waves, we propose that all the winking filaments and the jet were triggered by the weak (or invisible) lateral surface EUV wave. The oscillation of the filaments last for two or three cycles, and their periods, Doppler velocity amplitudes, and damping times are 11 - 22 minutes, 6 - 14 km/s, and 25 - 60 minutes, respectively. We further estimate the radial component magnetic field and the maximum kinetic energy of the filaments, and they are 5 - 10 Gauss and ~ 10^19 J, respectively. The estimated maximum kinetic energy is comparable to the minimum energy of ordinary EUV waves, suggesting that EUV waves can efficiently launch filament oscillations on their path. Based on our analysis results, we conclude that the EUV wave is a good agent for triggering and connecting successive but separated solar activities in the solar atmosphere, and it is also important for producing solar sympathetic eruptions.

This whitepaper describes how the VLASS could be designed in a manner to allow the identification of candidate dual active galactic nuclei (AGN) at separations <7 kpc. Dual AGN represent a clear marker of two supermassive black holes within an ongoing merger. A dual AGN survey will provide a wealth of studies in structure growth and gravitational-wave science. Radio wavelengths are ideal for identifying close pairs, as disturbed stellar and gaseous material can obscure their presence in optical and shorter wavelengths. With sufficiently high resolution and sensitivity, a large-scale radio imaging survey like the VLASS will uncover many of these systems and provide the means to broadly study the radio properties of candidate dual systems revealed at other wavelengths. We determine that the ideal survey for our purposes will be at as high a resolution as possible, with significantly more science return in A array at L-band or higher, or B array at C-band or higher. We describe a range of potential survey parameters within this document. Based on the analysis outlined in this whitepaper, our ideal survey would create a catalogue of $\gtrsim$100 dual AGN in either: 1) a medium-sensitivity (~1 mJy detection threshold), wide-field (few thousand square degree) survey, or 2) a high-sensitivity (~10 uJy threshold) survey of several hundred square degrees.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

The magnetic breakout model has been widely used to explain solar eruptive activities. Here, we apply it to explain successive filament eruptions occurred in a quadrupolar magnetic source region. Based on the high temporal and spatial resolution, multi-wavelengths observations taken by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO), we find some signatures that support the occurrence of breakout-like external reconnection just before the start of the successive filament eruptions. Furthermore, the extrapolated three-dimensional coronal field also reveals that the magnetic topology above the quadrupolar source region resembles that of the breakout model. We propose a possible mechanism within the framework of the breakout model to interpret the successive filament eruptions, in which the so-called magnetic implosion mechanism is firstly introduced to be the physical linkage of successive filament eruptions. We conclude that the structural properties of coronal fields are important for producing successive filament eruptions.

We present a detailed investigation of the variability of 428 C IV and 235 Si IV Broad Absorption Line (BAL) troughs identified in multi-epoch observations of 291 quasars by the Sloan Digital Sky Survey-I/II/III. These observations primarily sample rest-frame timescales of 1-3 yr over which significant rearrangement of the BAL wind is expected. We derive a number of observational results on, e.g., the frequency of BAL variability, the velocity range over which BAL variability occurs, the primary observed form of BAL-trough variability, the dependence of BAL variability upon timescale, the frequency of BAL strengthening vs. weakening, correlations between BAL variability and BAL-trough profiles, relations between C IV and Si IV BAL variability, coordinated multi-trough variability, and BAL variations as a function of quasar properties. We assess implications of these observational results for quasar winds. Our results support models where most BAL absorption is formed within an order-of-magnitude of the wind-launching radius, although a significant minority of BAL troughs may arise on larger scales. We estimate an average lifetime for a BAL trough along our line-of-sight of a few thousand years. BAL disappearance and emergence events appear to be extremes of general BAL variability, rather than being qualitatively distinct phenomena. We derive the parameters of a random-walk model for BAL EW variability, finding that this model can acceptably describe some key aspects of EW variability. The coordinated trough variability of BAL quasars with multiple troughs suggests that changes in "shielding gas" may play a significant role in driving general BAL variability.

We perform a systematic search for sub-parsec binary supermassive black holes (BHs) in normal broad line quasars at z<0.8, using multi-epoch SDSS spectroscopy of the broad Hbeta line. Our working model is that: only one of the two BHs in the binary is active, and dynamically dominates its own broad line region (BLR); the inactive companion BH is orbiting at a distance of a few R_BLR, where R_BLR~0.01-0.1 pc is the BLR size. We search for the expected line-of-sight acceleration of the broad line velocity from binary orbital motion by cross-correlating SDSS spectra from two epochs separated by up to several years in the quasar restframe. Out of ~700 pairs of spectra for which we have good measurements of the velocity shift between two epochs (1-sigma error~40 km/s), we detect 28 systems with significant velocity shifts in broad Hbeta, among which seven are the best candidates for the hypothesized binaries. We use the distribution of the observed accelerations (mostly non-detections) to place constraints on the abundance of such binary systems among the general quasar population. Excess variance in the velocity shift is inferred for observations separated by longer than 0.4 yr (quasar restframe). Attributing all the excess to binary motion would imply that most of the quasars in this sample must be in binaries; that the inactive BH must be on average more massive than the active one; and that the binary separation is at most a few times the size of the BLR. However, if this excess variance is partly or largely due to long-term broad line variability, the requirement of a large population of close binaries is much weakened or even disfavored for massive companions. [Abridged]

We present an observational study of a quasi-periodic fast propagating (QFP) magnetosonic wave on 2012, April 23. The multiple wave trains were observed along an active region open loop system which has a divergence geometry. The wave trains were first observed in 171 A observations at a distance of 150 Mm from the footpoint of the guiding loop system and with a speed of 689 km/s, then they appeared in 193 A observations after their interaction with a perpendicular underlaying loop system on the path, in the meantime, the wave speed decelerated to 343 km/s quickly within a short timescale. The sudden deceleration of the wave trains and their appearance in 193 A observations caused by the interaction are interpreted through geometric effect and the density increase of the guiding loop system, respectively. On the other hand, with Wavelet and Fourier analysis methods we find that the wave trains has a common period of 80 s with the associated flare. In addition, a few low frequencies are also identified in the QFP wave. We propose that the generation of the period of 80 s was caused by the periodic releasing of energy busts through some nonlinear processes in magnetic reconnection or the so-called oscillatory reconnection mechanism, while the low frequencies detected in the QFP wave were possibly the manifestations of the leakage of pressure-driven oscillations from the photosphere or chromosphere, which could be an important source for driving QFP waves in the low corona. Our observational results also indicate that the properties of the guiding magnetic structure such as the distributions of magnetic field and density as well as geometry are crucial for modulating the propagation behaviors of QFP waves. Therefore, QFP waves could be used for remote diagnostics of the local physical properties of the solar corona.

We present observations of the diffraction, refraction, and reflection of a global extreme-ultraviolet (EUV) wave propagating in the solar corona. These intriguing phenomena are observed when the wave interacts with two remote active regions, and they together exhibit the wave property of this EUV wave. When the wave approached AR11465, it became weaker and finally disappeared in the active region, but a few minutes latter a new wavefront appeared behind the active region, and it was not concentric with the incoming wave. In addition, a reflected wave was also observed simultaneously on the wave incoming side. When the wave approached AR11459, it transmitted through the active region directly and without reflection. The formation of the new wavefront and the transmission could be explained with diffraction and refraction effects, respectively. We propose that the different behaviors observed during the interactions may caused by different speed gradients at the boundaries of the two active regions. For the origin of the EUV wave, we find that it formed ahead of a group of expanding loops a few minutes after the start of the loops' expansion, which represents the initiation of the associated coronal mass ejection (CME). Based on these results, we conclude that the EUV wave should be a nonlinear magnetosonic wave or shock driven by the associated CME, which propagated faster than the ambient fast-mode speed and gradually slowed down to an ordinary linear wave. Our observations support the hybrid model that includes both fast wave and slow non-wave components.

We present detailed analysis of an extreme ultraviolet (EUV) wave and its interaction with active region (AR) loops observed by the Solar Dynamics Observatory/Atmospheric Imaging Assembly and the Hinode EUV Imaging Spectrometer (EIS). This wave was initiated from AR 11261 on 2011 August 4 and propagated at velocities of 430--910 km\s. It was observed to traverse another AR and cross over a filament channel on its path. The EUV wave perturbed neighboring AR loops and excited a disturbance that propagated toward the footpoints of these loops. EIS observations of AR loops revealed that at the time of the wave transit, the original redshift increased by about 3 km\s, while the original blueshift decreased slightly. After the wave transit, these changes were reversed. When the EUV wave arrived at the boundary of a polar coronal hole, two reflected waves were successively produced and part of them propagated above the solar limb. The first reflected wave above the solar limb encountered a large-scale loop system on its path, and a secondary wave rapidly emerged 144 Mm ahead of it at a higher speed. These findings can be explained in the framework of a fast-mode magnetosonic wave interpretation for EUV waves, in which observed EUV waves are generated by expanding coronal mass ejections.

We report the discovery in the Sloan Digital Sky Survey and the SDSS-III Baryon Oscillation Spectroscopic Survey of seventeen broad absorption line (BAL) quasars with high-ionization troughs that include absorption redshifted relative to the quasar rest frame. The redshifted troughs extend to velocities up to v=12,000 km/s and the trough widths exceed 3000 km/s in all but one case. Approximately 1 in 1000 BAL quasars with blueshifted C IV absorption also has redshifted C IV absorption; objects with C IV absorption present only at redshifted velocities are roughly four times rarer. In more than half of our objects, redshifted absorption is seen in C II or Al III as well as C IV, making low-ionization absorption at least ten times more common among BAL quasars with redshifted troughs than among standard BAL quasars. However, the C IV absorption equivalent widths in our objects are on average smaller than those of standard BAL quasars with low-ionization absorption. We consider several possible ways of generating redshifted absorption. The two most likely possibilities may be at work simultaneously, in the same objects or in different ones. Rotationally dominated outflows seen against a quasar's extended continuum source can produce redshifted and blueshifted absorption, but variability consistent with this scenario is seen in only one of the four objects with multiple spectra. The infall of relatively dense and low-ionization gas to radii as small as 400 Schwarzschild radii can in principle explain the observed range of trough profiles, but current models do not easily explain the origin and survival of such gas. Whatever the origin(s) of the absorbing gas in these objects, it must be located at small radii to explain its large redshifted velocities, and thus offers a novel probe of the inner regions of quasars.

We present a measurement of the Type I quasar luminosity function at z=5 using a large sample of spectroscopically confirmed quasars selected from optical imaging data. We measure the bright end (M_1450<-26) with Sloan Digital Sky Survey (SDSS) data covering ~6000 deg^2, then extend to lower luminosities (M_1450<-24) with newly discovered, faint z~5 quasars selected from 235 deg^2 of deep, coadded imaging in the SDSS Stripe 82 region (the celestial equator in the Southern Galactic Cap). The faint sample includes 14 quasars with spectra obtained as ancillary science targets in the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), and 59 quasars observed at the MMT and Magellan telescopes. We construct a well-defined sample of 4.7<z<5.1 quasars that is highly complete, with 73 spectroscopic identifications out of 92 candidates. Our color selection method is also highly efficient: of the 73 spectra obtained, 71 are high redshift quasars. These observations reach below the break in the luminosity function (M_1450* ~ -27). The bright end slope is steep (beta <~ -4), with a constraint of beta < -3.1 at 95% confidence. The break luminosity appears to evolve strongly at high redshift, providing an explanation for the flattening of the bright end slope reported previously. We find a factor of ~2 greater decrease in the number density of luminous quasars (M_1450<-26) from z=5 to z=6 than from z=4 to z=5, suggesting a more rapid decline in quasar activity at high redshift than found in previous surveys. Our model for the quasar luminosity function predicts that quasars generate ~30% of the ionizing photons required to keep the universe ionized at z=5.

We review our knowledge of the most basic properties of the AGN obscuring region - its location, scale, symmetry, and mean covering factor - and discuss new evidence on the distribution of covering factors in a sample of ~9000 quasars with WISE, UKIDSS, and SDSS photometry. The obscuring regions of AGN may be in some ways more complex than we thought - multi-scale, not symmetric, chaotic - and in some ways simpler - with no dependence on luminosity, and a covering factor distribution that may be determined by the simplest of considerations - e.g. random misalignments.

I review the current status of quasar black hole (BH) mass estimations. Spectroscopic methods have been developed to estimate BH mass in broad line quasars to an accuracy of ~0.5 dex. Despite their popularity, significant issues and confusion remain regarding these mass estimators. I provide an in-depth discussion on the merits and caveats of the single-epoch (SE) virial BH mass estimators, and a detailed derivation of the statistical biases of these SE mass estimates resulting from their errors. I show that error-induced sample biases on the order of a factor of several are likely present in the SE mass estimates for flux-limited, statistical quasar samples, and the distribution of SE masses in finite luminosity bins can be narrower than the nominal uncertainty of these mass estimates. I then discuss the latest applications of SE virial masses in quasar studies, including the early growth of supermassive black holes, quasar demography in the mass-luminosity plane, and the evolution of the BH-host scaling relations, with specific emphases on selection effects and sample biases in the SE masses. I conclude that there is a pressing need to understand and deal with the errors in these BH mass estimates, and to improve these BH weighing methods with substantially more and better reverberation mapping data.

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.

We present the measurement of the two-point cross-correlation function (CCF) of 8,198 Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) quasars and 349,608 DR10 CMASS galaxies from the Baryonic Oscillation Spectroscopic Survey (BOSS) at redshift <z>~0.5 (0.3<z<0.9). The cross-correlation function can be reasonably well fit by a power-law model xi_QG(r)=(r/r_0)^(-gamma) on projected scales of r_p=2-25 Mpc/h with r_0=6.61+-0.25 Mpc/h and gamma=1.69+-0.07. We estimate a quasar linear bias of b_Q=1.38+-0.10 at <z>=0.53 from the CCF measurements. This linear bias corresponds to a characteristic host halo mass of ~4x10^12 M_sun/h, compared to ~10^13 M_sun/h characteristic host halo mass for CMASS galaxies. We divide the quasar sample in luminosity and constrain the luminosity dependence of quasar bias to be db_Q/dlogL=0.20+-0.34 or 0.11+-0.32 (depending on different luminosity divisions) for quasar luminosities -23.5>M_i(z=2)>-25.5, implying a weak luminosity dependence of quasar clustering for the bright end of the quasar population at <z>~0.5. We compare our measurements with theoretical predictions, Halo Occupation Distribution (HOD) models and mock catalogs. These comparisons suggest quasars reside in a broad range of host halos, and the host halo mass distributions significantly overlap with each other for quasars at different luminosities, implying a poor correlation between halo mass and instantaneous quasar luminosity. We also find that the quasar HOD parameterization is largely degenerate such that different HODs can reproduce the CCF equally well, but with different outcomes such as the satellite fraction and host halo mass distribution. These results highlight the limitations and ambiguities in modeling the distribution of quasars with the standard HOD approach and the need for additional information in populating quasars in dark matter halos with HOD. [Abridged]

Kiloparsec-scale binary active galactic nuclei (AGNs) signal active supermassive black hole (SMBH) pairs in merging galaxies. Despite their significance, unambiguously confirmed cases remain scarce and most have been discovered serendipitously. In a previous systematic search, we optically identified four kpc-scale binary AGNs from candidates selected with double-peaked narrow emission lines at redshifts of 0.1--0.2. Here we present Chandra and Hubble Space Telescope Wide Field Camera 3 (WFC3) imaging of these four systems. We critically examine and confirm the binary-AGN scenario for two of the four targets, by combining high angular resolution X-ray imaging spectroscopy with Chandra ACIS-S, better nuclear position constraints from WFC3 F105W imaging, and direct starburst estimates from WFC3 F336W imaging; for the other two targets, the existing data are still consistent with the binary-AGN scenario, but we cannot rule out the possibility of only one AGN ionizing gas in both merging galaxies. We find tentative evidence for a systematically smaller X-ray-to-[O III] luminosity ratio and/or higher Compton-thick fraction in optically selected kpc-scale binary AGNs than in single AGNs, possibly caused by a higher nuclear gas column due to mergers and/or a viewing angle bias related to the double-peak narrow line selection. While our result lends some further support to the general approach of optically identifying kpc-scale binary AGNs, it also highlights the challenge and ambiguity of X-ray confirmation.

We investigate the range of covering factors (determined from the ratio of IR to UV/optical luminosity) seen in luminous type 1 quasars using a combination of data from the WISE, UKIDSS and SDSS surveys. Accretion disk (UV/optical) and obscuring dust (IR) luminosities are measured via the use of a simple three component SED model. We use these estimates to investigate the distribution of covering factors and its relationship to both accretion luminosity and IR SED shape. The distribution of covering factors (f_C) is observed to be log-normal, with a bias-corrected mean of <log10 f_C>=-0.41 and standard deviation of 0.2. The fraction of IR luminosity emitted in the near-IR (1--5 micron) is found to be high (~40 per cent), and strongly dependant on covering factor.

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

We present a new measurement of the optical Quasar Luminosity Function (QLF), using data from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III: BOSS). From the SDSS-III Data Release Nine (DR9), we select a uniform sample of 22,301 i<=21.8 quasars over an area of 2236 sq. deg with confirmed spectroscopic redshifts between 2.2<z<3.5, filling in a key part of the luminosity-redshift plane for optical quasar studies. We derive the completeness of the survey through simulated quasar photometry, and check this completeness estimate using a sample of quasars selected by their photometric variability within the BOSS footprint. We investigate the level of systematics associated with our quasar sample using the simulations, in the process generating color-redshift relations and a new quasar k-correction. We probe the faint end of the QLF to M_i(z=2.2) = -24.5 and see a clear break in the QLF at all redshifts up to z=3.5. We find that a log-linear relation (in log[Phi*] - M*) for a luminosity and density evolution (LEDE) model adequately describes our data within the range 2.2<z<3.5; across this interval the break luminosity increases by a factor of ~2.3 while Phi* declines by a factor of ~6. At z<2.2 our data is reasonably well fit by a pure luminosity evolution (PLE) model. We see only a weak signature of "AGN downsizing", in line with recent studies of the hard X-ray luminosity function. We compare our measured QLF to a number of theoretical models and find that models making a variety of assumptions about quasar triggering and halo occupation can fit our data over a wide range of redshifts and luminosities.

We present the Data Release 9 Quasar (DR9Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the survey, are spectrocopically confirmed as quasars via visual inspection, have luminosities Mi[z=2]<-20.5 (in a $\Lambda$CDM cosmology with H0 = 70 km/s/Mpc, $\Omega_{\rm M}$ = 0.3, and $\Omega_{\Lambda}$ = 0.7) and either display at least one emission line with full width at half maximum (FWHM) larger than 500 km/s or, if not, have interesting/complex absorption features. It includes as well, known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. This catalog contains 87,822 quasars (78,086 are new discoveries) detected over 3,275 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra newly derived from a training set of 8,632 spectra from SDSS-DR7. The number of quasars with $z>2.15$ (61,931) is ~2.8 times larger than the number of z>2.15 quasars previously known. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 7,533 broad absorption line quasars and gives their characteristics. For each object the catalog presents five-band (u,g,r,i,z) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

The coronal blowout jet is a peculiar category among various jet phenomena, of which the sheared base arch, often carrying a small filament, experiences a miniature version of blowout eruption that produces large-scale coronal mass ejection (CME). In this paper, we report such a coronal blowout jet with high-resolution multi-wavelength and multi-angle observations taken from Solar Dynamics Observatory, Solar Terrestrial Relations Observatory, and Big Bear Solar Observatory. For the first time, we find that a simultaneous bubble-like and a jet-like CMEs were dynamically related to the blowout jet that showed cool and hot components next to each other. Our observational results indicate that: (1) the cool component was resulted from the eruption of the filament contained within the jet's base arch, and it further caused the bubble-like CME; (2) the jet-like CME was associated with the hot component, which was the outward moving heated plasmas generated by the reconnection of the base arch and its ambient open field lines. On the other hand, bifurcation of the jet's cool component was also observed, which was resulted from the uncoupling of the erupting filament's two legs that were highly twisted at the very beginning. Based on these results, we propose a model to interpret the coronal blowout jet, of which the external reconnection not only produces the jet-like CME but also leads to the rising of the filament. Subsequently, internal reconnection starts underneath the rising filament and thereby causes the bubble-like CME.

Extreme Ultraviolet (EUV) waves have been found for about 15 years. However, significant controversy remains over their physical natures and origins. In this paper, we report an EUV wave that was accompanied by an X1.9 flare and a partial halo coronal mass ejection. Using high temporal and spatial resolution observations taken by the {\em Solar Dynamics Observatory} and the Solar-TErrestrial RElations Observatory, we are able to investigate the detailed kinematics of the EUV wave. We find several arguments that support the fast-mode wave scenario: (1) The speed of the EUV wave (570 km/s) is higher than the sound speed of quiet-Sun corona. (2) Significant deceleration of the EUV wave (-130 m/s2) is found during its propagation. (3) The EUV wave resulted in the oscillations of a loop and a filament along its propagation path, and a reflected wave from the polar coronal hole is also detected. (4) Refraction or reflection effect is observed when the EUV wave was passing through two coronal bright points. (5) The dimming region behind the wavefront stopped to expand when the wavefront started to become diffuse. (6) The profiles of the wavefront exhibited a dispersive nature, and the magnetosonic Mach number of the EUV wave derived from the highest intensity jump is about 1.4. In addition, triangulation indicates that the EUV wave propagated within a height range of about 60-100 Mm above the photosphere. We propose that the EUV wave observed should be a nonlinear fast-mode magnetosonic wave that propagated freely in the corona after it was driven by the CME expanding flanks during the initial period.

We report two sympathetic solar eruptions, including a partial and a full flux rope eruption in a quadrupolar magnetic region, where a large and a small filament resided above the middle and the east neutral lines respectively. The large filament first rose slowly at a speed of 8 km/s for 23 minutes and then it was accelerated to 102 km/s. Finally, this filament erupted successfully and caused a coronal mass ejection. During the slow rising phase, various evidence for breakout-like external reconnection has been identified at high and low temperature lines. The eruption of the small filament started around the end of the large filament's slow rising. This filament erupted partially and no associating coronal mass ejection could be detected. Based on a potential field extrapolation, we find that the topology of the three-dimensional coronal field above the source region is composed of three low-lying lobes and a large overlying flux system, and a null point located between the middle lobe and the overlying antiparallel flux system. We propose a possible mechanism within the framework of the magnetic breakout model to interpret the sympathetic filament eruptions, in which the magnetic implosion mechanism is thought to be a possible linkage between the sympathetic eruptions, and the external reconnection at the null point transfers field lines from the middle lobe to the lateral lobes and thereby leads to the full (partial) eruption of the observed large (small) filament. Other possible mechanisms are also discussed briefly. We conclude that the structural properties of coronal fields are important for producing sympathetic eruptions.

For the first time, we report a large-scale wave that was observed simultaneously in the photosphere, chromosphere, transition region and low corona layers of the solar atmosphere. Using the high temporal and high spatial resolution observations taken by the Solar Magnetic Activity Research Telescope at Hida Observatory and the Atmospheric Imaging Assembly (AIA) onboard Solar Dynamic Observatory, we find that the wave evolved synchronously at different heights of the solar atmosphere, and it propagated at a speed of 605 km/s and showed a significant deceleration (-424 m/s2) in the extreme-ultraviolet (EUV) observations. During the initial stage, the wave speed in the EUV observations was 1000 km/s, similar to those measured from the AIA 1700 {\AA} (967 km/s) and 1600 {\AA} (893 km/s) observations. The wave was reflected by a remote region with open fields, and a slower wave-like feature at a speed of 220 km/s was also identified following the primary fast wave. In addition, a type-II radio burst was observed to be associated with the wave. We conclude that this wave should be a fast magnetosonic shock wave, which was firstly driven by the associated coronal mass ejection and then propagated freely in the corona. As the shock wave propagated, its legs swept the solar surface and thereby resulted in the wave signatures observed in the lower layers of the solar atmosphere. The slower wave-like structure following the primary wave was probably caused by the reconfiguration of the low coronal magnetic fields, as predicted in the field-line stretching model.

For the future development of Chinese Giant Solar Telescope (CGST) in Western China, a new sky brightness monitor (SBM) has been produced for the site survey for CGST. To critically examine the performance and sensitivity of SBM, we used it in the observation of the annular solar eclipse in Dali City, Yunnan, on 15 January 2010. The observation met good weather condition with almost clear sky during the eclipse. The SBM measurement translates into the solar illuminance changes at a level of 2.4 \times 10-4I s-1 during the eclipse. The time of the minimal sky brightness in the field of view (FOV) is found consistent with the time of maximum eclipse. Two local sky regions in the FOV are chosen to make time series of calibrated skylight profiles. The evolution of the sky brightness thus calibrated also shows good consistency with the eclipse, particularly between the second and the third contacts. The minimal sky brightness in each local sky region took place within half a minute from the corresponding predicted contact time. Such small time delays were mainly caused by occasional cirri. The minimal sky brightness measured during the eclipse is a few millionths of I\odot with standard deviation of 0.11 millionths of I\odot. The observation supports that the singlescattering process (optically thin conditions) is the main contributor to the atmospheric scattering. We have demonstrated that many important aerosol optical parameters can be deduced from our data.We conclude that the new SBM is a sensitive sky photometer that can be used for our CGST and coronagraph site surveys.

On 2011 May 30, quasi-periodic fast propagating (QFP) magnetosonic waves accompanied by a C2.8 flare were directly imaged by the Atomospheric Imaging Assembly instrument on board the Solar Dynamics Observatory. The QFP waves successively emanated from the flare kernel, they propagated along a cluster of open coronal loops with a phase speed of 834 km/s during the flare's rising phase, and the multiple arc-shaped wave trains can be fitted with a series of concentric circles. We generate the k-omega diagram of the Fourier power and find a straight ridge that represents the dispersion relation of the waves. Along the ridge, we find a lot of prominent nodes which represent the available frequencies of the QFP waves. On the other hand, the frequencies of the flare are also obtained by analyzing the flare light curves using the wavelet technique. The results indicate that almost all the main frequencies of the flare are consistent with those of the QFP waves. This suggests that the flare and the QFP waves were possibly excited by a common physical origin. On the other hand, a few low frequencies revealed by the k-omega diagram can not be found in the accompanying flare. We propose that these low frequencies were possibly due to the leakage of the pressure-driven p-mode oscillations from the photosphere into the low corona, which should be a noticeable mechanism for driving the QFP waves observed in the corona.

Spatially extended emission regions of Active Galactic Nuclei (AGN) respond to continuum variations, if such emission regions are powered by energy reprocessing of the continuum. The response from different parts of the reverberating region arrives at different times lagging behind the continuum variation. The lags can be used to map the geometry and kinematics of the emission region (i.e., reverberation mapping, RM). If the extended emission region is not spherically symmetric in configuration and velocity space, reverberation may produce astrometric offsets in the emission region photocenter as a function of time delay and velocity, detectable with future micro-arcsec to tens of micro-arcsec astrometry. Such astrometric responses provide independent constraints on the geometric and kinematic structure of the extended emission region, complementary to traditional reverberation mapping. In addition, astrometric RM is more sensitive to infer the inclination of a flattened geometry and the rotation angle of the extended emission region.

We present 21 examples of C IV Broad Absorption Line (BAL) trough disappearance in 19 quasars selected from systematic multi-epoch observations of 582 bright BAL quasars (1.9 < z < 4.5) by the Sloan Digital Sky Survey-I/II (SDSS-I/II) and SDSS-III. The observations span 1.1-3.9 yr rest-frame timescales, longer than have been sampled in many previous BAL variability studies. On these timescales, ~2.3% of C IV BAL troughs disappear and ~3.3% of BAL quasars show a disappearing trough. These observed frequencies suggest that many C IV BAL absorbers spend on average at most a century along our line of sight to their quasar. Ten of the 19 BAL quasars showing C IV BAL disappearance have apparently transformed from BAL to non-BAL quasars; these are the first reported examples of such transformations. The BAL troughs that disappear tend to be those with small-to-moderate equivalent widths, relatively shallow depths, and high outflow velocities. Other non-disappearing C IV BALs in those nine objects having multiple troughs tend to weaken when one of them disappears, indicating a connection between the disappearing and non-disappearing troughs, even for velocity separations as large as 10000-15000 km/s. We discuss possible origins of this connection including disk-wind rotation and changes in shielding gas.

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

We present an estimate of the projected two-point correlation function (2PCF) of quasars in the Sloan Digital Sky Survey (SDSS) over the full range of one- and two-halo scales, 0.02-120 Mpc/h. This was achieved by combining data from SDSS DR7 on large scales and Hennawi et al. (2006; with appropriate statistical corrections) on small scales. Our combined clustering sample is the largest spectroscopic quasar clustering sample to date, containing ~48,000 quasars in the redshift range 0.4<z<2.5 with median redshift 1.4. We interpret these precise 2PCF measurements within the halo occupation distribution (HOD) framework and constrain the occupation functions of central and satellite quasars in dark matter halos. In order to explain the small-scale clustering, the HOD modeling requires that a small fraction of z~1.4 quasars, fsat=(7.4+/-1.4) 10^(-4), be satellites in dark matter halos. At z~1.4, the median masses of the host halos of central and satellite quasars are constrained to be Mcen=(4.1+0.3/-0.4) 10^12 Msun/h and Msat=(3.6+0.8/-1.0) 10^14 Msun/h, respectively. To investigate the redshift evolution of the quasar-halo relationship, we also perform HOD modeling of the projected 2PCF measured by Shen et al. (2007) for SDSS quasars with median redshift 3.2. We find tentative evidence for an increase in the mass scale of quasar host halos---the inferred median mass of halos hosting central quasars at z~3.2 is Mcen=(14.1+5.8/-6.9) 10^12 Msun/h. The cutoff profiles of the mean occupation functions of central quasars reveal that quasar luminosity is more tightly correlated with halo mass at higher redshifts. The average quasar duty cycle around the median host halo mass is inferred to be fq=(7.3+0.6/-1.5) 10^(-4) at z~1.4 and fq=(8.6+20.4/-7.2) 10^(-2) at z~3.2. We discuss the implications of our results for quasar evolution and quasar-galaxy co-evolution.

We measure the quasar two-point correlation function over the redshift range 2.2<z<2.8 using data from the Baryon Oscillation Spectroscopic Survey. We use a homogeneous subset of the data consisting of 27,129 quasars with spectroscopic redshifts---by far the largest such sample used for clustering measurements at these redshifts to date. The sample covers 3,600 square degrees, corresponding to a comoving volume of 9.7(Gpc/h)^3 assuming a fiducial LambdaCDM cosmology, and it has a median absolute i-band magnitude of -26, k-corrected to z=2. After accounting for redshift errors we find that the redshift space correlation function is fit well by a power-law of slope -2 and amplitude s_0=(9.7\pm 0.5)Mpc/h over the range 3<s<25Mpc/h. The projected correlation function, which integrates out the effects of peculiar velocities and redshift errors, is fit well by a power-law of slope -1 and r_0=(8.4\pm 0.6)Mpc/h over the range 4<R<16Mpc/h. There is no evidence for strong luminosity or redshift dependence to the clustering amplitude, in part because of the limited dynamic range in our sample. Our results are consistent with, but more precise than, previous measurements at similar redshifts. Our measurement of the quasar clustering amplitude implies a bias factor of b~3.5 for our quasar sample. We compare the data to models to constrain the manner in which quasars occupy dark matter halos at z~2.4 and infer that such quasars inhabit halos with a characteristic mass of <M>~10^{12}Msun/h with a duty cycle for the quasar activity of 1 per cent.

Single-epoch virial black hole (BH) mass estimators utilizing broad emission lines have been routinely applied to high-redshift quasars to estimate their BH masses. Depending on the redshift, different line estimators (Halpha, Hbeta, MgII, CIV) are often used with optical/near-infrared spectroscopy. Here we use a homogeneous sample of 60 intermediate-redshift (z~1.5-2.2) SDSS quasars with optical and near-infrared spectra covering CIV through Halpha to investigate the consistency between different line estimators. We critically compare restframe UV line estimators (CIV, CIII], and MgII) with optical estimators (Hbeta and Halpha) in terms of correlations between line widths and between continuum/line luminosities, for the high-luminosity regime (L_5100>10^45.4 erg/s) probed by our sample. The continuum luminosities of L_1350 and L_3000, and the broad line luminosities are well correlated with L_5100. We found that the MgII FWHM correlates well with the FWHMs of the Balmer lines, and that the MgII line estimator can be calibrated to yield consistent virial mass estimates with those based on the Hbeta/Halpha estimators, thus extending earlier results on less luminous objects. The CIV FWHM is poorly correlated with the Balmer line FWHMs, and the scatter between the CIV and Hbeta FWHMs consists of an irreducible part (~0.12 dex), and a part that correlates with the blueshift of the CIV centroid relative to that of Hbeta. The CIII] FWHM is found to correlate with the CIV FWHM, and hence is also poorly correlated with the Hbeta FWHM. While the CIV and CIII] lines can be calibrated to yield consistent virial mass estimates as Hbeta on average, the scatter is substantially larger than MgII, and the usage of CIV/CIII] FWHM in the mass estimators does not improve the agreement with the Hbeta estimator. (Abridged)

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.

A few percent of quasars show strong associated MgII absorption, with velocities (v_off) lying within a few thousand km/s from the quasar systemic redshift. These associated absorption line systems (AALs) are usually interpreted as absorbers that are either intrinsic to the quasar and its host, or arising from external galaxies clustering around the quasar. Using composite spectra of ~1,800 MgII AAL quasars selected from SDSS DR7 at 0.4<~ z<~2, we show that quasars with AALs with v_off<1500 km/s have a prominent excess in [OII]3727 emission (detected at >7sigma) at rest to the quasar host, compared to unabsorbed quasars. We interpret this [OII] excess as due to enhanced star formation in the quasar host. Our results suggest that a significant fraction of AALs with v_off<1500 km/s are physically associated with the quasar and its host, and are most likely large-scale flows with typical velocities of a few hundred km/s. AAL quasars also have dust reddening lying between normal quasars and the so-called dust-reddened quasars. We suggest that the unique properties of AAL quasars can be explained if they are the transitional population from heavily dust-reddened quasar to normal quasars in the formation process of quasars and their hosts. This scenario predicts a larger fraction of young bulges, disturbed morphologies and interactions of AAL quasar hosts compared to normal quasars. The intrinsic link between associated absorbers and quasar hosts opens a new window to probe massive galaxy formation and galactic-scale feedback processes, and provides a crucial test of the evolutionary picture of quasars.

Active galactic nuclei (AGNs) are occasionally seen in pairs, suggesting that tidal encounters are responsible for the accretion of material by both central supermassive black holes (BHs). In Paper I of this series, we selected a sample of AGN pairs with projected separations r_p < 100 kpc and velocity offsets < 600 km/s from the SDSS DR7 and quantified their frequency. In this paper, we address the BH accretion and recent star formation properties in their host galaxies. AGN pairs experience stronger BH accretion, as measured by their [O III]5007 luminosities (corrected for contribution from star formation) and Eddington ratios, than do control samples of single AGNs matched in redshift and host stellar mass. Their host galaxies have stronger post-starburst activity and younger mean stellar ages, as indicated by stronger H-delta absorption and smaller 4000 A break in their spectra. The BH accretion and recent star formation in the host galaxies both increase with decreasing projected separation in AGN pairs, for r_p ~< 10-30 kpc. The intensity of BH accretion, the post-starburst strength, and the mean stellar ages are correlated between the two AGNs in a pair. The luminosities and Eddington ratios of AGN pairs are correlated with recent star formation in their host galaxies, with a scaling relation consistent with that observed in single AGNs. Our results suggest that galaxy tidal interactions enhance both BH accretion and star formation in AGN pairs, even though the majority of low redshift AGNs is not coincident with on-going interactions.

Galaxy-galaxy mergers and close interactions have long been regarded as a viable mechanism for channeling gas toward the central supermassive black holes (SMBHs) of galaxies which are triggered as active galactic nuclei (AGNs). AGN pairs, in which the central SMBHs of a galaxy merger are both active, are expected to be common from such events. We conduct a systematic study of 1286 AGN pairs at z ~ 0.1 with line-of-sight velocity offsets dv < 600 km/s and projected separations r_p < 100 kpc, selected from the Seventh Data Release of the Sloan Digital Sky Survey (SDSS). This AGN pair sample was drawn from 138,070 AGNs optically identified based on diagnostic emission line ratios and/or line widths. The fraction of AGN pairs with 5 kpc ~< r_p < 100 kpc among all spectroscopically selected AGNs at 0.02 < z < 0.16 is 3.6% after correcting for SDSS spectroscopic incompleteness; ~ 30% of these pairs show morphological tidal features in their SDSS images, and the fraction becomes >~ 80% for pairs with the brightest nuclei. Our sample increases the number of known AGN pairs on these scales by more than an order of magnitude. We present their AGN and host-galaxy star formation properties in a companion paper.

We present an observational study of the kinematics and fine structure of an unwinding polar jet, with high temporal and spatial observations taken by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO) and the Solar Magnetic Activity Research Telescope (SMART). During the rising period, the shape of the jet resembled a cylinder with helical structures on the surface, while the mass of the jet was mainly distributed on the cylinder's shell. In the radial direction, the jet expanded successively at its western side. The radial expansion presented three distinct phases: the gradually expanding phase, the fast expanding phase, and the steady phase. Each phase lasted for about 12 minutes. The angular speed of the unwinding jet and the twist transferred into the outer corona during the eruption are estimated to be 11.1 \times 10{-3} rad/s (period = 564 s) and 1.17 to 2.55 turns (or 2.34 to 5.1{\pi}) respectively. On the other hand, by calculating the azimuthal component of the magnetic field in the jet and comparing the free energy stored in the non-potential magnetic field with the jet's total energy, we find that the non-potential magnetic field in the jet is enough to supply the energy for the ejection. These new observational results strongly support the scenario that the jets are driven by the magnetic twist, which is stored in the twisted closed field of a bipole, and released through magnetic reconnection between the bipole and its ambient open field.

Hierarchical galaxy mergers will lead to the formation of binary and, in the case of a subsequent merger before a binary coalesce, triple supermassive black holes (SMBHs), given that most massive galaxies harbor SMBHs. A triple of SMBHs becomes visible as a triple active galactic nucleus (AGN) when the BHs accrete large amounts of gas at the same time. Here we report the discovery of a kpc-scale triple AGN, SDSSJ1027+1749 at z = 0.066, from our systematic search for hierarchical mergers of AGNs. The galaxy contains three emission-line nuclei, two of which are offset by ~ 450 and 110 km/s in velocity and by 2.4 and 3.0 kpc in projected separation from the central nucleus. All three nuclei are classified as obscured AGNs based on optical diagnostic emission line ratios, with black hole mass estimates M_BH ~> 10^8 M_sun from stellar velocity dispersions measured in the associated stellar components. Based on dynamical friction timescale estimates, the three stellar components in SDSSJ1027+1749 will merge in ~ 40 Myr, and their associated SMBHs may evolve into a gravitationally interacting triple system in ~< 200 Myr. Our result sets a lower limit of ~ 5 x 10^(-5) for the fraction of kpc-scale triples in optically selected AGNs at z ~ 0.1.

We present stereoscopic observations of six sequent eruptions of a filament in the active region NOAA 11045 on 2010 Feb 8, with the advantage of the STEREO twin viewpoints in combination with the earth viewpoint from SOHO instruments and ground-based telescopes. The last one of the six eruptions is with a coronal mass ejection, while the others are not. The flare in this successful one is more intensive than in the others. Moreover, the filament material velocity of the successful one is also the largest among them. Interestingly, all the filament velocities are found proportional to their flare powers. We calculate magnetic field intensity at low altitude, the decay indexes of the external field above the filament, and the asymmetry properties of the overlying fields before and after the failed eruptions and find little difference between them, indicating the same coronal confinement for the failed and the successful eruptions. The results suggest that, besides the confinement of coronal magnetic field, the energy released in low corona should be another crucial element for production of a failed or successful filament eruption. That is, only a coronal mass ejection can be launched away if the energy released exceeds some critical value, given the same coronal conditions.

(Abridged) Approximately 1% of low redshift (z<0.3) optically-selected type 2 AGNs show a double-peaked [OIII] narrow emission line profile in their spatially-integrated spectra. Such features are usually interpreted as due either to kinematics, such as biconical outflows and/or disk rotation of the narrow line region (NLR) around single black holes, or to the relative motion of two distinct NLRs in a merging pair of AGNs. Here we report follow-up near infrared (NIR) imaging and optical slit spectroscopy of 31 double-peaked [OIII] type 2 AGNs drawn from the SDSS parent sample presented in Liu et al (2010). These data reveal a mixture of origins for the double-peaked feature. Roughly 10% of our objects are best explained by binary AGNs at (projected) kpc-scale separations, where two stellar components with spatially coincident NLRs are seen. ~ 50% of our objects have [OIII] emission offset by a few kpc, corresponding to the two velocity components seen in the SDSS spectra, but there are no corresponding double stellar components seen in the NIR imaging. For those objects with sufficiently high quality slit spectra, we see velocity and/or velocity dispersion gradients in [OIII] emission, suggestive of the kinematic signatures of a single NLR. The remaining ~40% of our objects are ambiguous, and will need higher spatial resolution observations to distinguish between the two scenarios. Our observations therefore favor the kinematics scenario with a single AGN for the majority of these double-peaked [OIII] type 2 AGNs. We emphasize the importance of combining imaging and slit spectroscopy in identifying kpc binary AGNs, i.e., in no cases does one of these alone allow an unambiguous identification. We estimate that ~ 0.5-2.5% of the z<0.3 type 2 AGNs are kpc-scale binary AGNs of comparable luminosities, with a relative orbital velocity >~150 km/s.

We present a compilation of properties of the 105,783 quasars in the SDSS Data Release 7 (DR7) quasar catalog. In this value-added product, we compile continuum and emission line measurements around the Halpha, Hbeta, MgII and CIV regions, as well as other quantities such as radio properties, broad absorption line quasar (BALQSO) flags, and disk emitters. We also compile virial black hole mass estimates based on various calibrations. For the fiducial virial mass estimates we use the Vestergaard & Peterson (VP06) calibrations for Hbeta and CIV, and our own calibration for MgII which matches the VP06 Hbeta masses on average. We describe the construction of this catalog, and discuss its limitations. The catalog and its future updates will be made publicly available online.

Using a sample of 30,000 quasars from SDSS-DR7, we explore the range of properties exhibited by high-ionization, broad emission lines, such as CIV 1549. Specifically we investigate the anti-correlation between L_UV and emission line EQW (the Baldwin Effect) and the "blueshifting" of high-ionization emission lines. The blueshift of the CIV emission line is nearly ubiquitous, with a mean shift of 810 km/s for radio-quiet (RQ) quasars and 360 km/s for radio-loud (RL) quasars, and the Baldwin Effect is present in both RQ and RL samples. Composite spectra are constructed as a function of CIV emission line properties in attempt to reveal empirical relationships between different line species and the SED. Within a two-component disk+wind model of the broad emission line region (BELR), where the wind filters the continuum seen by the disk component, we find that RL quasars are consistent with being dominated by the disk component, while BALQSOs are consistent with being dominated by the wind component. Some RQ objects have emission line features similar to RL quasars; they may simply have insufficient black hole (BH) spin to form radio jets. Our results suggest that there could be significant systematic errors in the determination of L_bol and BH mass that make it difficult to place these findings in a more physical context. However, it is possible to classify quasars in a paradigm where the diversity of BELR parameters are due to differences in an accretion disk wind between quasars (and over time); these differences are underlain primarily by the SED, which ultimately must be tied to BH mass and accretion rate.

Double-peaked broad emission lines in Active Galactic Nuclei (AGNs) may indicate the existence of a bound supermassive black hole (SMBH) binary where two distinct broad line regions (BLRs) contribute together to the line profile. An alternative interpretation is a disk emitter origin for the double-peaked line profile. Using simple BLR models, we calculate the expected broad line profile for a SMBH binary at different separations. Under reasonable assumptions that both BLRs are illuminated by the two active SMBHs and that the ionizing flux at the BLR location is roughly constant, we confirm the emergence of double-peaked features and radial velocity drifts of the two peaks due to the binary orbital motion. However, such a clear double-peaked feature only arises in a particular stage of the binary evolution when the two BHs are close enough such that the line-of-sight orbital velocity difference is larger than the FWHM of the individual broad components, while the two BLRs are still mostly distinct. Prior to this stage, the velocity splitting due to the orbit motion of the binary is too small to separate the emission from the two BLRs, leading to asymmetric broad line profiles in general. When the two BHs are even closer such that the two BLRs can no longer be distinct, the line profile becomes more complex and the splitting of the peaks does not correspond to the orbital motion of the binary. In this regime there are no coherent radial velocity drifts in the peaks with time. Asymmetric line profiles are probably a far more common signature of binary SMBHs than are double-peaked profiles. We discuss the temporal variations of the broad line profile for binary SMBHs and highlight the different behaviors of reverberation mapping in the binary and disk emitter cases, which may serve as a feasible tool to disentangle these two scenarios.

Extended extragalactic radio sources have traditionally been classified into FR I and II types, based on the ratio r of the separation S between the brightest regions on either sides of the host galaxy and the total size T of the radio source. Here we examine the distribution of various physical properties as a function of r of 1040 luminous radio galaxies (RGs) at z<0.3 from the SDSS, NVSS, and FIRST. About 2/3 of the RGs are lobe-dominated (LD), and 1/3 have prominent jets. If we follow the original definition of the FR types (a division based solely on r), FR I and FR II RGs overlap in their host galaxy properties. However, the rare, LD sources with r>0.8 AND OIII5007 line luminosity >10^6 Lsun are markedly different on average from the rest of the RGs, in the sense that they are hosted in lower mass galaxies, live in relatively sparse environments, and have higher accretion rates onto the central SMBH. Thus these objects and the rest of RGs form a well-defined dichotomy. Motivated by the stark differences in the nuclear emission line properties of the RG subsamples, we suggest that the accretion rate onto the SMBH may play the primary role in creating the different morphologies. At relatively high accretion rates, the accretion system may produce powerful jets that create the "classical double" morphology (roughly corresponding to the LD sources with r>0.8 and emission lines); at lower accretion rates the jets from a radiatively inefficient accretion flow generate radio lobes without apparent "hot spots" at the edge (corresponding to the majority of LD sources). At slightly lower accretion rates AND in galaxies with dense galactic structure, sources with prominent jets result. It is possible that while the high accretion rate systems could affect sub-Mpc scale environments, the jets from lower accretion rate systems may efficiently suppress activity within the host galaxies.

We report the discovery of four kpc-scale binary AGNs. These objects were originally selected from the Sloan Digital Sky Survey based on double-peaked [O III] 4959,5007 emission lines in their fiber spectra. The double peaks could result from pairing active supermassive black holes (SMBHs) in a galaxy merger, or could be due to bulk motions of narrow-line region gas around a single SMBH. Deep near-infrared (NIR) images and optical slit spectra obtained from the Magellan 6.5 m and the APO 3.5 m telescopes strongly support the binary SMBH scenario for the four objects. In each system, the NIR images reveal tidal features and double stellar components with a projected separation of several kpc, while optical slit spectra show two Seyfert 2 nuclei spatially coincident with the stellar components, with line-of-sight velocity offsets of a few hundred km/s. These objects were drawn from a sample of only 43 objects, demonstrating the efficiency of this technique to find kpc-scale binary AGNs.

We present the fifth edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog, which is based upon the SDSS Seventh Data Release. The catalog, which contains 105,783 spectroscopically confirmed quasars, represents the conclusion of the SDSS-I and SDSS-II quasar survey. The catalog consists of the SDSS objects that have luminosities larger than M_i = -22.0 (in a cosmology with H_0 = 70 km/s/Mpc Omega_M = 0.3, and Omega_Lambda = 0.7) have at least one emission line with FWHM larger than 1000 km/s or have interesting/complex absorption features, are fainter than i > 15.0 and have highly reliable redshifts. The catalog covers an area of 9380 deg^2. The quasar redshifts range from 0.065 to 5.46, with a median value of 1.49; the catalog includes 1248 quasars at redshifts greater than four, of which 56 are at redshifts greater than five. The catalog contains 9210 quasars with i < 18; slightly over half of the entries have i< 19. For each object the catalog presents positions accurate to better than 0.1" rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800-9200 Ang. at a spectral resolution R = 2000 the spectra can be retrieved from the SDSS public database using the information provided in the catalog. Over 96% of the objects in the catalog were discovered by the SDSS. We also include a supplemental list of an additional 207 quasars with SDSS spectra whose archive photometric information is incomplete.

We use Shen et al.'s (2009) measurements of luminosity-dependent clustering in the SDSS Data Release 5 Quasar Catalog, at redshifts 0.4 < z < 2.5, to constrain the relation between quasar luminosity and host halo mass and to infer the duty cycle f_opt, the fraction of black holes that shine as optically luminous quasars at a given time. We assume a monotonic mean relation between quasar luminosity and host halo mass, with log-normal scatter \Sigma. For specified f_opt and \Sigma, matching the observed quasar space density determines the normalization of the luminosity-halo mass relation, from which we predict the clustering bias. The data show no change of bias between the faint and bright halves of the quasar sample but a modest increase in bias for the brightest 10%. At the mean redshift z=1.45 of the sample, the data can be well described either by models with small intrinsic scatter (\Sigma=0.1 dex) and a duty cycle f_opt=6*10^(-4) or by models with much larger duty cycles and larger values of the scatter. "Continuity equation" models of the black hole mass population imply f_opt > 2*10^(-3) in this range of masses and redshifts, and the combination of this constraint with the clustering measurements implies scatter \Sigma > 0.4 dex. These findings contrast with those inferred from the much stronger clustering of high-luminosity quasars at z~4, which require minimal scatter between luminosity and halo mass and duty cycles close to one.

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.

The most distant quasars known, at redshifts z=6, generally have properties indistinguishable from those of lower-redshift quasars in the rest-frame ultraviolet/optical and X-ray bands. This puzzling result suggests that these distant quasars are evolved objects even though the Universe was only seven per cent of its current age at these redshifts. Recently one z=6 quasar was shown not to have any detectable emission from hot dust, but it was unclear whether that indicated different hot-dust properties at high redshift or if it is simply an outlier. Here we report the discovery of a second quasar without hot-dust emission in a sample of 21 z=6 quasars. Such apparently hot-dust-free quasars have no counterparts at low redshift. Moreover, we demonstrate that the hot-dust abundance in the 21 quasars builds up in tandem with the growth of the central black hole, whereas at low redshift it is almost independent of the black hole mass. Thus z=6 quasars are indeed at an early evolutionary stage, with rapid mass accretion and dust formation. The two hot-dust-free quasars are likely to be first-generation quasars born in dust-free environments and are too young to have formed a detectable amount of hot dust around them.

Type 2 quasars are luminous active galactic nuclei (AGN) whose central regions are obscured by large amounts of gas and dust. In this paper, we present a catalog of type 2 quasars from the Sloan Digital Sky Survey (SDSS), selected based on their optical emission lines. The catalog contains 887 objects with redshifts z < 0.83; this is six times larger than the previous version and is by far the largest sample of type 2 quasars in the literature. We derive the [OIII]5008 luminosity function for 10^8.3 Lsun < L[OIII] < 10^10 Lsun (corresponding to intrinsic luminosities up to M[2400A]-28 mag or bolometric luminosities up to 4x10^47 erg/sec). This luminosity function provides strong lower limits to the actual space density of obscured quasars, due to our selection criteria, the details of the spectroscopic target selection, as well as other effects. We derive the equivalent luminosity function for the complete sample of type 1 (unobscured) quasars; then, we determine the ratio of type 2/type 1 quasar number densities. Our best data constrain this ratio to be at least 1.5:1 for 10^8.3 Lsun < L[OIII] < 10^9.5 Lsun at z < 0.3, and at least 1.2:1 for L[OIII]=10^10 Lsun at 0.3 < z < 0.83. Type 2 quasars are at least as abundant as type 1 quasars in the relatively nearby Universe (z < 0.8) for the highest luminosities.

We present a sample of 167 type 2 AGNs with double-peaked [O III] 4959,5007 narrow emission lines, selected from the Seventh Data Release of the Sloan Digital Sky Survey. The double-peaked profiles can be well modeled by two velocity components, blueshifted and redshifted from the systemic velocity. Half of these objects have a more prominent redshifted component. In cases where the H-beta emission line is strong, it also shows two velocity components whose line-of-sight (LOS) velocity offsets are consistent with those of [O III]. The relative LOS velocity offset between the two components is typically a few hundred km/s, larger by a factor of ~ 1.5 than the line full width at half maximum of each component. The offset correlates with the host stellar velocity dispersion sigma*. The host galaxies of this sample show systematically larger sigma*, stellar masses, and concentrations, and older luminosity-weighted mean stellar ages than a regular type 2 AGN sample matched in redshift, [O III] 5007 equivalent width and luminosity; they show no significant difference in radio properties. These double-peaked features could be due to narrow-line region kinematics, or binary black holes. The statistical properties do not show strong preference for or against either scenario, and spatially resolved optical imaging, spectroscopy, radio or X-ray followup are needed to draw firm conclusions.

We present measurements of the small-scale (0.1<~ r <~ 1 Mpc/h) quasar two-point correlation function at z>2.9, for a flux-limited (i<21) sample of 15 binary quasars compiled by Hennawi et al. (2009). The amplitude of the small-scale clustering increases from z ~ 3 to z ~ 4. The small-scale clustering amplitude is comparable to or lower than power-law extrapolations (with slope gamma=2) from the large-scale correlation function of the i<20.2 quasar sample from the Sloan Digital Sky Survey. Using simple prescriptions relating quasars to dark matter halos, we model the observed small-scale clustering with halo occupation models. Reproducing the large-scale clustering amplitude requires that the active fraction of the black holes in the central galaxies of halos is near unity, but the level of small-scale clustering favors an active fraction of black holes in satellite galaxies 0.1 <~ f_s <~ 0.5 at z >~ 3.

The clustering of quasars on small scales yields fundamental constraints on models of quasar evolution and the buildup of supermassive black holes. This paper describes the first systematic survey to discover high redshift binary quasars. Using color-selection and photometric redshift techniques, we searched 8142 deg^2 of SDSS imaging data for binary quasar candidates, and confirmed them with follow-up spectroscopy. Our sample of 27 high redshift binaries (24 of them new discoveries) at redshifts 2.9 < z < 4.3 with proper transverse separations 10 kpc < R_{\perp} < 650 kpc increases the number of such objects known by an order of magnitude. Eight members of this sample are very close pairs with R_{\perp} < 100 kpc, and of these close systems four are at z > 3.5. The completeness and efficiency of our well-defined selection algorithm are quantified using simulated photometry and we find that our sample is ~ 50% complete. Our companion paper uses this knowledge to make the first measurement of the small scale clustering (R < 1 Mpc/h comoving) of high-redshift quasars. High redshift binaries constitute exponentially rare coincidences of two extreme (M >~ 10^9 Msun) supermassive black holes. At z ~ 4 there is about one close binary per 10 Gpc^3, thus these could be the highest sigma peaks, the analogs of superclusters, in the early Universe.

(Abridged) We present a simple framework for the growth and evolution of supermassive black holes (SMBHs) in the hierarchical structure formation paradigm. In our model, black hole accretion is triggered during major mergers (mass ratio>~0.3) between host dark matter halos. The successive evolution of quasar luminosities follows a universal light curve form: an initial exponential growth at constant Eddington ratio of order unity until it reaches the peak luminosity, followed by a power-law decay. Assuming that the peak luminosity correlates with the post-merger halo mass, we convolve the light curve with the triggering rate of quasar activity to predict the quasar luminosity function (LF). Our model reproduces the observed LF at 0.5<z<4.5 for the full luminosity ranges probed by current optical and X-ray surveys. Our model reproduces both the observed redshift evolution and luminosity dependence of the linear bias of quasar/AGN clustering. Due to the scatter between instantaneous luminosity and halo mass, quasar/AGN clustering weakly depends on luminosity at low to intermediate luminosities; but the linear bias rises rapidly with luminosity at the high luminosity end and at high redshift. In our model, the Eddington ratio distribution is roughly log-normal, which broadens and shifts to lower mean values from high luminosity quasars (L_bol>~10^46 erg s^-1) to low luminosity AGNs (L_bol<~ 10^45 erg s^-1), in good agreement with observations.

We identify a sample of 74 high-redshift quasars (z>3) with weak emission lines from the Fifth Data Release of the Sloan Digital Sky Survey and present infrared, optical, and radio observations of a subsample of four objects at z>4. These weak emission-line quasars (WLQs) constitute a prominent tail of the Lya+NV equivalent width distribution, and we compare them to quasars with more typical emission-line properties and to low-redshift active galactic nuclei with weak/absent emission lines, namely BL Lac objects. We find that WLQs exhibit hot (T~1000 K) thermal dust emission and have rest-frame 0.1-5 micron spectral energy distributions that are quite similar to those of normal quasars. The variability, polarization, and radio properties of WLQs are also different from those of BL Lacs, making continuum boosting by a relativistic jet an unlikely physical interpretation. The most probable scenario for WLQs involves broad-line region properties that are physically distinct from those of normal quasars.

We present measurements of the quasar two-point correlation function, \xi_{Q}, over the redshift range z=0.3-2.2 based upon data from the SDSS. Using a homogeneous sample of 30,239 quasars with spectroscopic redshifts from the DR5 Quasar Catalogue, our study represents the largest sample used for this type of investigation to date. With this redshift range and an areal coverage of approx 4,000 deg^2, we sample over 25 h^-3 Gpc^3 (comoving) assuming the current LCDM cosmology. Over this redshift range, we find that the redshift-space correlation function, xi(s), is adequately fit by a single power-law, with s_{0}=5.95+/-0.45 h^-1 Mpc and \gamma_{s}=1.16+0.11-0.16 when fit over s=1-25 h^-1 Mpc. Using the projected correlation function we calculate the real-space correlation length, r_{0}=5.45+0.35-0.45 h^-1 Mpc and \gamma=1.90+0.04-0.03, over scales of rp=1-130 h^-1 Mpc. Dividing the sample into redshift slices, we find very little, if any, evidence for the evolution of quasar clustering, with the redshift-space correlation length staying roughly constant at s_{0} ~ 6-7 h^-1 Mpc at z<2.2 (and only increasing at redshifts greater than this). Comparing our clustering measurements to those reported for X-ray selected AGN at z=0.5-1, we find reasonable agreement in some cases but significantly lower correlation lengths in others. We find that the linear bias evolves from b~1.4 at z=0.5 to b~3 at z=2.2, with b(z=1.27)=2.06+/-0.03 for the full sample. We compare our data to analytical models and infer that quasars inhabit dark matter haloes of constant mass M ~2 x 10^12 h^-1 M_Sol from redshifts z~2.5 (the peak of quasar activity) to z~0. [ABRIDGED]

We study the light scattering properties of random ballistic aggregates constructed in Shen et al. (Paper I). Using the discrete-dipole-approximation, we compute the scattering phase function and linear polarization for random aggregates with various sizes and porosities, and with two different compositions: 100% silicate and 50% silicate-50% graphite. We investigate the dependence of light scattering properties on wavelength, cluster size and porosity using these aggregate models. We find that while the shape of the phase function depends mainly on the size parameter of the aggregates, the linear polarization depends on both the size parameter and the porosity of the aggregates, with increasing degree of polarization as the porosity increases. Contrary to previous studies, we argue that monomer size has negligible effects on the light scattering properties of ballistic aggregates, as long as the constituent monomer is smaller than the incident wavelength up to 2*pi*a_0/lambda\sim 1.6 where a_0 is the monomer radius. Previous claims for such monomer size effects are in fact the combined effects of size parameter and porosity. Finally, we present aggregate models that can reproduce the phase function and polarization of scattered light from the AU Mic debris disk and from cometary dust, including the negative polarization observed for comets at scattering angles 160<theta<180 deg. These aggregates have moderate porosities, P\sim 0.6, and are of sub-micron-size for the debris disk case, or micron-size for the comet case.

Motivated by the holographic principle, it has been suggested that the dark energy density may be inversely proportional to the area of the event horizon of the Universe. However, such a model would have a causality problem. In this paper, we propose to replace the future event horizon area with the inverse of the Ricci scalar curvature. We show that this model does not only avoid the causality problem and is phenomenologically viable, but also naturally solves the coincidence problem of dark energy. Our analysis of the evolution of density perturbations show that the matter power spectra and CMB temperature anisotropy is only slightly affected by such modification.

Using a homogenous sample of 38,208 quasars with a sky coverage of $4000 {\rm deg^2}$ drawn from the SDSS Data Release Five quasar catalog, we study the dependence of quasar clustering on luminosity, virial black hole mass, quasar color, and radio loudness. At $z<2.5$, quasar clustering depends weakly on luminosity and virial black hole mass, with typical uncertainty levels $\sim 10%$ for the measured correlation lengths. These weak dependences are consistent with models in which substantial scatter between quasar luminosity, virial black hole mass and the host dark matter halo mass has diluted any clustering difference, where halo mass is assumed to be the relevant quantity that best correlates with clustering strength. However, the most luminous and most massive quasars are more strongly clustered (at the $\sim 2\sigma$ level) than the remainder of the sample, which we attribute to the rapid increase of the bias factor at the high-mass end of host halos. We do not observe a strong dependence of clustering strength on quasar colors within our sample. On the other hand, radio-loud quasars are more strongly clustered than are radio-quiet quasars matched in redshift and optical luminosity (or virial black hole mass), consistent with local observations of radio galaxies and radio-loud type 2 AGN. Thus radio-loud quasars reside in more massive and denser environments in the biased halo clustering picture. Using the Sheth et al.(2001) formula for the linear halo bias, the estimated host halo mass for radio-loud quasars is $\sim 10^{13} h^{-1}M_\odot$, compared to $\sim 2\times 10^{12} h^{-1}M_\odot$ for radio-quiet quasar hosts at $z\sim 1.5$.

We present a catalog of 5039 broad absorption line (BAL) quasars (QSOs) in the Sloan Digital Sky Survey (SDSS) Data Release 5 (DR5) QSO catalog that have absorption troughs covering a continuous velocity range >= 2000 km/s. We have fit ultraviolet (UV) continua and line emission in each case, enabling us to report common diagnostics of BAL strengths and velocities in the range -25,000 to 0 km/s for SiIV $\lambda$1400, CIV $\lambda$1549, AlIII $\lambda$1857, and MgII $\lambda$2799. We calculate these diagnostics using the spectrum listed in the DR5 QSO catalog, and also for spectra from additional SDSS observing epochs when available. In cases where BAL QSOs have been observed with Chandra or XMM-Newton, we report the X-ray monochromatic luminosities of these sources. We confirm and extend previous findings that BAL QSOs are more strongly reddened in the rest-frame UV than non-BAL QSOs and that BAL QSOs are relatively X-ray weak compared to non-BAL QSOs. The observed BAL fraction is dependent on the spectral signal-to-noise (S/N); for higher-S/N sources, we find an observed BAL fraction of approximately 15%. BAL QSOs show a similar Baldwin effect as for non-BAL QSOs, in that their CIV emission equivalent widths decrease with increasing continuum luminosity. However, BAL QSOs have weaker CIV emission in general than do non-BAL QSOs. Sources with higher UV luminosities are more likely to have higher-velocity outflows, and the BAL outflow velocity and UV absorption strength are correlated with relative X-ray weakness. These results are in qualitative agreement with models that depend on strong X-ray absorption to shield the outflow from over-ionization and enable radiative acceleration....

We conduct a systematic survey of the regions in which distant satellites can orbit stably around the four giant planets in the solar system, using orbital integrations of up to $10^9$ yr. In contrast to previous investigations, we use a grid of initial conditions on a surface of section to explore phase space uniformly inside and outside the planet's Hill sphere (radius $r_{\rm H}$; satellites outside the Hill sphere sometimes are also known as quasi-satellites). Our confirmations and extensions of old results and new findings include the following: (i) many prograde and retrograde satellites can survive out to radii $\sim 0.5r_{\rm H}$ and $\sim 0.7r_{\rm H}$, respectively, while some coplanar retrograde satellites of Jupiter and Neptune can survive out to $\sim r_{\rm H}$; (ii) stable orbits do not exist within the Hill sphere at high ecliptic inclinations when the semi-major axis is large enough that the solar tide is the dominant non-Keplerian perturbation; (iii) there is a gap between $\sim r_{\rm H}$ and $2r_{\rm H}$ in which no stable orbits exist; (iv) at distances $\gtrsim 2r_{\rm H}$ stable satellite orbits exist around Jupiter, Uranus and Neptune (but not Saturn). For Uranus and Neptune, in particular, stable orbits are found at distances as large as $\sim 10r_{\rm H}$; (v) the differences in the stable zones beyond the Hill sphere arise mainly from differences in the planet/Sun mass ratio and perturbations from other planets; in particular, the absence of stable satellites around Saturn is mainly due to perturbations from Jupiter. It is therefore likely that satellites at distances $\gtrsim 2r_{\rm H}$ could survive for the lifetime of the solar system around Uranus, Neptune, and perhaps Jupiter.

We investigate the scattering and absorption of light by random ballistic aggregates of spherical monomers. We present a general measure for the porosity of an irregular particle. Three different classes of ballistic aggregates are considered, with different degrees of porosity. Scattering and absorption cross sections are calculated, using the discrete dipole approximation (DDA), for grains of three compositions (50% silicate and 50% graphite; 50% silicate and 50% amorphous carbon; and 100% silicate), for wavelengths from 0.1 micron to 4 micron. For fixed particle mass, increased porosity increases the extinction at short wavelengths, but decreases the extinction at wavelengths long compared to the overall aggregate size. Scattering and absorption cross sections are insensitive to monomer size as long as the constituent monomers are small compared with the incident wavelength. We compare our accurate DDA results with two other approximations: the analytical multi-layer sphere (MLS) model and effective medium theory (EMT). For high porosity and/or absorptive materials, the MLS model does not provide a good approximation for scattering and absorption by ballistic aggregates. The EMT method provides a much better approximation than the MLS model for these aggregates, with a typical difference less than 20% in extinction and scattering cross sections compared with DDA results, for all types, compositions and wavelengths probed in this study.

According to Babichev et al., the accretion of a phantom test fluid onto a Schwarzschild black hole will induce the mass of the black hole to decrease, however the backreaction was ignored in their calculation. Using new exact solutions describing black holes in a background Friedmann-Robertson-Walker universe, we find that the physical black hole mass may instead increase due to the accretion of phantom energy. If this is the case, and the future universe is dominated by phantom dark energy, the black hole apparent horizon and the cosmic apparent horizon will eventually coincide and, after that, the black hole singularity will become naked in finite comoving time before the Big Rip occurs, violating the Cosmic Censorship Conjecture.

We investigate the estimation of orbital parameters by least-$\chi^2$ Keplerian fits to radial velocity (RV) data using synthetic data sets. We find that while the fitted period is fairly accurate, the best-fit eccentricity and $M_p\sin i$ are systematically biased upward from the true values for low signal-to-noise ratio $K/\sigma\lesssim 3$ and moderate number of observations $N_{\rm obs}\lesssim 60$, leading to a suppression of the number of nearly circular orbits. Assuming intrinsic distributions of orbital parameters, we generate a large number of mock RV data sets and study the selection effect on the eccentricity distribution. We find the overall detection efficiency only mildly decreases with eccentricity. This is because although high eccentricity orbits are more difficult to sample, they also have larger RV amplitudes for fixed planet mass and orbital semi-major axis. Thus the primary source of uncertainties in the eccentricity distribution comes from biases in Keplerian fits to detections with low-amplitude and/or small $N_{\rm obs}$, rather than from selection effects. Our results suggest that the abundance of low-eccentricity exoplanets may be underestimated in the current sample and we urge caution in interpreting the eccentricity distributions of low-amplitude detections in future RV samples.

We compile black hole (BH) masses for $\sim 60,000$ quasars in the redshift range $0.1 \lesssim z \lesssim 4.5$ included in the Fifth Data Release of the Sloan Digital Sky Survey (SDSS), using virial BH mass estimators based on the \hbeta, \MgII, and \CIV emission lines. We find that: (1) within our sample, the widths of the three lines follow log-normal distributions, with means and dispersions that do not depend strongly on luminosity or redshift;(2) the \MgII- and \hbeta-estimated BH masses are consistent with one another; and (3) the \CIV BH mass estimator may be more severely affected by a disk wind component than the \MgII and \hbeta estimators, giving a positive bias in mass correlated with the \CIV-\MgII blueshift. Most SDSS quasars have virial BH masses in the range $10^8-10^9 M_\odot$. There is a clear upper mass limit of $\sim 10^{10} M_\odot$ for active BHs at $z \gtrsim 2$, decreasing at lower redshifts. Making the reasonable assumptions that the underlying BH mass distribution decreases with mass and that the Eddington ratio distribution at fixed BH mass has non-zero width, we show that the measured virial BH mass distribution and Eddington ratio distribution are subject to Malmquist bias. A radio quasar subsample (with $1.5\lesssim z\lesssim 2.3$) has mean virial BH mass larger by $\sim 0.12$ dex than the whole sample. A broad absorption line (BAL) quasar subsample (with $1.7\lesssim z\lesssim 2.2$) has identical virial mass distribution as the nonBAL sample, with no mean offset. (Abridged)

We select a sample of $\sim 4200$ traditionally defined broad absorption line quasars (BALQs) from the Fifth Data Release quasar catalog of the Sloan Digital Sky Survey. For a statistically homogeneous quasar sample with $1.7\le z\le 4.2$, the BAL quasar fraction is $\sim 14%$ and is almost constant with redshift. We measure the auto-correlation of non-BAL quasars (nonBALQs) and the cross-correlation of BALQs with nonBALQs using this statistically homogeneous sample, both in redshift space and using the projected correlation function. We find no significant difference between the clustering strengths of BALQs and nonBALQs. Assuming a power-law model for the real space correlation function $\xi(r)=(r/r_0)^{-1.8}$, the correlation length for nonBALQs is $r_0=7.6\pm 0.8 h^{-1}{\rm Mpc}$; for BALQs, the cross-correlation length is $r_0=7.4\pm 1.1 h^{-1}{\rm Mpc}$. Our clustering results suggest that BALQs live in similar large-scale environments as do nonBALQs.

We present the fourth edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog. The catalog contains 77,429 objects; this is an increase of over 30,000 entries since the previous edition. The catalog consists of the objects in the SDSS Fifth Data Release that have luminosities larger than M_i = -22.0 (in a cosmology with H_0 = 70 km/s/Mpc, Omega_M = 0.3, and Omega_Lambda = 0.7) have at least one emission line with FWHM larger than 1000 km/s, or have interesting/complex absorption features, are fainter than i=15.0, and have highly reliable redshifts. The area covered by the catalog is 5740 sq. deg. The quasar redshifts range from 0.08 to 5.41, with a median value of 1.48; the catalog includes 891 quasars at redshifts greater than four, of which 36 are at redshifts greater than five. Approximately half of the catalog quasars have i < 19; nearly all have i < 21. For each object the catalog presents positions accurate to better than 0.2 arcsec. rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains basic radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800--9200A at a spectral resolution of ~2000. The spectra can be retrieved from the public database using the information provided in the catalog. The average SDSS colors of quasars as a function of redshift, derived from the catalog entries, are presented in tabular form. Approximately 96% of the objects in the catalog were discovered by the SDSS.

(Abridged) We study the two-point correlation function of a uniformly selected sample of 4,426 luminous optical quasars with redshift $2.9 \le z\le 5.4$ selected over 4041 deg$^2$ from the Fifth Data Release of the Sloan Digital Sky Survey. For a real-space correlation function of the form $\xi(r)=(r/r_0)^{-\gamma}$, the fitted parameters in comoving coordinates are $r_0 = 15.2 \pm 2.7 h^{-1}$ Mpc and $\gamma = 2.0 \pm 0.3$, over a scale range $4\le r_p\le 150 h^{-1}$ Mpc. Thus high-redshift quasars are appreciably more strongly clustered than their $z \approx 1.5$ counterparts, which have a comoving clustering length $r_0 \approx 6.5 h^{-1}$ Mpc. Dividing our sample into two redshift bins: $2.9\le z\le 3.5$ and $z\ge 3.5$, and assuming a power-law index $\gamma=2.0$, we find a correlation length of $r_0 = 16.9 \pm 1.7 h^{-1}$ Mpc for the former, and $r_0 = 24.3 \pm 2.4 h^{-1}$ Mpc for the latter. Following Martini & Weinberg, we relate the clustering strength and quasar number density to the quasar lifetimes and duty cycle. Using the Sheth & Tormen halo mass function, the quasar lifetime is estimated to lie in the range $4\sim 50$ Myr for quasars with $2.9\le z\le 3.5$; and $30\sim 600$ Myr for quasars with $z\ge 3.5$. The corresponding duty cycles are $0.004\sim 0.05$ for the lower redshift bin and $0.03\sim 0.6$ for the higher redshift bin. The minimum mass of halos in which these quasars reside is $2-3\times 10^{12}\ h^{-1}M_\odot$ for quasars with $2.9\le z\le 3.5$ and $4-6\times 10^{12}\ h^{-1}M_\odot$ for quasars with $z\ge 3.5$.

We combine optical and X-ray data for eight low redshift ($z\sim 0.06$) poor groups of galaxies from the {\it XI}({\it XMM/IMACS}) Groups Project to study the AGN population in the group environment. Among $\sim 140$ group members, we identify five AGN based on their optical emission lines. None of these optically-selected AGN are detected by {\it XMM-Newton}. One additional AGN is discovered in the {\it XMM-Newton} observations. This X-ray detected AGN, which has no obvious AGN emission line signatures in its optical spectrum, is a member of the only X-ray luminous group in our sample. The lack of a significant population of X-ray bright, but optically dull AGN among less dynamically evolved groups is in stark contrast to the large fraction of such objects in rich clusters of galaxies (Martini et al. 2006). We suggest this result can be explained by a physical scenario for AGN accretion evolution: AGN activity is initially triggered by galaxy merging, leading to a high accretion rate and an optically dominant phase (via thin disk accretion). As the accretion rate drops in time, the AGN gradually enters an X-ray dominant low-accretion phase (via a radiative inefficient accretion flow). In this picture, optical- and X-ray-selected AGN are the same population of supermassive black holes observed at different epochs. Within the context of this scenario, the majority of AGN in poor groups are in the high-accretion optically dominant phase, while the AGN population in rich clusters is mostly in the low-accretion X-ray dominant phase.

We carry out three-dimensional, high resolution (up to $1024^2\times 256$) hydrodynamic simulations of the evolution of vortices in vertically unstratified Keplerian disks using the shearing sheet approximation. The transient amplification of incompressible, linear amplitude leading waves (which has been proposed as a possible route to nonlinear hydrodynamical turbulence in disks) is used as one test of our algorithms; our methods accurately capture the predicted amplification, converges at second-order, and is free from aliasing. Waves expected to reach nonlinear amplitude at peak amplification become unstable to Kelvin-Helmholtz modes when $\mid W_{\rm max}\mid\gtrsim \Omega$ (where $W_{\rm max}$ is the local maximum of vorticity and $\Omega$ the angular velocity). We study the evolution of a power-law distribution of vorticity consistent with Kolmogorov turbulence; in two-dimensions long-lived vortices emerge and decay slowly, similar to previous studies. In three-dimensions, however, vortices are unstable to bending modes, leading to rapid decay. Only vortices with a length to width ratio smaller than one survive; in three-dimensions the residual kinetic energy and shear stress is at least one order of magnitude smaller than in two-dimensions. No evidence for sustained hydrodynamical turbulence and transport is observed in three-dimensions. Instead, at late times the residual transport is determined by the amplitude of slowly decaying, large-scale vortices (with horizontal extent comparable to the scale height of the disk), with additional contributions from nearly incompressible inertial waves possible. Evaluating the role that large-scale vortices play in astrophysical accretion disks will require understanding the mechanisms that generate and destroy them.

We study the dispersal of a gaseous disc surrounding a central high-mass stellar core once this circumstellar disc becomes fully ionized. If the stellar and surrounding EUV and X-ray radiations are so strong as to rapidly heat up and ionize the entire circumstellar disc as further facilitated by disc magnetohydrodynamic (MHD) turbulence, a shock can be driven to travel outward in the fully ionized disc, behind which the disc expands and thins. For an extremely massive and powerful stellar core, the ionized gas pressure overwhelms the centrifugal and gravitational forces in the disc. In this limit, we construct self-similar shock solutions for such an expansion and depletion phase. As a significant amount of circumstellar gas being removed, the relic disc becomes vulnerable to strong stellar winds and fragments into clumps. We speculate that disc disappearance happens rapidly, perhaps on a timescale of $\sim 10^3-10^4\hbox{yr}$ once the disc becomes entirely ionized sometime after the onset of thermal nuclear burning in a high-mass stellar core.

Given a large-scale mixture of self-interacting dark matter (SIDM) particles and baryon matter distributed in the early Universe, we advance here a two-phase accretion scenario for forming supermassive black holes (SMBHs) with masses around $\sim 10^9 M_{\odot}$ at high redshifts $z (\gsim 6)$. The first phase is conceived to involve a rapid quasi-spherical and quasi-steady Bondi accretion of mainly SIDM particles embedded with baryon matter onto seed black holes (BHs) created at redshifts $z\lsim 30$ by the first generation of massive Population III stars; this earlier phase rapidly gives birth to significantly enlarged seed BH masses of $M_{\hbox{\tiny BH},t_1}\backsimeq 1.4\times 10^6\ M_\odot \sigma_0/(1\hbox{cm}^2\hbox{g}^{-1})(C_s/30\hbox{km s}^{-1})^4$ during $z\sim 20-15$, where $\sigma_0$ is the cross section per unit mass of SIDM particles and $C_s$ is the velocity dispersion in the SIDM halo referred to as an effective "sound speed". The second phase of BH mass growth is envisaged to proceed primarily via baryon accretion, eventually leading to SMBH masses of $M_{\hbox{\tiny BH}}\sim 10^9 M_\odot$; such SMBHs may form either by $z\sim 6$ for a sustained accretion at the Eddington limit or later at lower $z$ for sub-Eddington mean accretion rates. We intend to account for the reported detections of a few SMBHs at early epochs, e.g., SDSS 1148+5251 and so forth, without necessarily resorting to either super-Eddington baryon accretion or very frequent BH merging processes. Only extremely massive dark SIDM halos associated with rare peaks of density fluctuations in the early Universe may harbour such early SMBHs or quasars. Observational consequences are discussed.

The role of the auxiliary scalar field $\phi$ of Brans-Dicke theory played in baryon number asymmetry is discussed in this paper. We consider a derivative coupling of this gravitational scalar field to the baryon current ${J^{\mu}}_B$ or the current of the baryon number minus lepton number ${J^{\mu}}_{B-L}$ based on a series of works of R. Morganstern about the Brans-Dicke cosmology. We find that the spontaneous baryogenesis by this coupling is capable to yield a sufficient baryon asymmetry $n_B/s\sim 10^{10} $ for the time of the grand unification is in a little advanced. In addition, Davoudiasl et al have recently introduced a new type of interaction between the Ricci scalar $R$ and the bayon current $J^{\mu}$, $\partial_{\mu}R J^{\mu}$ and also proposed a mechanism for baryogenesis, the gravitational baryogenesis. However, the Einstein equation tell us that $\dot{R}=0$ in the radiation-dominated epoch of the standard FRW cosmology. In this paper we reconsider the feasibility of having gravitational baryongenesis with such a form of interaction in radiation-filled Brans-Dicke cosmology. We will show that $\dot{R}$ does not vanish in this case and the required baryon number asymmetry can also be achieved.

We construct analytically stationary global configurations for both aligned and logarithmic spiral coplanar magnetohydrodynamic (MHD) perturbations in an axisymmetric background MHD disc with a power-law surface mass density $\Sigma_0\propto r^{-\alpha}$, a coplanar azimuthal magnetic field $B_0\propto r^{-\gamma}$, a consistent self-gravity and a power-law rotation curve $v_0\propto r^{-\beta}$ where $v_0$ is the linear azimuthal gas rotation speed. The barotropic equation of state $\Pi\propto\Sigma^{n}$ is adopted for both MHD background equilibrium and coplanar MHD perturbations where $\Pi$ is the vertically integrated pressure and $n$ is the barotropic index. For a scale-free background MHD equilibrium, a relation exists among $\alpha$, $\beta$, $\gamma$ and $n$ such that only one parameter (e.g., $\beta$) is independent. For a linear axisymmetric stability analysis, we provide global criteria in various parameter regimes. For nonaxisymmetric aligned and logarithmic spiral cases, two branches of perturbation modes (i.e., fast and slow MHD density waves) can be derived once $\beta$ is specified. To complement the magnetized singular isothermal disc (MSID) analysis of Lou, we extend the analysis to a wider range of $-1/4<\beta<1/2$. As an example of illustration, we discuss specifically the $\beta=1/4$ case when the background magnetic field is force-free. Angular momentum conservation for coplanar MHD perturbations and other relevant aspects of our approach are discussed.

We seek for self-similar solutions describing the time-dependent evolution of self-gravity systems with either spherical symmetry or axisymmetric disk geometry. By assuming self-similar variable $x\equiv r/at$ where $a$ is isothermal sound speed we find self-similar solutions extending from the initial instant $t=0$ to the final stage $t\to \infty$ using standard semi-analytical methods. Different types of solutions are constructed, which describe overall expansion or collapse, envelope expansion with core collapse (EECC), the formation of central rotationally supported quasi-equilibrium disk as well as shocks. Though infinitely many, these self-similarity solutions have similar asymptotic behaviors which may impose diagnosis on the velocity and density structures in astrophysical systems.

We propose self-similar shocked flow models for certain dynamical evolution phases of young stellar objects (YSOs), `champagne flows' of H {\sevenrm II} regions surrounding OB stars and shaping processes of planetary nebulae (PNe). We analyze an isothermal fluid of spherical symmetry and construct families of similarity shocked flow solutions featured by: 1. either a core expansion with a finite central density or a core accretion at constant rate with a density scaling $\propto r^{-3/2}$; 2. a shock moving outward at a constant speed; 3. a preshock gas approaching a constant speed at large $r$ with a density scaling $\propto r^{-2}$. In addition to testing numerical codes, our models can accommodate diverse shocked flows with or without a core collapse or outflow and an envelope expansion or contraction. As an application, we introduce our model analysis to observations of Bok globule B335.

In a composite system of gravitationally coupled stellar and gaseous discs, we perform linear stability analysis for axisymmetric coplanar perturbations using the two-fluid formalism. The background stellar and gaseous discs are taken to be scale-free with all physical variables varying as powers of cylindrical radius $r$ with compatible exponents. The unstable modes set in as neutral modes or stationary perturbation configurations with angular frequency $\omega=0$.

In a composite fluid system of two gravitationally coupled barotropic scale-free discs bearing a rotation curve $v\propto r^{-\beta}$ and a power-law surface mass density $\Sigma_0\propto r^{-\alpha}$ with $\alpha=1+2\beta$, we construct coplanar stationary aligned and spiral perturbation configurations in the two discs. Due to the mutual gravitational interaction, there are two independent classes of perturbation modes with surface mass density disturbances in the two coupled discs being either in-phase or out-of-phase. We derive analytical criteria for such perturbation modes to exist and show numerical examples. We compute the aligned and spiral perturbation modes systematically to explore the entire parameter regime. For the axisymmetric $m=0$ case with radial oscillations, there are two unstable regimes of ring-fragmentation and collapse corresponding to short and long radial wavelengths, respectively. Only within a certain range of the rotation parameter $D_s^2$ (square of the effective Mach number for the stellar disc), can a composite disc system be stable against all axisymmetric perturbations. Compared with a single-disc system, the coupled two-disc system becomes less stable against such axisymmetric instabilities.

The effects of the generalized uncertainty principle({\bf GUP}) on the cosmological constant problem are discussed in the Schwarzchild-deSitter spacetime, through studying the corrections to its thermodynamics. We derive the correction to the Hawking temperature of the cosmological horizon, by a heuristic method enlighten by gr-qc/0106080 . The logarithmic correction to the Bekenstein-Hawking entropy is also obtained. For an ordinary star (not a black hole), the probable ratio of the vacuum energy to the total energy within the cosmological horizon is 2/3, which roughly coincides with the evidences from the astronomical observations. For a black hole, the ratio tends to decrease. An inequality associating the energy density with the length of system is put forward for understanding the smallness of the cosmological constant, and the relation between the Bekenstein entropy bound and the Bekenstein-Hawking entropy is also briefly discussed.

We investigate self-similar dynamical processes in an isothermal self-gravitational fluid with spherical symmetry. In reference to earlier complementary solution results of Larson, Penston, Shu, Hunter and Whitworth & Summers, we further explore the `semi-complete solution space' from an initial instant $t\to 0^{+}$ to a final stage $t\to +\infty$. These similarity solutions can describe and accommodate physical processes of radial inflow, core collapse, oscillations and envelope expansion (namely, outflow or wind) or contraction as well as shocks. In particular, we present new classes of self-similar solutions, referred to as `envelope expansion with core collapse' (EECC) solutions, that are featured by an interior core collapse and an exterior envelope expansion concurrently. The interior collapse towards the central core approaches a free-fall state as the radius $r\to 0$, while the exterior envelope expansion gradually approaches a constant radial flow speed as $r\to +\infty$. There exists at least one spherical stagnation surface of zero flow speed that separates the core collapse and the envelope outflow and that travels outward at constant speed, either subsonically or supersonically, in a self-similar manner.

We investigate properties of stationary aligned and unaligned spiral perturbation configurations in a composite system of gravitationally coupled stellar and gaseous singular isothermal disks (SIDs) using the two-fluid formalism. We obtain stationary perturbation solutions for aligned and unaligned spiral logarithmic configurations in such a composite SID system and derive analytically existence criteria for these solutions. In comparison with the similar problem of a single SID studied by Shu et al. (2000), there are now two possible sets of solutions owing to an additional SID.

Using the two-fluid formalism with the polytropic approximation, we examine the axisymmetric stability criterion for a composite system of gravitationally coupled stellar and gaseous singular isothermal discs (SIDs).