Correlations between the mass of a supermassive black hole and the properties of its host galaxy (e.g., total stellar mass (M*), luminosity (Lhost)) suggest an evolutionary connection. A powerful test of a co-evolution scenario is to measure the relations MBH-Lhost and MBH-M* at high redshift and compare with local estimates. For this purpose, we acquired HST imaging with WFC3 of 32 X-ray-selected broad-line AGN at 1.2<z<1.7 in deep survey fields. By applying state-of-the-art tools to decompose the HST images including available ACS data, we measured the host galaxy luminosity and stellar mass along with other properties through the 2D model fitting. The black hole mass was determined using the broad Halpha line, detected in the near-infrared with Subaru/FMOS, which potentially minimizes systematic effects using other indicators. We find that the observed ratio of MBH to total M* is 2.7 times larger at z~1.5 than in the local universe, while the scatter is equivalent between the two epochs. A non-evolving mass ratio is consistent with the data at the 2-3 sigma confidence level when accounting for selection effects and their uncertainties. The relationship between MBH-Lhost paints a similar picture. Therefore, our results cannot distinguish whether SMBHs and their total M* and Lhost proceed in lockstep or whether the growth of the former somewhat overshoots the latter, given the uncertainties. Based on a statistical estimate of the bulge-to-total mass fraction, the ratio MBH/M* is offset from the local value by a factor of ~7 which is significant even accounting for selection effects. Taken together, these observations are consistent with a scenario in which stellar mass is subsequently transferred from an angular momentum supported component of the galaxy to the pressure supported one through secular processes or minor mergers at a faster rate than mass accretion onto the SMBH.

Active galactic nuclei (AGNs) show a correlation between the size of the broad line region (BLR) and the monochromatic continuum luminosity at 5100 \AA, allowing black hole mass estimation based on single-epoch spectra. However, the validity of the correlation is yet to be clearly tested for high-luminosity AGNs. We present the first reverberation-mapping results of the Seoul National University AGN monitoring program (SAMP), which is designed to focus on luminous AGNs for probing the high end of the size-luminosity relation. We report time lag measurements of two AGNs, namely, 2MASS J10261389+5237510 and SDSS J161911.24+501109.2, using the light curves obtained over a $\sim$1000 day period with an average cadence of $\sim$10 and $\sim$20 days, respectively for photometry and spectroscopy monitoring. Based on a cross-correlation analysis and H$\beta$ line width measurements, we determine the H$\beta$ lag as $41.8^{+4.9}_{-6.0}$ and $52.6^{+17.6}_{-14.7}$ days in the observed-frame, and black hole mass as $3.65^{+0.49}_{-0.57} \times 10^7 M_{\odot}$ and $23.02^{+7.81}_{-6.56} \times 10^7 M_{\odot}$, respectively for 2MASS J1026 and SDSS J1619.

While the reverberation mapping technique is the best available method for measuring black hole mass in active galactic nuclei (AGN) beyond the local volume, this method has been mainly applied to relatively low-to-moderate luminosity AGNs at low redshift. We present the strategy of the Seoul National University AGN Monitoring Project, which aims at measuring the time delay of the \Hb\ line emission with respect to AGN continuum, using a sample of relatively high luminosity AGNs out to z $\sim$0.5. We present simulated cross correlation results based on a number of mock light curves, in order to optimally determine monitoring duration and cadence. We describe our campaign strategy based on the simulation results and the availability of observing facilities. We present the sample selection, and the properties of the selected 100 AGNs, including the optical luminosity, expected time lag, black hole mass, and Eddington ratio.

For a sample of $\sim$80 local ($0.02 \leq z \leq 0.1$) Seyfert-1 galaxies with high-quality long-slit Keck spectra and spatially-resolved stellar-velocity dispersion ($\sigma_{\star}$) measurements, we study the profile of the [OIII]$\lambda$5007A emission line to test the validity of using its width as a surrogate for $\sigma_{\star}$. Such an approach has often been used in the literature, since it is difficult to measure $\sigma_{\star}$ for type-1 active galactic nuclei (AGNs) due to the AGN continuum outshining the stellar-absorption lines. Fitting the [OIII] line with a single Gaussian or Gauss-Hermite polynomials overestimates $\sigma_{\star}$ by 50-100%. When line asymmetries from non-gravitational gas motion are excluded in a double Gaussian fit, the average ratio between the core [OIII] width ($\sigma_{\rm {[OIII],D}}$) and $\sigma_{\star}$ is $\sim$1, but with individual data points off by up to a factor of two. The resulting black-hole-mass-$\sigma_{\rm {[OIII],D}}$ relation scatters around that of quiescent galaxies and reverberation-mapped AGNs. However, a direct comparison between $\sigma_{\star}$ and $\sigma_{\rm {[OIII],D}}$ shows no close correlation, only that both quantities have the same range, average and standard deviation, probably because they feel the same gravitational potential. The large scatter is likely due to the fact that line profiles are a luminosity-weighted average, dependent on the light distribution and underlying kinematic field. Within the range probed by our sample (80-260 km s$^{-1}$), our results strongly caution against the use of [OIII] width as a surrogate for $\sigma_{\star}$ on an individual basis. Even though our sample consists of radio-quiet AGNs, FIRST radio-detected objects have, on average, a $\sim$10% larger [OIII] core width.

We present the single-epoch black hole mass (M$_{\rm BH}$) calibrations based on the rest-frame UV and optical measurements of Mg II 2798\AA\ and H$\beta$ 4861\AA\ lines and AGN continuum, using a sample of 52 moderate-luminosity AGNs at z$\sim$0.4 and z$\sim$0.6 with high-quality Keck spectra. We combine this sample with a large number of luminous AGNs from the Sloan Digital Sky Survey to increase the dynamic range for a better comparison of UV and optical velocity and luminosity measurements. With respect to the reference M$_{\rm BH}$ based on the line dispersion of H$\beta$ and continuum luminosity at 5100\AA, we calibrate the UV and optical mass estimators, by determining the best-fit values of the coefficients in the mass equation. By investigating whether the UV estimators show systematic trend with Eddington ratio, FWHM of H$\beta$, the Fe II strength, and the UV/optical slope, we find no significant bias except for the slope. By fitting the systematic difference of Mg II-based and H$\beta$-based masses with the L$_{3000}$/L$_{5100}$ ratio, we provide a correction term as a function of the spectral index as $\Delta$C = 0.24 (1+$\alpha_{\lambda}$) + 0.17, which can be added to the Mg II-based mass estimators if the spectral slope can be well determined. The derived UV mass estimators typically show $>$$\sim$0.2 dex intrinsic scatter with respect to H$\beta$-based M$_{\rm BH}$, suggesting that the UV-based mass has an additional uncertainty of $\sim$0.2 dex, even if high quality rest-frame UV spectra are available.

We present the single-epoch black hole mass (M$_{\rm BH}$) calibrations based on the rest-frame UV and optical measurements of Mg II 2798\AA\ and H$\beta$ 4861\AA\ lines and AGN continuum, using a sample of 52 moderate-luminosity AGNs at z$\sim$0.4 and z$\sim$0.6 with high-quality Keck spectra. We combine this sample with a large number of luminous AGNs from the Sloan Digital Sky Survey to increase the dynamic range for a better comparison of UV and optical velocity and luminosity measurements. With respect to the reference M$_{\rm BH}$ based on the line dispersion of H$\beta$ and continuum luminosity at 5100\AA, we calibrate the UV and optical mass estimators, by determining the best-fit values of the coefficients in the mass equation. By investigating whether the UV estimators show systematic trend with Eddington ratio, FWHM of H$\beta$, the Fe II strength, and the UV/optical slope, we find no significant bias except for the slope. By fitting the systematic difference of Mg II-based and H$\beta$-based masses with the L$_{3000}$/L$_{5100}$ ratio, we provide a correction term as a function of the spectral index as $\Delta$C = 0.24 (1+$\alpha_{\lambda}$) + 0.17, which can be added to the Mg II-based mass estimators if the spectral slope can be well determined. The derived UV mass estimators typically show $>$$\sim$0.2 dex intrinsic scatter with respect to H$\beta$-based M$_{\rm BH}$, suggesting that the UV-based mass has an additional uncertainty of $\sim$0.2 dex, even if high quality rest-frame UV spectra are available.

We present the variability and time lag measurements of PG 0934+013 based on a photometric and spectroscopic monitoring campaign over a two year period. We obtained 46 epochs of data from the spectroscopic campaign, which was carried out using the Southern African Large Telescope with $\sim$1 week cadence over two sets of 4 month-long observing period, while we obtained 80 epochs of \textit{B}-band imaging data using a few 1-m class telescopes. Due to the seven month gap between the two observing periods, we separately measured the time lags of broad emission lines including H$\beta$, by comparing the emission line light curve with the \textit{B}-band continuum light curve using the cross-correlation function techniques. We determined the H$\beta$ lag, $\tau_{\rm cent} = 8.46^{+2.08}_{-2.14}$ days in the observed-frame based on Year 2 data, while the time lag from Year 1 data was not reliably determined. Using the rms spectrum of Year 2 data, we measured the \Hb\ line dispersion \sigmaline = 668 $\pm$ 44 \kms\ after correcting for the spectral resolution. Adopting a virial factor f = 4.47 from Woo et al. 2015, we determined the black hole mass M$_{BH}$ = $3.13 ^{+0.91} _{-0.93} \times 10^{6}$ \msun based on the \Hb\ time lag and velocity.

We present rest-frame optical spectra of 12 nitrogen-loud quasars at z ~ 2.2, whose rest-frame ultraviolet (UV) spectra show strong nitrogen broad emission lines. To investigate their narrow-line region (NLR) metallicities, we measure the equivalent width (EW) of the [OIII]5007 emission line: if the NLR metallicity is remarkably high as suggested by strong UV nitrogen lines, the [OIII]5007 line flux should be very week due to the low equilibrium temperature of the ionized gas owing to significant metal cooling. In the result, we found that our spectra show moderate EW of the [OIII]5007 line similar to general quasars. This indicates nitrogen-loud quasars do not have extremely metal-rich gas clouds in NLRs. This suggests that strong nitrogen lines from broad-line regions (BLRs) are originated by exceptionally high abundances of nitrogen relative to oxygen without very high BLR metallicities. This result indicates that broad-emission lines of nitrogen are not good indicators of the BLR metallicity in some cases. On the other hand, we also investigate virial black-hole masses and Eddington ratios by using the Hbeta and CIV1549 lines for our sample. As a result, we found that black-hole masses and Eddington ratios of nitrogen-loud quasars tend to be low and high relative to normal quasars, suggesting that nitrogen-loud quasars seem to be in a rapidly-accreting phase. This can be explained in terms of a positive correlation between Eddington ratios and nitrogen abundances of quasars, that is probably caused by the connection between the mass accretion onto black holes and nuclear star formation.

Strongly lensed active galactic nuclei (AGN) provide a unique opportunity to make progress in the study of the evolution of the correlation between the mass of supermassive black holes ($\mathcal M_{BH}$) and their host galaxy luminosity ($L_{host}$). We demonstrate the power of lensing by analyzing two systems for which state-of-the-art lens modelling techniques have been applied to Hubble Space Telescope imaging data. We use i) the reconstructed images to infer the total and bulge luminosity of the host and ii) published broad-line spectroscopy to estimate $\mathcal M_{BH}$ using the so-called virial method. We then enlarge our sample with new calibration of previously published measurements to study the evolution of the correlation out to z~4.5. Consistent with previous work, we find that without taking into account passive luminosity evolution, the data points lie on the local relation. Once passive luminosity evolution is taken into account, we find that BHs in the more distant Universe reside in less luminous galaxies than today. Fitting this offset as $\mathcal M_{BH}$/$L_{host}$ $\propto$ (1+z)$^{\gamma}$, and taking into account selection effects, we obtain $\gamma$ = 0.6 $\pm$ 0.1 and 0.8$\pm$ 0.1 for the case of $\mathcal M_{BH}$-$L_{bulge}$ and $\mathcal M_{BH}$-$L_{total}$, respectively. To test for systematic uncertainties and selection effects we also consider a reduced sample that is homogeneous in data quality. We find consistent results but with considerably larger uncertainty due to the more limited sample size and redshift coverage ($\gamma$ = 0.7 $\pm$ 0.4 and 0.2$\pm$ 0.5 for $\mathcal M_{BH}$-$L_{bulge}$ and $\mathcal M_{BH}$-$L_{total}$, respectively), highlighting the need to gather more high-quality data for high-redshift lensed quasar hosts. Our result is consistent with a scenario where the growth of the black hole predates that of the host galaxy.

We investigate properties of black hole (BH) binaries formed in globular clusters via dynamical processes, using direct N-body simulations. We pay attention to effects of BH mass function on the total mass and mass ratio distributions of BH binaries ejected from clusters. Firstly, we consider BH populations with two different masses in order to learn basic differences from models with single-mass BHs only. Secondly, we consider continuous BH mass functions adapted from recent studies on massive star evolution in a low metallicity environment, where globular clusters are formed. In this work, we consider only binaries that are formed by three-body processes and ignore stellar evolution and primordial binaries for simplicity. Our results imply that most BH binary mergers take place after they get ejected from the cluster. Also, mass ratios of dynamically formed binaries should be close to one or likely to be less than 2:1. Since the binary formation efficiency is larger for higher-mass BHs, it is likely that a BH mass function sampled by gravitational-wave observations would be weighed toward higher masses than the mass function of single BHs for a dynamically formed population. Applying conservative assumptions regarding globular cluster populations such as small BH mass fraction and no primordial binaries, the merger rate of BH binaries originated from globular clusters is estimated to be at least 6.5 per yr per Gpc^3. Actual rate can be up to more than several times of our conservative estimate.

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.

A sample of 102 local (0.02 < z < 0.1) Seyfert galaxies with black hole masses MBH > 10^7 M_sun was selected from the Sloan Digital Sky Survey (SDSS) and observed using the Keck 10-m telescope to study the scaling relations between MBH and host galaxy properties. We study profile changes of the broad Hbeta emission line within the ~3-9 year time-frame between the two sets of spectra. The variability of the broad Hbeta emission line is of particular interest, not only since it is used to estimate MBH, but also since its strength and width is used to classify Seyfert galaxies into different types. At least some form of broad-line variability (in either width or flux) is observed in the majority (~66%) of the objects, resulting in a Seyfert-type change for ~38% of the objects, likely driven by variable accretion and/or obscuration. The broad Hbeta line virtually disappears in 3/102 (~3%) extreme cases. We discuss potential causes for these changing-look AGNs. While similar dramatic transitions have previously been reported in the literature, either on a case-by-case basis or in larger samples focusing on quasars at higher redshifts, our study provides statistical information on the frequency of H$\beta$ line variability in a sample of low-redshift Seyfert galaxies.

The ability to accurately derive black hole (BH) masses at progressively higher redshifts and over a wide range of continuum luminosities has become indispensable in the era of large-area extragalactic spectroscopic surveys. In this paper we present an extension of existing comparisons between rest-frame UV and optical virial BH mass estimators to intermediate redshifts and luminosities comparable to the local H$\beta$ reverberation mapped active galactic nuclei (AGN). We focus on the MgII, CIV, and CIII] broad emission lines and compare them to both H$\alpha$ and H$\beta$. We use newly acquired near-infrared spectra from the FMOS instrument on the Subaru telescope for 89 broad-lined AGN at redshifts between 0.3 and 3.5, complemented by data from the AGES survey. We employ two different prescriptions for measuring the emission line widths and compare the results. We confirm that MgII shows a tight correlation with H$\alpha$ and H$\beta$, with a scatter of ~0.25 dex. The CIV and CIII] estimators, while showing larger scatter, are viable virial mass estimators after accounting for a trend with the UV-to-optical luminosity ratio. We find an intrinsic scatter of ~0.37 dex between Balmer and carbon virial estimators by combining our dataset with previous high redshift measurements. This updated comparison spans a total of 3 decades in BH mass. We calculate a virial factor for CIV/CIII] logf(CIV/CIII])=0.87 with an estimated systematic uncertainty of ~0.4 dex and find excellent agreement between the local reverberation mapped AGN sample and our high-z sample.

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.

Intra Cluster Media (ICMs) located at galaxy clusters is in the state of hot, tenuous, magnetized, and highly ionized X-ray emitting plasmas. This overall collisionless, viscous, and conductive magnetohydrodynamic (MHD) turbulence in ICM is simulated using hyper and physical magnetic diffusivity. The results show that fluctuating random plasma motion amplifies the magnetic field, which cascades toward the diffusivity scale passing through the viscous scale. The kinetic eddies in the subviscous scale are driven and constrained by the magnetic tension which finally gets balanced with the highly damping effect of the kinetic eddies. However, the saturated kinetic energy spectrum is deeper than that of the incompressible or compressible hydrodynamics fluid. To explain this unusual field profile we set up two simultaneous differential equations for the kinetic and magnetic energy spectrum using an Eddy Damped Quasi Normal Markovianized (EDQNM) approximation. The analytic solution tells us that the magnetic energy in addition to the viscous damping effect constrains the plasma motion leading to the power spectra: kinetic energy spectrum $E_V^k\sim k^{-3}$ and corresponding representative magnetic energy spectrum $E_M^k\sim k^{-1/2}$. Also the comparison of simulation results with different resolutions and magnetic diffusivities implies the role of small scale magnetic energy in dynamo.

We create a baseline of the black hole (BH) mass (MBH) - stellar-velocity dispersion (sigma) relation for active galaxies, using a sample of 66 local (0.02<z<0.09) Seyfert-1 galaxies, selected from the Sloan Digital Sky Survey (SDSS). Analysis of SDSS images yields AGN luminosities free of host-galaxy contamination and morphological classification. 51/66 galaxies have spiral morphology. 28 bulges have Sersic index n<2 and are considered candidate pseudo bulges, with eight being definite pseudo bulges based on multiple classification criteria met. Only 4/66 galaxies show sign of interaction/merging. High signal-to-noise ratio Keck spectra provide the width of the broad Hbeta emission line free of FeII emission and stellar absorption. AGN luminosity and Hbeta line widths are used to estimate MBH. The Keck-based spatially-resolved kinematics is used to determine stellar-velocity dispersion within the spheroid effective radius. We find that sigma can vary on average by up to 40% across definitions commonly used in the literature, emphasizing the importance of using self-consistent definitions in comparisons and evolutionary studies. The MBH-sigma relation for our Seyfert-1 galaxies has the same intercept and scatter as that of reverberation-mapped AGNs as well as quiescent galaxies, consistent with the hypothesis that our single epoch MBH estimator and sample selection do not introduce significant biases. Barred galaxies, merging galaxies, and those hosting pseudo bulges do not represent outliers in the MBH-sigma relation. This is in contrast with previous work, although no firm conclusion can be drawn due to the small sample size and limited resolution of the SDSS images.

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.

Narrow-line Seyfert 1 galaxies (NLS1s) are arguably one of the key AGN subclasses in investigating the origin of the black hole mass - stellar velocity dispersion (M-sigma) relation because of their high accretion rate and significantly low black hole mass. Currently, it is under discussion whether present-day NLS1s offset from the M-sigma relation. Using the directly measured stellar velocity dispersion of 93 NLS1s at z<0.1, and black hole mass estimates based on the updated mass estimators, we investigate the M-sigma relation of NLS1s in comparison with broad-line AGNs. We find no strong evidence that the NLS1s deviates from the M-sigma relation, which is defined by reverberation-mapped type 1 AGNs and quiescent galaxies. However, there is a clear trend of the offset with the host galaxy morphology, i.e., more inclined galaxies toward the line-of-sight have higher stellar velocity dispersion, suggesting that the rotational broadening plays a role in measuring stellar velocity dispersion based on the single-aperture spectra from the Sloan Digital Sky Survey. In addition, we provide the virial factor log f = $0.05 \pm 0.12$ (f = 1.12), for black hole mass estimators based on the FWHM of H$\beta$, by jointly fitting the M-sigma relation using quiescent galaxies and reverberation-mapped AGNs.

We use the high-resolution simulation MassiveBlackII to examine scaling relations between black hole mass (MBH) and host galaxy properties (sigma, M*, and LV), finding good agreement with observational data, especially at the high-mass end. The simulations have less intrinsic scatter than observations, and the MBH-LV correlation has the largest scatter, suggesting it may the the least fundamental of the three relations. We find Gaussian scatter about all three relations, except among the highest mass galaxies, which host more massive black holes. Below z~2 the slopes for the full population remain roughly z-independent, and only steepen by 50% by z~4. The normalization of the sigma, LV relations evolve by 0.3, 0.43 dex, while the MBH correlation does not evolve to at least z~2. Testing for selection biases, we find samples selected by MBH or M* have steeper slopes than randomly selected samples. If unaccounted for, such a selection function would find faster evolution than inferred from a randomly selected sample, as objects at the highend of the relation tend to evolve more rapidly. We find a potential bias among high-LBH subsamples (tending to reside in higher mass galaxies), but these bright-AGN exhibit no intrinsic bias relative to fainter ones in equivalent-mass hosts, nor is there a significant difference between active- and inactive-samples. Finally we characterize the evolution of individual black holes along the scaling planes. Below the local relation, black holes grow faster than their host (72% of black holes 0.3 dex below the mean relation have a MBH-M* trajectory steeper than the local relation), while those above have shallower trajectories (only 14% are steeper than local). Black holes tend to grow faster than their hosts until surpassing the local relation, at which point their growth is suppressed while their hosts continue to grow, returning them to the mean relation.

We investigate the cosmic evolution of the black hole (BH) mass -- bulge luminosity relation using a sample of 52 active galaxies at $z \sim 0.36$ and $z \sim 0.57$ in the BH mass range of $10^{7.4-9.1} M_{\odot}$. By consistently applying multi-component spectral and structural decomposition to high-quality Keck spectra and high-resolution HST images, BH masses ($M_{\rm BH}$) are estimated using the H$\beta$ broad emission line combined with the 5100 \AA\ nuclear luminosity, and bulge luminosities ($L_{\rm bul}$) are derived from surface photometry. Comparing the resulting $M_{\rm BH}-L_{\rm bul}$ relation to local active galaxies and taking into account selection effects, we find evolution of the form $M_{\rm BH} / L_{\rm bul} \propto (1+z)^{\gamma}$ with $\gamma=1.8\pm0.7$, consistent with BH growth preceding that of the host galaxies. Including an additional sample of 27 active galaxies with $0.5<z<1.9$ taken from the literature and measured in a consistent way, we obtain $\gamma=0.9\pm0.7$ for the $M_{\rm BH}-L_{\rm bul}$ relation and $\gamma=0.4\pm0.5$ for the $M_{\rm BH}$--total host galaxy luminosity ($L_{\rm host}$) relation. The results strengthen the findings from our previous studies and provide additional evidence for host-galaxy bulge growth being dominated by disk-to-bulge transformation via minor mergers and/or disk instabilities.

We search for misclassified type 1 AGNs among type 2 AGNs identified with emission line flux ratios, and investigate the properties of the sample. Using 4\,113 local type 2 AGNs at $0.02<z<0.05$ selected from Sloan Digital Sky Survey Data Release 7, we detected a broad component of the \Ha\ line with a Full-Width at Half-Maximum (FWHM) ranging from 1\,700 to 19\,090 \kms\ for 142 objects, based on the spectral decomposition and visual inspection. The fraction of the misclassified type 1 AGNs among type 2 AGN sample is $\sim$3.5%, implying that a large number of missing type 1 AGN population may exist. The misclassified type 1 AGNs have relatively low luminosity with a mean broad \Ha\ luminosity, log L$_{H\alpha} = 40.50\pm0.35$ \ergs, while black hole mass of the sample is comparable to that of the local black hole population, with a mean black hole mass, log M$_{\rm BH} = 6.94\pm0.51$ M$_{\odot}$. The mean Eddington ratio of the sample is log L$_{\rm bol}$/L$_{\rm Edd}$ = $-2.00\pm0.40$, indicating that black hole activity is relatively weak, hence, AGN continuum is too weak to change the host galaxy color. We find that the \OIII\ lines show significant velocity offsets, presumably due to outflows in the narrow-line region, while the velocity offset of the narrow component of the \Ha\ line is not prominent, consistent with the ionized gas kinematics of general type 1 AGN population.

We present the kinematic properties of a type-2 QSO, SDSS J132323.33-015941.9 at z~0.35, based on the analysis of Very Large Telescope integral field spectroscopy and Hubble Space Telescope (HST) imaging, which suggest that the target is a binary active galactic nucleus (AGN) with double narrow-line regions. The QSO features double-peaked emission lines ([OIII] and Hb) which can be decomposed into two kinematic components. The flux-weighted centroids of the blue and red components are separated by ~0.2" (0.8 kpc in projection) and coincide with the location of the two stellar cores detected in the HST broad-band images, implying that both stellar cores host an active black hole. The line-of-sight velocity of the blue component is comparable to the luminosity-weighted velocity of stars in the host galaxy while the red component is redshifted by ~240 km/s, consistent with typical velocity offsets of two cores in a late stage of a galaxy merger. If confirmed, the target is one of the rare cases of sub-kpc scale binary AGNs, providing a test-bed for understanding the binary AGN population.

We present the single-epoch black hole mass estimators based on the CIV (1549 A) broad emission line, using the updated sample of the reverberation-mapped AGNs and high-quality UV spectra. By performing multi-component spectral fitting analysis, we measure the CIV line widths (FWHM_CIV) and line dispersion (sigma_CIV) and the continuum luminosity at 1350 A (L_1350) to calibrate the CIV-based mass estimators. By comparing with the Hbeta reverberation-based masses, we provide new mass estimators with the best-fit relationships, i.e., M_BH \propto L_1350 ^ (0.50+-0.07) sigma_CIV ^2 and M_BH \propto L_1350 ^ (0.52+-0.09) FWHM_CIV ^ (0.56+-0.48). The new CIV-based mass estimators show significant mass-dependent systematic difference compared to the estimators commonly used in the literature. Using the published Sloan Digital Sky Survey QSO catalog, we show that the black hole mass of high-redshift QSOs decreases on average by ~0.25 dex if our recipe is adopted.

To investigate the validity of the assumption that quiescent galaxies and active galaxies follow the same black hole mass - stellar velocity dispersion relation, as required for the calibration of M_BH estimators for broad line AGNs, we determine and compare the M-sigma relations, respectively, for quiescent and active galaxies. For the quiescent galaxy sample, composed of 72 dynamical M_BH measurements, we update sigma for 28 galaxies using homogeneous H-band measurements that are corrected for galaxy rotation. For active galaxies, we collect 25 reverberation-mapped AGNs and improve sigma measurement for two objects . Combining the two samples, we determine the virial factor $f$, first by scaling the active galaxy sample to the M-sigma relation of quiescent galaxies, and second by simultaneously fitting the quiescent and active galaxy samples, as f=5.1_{-1.1}^{+1.5} and f=5.9_{-1.5}^{+2.1}, respectively. The M-sigma relation of active galaxies appears to be shallower than that of quiescent galaxies. However, the discrepancy is caused by a difference in the accessible M_BH distribution at given sigma, primarily due to the difficulty of measuring reliable stellar velocity dispersion for the host galaxies of luminous AGNs. Accounting for the selection effects, we find that active and quiescent galaxies are consistent with following intrinsically the same M-sigma relation.

To calibrate stellar velocity dispersion measurements from optical and near-IR stellar lines, and to improve the black hole mass (MBH)- stellar velocity dispersion (sigma) relation, we measure stellar velocity dispersions based on high quality H-band spectra for a sample of 31 nearby galaxies, for which dynamical MBH is available in the literature. By comparing velocity dispersions measured from stellar lines in the H-band with those measured from optical stellar lines, we find no significant difference, suggesting that optical and near-IR stellar lines represent the same kinematics and that dust effect is negligible for early-type galaxies. Based on the spatially-resolved rotation and velocity dispersion measurements along the major axis of each galaxy, we find that a rotating stellar disk is present for 80% of galaxies in the sample. For galaxies with a rotation component, velocity dispersions measured from a single aperture spectrum can vary by up to ~20%, depending on the size of the adopted extraction aperture. To correct for the rotational broadening, we derive luminosity-weighted velocity dispersion within the effective radius of each galaxy, providing uniformly measured velocity dispersions to improve the MBH-sigma relation.

Determining the virial factor of the broad-line region (BLR) gas is crucial for calibrating AGN black hole mass estimators, since the measured line-of-sight velocity needs to be converted into the intrinsic virial velocity. The average virial factor has been empirically calibrated based on the M-sigma relation of quiescent galaxies, but the claimed values differ by a factor of two in recent studies. We investigate the origin of the difference by measuring the M-sigma relation using an updated galaxy sample from the literature, and explore the dependence of the virial factor on various fitting methods. We find that the discrepancy is primarily caused by the sample selection, while the difference stemming from the various regression methods is marginal. However, we generally prefer the FITEXY and Bayesian estimators based on Monte Carlo simulations for the M-sigma relation. In addition, the choice of independent variable in the regression leads to ~0.2 dex variation in the virial factor inferred from the calibration process. Based on the determined virial factor, we present the updated M-sigma relation of local active galaxies.

Young radio galaxies (YRGs) provide an ideal laboratory to explore the connection between accretion disk and radio jet thanks to their recent jet formation. We investigate the relationship between the emission-line properties, the black hole accretion rate, and the radio properties using a sample of 34 low-redshift (z < 0.4) YRGs. We classify YRGs as high-excitation galaxies (HEGs) and low-excitation galaxies (LEGs) based on the flux ratio of high-ionization to low-ionization emission lines. Using the H{\alpha} luminosities as a proxy of accretion rate, we find that HEGs in YRGs have \sim1 dex higher Eddington ratios than LEGs in YRGs, suggesting that HEGs have higher mass accretion rate or higher radiative efficiency than LEGs. In agreement with previous studies, we find that the luminosities of emission lines, in particular H{\alpha}, are correlated with radio core luminosity, suggesting that accretion and young radio activities are fundamentally connected.

PreCam, a precursor observational campaign supporting the Dark Energy Survey (DES), is designed to produce a photometric and astrometric catalog of nearly a hundred thousand standard stars within the DES footprint, while the PreCam instrument also serves as a prototype testbed for the Dark Energy Camera (DECam)'s hardware and software. This catalog represents a potential 100-fold increase in Southern Hemisphere photometric standard stars, and therefore will be an important component in the calibration of the Dark Energy Survey. We provide details on the PreCam instrument's design, construction and testing, as well as results from a subset of the 51 nights of PreCam survey observations on the University of Michigan Department of Astronomy's Curtis-Schmidt telescope at Cerro Tololo Inter-American Observatory. We briefly describe the preliminary data processing pipeline that has been developed for PreCam data and the preliminary results of the instrument performance, as well as astrometry and photometry of a sample of stars previously included in other southern sky surveys.

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.

We investigate the connection between starburst and AGN activity by comparing the 3.3 {\mu}m PAH eimission with AGN properties. Utilizing the slit-less spectroscopic capability of the AKARI space telescope, we observe moderate-luminosity Type I AGN at z~0.4 to measure global starburst activity. The 3.3 {\mu}m PAH emissions are detected for 7 out of 26 target galaxies. We find no strong correlation between the 3.3 {\mu}m PAH emission and AGN luminosity in the limted range of the observed AGN luminosity, suggesting that global star formation may not be tightly related with AGN activity. Combining our measurements with the previous 3.3 {\mu}m measurements of the low redshift Type I AGN in the literature, we investigate the connection between nuclear starburst and AGN activity. In contrast to global star formation, the 3.3 {\mu}m PAH luminosity measured from the central part of galaxies correlates with AGN luminosity, implying that starburst and AGN activity are directly connected at the nuclear region.

We investigate the calibration and uncertainties of black hole mass estimates based on the single-epoch (SE) method, using homogeneous and high-quality multi-epoch spectra obtained by the Lick Active Galactic Nucleus (AGN) Monitoring Project for 9 local Seyfert 1 galaxies with black hole masses < 10^8 M_sun. By decomposing the spectra into their AGN and stellar components, we study the variability of the single-epoch Hbeta line width (full width at half-maximum intensity, FWHM_Hbeta; or dispersion, sigma_Hbeta) and of the AGN continuum luminosity at 5100A (L_5100). From the distribution of the "virial products" (~ FWHM_Hbeta^2 L_5100^0.5 or sigma_Hbeta^2 L_5100^0.5) measured from SE spectra, we estimate the uncertainty due to the combined variability as ~ 0.05 dex (12%). This is subdominant with respect to the total uncertainty in SE mass estimates, which is dominated by uncertainties in the size-luminosity relation and virial coefficient, and is estimated to be ~ 0.46 dex (factor of ~ 3). By comparing the Hbeta line profile of the SE, mean, and root-mean-square (rms) spectra, we find that the Hbeta line is broader in the mean (and SE) spectra than in the rms spectra by ~ 0.1 dex (25%) for our sample with FWHM_Hbeta < 3000 km/s. This result is at variance with larger mass black holes where the difference is typically found to be much less than 0.1 dex. To correct for this systematic difference of the Hbeta line profile, we introduce a line-width dependent virial factor, resulting in a recalibration of SE black hole mass estimators for low-mass AGNs.

The Lick AGN Monitoring Project 2011 observing campaign was carried out over the course of 11 weeks in Spring 2011. Here we present the first results from this program, a measurement of the broad-line reverberation lag in the Seyfert 1 galaxy Mrk 50. Combining our data with supplemental observations obtained prior to the start of the main observing campaign, our dataset covers a total duration of 4.5 months. During this time, Mrk 50 was highly variable, exhibiting a maximum variability amplitude of a factor of 4 in the U-band continuum and a factor of 2 in the H-beta line. Using standard cross-correlation techniques, we find that H-beta and H-gamma lag the V-band continuum by tau_cen = 10.64(-0.93,+0.82) and 8.43(-1.28,+1.30) days, respectively, while the lag of He II 4686 is unresolved. The H-beta line exhibits a symmetric velocity-resolved reverberation signature with shorter lags in the high-velocity wings than in the line core, consistent with an origin in a broad-line region dominated by orbital motion rather than infall or outflow. Assuming a virial normalization factor of f=5.25, the virial estimate of the black hole mass is (3.2+-0.5)*10^7 solar masses. These observations demonstrate that Mrk 50 is among the most promising nearby active galaxies for detailed investigations of broad-line region structure and dynamics.

Cosmic filaments play a role of bridges along which matter and gas accrete onto galaxies to trigger star formation and feed central black holes. Here we explore the correlations between the intrinsic properties of void galaxies and the linearity R_L of void filaments (degree of filament's straightness). We focus on void regions since the bridge effect of filaments should be most conspicuous in the pristine underdense regions like voids. Analyzing the Millennium-Run semi-analytic galaxy catalogue, we identify void filaments consisting of more than four galaxies (three edges) and calculate the means of central black hole mass, star formation rate, and stellar mass as a function of R_L. It is shown that the void galaxies constituting more straight filaments tend to have higher luminosity, more massive central black holes and higher star formation rate. Among the three properties, the central black hole mass is most strongly correlated with R_L. It is also shown that the dark halos constituting straight filaments tend to have similar masses. Our results suggest that the fuel-supply for central black holes and star formation of void galaxies occurs most efficiently along straight void filaments whose potential wells are generated by similar-mass dark halos.

The size distribution of mini-filaments in voids has been derived from the Millennium Run halo catalogs at redshifts z=0,0.5,1 and 2. It is assumed that the primordial tidal field originated the presence of filamentary substructures in voids and that the void filaments have evolved only little, keeping the initial memory of the primordial tidal field. Applying the filament-finding algorithm based on the minimal spanning tree (MST) technique to the Millennium voids, we identify the mini-filaments running through voids and measure their sizes at each redshift. Then, we calculate the comoving number density of void filaments as a function of their sizes in the logarithmic interval and determine an analytic fitting function for it. It is found that the size distribution of void mini-filaments in the logarithmic interval has an almost universal shape, insensitive to the redshift: In the short-size section it is well approximated as a power-law, while in the long-size section it decreases exponentially. We expect that the universal size distribution of void filaments may provide a useful cosmological probe without resorting to the rms density fluctuations.

Our universe is observed to be accelerating due to the dominant dark energy with negative pressure. The dark energy equation of state (w) holds a key to understanding the ultimate fate of the universe. The cosmic voids behave like bubbles in the universe so that their shapes must be quite sensitive to the background cosmology. Assuming a flat universe and using the priors on the matter density parameter (Omega_m) and the dimensionless Hubble parameter (h), we demonstrate analytically that the ellipticity evolution of cosmic voids may be a sensitive probe of the dark energy equation of state. We also discuss the parameter degeneracy between w and Omega_m.

We investigate alignments between the spin axes of cosmic voids and the principal axes of nearby superclusters using the Millennium Run simulation of a LCDM cosmology. The concept of void spin was first introduced by Lee and Park in 2006 to quantify the tidal effect on voids from the surrounding matter distribution. Our numerical analysis reveals that the void spin axes are strongly aligned with the supercluster minor axes, but anti-aligned with the major axes, and have no correlations with the intermediate axes. We provide physical explanations to this numerical results on the basis of tidally induced correlations. It is expected that our work will provide a new insight into the characterization of the cosmic web on the largest scale.

Cosmic voids refer to the large empty regions in the universe with a very low number density of galaxies. Voids are likely to be severely disturbed by the tidal effect from the surrounding dark matter. We derive a completely analytic model for the void ellipticity distribution from physical principles. We use the spatial distribution of galaxies in a void as a measure of its shape, tracking the trajectory of the void galaxies under the influence of the tidal field using the Lagrangian perturbation theory. Our model implies that the void ellipticity distribution depends sensitively on the cosmological parameters. Testing our model against the high-resolution Millennium Run simulation, we find excellent quantitative agreement of the analytic predictions with the numerical results.

We present a theoretical study of void spins and their correlation properties. The concept of the spin angular momentum for an unbound void is introduced to quantify the effect of the tidal field on the distribution of matter that make up the void. Both the analytical and numerical approaches are used for our study. Analytically, we adopt the linear tidal torque model to evaluate the void spin-spin and spin-density correlations, assuming that a void forms in the initial region where the inertia momentum and the tidal shear tensors are maximally uncorrelated with each other. Numerically, we use the Millennium run galaxy catalog to find voids and calculate their spin statistics. The numerical results turn out to be in excellent agreement with the analytic predictions, both of which consistently show that there are strong spatial alignments between the spin axes of neighbor voids and strong anti-alignments between the void spin axes and the directions to the nearest voids. We expect that our work will provide a deeper insight into the origin and properties of voids and the large scale structure.