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.

Black hole accretion is one of nature's most efficient energy extraction processes. When gas falls in, a significant fraction of its gravitational binding energy is either converted into radiation or flows outwards in the form of black hole-driven jets and disk-driven winds. Recently, the Event Horizon Telescope (EHT), an Earth-size sub-millimetre radio interferometer, captured the first images of M87's black hole (or M87*). These images were analysed and interpreted using general-relativistic magnetohydrodynamics (GRMHD) models of accretion disks with rotation axes aligned with the black hole spin axis. However, since infalling gas is often insensitive to the black hole spin direction, misalignment between accretion disk and black hole spin may be a common occurrence in nature. In this work, we use the general-relativistic radiative transfer (GRRT) code \texttt{BHOSS} to calculate the first synthetic radio images of (highly) tilted disk/jet models generated by our GPU-accelerated GRMHD code \texttt{HAMR}. While the tilt does not have a noticeable effect on the system dynamics beyond a few tens of gravitational radii from the black hole, the warping of the disk and jet can imprint observable signatures in EHT images on smaller scales. Comparing the images from our GRMHD models to the 43 GHz and 230 GHz EHT images of M87, we find that M87 may feature a tilted disk/jet system. Further, tilted disks and jets display significant time variability in the 230 GHz flux that can be further tested by longer-duration EHT observations of M87.

General-relativistic magnetohydrodynamic (GRMHD) simulations have revolutionized our understanding of black-hole accretion. Here, we present a GPU-accelerated GRMHD code H-AMR with multi-faceted optimizations that, collectively, accelerate computation by 2-5 orders of magnitude for a wide range of applications. Firstly, it involves a novel implementation of a spherical-polar grid with 3D adaptive mesh refinement that operates in each of the 3 dimensions independently. This allows us to circumvent the Courant condition near the polar singularity, which otherwise cripples high-res computational performance. Secondly, we demonstrate that local adaptive time-stepping (LAT) on a logarithmic spherical-polar grid accelerates computation by a factor of $\lesssim10$ compared to traditional hierarchical time-stepping approaches. Jointly, these unique features lead to an effective speed of $\sim10^9$ zone-cycles-per-second-per-node on 5,400 NVIDIA V100 GPUs (i.e., 900 nodes of the OLCF Summit supercomputer). We demonstrate its computational performance by presenting the first GRMHD simulation of a tilted thin accretion disk threaded by a toroidal magnetic field around a rapidly spinning black hole. With an effective resolution of $13$,$440\times4$,$608\times8$,$092$ cells, and a total of $\lesssim22$ billion cells and $\sim0.65\times10^8$ timesteps, it is among the largest astrophysical simulations ever performed. We find that frame-dragging by the black hole tears up the disk into two independently precessing sub-disks. The innermost sub-disk rotation axis intermittently aligns with the black hole spin, demonstrating for the first time that such long-sought alignment is possible in the absence of large-scale poloidal magnetic fields.

Depending on the value of the accretion rate, black hole accretion is divided into cold and hot modes. The two modes have distinctly different physics and correspond to two feedback modes. Most previous feedback works either focus only on one mode, or the accretion physics is not always properly adopted. Here, by performing hydrodynamical numerical simulations of AGN feedback in an elliptical galaxy, we show that including both is important, and gives different results from including just one or the other. We specifically focus on the wind and radiation (but neglecting the jet) feedback in the hot mode, to explore whether this particular mode of feedback can play any important feedback role. For this aim, we have run two test models. In one model, we always adopt the cold mode no matter what value the accretion rate is; in another model, we turn off the AGN once it enters into the hot mode. We have calculated the AGN light curves, black hole growth, star formation, and AGN energy duty-cycle; and compared the results with the model in which the two modes are correctly included. Important differences are found. For example, if we were to adopt only the cold mode, the total mass of newly formed stars would become two orders of magnitude smaller, and the fraction of energy ejected within the high accretion regime (i.e., $L_{\rm BH} > 2\% L_{\rm Edd}$) would be too small to be consistent with observations. We have also investigated the respective roles of wind and radiation in the hot mode.

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 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.

Understanding how galaxies maintain the inefficiency of star formation with physically self-consistent models is a central problem for galaxy evolution. Although numerous theoretical models have been proposed in recent decades, the debate still exists. By means of high-resolution two-dimensional hydrodynamical simulations, we study the three feedback effects (the stellar wind heating, SNe feedback, and AGN feedback) in suppressing star formation activities on the evolution of early-type galaxies with different stellar masses. AGN feedback models are updated based on \citet{Yuan2018}. The gas sources comes exclusively from the mass losses of dying low-mass stars for most of our models. We find that SNe feedback can keep star formation at a significantly low level for low mass elliptical galaxies for a cosmological evolution time. For the high mass galaxies, AGN feedback can efficiently offset the radiative cooling and thus regulate the star formation activities. Such a suppression of star formation is extremely efficient in the inner region of the galaxies. AGB heating cannot account for this suppression for low and high mass galaxies. The X-ray temperature $T_{\rm X}$ and luminosity $L_{\rm X}$ of hot plasma can be in agreement with the observed data with the inclusion of effective feedback processes. These results thus suggest that we can use $T_{\rm X}$ and $L_{\rm X}$ to probe the role of different feedback processes. The inclusion of additional gas sources can make the mass scale between SNe and AGN feedback dominating in suppressing star formation decrease to an observationally inferred value of a few $10^{10}~M_{\odot}$.

This is the second paper of our series of works of studying the effects of active galactic nuclei (AGN) feedback on the cosmological evolution of an isolated elliptical galaxy by performing two-dimensional high-resolution hydrodynamical numerical simulations. In these simulations, the inner boundary is chosen so that the Bondi radius is resolved. Physical processes like star formation, SNe Ia and II are taken into account. Compared to previous works, the main improvements is that we adopt the most updated AGN physics, which is described in detail in the first paper of this series (Yuan et al. 2018, Paper I). These improvements include the discrimination of the two accretion modes of the central AGN and the most updated descriptions of the wind and radiation in the two modes. In Paper I, we consider the case that the specific angular momentum of the gas in the galaxy is very low. In this paper, we consider the case that the specific angular momentum of the gas is high. In the galactic scale, we adopt the gravitational torques raised due to non-axisymmetric structure in the galaxy as the mechanism of the transfer of angular momentum of gas, as proposed in some recent works. Since our simulations are axisymmetric, we make use of a parameterized prescription to mimic this mechanism. Same as Paper I, we investigate the AGN light curve, typical AGN lifetime, growth of the black hole mass, AGN duty-cycle, star formation, and the X-ray surface brightness of the galaxy. Special attention is paid to the effects of specific angular momentum of the galaxy on these properties. We find that some results are qualitatively similar to those shown in Paper I, while some results such as star formation and black hole growth do show a significant difference due to the mass concentration in the galactic disk as a consequence of galactic rotation.

This paper summarizes our recent works of studying AGN feedback in an isolated elliptical galaxy by performing high-resolution hydrodynamical numerical simulations. Bondi radius is resolved and the mass accretion rate of the black hole is calculated. The most updated AGN physics, namely the discrimination of cold and hot accretion modes and the exact descriptions of the AGN radiation and wind for a given accretion rate are adopted and their interaction with the gas in the host galaxy is calculated. Physical processes such as star formation and SNe feedback are taken into account. Consistent with observation, we find the AGN spends most of the time in the low-luminosity regime. AGN feedback overall suppresses the star formation; but depending on location in the galaxy and time, it can also enhance it. The light curve of specific star formation rate is not synchronous with the AGN light curve. We find that wind usually plays a dominant role in controlling the AGN luminosity and star formation, but radiation also cannot be neglected.

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 investigate the effects of AGN feedback on the cosmological evolution of an isolated elliptical galaxy by performing two-dimensional high-resolution hydrodynamical numerical simulations. The inner boundary of the simulation is chosen so that the Bondi radius is resolved. Compared to previous works, the two accretion modes, namely hot and cold, which correspond to different accretion rates and have different radiation and wind outputs, are carefully discriminated and the feedback effects by radiation and wind in each mode are taken into account. The most updated AGN physics, including the descriptions of radiation and wind from the hot accretion flows and wind from cold accretion disks, are adopted. Physical processes like star formation, Type Ia and Type II supernovae are taken into account. We study the AGN light curve, typical AGN lifetime, growth of the black hole mass, AGN duty-cycle, star formation, and the X-ray surface brightness of the galaxy. We compare our simulation results with observations and find general consistency. Comparisons with previous simulation works find significant differences, indicating the importance of AGN physics. The respective roles of radiation and wind feedbacks are examined and it is found that they are different for different problems of interest such as AGN luminosity and star formation. We find that it is hard to neglect any of them, so we suggest to use the names of "cold feedback mode" and "hot feedback mode" to replace the currently used ones.

We carry out a series of simulations of G2-type clouds interacting with the black hole at the galactic center, to determine why no large changes in the luminosity of Sgr A* were seen, and to determine the nature of G2. We measure the accretion rate from the gas cloud onto Sgr A* for a range of simulation parameters, such as cloud structure, background structure, background density, grid resolution, and accretion radius. For a broad range of parameters, the amount of cloud material accreted is small relative to the amount of background material accreted. The total accretion rate is not significantly effected for at least 30 years after periapsis. We find that reproducing observations of G2 likely requires two components for the object: an extended, cold gas cloud responsible for the Br-gamma emission, and a compact core or dusty stellar object dominating the bolometric luminosity. In simulations, the bolometric and X-ray luminosity have a peak lasting from about 1 year before to 1 year after periapsis, a feature not detected in observations. Our simulated Br-gamma emission is largely consistent with observations leading up to periapsis, with a slight increase in luminosity and a large increase in the FWHM of the line velocity. All emission from a gaseous component of G2 should fade rapidly after periapsis and be undetectable after 1 year, due to shock heating and expansion of the cloud. Any remaining emission should be from the compact component of G2.

Bow-shock pulsar wind nebulae are a subset of pulsar wind nebulae that form when the pulsar has high velocity due to the natal kick during the supernova explosion. The interaction between the relativistic wind from the fast-moving pulsar and the interstellar medium produces a bow-shock and a trail, which are detectable in H$_{\alpha}$ emission. Among such bow-shock pulsar wind nebulae, the Guitar Nebula stands out for its peculiar morphology, which consists of a prominent bow-shock head and a series of bubbles further behind. We present a scenario in which multiple bubbles can be produced when the pulsar encounters a series of density discontinuities in the ISM. We tested the scenario using 2-D and 3-D hydrodynamic simulations. The shape of the guitar nebula can be reproduced if the pulsar traversed a region of declining low density. We also show that if a pulsar encounters an inclined density discontinuity, it produces an asymmetric bow-shock head, consistent with observations of the bow-shock of the millisecond pulsar J2124-3358.

We study conditions for formation of recollimation shocks in jets interacting with stellar winds in high-mass X-ray binaries. We show the existence of a critical jet power, dependent on the wind rate and velocity and the jet velocity, above which a recollimation shock is not formed. For jet powers below critical, we derive the location of the shock. We show that surface shocks may still exist above the critical power, but only occupy a small volume of the jet and do not significantly alter the jet opening angle. We test these prediction by 3-D numerical simulations, which confirm the existence and the value of the critical power. We apply our results to Cyg X-1 and Cyg X-3.

Jets powered by high-mass X-ray binaries must traverse the powerful wind of the companion star. We present the first global 3D simulations of jet-wind interaction in high-mass X-ray binaries. We show that the wind momentum flux intercepted by the jet can lead to significant bending of the jet and that jets propagating through a spherical wind will be bent to an asymptotic angle $\psi_{\infty}$. We derive simple expressions for $\psi_{\infty}$ as a function of jet power and wind thrust. For known wind parameters, measurements of $\psi_{\infty}$ can be used to constrain the jet power. In the case of Cygnus X-1, the lack of jet precession as a function of orbital phase observed by the VLBA can be used to put a lower limit on the jet power of $L_{\rm jet} \gtrsim 10^{36}\,{\rm ergs\,s^{-1}}$. We further discuss the case where the initial jet is inclined relative to the binary orbital axis. We also analyze the case of Cygnus X-3 and show that jet bending is likely negligible unless the jet is significantly less powerful or much wider than currently thought. Our numerical investigation is limited to isotropic stellar winds. We discuss the possible effect of wind clumping on jet-wind interaction, which are likely significant, but argue that our limits on jet power for Cygnus X-1 are likely unaffected by clumping unless the global wind mass loss rate is orders of magnitude below the commonly assumed range for Cyg X-1.

The inner regions of barred galaxies contain substructures such as off-axis shocks, nuclear rings, and nuclear spirals. These substructure may affect star formation, and control the activity of a central black hole (BH) by determining the mass inflow rate. We investigate the formation and properties of such substructures using high-resolution, grid-based hydrodynamic simulations. The gaseous medium is assumed to be infinitesimally-thin, isothermal, and non-self-gravitating. The stars and dark matter are represented by a static gravitational potential with four components: a stellar disk, the bulge, a central BH, and the bar. To investigate various galactic environments, we vary the gas sound speed c_s as well as the mass of the central BH M_BH. Once the flow has reached a quasi-steady state, off-axis shocks tend to move closer to the bar major axis as c_s increases. Nuclear rings shrink in size with increasing c_s, but are independent of M_BH, suggesting that ring position is not determined by the Lindblad resonances. Rings in low-c_s models are narrow since they are occupied largely by gas on x2-orbits and well decoupled from nuclear spirals, while they become broad because of large thermal perturbations in high-c_s models. Nuclear spirals persist only when either c_s is small or M_BH is large; they would otherwise be destroyed completely by the ring material on eccentric orbits. The shape and strength of nuclear spirals depend sensitively on c_s and M_BH such that they are leading if both c_s and M_BH are small, weak trailing if c_s is small and M_BH is large, and strong trailing if both c_s and M_BH are large. While the mass inflow rate toward the nucleus is quite small in low-c_s models because of the presence of a narrow nuclear ring, it becomes larger than 0.01 Msun/yr when c_s is large, providing a potential explanation of nuclear activity in Seyfert galaxies.

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.

A sub-set of microquasars exhibit high peculiar velocity with respect to the local standard of rest due to the kicks they receive when being born in supernovae. The interaction between the radio plasma released by microquasar jets from such high-velocity binaries with the ISM must lead to the production of trails and bow shocks similar to what is observed in narrow-angle tailed radio galaxies and pulsar wind nebulae. We present a set of numerical simulations of this interaction that illuminate the long term dynamical evolution and the observational properties of these microquasar bow shock nebulae and trails. We find that this interaction always produces a structure that consists of a bow shock, a trailing neck, and an expanding bubble. Using our simulations to model emission, we predict that the shock surrounding the bubble and the neck should be visible in H{\alpha} emission, the interior of the bubble should be visible in synchrotron radio emission, and only the bow shock is likely to be detectable in X-ray emission. We construct an analytic model for the evolution of the neck and bubble shape and compare this model with observations of X-ray binary SAX J1712.6-3739.