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 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., Comment: 18 pages, 17 Figures, 3 Tables, submitted to 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) for a sample of five Seyfert 1 galaxies using reverberation mapping data taken by the Lick AGN Monitoring Project in 2008. By modeling the AGN continuum light curve and H$\beta$ line profiles directly we are able to constrain the geometry and kinematics of the BLR and make a measurement of the black hole mass that does not depend upon the virial factor, $f$, needed in traditional reverberation mapping analysis. We find that the geometry of the BLR is generally a thick disk viewed close to face-on. While the H$\beta$ emission is found to come preferentially from the far side of the BLR, the mean size of the BLR is consistent with the lags measured with cross-correlation analysis. The BLR kinematics are found to be consistent with either inflowing motions or elliptical orbits, often with some combination of the two. We measure black hole masses of $\log_{10}(M_{\rm\,BH}/M_\odot)=6.62^{+0.10}_{-0.13}$ for Arp 151, $7.42^{+0.26}_{-0.27}$ for Mrk 1310, $7.51^{+0.23}_{-0.14}$ for NGC 5548, $6.42^{+0.24}_{-0.18}$ for NGC 6814, and $6.99^{+0.32}_{-0.25}$ for SBS 1116+583A. The $f$ factors measured individually for each AGN are found to correlate with inclination angle, although not with $M_{\rm\,BH}$, $L_{5100}$, or FWHM/$\sigma$ of the emission line profile., Comment: 21 pages, 24 figures, 3 tables, Accepted for publication in MNRAS, corrected masses for NGC 5548 and NGC 6814 in the abstract

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

Mapping Mg II resonance emission scattered by galactic winds offers a means to determine the spatial extent and density of the warm outflow. Using Keck/LRIS spectroscopy, we have resolved scattered Mg II emission to the east of 32016857, a star-forming galaxy at z =0.9392 with an outflow. The Mg II emission from this galaxy exhibits a P-Cygni profile, extends further than both the continuum and [O II] emission along the eastern side of the slit, and has a constant Doppler shift along the slit which does not follow the velocity gradient of the nebular [O II] emission. Using the Sobolev approximation, we derive the density of Mg+ ions at a radius of 12 to 18 kpc in the outflow. We model the ionization correction and find that much of the outflowing Mg is in Mg++. We estimate that the total mass flux could be as large as 330 - 500 solar masses per year, with the largest uncertainties coming from the depletion of Mg onto grains and the clumpiness of the warm outflow. We show that confining the warm clouds with a hot wind reduces the estimated mass flux of the warm outflow and indicates amass-loading factor near unity in the warm phase alone. Based on the high blue luminosities that distinguish 32016857 and TKRS 4389, described by Rubin et al. 2011, from other galaxies with P-Cygni emission, we suggest that, as sensitivity to diffuse emission improves, scattering halos may prove to be a generic property of star-forming galaxies at intermediate redshifts., Comment: Manuscript accepted for publication in the Astrophysical Journal. The 22 pages include 2 tables, 10 figures, and an Appendix. Several typos corrected in version 2

The prominent broad Fe II emission blends in the spectra of active galactic nuclei have been shown to vary in response to continuum variations, but past attempts to measure the reverberation lag time of the optical Fe II lines have met with only limited success. Here we report the detection of Fe II reverberation in two Seyfert 1 galaxies, NGC 4593 and Mrk 1511, based on data from a program carried out at Lick Observatory in Spring 2011. Light curves for emission lines including H-beta and Fe II were measured by applying a fitting routine to decompose the spectra into several continuum and emission-line components, and we use cross-correlation techniques to determine the reverberation lags of the emission lines relative to V-band light curves. In both cases the measured lag (t_cen) of Fe II is longer than that of H-beta, although the inferred lags are somewhat sensitive to the choice of Fe II template used in the fit. For spectral decompositions done using the Fe II template of Veron-Cetty et al. (2004), we find t_cen(Fe II)/t_cen(H-beta) = 1.9+-0.6 in NGC 4593 and 1.5+-0.3 in Mrk 1511. The detection of highly correlated variations between Fe II and continuum emission demonstrates that the Fe II emission in these galaxies originates in photoionized gas, located predominantly in the outer portion of the broad-line region., Comment: 14 pages, 6 figures. Accepted for publication in ApJ

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

Some active galactic nuclei (AGN) are surrounded by extended emission-line regions (EELRs), which trace both the illumination pattern of escaping radiation and its history over the light-travel time from the AGN to the gas. From a new set of such EELRs, we present evidence that the AGN in many Seyfert galaxies undergo luminous episodes 20,000-200,000 years in duration. Motivated by the discovery of the spectacular nebula known as Hanny's Voorwerp, ionized by a powerful AGN which has apparently faded dramatically within ~ 100,000 years, Galaxy Zoo volunteers have carried out both targeted and serendipitous searches for similar emission-line clouds around low-redshift galaxies.We present the resulting list of candidates and describe spectroscopy identifying 19 galaxies with AGN-ionized regions at projected radii > 10 kpc. This search recovered known EELRs and identified additional previously unknown cases, one with detected emission to r = 37 kpc. At least 14/19 are in interacting or merging systems; tidal tails are a prime source of extraplanar ionized gas. We see a mix of one- and two-sided structures, with observed cone angles from 23-112 degrees. We consider the energy balance in the ionized clouds, with lower and upper bounds on ionizing luminosity from recombination and ionization-parameter arguments, and estimate the luminosity of the core from the far-infrared data. The implied ratio of ionizing radiation seen by the clouds to that emitted by the nucleus, for a constant nuclear source, ranges from 0.02 to > 12; 7/19 exceed unity. Small values imply heavily obscured AGN. However, large values may require that the AGN has faded over tens of thousands of years, giving us several examples of systems in which such dramatic long-period variation has occurred; this is the only current technique for addressing these timescales in AGN history. (Abridged), Comment: Accepted by MNRAS. Reduced-resolution figures; full-resolution PDF file available until publication at http://astronomy.ua.edu/keel/research/agncloud-rev.pdf Updated to match published version

Supermassive black holes are believed to be ubiquitous at the centers of galaxies. Measuring their masses is extremely challenging yet essential for understanding their role in the formation and evolution of cosmic structure. We present a direct measurement of the mass of a black hole in an active galactic nucleus (Arp 151) based on the motion of the gas responsible for the broad emission lines. By analyzing and modeling spectroscopic and photometric time series, we find that the gas is well described by a disk or torus with an average radius of 3.99 +- 1.25 light days and an opening angle of 68.9 (+21.4, -17.2) degrees, viewed at an inclination angle of 67.8 +- 7.8 degrees (that is, closer to face-on than edge-on). The black hole mass is inferred to be 10^(6.51 +- 0.28) solar masses. The method is fully general and can be used to determine the masses of black holes at arbitrary distances, enabling studies of their evolution over cosmic time., Comment: Accepted for publication in ApJ Letters

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

In the Spring of 2011 we carried out a 2.5 month reverberation mapping campaign using the 3 m Shane telescope at Lick Observatory, monitoring 15 low-redshift Seyfert 1 galaxies. This paper describes the observations, reductions and measurements, and data products from the spectroscopic campaign. The reduced spectra were fitted with a multicomponent model in order to isolate the contributions of various continuum and emission-line components. We present light curves of broad emission lines and the active galactic nucleus (AGN) continuum, and measurements of the broad Hβ line widths in mean and rms spectra. For the most highly variable AGNs we also measured broad H beta line widths and velocity centroids from the nightly spectra. In four AGNs exhibiting the highest variability amplitudes, we detect anticorrelations between broad H beta width and luminosity, demonstrating that the broad-line region "breathes" on short timescales of days to weeks in response to continuum variations. We also find that broad H beta velocity centroids can undergo substantial changes in response to continuum variations; in NGC 4593, the broad H beta velocity shifted by approximately 250 km s(exp -1) over a 1 month period. This reverberation-induced velocity shift effect is likely to contribute a significant source of confusion noise to binary black hole searches that use multi-epoch quasar spectroscopy to detect binary orbital motion. We also present results from simulations that examine biases that can occur in measurement of broad-line widths from rms spectra due to the contributions of continuum variations and photon-counting noise.