We performed simultaneous observations of the H2O 6(1,6) - 5(2,3) (22.235080 GHz) and SiO v= 1, 2, J = 1 - 0, SiO v = 1, J = 2 - 1, 3 - 2 (43.122080, 42.820587, 86.243442, and 129.363359 GHz) masers towards the suspected D-type symbiotic star, V627 Cas, using the Korean VLBI Network. Here, we present astrometrically registered maps of the H2O and SiO v = 1, 2, J = 1 - 0, SiO v = 1, J = 2 - 1 masers for five epochs from January 2016 to June 2018. Distributions of the SiO maser spots do not show clear ring-like structures, and those of the H2O maser are biased towards the north-north-west to west with respect to the SiO maser features according to observational epochs. These asymmetric distributions of H2O and SiO masers are discussed based on two scenarios of a bipolar outflow and the presence of the hot companion, a white dwarf, in V627 Cas. We carried out ring fitting of SiO v = 1, and v = 2 masers and estimated the expected position of the cool red giant. The ring radii of the SiO v = 1 maser are slightly larger than those of the SiO v = 2 maser, as previously known. Our assumption for the physical size of the SiO maser ring of V627 Cas to be the typical size of a SiO maser ring radius (\sim4 au) of red giants yields the distance of V627 Cas to be \sim1 kpc.

H2O (22 GHz) and SiO masers (43, 86, 129 GHz) in the bipolar proto-planetary nebula OH231.8+4.2 were simultaneously monitored using the 21-m antennas of the Korean VLBI Network in 2009-2015. Both species exhibit periodic flux variations that correlate with the central star's optical light curve, with a phase delay of up to 0.15 for the maser flux variations with respect to the optical light curve. The flux densities of SiO v = 2, J = 1-0 and H2O masers decrease with time, implying that they may disappear in 10-20 years. However, there seems to have been a transient episode of intense H2O maser emission around 2010. We also found a systematic behaviour in the velocity profiles of these masers. The velocities of the H2O maser components appear to be remarkably constant, suggesting ballistic motion for the bipolar outflow in this nebula. On the other hand, those of the SiO maser clumps show a systematic radial acceleration of the individual clumps, converging to the outflow velocity of the H2O maser clumps. Measuring the full widths at zero power of the detected lines, we estimated the expansion velocities of the compact bipolar outflow traced by H2O maser and SiO thermal line, and discussed the possibility of the expanding SiO maser region in the equatorial direction. All of our analyses support that the central host star of OH231.8 is close to the tip of the AGB phase, and that the mass-loss rate recently started to decrease because of incipient post-AGB evolution.

We present the jet kinematics of the flat spectrum radio quasar (FSRQ) 4C +21.35 using time-resolved KaVA very long baseline interferometry array radio maps obtained from September 2014 to July 2016. During two out of three observing campaigns, observations were performed bi-weekly at 22 and 43 GHz quasi-simultaneously. At 22 GHz, we identified three jet components near the core with apparent speeds up to (14.4+/-2.1)c. The timing of the ejection of a new component detected in 2016 is consistent with a gamma-ray flare in November 2014. At 43 GHz, we found four inner jet (<3 mas) components with speeds from (3.5+/-1.4)c to (6.8+/-1.5)c. Jet component speeds tend to be higher with increasing distances from the core. We compared our data with archival Very Long Baseline Array (VLBA) data from the Boston University (BU) 43 GHz and the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) 15.4 GHz monitoring programs. Whereas MOJAVE data and our data are in good agreement, jet speeds obtained from the BU Program data in the same time period are about twice as high as the ones we obtain from the KaVA data. The discrepancy at 43 GHz indicates that radio arrays with different angular resolution identify and trace different jet features even when the data are obtained at the same frequency and at the same time. The flux densities of jet components decay exponentially, in agreement with a synchrotron cooling time scale of about 1 year. Using known electron Lorentz factor values (about 9,000), we estimate the magnetic field strength to be around 1-3 micro-Tesla. When adopting a jet viewing angle of 5 degrees, the intrinsic jet speed is of order 0.99c.

We obtained, for the first time, astrometrically registered maps of the 22.2 GHz H2O and 42.8, 43.1, and 86.2 GHz SiO maser emission toward the semiregular b-type variable (SRb) R Crateris, at three epochs (2015 May 21, and 2016 January 7 and 26) using the Korean Very-long-baseline Interferometry Network. The SiO masers show a ring-like spatial structure, while the H2O maser shows a very asymmetric one-side outflow structure, which is located at the southern part of the ring-like SiO maser feature. We also found that the 86.2 GHz SiO maser spots are distributed in an inner region, compared to those of the 43.1 GHz SiO maser, which is different from all previously known distributions of the 86.2 GHz SiO masers in variable stars. The different distribution of the 86.2 GHz SiO maser seems to be related to the complex dynamics caused by the overtone pulsation mode of the SRb R Crateris. Furthermore, we estimated the position of the central star based on the ring fitting of the SiO masers, which is essential for interpreting the morphology and kinematics of a circumstellar envelope. The estimated stellar coordinate corresponds well to the position measured by Gaia.

The supergiant VX Sagittarii is a strong emitter of both H2O and SiO masers. However, previous VLBI observations have been performed separately, which makes it difficult to spatially trace the outward transfer of the material consecutively. Here we present the astrometrically registered, simultaneous maps of 22.2 GHz H2O and 43.1/42.8/86.2/129.3 GHz SiO masers toward VX Sagittarii. The H2O masers detected above the dust-forming layers have an asymmetric distribution. The multi-transition SiO masers are nearly circular ring, suggesting spherically symmetric wind within a few stellar radii. These results provide the clear evidence that the asymmetry in the outflow is enhanced after the smaller molecular gas clump transform into the inhomogeneous dust layers. The 129.3 GHz maser arises from the outermost region compared to that of 43.1/42.8/86.2 GHz SiO masers. The ring size of the 129.3 GHz maser is maximized around the optical maximum, suggesting that radiative pumping is dominant.

We report on the astrometric registration of VLBI images of the SiO and water masers in OH231.8+4.2, the iconic Proto-Planetary Nebula also known as the Calabash nebula, using the KVN and Source/Frequency Phase Referencing. This, for the first time, robustly confirms the alignment of the SiO masers, close to the AGB star, driving the bi-lobe structure with the water masers in the out-flow. We are able to trace the bulk motions for the water masers over the last few decades to be 19 km/s and deduce that the age of this expansion stage is 38$\pm$2 years. The combination of this result with the distance allows a full 3D reconstruction, and confirms that the water masers lie on and expand along the known large-scale symmetry axis and that the outflow is only a few decades old, so mass loss is almost certainly on-going. Therefore we conclude that the SiO emission marks the stellar core of the nebular, the water emission traces the expansion, and that there must be multiple epochs of ejection to drive the macro-scale structure.

We report the initial results of our high-cadence monitoring program on the radio jet in the active galaxy M87, obtained by the KVN and VERA Array (KaVA) at 22 GHz. This is a pilot study that preceded a larger KaVA-M87 monitoring program, which is currently ongoing. The pilot monitoring was mostly performed every two to three weeks from December 2013 to June 2014, at a recording rate of 1 Gbps, obtaining the data for a total of 10 epochs. We successfully obtained a sequence of good quality radio maps that revealed the rich structure of this jet from <~1 mas to 20 mas, corresponding to physical scales (projected) of ~0.1-2 pc (or ~140-2800 Schwarzschild radii). We detected superluminal motions at these scales, together with a trend of gradual acceleration. The first evidence for such fast motions and acceleration near the jet base were obtained from recent VLBA studies at 43 GHz, and the fact that very similar kinematics are seen at a different frequency and time with a different instrument suggests these properties are fundamental characteristics of this jet. This pilot program demonstrates that KaVA is a powerful VLBI array for studying the detailed structural evolution of the M87 jet and also other relativistic jets.

Most bipolar outflows are associated with individual young stellar objects and have small opening angles. Here we report the discovery of an extremely wide-angle ($\sim$180$\arcdeg$) bipolar outflow ("EWBO") in a cluster forming region AFGL 5142 from low-velocity emission of the HCN (3-2) and HCO$^{+}$ (3-2) lines. This bipolar outflow is along a north-west to south-east direction with a line-of-sight flow velocity of about 3 km~s$^{-1}$ and is spatially connected to the high-velocity jet-like outflows. It seems to be a collection of low-velocity material entrained by the high-velocity outflows due to momentum feedback. The total ejected mass and mass loss rate due to both high velocity jet-like outflows and the "EWBO" are $\sim$24.5 M$_{\sun}$ and $\sim1.7\times10^{-3}$ M$_{\sun}$~yr$^{-1}$, respectively. Global collapse of the clump is revealed by the "blue profile" in the HCO$^{+}$ (1-0) line. A hierarchical network of filaments was identified in NH$_{3}$ (1,1) emission. Clear velocity gradients of the order of 10 km~s$^{-1}$~pc$^{-1}$ are found along filaments, indicating gas inflow along the filaments. The sum of the accretion rate along filaments and mass infall rate along the line of sight is $\sim$3.1$\times10^{-3}$ M$_{\sun}$~yr$^{-1}$, which exceeds the total mass loss rate, indicating that the central cluster is probably still gaining mass. The central cluster is highly fragmented and 22 condensations are identified in 1.1 mm continuum emission. The fragmentation process seems to be determined by thermal pressure and turbulence. The magnetic field may not play an important role in fragmentation.

Simultaneous time monitoring observations of H$_{2}$O and SiO maser lines were performed toward the D-type symbiotic binary system V407 Cyg with the Korean VLBI Network single dish radio telescope. These monitoring observations were carried out from March 2, 2010 (optical phase $\phi$ = 0.0), 8 days before the nova outburst on March 10, 2010 to June 5, 2014 ($\phi$ = 2.13). Eight days before the nova outburst, we detected the SiO $v$ = 1, 2, $J$ = 1--0 maser lines which exhibited values of 0.51 K ($\sim$ 6.70 Jy) and 0.71 K ($\sim$ 9.30 Jy), respectively, while after the outburst we could not detect them on April 2 ($\phi$ = 0.04), May 5 ($\phi$ = 0.09), May 8 ($\phi$ = 0.09), or on June 5, 2010 ($\phi$ = 0.13) within the upper limits of our KVN observations. After restarting our monitoring observations, we detected SiO $v$ = 2, $J$ = 1--0 masers starting on October 20, 2011 ($\phi$ = 0.83) and detected SiO $v$ = 1, $J$ = 1--0 masers starting on December 22, 2011 ($\phi$ = 0.92). These results provide clear evidence of the interaction between the shock from the nova outburst and the SiO maser regions of the Mira envelope. The peak emission of SiO $v$ = 1, 2, $J$ = 1--0 masers always occurred at blueshifted velocities with respect to the stellar velocity except for that of SiO $v$ = 1 at one epoch. These phenomena may be related to the redistribution of SiO maser regions after the outburst. The peak velocity variations of SiO masers associated with stellar pulsation phases show an increasing blueshifted trend during our monitoring interval after the outburst.

W49N is a mini-starburst in the Milky Way and thus an ideal laboratory for high-mass star formation studies. Due to its large distance (11.1$_{-0.7}^{+0.9}$ kpc), the kinematics inside and between the dense molecular clumps in W49N are far from well understood. The SMA observations resolved the continuum emission into two clumps. The molecular line observation of SO$_{2}$ (28$_{4,24}$-28$_{3,25}$) suggests that the two clumps have a velocity difference of $\sim$7 km~s$^{-1}$. The eastern clump is very close to two radio sources "G1" and "G2", and the western clump coincides with a radio source "B". The HCN (3-2) line reveals an extremely energetic outflow, which is among the most energetic molecular outflows in the Milky Way. This is the first report of high-velocity molecular outflow detection in W49N. The outflow jet might be in precession, which could account for the distribution, velocity and rotation of water maser spots. Three absorption systems are identified in HCO$^{+}$ (3-2) spectra. The absorption features are blueshifted with respect to the emission of SO$_{2}$ (28$_{4,24}$-28$_{3,25}$) lines, indicating that a cold layer is expanding in front of the warm gas. Further analysis indicates that the expansion is decelerated from the geometric expansion centers.

Simultaneous time monitoring observations of H$_{2}$O $6_{16}-5_{23}$, SiO $J$ = 1--0, 2--1, 3--2, and $^{29}$SiO $v$ = 0, $J$ = 1--0 lines were carried out in the direction of the Mira variable star TX Cam with the Korean VLBI Network single dish radio telescopes. For the first time, the H$_{2}$O maser emission from TX Cam was detected near the stellar velocity at five epochs from April 10, 2013 ($\phi$ = 3.13) to June 4, 2014 ($\phi$ = 3.89) including minimum optical phases. The intensities of H$_{2}$O masers are very weak compared to SiO masers. The variation of peak antenna temperature ratios among SiO $v$ = 1, $J$ = 1--0, $J$ = 2--1, and $J$ = 3--2 masers is investigated according to their phases. The shift of peak velocities of H$_{2}$O and SiO masers with respect to the stellar velocity is also investigated according to observed optical phases. The H$_{2}$O maser emission occurs around the stellar velocity during our monitoring interval. On the other hand, the peak velocities of SiO masers show a spread compared to the stellar velocity. The peak velocities of SiO $J$ = 2--1, and $J$ = 3--2 masers show a smaller spread with respect to the stellar velocity than those of SiO $J$ = 1--0 masers. These simultaneous observations of multi-frequencies will provide a good constraint for maser pumping models and a good probe for investigating the stellar atmosphere and envelope according to their different excitation conditions.

Oxygen-rich Asymptotic Giant Branch (AGB) stars can be intense emitters of SiO ($v$=1 and 2, J=1$\rightarrow$0) and H$_2$O maser lines at 43 and 22 GHz, respectively. VLBI observations of the maser emission provide a unique tool to probe the innermost layers of the circumstellar envelopes in AGB stars. Nevertheless, the difficulties in achieving astrometrically aligned \water\ and $v$=1 and $v$=2 SiO maser maps have traditionally limited the physical constraints that can be placed on the SiO maser pumping mechanism. We present phase referenced simultaneous spectral-line VLBI images for the SiO $v$=1 and $v$=2, J=1$\rightarrow$0, and H$_2$O maser emission around the AGB star R\,LMi, obtained from the Korean VLBI Network (KVN). The simultaneous multi-channel receivers of the KVN offer great possibilities for astrometry in the frequency domain. With this facility we have produced images with bona-fide absolute astrometric registration between high frequency maser transitions of different species to provide the positions of the \water\ maser emission, and the centre of the SiO maser emission, and hence reducing the uncertainty in the proper motion for R\,LMi by an order of magnitude over that from Hipparcos. This is the first successful demonstration of source frequency phase referencing for mm-VLBI spectral-line observations and also where the ratio between the frequencies is not an integer.

The Korean VLBI Network (KVN) is a new mm-VLBI dedicated array with capability for simultaneous observations at multiple frequencies, up to 129 GHz. The innovative multi-channel receivers present significant benefits for astrometric measurements in the frequency domain. The aim of this work is to verify the astrometric performance of the KVN using a comparative study with the VLBA, a well established instrument. For that purpose, we carried out nearly contemporaneous observations with the KVN and the VLBA, at 14/7 mm, in April 2013. The KVN observations consisted of simultaneous dual frequency observations, while the VLBA used fast frequency switching observations. We used the Source Frequency Phase Referencing technique for the observational and analysis strategy. We find that having simultaneous observations results in a superior performance for compensation of all atmospheric terms in the observables, in addition to offering other significant benefits for astrometric analysis. We have compared the KVN astrometry measurements to those from the VLBA. We find that the structure blending effects introduce dominant systematic astrometric shifts and these need to be taken into account. We have tested multiple analytical routes to characterize the impact of the low resolution effects for extended sources in the astrometric measurements. The results from the analysis of KVN and full VLBA datasets agree within 2-$\sigma$ of the thermal error estimate. We interpret the discrepancy as arising from the different resolutions. We find that the KVN provides astrometric results with excellent agreement, within 1-$\sigma$, when compared to a VLBA configuration which has a similar resolution. Therefore this comparative study verifies the astrometric performance of KVN using SFPR at 14/7 mm, and validates the KVN as an astrometric instrument.

The Korean very-long-baseline interferometry (VLBI) network (KVN) and VLBI Exploration of Radio Astrometry (VERA) Array (KaVA) is the first international VLBI array dedicated to high-frequency (23 and 43 GHz bands) observations in East Asia. Here, we report the first imaging observations of three bright active galactic nuclei (AGNs) known for their complex morphologies: 4C 39.25, 3C 273, and M 87. This is one of the initial result of KaVA early science. Our KaVA images reveal extended outflows with complex substructure such as knots and limb brightening, in agreement with previous Very Long Baseline Array (VLBA) observations. Angular resolutions are better than 1.4 and 0.8 milliarcsecond at 23 GHz and 43 GHz, respectively. KaVA achieves a high dynamic range of ~1000, more than three times the value achieved by VERA. We conclude that KaVA is a powerful array with a great potential for the study of AGN outflows, at least comparable to the best existing radio interferometric arrays.

We have carried out the first very long baseline interferometry (VLBI) imaging of 44 GHz class I methanol maser (7_{0}-6_{1}A^{+}) associated with a millimeter core MM2 in a massive star-forming region IRAS 18151-1208 with KaVA (KVN and VERA Array), which is a newly combined array of KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). We have succeeded in imaging compact maser features with a synthesized beam size of 2.7 milliarcseconds x 1.5 milliarcseconds (mas). These features are detected at a limited number of baselines within the length of shorter than approximately 650 km corresponding to 100 Mlambda in the uv-coverage. The central velocity and the velocity width of the 44 GHz methanol maser are consistent with those of the quiescent gas rather than the outflow traced by the SiO thermal line. The minimum component size among the maser features is ~ 5 mas x 2 mas, which corresponds to the linear size of ~ 15 AU x 6 AU assuming a distance of 3 kpc. The brightness temperatures of these features range from ~ 3.5 x 10^{8} to 1.0 x 10^{10} K, which are higher than estimated lower limit from a previous Very Large Array observation with the highest spatial resolution of ~ 50 mas. The 44 GHz class I methanol maser in IRAS 18151-1208 is found to be associated with the MM2 core, which is thought to be less evolved than another millimeter core MM1 associated with the 6.7 GHz class II methanol maser.

Multi-epoch phase-referencing VLBI (Very Long Baseline Interferometry) observations with VERA (VLBI Exploration of Radio Astrometry) were performed for the symbiotic star R Aquarii (R Aqr) from September 2005 to Oct 2006. Tracing one of the $v=2$, $J=1-0$ SiO maser spots, we measured an annual parallax of $\pi = 4.59\pm0.24$ mas, corresponding to a distance of $218_{-11}^{+12}$ pc. Our result is consistent with earlier distance measurements, but yields the highest accuracy of about $5\%$ level. Applying our distance, we derived an absolute K-band magnitude of $M_{\mathrm{K}} = -7.71 \pm 0.11$, which is consistent with the recent Period-Luminosity relation by VLBI parallax measurements for 5 OH-Mira variables. In addition, the expansion age of an inner nebulae around R Aqr is found to be about 240 years, corresponds to about the year 1773.

In this Letter, we report detections of SiO v=3 J=1--0 maser emission in very long baseline interferometric (VLBI) observations towards 4 out of 12 long-period variable stars: WX Psc, R Leo, W Hya, and T Cep. The detections towards WX Psc and T Cep are new ones. We also present successful astrometric observations of SiO v=2 and v=3 J=1--0 maser emissions associated with two stars: WX Psc and W Hya and their position-reference continuum sources: J010746.0+131205 and J135146.8-291218 with the VLBI Exploration of Radio Astrometry (VERA). The relative coordinates of the position-reference continuum source and SiO v=3 maser spots were measured with respect to those of an SiO v=2 maser spot adopted as fringe-phase reference. Thus the faint continuum sources were inversely phase-referenced to the bright maser sources. It implies possible registration of multiple SiO maser line maps onto a common coordinate system with 10 microarcsecond-level accuracy.

We report simultaneous multi-frequency observing performance at 22 and 43 GHz of the 21-m shaped-Cassegrain radio telescopes of the Korean VLBI Network (KVN). KVN is the first millimeter-dedicated VLBI network in Korea having a maximum baseline length of 480 km. It currently operates at 22 and 43 GHz and planed to operate in four frequency bands, 22, 43, 86, and 129 GHz. The unique quasioptics of KVN enable simultaneous multi-frequency observations based on efficient beam filtering and accuarate antenna-beam alignment at 22 and 43 GHz. We found that the offset of the beams is within <5 arcseconds over all pointing directions of antenna. The dual polarization, cooled HEMT receivers at 22 and 43 GHz result in receiver noise temperatures less than 40 K at 21.25-23.25 GHz and 80 K at 42.11-44.11 GHz. The pointing accuracies have been measured to be 3 arcseconds in azimuth and elevation for all antennas. The measured aperture efficiencies are 65%(K)/67%(Q), 62%(K)/59%(Q), and 66%(K)/60%(Q) for the three KVN antennas, KVNYS, KVNUS, and KVNTN, respectively. The main-beam efficiencies are measured to be 50%(K)/52%(Q), 48%(K)/50%(Q), and 50%(K)/47%(Q) for KVNYS, KVNUS, and KVNTN, respectively. The estimated Moon efficiencies are 77%(K)/90%(Q), 74%(K)/79%(Q), and 80%(K)/86%(Q) for KVNYS, KVNUS, KVNTN, respectively. The elevation dependence of the aperture efficiencies is quite flat for elevations > 20 degrees.

In this paper, we report results of our near-infrared (NIR) photometric variability studies of the BL Lacertae object S5 0716+714. NIR photometric observations spread over 7 nights during our observing run April 2-9, 2007 at 1.8 meter telescope equipped with KASINICS (Korea Astronomy and Space Science Institute Near Infrared Camera System) and J, H, and Ks filters at Bohyunsan Optical Astronomy Observatory (BOAO), South Korea. We searched for intra-day variability, short term variability and color variability in the BL Lac object. We have not detected any genuine intra-day variability in any of J, H, and Ks passbands in our observing run. Significant short term variability ~ 32.6%, 20.5% and 18.2% have been detected in J, H, Ks passbands, respectively, and ~ 11.9% in (J-H) color.

R Aqr is known to be a symbiotic binary system with an associated extended emission nebula, possibly produced by a historic outburst. To find the associated historic records, we searched for and compiled all 'Guest Star' and 'Peculiar Star' records in three Korean 'official' history books that cover almost two thousand years, Samguksagi, Goryeosa, Joseonwangjosillok. In addition to the record of A.D. 1073, previously noted by Li (1985), we have found in Goryeosa another candidate record of A.D. 1074, which has the same positional description as that of A.D. 1073 with an additional brightness description. We examined various aspects of the two records and conclude that they both are likely to be the records of outburst of R Aqr. This means that there were two successive outbursts in A.D. 1073 and in A.D. 1074, separated by approximately one year. Based on these records, we estimate the distance to R Aqr to be 273 pc if the expansion of the nebula has been at a constant rate. The brightness record of A.D. 1074 corresponds to the absolute magnitude at outburst of M_(outburst) = -6.2 mag. ~ -5.2 mag. at this distance. The two Korean records associated with R Aqr may provide astronomically meaningful constraints to the outburst model of R Aqr and the formative process of its nebulosity.

We present the results of a spectral line survey from 138.3 to 150.7 GHz toward Orion-KL. The observations were made using the 14 m radio telescope of Taeduk Radio Astronomy Observatory. Typical system temperatures were between 500 and 700 K, with the sensitivity between $0.02 - 0.06$ K in units of $\rm T_A^*$. A total of 149 line spectra are detected in this survey. Fifty lines have been previously reported, however we find 99 new detections. Among these new lines, 32 are `unidentified', while 67 are from molecular transitions with known identifications. There is no detection of H or He recombination lines. The identified spectra are from a total of 16 molecular species and their isotopic variants. In the range from 138.3 to 150.7 GHz, the strongest spectral line is the J=3-2 transition of CS molecule, followed by transitions of the $\rm H_2CO$, $\rm CH_3OH$, $\rm CH_3CN$, and $\rm SO_2$. Spectral lines from the large organic molecules such as $\rm CH_3OH$, $\rm CH_3OCH_3$, $\rm HCOOCH_3$, $\rm C_2H_5CN$ and $\rm CH_3CN$ are prominent; with 80 % of the identified lines arising from transitions of these molecules. The rotational temperatures and column densities are derived using the standard rotation diagram analysis for $\rm CH_3OH$ ($\rm ^{13}CH_3OH$), $\rm HCOOCH_3$, $\rm CH_3CN$ and $\rm SO_2$ with $\rm 10\sim 270 K$ and $\rm 0.2\sim 20\times 10^{15} cm^{-2}$. These estimates are fairly comparable to the values for the same molecule in other frequency regions by other studies.