Sgr A*, the supermassive black hole (SMBH) in our Galaxy, is dormant today, but it should have gone through multiple gas-accretion episodes in the past to grow to its current mass of $4\times10^6\,M_\odot$. Each episode temporarily ignites the SMBH and turns the Galactic Center into an active galactic nucleus (AGN). Recently, we showed that the AGN could produce large amount of hard X-rays that can penetrate the dense interstellar medium in the Galactic plane. Here we further study the impact of the X-rays on the molecular chemistry in our Galaxy. We use a chemical reaction network to simulate the evolution of several molecular species including $\rm{H_2O}$, $\rm{CH_3OH}$, and $\rm{H_2CO}$, both in the gas phase and on the surface of dust grains. We find that the X-ray irradiation could significantly enhance the abundances of these species. The effect is the most significant in those young, high-density molecular clouds, and could be prominent at a Galactic distance of $8$ kpc or smaller. The imprint in the chemical abundance is visible even several million years after the AGN turns off. The potential impact on the origin and evolution of organic and prebiotic molecules in the Milky Way deserves further investigation.

We performed high-sensitivity mapping observations of several complex organic molecules around Sagittarius B2 with ARO 12m telescope at 3-mm wavelength. Based on their spatial distribution, molecules can be classified as either "extended" that detected not only in Sgr B2(N) and Sgr B2(M), or "compact" that only detected toward or near to Sgr B2(N) and Sgr B2(M). The "extended" molecules including glycolaldehyde (CH2OHCHO), methyl formate (CH3OCHO), formic acid (t-HCOOH), ethanol (C2H5OH) and methyl amine (CH3NH2), while the "compact" molecules including dimethyl ether (CH3OCH3), ethyl cyanide (C2H5CN), and amino acetonitrile (H2NCH2CN). These "compact" molecules are likely produced under strong UV radiation, while "extended" molecules are likely formed under low-temperature, via gas-phase or grain surface reactions. The spatial distribution of "warm" CH2OHCHO at 89 GHz differ from the spatial distribution of "cold" CH2OHCHO observed at 13 GHz. We found evidence for an overabundance of CH2OHCHO compared to that expected from the gas-phase model, which indicates that grain-surface reactions are necessary to explain the origin of CH2OHCHO in Sagittarius B2. Grain-surface reactions are also needed to explain the correlation between the abundances of "cold" CH2OHCHO and C2H5OH. These results demonstrate the importance of grain-surface chemistry in the production of complex organic molecules.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

We report first results from the CHinese Exoplanet Searching Program from Antarctica (CHESPA)---a wide-field high-resolution photometric survey for transiting exoplanets carried out using telescopes of the AST3 (Antarctic Survey Telescopes times 3) project. There are now three telescopes (AST3-I, AST3-II, and CSTAR-II) operating at Dome A---the highest point on the Antarctic Plateau---in a fully automatic and remote mode to exploit the superb observing conditions of the site, and its long and uninterrupted polar nights. The search for transiting exoplanets is one of the key projects for AST3. During the Austral winters of 2016 and 2017 we used the AST3-II telescope to survey a set of target fields near the southern ecliptic pole, falling within the continuous viewing zone of the TESS mission \citep{Ricker10}. The first data release of the 2016 data, including images, catalogs and lightcurves of 26578 bright stars ($7.5\le i \le15$) was presented in \citet{Zhang18}. The best precision, as measured by the RMS of the lightcurves at the optimum magnitude of the survey ($i=10$), is around 2\,mmag. We detect 222 objects with plausible transit signals from these data, 116 of which are plausible transiting exoplanet candidates according to their stellar properties as given by the TESS Input Catalog \citep{Stassun17}, Gaia DR2 \citep{Gaia18} and TESS-HERMES spectroscopy \citep{Sharma18}. With the first data release from TESS expected in late 2018, this candidate list will be a timely for improving the rejection of potential false positives.

Located at Dome A, the highest point of the Antarctic plateau, the Chinese Kunlun station is considered to be one of the best ground-based photometric sites because of its extremely cold, dry, and stable atmosphere(Saunders et al. 2009). A target can be monitored from there for over 40 days without diurnal interruption during a polar winter. This makes Kunlun station a perfect site to search for short-period transiting exoplanets. Since 2008, an observatory has been built at Kunlun station and three telescopes are working there. Using these telescopes, the AST3 project has been carried out over the last six years with a search for transiting exoplanets as one of its key programs (CHESPA). In the austral winters of 2016 and 2017, a set of target fields in the Southern CVZ of TESS (Ricker et al. 2009) were monitored by the AST3-II telescope. In this paper, we introduce the CHESPA and present the first data release containing photometry of 26,578 bright stars (m_i < 15). The best photometric precision at the optimum magnitude for the survey is around 2 mmag. To demonstrate the data quality, we also present a catalog of 221 variables with a brightness variation greater than 5 mmag from the 2016 data. Among these variables, 179 are newly identified periodic variables not listed in the AAVSO databasea), and 67 are listed in the Candidate Target List(Stassun et al. 2017). These variables will require careful attention to avoid false-positive signals when searching for transiting exoplanets. Dozens of new transiting exoplanet candidates will be also released in a subsequent paper(Zhang et al. 2018b).

We present an analysis of meteorological data from the second generation of the Kunlun Automated Weather Station (KLAWS-2G) at Dome A, Antarctica during 2015 and 2016. We find that a strong temperature inversion exists for all the elevations up to 14 m that KLAWS-2G can reach, and lasts for more than 10 hours for 50% or more of the time when temperature inversion occurs. The average wind speeds at 4 m elevation are 4.2 m/s and 3.8 m/s during 2015 and 2016, respectively. The strong temperature inversion and moderate wind speed lead to a shallow turbulent boundary layer height at Dome A. By analyzing the temperature and wind shear profiles, we note telescopes should be elevated by at least 8 m above the ice. We also find that the duration of temperature inversions, and the wind speed, vary considerably from year to year. Therefore, long-term and continuous data are still needed for the site survey at Dome A.

The three Antarctic Survey Telescopes (AST3) aim to carry out time domain imaging survey at Dome A, Antarctica. The first of the three telescopes (AST3-1) was successfully deployed on January 2012. AST3-1 is a 500\,mm aperture modified Schmidt telescope with a 680\,mm diameter primary mirror. AST3-1 is equipped with a SDSS $i$ filter and a 10k $\times$ 10k frame transfer CCD camera, reduced to 5k $\times$ 10k by electronic shuttering, resulting in a 4.3 deg$^2$ field-of-view. To verify the capability of AST3-1 for a variety of science goals, extensive commissioning was carried out between March and May 2012. The commissioning included a survey covering 2000 deg$^2$ as well as the entire Large and Small Magellanic Clouds. Frequent repeated images were made of the center of the Large Magellanic Cloud, a selected exoplanet transit field, and fields including some Wolf-Rayet stars. Here we present the data reduction and photometric measurements of the point sources observed by AST3-1. We have achieved a survey depth of 19.3\,mag in 60 s exposures with 5\,mmag precision in the light curves of bright stars. The facility achieves sub-mmag photometric precision under stable survey conditions, approaching its photon noise limit. These results demonstrate that AST3-1 at Dome A is extraordinarily competitive in time-domain astronomy, including both quick searches for faint transients and the detection of tiny transit signals.

Proxima Centauri is known as the closest star from the Sun. Recently, radial velocity observations revealed the existence of an Earth-mass planet around it. With an orbital period of ~11 days, the surface of Proxima Centauri b is temperate and might be habitable. We took a photometric monitoring campaign to search for its transit, using the Bright Star Survey Telescope at the Zhongshan Station in Antarctica. A transit-like signal appearing on 2016 September 8th, is identified tentatively. Its midtime, $T_{C}=2,457,640.1990\pm0.0017$ HJD, is consistent with the predicted ephemeris based on RV orbit in a 1$\sigma$ confidence interval. Time-correlated noise is pronounced in the light curve of Proxima Centauri, affecting detection of transits. We develop a technique, in a Gaussian process framework, to gauge the statistical significance of potential transit detection. The tentative transit signal reported here, has a confidence level of $2.5\sigma$. Further detection of its periodic signals is necessary to confirm the planetary transit of Proxima Centauri b. We plan to monitor Proxima Centauri in next Polar night at Dome A in Antarctica, taking the advantage of continuous darkness. \citet{Kipping17} reported two tentative transit-like signals of Proxima Centauri b, observed by the Microvariability and Oscillation of Stars space Telescope in 2014 and 2015, respectively. The midtransit time of our detection is 138 minutes later than that predicted by their transit ephemeris. If all the signals are real transits, the misalignment of the epochs plausibly suggests transit timing variations of Proxima Centauri b induced by an outer planet in this system.

The LIGO detection of gravitational waves (GW) from merging black holes in 2015 marked the beginning of a new era in observational astronomy. The detection of an electromagnetic signal from a GW source is the critical next step to explore in detail the physics involved. The Antarctic Survey Telescopes (AST3), located at Dome A, Antarctica, is uniquely situated for rapid response time-domain astronomy with its continuous night-time coverage during the austral winter. We report optical observations of the GW source (GW~170817) in the nearby galaxy NGC 4993 using AST3. The data show a rapidly fading transient at around 1 day after the GW trigger, with the $i$-band magnitude declining from $17.23\pm0.13$ magnitude to $17.72\pm0.09$ magnitude in $\sim 1.8$ hour. The brightness and time evolution of the optical transient associated with GW~170817 are broadly consistent with the predictions of models involving merging binary neutron stars. We infer from our data that the merging process ejected about $\sim 10^{-2}$ solar mass of radioactive material at a speed of up to $30\%$ the speed of light.

We performed very deep searches for 2 ground-state water transitions in 13 protoplanetary disks with the HIFI instrument on-board the Herschel Space Observatory, with integration times up to 12 hours per line. Two other water transitions that sample warmer gas were also searched for with shallower integrations. The detection rate is low, and the upper limits provided by the observations are generally much lower than predictions of thermo-chemical models with canonical inputs. One ground-state transition is newly detected in the stacked spectrum of AA Tau, DM Tau, LkCa 15, and MWC 480. We run a grid of models to show that the abundance of gas-phase oxygen needs to be reduced by a factor of at least ~100 to be consistent with the observational upper limits (and positive detections) if a dust-to-gas mass ratio of 0.01 were to be assumed. As a continuation of previous ideas, we propose that the underlying reason for the depletion of oxygen (hence the low detection rate) is the freeze-out of volatiles such as water and CO onto dust grains followed by grain growth and settling/migration, which permanently removes these gas-phase molecules from the emissive upper layers of the outer disk. Such depletion of volatiles is likely ubiquitous among different disks, though not necessarily to the same degree. The volatiles might be returned back to the gas phase in the inner disk (within about 15 AU), which is consistent with current constraints. Comparison with studies on disk dispersal due to photoevaporation indicates that the timescale for volatile depletion is shorter than that of photoevaporation.

AST3-1 is the second-generation wide-field optical photometric telescope dedicated to time domain astronomy at Dome A, Antarctica. Here we present the results of $i$ band images survey from AST3-1 towards one Galactic disk field. Based on time-series photometry of 92,583 stars, 560 variable stars were detected with $i$ magnitude $\leq$ 16.5 mag during eight days of observations; 339 of these are previously unknown variables. We tentatively classify the 560 variables as 285 eclipsing binaries (EW, EB, EA), 27 pulsating variable stars ($\delta$~Scuti, $\gamma$~Doradus, $\delta$~Cephei variable and RR Lyrae stars) and 248 other types of variables (unclassified periodic, multi-periodic and aperiodic variable stars). Among the eclipsing binaries, 34 show O'Connell effects. One of the aperiodic variables shows a plateau light curve and another one shows a secondary maximum after peak brightness. We also detected a complex binary system with RS CVn-like light curve morphology; this object is being followed-up spectroscopically using the Gemini South telescope.

We report observations of resolved C2H emission rings within the gas-rich protoplanetary disks of TW Hya and DM Tau using the Atacama Large Millimeter Array (ALMA). In each case the emission ring is found to arise at the edge of the observable disk of mm-sized grains (pebbles) traced by (sub)mm-wave continuum emission. In addition, we detect a C3H2 emission ring with an identical spatial distribution to C2H in the TW Hya disk. This suggests that these are hydrocarbon rings (i.e. not limited to C2H). Using a detailed thermo-chemical model we show that reproducing the emission from C2H requires a strong UV field and C/O > 1 in the upper disk atmosphere and outer disk, beyond the edge of the pebble disk. This naturally arises in a disk where the ice-coated dust mass is spatially stratified due to the combined effects of coagulation, gravitational settling and drift. This stratification causes the disk surface and outer disk to have a greater permeability to UV photons. Furthermore the concentration of ices that transport key volatile carriers of oxygen and carbon in the midplane, along with photochemical erosion of CO, leads to an elemental C/O ratio that exceeds unity in the UV-dominated disk. Thus the motions of the grains, and not the gas, lead to a rich hydrocarbon chemistry in disk surface layers and in the outer disk midplane.

Deuterium-to-hydrogen (D/H) enrichments in molecular species provide clues about their original formation environment. The organic materials in primitive solar system bodies have generally higher D/H ratios and show greater D/H variation when compared to D/H in solar system water. We propose this difference arises at least in part due to 1) the availability of additional chemical fractionation pathways for organics beyond that for water, and 2) the higher volatility of key carbon reservoirs compared to oxygen. We test this hypothesis using detailed disk models, including a sophisticated, new disk ionization treatment with a low cosmic ray ionization rate, and find that disk chemistry leads to higher deuterium enrichment in organics compared to water, helped especially by fractionation via the precursors CH$_2$D$^+$/CH$_3^+$. We also find that the D/H ratio in individual species varies significantly depending on their particular formation pathways. For example, from $\sim20-40$ AU, CH$_4$ can reach $\rm{D/H\sim2\times10^{-3}}$, while D/H in CH$_3$OH remains locally unaltered. Finally, while the global organic D/H in our models can reproduce intermediately elevated D/H in the bulk hydrocarbon reservoir, our models are unable to reproduce the most deuterium-enriched organic materials in the solar system, and thus our model requires some inheritance from the cold interstellar medium from which the Sun formed.

This paper discusses compelling science cases for a future long-baseline interferometer operating at millimeter and centimeter wavelengths, like the proposed Next Generation Vary Large Array (ngVLA). We report on the activities of the Cradle of Life science working group, which focused on the formation of low- and high-mass stars, the formation of planets and evolution of protoplanetary disks, the physical and compositional study of Solar System bodies, and the possible detection of radio signals from extraterrestrial civilizations. We propose 19 scientific projects based on the current specification of the ngVLA. Five of them are highlighted as possible Key Science Projects: (1) Resolving the density structure and dynamics of the youngest HII regions and high-mass protostellar jets, (2) Unveiling binary/multiple protostars at higher resolution, (3) Mapping planet formation regions in nearby disks on scales down to 1 AU, (4) Studying the formation of complex molecules, and (5) Deep atmospheric mapping of giant planets in the Solar System. For each of these projects, we discuss the scientific importance and feasibility. The results presented here should be considered as the beginning of a more in-depth analysis of the science enabled by such a facility, and are by no means complete or exhaustive.

An abundance decrease in carbon- and oxygen-bearing species relative to dust has been frequently found in planet-forming disks, which can be attributed to an overall reduction of gas mass. However, in the case of TW Hya, the only disk with gas mass measured directly with HD rotational lines, the inferred gas mass ($\lesssim$0.005 solar mass) is significantly below the directly measured value ($\gtrsim$0.05 solar mass). We show that this apparent conflict can be resolved if the elemental abundances of carbon and oxygen are reduced in the upper layers of the outer disk but are normal elsewhere (except for a possible enhancement of their abundances in the inner disk). The implication is that in the outer disk, the main reservoir of the volatiles (CO, water, ...) resides close to the midplane, locked up inside solid bodies that are too heavy to be transported back to the atmosphere by turbulence. An enhancement in the carbon and oxygen abundances in the inner disk can be caused by inward migration of these solid bodies. This is consistent with estimates based on previous models of dust grain dynamics. Indirect measurements of the disk gas mass and disk structure from species such as CO will thus be intertwined with the evolution of dust grains, and possibly also with the formation of planetesimals.

Identifying the source of Earth's water is central to understanding the origins of life-fostering environments and to assessing the prevalence of such environments in space. Water throughout the solar system exhibits deuterium-to-hydrogen enrichments, a fossil relic of low-temperature, ion-derived chemistry within either (i) the parent molecular cloud or (ii) the solar nebula protoplanetary disk. Utilizing a comprehensive treatment of disk ionization, we find that ion-driven deuterium pathways are inefficient, curtailing the disk's deuterated water formation and its viability as the sole source for the solar system's water. This finding implies that if the solar system's formation was typical, abundant interstellar ices are available to all nascent planetary systems.

Water vapor has been detected in protoplanetary disks. In this work we model the distribution of water vapor in protoplanetary disks with a thermo-chemical code. For a set of parameterized disk models, we calculate the distribution of dust temperature and radiation field of the disk with a Monte Carlo method, and then solve the gas temperature distribution and chemical composition. The radiative transfer includes detailed treatment of scattering by atomic hydrogen and absorption by water of Lyman alpha photons, since the Lyman alpha line dominates the UV spectrum of accreting young stars. In a fiducial model, we find that warm water vapor with temperature around 300 K is mainly distributed in a small and well-confined region in the inner disk. The inner boundary of the warm water region is where the shielding of UV field due to dust and water itself become significant. The outer boundary is where the dust temperature drops below the water condensation temperature. A more luminous central star leads to a more extended distribution of warm water vapor, while dust growth and settling tends to reduce the amount of warm water vapor. Based on typical assumptions regarding the elemental oxygen abundance and the water chemistry, the column density of warm water vapor can be as high as 1E22 cm^{-2}. A small amount of hot water vapor with temperature higher than ~300 K exists in a more extended region in the upper atmosphere of the disk. Cold water vapor with temperature lower than 100 K is distributed over the entire disk, produced by photodesorption of the water ice.

Context: Hydrogen peroxide (HOOH) was recently detected toward \rho Oph A. Subsequent astrochemical modeling that included reactions in the gas phase and on the surface of dust grains was able to explain the observed abundance, and highlighted the importance of grain chemistry in the formation of HOOH as an intermediate product in water formation. This study also predicted that the hydroperoxyl radical HO2, the precursor of HOOH, should be detectable. Aims: We aim at detecting the hydroperoxyl radical HO2 in \rho Oph A. Methods: We used the IRAM 30m and the APEX telescopes to target the brightest HO2 lines at about 130 and 260 GHz. Results: We detect five lines of HO2 (comprising seven individual molecular transitions). The fractional abundance of HO2 is found to be about 1e-10, a value similar to the abundance of HOOH. This observational result is consistent with the prediction of the above mentioned astrochemical model, and thereby validates our current understanding of the water formation on dust grains. Conclusions: This detection, anticipated by a sophisticated gas-grain chemical model, demonstrates that models of grain chemistry have improved tremendously and that grain surface reactions now form a crucial part of the overall astrochemical network.

Although water is an essential and widespread molecule in star-forming regions, its chemical formation pathways are still not very well constrained. Observing the level of deuterium fractionation of OH, a radical involved in the water chemical network, is a promising way to infer its chemical origin. We aim at understanding the formation mechanisms of water by investigating the origin of its deuterium fractionation. This can be achieved by observing the abundance of OD towards the low-mass protostar IRAS16293-2422, where the HDO distribution is already known. Using the GREAT receiver on board SOFIA, we observed the ground-state OD transition at 1391.5 GHz towards the low-mass protostar IRAS16293-2422. We also present the detection of the HDO 111-000 line using the APEX telescope. We compare the OD/HDO abundance ratio inferred from these observations with the predictions of chemical models. The OD line is detected in absorption towards the source continuum. This is the first detection of OD outside the solar system. The SOFIA observation, coupled to the observation of the HDO 111-000 line, provides an estimate of the abundance ratio OD/HDO ~ 17-90 in the gas where the absorption takes place. This value is fairly high compared with model predictions. This may be reconciled if reprocessing in the gas by means of the dissociative recombination of H2DO+ further fractionates OH with respect to water. The present observation demonstrates the capability of the SOFIA/GREAT instrument to detect the ground transition of OD towards star-forming regions in a frequency range that was not accessible before. Dissociative recombination of H2DO+ may play an important role in setting a high OD abundance. Measuring the branching ratios of this reaction in the laboratory will be of great value for chemical models.

Context. The formation of water on the dust grains in the interstellar medium may proceed with hydrogen peroxide (H2O2) as an intermediate. Recently gas-phase H2O2 has been detected in {\rho} Oph A with an abundance of ~1E-10 relative to H2. Aims. We aim to reproduce the observed abundance of H2O2 and other species detected in {\rho} Oph A quantitatively. Methods. We make use of a chemical network which includes gas phase reactions as well as processes on the grains; desorption from the grain surface through chemical reaction is also included. We run the model for a range of physical parameters. Results. The abundance of H2O2 can be best reproduced at ~6E5 yr, which is close to the dynamical age of {\rho} Oph A. The abundances of other species such as H2CO, CH3OH, and O2 can be reasonably reproduced also at this time. In the early time the gas-phase abundance of H2O2 can be much higher than the current detected value. We predict a gas phase abundance of O2H at the same order of magnitude as H2O2, and an abundance of the order 1E-8 for gas phase water in {\rho} Oph A. A few other species of interest are also discussed. Conclusions. We demonstrate that H2O2 can be produced on the dust grains and released into the gas phase through non-thermal desorption via surface exothermic reactions. The H2O2 molecule on the grain is an important intermediate in the formation of water. The fact that H2O2 is over-produced in the gas phase for a range of physical conditions suggests that its destruction channel in the current gas phase network may be incomplete.

The structure in density and temperature of protoplanetary disks surrounding low-mass stars is not yet well known. The protoplanetary disks mid-planes are expected to be very cold and thus depleted in molecules in gas phase, especially CO. Recent observations of molecules at very low apparent temperature (~ 6 K) challenge this current picture of the protoplanetary disk structures. We aim at constraining the physical conditions, and in particular the gas-phase CO abundance in the mid-plane of protoplanetary disks. The light molecule H2D+,is a tracer of cold and CO-depleted environment. It is therefore a good candidate to explore the disks mid-planes. We performed a deep search for H2D+ in the two well-known disks surrounding TW Hya and DM Tau using the APEX and JCMT telescopes. The analysis of the observations are done with DISKFIT, a radiative transfer code dedicated to disks. In addition, we used a chemical model describing deuterium chemistry to infer the implications of our observations on the level of CO depletion and on the ionization rate in the disks mid-plane. The ortho-H2D+(1(1,0))-1(1,1)) line at 372 GHz was not detected. Although our limit is three times better than previous observations, comparison with the chemical modeling indicates that it remains insufficient to put valuable constraints on the CO abundance in the disk mid-plane. Even with ALMA, the detection of H2D+ may not be straightforward, and H2D+ may not be a sufficiently sensitive tracer of the protoplanetary disks mid-plane

[Context] The stochasticity of grain chemistry requires special care in modeling. Previously methods based on the modified rate equation, the master equation, the moment equation, and Monte Carlo simulations have been used. [Aims] We attempt to develop a systematic and efficient way to model the gas-grain chemistry with a large reaction network as accurately as possible. [Methods] We present a hybrid moment equation approach which is a general and automatic method where the generating function is used to generate the moment equations. For large reaction networks, the moment equation is cut off at the second order, and a switch scheme is used when the average population of certain species reaches 1. For small networks, the third order moments can also be utilized to achieve a higher accuracy. [Results] For physical conditions in which the surface reactions are important, our method provides a major improvement over the rate equation approach, when benchmarked against the rigorous Monte Carlo results. For either very low or very high temperatures, or large grain radii, results from the rate equation are similar to those from our new approach. Our method is faster than the Monte Carlo approach, but slower than the rate equation approach. [Conclusions] The hybrid moment equation approach with a cutoff and switch scheme is applicable to large gas-grain networks, and is accurate enough to be used for astrochemistry studies. The layered structure of the grain mantle could also be incorporated into this approach, although a full implementation of the grain micro-physics appears to be difficult.

Context: In the last years, the H2D+ and D2H+ molecules have gained great attention as probes of cold and depleted dense molecular cloud cores. These ions are at the basis of molecular deuterium fractionation, a common characteristic observed in star forming regions. H2D+ is now routinely observed, but the search for its isotopologue D2H+ is still difficult because of the high frequency of its ground para transition (692 GHz). Aims: We have observed molecular transitions of H2D+ and D2H+ in a cold prestellar core to characterize the roots of deuterium chemistry. Methods: Thanks to the sensitive multi-pixel CHAMP+ receiver on the APEX telescope where the required excellent weather conditions are met, we not only successfully detect D2H+ in the H-MM1 prestellar core located in the L1688 cloud, but also obtain information on the spatial extent of its emission. We also detect H2D+ at 372 GHz in the same source. We analyse these detections using a non-LTE radiative transfer code and a state-of-the-art spin-dependent chemical model. Results: This observation is the first secure detection of D2H+ in space. The emission is moreover extended over several pixels of the CHAMP+ array, i.e. on a scale of at least 40'', corresponding to ~ 4800 AU. We derive column densities on the order of 1e12-1e13 cm-2 for both molecules in the LTE approximation depending on the assumed temperature, and up to two orders of magnitude higher based on a non-LTE analysis. Conclusions: Our modeling suggests that the level of CO depletion must be extremely high (>10, and even >100 if the temperature of the core is around 10 K) at the core center, in contradiction with CO depletion levels directly measured in other cores. Observation of the H2D+ spatial distribution and direct measurement of the CO depletion in H-MM1 will be essential to confirm if present chemical models investigating the basis of deuterium [...].

InfraRed Dark Clouds (IRDCs) are extinction features against the Galactic infrared background, mainly in mid-infrared band. Recently they were proposed to be potential sites of massive star formation. In this work we have made 12CO, 13CO, and C18O (J=1-0) survey of 61 IRDCs, 52 of which are in the first Galactic quadrant, selected from a catalog given by Simon et al. (2006), while the others are in the outer Galaxy, selected by visually inspecting the MSX images. Detection rates in the three CO lines are 90%, 71%, and 62% respectively. The distribution of IRDCs in the first Galactic quadrant is consistent with the 5 kpc molecular ring picture, while slight trace of spiral pattern can also be noticed, which needs to be further examined. The IRDCs have typical excitation temperature of 10 K and typical column density of several 10^{22} cm^{-2}. Their typical physical size is estimated to be several pcs using angular sizes from the Simon catalog. Typical volume density and typical LTE mass are ~5000 cm^{-3} and ~5000 M_{sun} respectively. The IRDCs are in or near virial equilibrium. The properties of IRDCs are similar to those of star forming molecular clumps, and they seem to be intermediate between giant molecular clouds and Bok globules, thus they may represent early stages of massive star formation.