Search

We carried out spectroscopic monitoring of 21 low-redshift Seyfert 1 galaxies using the Kast double spectrograph on the 3-m Shane telescope at Lick Observatory from April 2016 to May 2017. Targeting active galactic nuclei (AGN) with luminosities of {\lambda}L{\lambda} (5100 {\AA}) = 10^44 erg/s and predicted H{\beta} lags of 20-30 days or black hole masses of 10^7-10^8.5 Msun, our campaign probes luminosity-dependent trends in broad-line region (BLR) structure and dynamics as well as to improve calibrations for single-epoch estimates of quasar black hole masses. Here we present the first results from the campaign, including H{\beta} emission-line light curves, integrated H{\beta} lag times (8-30 days) measured against V-band continuum light curves, velocity-resolved reverberation lags, line widths of the broad H{\beta} components, and virial black hole mass estimates (10^7.1-10^8.1 Msun). Our results add significantly to the number of existing velocity-resolved lag measurements and reveal a diversity of BLR gas kinematics at moderately high AGN luminosities. AGN continuum luminosity appears not to be correlated with the type of kinematics that its BLR gas may exhibit. Follow-up direct modeling of this dataset will elucidate the detailed kinematics and provide robust dynamical black hole masses for several objects in this sample., Comment: 34 pages, 14 figures, accepted for publication in ApJ

A supermassive black hole SgrA* with the mass ~4x10^6 Msun resides at the centre of our galaxy. Building up such a massive black hole within the ~10^10 year lifetime of our galaxy would require a mean accretion rate of ~4x10^-4 Msun/yr. At present, X-ray observations constrain the rate of hot gas accretion at the Bondi radius (10^5 R_Sch = 0.04 pc at 8kpc) to \dot{M}_Bondi ~ 3x10^-6 Msun/yr, and polarization measurements constrain it near the event horizon to \dot{M}_horizon ~ 10^{-8} Msun/yr. A range of models was developed to describe the accretion gas onto an underfed black hole. However, the exact physics still remains to be understood. One challenge with the radiation inefficient accretion flows is that even if one understands the dynamics there is no accepted prescription for associating emissivity (and absorption) with the flow. The other issue is the lack of model-independent probes of accretion flow at intermediate radii (between few and ~ 10^5 R_Sch), i.e. the constraints that do not assume a model of accretion flow as an input parameter. We report detection and imaging of the 10^4 K ionized gas disk within 2x10^4 R_Sch in a mm hydrogen recombination line H30alpha: n = 31 -> 30 at 231.9 GHz using the ALMA. The emission was detected with a double-peaked line profile spanning full width of 2,200 km/s with the approaching and the receding components straddling Sgr A*, each offset from it by 0.11arcsec= 0.004pc. The red-shifted side is displaced to the north-east, while the blue-shifted side is displaced to the south-west. The limit on the total mass of ionized gas estimated from the emission is 10^-4 - 10^-5 Sun at a mean hydrogen density 10^5-10^6 cm^-3, depending upon whether or not we assume the presence of a uniform density disk or an ensemble of orbiting clouds, and the amplification factor of the mm radiation due to the strong background source which is Sgr A* continuum., Comment: Authors' own extended version

(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., Comment: 9 chapters, 301 pages, 100 figures. This version of the DSC is a comprehensive update of the original version, released in 2016, which can be downloaded at arXiv:1606.00043. A detailed summary of the design of MSE is available in the MSE Book 2018, available at arXiv:1810.08695

We present models of the H$\beta$-emitting broad-line region (BLR) in seven Seyfert 1 galaxies from the Lick AGN (Active Galactic Nucleus) Monitoring Project 2011 sample, drawing inferences on the BLR structure and dynamics as well as the mass of the central supermassive black hole. We find that the BLR is generally a thick disk, viewed close to face-on, with preferential emission back toward the ionizing source. The dynamics in our sample range from near-circular elliptical orbits to inflowing or outflowing trajectories. We measure black hole masses of $\log_{10}(M_{\rm BH}/M_\odot) = 6.48^{+0.21}_{-0.18}$ for PG 1310$-$108, $7.50^{+0.25}_{-0.18}$ for Mrk 50, $7.46^{+0.15}_{-0.21}$ for Mrk 141, $7.58^{+0.08}_{-0.08}$ for Mrk 279, $7.11^{+0.20}_{-0.17}$ for Mrk 1511, $6.65^{+0.27}_{-0.15}$ for NGC 4593, and $6.94^{+0.14}_{-0.14}$ for Zw 229$-$015. We use these black hole mass measurements along with cross-correlation time lags and line widths to recover the scale factor $f$ used in traditional reverberation mapping measurements. Combining our results with other studies that use this modeling technique, bringing our sample size to 16, we calculate a scale factor that can be used for measuring black hole masses in other reverberation mapping campaigns. When using the root-mean-square (rms) spectrum and using the line dispersion to measure the line width, we find $\log_{10}(f_{{\rm rms},\sigma})_{\rm pred} = 0.57 \pm 0.19$. Finally, we search for correlations between $f$ and other AGN and BLR parameters and find marginal evidence that $f$ is correlated with $M_{\rm BH}$ and the BLR inclination angle, but no significant evidence of a correlation with the AGN luminosity or Eddington ratio., Comment: 26 pages, 14 figures. Accepted for publication in ApJ

We consider AGN photoionization of gas in companion galaxies (cross-ionization) as a way to sample the intensity of AGN radiation in both direction and time, independent of the gas properties of the AGN host galaxies. From an initial set of 212 AGN+companion systems, identified with the help of Galaxy Zoo participants, we obtained long-slit optical spectra of 32 pairs which were a priori likely to show cross-ionization based on projected separation or angular extent of the companion. From emission-line ratios, 10 of these systems are candidates for cross-ionization, roughly the fraction expected if most AGN have ionization cones with 70-degree opening angles. Among these, Was 49 remains the strongest nearby candidate. NGC 5278/9 and UGC 6081 are dual-AGN systems with tidal debris, complicating identification of cross-ionization. The two weak AGN in the NGC 5278/9 system ionize gas filaments to a projected radius 14 kpc from each galaxy. In UGC 6081, an irregular high-ionization emission region encompasses both AGN, extending more than 15 kpc from each. The observed AGN companion galaxies with and without signs of external AGN photoionization have similar distributions in estimated incident AGN flux, suggesting that geometry of escaping radiation or long-term variability control this facet of the AGN environment. This parallels conclusions for luminous QSOs based on the proximity effect among Lyman-alpha absorbers. In some galaxies, mismatch between spectroscopic classifications in the common BPT diagram and the intensity of weaker He II and [Ne V] emission lines highlights the limits of common classifications in low-metallicity environments., Comment: Accepted version to appear in MNRAS. New version has spectra showing AGN ionization in NGC 5278/9 filaments and tunable-filter mapping of clouds around UGC 6081 system

We provide an updated calibration of CIV $\lambda1549$ broad emission line-based single-epoch (SE) black hole (BH) mass estimators for active galactic nuclei (AGNs) using new data for six reverberation-mapped AGNs at redshift $z=0.005-0.028$ with BH masses (bolometric luminosities) in the range $10^{6.5}-10^{7.5}$ $M_{\odot}$ ($10^{41.7}-10^{43.8}$ erg s$^{\rm -1}$). New rest-frame UV-to-optical spectra covering 1150-5700 \AA\ for the six AGNs were obtained with the Hubble Space Telescope (HST). Multi-component spectral decompositions of the HST spectra were used to measure SE emission-line widths for the CIV, MgII, and H$\beta$ lines as well as continuum luminosities in the spectral region around each line. We combine the new data with similar measurements for a previous archival sample of 25 AGNs to derive the most consistent and accurate calibrations of the CIV-based SE BH mass estimators against the H$\beta$ reverberation-based masses, using three different measures of broad-line width: full-width at half maximum (FWHM), line dispersion ($\sigma_{\rm line}$) and mean absolute deviation (MAD). The newly expanded sample at redshift $z=0.005-0.234$ covers a dynamic range in BH mass (bolometric luminosity) of $\log\ M_{\rm BH}/M_{\odot} = 6.5-9.1$ ($\log\ L_{\rm bol}/$erg s$^{\rm -1}=41.7-46.9$), and we derive the new CIV-based mass estimators using a Bayesian linear regression analysis over this range. We generally recommend the use of $\sigma_{\rm line}$ or MAD rather than FWHM to obtain a less biased velocity measurement of the CIV emission line, because its narrow-line component contribution is difficult to decompose from the broad-line profile., Comment: Accepted for publication in ApJ

We consider energy budgets and radiative history of 8 fading AGN, identified from mismatch between the ionizion of very extended (>10 kpc) ionized clouds and the luminosity of the nucleus viewed directly. All show significant fading on ~50,000-year timescales. We explore the use of minimum ionizing luminosity Q derived from photoionization balance in the brightest pixels in H-alpha at each projected radius. Tests using PG QSOs, and one target with detailed photoionization modeling, suggest that we can derive useful histories of individual AGN; the minimum ionizing luminosity is always an underestimate and subject to fine structure in the ionized material. These tests suggest that the underestimation from the upper envelope of Q values is roughly constant for a given object. These AGN show rapid drops and standstills; the common feature is a rapid drop in the last 20,000 years before our view of the nucleus. E-folding timescales are mostly thousands of years, with a few episodes as short as 400. In the limit of largely obscured AGN, we find additional evidence for fading, comparing lower limits from recombination balance and the maximum luminosities derived from from infrared fluxes. We compare these long-term light curves to simulations of AGN accretion; the strongest variations on these timespans are seen in models with strong and local feedback. Gemini integral-field optical spectroscopy shows a very limited role for outflows in these structures. While rings and loops of emission are common, their kinematic structure shows some to be in regular rotation. UGC 7342 exhibits local signatures of outflows <300 km/s, largely associated with very diffuse emission. Only in the Teacup AGN do we see outflow signatures of order 1000 km/s. Clouds around these fading AGN consist largely of tidal debris being externally illuminated but not displaced by AGN outflows. (Abridged), Comment: Accepted by Astrophysical Journal; text is accepted version. Some figures have reduced resolution

We report the discovery of 14 Lyman-alpha blobs (LABs) at z~0.3, existing at least 4-7 billion years later in the Universe than all other LABs known. Their optical diameters are 20-70 kpc, and GALEX data imply Ly-alpha luminosities of (0.4-6.3)x10^43 erg/s. Contrary to high-z LABs, they live in low-density areas. They are ionized by AGN, suggesting that cold accretion streams as a power source must deplete between z=2 and z=0.3. We also show that transient AGN naturally explain the ionization deficits observed in many LABs: Their Ly-alpha and X-ray fluxes decorrelate below 10^6 years because of the delayed escape of resonantly scattering Ly-alpha photons. High Ly-alpha luminosities do not require currently powerful AGN, independent of obscuration. Chandra X-ray data reveal intrinsically weak AGN, confirming the luminous optical nebulae as impressive ionization echoes. For the first time, we also report mid-infrared thermal echoes from the dusty tori. We conclude that the AGN have faded by 3-4 orders of magnitude within the last 10^(4-5) years, leaving fossil UV, optical and thermal radiation behind. The host galaxies belong to the group of previously discovered Green Bean galaxies (GBs). Gemini optical imaging reveals smooth spheres, mergers, spectacular outflows and ionization cones. Because of their proximity and high flux densities, GBs are perfect targets to study AGN feedback, mode switching and the Ly-alpha escape. The fully calibrated, coadded optical FITS images are publicly available., Comment: 36 pages, 27 figures. Accepted for publication in MNRAS 2016 July 21. The coadded gri FITS images are publicly available at https://zenodo.org/record/56059

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), Comment: 210 pages, 91 figures. Exposure draft. Appendices to the Detailed Science Case can be found at http://mse.cfht.hawaii.edu/docs/

Reverberation mapping probes the structure of the broad emission-line region (BLR) in active galactic nuclei (AGN). The kinematics of the BLR gas can be used to measure the mass of the central supermassive black hole. The main uncertainty affecting black hole mass determinations is the structure of the BLR. We present a new method for reverberation mapping based on regularized linear inversion (RLI) that includes modelling of the AGN continuum light curves. This enables fast calculation of velocity-resolved response maps to constrain BLR structure. RLI allows for negative response, such as when some areas of the BLR respond in inverse proportion to a change in ionizing continuum luminosity. We present time delays, integrated response functions, and velocity-delay maps for the $\rm{H}\,\beta$ broad emission line in five nearby AGN, as well as for $\rm{H}\,\alpha$ and $\rm{H}\,\gamma$ in Arp 151, using data from the Lick AGN Monitoring Project 2008. We find indications of prompt response in three of the objects (Arp 151, NGC 5548 and SBS 1116+583A) with additional prompt response in the red wing of $\rm{H}\,\beta$. In SBS 1116+583A we find evidence for a multimodal broad prompt response followed by a second narrow response at 10 days. We find no clear indications of negative response. The results are complementary to, and consistent with, other methods such as cross correlation, maximum entropy and dynamical modelling. Regularized linear inversion with continuum light curve modelling provides a fast, complementary method for velocity-resolved reverberation mapping and is suitable for use on large datasets., Comment: 20 pages, 13 figures, accepted to MNRAS

We present narrow- and medium-band HST imaging, with additional supporting ground-based data, for 8 galaxies identified as hosting fading AGN. These have AGN-ionized gas projected >10 kpc from the nucleus, and significant shortfall of ionizing radiation between the distant gas and the AGN, indicating fading AGN on ~50,000-year timescales. Every system shows evidence of ongoing or past interactions; a similar sample of obscured AGN with extended ionized clouds shares this incidence of disturbances. Several systems show multiple dust lanes in different orientations, broadly fit by differentially precessing disks of accreted material ~1.5 Gyr after initial arrival. The gas has lower metallicity than the nuclei; three systems have abundances uniformly well below solar, consistent with an origin in tidally disrupted low-luminosity galaxies, while some systems have more nearly solar abundances (accompanied by such signatures as multiple Doppler components), which may suggest redistribution of gas by outflows within the host galaxies themselves. These aspects are consistent with a tidal origin for the extended gas in most systems, although the ionized gas and stellar tidal features do not always match closely. In contrast to clouds near radio-loud AGN, these are dominated by rotation, in some cases in warped disks. Outflows are important only in localized regions near some of the AGN. In UGC 7342 and UGC 11185, luminous star clusters are seen within projected ionization cones, potentially triggered by outflows. As in the discovery example Hanny's Voorwerp/IC 2497, some clouds lack a strong correlation between H-alpha surface brightness and ionization parameter, indicating unresolved fine structure. Together with thin coherent filaments spanning several kpc, persistence of these structures over their orbital lifetimes may require a role for magnetic confinement. (Abridged), Comment: Accepted by the Astronomical Journal. This revised version adds spectroscopic estimates of abundances and new kinematic data. Figures compressed for space

We present an improved and expanded simply parameterized phenomenological model of the broad line region (BLR) in active galactic nuclei (AGN) for modeling reverberation mapping data. By modeling reverberation mapping data directly, we can constrain the geometry and dynamics of the BLR and measure the black hole mass without relying on the normalization factor needed in the traditional analysis. For realistic simulated reverberation mapping datasets of high-quality, we can recover the black hole mass to $0.05-0.25$ dex uncertainty and distinguish between dynamics dominated by elliptical orbits and inflowing gas. While direct modeling of the integrated emission line light curve allows for measurement of the mean time lag, other details of the geometry of the BLR are better constrained by the full spectroscopic dataset of emission line profiles. We use this improved model of the BLR to explore possible sources of uncertainty in measurements of the time lag using cross-correlation function (CCF) analysis and in measurements of the black hole mass using the virial product. Sampling the range of geometries and dynamics in our model of the BLR suggests that the theoretical uncertainty in black hole masses measured using the virial product is on the order of 0.25 dex. These results support the use of the CCF to measure time lags and the virial product to measure black hole masses when direct modeling techniques cannot be applied, provided the uncertainties associated with the interpretation of the results are taken into account., Comment: 20 pages, 14 figures, 4 tables, published in MNRAS

We present dynamical modeling of the broad line region (BLR) in the Seyfert 1 galaxy Mrk 50 using reverberation mapping data taken as part of the Lick AGN Monitoring Project (LAMP) 2011. We model the reverberation mapping data directly, constraining the geometry and kinematics of the BLR, as well as deriving a black hole mass estimate that does not depend on a normalizing factor or virial coefficient. We find that the geometry of the BLR in Mrk 50 is a nearly face-on thick disk, with a mean radius of 9.6(+1.2,-0.9) light days, a width of the BLR of 6.9(+1.2,-1.1) light days, and a disk opening angle of 25\pm10 degrees above the plane. We also constrain the inclination angle to be 9(+7,-5) degrees, close to face-on. Finally, the black hole mass of Mrk 50 is inferred to be log10(M(BH)/Msun) = 7.57(+0.44,-0.27). By comparison to the virial black hole mass estimate from traditional reverberation mapping analysis, we find the normalizing constant (virial coefficient) to be log10(f) = 0.78(+0.44,-0.27), consistent with the commonly adopted mean value of 0.74 based on aligning the M(BH)-{\sigma}* relation for AGN and quiescent galaxies. While our dynamical model includes the possibility of a net inflow or outflow in the BLR, we cannot distinguish between these two scenarios., Comment: Accepted for publication in ApJ. 8 pages, 6 figures

We present a general method to analyze reverberation mapping data that provides both estimates for the black hole mass and for the geometry and dynamics of the broad line region (BLR) in active galactic nuclei (AGN). Our method directly infers the spatial and velocity distribution of the BLR from the data, allowing us to easily derive a velocity-resolved transfer function and allowing for a self-consistent estimate of the black hole mass without a virial coefficient. We obtain estimates and reasonable uncertainties of the BLR model parameters by implementing a Markov Chain Monte Carlo algorithm using the formalism of Bayesian probability theory. We use Gaussian Processes to interpolate the the continuum light curve data and create mock light curves that can be fitted to the data. We test our method by creating simulated reverberation mapping data-sets with known true parameter values and by trying to recover these parameter values using our models. We are able to recover the parameters with realistic uncertainties that depend upon the variability of the AGN and the quality of the reverberation mapping campaign. With a geometry model we can recover the mean radius of the BLR to within ~0.1dex random uncertainty for simulated data with an integrated line flux uncertainty of 1.5%, while with a dynamical model we can recover the black hole mass and the mean radius to within ~0.05dex random uncertainty, for simulated data with a line profile average signal to noise ratio of 4 per spectral pixel. These uncertainties do not include modeling errors, which are likely to be present in the analysis of real data, and should therefore be considered as lower limits to the accuracy of the method., Comment: 16 pages, 16 Figures, Accepted to ApJ

Abridged: We present a detailed investigation of the morphological distribution and level of star formation and dust obscuration in the nearby tidally distorted galaxy NGC2442. Spitzer images in the IR at 3.6, 4.5, 5.8, 8.0um, and 24um and GALEX images at 1500\AA{} and 2300\AA{} allow us to resolve the galaxy on scales between 240-600pc. We supplement these with archival data in the B, J, H, and K bands. We use the 8um, 24um and FUV (1500\AA) emission to study the star formation rate (SFR). We find that globally, these tracers of star formation give a range of results of ~6-11\msun/yr, with the dust-corrected FUV giving the highest value of SFR. We can reconcile the UV and IR-based estimates by adopting a steeper UV extinction curve that lies in between the starburst (Calzetti) and SMC extinction curves. However, the regions of highest SFR intensity along the spiral arms are consistent with a starburst-like extinction. Overall, the level of star-formation we find is higher than previously published for this galaxy, by about a factor of two, which, contrary to previous conclusions, implies that the interaction that caused the distorted morphology of NGC2442 likely also triggered increased levels of star-formation activity. Outside of the spiral arms, we discover what appears to be a superbubble, ~1.7kpc across in the IRAC images. Significant H{\alpha}, UV and IR emission in the area also suggest vigorous ongoing star-formation. A known, recent supernova (SN1999ga) is located at the edge of this superbubble. Although speculative at this stage, this area suggests a large star-forming region with a morphology shaped by generations of supernovae. Lastly, we discover an 8um (PAH) circumnuclear ring with an ~0.8kpc radius. The H{\alpha} emission is largely concentrated inside that ring and shows a vague spiral structure in the rest of the galaxy., Comment: 16 pages, 13 figures, accepted for publication in ApJ