Offsets of molecular line emission peaks from continuum peaks are very common but frequently difficult to explain with a single spherical cloud chemical model. We propose that the spatial projection effects of an irregular three dimensional (3D) cloud structure can be a solution. This work shows that the idea can be successfully applied to the Planck cold clump G224.4-0.6 by approximating it with four individual spherically symmetric cloud cores whose chemical patterns overlap with each other to produce observable line maps. With the empirical physical structures inferred from the observation data of this clump and a gas-grain chemical model, the four cores can satisfactorily reproduce its 850 $\mu$m continuum map and the diverse peak offsets of CCS, HC$_3$N and N$_2$H$^+$ simultaneously at chemical ages of about $8\times 10^5\sim 3\times 10^6$ yrs. The 3D projection effects on chemistry has the potential to explain such asymmetrical distributions of chemicals in many other molecular clouds.

Abundant carbon-chain species have been observed towards lukewarm corinos L1527, B228, and L483. These carbon-chain species are believed to be synthesized in the gas phase after CH$_4$ desorbs from the dust grain surface at the temperature around 30 K. We investigate carbon-chain species formation in lukewarm corinos using a more rigorous numerical method and advanced surface chemical models. We use the macroscopic Monte Carlo method in simulations. In addition to the two-phase model, the basic multiphase model and the new multiphase models are used for modeling surface chemistry on dust grains. All volatile species can sublime at their sublimation temperatures in the two-phase model while most volatile species are frozen in the ice mantle before water ice sublimes in the basic and the new multiphase models. The new multiphase models allow more volatile species to sublime at their sublimation temperatures than the basic multiphase model does. When T $\sim$ 30 K, the abundances of gaseous CH$_4$ and CO in the two-phase model are the highest while the basic multiphase model predicts the lowest CO and CH$_4$ abundances among all models. The abundances of carbon-chain species in the basic and the new multiphase models are lower than that in the two-phase model when T $\sim$ 30 K because CH$_4$ is crucial for the synthesis of carbon-chain species. The two-phase model performs the best to predict carbon-chain species abundances to fit observations while the basic multiphase model works the worst. The abundances of carbon-chain species predicted by the new multiphase models agree reasonably well with observations. The amount of CH$_4$ that can diffuse inside the ice mantle, thus sublime upon warm-up plays a crucial role in the synthesis of carbon-chain species in the gas phase. The carbon-chain species observed in lukewarm corinos may be able to gauge surface chemical models.

We investigate the chemical evolution of a collapsing core that starts from a hydrostatic core and finally form a low-mass protostar. New multiphase gas-grain models that include bulk diffusion and photon penetration are simulated by the macroscopic Monte Carlo method in order to derive the chemical evolution. There are two types of species in the ice bulk in the new multiphase models. Interstitial species can diffuse and sublime at their own sublimation temperatures while normal species are locked in the ice bulk. Photodissociation rates of icy species are reduced by the exponential decay of UV flux within the ice mantle. Two-phase models and basic multiphase models without bulk diffusion and photon penetration are also simulated for comparison. Our physical model for the collapsing core is base on a one-dimensional radiation hydrodynamics model. Abundant icy radicals are produced at around 10 K in the new multiphase models. Interstitial radicals can diffuse inside ice mantles to form complex organic molecules (COMS) upon warming-up. Thus, COMs produced by radical recombination at higher temperatures in the new multiphase models are more than one order of magnitude higher than those in the two-phase and basic multiphase models. Moreover, COMs produced at around 10 K in the new multiphase models are about one order of magnitude higher than those in the two-phase model. Our model shows a reasonable agreement with observations toward low-mass protostars. Moreover, molecular oxygen abundances predicted by our new multiphase models agree reasonably well with that found in cometary materials.

The majority of hydrogen in the interstellar medium (ISM) is in atomic form. The transition from atoms to molecules and, in particular, the formation of the H$_2$ molecule, is a key step in cosmic structure formation en route to stars. Quantifying H$_2$ formation in space is difficult, due to the confusion in the emission of atomic hydrogen (HI) and the lack of a H$_2$ signal from the cold ISM. Here we present the discovery of a rare, isolated dark cloud currently undergoing H$_2$ formation, as evidenced by a prominent "ring" of HI self-absorption. Through a combined analysis of HI narrow self-absorption, CO emission, dust emission, and extinction, we directly measured, for the first time, the [HI]/[H$_2$] abundance varying from 2% to 0.2%, within one region. These measured HI abundances are orders of magnitude higher than usually assumed initial conditions for protoplanetary disk models. None of the fast cloud formation model could produce such low atomic hydrogen abundance. We derived a cloud formation timescale of 6$\times$10$^6$ years, consistent with the global Galactic star formation rate, and favoring the classical star formation picture over fast star formation models. Our measurements also help constrain the H$_2$ formation rate, under various ISM conditions.

The temperatures of dust grains play important roles in the chemical evolution of molecular clouds. Unlike large grains, the temperature fluctuations of small grains induced by photons may be significant. Therefore, if the grain size distribution is included in astrochemical models, the temperatures of small dust grains may not be assumed to be constant. We simulate a full gas-grain reaction network with a set of dust grain radii using the classical MRN grain size distribution and include the temperature fluctuations of small dust grains. Monte Carlo method is used to simulate the real-time dust grain's temperature fluctuations which is caused by the external low energy photons and the internal cosmic ray induced secondary photons. The increase of dust grains radii as ice mantles accumulate on grain surfaces is also included in our models. We found that surface CO$_2$ abundances in models with grain size distribution and temperature fluctuations are more than one order of magnitude larger than those with single grain size. Small amounts of terrestrial complex organic molecules (COMs) can also be formed on small grains due to the temperature spikes induced by external low energy photons. However, cosmic ray induced secondary photons overheat small grains so that surface CO sublime and less radicals are formed on grains surfaces, thus the production of surface CO$_2$ and COMs decreases by about one order of magnitude. The overheating of small grains can be offset by grain growth so that the formation of surface CO$_2$ and COMs becomes more efficient.

We present 2.5-square-degree C$_{2}$H N=1-0 and N$_2$H$^+$ J=1-0 maps of the $\rho$ Ophiuchi molecular cloud complex. These are the first large-scale maps of the $\rho$ Ophiuchi molecular cloud complex with these two tracers. The C$_{2}$H emission is spatially more extended than the N$_2$H$^+$ emission. One faint N$_2$H$^+$ clump Oph-M and one C$_{2}$H ring Oph-RingSW are identified for the first time. The observed C$_{2}$H to N$_{2}$H$^{+}$ abundance ratio ([C$_{2}$H]/[N$_{2}$H$^{+}$]) varies between 5 and 110. We modeled the C$_{2}$H and N$_2$H$^+$ abundances with 1-D chemical models which show a clear decline of [C$_2$H]/[N$_2$H$^+$] with chemical age. Such an evolutionary trend is little affected by temperatures when they are below 40 K. At high density (n$_H$ $>$ 10$^5$ cm$^{-3}$), however, the time it takes for the abundance ratio to drop at least one order of magnitude becomes less than the dynamical time (e.g., turbulence crossing time $\rm \sim$10$^5$ years). The observed [C$_2$H]/[N$_2$H$^+$] difference between L1688 and L1689 can be explained by L1688 having chemically younger gas in relatively less dense regions. The observed [C$_{2}$H]/[N$_{2}$H$^{+}$] values are the results of time evolution, accelerated at higher densities. For the relative low density regions in L1688 where only C$_2$H emission was detected, the gas should be chemically younger.

The recent discovery of methyl formate and dimethyl ether in the gas phase of cold cores with temperatures as cold as 10 K challenges our previous astrochemical models concerning the formation of complex organic molecules. The strong correlation between the abundances and distributions of methyl formate and dimethyl ether further shows that current astrochemical models may be missing important chemical processes in cold astronomical sources. We investigate a scenario in which complex organic molecules and the methoxy radical can be formed on dust grains via a so-called "chain reaction" mechanism, in a similar manner to CO$_2$. A unified gas-grain microscopic-macroscopic Monte Carlo approach with both normal and interstitial sites for icy grain mantles is used to perform the chemical simulations. Reactive desorption with varying degrees of efficiency is included to enhance the non-thermal desorption of species formed on cold dust grains. In addition, varying degrees of efficiency for the surface formation of methoxy are also included. The observed abundances of a variety of organic molecules in cold cores can be reproduced in our models. The strong correlation between the abundances of methyl formate and dimethyl ether in cold cores can also be explained. Non-diffusive chemical reactions on dust grain surfaces may play a key role in the formation of some complex organic molecules.

Water is usually the main component of ice mantles, which cover the cores of dust grains in cold portions of dense interstellar clouds. When molecular hydrogen is adsorbed onto an icy mantle through physisorption, a common assumption in gas-grain rate equation models is to use an adsorption energy for molecular hydrogen on a pure water substrate. However, at high density and low temperature, when H2 is efficiently adsorbed onto the mantle, its surface abundance can be strongly overestimated if this assumption is still used. Unfortunately, the more detailed microscopic Monte Carlo treatment cannot be used to study the abundance of H2 in ice mantles if a full gas-grain network is utilized. We present a numerical method adapted for rate-equation models that takes into account the possibility that an H2 molecule can, while diffusing on the surface, find itself bound to another hydrogen molecule, with a far weaker bond than the H2-water bond, which can lead to more efficient desorption. We label the ensuing desorption "encounter desorption". The method is implemented first in a simple system consisting only of hydrogen molecules at steady state between gas and dust using the rate-equation approach and comparing the results with the results of a microscopic Monte Carlo calculation. We then discuss the use of the rate-equation approach with encounter desorption embedded in a complete gas-grain chemical network. For both systems, the rate-equation model with encounter desorption reproduces the H2 granular coverage computed by the microscopic Monte Carlo model. The method is especially useful for dense and cold environments, and for time-dependent physical conditions, such as occur in the collapse of dense cores and the formation of protoplanetary disks. It is not significantly CPU time consuming, so can be used for example with complex 3D chemical-hydrodynamical simulations.

We have designed an improved algorithm that enables us to simulate the chemistry of cold dense interstellar clouds with a full gas-grain reaction network. The chemistry is treated by a unified microscopic-macroscopic Monte Carlo approach that includes photon penetration and bulk diffusion. To determine the significance of these two processes, we simulate the chemistry with three different models. In Model 1, we use an exponential treatment to follow how photons penetrate and photodissociate ice species throughout the grain mantle. Moreover, the products of photodissociation are allowed to diffuse via bulk diffusion and react within the ice mantle. Model 2 is similar to Model 1 but with a slower bulk diffusion rate. A reference Model 0, which only allows photodissociation reactions to occur on the top two layers, is also simulated. Photodesorption is assumed to occur from the top two layers in all three models. We found that the abundances of major stable species in grain mantles do not differ much among these three models, and the results of our simulation for the abundances of these species agree well with observations. Likewise, the abundances of gas-phase species in the three models do not vary. However, the abundances of radicals in grain mantles can differ by up to two orders of magnitude depending upon the degree of photon penetration and the bulk diffusion of photodissociation products. We also found that complex molecules can be formed at temperatures as low as 10 K in all three models.

For the first time, we report a unified microscopic-macroscopic Monte Carlo simulation of gas-grain chemistry in cold interstellar clouds in which both the gas-phase and the grain surface chemistry are simulated by a stochastic technique. The surface chemistry is simulated with a microscopic Monte Carlo method in which the chemistry occurs on an initially flat surface. The surface chemical network consists of 29 reactions initiated by the accreting species H, O, C, and CO. Four different models are run with diverse but homogeneous physical conditions including temperature, gas density, and diffusion-barrier-to-desorption energy ratio. As time increases, icy interstellar mantles begin to grow. Our approach allows us to determine the morphology of the ice, layer by layer, as a function of time, and to ascertain the environment or environments for individual molecules. Our calculated abundances can be compared with observations of ices and gas-phase species, as well as the results of other models.

AIM: We have recently developed a microscopic Monte Carlo approach to study surface chemistry on interstellar grains and the morphology of ice mantles. The method is designed to eliminate the problems inherent in the rate-equation formalism to surface chemistry. Here we report the first use of this method in a chemical model of cold interstellar cloud cores that includes both gas-phase and surface chemistry. The surface chemical network consists of a small number of diffusive reactions that can produce molecular oxygen, water, carbon dioxide, formaldehyde, methanol and assorted radicals. METHOD: The simulation is started by running a gas-phase model including accretion onto grains but no surface chemistry or evaporation. The starting surface consists of either flat or rough olivine. We introduce the surface chemistry of the three species H, O and CO in an iterative manner using our stochastic technique. Under the conditions of the simulation, only atomic hydrogen can evaporate to a significant extent. Although it has little effect on other gas-phase species, the evaporation of atomic hydrogen changes its gas-phase abundance, which in turn changes the flux of atomic hydrogen onto grains. The effect on the surface chemistry is treated until convergence occurs. We neglect all non-thermal desorptive processes. RESULTS: We determine the mantle abundances of assorted molecules as a function of time through 2x10^5 yr. Our method also allows determination of the abundance of each molecule in specific monolayers. The mantle results can be compared with observations of water, carbon dioxide, carbon monoxide, and methanol ices in the sources W33A and Elias 16. Other than a slight underproduction of mantle CO, our results are in very good agreement with observations.