Velocity-resolved Reverberation Mapping of Five Bright Seyfert 1 Galaxies

Variability of the broad emission-line fluxes and profiles is commonly seen in the spectra of Type 1 active galactic nuclei (AGNs). A number of isolated cases of dramatic emission-line changes were reported based on photographic spectrograms by the late 1960s and early 1970s (e.g., Andrillat 1968; Pastoriza & Gerola 1970; see the reviews by Pronik 1980 and Collin-Souffrin 1980). Additional and more convincing instances of emission-line changes were found on surprisingly short timescales with the advent of linear detectors for spectrometers on ground-based telescopes (e.g., Tohline & Osterbrock 1976; Boksenberg & Netzer 1977; Foltz et al. 1981; Kollatschny et al. 1981; Schulz & Rafanelli 1981; Peterson et al. 1982; Antonucci & Cohen 1983) and in the UV with the International Ultraviolet Explorer (e.g., Ulrich et al. 1984). The interested reader is referred to Peterson (1988) for a review of the early studies of emission-line variability in AGNs.

That correlated variability of continuum and emission-line fluxes could be used to probe the structure of the broad-line region (BLR) in AGNs was recognized in the first decade of quasar research (Bahcall et al. 1972). The concept was refined in the early 1980s and has been known since as "reverberation mapping" (Blandford & McKee 1982) because the emission lines "reverberate" in response to continuum variations. Reverberation mapping has since become a standard tool for studying the structure and dynamics of the BLR (Peterson 1993, 2014). Many programs were undertaken in the 1990s, largely enabled by the proliferation of high-quality detectors on small- to medium-sized telescopes where groups of observers could obtain enough telescope time for long-term monitoring campaigns.

In its simplest form, reverberation mapping is used to measure the mean response time τ of emission lines to continuum variations, and this is interpreted as the light travel-time across the BLR radius R = cτ. By combining the measured time delay, or lag, between continuum and emission-line flux variations with some suitable measure of the emission-line width ΔV, it is possible to estimate the mass of the supermassive black hole that is the central engine of the AGN. The mass is usually expressed as

$$\begin{eqnarray}&&{M}_{\mathrm{BH}}=f\left(\displaystyle \frac{c\tau {\rm{\Delta }}{V}^{2}}{G}\right),\end{eqnarray} \tag{ 1 }$$

where G is the gravitational constant and f is an unknown scaling constant. The quantity in parentheses is often referred to as the "virial product" (VP), which has units of mass and contains only the two observables (τ and ΔV) and physical constants. All complicating factors, such as the inclination of the system or the effects of anisotropic line emission, are subsumed into the constant f. Thus, f is expected to be different for every individual AGN, but should be approximately constant for every emission line in a given AGN assuming similar geometries and dynamics of the line-emitting gas. In every case where the lags from multiple emission lines can be measured in a single object, it is found that $\tau \propto {\rm{\Delta }}{V}^{-2}$, as expected from Equation (1), suggesting that this is the case (Peterson & Wandel 1999, 2000; Kollatschny 2003; Peterson et al. 2004; Bentz et al. 2010b). It is worth reminding the reader that the scale factor f depends on which parameter is used to characterize the emission-line width, as we discuss in Section 4.

The scaling factor f can be determined for an AGN if there is an independent measurement of the black hole mass. Unfortunately, there are few cases where the black hole radius of influence is large enough that either stellar or gas dynamical modeling can also be used. At the present time, there are stellar dynamical masses (Davies et al. 2006; Onken et al. 2014) and gas dynamical masses (Hicks & Malkan 2008) for NGC 3227 and NGC 4151; these are useful for comparison purposes, but it would be unwise to attempt to calibrate the entire reverberation-based mass scale on only two objects. Instead, one of the well-known correlations between central black hole mass and properties of the host galaxies can be used. The first of these to be used to calibrate the AGN black hole mass scale was the correlation between black hole mass and host galaxy luminosity (Magorrian et al. 1998) by Laor (1998). More recently, calibration of the AGN black hole mass scale has been based on the strong correlation between the black hole mass and the velocity dispersion of host-galaxy bulge, the M_BH–σ_* relationship, which applies to both quiescent (Ferrarese & Merritt 2000; Gebhardt et al. 2000a; Tremaine et al. 2002; Gültekin et al. 2009; McConnell et al. 2011; McConnell & Ma 2013) and active galaxies (Gebhardt et al. 2000b; Ferrarese et al. 2001; Nelson et al. 2004; Onken et al. 2004; Dasyra et al. 2007; Woo et al. 2010, 2015; Graham et al. 2011; Park et al. 2012; Grier et al. 2013a; Batiste et al. 2017). With σ_* measurements now available for ∼30 AGNs from the reverberation-mapping database (Woo et al. 2015), an ensemble average $\langle f\rangle =4.47\pm 1.25$ can be computed by comparing the predicted masses from the M_BH–σ_* relationship with the observed VPs, using the line dispersion to characterize the line width. Using this prescription, black hole masses have been measured for ∼60 AGNs using reverberation mapping (see Bentz & Katz 2015 for an up-to-date compilation).

An important result from reverberation mapping is the observed "R–L" relationship between the size of the BLR and the AGN luminosity (Wandel et al. 1999; Kaspi et al. 2000, 2005; Bentz et al. 2006b, 2009a, 2013). This R–L relationship allows us to bypass resource-intensive reverberation mapping by using the luminosity to infer the BLR radius. By combining the estimate of the BLR radius with the emission-line width, we can apply Equation (1) to estimate the black hole mass (see Vestergaard et al. 2011 for a review on single-epoch M_BH estimates).

The frontier of reverberation mapping is determination of the kinematics and geometry of the BLR by examination of the emission-line response as a function of line-of-sight velocity. The ultimate goal is to either model the BLR geometry and kinematics directly (e.g., Pancoast et al. 2012, 2014; Waters et al. 2016) or to recover velocity–delay maps and model the BLR indirectly (e.g., Bentz et al. 2010a; Grier et al. 2013b). Observational results are only now beginning to appear as the technical requirements for detailed reverberation mapping are quite demanding (Horne et al. 2004).

Over the last decade, we have undertaken a new series of reverberation programs with several specific goals in mind.

• 1.  
  To increase the number of AGNs for which reverberation lags are measured for the Hβ emission line. Additional data can better constrain the R–L relationship (Bentz et al. 2013) and the AGN M_BH–σ_* relationship that underlies the reverberation-based black hole mass calibration scale (Grier et al. 2013a).
• 2.  
  To improve upon previous reverberation results. Our re-analysis of nearly all the reverberation data that existed a decade ago revealed that many of the sources would benefit from a higher sampling rate (Peterson et al. 2004).
• 3.  
  To obtain higher-quality, higher-time-resolution spectra that would enable recovery of velocity–delay maps (e.g., Grier et al. 2013b).

These programs were designed to meet the criteria described by Horne et al. (2004) to enable recovery of velocity–delay maps. They were carried out at MDM and partner observatories in 2005 (Bentz et al. 2006a, 2007; Denney et al. 2006), 2007 (Denney et al. 2009a, 2009b, 2010), and 2010 (Grier et al. 2012a, 2012b, 2013b). In addition to these ground-based programs, we carried out an intensive multiwavelength campaign on NGC 5548 known as the AGN Space Telescope and Optical Reverberation Mapping (AGN STORM) project (De Rosa et al. 2015; Edelson et al. 2015; Fausnaugh et al. 2016; Goad et al. 2016; Mathur et al. 2017; Pei et al. 2017; Starkey et al. 2017) and a concurrent optical monitoring program on additional AGNs (Fausnaugh et al. 2017). The amount of AGN reverberation data has increased dramatically over the last few years, with several other groups carrying out campaigns similar to ours (Bentz et al. 2009b, 2010a, 2014, 2016a, 2016b; Barth et al. 2011a, 2011b, 2013, 2015; Du et al. 2014, 2015, 2016; Pei et al. 2014; Wang et al. 2014). Ongoing large multi-object reverberation-mapping campaigns (King et al. 2015; Shen et al. 2015) are expected to significantly increase the number of reverberation-mapped AGNs, as well as increase redshift and luminosity ranges of the sample, especially for emission lines other than the Balmer series.

Here we report results from a campaign undertaken in early 2012. We also include additional results on NGC 3227 from 2014. We describe the observations and data analysis in Section 2. Our time-series analysis is presented in Section 3 and our black hole mass measurement is explained in Section 4. We briefly discuss and summarize our results in Section 5. When needed, we adopt a cosmological model with Ω_m = 0.30, Ω_Λ = 0.70, and H₀ = 70 km s⁻¹ Mpc⁻¹.

2.1. Target Selection

The primary objective of this campaign is to determine the kinematics and structure of the BLR in a few well-studied bright AGNs. In particular, we are re-examining NGC 3227, NGC 3516, and NGC 5548 from Denney et al. (2009a), for which cross-correlation of individual velocity bins suggested gross kinematics of outflow, infall, and rotation/virialization, respectively. As we discuss here and elsewhere, these results need to be checked and more thoroughly characterized. We also included in our observing program NGC 4151, for which the best reverberation data are from a weather-abbreviated campaign in 2005 (Bentz et al. 2006a). In addition to these primary targets, we added a few sources that could only be observed for part of our campaign on account of their location in the sky. Sources included in the 2012 campaign were Mrk 374, Mrk 382, Mrk 478, Mrk 618, and Mrk 704. Because of the shorter monitoring period, the failure rate for these secondary sources was high, with only Mrk 704 yielding data useful for reverberation purposes.

The properties of the sources studied in this paper are summarized in Table 1. Both NGC 3227 and NGC 4151 are too close for redshift-based distances to be reliable. NGC 3227 is interacting with an elliptical companion, NGC 3226, which has a surface-brightness fluctuation distance of 23.5 Mpc (Tonry et al. 2001). We therefore adopt this as the distance to NGC 3227. In the case of NGC 4151, we are currently working on a Cepheid-based distance, but here we use the distance of 13.9 Mpc adopted by Onken et al. (2014) in their recent stellar dynamical study (although this distance is derived from Hubble's law).

Table 1.  Source Properties

Object	z	DL	$\mathrm{log}\,\lambda {L}_{5100\mathring{\rm A} }$	$\mathrm{log}\,\lambda {L}_{\mathrm{host}}$	E(B−V)
		(Mpc)	(erg s−1)	(erg s−1)	(mag)
(1)	(2)	(3)	(4)	(5)	(6)
Mrk 704	0.0292	128.0	43.72	43.27	0.03
NGC 3227	0.0038	23.5	42.74	42.48	0.02
NGC 3516	0.0088	38.1	43.29	43.21	0.04
NGC 4151	0.0033	13.9	42.61	42.37	0.02
NGC 5548	0.0171	74.5	43.45	43.20	0.02

Note. Column 2 is taken from the NASA Extragalactic Database. Column 3 gives the luminosity distance in a consensus cosmology, except for NGC 3227 and NGC 4151 as explained in the text (see Section 2.1). Column 4 gives the observed luminosity (corrected for Galactic extinction), calculated from the observed 5100 Å rest-frame light curve and Column 3 and is corrected for the starlight contribution which is given in Column 5. Column 6 gives the Galactic reddening value from Schlafly & Finkbeiner (2011).

Download table as:  ASCIITypeset image

2.2. Observations

2.3. Data Processing and Light Curves

The procedures we followed for reducing the data, producing calibrated light curves, and assessing uncertainties are described in detail by Fausnaugh et al. (2017). We provide a brief recapitulation here.

2.3.1. Spectroscopy

Two-dimensional spectra were reduced using standard IRAF tasks to deal with bias, flat field, sky subtraction, and wavelength calibration. An extraction window of 12'' was used throughout. Cosmic ray removal was done using LA Cosmic (van Dokkum 2001). Flux calibration relied on observations of standard stars, usually Feige 34 and/or BD +33°2642 (Oke 1990).

We used the narrow [O iii] λ5007 emission line as an internal flux standard for both relative and absolute calibration. While narrow emission lines have been found to vary on long timescales (years to decades; e.g., Peterson et al. 2013), they are effectively constant in flux on BLR reverberation timescales (days to months). We identified all the individual spectra where the observer reported "clear" or "photometric" observing conditions. The [O iii] λ5007 flux was measured, and from these a mean and standard deviation was computed. Outliers greater than 3σ from the mean were rejected and the mean and standard deviation were recomputed. The number of observations used for the calibration is given in Column (2) of Table 2 and the adopted [O iii] λ5007 fluxes appear in Column (3) of the same table. This provides the absolute flux calibration for the spectrophotometric observations. We note that the [O iii] λ5007 flux in NGC 5548 is in good agreement with the preliminary measurement we presented earlier (Peterson et al. 2013).

Table 2.  [O iii] λ 5007 Flux Calibration

Object	No.	F([O iii]λ5007)	Percent
	Photometric	(10−13 erg s−1 cm−2)	Scatter
(1)	(2)	(3)	(4)
Mrk 704	21	1.31 ± 0.03	0.40
NGC 3227	24	7.81 ± 0.16	0.19
NGC 3516	21	4.58 ± 0.07	0.27
NGC 4151	20	107 ± 2	0.15
NGC 5548	21	4.91 ± 0.08	0.18

Note. Column 2 gives the number of nights with clear and stable conditions and judged to be photometric. Each object had three observations per night, which were used to calculate the narrow [O iii]λ5007 line flux. The line flux and its uncertainty are given in Column 3. Column 4 gives the fractional variation of the [O iii]λ5007 line light curve, which serves as an estimate of the night-to-night calibration error.

Download table as:  ASCIITypeset image

For each AGN, the spectra with the highest signal-to-noise ratios (S/Ns) and no obvious flaws are combined to form a reference spectrum, which is scaled to have the adopted [O iii] λ5007 flux. We then place all the individual spectra on the same relative flux by scaling each spectrum to the adopted [O iii] λ5007 flux. This is done using a Markov chain Monte Carlo code called mapspec (Fausnaugh 2017) which adjusts the flux, wavelength shift, and resolution of each individual spectrum to match that of the reference spectrum, as measured by the [O iii] λ5007 emission line profiles. This affords a factor of several improvement over the long-used method of van Groningen & Wanders (1992), as assessed by the root-mean-square (rms) scatter of the [O iii] λ5007 flux across the full time series.

Once flux calibration is complete, we combine all N spectra for each object to form a weighted mean spectrum

$$\begin{eqnarray}&&\langle F(\lambda )\rangle =\displaystyle \frac{{\displaystyle \sum }_{i=1}^{N}F(\lambda ,{t}_{i})/{\sigma }^{2}(\lambda ,{t}_{i})}{{\displaystyle \sum }_{i=1}^{N}1/{\sigma }^{2}(\lambda ,{t}_{i})},\end{eqnarray} \tag{ 2 }$$

where F(λ, t_i) is the flux at epoch t_i and σ(λ, t_i) is the associated uncertainty. We also form an rms residual spectrum

$$\begin{eqnarray}&&{\sigma }_{\mathrm{rms}}(\lambda )={\left\{\displaystyle \frac{1}{N-1}\displaystyle \sum _{i=1}^{N}{[F(\lambda ,{t}_{i})-\langle F(\lambda )\rangle ]}^{2}\right\}}^{1/2}.\end{eqnarray} \tag{ 3 }$$

The mean and rms spectra for our sources are shown in Figures 1–6. The rms spectrum is of special value in this context, since the constant components (e.g., host-galaxy starlight, narrow emission lines) vanish, isolating the variable part of the spectrum. However, the total variability power contains contributions not only from intrinsic variability, but also from statistical fluctuations and/or measurement errors. We therefore attempt to isolate the intrinsic variability by minimizing the negative log-likelihood

$$\begin{eqnarray}\begin{array}{rcl}-2\mathrm{ln}{ \mathcal L } & = & \displaystyle \sum _{i=1}^{N}\displaystyle \frac{{[F(\lambda ,{t}_{i})-\hat{F}(\lambda )]}^{2}}{{\sigma }^{2}(\lambda ,{t}_{i})+{\sigma }_{\mathrm{var}}^{2}(\lambda )}\\ & & +\ \displaystyle \sum _{i=1}^{N}\mathrm{ln}[{\sigma }^{2}(\lambda ,{t}_{1})+{\sigma }_{\mathrm{var}}^{2}(\lambda )],\end{array}\end{eqnarray} \tag{ 4 }$$

where $\hat{F}(\lambda )$ is the optimal average weighted by ${\sigma }^{2}({t}_{i})+{\sigma }_{\mathrm{var}}^{2}$ and σ_var(λ) is the intrinsic variability. We solve simultaneously for $\hat{F}(\lambda )$ and σ_var(λ). Our estimator of the intrinsic variability σ_var(λ) is also shown in Figures 1–6.

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

2.3.2. Imaging

Independent continuum light curves were generated for each bandpass for each set of imaging data using the image subtraction software package ISIS (Alard & Lupton 1998; Alard 2000). We followed the procedures as described by Shappee & Stanek (2011). First, we aligned the images using Sexterp (Siverd et al. 2012). We then created a reference image with ISIS for each AGN field by combining the images with the best seeing and lowest background counts; typically we used 5–15 images to construct the reference image. ISIS convolves the images of each AGN with a convolution kernel that is allowed to vary across the field in order to transform all the images to the same point-spread function (PSF) and background level. The reference image was convolved to match each individual frame and ISIS then subtracted each image from the convolved reference frame. The fluxes of the AGN and control stars to estimate errors were determined by fitting a PSF-weighted aperture over each source, thus producing a differential light curve.

2.3.3. Construction of Light Curves

A spectroscopic continuum light curve, nominally at ∼5100 Å in the rest frame of each AGN, is formed by averaging the flux densities over the wavelength ranges given in Table 3 and shown as a shaded region in Figures 1–6. Our final continuum light curves are constructed by merging the differential V-band light curves with the 5100 Å spectroscopic light curve by scaling multiplicatively (to match the variations) and shifting additively (to account for the different mean flux levels in each reference image) the differential continuum light curves. We found that the uncertainties on the differential light curves are systematically too small because ISIS takes into account only Poisson errors. To bypass this problem, we rescaled the errors based on measurements of other stars in the field of view, as described in detail by Fausnaugh et al. (2016).

Table 3.  Observed-frame Integration Windows

Object	5100 Å	Hβ	[O iii] λ5007
	(Å)	(Å)	(Å)
Mrk 704	5250–5270	4910–5122	5138–5168
NGC 3227	5110–5130	4812–4942	5005–5047
NGC 3516	5128–5170	4775–4960	5032–5066
NGC 4151	5110–5140	4815–4955	4998–5055
NGC 5548	5179–5210	4830–5052	5070–5110

Download table as:  ASCIITypeset image

The emission-line light curves are generated by interpolating a simple linear continuum underneath the emission lines using the windows given in Table 4 and integrating the flux above this continuum between the limits given in Table 3 and illustrated in Figures 1–6. These measurements are fairly crude, but are intended to capture the emission-line variations as opposed to all the emission-line flux. A more sophisticated treatment is deferred to a future paper. We estimate the uncertainties using a local linear interpolation method described in detail by Fausnaugh et al. (2017), which rescales the statistical uncertainties of the light curves so that they are consistent with the observed night-to-night scatter.

Table 4.  Observed-frame Continuum Fitting Windows

Object	Line Side	Hβ	[O iii] λ4959	[O iii] λ5007
		(Å)	(Å)	(Å)
Mrk 704	Blue	4890–4910	5080–5090	5128–5137
	Red	5235–5245	5113–5120	5169–5175
NGC 3227	Blue	4801–4811	4950–4961	5000–5005
	Red	4941–4950	4992–5002	5047–5055
NGC 3516	Blue	4743–4752	4970–4980	5022–5032
	Red	5128–5170	5014–5030	5067–5075
NGC 4151	Blue	4510–4520	4950–4955	4993–4998
	Red	5110–5140	4990–5000	5055–5070
NGC 5548	Blue	4535–4545	5017–5027	5060–5070
	Red	5136–5159	5058–5070	5110–5120

Download table as:  ASCIITypeset image

The final continuum and emission line light curves are given in Tables 5 and 6 respectively. All the light curves are shown in Figures 1–6. The statistical properties of the light curves are summarized in Table 7, including the number of observations N_obs, the uncertainty rescaling factor, median cadence Δt_med, mean flux $\langle F\rangle$, $\langle {\rm{S}}/{\rm{N}}\rangle$, excess variance ${F}_{\mathrm{var}}={\sigma }_{\mathrm{var}}/\hat{F}$, where ${\sigma }_{\mathrm{var}}$ and $\hat{F}$ are determined in the same way as in Equation (4) (after integrating over λ to produce the light curves), and the significance

$$\begin{eqnarray}&&{({\rm{S}}/{\rm{N}})}_{\mathrm{var}}=\displaystyle \frac{{\sigma }_{\mathrm{var}}}{\overline{\sigma }\sqrt{2/{N}_{\mathrm{obs}}}}\end{eqnarray} \tag{ 5 }$$

at which variability is detected, where $\overline{\sigma }$ is the mean measurement uncertainty. Further details can be found in Fausnaugh et al. (2017).

Table 5.  Continuum Light Curves

Mrk 704		NGC 3227		NGC 3227 (2014)		NGC 3516		NGC 4151		NGC 5548
HJDa	Fλb	HJDa	Fλb	HJDa	Fλb	HJDa	Fλb	HJDa	Fλb	HJDa	Fλb
5932.26	4.29 ± 0.05 M	5933.82	14.74 ± 0.22 M	6645.61	10.81 ± 0.22 W1	5932.39	19.16 ± 0.20 M	5931.51	35.52 ± 1.03 C	5931.55	8.46 ± 0.22 C
5933.57	4.27 ± 0.02 W1	5935.32	15.25 ± 0.32 M	6646.60	10.89 ± 0.23 W1	5933.59	18.57 ± 0.05 W1	5932.38	37.03 ± 0.55 M	5932.44	8.83 ± 0.15 M
5933.75	4.26 ± 0.11 M	5935.55	15.67 ± 0.02 W1	6647.62	11.16 ± 0.20 W1	5933.86	18.71 ± 0.29 M	5933.91	36.82 ± 0.42 M	5933.62	8.60 ± 0.05 W1
5935.26	4.30 ± 0.08 M	5936.81	15.08 ± 0.26 M	6648.59	10.98 ± 0.14 W1	5934.89	19.15 ± 0.28 M	5935.38	34.92 ± 0.30 M	5933.94	8.50 ± 0.10 M
5939.50	4.44 ± 0.02 W1	5937.60	16.01 ± 0.15 A1	6650.51	11.06 ± 0.16 W1	5935.85	18.43 ± 0.31 M	5935.90	34.44 ± 1.32 M	5935.62	8.81 ± 0.03 W1
5940.75	4.32 ± 0.06 M	5937.80	14.34 ± 0.63 M	6653.59	12.52 ± 0.34 W1	5936.47	19.03 ± 0.04 W1	5936.51	34.22 ± 0.19 W1	5935.93	8.57 ± 0.11 M
5942.35	4.30 ± 0.03 W1	5938.80	15.46 ± 0.47 M	6655.50	13.01 ± 0.37 W1	5936.86	18.68 ± 0.31 M	5937.40	33.10 ± 0.26 M	5936.54	8.75 ± 0.03 W1
5943.56	4.41 ± 0.02 W1	5940.55	15.64 ± 0.14 A1	6656.51	13.06 ± 0.37 W1	5937.85	18.30 ± 0.31 M	5937.56	33.65 ± 1.33 A1	5936.93	8.55 ± 0.12 M
5944.76	4.27 ± 0.06 M	5940.80	15.60 ± 0.25 M	6661.93	13.70 ± 0.16 M	5938.14	20.18 ± 0.08 A2	5938.65	33.64 ± 1.33 A1	5937.93	8.49 ± 0.11 M
5945.74	4.30 ± 0.06 M	5941.48	15.65 ± 0.14 A1	6662.48	13.07 ± 0.26 W1	5939.63	19.37 ± 0.05 W1	5939.62	33.30 ± 1.32 A1	5938.59	8.34 ± 0.24 A1
5946.75	4.34 ± 0.05 M	5942.47	15.88 ± 0.15 A1	6663.50	13.00 ± 0.41 W1	5939.82	19.46 ± 0.32 M	5939.65	33.52 ± 0.21 W1	5939.69	8.88 ± 0.26 A1
5947.75	4.30 ± 0.04 M	5942.54	15.93 ± 0.02 W1	6663.90	13.75 ± 0.13 M	5942.59	19.69 ± 0.07 W1	5940.59	33.46 ± 1.32 A1	5940.43	8.03 ± 0.05 M
5949.33	4.55 ± 0.10 M	5943.45	16.01 ± 0.15 A1	6664.88	13.94 ± 0.14 M	5943.60	20.06 ± 0.08 W1	5940.89	33.57 ± 0.35 M	5940.61	8.06 ± 0.23 A1
5951.75	4.63 ± 0.04 M	5943.60	16.39 ± 0.03 W1	6665.42	14.26 ± 0.21 W1	5944.87	19.93 ± 0.30 M	5941.58	33.46 ± 1.32 A1	5941.60	8.00 ± 0.23 A1
5952.45	4.48 ± 0.03 W1	5944.01	16.14 ± 0.69 M	6666.93	14.02 ± 0.17 M	5945.84	19.77 ± 0.32 M	5942.57	34.30 ± 1.35 A1	5942.59	7.91 ± 0.23 A1
5952.76	4.79 ± 0.08 M	5944.44	15.90 ± 0.15 A1	6666.94	13.95 ± 0.14 M	5946.89	19.50 ± 0.30 M	5942.62	33.43 ± 0.23 W1	5942.62	7.92 ± 0.04 W1
5953.45	4.52 ± 0.03 W1	5945.31	16.14 ± 0.14 M	6670.45	13.35 ± 0.17 W1	5947.11	19.43 ± 0.08 A2	5943.56	34.61 ± 1.36 A1	5943.58	7.72 ± 0.22 A1
5953.75	4.46 ± 0.05 M	5945.44	15.86 ± 0.14 A1	6671.63	13.66 ± 0.21 C	5948.65	20.10 ± 0.07 W1	5943.62	33.74 ± 0.17 W1	5943.59	8.00 ± 0.04 W1
5955.40	4.51 ± 0.03 W1	5946.49	16.39 ± 0.16 A1	6671.89	12.96 ± 0.30 FWO	5948.94	20.31 ± 0.31 M	5944.53	32.46 ± 1.29 A1	5945.43	7.95 ± 0.05 M
5955.75	4.54 ± 0.08 M	5947.33	15.84 ± 0.14 M	6672.45	13.09 ± 0.17 W1	5949.86	20.08 ± 0.32 M	5945.40	31.77 ± 0.48 M	5947.45	8.04 ± 0.05 M
5956.24	4.57 ± 0.01 W1	5948.52	16.38 ± 0.16 A1	6672.91	13.19 ± 0.28 FWO	5950.63	20.60 ± 0.08 W1	5945.54	30.85 ± 1.25 A1	5948.64	8.21 ± 0.24 A1
5956.75	4.54 ± 0.05 M	5948.65	16.36 ± 0.03 W1	6674.00	13.27 ± 0.27 FWO	5951.86	20.28 ± 0.31 M	5946.93	31.13 ± 0.38 M	5948.65	8.26 ± 0.04 W1
5957.76	4.58 ± 0.09 M	5949.46	16.20 ± 0.15 A1	6679.44	13.44 ± 0.13 W1	5952.61	20.33 ± 0.09 W1	5947.90	32.61 ± 0.71 M	5949.94	8.14 ± 0.09 M
5958.60	4.57 ± 0.03 W1	5949.82	16.58 ± 0.28 M	6679.91	13.83 ± 0.32 FWO	5952.87	21.06 ± 0.33 M	5948.61	32.28 ± 1.29 A1	5950.62	8.01 ± 0.23 A1
5958.77	4.52 ± 0.04 M	5950.63	16.22 ± 0.03 W1	6682.40	14.58 ± 0.27 W1	5953.57	20.52 ± 0.10 W1	5948.65	32.69 ± 0.24 W1	5950.64	8.11 ± 0.05 W1
...	...	...	...	...	...	...	...	...	...	...	...

Notes. The alpha-numeric codes in the F_λ columns indicate the contributing observatory: M is MDM, W1 is Wise, W2 is West Mountain, A1 is Asiago, A2 is ASAS, C is Crimean Astrophysical, and F is Fountainwood.

^aDays − 2,450,000. ^b10⁻¹⁵ erg s⁻¹ cm⁻² Å⁻¹.

Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.

Download table as:  DataTypeset image

Table 6.  Hβ Light Curves

Mrk 704		NGC 3227		NGC 3227 (2014)		NGC 3516		NGC 4151		NGC 5548
HJDa	Fb	HJDa	Fb	HJDa	Fb	HJDa	Fb	HJDa	Fb	HJDa	Fb
5932.26	2.92 ± 0.07 M	5933.82	4.90 ± 0.20 M	6661.93	3.48 ± 0.15 M	5932.39	5.58 ± 0.08 M	5931.51	48.24 ± 1.20 C	5931.55	7.08 ± 0.15 C
5933.75	2.94 ± 0.12 M	5935.32	4.89 ± 0.08 M	6663.90	3.61 ± 0.05 M	5933.86	5.62 ± 0.07 M	5932.38	47.52 ± 0.55 M	5932.44	7.19 ± 0.10 M
5935.26	2.93 ± 0.08 M	5936.81	4.93 ± 0.04 M	6664.88	3.74 ± 0.06 M	5934.89	5.37 ± 0.21 M	5933.91	48.33 ± 0.95 M	5933.94	7.12 ± 0.05 M
5940.75	3.02 ± 0.05 M	5937.60	5.38 ± 0.16 A1	6666.93	3.90 ± 0.09 M	5935.85	5.55 ± 0.27 M	5935.38	46.72 ± 0.33 M	5935.93	6.92 ± 0.10 M
5944.76	2.99 ± 0.03 M	5937.80	5.05 ± 0.11 M	6666.94	3.97 ± 0.09 M	5936.86	5.48 ± 0.12 M	5935.90	46.44 ± 0.24 M	5936.93	6.89 ± 0.15 M
5945.74	2.96 ± 0.02 M	5938.80	5.39 ± 0.55 M	6777.72	4.30 ± 0.05 M	5937.85	5.51 ± 0.26 M	5937.40	45.51 ± 0.42 M	5937.93	6.89 ± 0.32 M
5946.75	2.96 ± 0.02 M	5940.55	4.99 ± 0.15 A1	6778.70	4.13 ± 0.07 M	5939.82	5.62 ± 0.20 M	5937.56	46.44 ± 1.12 A1	5938.59	6.90 ± 0.16 A1
5947.75	2.97 ± 0.02 M	5940.80	4.95 ± 0.17 M	6779.71	4.19 ± 0.05 M	5944.87	5.75 ± 0.07 M	5938.65	46.52 ± 1.13 A1	5939.69	6.36 ± 0.14 A1
5949.33	2.88 ± 0.03 M	5941.48	5.00 ± 0.15 A1	6780.71	4.14 ± 0.05 M	5945.84	5.68 ± 0.07 M	5939.62	43.57 ± 1.04 A1	5940.43	6.60 ± 0.14 M
5951.75	2.95 ± 0.05 M	5942.47	5.24 ± 0.15 A1	6781.70	4.08 ± 0.07 M	5946.89	5.95 ± 0.11 M	5940.59	44.43 ± 1.06 A1	5940.61	6.66 ± 0.15 A1
5952.76	2.89 ± 0.06 M	5943.45	5.38 ± 0.16 A1	6782.71	4.04 ± 0.06 M	5948.94	5.48 ± 0.09 M	5940.77	45.29 ± 0.74 M	5941.60	6.52 ± 0.15 A1
5953.75	2.93 ± 0.04 M	5944.01	5.27 ± 0.19 M	6783.70	4.00 ± 0.07 M	5949.86	5.80 ± 0.12 M	5941.58	45.53 ± 1.10 A1	5942.59	6.53 ± 0.15 A1
5955.75	2.82 ± 0.05 M	5944.44	5.26 ± 0.15 A1	6784.77	4.05 ± 0.09 M	5951.86	6.09 ± 0.08 M	5942.57	45.37 ± 1.09 A1	5943.58	6.49 ± 0.15 A1
5956.75	3.00 ± 0.05 M	5945.31	5.52 ± 0.06 M	6785.70	3.86 ± 0.06 M	5952.87	6.03 ± 0.17 M	5943.56	46.06 ± 1.11 A1	5945.43	6.17 ± 0.07 M
5957.76	2.94 ± 0.07 M	5945.44	5.23 ± 0.15 A1	6786.71	3.96 ± 0.06 M	5953.85	5.97 ± 0.16 M	5944.53	43.51 ± 1.04 A1	5947.45	6.19 ± 0.05 M
5958.77	2.88 ± 0.03 M	5946.49	5.79 ± 0.17 A1	6787.69	3.96 ± 0.06 M	5955.01	6.41 ± 0.18 M	5945.40	44.06 ± 0.32 M	5948.64	5.94 ± 0.13 A1
5960.27	2.97 ± 0.02 M	5947.33	5.58 ± 0.03 M	6788.70	4.03 ± 0.07 M	5955.86	6.13 ± 0.11 M	5945.54	43.99 ± 1.05 A1	5949.94	5.86 ± 0.10 M
5961.75	2.97 ± 0.14 M	5948.52	5.78 ± 0.17 A1	6789.71	4.06 ± 0.07 M	5956.86	6.15 ± 0.08 M	5946.93	43.56 ± 0.29 M	5950.62	6.14 ± 0.13 A1
5962.77	2.96 ± 0.05 M	5949.46	5.58 ± 0.16 A1	6790.68	4.03 ± 0.08 M	5957.87	6.53 ± 0.11 M	5947.90	44.77 ± 1.17 M	5951.54	5.97 ± 0.13 A1
5963.50	2.98 ± 0.08 M	5949.82	5.57 ± 0.09 M	6791.69	3.81 ± 0.07 M	5958.88	6.05 ± 0.23 M	5948.61	42.95 ± 1.02 A1	5951.95	6.21 ± 0.08 M
5968.27	2.98 ± 0.07 M	5951.45	5.62 ± 0.16 A1	6792.69	3.90 ± 0.07 M	5959.88	6.30 ± 0.07 M	5949.03	43.95 ± 0.32 M	5953.45	6.08 ± 0.07 M
5973.69	3.10 ± 0.05 M	5951.80	5.75 ± 0.04 M	6793.70	3.94 ± 0.07 M	5960.86	6.10 ± 0.09 M	5949.55	42.68 ± 1.01 A1	5953.69	6.04 ± 0.13 A1
5978.76	3.16 ± 0.05 M	5952.70	5.78 ± 0.17 A1	6795.69	4.00 ± 0.10 M	5961.86	6.48 ± 0.18 M	5949.91	42.63 ± 0.23 M	5955.52	6.28 ± 0.12 C
5979.76	3.26 ± 0.08 M	5953.32	5.68 ± 0.08 M	6798.71	4.19 ± 0.06 M	5962.88	6.17 ± 0.16 M	5950.59	44.28 ± 1.06 A1	5956.54	5.99 ± 0.07 M
5980.73	3.24 ± 0.06 M	5955.81	5.60 ± 0.05 M	6799.70	4.07 ± 0.06 M	5964.24	6.07 ± 0.46 M	5951.52	43.22 ± 1.03 A1	5957.92	6.20 ± 0.14 M
...	...	...	...	...	...	...	...	...	...	...	...

Notes. The alpha-numeric codes in the F_λ columns indicate the contributing observatory: M is MDM, A1 is Asiago, and C is Crimean Astrophysical.

^aDays − 2,450,000. ^b10⁻¹³ erg s⁻¹ cm⁻².

Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.

Download table as:  DataTypeset image

Table 7.  Light-curve Properties

Object	Light Curve	Nobs	Δtmed	Uncertainty	$\langle F\rangle$	$\langle {\rm{S}}/{\rm{N}}\rangle$	σvar	(S/N)var	rmax
			(days)	Rescaling Factor
(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
Mrk 704	5100 Å	97	1.01	1.53	4.81	205.5	0.06	88.1	...
	Hβ	72	1.01	1.50	3.30	89.7	0.07	40.0	0.92 ± 0.02
NGC 3227	5100 Å	94	1.01	2.53	15.60	175.0	0.02	27.5	...
	Hβ	84	1.02	1.97	4.94	74.8	0.08	40.3	0.69 ± 0.06
NGC 3227 (2014)	5100 Å	56	1.00	1.28	13.60	99.8	0.05	24.2	...
	Hβ	34	1.00	1.41	4.14	55.1	0.04	9.0	0.77 ± 0.05
NGC 3516	5100 Å	109	1.03	1.64	19.70	226.5	0.03	54.9	...
	Hβ	77	1.01	2.54	5.96	46.8	0.07	20.1	0.76 ± 0.04
NGC 4151	5100 Å	119	1.01	3.09	32.20	159.3	0.11	131.7	...
	Hβ	97	1.01	3.59	40.10	111.7	0.13	103.7	0.94 ± 0.01
NGC 5548	5100 Å	112	1.03	1.84	7.94	156.3	0.07	87.6	...
	Hβ	91	1.03	2.61	5.77	64.2	0.10	44.7	0.74 ± 0.04

Note. Column 3 gives the number of observations in each light curve. Column 4 gives the median cadence. Column 5 gives the rescaling factor by which the statistical uncertainties are multiplied to account for additional systematic errors (see Section 2.3.3). Column 6 gives the mean flux level of each light curve. The rest-frame 5100 Å continuum light curves are in units of 10⁻¹⁵ erg cm⁻² s⁻¹ Å⁻¹, and the emission line light curves are in units of 10⁻¹³ erg cm⁻² s⁻¹. Column 7 gives the mean signal-to-noise ratio $\langle {\rm{S}}/{\rm{N}}\rangle$. Column 8 gives the rms fractional variability defined in Equation (4). Column 9 gives the approximate S/N at which we detect variability (see Section 2.2.2). Column 10 gives the maximum value of the interpolated cross-correlation function (ICCF) (see Section 3).

Download table as:  ASCIITypeset image

3.1. Mean Emission-line Lags

Our initial goal is to determine the mean timescale for the response of the Hβ emission line to continuum variations, which we later use to determine the mass of the central black hole.

The time series analysis is carried out using two common methodologies, interpolated cross-correlation (Gaskell & Sparke 1986; Gaskell & Peterson 1987; White & Peterson 1994; Peterson et al. 1998, 2004) and the stochastic process modeling algorithm JAVELIN⁶² (Zu et al. 2011). A more complete description of how we have employed these methods for such analysis is provided by Fausnaugh et al. (2017).

Results of the time-series analysis are given in Table 8 and shown graphically in the right-hand panels of Figures 7–11. It is interesting to notice that the three AGNs from Denney et al. (2009a), re-observed in this program, all have shorter lags than they did in 2007. In the case of NGC 3516, the factor-of-two decrease in the Hβ lag is consistent with the factor-of-four decrease in the AGN luminosity and the expected scaling relation ${R}_{\mathrm{BLR}}\propto {L}_{\mathrm{AGN}}^{1/2}$. In the case of NGC 3227, the Hβ lag also decreased by a factor of two from 2007, but the AGN luminosity is in fact slightly higher in 2012 and 2014. In 2012, NGC 5548 had been in a prolonged faint state for a few years (Peterson et al. 2013) and by 2013–14 heavy internal absorption became an important factor (Kaastra et al. 2014; De Rosa et al. 2015). In both 2012 and 2014 (Pei et al. 2017), the Hβ lag is found to be surprisingly short given the AGN luminosity at the time. The implications of this are not yet clear, although it appears that increased absorption within the BLR plays some role. The 2014 data on NGC 3227 are quite marginal, and JAVELIN was unable to converge on a solution for the lag. The ICCF analysis, however, shows consistency with the 2012 results. In the case of NGC 4151, the Hβ lag is in good agreement with that obtained by Bentz et al. (2006a).

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

Table 8.  Rest-frame Hβ Lags

Object	τcent	τpeak	τJAV
	(days)	(days)	(days)
(1)	(2)	(3)	(4)
Mrk 704	${12.65}_{-2.14}^{+1.49}$	${14.87}_{-2.45}^{+5.85}$	${14.32}_{-1.06}^{+0.87}$
NGC 3227 (2012)	${1.29}_{-1.27}^{+1.56}$	${1.74}_{-1.59}^{+1.69}$	${2.29}_{-0.20}^{+0.23}$
NGC 3227 (2014)	${2.58}_{-1.31}^{+1.20}$	${2.80}_{-1.60}^{+1.0}$	...
NGC 3516	${5.74}_{-2.04}^{+2.26}$	${4.24}_{-3.93}^{+2.16}$	${8.27}_{-0.64}^{+1.12}$
NGC 4151	${6.82}_{-0.57}^{+0.48}$	${6.50}_{-1.39}^{+0.99}$	${6.58}_{-0.22}^{+0.19}$
NGC 5548	${2.83}_{-0.96}^{+0.88}$	${2.66}_{-1.55}^{+1.06}$	${3.66}_{-0.52}^{+0.53}$

Note. Columns 2 and 3 give the centroids and peaks, respectively, of the interpolated cross correlation functions. The uncertainties give the central 68% confidence intervals of the cross-correlation centroid distribution (Peterson et al. 1998). Column 4 gives the lag fit by JAVELIN. The uncertainties give the central 68% confidence intervals of the JAVELIN posterior lag distributions. All lags are relative to the 5100 Å continuum light curve and corrected to the rest frame.

Download table as:  ASCIITypeset image

We note in passing that we also attempted to measure the variations of the He ii λ4686 line, which is clearly seen in the rms residual spectra of each source (Figures 1–6). Unfortunately this is a weak, low-contrast feature, and the measurements are very uncertain on account of the difficulties in defining the underlying continuum. Contamination of the spectra by the host-galaxy starlight is a significant problem in low-luminosity AGNs, and it needs to be modeled and subtracted for a reliable He ii measurement. We defer this to a future paper.

3.2. Velocity-resolved Lags

The individual spectra are of high enough S/N and sufficiently well-sampled in time that we can also divide each emission line into line-of-sight velocity bins to see if there are any indications of gross kinematic signature and, in the cases of NGC 3227, NGC 3516, and NGC 5548, compare these results with those obtained by Denney et al. (2009a). This is not a foolproof method of discerning the velocity field of the BLR as experience has shown that reverberation effects are quite subtle, and attempting to characterize an individual velocity bin with a single average lag could be misleading. While we must interpret the results cautiously, detection of a velocity-dependent lag signature identifies good candidates for more ambitious attempts to determine the BLR structure and velocity field by either forward modeling (Pancoast et al. 2012, 2014; Grier et al. 2017) or regularization (Horne et al. 2004; Bentz et al. 2010a; Grier et al. 2013b; Skielboe et al. 2015). The results of measuring velocity-dependent lags are shown in the lower panels of Figures 7–11 in a format that can be easily compared with Figure 3 of Denney et al. (2009a) in Figures 12–16. We comment on each source individually:

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

Mrk 704. Figure 7 shows that the highest-velocity blueshifted and redshifted bins have large lag uncertainties, so we will disregard these. The remaining bins (Figure 12) show a local lag minimum around line center (V = 0). At higher redshifted velocities, the lags increase to a maximum at ∼2000 km s⁻¹, then become smaller in the far wings. On the blueshifted side, the lags also increase, but we do not see a turnover toward smaller lags at higher velocity. A similar pattern with relatively small lags at the line center compared to the wings is seen in NGC 5548 in 2014 in Hβ (Pei et al. 2017) and also probably in Lyα (De Rosa et al. 2015). The BLR velocity field in Mrk 704 may well have multiple components, and requires more sophisticated modeling.

NGC 3227. As with Mrk 704, the highest-velocity bins have large errors and should be disregarded (Figure 8). The remaining bins show a pattern that suggests a virialized BLR (Figure 13), with large lags at the line center and shorter lags at higher positive and negative velocities. These can be compared with results from Denney et al. (2009a), which do not show a decrease in the lag at higher positive velocities. Lower lags at high negative velocity might be interpreted as evidence for outflow. There is no strong evidence for outflow in the 2012 data. Again, more sophisticated modeling will clarify the situation. As noted earlier, the 2014 data on this source are marginal and are not included in this analysis. We note that a very similar dependence of lag on velocity bin is seen from an independent RM campaign from 2017 (M.S. Brotherton 2018, private communication).

NGC 3516. In 2007 (Denney et al. 2009a), the highest positive velocities in the Hβ emission line showed the shortest lags, with the lags steadily increasing toward the line center and continuing to increase slightly to higher negative velocity (Figure 14). This behavior could be interpreted as an infall signature. In 2012, at least in the core of the line, this trend seems to be reversed.

NGC 4151. On account of the brightness and favorable variability characteristics of NGC 4151 during this campaign, the results for this AGN are superb. The uncertainties in the lag for each velocity bin are quite small (Figure 10) and there is a very clear virial-like pattern where the largest lags are seen at the lowest velocities (Figure 15).

NGC 5548. Due to less favorable variability characteristics in 2012, the NGC 5548 results are not as clear as they were in either 2007 (Denney et al. 2009a) or in 2014 (Pei et al. 2017); the uncertainties in each velocity bin are comparatively large (Figure 11). The pattern as a function of wavelength seems to be quite similar to the complex pattern observed in 2014 (see Figure 10 of Pei et al. 2017) as well as in 2015 Lu et al. (2016). This is also similar to Mrk 704 (Figure 12), and possibly indicates a multicomponent BLR.

In order to compute the mass of the central black hole from Equation (1), we need to characterize the line width ΔV in addition to the mean emission-line lag τ. The two line width measures commonly used for this are full width at half maximum (FWHM) and the line dispersion

$$\begin{eqnarray}&&{\sigma }_{\mathrm{line}}={\left[\displaystyle \frac{\int {v}^{2}P(v){dv}}{\int P(v){dv}}\right]}^{1/2},\end{eqnarray} \tag{ 6 }$$

which is the square root of the second moment of the line. The integral is over the line profile P(v) as a function of line-of-sight (Doppler) velocity. There are practical advantages and disadvantages to each of these. The FWHM is usually trivial to measure, but presents problems when the data are noisy or the profiles are complex. The line dispersion, on the other hand, requires attention to blending with other features, but is computationally well-defined for any profile. The two measures are not interchangeable, as their ratio varies with line shape, which is correlated with line width. There are compelling, but not conclusive, arguments that line dispersion is the better parameter for computing masses (Denney et al. 2013; Peterson 2014), so we use σ_line in our mass calculations, but report both measures for both the mean and rms spectra in Table 9.

Table 9.  Rest-frame Hβ Velocity Measurements

		rms Spectrum		Mean Spectrum
Object	Line	σline	FWHM	σline	FWHM	Smoothing Width
		(km s−1)	(km s−1)	(km s−1)	(km s−1)	(km s−1)
(1)	(2)	(3)	(4)	(5)	(6)	(7)
Mrk 704	Hβ	${1860}_{-130}^{+108}$	${3406}_{-240}^{+310}$	${2650}_{-3}^{+4}$	${3502}_{-30}^{+32}$	294
NGC 3227 (2012)	Hβ	${1368}_{-37}^{+38}$	${3837}_{-107}^{+81}$	${1402}_{-2}^{+2}$	${1602}_{-17}^{+18}$	313
NGC 3227 (2014)	Hβ	${1428}_{-106}^{+97}$	${2236}_{-387}^{+487}$	${1301}_{-3}^{+4}$	${1324}_{-17}^{+17}$	226
NGC 3516	Hβ	${2448}_{-74}^{+63}$	${3488}_{-146}^{+219}$	${2633}_{-3}^{+3}$	${3231}_{-15}^{+13}$	339
NGC 4151	Hβ	${1940}_{-22}^{+22}$	${4393}_{-110}^{+110}$	${2078}_{-2}^{+2}$	${5174}_{-32}^{+32}$	369
NGC 5548	Hβ	${2772}_{-34}^{+33}$	${7038}_{-110}^{+133}$	${3056}_{-4}^{+3}$	${1094}_{-9}^{+10}$	329

Note. Columns 3 and 4 give the rms line width and FWHM in the rms spectrum. Columns 5 and 6 give the same but in the mean spectrum. All values are corrected for instrumental broadening and the smoothing introduced by the scaling algorithm (see Section 2.3.1); the FWHM of the Gaussian smoothing kernel is given in Column 7. Apart from Column 7, all values are reported in the rest frame.

Download table as:  ASCIITypeset image

For the mass calculation, we use σ_line from the rms spectra as the line width measure because the rms profile reflects the motions of the gas that is actually responding to the continuum flux variations. For the time delay, we use τ_JAV, though the uncertainties in this quantity depend strongly on the assumption that all errors are Gaussian. In Table 10, we list the time lags and line widths adopted for each data set, and combine these to form the VP. To put the VPs on a calibrated mass scale, we adopt a mean scale factor of f = 4.47 ± 1.25 (Woo et al. 2015). The uncertainty in the scale factor is propagated into the masses given in Table 10.

Table 10.  Black Hole Masses

Object	τJAV (days)	σline(km s−1)	$\mathrm{log}\,\mathrm{VP}({M}_{\odot })$ (current)	$\mathrm{log}\,\mathrm{VP}({M}_{\odot })$ (previous)	$\mathrm{log}M({M}_{\odot })$
(1)	(2)	(3)	(4)	(5)	(6)
Mrk 704	${14.19}_{-0.79}^{+0.87}$	${1860}_{-130}^{+108}$	6.98 ± 0.06	...	7.63 ± 0.14
NGC 3227 (2012)	${2.30}_{-0.20}^{+0.22}$	${1368}_{-37}^{+38}$	5.92 ± 0.05	6.21 ± 0.04	6.57 ± 0.13
NGC 3227 (2014)	2.6 ± 1.0	${1428}_{-106}^{+97}$	6.01 ± 0.19	...	6.66 ± 0.24
NGC 3516	${8.11}_{-0.58}^{+0.75}$	${2448}_{-74}^{+63}$	6.99 ± 0.05	6.86 ± 0.04	7.63 ± 0.13
NGC 4151	${6.59}_{-0.21}^{+0.19}$	${1940}_{-22}^{+22}$	6.68 ± 0.01	6.93 ± 0.04	7.33 ± 0.13
NGC 5548	${3.68}_{-0.52}^{+0.43}$	${2772}_{-34}^{+33}$	6.74 ± 0.06	7.08 ± 0.16	7.39 ± 0.14

Note. Columns 2 and 3 give the time delays measured by JAVELIN and line dispersion in the rms spectra. The log of the virial product (Equation (1)) is in Column 4, and previous determinations of the virial product are in Column 5; the previous NGC 3227 and NGC 3516 measurements are from Denney et al. (2010), NGC 4151 is from Bentz et al. (2006a), and the value for NGC 5548 is the mean and standard deviation of 16 reverberation results drawn from the literature. Column 6 gives the black hole mass based on the data from this campaign and assuming $\langle f\rangle =4.47\pm 1.25$ (Woo et al. 2015).

Download table as:  ASCIITypeset image

The virial product, VP = cτΔV²/G, is useful for comparing the masses derived in different reverberation programs because it involves only the two observables and physical constants. While the VPs obtained here are in reasonable agreement with earlier measurements (Columns 4 and 5 of Table 10), it is also clear that the formal uncertainties derived from the time delay and line width are too small. There is clearly some additional intrinsic scatter in the VP values beyond these formal estimates, indicating additional systematic uncertainties (perhaps due to the choice of integration windows or blended spectral components) and/or underestimated measure uncertainties. The previous value for NGC 5548 in the last row of Table 10 underscores this point: for this entry, we used the mean and standard deviation from 16 previous measurements of the VP based on Hβ reverberation, spanning the range $6.74\lt \mathrm{log}\,\mathrm{VP}\lt 7.38$. The standard deviation of this distribution is ${\rm{\Delta }}\mathrm{log}\,\mathrm{VP}\approx 0.15$, which is probably a good indicator of the true uncertainties in typical measurements. If this is true generally, then the VP values measured here are all in agreement with previous determinations.

We have presented new reverberation results for five bright local Seyfert galaxies. All five have been targets in previous reverberation campaigns. In two cases, Mrk 704 and NGC 4151, previous campaigns did not provide good measurements of the emission-line lags or black hole masses. Mrk 704 did not vary in a fashion conducive to reverberation (Barth et al. 2015), showing only monotonically decreasing light curves. Our new data on NGC 4151 are far more extensive than those from our 2005 campaign (Bentz et al. 2006a), which was abbreviated by poor weather. The other three AGNs—NGC 3227, NGC 3516, and NGC 5548—have been targets in multiple previous reverberation campaigns, and were specifically included in this campaign to compare the velocity-dependent lags, which might be interpreted as indicators of the gross kinematics of the BLR, with previous results from our 2007 campaign (Denney et al. 2009a). In all three cases, the pattern of the lags as a function of velocity has changed. The most likely reason for this is that the BLR structure is probably complex and consists of multiple components—a disk and a wind, for example (e.g., Storchi-Bergmann et al. 2017)—and characterizing any particular velocity bin by a single lag is simply inadequate to describe the BLR structure and kinematics. The important point is that the apparent differences between the 2007 and 2012 results suggest that changes may occur over a BLR dynamical timescale. In a future paper, we will undertake a more detailed analysis of these data with the aim of determining the structure and kinematics of the BLR in these sources and determine whether or not the apparent changes are significant.

G.D.R., C.J.G., B.M.P., and R.W.P. are grateful for the support of the National Science Foundation through grant AST-1008882 to The Ohio State University. K.D.D., B.J.S., C.B.H., and J.L.V. acknowledge support by NSF Fellowships. M.C.B. gratefully acknowledges support from the NSF through CAREER grant AST-1253702. A.M.M. and D.M.S. acknowledge the support of NSF grants AST-1004756 and AST-1009756. C.S.K. is supported by NSF grant AST-1515876. S.K. is supported at the Technion by the Kitzman Fellowship and by a grant from the Israel-Niedersachsen collaboration program. S.R. is supported at Technion by the Zeff Fellowship. S.G.S. acknowledges the support to CrAO in the frame of the "CosmoMicroPhysics" Target Scientific Research Complex Programme of the National Academy of Sciences of Ukraine (2007–2012). M.V. gratefully acknowledges support from the Danish Council for Independent Research via grant no. DFF 4002-00275. V.T.D. acknowledges the support of the Russian Foundation of Research (RF project no. 12-02-01237-a). The CrAO CCD cameras were purchased through the US Civilian Research and Development for Independent States of the Former Soviet Union (CRDF) awards UP1-2116 and UP1-2549-CR-03. This research has been partly supported by the Grants-in-Aid of Scientific Research (17104002, 20041003, 21018003, 21018005, 22253002, and 22540247) of the Ministry of Education, Science, Culture and Sports of Japan. This research has made use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.

Software: IRAF (Tody 1986, 1993), LA Cosmic (van Dokkum 2001), mapspec (Fausnaugh 2017), ISIS (Alard & Lupton 1998; Alard 2000), Sexterp (Siverd et al. 2012), JAVELIN (Zu et al. 2011).