A two dimensional hydrochemical hybrid code, KM2, is constructed to deal with astrophysical problems that would require coupled hydrodynamical and chemical evolution. The code assumes axisymmetry in cylindrical coordinate system, and consists of two modules: a hydrodynamics module and a chemistry module. The hydrodynamics module solves hydrodynamics using a Godunov-type finite volume scheme and treats included chemical species as passively advected scalars. The chemistry module implicitly solves non-equilibrium chemistry and change of the energy due to thermal processes with transfer of external ultraviolet radiation. Self-shielding effects on photodissociation of CO and H$_2$ are included. In this introductory paper, the adopted numerical method is presented, along with code verifications using the hydrodynamics modules, and a benchmark on the chemistry module with reactions specific to a photon-dominated region (PDR). Finally, as an example of the expected capability, the hydrochemical evolution of a PDR is presented based on the PDR benchmark.

We present an effective system for simulating dust continuum observations by radiative transfer simulations. By using workflow management system RENKEI-WFT, we utilized distributed computing resources and automated a sequence of computational tasks required for radiative transfer modeling, namely, main radiative transfer simulations, pre-/post-processes, and data transfer between computing resources. Our system simultaneously executes a lot of radiative transfer simulations with different input parameters on distributed computing resources. This capability of our system enables us to conduct effective research by radiative transfer simulation. As a demonstration of our system, we simulated dust continuum observations of protoplanetary disk. We performed hydrodynamic simulation modeling photoevaporating protoplanetary disk irradiated by ultra violet radiation from nearby massive stars. Results of this hydrodynamic simulation were used as input data for radiative transfer simulations. Expected spectral energy distributions and intensity maps were obtained by our system.

The bright-rimmed cloud SFO 22 was observed with the 45 m telescope of Nobeyama Radio Observatory in the ^{12}CO (J = 1-0), ^{13}CO (J = 1-0), and C^{18}O (J = 1-0) lines, where well-developed head-tail structure and small line widths were found. Such features were predicted by radiation-driven implosion models, suggesting that SFO 22 may be in a quasi-stationary equilibrium state. We compare the observed properties with those from numerical models of a photo-evaporating cloud, which include effects of magnetic pressure and heating due to strong far-ultraviolet (FUV) radiation from an exciting star. The magnetic pressure may play a more important role in the density structures of bright-rimmed clouds, than the thermal pressure that is enhanced by the FUV radiation. The FUV radiation can heat the cloud surface to near 30 K, however, its effect is not enough to reproduce the observed density structure of SFO 22. An initial magnetic field of 5 \mu G in our numerical models produces the best agreement with the observations, and its direction can affect the structures of bright-rimmed clouds.

Context. Molecular clouds near the H II regions tend to harbor more luminous protostars. Aims. Our aim in this paper is to investigate whether or not radiation-driven implosion mechanism enhances luminosity of protostars near regions of high-ionizing fluxes. Methods. We performed numerical simulations to model collapse of cores exposed to UV radiation from O stars. We investigated dependence of mass loss rates on the initial density profiles of cores and variation of UV fluxes. We derived simple analytic estimates of accretion rates and final masses of protostars. Results. Radiation-driven implosion mechanism can increase accretion rates of protostars by 1-2 orders of magnitude. On the other hand, mass loss due to photo-evaporation is not large enough to have a significant impact on the luminosity. The increase of accretion rate makes luminosity 1-2 orders higher than those of protostars that form without external triggering. Conclusions. Radiation-driven implosion can help explain the observed higher luminosity of protostars in molecular clouds near H II regions.

Recent observations indicate that some class 0 sources have orders of magnitude higher accretion rates than those of class I. We investigated the conditions for the high accretion rates of some class 0 sources by numerical calculations, modelling an external trigger. For no external trigger, we find that the maximum value of the accretion rate is determined by the ratio $\alpha$ of the gravitational energy to the thermal one within a flat inner region of the cloud core. The accretion rate reaches $\sim 10^{-4} M_{\sun} yr^{-1}$ if the cloud core has $ \alpha > 2$. For an external trigger we find that the maximum value of the accretion rate is proportional to the momentum given to the cloud core. The accretion rate reaches $ > 10^{-4} M_{\sun} yr^{-1}$ with a momentum of $\sim 0.1 M_{\sun} km s^{-1}$ when the initial central density of the cloud core is $\sim 10^{-18} g cm^{-3}$. A comparison between recent observational results for prestellar cores and our no triggered collapse model indicates that the flat inner regions of typical prestellar cores are not large enough to cause accretion rates of $\sim 10^{-4} M_{\sun} yr^{-1}$. Our results show that the triggered collapse of the cloud core is more preferable for the origin of the high accretion rates of class 0 sources than no triggered collapse.