We present the results of high-resolution (R $\ge$ 30,000) optical and near-infrared spectroscopic monitoring observations of a FU Orionis-type object, V960 Mon, which underwent an outburst in 2014 November. We have monitored this object with the Bohyunsan Optical Echelle Spectrograph (BOES) and the Immersion GRating INfrared Spectrograph (IGRINS) since 2014 December. Various features produced by a wind, disk, and outflow/jet were detected. The wind features varied over time and continually weakened after the outburst. We detected double-peaked line profiles in the optical and near-infrared, and the line widths tend to decrease with increasing wavelength, indicative of Keplerian disk rotation. The disk features in the optical and near-infrared spectra fit well with G-type and K-type stellar spectra convolved with a kernel to account for the maximum projected disk rotation velocity of about 40.3$\pm$3.8 km s$^{-1}$ and 36.3$\pm$3.9 km s$^{-1}$, respectively. We also report the detection of [S II] and H$_{2}$ emission lines, which are jet/outflow tracers and rarely found in FUors.

HD 81817 is known as a hybrid star. Hybrid stars have both cool stellar wind properties and Ultraviolet (UV) or even X-ray emission features of highly ionized atoms in their spectra. A white dwarf companion has been suggested as the source of UV or X-ray features. HD 81817 has been observed since 2004 as a part of a radial velocity (RV) survey program to search for exoplanets around K giant stars using the Bohyunsan Observatory Echelle Spectrograph at the 1.8 m telescope of Bohyunsan Optical Astronomy Observatory in Korea. We obtained 85 RV measurements between 2004 and 2019 for HD 81817 and found two periodic RV variations. The amplitudes of RV variations are around 200 m s^-1, which are significantly lower than that expected from a closely orbiting white dwarf companion. Photometric data and relevant spectral lines were also analyzed to help determine the origin of the periodic RV variations. We conclude that 627.4-day RV variations are caused by intrinsic stellar activities such as long-term pulsations or rotational modulations of surface activities based on H{\alpha} equivalent width (EW) variations of a similar period. On the other hand, 1047.1-day periodic RV variations are likely to be caused by a brown dwarf or substellar companion, which is corroborated by a recent GAIA proper motion anomaly for HD 81817. The Keplerian fit yields a minimum mass of 27.1 M_Jup, a semimajor axis of 3.3 AU, and an eccentricity of 0.17 for the stellar mass of 4.3 M_sun for HD 81817. The inferred mass puts HD 81817 b in the brown dwarf desert.

The cosmic origin of fluorine is still not well constrained. Several nucleosynthetic channels at different phases of stellar evolution have been suggested, but these must be constrained by observations. For this, the fluorine abundance trend with metallicity spanning a wide range is required. Our aim is to determine stellar abundances of fluorine for -1.1<[Fe/H]<+0.4. We determine the abundances from HF lines in infrared K-band spectra (at 2.3 microns) of cool giants, observed with the IGRINS and Phoenix high-resolution spectrographs. We derive accurate stellar parameters for all our observed K giants, which is important since the HF lines are very temperature sensitive. We find that [F/Fe] is flat as a function of metallicity at [F/Fe]~0, but increases as the metallicity increases. The fluorine slope shows a clear secondary behavior in this metallicity range. We also find that the [F/Ce] ratio is relatively flat for -0.6<[Fe/H]<0, and that for two metal-poor ([Fe/H]<-0.8), s-process element enhanced giants, we do not detect an elevated fluorine abundance. We interpret all these observational constraints to indicate that several major processes are at play for the cosmic budget of fluorine over time; from those in massive stars at low metallicities, through the asymptotic giant branch-star contribution at -0.6<[Fe/H]<0, to processes with increasing yields with metallicity at super-solar metallicities. The origins of the latter, and whether or not Wolf-Rayet stars and/or novae could contribute at super-solar metallicities, is currently not known. To quantify these observational results, theoretical modelling is required. More observations in the metal-poor region are required to clarify the processes there.

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.

Cassiopeia A is a nearby young supernova remnant that provides a unique laboratory for the study of core-collapse supernova explosions. Cassiopeia A is known to be a Type IIb supernova from the optical spectrum of its light echo, but the immediate progenitor of the supernova remains uncertain. Here we report results of near-infrared, high-resolution spectroscopic observations of Cassiopeia A where we detected the pristine circumstellar material of the supernova progenitor. Our observations revealed a strong emission line of iron (Fe) from a circumstellar clump that has not yet been processed by the supernova shock wave. A comprehensive analysis of the observed spectra, together with an HST image, indicates that the majority of Fe in this unprocessed circumstellar material is in the gas phase, not depleted onto dust grains as in the general interstellar medium. This result is consistent with a theoretical model of dust condensation in material that is heavily enriched with CNO-cycle products, supporting the idea that the clump originated near the He core of the progenitor. It has been recently found that Type IIb supernovae can result from the explosion of a blue supergiant with a thin hydrogen envelope, and our results support such a scenario for Cassiopeia A.

We describe the design and assembly of the LUX-ZEPLIN experiment, a direct detection search for cosmic WIMP dark matter particles. The centerpiece of the experiment is a large liquid xenon time projection chamber sensitive to low energy nuclear recoils. Rejection of backgrounds is enhanced by a Xe skin veto detector and by a liquid scintillator Outer Detector loaded with gadolinium for efficient neutron capture and tagging. LZ is located in the Davis Cavern at the 4850' level of the Sanford Underground Research Facility in Lead, South Dakota, USA. We describe the major subsystems of the experiment and its key design features and requirements.

High-resolution, near-infrared spectra will be the primary tool for finding and characterizing Earth-like planets around low-mass stars. Yet, the properties of exoplanets can not be precisely determined without accurate and precise measurements of the host star. Spectra obtained with the Immersion GRating INfrared Spectrometer (IGRINS) simultaneously provide diagnostics for most stellar parameters, but the first step in any analysis is the determination of the effective temperature. Here we report the calibration of high-resolution H-band spectra to accurately determine effective temperature for stars between 4000-3000 K ($\sim$K8--M5) using absorption line depths of Fe I, OH, and Al I. The field star sample used here contains 254 K and M stars with temperatures derived using BT-Settl synthetic spectra. We use 106 stars with precise temperatures in the literature to calibrate our method with typical errors of about 140 K, and systematic uncertainties less than $\sim$120 K. For the broadest applicability, we present T$_{\rm eff}$--line-depth-ratio relationships, which we test on 12 members of the TW Hydrae Association and at spectral resolving powers between $\sim$10,000--120,000. These ratios offer a simple but accurate measure of effective temperature in cool stars that is distance and reddening independent.

We present a high-resolution, near-IR spectroscopic study of multiple outflows in the LkH$\alpha$ 234 star formation region using the Immersion GRating INfrared Spectrometer (IGRINS). Spectral mapping over the blueshifted emission of HH 167 allowed us to distinguish at least three separate, spatially overlapped, outflows in H${_2}$ and [Fe II] emission. We show that the H${_2}$ emission represents not a single jet, but complex multiple outflows driven by three known embedded sources: MM1, VLA 2, and VLA 3. There is a redshifted H${_2}$ outflow at a low velocity, $\VLSR$ $<$ $+$50 {\kms}, with respect to the systemic velocity of $\VLSR$ $=$ $-$11.5 {\kms}, that coincides with the H${_2}$O masers seen in earlier radio observations two arcseconds southwest of VLA 2. We found that the previously detected [Fe II] jet with $|$$\VLSR$$|$ $>$ 100 {\kms} driven by VLA 3B is also detected in H${_2}$ emission, and confirm that this jet has a position angle about 240$\degree$. Spectra of the redshifted knots at 14$\arcsec$$-$65$\arcsec$ northeast of LkH$\alpha$ 234 are presented for the first time. These spectra also provide clues to the existence of multiple outflows. We detected high-velocity (50$-$120 {\kms}) H${_2}$ gas in the multiple outflows around LkH$\alpha$ 234. Since these gases move at speeds well over the dissociation velocity ($>$ 40 {\kms}), the emission must originate from the jet itself rather than H${_2}$ gas in the ambient medium. Also, position-velocity diagrams and excitation diagram indicate that emission from knot C in HH 167 come from two different phenomena, shocks and photodissociation.

OB associations are the prevailing star forming sites in the Galaxy. Up to now, the process of how OB associations were formed remained a mystery. A possible process is self-regulating star formation driven by feedback from massive stars. However, although a number of observational studies uncovered various signposts of feedback-driven star formation, the effectiveness of such feedback has been questioned. Stellar and gas kinematics is a promising tool to capture the relative motion of newborn stars and gas away from ionizing sources. We present high-resolution spectroscopy of stars and gas in the young open cluster NGC 1893. Our findings show that newborn stars and the tadpole nebula Sim 130 are moving away from the central cluster containing two O-type stars, and that the timescale of sequential star formation is about 1 Myr within a 9 parsec distance. The newborn stars formed by feedback from massive stars account for at least 18 per cent of the total stellar population in the cluster, suggesting that this process can play an important role in the formation of OB associations. These results support the self-regulating star formation model.

At 60 pc, TW Hydra (TW Hya) is the closest example of a star with a gas-rich protoplanetary disk, though TW Hya may be relatively old (3-15 Myr). As such, TW Hya is especially appealing to test our understanding of the interplay between stellar and disk evolution. We present a high-resolution near-infrared spectrum of TW Hya obtained with the Immersion GRating INfrared Spectrometer (IGRINS) to re-evaluate the stellar parameters of TW Hya. We compare these data to synthetic spectra of magnetic stars produced by MoogStokes, and use sensitive spectral line profiles to probe the effective temperature, surface gravity, and magnetic field. A model with T_eff= 3800 K, log g=4.2, and B=3.0 kG best fits the near-infrared spectrum of TW Hya. These results correspond to a spectral type of M0.5 and an age of 8 Myr, which is well past the median life of gaseous disks.

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 present the detection of long-period RV variations in HD 36384, HD 52030, and HD 208742 by using the high-resolution, fiber-fed Bohyunsan Observatory Echelle Spectrograph (BOES) for the precise radial velocity (RV) survey of about 200 northern circumpolar stars. Analyses of RV data, chromospheric activity indicators, and bisector variations spanning about five years suggest that the RV variations are compatible with planet or brown dwarf companions in Keplerian motion. However, HD 36384 shows photometric variations with a period very close to that of RV variations as well as amplitude variations in the weighted wavelet Z-transform (WWZ) analysis, which argues that the RV variations in HD~36384 are from the stellar pulsations. Assuming that the companion hypothesis is correct, HD~52030 hosts a companion with minimum mass 13.3 M_Jup$ orbiting in 484 days at a distance of 1.2 AU. HD~208742 hosts a companion of 14.0 M_Jup at 1.5 AU with a period of 602 days. All stars are located at the asymptotic giant branch (AGB) stage on the H-R diagram after undergone the helium flash and left the giant clump.With stellar radii of 53.0 R_Sun and 57.2 R_Sun for HD 52030 and HD 208742, respectively, these stars may be the largest yet, in terms of stellar radius, found to host sub-stellar companions. However, given possible RV amplitude variations and the fact that these are highly evolved stars the planet hypothesis is not yet certain.

We present a deep near-infrared spectrum of the Orion Bar Photodissociation Region (PDR) taken with the Immersion Grating INfrared Spectrometer (IGRINS) on the 2.7 m telescope at the McDonald Observatory. IGRINS has high spectral resolution (R~45000) and instantaneous broad wavelength coverage (1.45-2.45 microns), enabling us to detect 87 emission lines from rovibrationally excited molecular hydrogen (H_2) that arise from transitions out of 69 upper rovibration levels of the electronic ground state. These levels cover a large range of rotational and vibrational quantum numbers and excitation energies, making them an excellent probe of the excitation mechanisms of H_2 and physical conditions within the PDR. The Orion Bar PDR is thought to consist of cooler high density clumps or filaments (T=50-250 K, n_H = 10^5 - 10^7 cm^-3) embedded in a warmer lower density medium (T=250-1000 K, n_H=10^4 - 10^5 cm^-3). We fit a grid of simple constant-temperature and constant-density Cloudy models, which recreate the observed H_2 level populations well, to constrain the temperature to a range of 600 to 650 K and the density to n_H = 2.5 x 10^3 to 10^4 cm^-3. The best fit model gives T = 625 K and n_H = 5x10^3 cm^-3. This well constrained warm temperature is consistent with kinetic temperatures found by other studies for the Orion Bar's lower density medium. However, the range of densities well fit by the model grid is marginally lower than those reported by other studies. We could be observing lower density gas than the surrounding medium, or perhaps a density-sensitive parameter in our models is not properly estimated.

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.

From high resolution (R = 45,000), high signal-to-noise (S/N > 400) spectra gathered with the Immersion Grating Infrared Spectrograph (IGRINS) in the H and K photometric bands, we have derived elemental abundances of two bright, well-known metal-poor halo stars: the red giant HD 122563 and the subgiant HD 140283. Since these stars have metallicities approaching [Fe/H] = -3, their absorption features are generally very weak. Neutral-species lines of Mg, Si, S and Ca are detectable, as well as those of the light odd-Z elements Na and Al. The derived IR-based abundances agree with those obtained from optical-wavelength spectra. For Mg and Si the abundances from the infrared transitions are improvements to those derived from shorter wavelength data. Many useful OH and CO lines can be detected in the IGRINS HD 122563 spectrum, from which derived O and C abundances are consistent to those obtained from the traditional [O I] and CH features. IGRINS high resolutions H- and K-band spectroscopy offers promising ways to determine more reliable abundances for additional metal-poor stars whose optical features are either not detectable, or too weak, or are based on lines with analytical difficulties.

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.

Aims. This program originated as the north pole region extension of the established exoplanet survey using 1.8 m telescope at Bohyunsan Optical Astronomy Observatory (BOAO). The aim of our paper is to find exoplanets in northern circumpolar stars with a precise radial velocity (RV) survey. Methods. We have selected about 200 northern circumpolar stars with the following criteria: Dec. > 70 degree, 0.6 < B-V < 1.6, HIPPARCOS_scat < 0.05 magnitude, and 5.0 < mv < 7.0. The high-resolution, fiber-fed Bohyunsan Observatory Echelle Spectrograph (BOES) was used for the RV survey. Chromospheric activities, the HIPPARCOS photometry, and line bisectors were analyzed to exclude other causes for the RV variations. Results. In 2010, we started to monitor the candidates and have completed initial screening for all stars for the last five years. We present the detection of four new exoplanets. Stars HD 11755, HD 12648, HD 24064, and 8 UMi all show evidence for giant planets in Keplerian motion. The companion to HD 11755 has a minimum mass of 6.5 M_Jup in a 433-day orbit with an eccentricity of 0.19. HD 12648 is orbited by a companion of minimum mass of 2.9 M_Jup having a period of 133 days and an eccentricity of 0.04. Weak surface activity was suspected in HD 24064. However, no evidence was found to be associated with the RV variations. Its companion has a minimum mass of 9.4 M_Jup, a period of 535 days, and an eccentricity of 0.35. Finally, 8 UMi has a minimum mass of 1.5 M_Jup, a period of 93 days with an eccentricity of 0.06.

We present observational results from optical long-slit spectroscopy of parsec-scale jets of DG Tau. From HH 158 and HH 702, the optical emission lines of H_alpha, [O I] 6300, 6363, [N II] 6548, 6584, and [S II] 6716, 6731 are obtained. The kinematics and physical properties (i.e., electron density, electron temperature, ionization fraction, and mass-loss rate) are investigated along the blueshifted jet up to 650 arcseconds distance from the source. For HH 158, the radial velocity ranges from -50 to -250 km/s. The proper motion of the knots is 0.196 - 0.272 /yr. The electron density is ~ 10^{4}/cm^{3} close to the star, and decreases to ~ 10^{2} /cm^{3} at 14 arcseconds away from the star. Ionization fraction indicates that the gas is almost neutral in the vicinity of the source. It increases up to over 0.4 along the distance. HH 702 is located at ~ 650 arcseconds from the source. It shows ~ -80 km/s in the radial velocity. Its line ratios are similar to those at knot C of HH 158. The mass-loss rate is estimated to be about ~ 10^{-7} M_{sun}/yr, which is similar to values obtained from previous studies.

We present the results of extensive multi-waveband monitoring of the blazar 3C~279 between 1996 and 2007 at X-ray energies (2-10 keV), optical R band, and 14.5 GHz, as well as imaging with the Very Long Baseline Array (VLBA) at 43 GHz. In all bands the power spectral density corresponds to "red noise" that can be fit by a single power law over the sampled time scales. Variations in flux at all three wavebands are significantly correlated. The time delay between high and low frequency bands changes substantially on time scales of years. A major multi-frequency flare in 2001 coincided with a swing of the jet toward a more southerly direction, and in general the X-ray flux is modulated by changes in the position angle of the jet near the core. The flux density in the core at 43 GHz--increases in which indicate the appearance of new superluminal knots--is significantly correlated with the X-ray flux. We decompose the X-ray and optical light curves into individual flares, finding that X-ray leads optical variations (XO) in 6 flares, the reverse occurs in 3 flares (OX), and there is essentially zero lag in 4 flares. Upon comparing theoretical expectations with the data, we conclude that (1) XO flares can be explained by gradual acceleration of radiating electrons to the highest energies; (2) OX flares can result from either light-travel delays of the seed photons (synchrotron self-Compton scattering) or gradients in maximum electron energy behind shock fronts; and (3) events with similar X-ray and optical radiative energy output originate well upstream of the 43 GHz core, while those in which the optical radiative output dominates occur at or downstream of the core.