We present the variability study of the lowest-luminosity Seyfert 1 galaxy NGC 4395 based on the photometric monitoring campaigns in 2017 and 2018. Using 22 ground-based and space telescopes, we monitored NGC 4395 with a $\sim$5 minute cadence during a period of 10 days and obtained light curves in the UV, V, J, H, and K/Ks bands as well as the H$\alpha$ narrow-band. The RMS variability is $\sim$0.13 mag on \emph{Swift}-UVM2 and V filter light curves, decreasing down to $\sim$0.01 mag on K filter. After correcting for continuum contribution to the H$\alpha$ narrow-band, we measured the time lag of the H$\alpha$ emission line with respect to the V-band continuum as ${55}^{+27}_{-31}$ to ${122}^{+33}_{-67}$ min. in 2017 and ${49}^{+15}_{-14}$ to ${83}^{+13}_{-14}$ min. in 2018, depending on the assumption on the continuum variability amplitude in the H$\alpha$ narrow-band. We obtained no reliable measurements for the continuum-to-continuum lag between UV and V bands and among near-IR bands, due to the large flux uncertainty of UV observations and the limited time baseline. We determined the AGN monochromatic luminosity at 5100\AA\ $\lambda L_\lambda = \left(5.75\pm0.40\right)\times 10^{39}\,\mathrm{erg\,s^{-1}}$, after subtracting the contribution of the nuclear star cluster. While the optical luminosity of NGC 4395 is two orders of magnitude lower than that of other reverberation-mapped AGNs, NGC 4395 follows the size-luminosity relation, albeit with an offset of 0.48 dex ($\geq$2.5$\sigma$) from the previous best-fit relation of Bentz et al. (2013).

We have developed observation control software for the Immersion GRating INfrared Spectrometer (IGRINS) slit-viewing camera module, which maintains the position of an astronomical target on the spectroscopic slit. It is composed of several packages that monitor and control the system, acquire the images, and compensate for the tracking error by sending tracking feedback information to the telescope control system. For efficient development and maintenance of each software package, we have applied software engineering methods, i.e., a spiral software development with model-based design. It is not trivial to define the shape and center of astronomical object point spread functions (PSFs), which do not have symmetric Gaussian profiles in short exposure (<4 s) guiding images. Efforts to determine the PSF centroid are additionally complicated by the core saturation of bright guide stars. We have applied both a two-dimensional Gaussian fitting algorithm (2DGA) and center balancing algorithm (CBA) to identify an appropriate method for IGRINS in the near-infrared K-band. The CBA derives the expected center position along the slit-width by referencing the spillover flux ratio of the PSF wings on both sides of the slit. In this research, we have compared the accuracy and reliability of the CBA to the 2DGA by using data from IGRINS commissioning observations at McDonald Observatory. We find that the performance of each algorithm depends on the brightness of the targets and the seeing conditions, with the CBA performing better in typical observing scenarios. The algorithms and test results we present can be utilized with future spectroscopic slit observations in various observing conditions and for a variety of spectrograph designs.

We present the kinematical properties of the UV and optical \FeII\ emission gas based on the velocity shift and line width measurements of a sample of 223 Type 1 active galactic nuclei (AGNs) at 0.4 $<$ z $<$ 0.8. We find a strong correlation between the line widths of the UV and optical \FeII\ emission lines, indicating that both \FeII\ emission features arise from similar distances in the broad line region (BLR). However, in detail we find differing trends, depending on the width of \FeII. While the velocity shifts and dispersions of the UV Fe II (\FeUV) and optical Fe II (\FeOPT) emission lines are comparable to each other for AGNs with relatively narrow \FeOPT\ line widths (i.e., FWHM < 3200 \kms; Group A), \FeOPT\ is broader than \FeUV\ for AGNs with relatively broad \FeOPT\ (i.e., FWHM > 3200 \kms; Group B). \FeII\ emission lines are on average narrower than \Hb\ and \MgII\ for Group A, indicating the \FeII\ emission region is further out in the BLR, while for Group B AGNs \FeOPT\ is comparable to \Hb\ and broader than \MgII. While \FeII\ emission lines are on average redshifted ($40\pm141$ \kms\ and $182\pm95$, respectively for \FeUV\ and \FeOPT), indicating inflow, the sample as a whole shows a large range of velocity shifts, suggesting complex nature of gas kinematics.

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.

The spectral energy distribution of blazars around the synchrotron peak can be well described by the log-parabolic model that has three parameters: peak energy ($E_\textrm{p}$), peak luminosity ($L_\textrm{p}$) and the curvature parameter ($b$). It has been suggested that $E_\textrm{p}$ shows relations with $L_\textrm{p}$ and $b$ in several sources, which can be used to constrain the physical properties of the emitting region and/or acceleration processes of the emitting particles. We systematically study the $E_\textrm{p}$-$L_\textrm{p}$ and $E_\textrm{p}$-(1$/b$) relations for 14 BL Lac objects using the 3-25~keV $RXTE$/PCA and 0.3-10~keV $Swift$/XRT data. Most objects (9/14) exhibit positive $E_\textrm{p}$-$L_\textrm{p}$ correlations, three sources show no correlation, and two sources display negative correlations. In addition, most targets (7/14) present no correlation between $E_\textrm{p}$ and 1$/b$, five sources pose negative correlations, and two sources demonstrate positive correlations. 1ES~1959+650 displays two different $E_\textrm{p}$-$L_\textrm{p}$ relations in 2002 and 2016. We also analyze $E_\textrm{p}$-$L_\textrm{p}$ and $E_\textrm{p}$-(1$/b$) relations during flares lasting for several days. The $E_\textrm{p}$-$L_\textrm{p}$ relation does not exhibit significant differences between flares, while the $E_\textrm{p}$-(1$/b$) relation varies from flare to flare. For the total sample, when $L_\textrm{p}$ < $\textrm{10}^\textrm{45}\ \textrm{erg}\ \textrm{s}^\textrm{-1}$, there seems to be a positive $E_\textrm{p}$-$L_\textrm{p}$ correlation. $L_\textrm{p}$ and the slope of $E_\textrm{p}$-$L_\textrm{p}$ relation present an anti-correlation, which indicates that the causes of spectral variations might be different between luminous and faint sources. $E_\textrm{p}$ shows a positive correlation with the black hole mass. We discuss the implications of these results.

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.

The motions of gas and stars in the nuclei of nearby large galaxies have demonstrated that massive black holes are common and that their masses strongly correlate with the stellar velocity dispersion $\sigma_{\star}$ of the bulge. This correlation suggests that massive black holes and galaxies influence each other's growth. Dynamical measurements are less reliable when the sphere of influence is unresolved, thus it remains unknown whether this correlation exists in galaxies much smaller than the Milky Way, as well as what fraction of these galaxies have central black holes. Light echoes from photoionized clouds around accreting black holes, in combination with the velocity of these clouds, yield a direct mass measurement that circumvents this difficulty. Here we report an exceptionally low reverberation delay of $83\pm14$ minutes between variability in the accretion disk and high velocity H$\alpha$ emission from the nucleus of the bulgeless dwarf galaxy NGC~4395. Combined with the H$\alpha$ line-of-sight velocity dispersion $\sigma_{\rm line}=426\pm1$~km~s$^{-1}$, this lag determines a mass of about 10,000~$M_{\odot}$ for the black hole. This mass is among the smallest central black hole masses reported, near the low end of expected masses for heavy "seeds", and the best direct mass measurement for a galaxy of this size. Despite the lack of a bulge, NGC~4395 is consistent with the $M_{\rm BH} - \sigma_{\star}$ relation when $\sigma_{\star}$ is measured from the central region. This indicates that the relation need not originate from hierarchical galaxy assembly nor from black hole feedback.

The FeII/MgII line flux ratio has been used to investigate the chemical evolution of high redshift active galactic nuclei (AGNs). No strong evolution has been found out to z~6, implying that the type 1a supernova activity has been already occurred in early universe. However, the trend of no evolution can be caused by the sample selection bias since the previous studies have utilized mostly very luminous AGNs, which may be already chemically matured at the observed redshift. As motivated by the previously reported correlation between AGN luminosity and metallicity, we investigate the FeII/MgII flux ratio over a large dynamic range of luminosity, by adding a new sample of 12 quasars at z~3, of which the lower luminosity limit is more than 1 dex smaller than that of the previously studied high-z quasars. Based on the Gemini/GNIRS observations, we find that the seven low-luminosity quasars with a mean bolometric luminosity log L_bol~46.5 \pm 0.2 has an average FeII/MgII ratio of 0.68 \pm 0.11 dex. This ratio is comparable to that of high-luminosity quasars (log L_bol~47.3 \pm 0.3) in our sample (i.e., FeII/MgII ratio of 0.59 \pm 0.15 dex) and that of the previously studied high-luminosity quasars at higher redshifts. One possible scenario is that the low-luminosity quasars in our sample are still relatively luminous and already chemically matured. To search for chemically-young AGNs, and fully understand the chemical evolution based on the FeII/MgII flux ratio, it is required to investigate much lower-luminosity AGNs.

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.

To investigate the connection between radio activity and AGN outflows, we present a study of ionized gas kinematics based on [O III] $\lambda$5007 emission line along the large-scale radio jet for six radio AGNs. These AGNs are selected based on the radio activity (i.e., $\mathrm{L_{1.4GHz}}$ $\geqslant$ 10$^{39.8}$ erg s$^{-1}$) as well as optical emission line properties as type 2 AGNs. Using the Red Channel Cross Dispersed Echellette Spectrograph at the Multiple Mirror Telescope, we investigate in detail the [O III] and stellar kinematics. We spatially resolve and probe the central AGN-photoionization sizes, which is important in understanding the structures and evolutions of galaxies. We find that the typical central AGN-photoionization radius of our targets are in range of 0.9$-$1.6 kpc, consistent with the size-luminosity relation of [O III] in the previous studies. We investigate the [O III] kinematics along the large-scale radio jets to test whether there is a link between gas outflows in the narrow-line region and extended radio jet emissions. Contrary to our expectation, we find no evidence that the gas outflows are directly connected to the large scale radio jets.

We present analysis of near-infrared, high-resolution spectroscopy towards the Flat-spectrum YSO (Young Stellar Object) ESO H$\alpha$ 279a (1.5 solar mass) in the Serpens star forming region, at the distance of 429 pc. Using the Immersion GRating INfrared Spectrometer (IGRINS, R=45,000), we detect emission lines originating from the accretion channel flow, jet, and inner disk. Specifically, we identify hydrogen Brackett series recombination, [Fe II], [Fe III], [Fe IV], Ca I, Na I, H2, H2O and CO overtone emission lines. By modeling five bands of CO overtone emission lines, and the symmetric double-peaked line profile for Na I emission lines, we find that ESO H$\alpha$ 279a has an actively accreting Keplerian disk. From our Keplerian disk model, we find that Na I emission lines originate between 0.04 AU and 1.00 AU, while CO overtone emission lines are from the outer part of disk, in the range between 0.22 AU and 3.00 AU. It reveals that the neutral atomic Na gas is a good tracer of the innermost region of the actively accreting disk. We derive a mass accretion rate of 2-10x10^{-7} M_solar/yr from the measured Br_gamma emission luminosity of 1.78(+-0.31)x10^{31} erg/s.

We have analyzed the temperature, velocity and density of H2 gas in NGC 7023 with a high-resolution near-infrared spectrum of the northwestern filament of the reflection nebula. By observing NGC 7023 in the H and K bands at R ~ 45,000 with the Immersion GRating INfrared Spectrograph (IGRINS), we detected 68 H2 emission lines within the 1" x 15" slit. The diagnostic ratios of 2-1 S(1)/1-0 S(1) is 0.41-0.56. In addition, the estimated ortho-to-para ratios (OPR) is 1.63-1.82, indicating that the H2 emission transitions in the observed region arises mostly from gas excited by UV fluorescence. Gradients in the temperature, velocity, and OPR within the observed area imply motion of the photodissociation region (PDR) relative to the molecular cloud. In addition, we derive the column density of H2 from the observed emission lines and compare these results with PDR models in the literature covering a range of densities and incident UV field intensities. The notable difference between PDR model predictions and the observed data, in high rotational J levels of v = 1, is that the predicted formation temperature for newly-formed H2 should be lower than that of the model predictions. To investigate the density distribution, we combine pixels in 1" x 1" areas and derive the density distribution at the 0.002 pc scale. The derived gradient of density suggests that NGC 7023 has a clumpy structure, including a high clump density of ~10^5 cm^-3 with a size smaller than ~5 x 10^-3 pc embedded in lower density regions of 10^3-10^4 cm^-3.

We report a study of the three-dimensional (3D) outflow structure of a 15$\arcsec$ $\times$ 13$\arcsec$ area around H$_{2}$ peak 1 in Orion KL with slit-scan observations (13 slits) using the Immersion Grating Infrared Spectrograph. The datacubes, with high velocity-resolution ($\sim$ 7.5 {\kms}) provide high contrast imaging within ultra-narrow bands, and enable the detection of the main stream of the previously reported H$_{2}$ outflow fingers. We identified 31 distinct fingers in H$_{2}$ 1$-$0 S(1) $\lambda$2.122 $\micron$ emission. The line profile at each finger shows multiple-velocity peaks with a strong low-velocity component around the systemic velocity at ${\VLSR}$ = $+$8 {\kms} and high velocity emission ($|$${\VLSR}$$|$ = 45$-$135 {\kms}) indicating a typical bow-shock. The observed radial velocity gradients of $\sim$ 4 {\kms} arcsec$^{-1}$ agree well with the velocities inferred from large-scale proper motions, where the projected motion is proportional to distance from a common origin. We construct a conceptual 3D map of the fingers with the estimated inclination angles of 57$\degree$$-$74$\degree$. The extinction difference ($\Delta$$A_{\rm v}$ $>$ 10 mag) between blueshifted and redshifted fingers indicates high internal extinction. The extinction, the overall angular spread and scale of the flow argue for an ambient medium with very high density (10$^{5}$$-$10$^{6}$ cm$^{-3}$), consistent with molecular line observations of the OMC core. The radial velocity gradients and the 3D distributions of the fingers together support the hypothesis of simultaneous, radial explosion of the Orion KL outflow.

We present the results of high-resolution near-IR spectroscopy toward the multiple outflows around the Herbig Be star Lk{\Ha} 234 using the Immersion Grating Infrared Spectrograph (IGRINS). Previous studies indicate that the region around Lk{\Ha} 234 is complex, with several embedded YSOs and the outflows associated with them. In simultaneous H$-$ and K$-$band spectra from HH 167, we detected 5 {\FeII} and 14 H$_{2}$ emission lines. We revealed a new {\FeII} jet driven by radio continuum source VLA 3B. Position-velocity diagrams of H$_{2}$ 1$-$0 S(1) $\lambda$2.122 $\micron$ line show multiple velocity peaks. The kinematics may be explained by a geometrical bow shock model. We detected a component of H$_{2}$ emission at the systemic velocity (V$_{LSR}$ $=$ $-$10.2 {\kms}) along the whole slit in all slit positions, which may arise from the ambient photodissociation region. Low-velocity gas dominates the molecular hydrogen emission from knots A and B in HH 167, which is close to the systemic velocity, {\FeII} emission lines are detected at farther from the systemic velocity, at V$_{LSR}$ $=$ $-$100 $-$ $-$130 {\kms}. We infer that the H$_{2}$ emission arises from shocked gas entrained by a high-velocity outflow. Population diagrams of H$_{2}$ lines imply that the gas is thermalized at a temperature of 2,500 $-$ 3,000 K and the emission results from shock excitation.

The RGB Bayer filter system consists of a mosaic of R, G, and B filters on the grid of the photo sensors which typical commercial DSLR (Digital Single Lens Reflex) cameras and CCD cameras are equipped with. Lot of unique astronomical data obtained using an RGB Bayer filter system are available, including transient objects, e.g. supernovae, variable stars, and solar system bodies. The utilization of such data in scientific research requires that reliable photometric transformation methods are available between the systems. In this work, we develop a series of equations to convert the observed magnitudes in the RGB Bayer filter system ($R_B$, $G_B$, and $B_B$) into the Johnson-Cousins BVR filter system ($B_J$, $V_J$, and $R_C$). The new transformation equations derive the calculated magnitudes in the Johnson-Cousins filters ($B_{Jcal}$, $V_{Jcal}$, and $R_{Ccal}$) as functions of RGB magnitudes and colors. The mean differences between the transformed magnitudes and original magnitudes, i.e. the residuals, are $\Delta(B_J-B_{Jcal})$ = 0.064 mag, $\Delta(V_J-V_{Jcal})$ = 0.041 mag, and $\Delta(R_C-R_{Ccal})$ = 0.039 mag. The calculated Johnson-Cousins magnitudes from the transformation equations show a good linear correlation with the observed Johnson-Cousins magnitudes.

We present an exposure-time calculator (ETC) for the Immersion Grating Infrared Spectrograph (IGRINS). The signal and noise values are calculated by taking into account the telluric background emission and absorption, the emission and transmission of the telescope and instrument optics, and the dark current and read noise of the infrared detector arrays. For the atmospheric transmission, we apply models based on the amount of precipitable water vapor along the line of sight to the target. The ETC produces the expected signal-to-noise ratio (S/N) for each resolution element, given the exposure-time and number of exposures. In this paper, we compare the simulated continuum S/N for the early-type star HD 124683 and the late-type star GSS 32, and the simulated emission line S/N for the H2 rovibrational transitions from the Iris Nebula NGC 7023 with the observed IGRINS spectra. The simulated S/N from the ETC is overestimated by 40-50% for the sample continuum targets.

We present medium resolution near-infrared host galaxy spectra of low redshift quasars, PG 0844 + 349 (z=0.064), PG 1226 + 023 (z=0.158), and PG 1426+015 (z=0.086). The observations were done by using the Infrared Camera and Spectrograph (IRCS) at the Subaru 8.2 m telescope. The full width at half maximum of the point spread function was about 0.3 arcsec by operations of an adaptive optics system, which can effectively resolve the quasar spectra from the host galaxy spectra. We spent up to several hours per target and developed data reduction methods to reduce the systematic noises of the telluric emissions and absorptions. From the obtained spectra, we identified absorption features of Mg I (1.503 um), Si I (1.589 um) and CO (6-3) (1.619 um), and measured the velocity dispersions of PG 0844 + 349 to be 132+/-110 km s-1 and PG 1426 + 015 to be 264+/-215 km s-1. By using an M_BH-sigma relation of elliptical galaxies, we derived the black hole (BH) mass of PG 0844+349, log(M_BH/M_SUN) = 7.7+/-5.5 and PG 1426+015, log(M_BH/M_SUN) = 9.0+/-7.5. These values are consistent with the BH mass values from broad emission lines with an assumption of a virial factor of 5.5.

We present medium resolution (R = 5000-6000) spectra in the near-infrared band, 1.4-1.8 microns, for template stars in G, K, and M types observed by the echelle spectrometer, IRCS, at the SUBARU 8.2 m telescope. The identification of lines is based on the spectra of Arcturus (K2 III) in the literature. We measured the equivalent of widths and compared our results to those of Meyer et al. (1998). We conclude that our spectral resolution (R = 6000) data can investigate more accurately the properties of lines in stellar spectra. The library of the template stellar spectra in ASCII format are available for download on the World Wide Web.