We present a simple, novel method of efficiently rendering ray cast soft shadows on curved terrain by using dynamic programming and maximum mipmaps to rapidly find a global minimum shadow cost in constant runtime complexity. Additionally, we apply a new method of reducing view ray computation times that pre-displaces the terrain mesh to bootstrap ray starting positions. Combining these two methods, our ray casting engine runs in real-time with more than 200% speed up over uniform ray stepping with comparable image quality and without hardware ray tracing acceleration. To add support for accurate planetary ephemerides and interactive features, we integrated the engine into celestia.Sci, a general space simulation software. We demonstrate the ability of our engine to accurately handle a large range of distance scales by using it to generate videos of lunar landing trajectories. The numerical error when compared with real lunar mission imagery is small, demonstrating the accuracy and efficiency of our approach.

We study the kinematics of the M87 jet using the first year data of the KVN and VERA Array (KaVA) large program, which has densely monitored the jet at 22 and 43 GHz since 2016. We find that the apparent jet speeds generally increase from $\approx0.3c$ at $\approx0.5$ mas from the jet base to $\approx2.7c$ at $\approx20$ mas, indicating that the jet is accelerated from subluminal to superluminal speeds on these scales. We perform a complementary jet kinematic analysis by using archival Very Long Baseline Array monitoring data observed in $2005-2009$ at 1.7 GHz and find that the jet is moving at relativistic speeds up to $\approx5.8c$ at distances of $200-410$ mas. We combine the two kinematic results and find that the jet is gradually accelerated over a broad distance range that coincides with the jet collimation zone, implying that conversion of Poynting flux to kinetic energy flux takes place. If the jet emission consists of a single streamline, the observed trend of jet acceleration ($\Gamma\propto z^{0.16\pm0.01}$) is relatively slow compared to models of a highly magnetized jet. This indicates that Poynting flux conversion through the differential collimation of poloidal magnetic fields may be less efficient than expected. However, we find a non-negligible dispersion in the observed speeds for a given jet distance, making it difficult to describe the jet velocity field with a single power-law acceleration function. We discuss the possibility that the jet emission consists of multiple streamlines following different acceleration profiles, resulting in jet velocity stratification.

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.

In the circumnuclear region of the nearby radio galaxy NGC~1052. Our $\lambda$3mm VLBI observations with the Korean VLBI Network have spatially resolved the broad HCO$^{+}$ absorption at a velocity range of 1350--1850 km~s$^{-1}$ against a double-sided nuclear jet, and have revealed that the HCO$^{+}$ absorption is concentrated on the receding jet and the nuclear components. The distribution of the HCO$^{+}$ absorbing gas strongly supports the circumnuclear torus surrounding the supermassive black hole. From the estimations of the column density and the volume density of molecular hydrogen, the size of the molecular gas region in the torus is at least 1~pc. The broad spectral profile of HCO$^{+}$ is likely to be a blend of multiple gas clumps with various velocities. The HCO$^{+}$ absorption of NGC~1052 could trace complex kinematics in the vicinity of the active galactic nucleus, such as inflow, outflow, turbulence, and so on.

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.

We present the first VLBI detection of HCN molecular absorption in the nearby active galactic nucleus NGC 1052. Utilizing the 1 milliarcsecond resolution achieved by the Korean VLBI Network, we have spatially resolved the HCN absorption against a double-sided nuclear jet structure. Two velocity features of HCN absorption are detected significantly at the radial velocity of 1656 and 1719 km/s, redshifted by 149 and 212 km/s with respect to the systemic velocity of the galaxy. The column density of the HCN molecule is estimated to be 10^{15}-10^{16} cm^{-2}, assuming the excitation temperature of 100-230 K. The absorption features show high optical depth localized on the receding jet side, where the free-free absorption occurred due to the circumnuclear torus. The size of the foreground absorbing molecular gas is estimated to be on approximately one-parsec scales, which agrees well with the approximate size of the circumnuclear torus. HCN absorbing gas is likely to be several clumps smaller than 0.1 parsec inside the circumnuclear torus. The redshifted velocities of the HCN absorption features imply that HCN absorbing gas traces ongoing infall motion inside the circumnuclear torus onto the central engine.

The Disk Detective citizen science project aims to find new stars with 22 micron excess emission from circumstellar dust using data from NASA's WISE mission. Initial cuts on the AllWISE catalog provide an input catalog of 277,686 sources. Volunteers then view images of each source online in 10 different bands to identify false-positives (galaxies, background stars, interstellar matter, image artifacts, etc.). Sources that survive this online vetting are followed up with spectroscopy on the FLWO Tillinghast telescope. This approach should allow us to unleash the full potential of WISE for finding new debris disks and protoplanetary disks. We announce a first list of 37 new disk candidates discovered by the project, and we describe our vetting and follow-up process. One of these systems appears to contain the first debris disk discovered around a star with a white dwarf companion: HD 74389. We also report four newly discovered classical Be stars (HD 6612, HD 7406, HD 164137, and HD 218546) and a new detection of 22 micron excess around a previously known debris disk host star, HD 22128.

We report results of investigation of amplitude calibration for very long baseline interferometry (VLBI) observations with Korean VLBI Network (KVN). Amplitude correction factors are estimated based on comparison of KVN observations at 22~GHz correlated by Daejeon hardware correlator and DiFX software correlator in Korea Astronomy and Space Science Institute (KASI) with Very Long Baseline Array (VLBA) observations at 22~GHz by DiFX software correlator in National Radio Astronomy Observatory (NRAO). We used the observations for compact radio sources, 3C 454.3, NRAO 512, OJ 287, BL Lac, 3C 279, 1633+382, and 1510-089, which are almost unresolved for baselines in a range of 350-477 km. Visibility data of the sources obtained with similar baselines at KVN and VLBA are selected, fringe-fitted, calibrated, and compared for their amplitudes. We found that visibility amplitudes of KVN observations should be corrected by factors of 1.10 and 1.35 when correlated by DiFX and Daejeon correlators, respectively. These correction factors are attributed to the combination of two steps of 2-bit quantization in KVN observing systems and characteristics of Daejeon correlator.

We report results of the performance evaluation of a new hardware correlator in Korea, the Daejeon correlator, developed by the Korea Astronomy and Space Science Institute (KASI) and the National Astronomical Observatory of Japan (NAOJ). We conducted Very Long Baseline Interferometry (VLBI) observations at 22~GHz with the Korean VLBI Network (KVN) in Korea and the VLBI Exploration of Radio Astrometry (VERA) in Japan, and correlated the aquired data with the Daejeon correlator. For evaluating the performance of the new hardware correlator, we compared the correlation outputs from the Daejeon correlator for KVN observations with those from a software correlator, the Distributed FX (DiFX). We investigated the correlated flux densities and brightness distributions of extragalactic compact radio sources. The comparison of the two correlator outputs show that they are consistent with each other within $<8\%$, which is comparable with the amplitude calibration uncertainties of KVN observations at 22~GHz. We also found that the 8\% difference in flux density is caused mainly by (a) the difference in the way of fringe phase tracking between the DiFX software correlator and the Daejeon hardware correlator, and (b) an unusual pattern (a double-layer pattern) of the amplitude correlation output from the Daejeon correlator. The visibility amplitude loss by the double-layer pattern is as small as 3\%. We conclude that the new hardware correlator produces reasonable correlation outputs for continuum observations, which are consistent with the outputs from the DiFX software correlator.

We investigate the spatial density configuration of stars around four metal-poor globular clusters (NGC 6266, NGC 6626, NGC 6642 and NGC 6723) in the Galactic bulge region using wide-field deep J, H, and K imaging data obtained with the WFCAM near-infrared array on United Kingdom Infrared Telescope. Statistical weighted filtering algorithm for the stars on the color-magnitude diagram is applied in order to sort cluster member candidates from the field star contamination. In two-dimensional isodensity contour maps of the clusters, we find that all of the four globular clusters exhibit strong evidence of tidally stripping stellar features beyond tidal radius, in the form of tidal tail or small density lobes or chunk. The orientations of the extended stellar substructures are likely to be associated with the effect of the dynamic interaction with the Galaxy and the cluster space motion. The observed radial density profiles of the four globular clusters also describe the extended substructures; they depart from theoretical King and Wilson models and have an overdensity feature with a break in a slope of profile at the outer region of clusters. The observed results could imply that four globular clusters in the Galactic bulge region have experienced strong environmental effect such as tidal force or bulge or disk shock of the Galaxy in the dynamical evolution of the globular clusters. These observational results provide us further constraints to understand the evolution of clusters in the Galactic bulge region as well as the formation of the Galaxy.

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.