We have carried out mapping observations of molecular emission lines of HC$_{3}$N and CH$_{3}$OH toward two massive cluster-forming clumps, NGC2264-C and NGC2264-D, using the Nobeyama 45-m radio telescope. We derive an $I$(HC$_{3}$N)/$I$(CH$_{3}$OH) integrated intensity ratio map, showing a higher value at clumps including 2MASS point sources at the northern part of NGC2264-D. Possible interpretations of the $I$(HC$_{3}$N)/$I$(CH$_{3}$OH) ratio are discussed. We have also observed molecular emission lines from CCS and N$_{2}$H$^{+}$ toward five positions in each clump. We investigate the $N$(N$_{2}$H$^{+}$)/$N$(CCS) and $N$(N$_{2}$H$^{+}$)/$N$(HC$_{3}$N) column density ratios among the ten positions in order to test whether they can be used as chemical evolutionary indicators in these clumps. The $N$(N$_{2}$H$^{+}$)/$N$(CCS) ratio shows a very high value toward a bright embedded IR source (IRS1), whereas the $N$(N$_{2}$H$^{+}$)/$N$(HC$_{3}$N) ratio at IRS1 is comparable with those at the other positions. These results suggest that UV radiation affects the chemistry around IRS1. We find that there are positive correlations between these column density ratios and the excitation temperatures of N$_{2}$H$^{+}$, which implies the chemical evolution of clumps. These chemical evolutionary indicators likely reflect the combination of evolution along the filamentary structure and evolution of each clump.

In this work, we reexamine sulfur chemistry occurring on and in the ice mantles of interstellar dust grains, and report the effects of two new modifications to standard astrochemical models; namely, (a) the incorporation of cosmic ray-driven radiation chemistry and (b) the assumption of fast, non-diffusive reactions for key radicals in the bulk. Results from our models of dense molecular clouds show that these changes can have a profound influence on the abundances of sulfur-bearing species in ice mantles, including a reduction in the abundance of solid-phase H$_2$S and HS, and a significant increase in the abundances of OCS, SO$_2$, as well as pure allotropes of sulfur, especially S$_8$. These pure-sulfur species - though nearly impossible to observe directly - have long been speculated to be potential sulfur reservoirs and our results represent possibly the most accurate estimates yet of their abundances in the dense ISM. Moreover, the results of these updated models are found to be in good agreement with available observational data. Finally, we examine the implications of our findings with regard to the as-yet-unknown sulfur reservoir thought to exist in dense interstellar environments.

We have carried out observations of CCH and its two $^{13}$C isotopologues, $^{13}$CCH and C$^{13}$CH, in the 84 - 88 GHz band toward two starless cores, L1521B and L134N (L183), using the Nobeyama 45 m radio telescope. We have detected C$^{13}$CH with a signal-to-noise (S/N) ratio of 4, whereas no line of $^{13}$CCH was detected in either the dark clouds. The column densities of the normal species were derived to be ($1.66 \pm 0.18$)$\times 10^{14}$ cm$^{-2}$ and ($7.3 \pm 0.9$)$\times 10^{13}$ cm$^{-2}$ ($1 \sigma$) in L1521B and L134N, respectively. The column density ratios of $N$(C$^{13}$CH)/$N$($^{13}$CCH) were calculated to be $>1.1$ and $>1.4$ in L1521B and L134N, respectively. The characteristic that $^{13}$CCH is less abundant than C$^{13}$CH is likely common for dark clouds. Moreover, we find that the $^{12}$C/$^{13}$C ratios of CCH are much higher than those of HC$_{3}$N in L1521B by more than a factor of 2, as well as in Taurus Molecular Cloud-1 (TMC-1). In L134N, the differences in the $^{12}$C/$^{13}$C ratios between CCH and HC$_{3}$N seem to be smaller than those in L1521B and TMC-1. We discuss the origins of the $^{13}$C isotopic fractionation of CCH and investigate possible routes that cause the significantly high $^{12}$C/$^{13}$C ratio of CCH especially in young dark clouds, with the help of chemical simulations. The high $^{12}$C/$^{13}$C ratios of CCH seem to be caused by reactions between hydrocarbons (e.g., CCH, C$_{2}$H$_{2}$, $l,c$-C$_{3}$H) and C$^{+}$.

Vibrationally excited states of detected interstellar molecules have been shown to account for a large portion of unidentified spectral lines in observed interstellar spectra toward chemically rich sources. Here, we present the first interstellar detection of the first and second vibrationally excited torsional states of acetic acid ($v_\mathrm{t} = 1, 2$) toward the high-mass star-forming region NGC 6334I. The observations presented were taken with the Atacama Large Millimeter/submillimeter Array in bands 4, 6, and 7 covering a frequency range of 130 - 352 GHz. By comparing a single excitation temperature model to the observations, the best-fit excitation temperature and column density are obtained to be 142(25) K and $1.12(7) \times 10^{17} \mathrm{cm^{-2}}$ respectively. Based on the intensity maps of the vibrationally excited CH$_3$COOH transitions, we found that the CH$_3$COOH emissions are compact and concentrated toward the MM1 and MM2 regions with a source size smaller than 2 arcsec. After locating the emission from different CH$_3$COOH transitions, which cover a large range of excitation energies, we are able to explain the variation of the CH$_3$COOH emission peak within the MM2 core by invoking continuum absorption or outflows.

Recent radio astronomical observations have revealed that HC$_{5}$N, the second shortest cyanopolyyne (HC$_{2n+1}$N), is abundant around some massive young stellar objects (MYSOs), which is not predicted by classical carbon-chain chemistry. For example, the observed HC$_{5}$N abundance toward the G28.28$-$0.36 MYSO is higher than that in L1527, which is one of the warm carbon chain chemistry (WCCC) sources, by more than one order of magnitude (Taniguchi et al., 2017). In this paper, we present chemical simulations of hot-core models with a warm-up period using the astrochemical code Nautilus. We find that the cyanopolyynes are formed initially in the gas phase and accreted onto the bulk and surface of granular ice mantles during the lukewarm phase, which occurs at $25 < T < 100$ K. In slow warm-up period models, the peak abundances occur as the cyanopolyynes desorb from dust grains after the temperature rises above 100 K. The lower limits of the abundances of HC$_{5}$N, CH$_{3}$CCH, and CH$_{3}$OH observed in the G28.28$-$0.36 MYSO can be reproduced in our hot-core models, after their desorption from dust grains. Moreover, previous observations suggested chemical diversity in envelopes around different MYSOs. We discuss possible interpretations of relationships between stages of the star-formation process and such chemical diversity, such as the different warm-up timescales. This timescale depends not only on the mass of central stars but also on the relationship between the size of warm regions and their infall velocity.

Shocks are a crucial probe for understanding the ongoing chemistry within ices on interstellar dust grains where many complex organic molecules (COMs) are believed to be formed. However, previous work has been limited to the initial liberation into the gas phase through non-thermal desorption processes such as sputtering. Here, we present results from the adapted three-phase gas-grain chemical network code NAUTILUS, with the inclusion of additional high-temperature reactions, non-thermal desorption, collisional dust heating, and shock-physics parameters. This enhanced model is capable of reproducing many of the molecular distributions and abundance ratios seen in our prior observations of the prototypical shocked-outflow L1157. In addition, we find that, among others, NH$_2$CHO, HCOOCH$_3$, and CH$_3$CHO have significant post-shock chemistry formation routes that differ from those of many other COMs observed in shocks. Finally, a number of selected species and phenomena are studied here with respect to their usefulness as shock tracers in various astrophysical sources.

Many astrochemical models today explicitly consider the species that comprise the bulk of interstellar dust grain ice-mantles separately from those in the top few monolayers. Bombardment of these ices by ionizing radiation - whether in the form of cosmic rays, stellar winds, or radionuclide emission - represents an astrochemically viable means of driving a rich chemistry even in the bulk of the ice-mantle, now supported by a large body of work in laboratory astrophysics. In this study, using an existing rate equation-based astrochemical code modified to include a method of considering radiation chemistry recently developed by us, we attempted to simulate two such studies in which (a) pure O$_2$ ice at 5 K and, (b) pure H$_2$O ice at 16 K and 77 K, were bombarded by keV H$^+$ ions. Our aims are twofold: (1) to test the capability of our newly developed method to replicate the results of ice-irradiation experiments, and (2) to determine in such a well-constrained system how bulk chemistry is best handled using the same gas-grain codes that are used to model the interstellar medium (ISM). We find that our modified astrochemical model is able to reproduce both the abundance of O$_3$ in the 5 K pure O$_2$ ice, as well as both the abundance of H$_2$O$_2$ in the 16 K water ice and the previously noted decrease of hydrogen peroxide at higher temperatures. However, these results require the assumption that radicals and other reactive species produced via radiolysis react quickly and non-diffusively with neighbors in the ice.

Thioformaldehyde is an abundant molecule in various regions of the interstellar medium. However, available laboratory data limit the accuracies of calculated transition frequencies in the submillimeter region, in particular for minor isotopic species. We aim to determine spectroscopic parameters of isotopologs of H2CS that are accurate enough for predictions well into the submillimeter region. We investigated the laboratory rotational spectra of numerous isotopic species in natural isotopic composition almost continuously between 110 and 377 GHz. Individual lines were studied for most species in two frequency regions between 566 and 930 GHz. Further data were obtained for the three most abundant species in the 1290-1390 GHz region. New or improved spectroscopic parameters were determined for seven isotopic species. Quantum-chemical calculations were carried out to evaluate the differences between ground state and equilibrium rotational parameters to derive semi-empirical equilibrium structural parameters. The spectroscopic parameters are accurate enough for predictions well above 1 THz with the exception of H2(13)C(34)S where the predictions should be reliable to around 700 GHz.

The C$_2$H$_4$O$_2$ isomers have been previously investigated primarily via disparate sets of observations involving single dish and array measurements. The only attempt at using a uniform set of observations was performed with the IRAM 30 m observation in 2013 (Belloche et al. 2013). In this study, we present an intensive and rigorous spectral and morphological analysis of the C$_2$H$_4$O$_2$ isomers towards Sgr B2(N) with interferometers, ALMA Band 3 observations. We propose a quantitative selection method, which automates the determination of the most uncontaminated transitions and allows us to report the discovery of previously undetected transitions of the three isomers. With the least contaminated transitions, the high spatial-resolution millimeter maps of the C$_2$H$_4$O$_2$ isomers reveal that HCOOCH$_3$ and CH$_2$OHCHO each display two different velocity components, while only one velocity component of CH$_3$COOH is resolved. Moreover, the distribution of HCOOCH$_3$ is extended and offset from the continuum emission, unlike CH$_2$OHCHO and CH$_3$COOH, for which the low-velocity component is found to be compact and concentrated toward the continuum emission peak of Sgr B2(N). The distinct morphologies of these C$_2$H$_4$O$_2$ isomeric species indicate that HCOOCH$_3$ have significant differences in chemical processes than CH$_2$OHCHO and CH$_3$COOH, which display similar spatial distributions.

The emission line spectra of cyanoacetylene and methanol reveal chemical and physical heterogeneity on very small (< 0.1 pc) scales toward the peak in cyanopolyyne emission in the Taurus Molecular Cloud, TMC-1 (CP). We generate grids of homogeneous chemical models using a three-phase rate equation approach to obtain all time-dependent abundances spanning the physical conditions determined from molecular tracers of compact and extended regions of emission along this line of sight. Each time-dependent abundance is characterized by one of four features: a maximum/minimum, a monotonic increase/decrease, oscillatory behavior, or inertness. We similarly classify the time-dependent agreement between modeled and observed abundances by calculating both the root-mean-square logarithm difference and root-mean-square deviation between the modeled and observed abundances at every point in our grid models for three groups of molecules: (i) a composite group of all species present in both the observations and our chemical network G, (ii) the cyanopolyynes C = {HC$_3$N, HC$_5$N, HC$_7$N, HC$_9$N}, and (iii) the oxygen-containing organic species methanol and acetaldehyde S = {CH$_3$OH, CH$_3$CHO}. We discuss how the Bayesian uncertainties in the observed abundances constrain solutions within the grids of chemical models. The calculated best fit times at each grid point for each group are tabulated to reveal the minimum solution space of the grid models and the effects the Bayesian uncertainties have on the grid model solutions. The results of this approach separate the effect different physical conditions and model-free parameters have on reproducing accurately the abundances of different groups of observed molecular species.

In this paper, we present preliminary results illustrating the effect of cosmic rays on solid-phase chemistry in models of both TMC-1 and several sources with physical conditions identical to TMC-1 except for hypothetically enhanced ionization rates. Using a recent theory for the addition of cosmic ray-induced reactions to astrochemical models, we calculated the radiochemical yields, called $G$ values, for the primary dust grain ice-mantle constituents. We show that the inclusion of this non-thermal chemistry can lead to the formation of complex organic molecules from simpler ice-mantle constituents, even under cold core conditions. In addition to enriching ice-mantles, we find that these new radiation-chemical processes can lead to increased gas-phase abundances as well, particularly for HOCO, NO$_2$, HC$_2$O, methyl formate (HCOOCH$_3$), and ethanol (CH$_3$CH$_2$OH). These model results imply that HOCO - and perhaps NO$_2$ - might be observable in TMC-1. Future detections of either of these two species in cold interstellar environments could provide strong support for the importance of cosmic ray-driven radiation chemistry. The increased gas-phase abundance of methyl formate can be compared with abundances achieved through other formation mechanisms such as pure gas-phase chemistry and three-body surface reactions.

We investigate deuterium chemistry coupled with the nuclear spin-state chemistry of H$_2$ and H$_3^+$ in protoplanetary disks. Multiple paths of deuterium fractionation are found; exchange reactions with D atoms, such as HCO$^+$ + D, are effective in addition to those with HD. In a disk model with grain sizes appropriate for dark clouds, the freeze-out of molecules is severe in the outer midplane, while the disk surface is shielded from UV radiation. Gaseous molecules, including DCO$^+$, thus become abundant at the disk surface, which tends to make their column density distribution relatively flat. If the dust grains have grown to millimeter size, the freeze-out rate of neutral species is reduced, and the abundances of gaseous molecules, including DCO$^+$ and N$_2$D$^+$, are enhanced in the cold midplane. Turbulent diffusion transports D atoms and radicals at the disk surface to the midplane, and stable ice species in the midplane to the disk surface. The effects of turbulence on chemistry are thus multifold; while DCO$^+$ and N$_2$D$^+$ abundances increase or decrease depending on the regions, HCN and DCN in the gas and ice are much reduced at the innermost radii, compared with the model without turbulence. When cosmic rays penetrate the disk, the ortho-to-para ratio (OPR) of H$_2$ is found to be thermal in the disk, except in the cold ($\lesssim 10$ K) midplane. We also analyze the OPR of H$_3^+$ and H$_2$D$^+$, as well as the main reactions of H$_2$D$^+$, DCO$^+$, and N$_2$D$^+$ to analytically derive their abundances in the cold midplane.

We present high-spectral resolution observations from 12.96 - 13.33 microns towards Orion IRc2 using the mid-infrared spectrograph, EXES, on SOFIA. These observations probe the physical and chemical conditions of the Orion hot core, which is sampled by a bright, compact, mid-infrared background continuum source in the region, IRc2. All ten of the rovibrational C2H2 transitions expected in our spectral coverage, are detected with high S/N, yielding continuous coverage of the R-branch lines from J=9-8 to J=18-17, including both ortho and para species. Eight of these rovibrational transitions are newly reported detections. The isotopologue, 13CCH2, is clearly detected with high signal-to-noise. This enabled a direct measurement of the 12C/13C isotopic ratio for the Orion hot core of 14 +/- 1 and an estimated maximum value of 21. We also detected several HCN rovibrational lines. The ortho and para C2H2 ladders are clearly separate and tracing two different temperatures, 226 K and 164 K, respectively, with a non-equilibrium ortho to para ratio (OPR) of 1.7 +\- 0.1. Additionally, the ortho and para V_LSR values differ by about 1.8 +/- 0.2 km/s, while, the mean line widths differ by 0.7 +/- 0.2 km/s, suggesting that these species are not uniformly mixed along the line of sight to IRc2. We propose that the abnormally low C2H2 OPR could be a remnant from an earlier, colder phase, before the density enhancement (now the hot core) was impacted by shocks generated from an explosive event 500 yrs ago.

Polycyclic aromatic hydrocarbons and polycyclic aromatic nitrogen heterocycles are thought to be widespread throughout the Universe, because these classes of molecules are probably responsible for the unidentified infrared bands, a set of emission features seen in numerous Galactic and extragalactic sources. Despite their expected ubiquity, astronomical identification of specific aromatic molecules has proven elusive. We present the discovery of benzonitrile ($c$-C$_6$H$_5$CN), one of the simplest nitrogen-bearing aromatic molecules, in the interstellar medium. We observed hyperfine-resolved transitions of benzonitrile in emission from the molecular cloud TMC-1. Simple aromatic molecules such as benzonitrile may be precursors for polycyclic aromatic hydrocarbon formation, providing a chemical link to the carriers of the unidentified infrared bands.

It is important to investigate the relationships between the power sources and the chemical compositions of galaxies for understanding the scenario of galaxy evolution. We carried out an unbiased molecular line survey towards AGN host galaxy NGC1068, and prototypical starburst galaxies, NGC 253 and IC 342, with the Nobeyama 45-m telescope in the 3-mm band. The advantage of this line survey is that the obtained spectra have the highest angular resolution ever obtained with single-dish telescopes. In particular, the beam size of this telescope is ~15"--19", which is able to spatially separate the nuclear molecular emission from that of the starburst ring (d~30") in NGC 1068. We successfully detected approximately 23 molecular species in each galaxy, and calculated rotation temperatures and column densities. We estimate the molecular fractional abundances with respect to 13CO and CS molecules and compare them among three galaxies in order to investigate the chemical signatures of an AGN environment. As a result, we found clear trends on the abundances of molecules surrounding the AGN on 1 kpc scale. HCN, H13CN, CN, 13CN, and HC3N are more abundant, and CH3CCH is deficient in NGC 1068 compared with the starburst galaxies. High abundances of HCN, H13CN, and HC3N suggest that the circumnuclear disk in NGC 1068 is in a high-temperature environment. The reason for the non-detection of CH3CCH is likely to be dissociation by high energy radiation or less sublimation of a precursor of CH3CCH from grains.

We report the detection of interstellar methoxymethanol (CH$_3$OCH$_2$OH) in ALMA Bands 6 and 7 toward the MM1 core in the high-mass star-forming region NGC 6334I at ~0.1" - 1" spatial resolution. A column density of 4(2) x $10^{18}$ cm$^{-2}$ at $T_{ex}$ = 200 K is derived toward MM1, ~34 times less abundant than methanol (CH$_3$OH), and significantly higher than predicted by astrochemical models. Probable formation and destruction pathways are discussed, primarily through the reaction of the CH$_3$OH photodissociation products, the methoxy (CH$_3$O) and hydroxymethyl (CH$_2$OH) radicals. Finally, we comment on the implications of these mechanisms on gas-phase vs grain-surface routes operative in the region, and the possibility of electron-induced dissociation of CH$_3$OH rather than photodissociation.

We report the first interstellar detection of DC$_7$N and six $^{13}$C-bearing isotopologues of HC$_7$N toward the dark cloud TMC-1 through observations with the Green Bank Telescope, and confirm the recent detection of HC$_5$$^{15}$N. For the average of the $^{13}$C isotopomers, DC$_7$N, and HC$_5$$^{15}$N, we derive column densities of 1.9(2)$\times$10$^{11}$, 2.5(9)$\times$10$^{11}$, and 1.5(4)$\times$10$^{11}$ cm$^{-2}$, respectively. The resulting isotopic ratios are consistent with previous values derived from similar species in the source, and we discuss the implications for the formation chemistry of the observed cyanopolyynes. Within our uncertainties, no significant $^{13}$C isotopomer variation is found for HC$_7$N, limiting the significance CN could have in its production. The results further show that, for all observed isotopes, HC$_5$N may be isotopically depleted relative to HC$_3$N and HC$_7$N, suggesting that reactions starting from smaller cyanopolyynes may not be efficient to form HC$_{n}$N. This leads to the conclusion that the dominant production route may be the reaction between hydrocarbon ions and nitrogen atoms.

Multi-hydrogenated species with proper symmetry properties can present different spin configurations, and thus exist under different spin symmetry forms, labeled as para and ortho for two-hydrogen molecules. We investigated here the ortho-to-para ratio (OPR) of H$_2$Cl$^+$ in the light of new observations performed in the z=0.89 absorber toward the lensed quasar PKS 1830-211 with the Atacama Large Millimeter/submillimeter Array (ALMA). Two independent lines of sight were observed, to the southwest (SW) and northeast (NE) images of the quasar, with OPR values found to be $3.15 \pm 0.13$ and $3.1 \pm 0.5$ in each region, respectively, in agreement with a spin statistical weight of 3:1. An OPR of 3:1 for a molecule containing two identical hydrogen nuclei can refer to either a statistical result or a high-temperature limit depending on the reaction mechanism leading to its formation. It is thus crucial to identify rigorously how OPRs are produced in order to constrain the information that these probes can provide. To understand the production of the H$_2$Cl$^+$ OPR, we undertook a careful theoretical study of the reaction mechanisms involved with the aid of quasi-classical trajectory calculations on a new global potential energy surface fit to a large number of high-level ab initio data. Our study shows that the major formation reaction for H$_2$Cl$^+$ produces this ion via a hydrogen abstraction rather than a scrambling mechanism. Such a mechanism leads to a 3:1 OPR, which is not changed by destruction and possible thermalization reactions for H$_2$Cl$^+$ and is thus likely to be the cause of observed 3:1 OPR ratios, contrary to the normal assumption of scrambling.

The 1-50 GHz GBT PRIMOS data contains ~50 molecular absorption lines observed in diffuse and translucent clouds located in the Galactic Center, Bar, and spiral arms in the line-of-sight to Sgr B2(N). We measure the column densities and estimate abundances, relative to H2, of 11 molecules and additional isotopologues. We use absorption by optically thin transitions of c-C3H2 to estimate the N(H2), and argue that this method is preferable to more commonly used methods. We discuss the kinematic structure and abundance patterns of small molecules including the sulfur-bearing species CS, SO, CCS, H2CS, and HCS+; oxygen-bearing molecules OH, SiO, and H2CO; and simple hydrocarbon molecules c-C3H2, l-C3H, and l-C3H+. We discuss the implications of the observed chemistry for the structure of the gas and dust in the ISM. Highlighted results include the following. First, whereas gas in the disk has a molecular hydrogen fraction of 0.65, clouds on the outer edge of the Galactic Bar and in or near the Galactic Center have molecular fractions of 0.85 and >0.9, respectively. Second, we observe trends in isotope ratios with Galactocentric distance; while carbon and silicon show enhancement of the rare isotopes at low Galactocentric distances, sulfur exhibits no trend with Galactocentric distance; the ratio of c-C3H2/c-H13CCCH provides a good estimate of the 12C:13C ratio, whereas H2CO/H2^13CO exhibits fractionation. Third, we report the presence of l-C3H+ in diffuse clouds for the first time. Finally, we suggest that CS has an enhanced abundance within higher density clumps of material in the disk, and therefore may be diagnostic of cloud conditions. If this holds, the diffuse clouds in the Galactic disk contain multiple embedded hyperdensities in a clumpy structure, and the density profile is not a simple function of A_V.

A wide variety of molecules have recently been detected in the Horsehead nebula photodissociation region (PDR) suggesting that: (i) gas-phase and grain chemistries should both contribute to the formation of organic molecules, and (ii) far-ultraviolet (FUV) photodesorption may explain the release into the gas phase of grain surface species. In order to tackle these specific problems and more generally in order to better constrain the chemical structure of these types of environments we present a study of the Horsehead nebula gas-grain chemistry. To do so we used the 1D astrochemical gas-grain code Nautilus with an appropriate physical structure computed with the Meudon PDR Code and compared our modeled outcomes with published observations and with previously modeled results when available. The use of a large set of chemical reactions coupled with the time-dependent code Nautilus allows us to reproduce most of the observations well, including those of the first detections in a PDR of the organic molecules HCOOH, CH2CO, CH3CHO and CH3CCH, which are mostly associated with hot cores. We also provide some abundance predictions for other molecules of interest. Understanding the chemistry behind the detection of these organic molecules is crucial to better constrain the environments these molecules can probe.

We report the detection of the carbon-chain radical HC$_5$O for the first time in the interstellar medium toward the dark cloud TMC-1 using the 100 m Green Bank Telescope. We observe four hyperfine components of this radical in the $J = 17/2 \rightarrow 15/2$ rotational transition that originates from the $^2\Pi_{1/2}$ fine structure level of its ground state, and calculate an abundance of $n/n_{H_2}$ = $1.7\times 10^{-10}$, assuming an excitation temperature of $T_{ex} = 7$~K. No indication of HC$_3$O, HC$_4$O, HC$_6$O, is found in these or archival observations of the source, while we report tentative evidence for HC$_7$O. We compare calculated upper limits, and the abundance of HC$_5$O to predictions based on (1) the abundance trend of the analogous HC$_n$N family in TMC-1 and (2) a gas-grain chemical model. We find that the gas-grain chemical model well reproduces the observed abundance of HC$_5$O, as well as the upper limits of HC$_3$O, HC$_6$O, and HC$_7$O, but HC$_4$O is over produced. The prospects for astronomical detection of both shorter and longer HC$_n$O chains are discussed.

The collisions between high-energy ions and solids can result in significant physical and chemical changes to the material. These effects are potentially important for better understanding the chemistry of interstellar and planetary bodies, which are exposed to cosmic radiation and the solar wind, respectively, however, modeling such collisions on a detailed microscopic basis has thus far been largely unsuccessful. To that end, a new model, entitled CIRIS: the Chemistry of Ionizing Radiation in Solids, was created to calculate the physical and chemical effects of the irradiation of solid materials. With the new code, we simulate O2 ice irradiated with 100 keV protons. Our models are able to reproduce the measured ozone abundances of a previous experimental study, as well as independently predict the approximate thickness of the ice used in that work.

Propargyl alcohol (HC2CH2OH, PA) has yet to be observed in the interstellar medium (ISM) although one of its stable isomers, propenal (CH2CHCHO), has already been detected in Sagittarius B2(N) with the 100-meter Green Bank Telescope in the frequency range $18-26$ GHz. In this paper, we investigate the formation of propargyl alcohol along with one of its deuterated isotopomers, HC2CH2OD(OD-PA), in a dense molecular cloud. Various pathways for the formation of PA in the gas and on ice mantles surrounding dust particles are discussed. We use a large gas-grain chemical network to study the chemical evolution of PA and its deuterated isotopomer. Our results suggest that gaseous HC2CH2OH can most likely be detected in hot cores or in collections of hot cores such as the star-forming region Sgr B2(N). A simple LTE (Local thermodynamic equilibrium) radiative transfer model is employed to check the possibility of detecting PA and OD-PA in the millimeter-wave regime. In addition, we have carried out quantum chemical calculations to compute the vibrational transition frequencies and intensities of these species in the infrared for perhaps future use in studies with the James Webb Space Telescope (JWST).

Because of the very peculiar conditions of chemistry in many astrophysical gases (low densities, mostly low temperatures, kinetics-dominated chemical evolution), great efforts have been devoted to study molecular signatures and chemical evolution. While experiments are being performed in many laboratories, it appears that the efforts directed towards theoretical works are not as strong. This report deals with the present status of chemical physics/physical chemistry theory, for the qualitative and quantitative understanding of kinetics of molecular scattering, being it reactive or inelastic. By gathering several types of expertise, from applied mathematics to physical chemistry, dialog is made possible, as a step towards new and more adapted theoretical frameworks, capable of meeting the theoretical, methodological and numerical challenges of kinetics-dominated gas phase chemistry in astrophysical environments. A state of the art panorama is presented, alongside present-day strengths and shortcomings. However, coverage is not complete, being limited in this report to actual attendance of the workshop. Some paths towards relevant progress are proposed.

The chemistry of dense interstellar regions was analyzed using a time-dependent gas-grain astrochemical simulation and a new chemical network that incorporates deuterated chemistry taking into account nuclear spin-states for the hydrogen chemistry and its deuterated isotopologues. With this new network, the utility of the [HCO$^+$]/[DCO$^+$] abundance ratio as a probe of the cosmic ray ionization rate has been reexamined, with special attention paid to the effect of the initial value of the molecular hydrogen ortho-to-para ratio (OPR). After discussing the use of the probe for cold cores, we then compare our results with previous theoretical and observational results for a molecular cloud close to the supernova remnant W51C, which is thought to have an enhanced cosmic ray ionization rate $\zeta$ caused by the nearby $\gamma$-ray source. In addition, we attempt to use our approach to estimate the cosmic ray ionization rate for L1174, a dense core with an embedded star. Beyond the previously known sensitivity of [HCO$^+$]/[DCO$^+$] to $\zeta$, we demonstrate its additional dependence on the initial OPR and, secondarily, on the age of the source, its temperature, and its density. We conclude that the usefulness of the [HCO$^+$]/[DCO$^+$] abundance ratio to constrain the cosmic ray ionization rate in dense regions increases with source age and ionization rate as the ratio becomes far less sensitive to the initial value of the OPR.

Based on recent $Herschel$ results, the ortho-to-para ratio (OPR) of NH$_2$ has been measured towards the following high-mass star-forming regions: W31C (G10.6-0.4), W49N (G43.2-0.1), W51 (G49.5-0.4), and G34.3+0.1. The OPR at thermal equilibrium ranges from the statistical limit of three at high temperatures to infinity as the temperature tends toward zero, unlike the case of H$_{2}$. Depending on the position observed along the lines-of-sight, the OPR was found to lie either slightly below the high temperature limit of three (in the range $2.2-2.9$) or above this limit ($\sim3.5$, $\gtrsim 4.2$, and $\gtrsim 5.0$). In low temperature interstellar gas, where the H$_{2}$ is para-enriched, our nearly pure gas-phase astrochemical models with nuclear-spin chemistry can account for anomalously low observed NH$_2$-OPR values. We have tentatively explained OPR values larger than three by assuming that spin thermalization of NH$_2$ can proceed at least partially by H-atom exchange collisions with atomic hydrogen, thus increasing the OPR with decreasing temperature. In this paper, we present quasi-classical trajectory calculations of the H-exchange reaction NH$_2$ + H, which show the reaction to proceed without a barrier, confirming that the H-exchange will be efficient in the temperature range of interest. With the inclusion of this process, our models suggest both that OPR values below three arise in regions with temperatures $\gtrsim20-25$~K, depending on time, and values above three but lower than the thermal limit arise at still lower temperatures.

The first detection of gas-phase methanol in a protoplanetary disk (TW Hya) is presented. In addition to being one of the largest molecules detected in disks to date, methanol is also the first disk organic molecule with an unambiguous ice chemistry origin. The stacked methanol emission, as observed with ALMA, is spectrally resolved and detected across six velocity channels ($>3 \sigma$), reaching a peak signal-to-noise of $5.5\sigma$, with the kinematic pattern expected for TW~Hya. Using an appropriate disk model, a fractional abundance of $3\times 10^{-12} - 4 \times 10^{-11}$ (with respect to H$_2$) reproduces the stacked line profile and channel maps, with the favoured abundance dependent upon the assumed vertical location (midplane versus molecular layer). The peak emission is offset from the source position suggesting that the methanol emission has a ring-like morphology: the analysis here suggests it peaks at $\approx 30$~AU reaching a column density $\approx 3-6\times10^{12}$~cm$^{-2}$. In the case of TW Hya, the larger (up to mm-sized) grains, residing in the inner 50~AU, may thus host the bulk of the disk ice reservoir. The successful detection of cold gas-phase methanol in a protoplanetary disk implies that the products of ice chemistry can be explored in disks, opening a window to studying complex organic chemistry during planetary system formation.

L1157, a molecular dark cloud with an embedded Class 0 protostar possessing a bipolar outflow, is an excellent source for studying shock chemistry, including grain-surface chemistry prior to shocks, and post-shock, gas-phase processing. The L1157-B1 and B2 positions experienced shocks at an estimated ~2000 and 4000 years ago, respectively. Prior to these shock events, temperatures were too low for most complex organic molecules to undergo thermal desorption. Thus, the shocks should have liberated these molecules from the ice grain-surfaces en masse, evidenced by prior observations of SiO and multiple grain mantle species commonly associated with shocks. Grain species, such as OCS, CH3OH, and HNCO, all peak at different positions relative to species that are preferably formed in higher velocity shocks or repeatedly-shocked material, such as SiO and HCN. Here, we present high spatial resolution (~3") maps of CH3OH, HNCO, HCN, and HCO+ in the southern portion of the outflow containing B1 and B2, as observed with CARMA. The HNCO maps are the first interferometric observations of this species in L1157. The maps show distinct differences in the chemistry within the various shocked regions in L1157B. This is further supported through constraints of the molecular abundances using the non-LTE code RADEX (Van der Tak et al. 2007). We find the east/west chemical differentiation in C2 may be explained by the contrast of the shock's interaction with either cold, pristine material or warm, previously-shocked gas, as seen in enhanced HCN abundances. In addition, the enhancement of the HNCO abundance toward the the older shock, B2, suggests the importance of high-temperature O-chemistry in shocked regions.

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.

The first hydrostatic core, also called the first Larson core, is one of the first steps in low-mass star formation, as predicted by theory. With recent and future high performance telescopes, details of these first phases become accessible, and observations may confirm theory and even bring new challenges for theoreticians. In this context, we study from a theoretical point of view the chemical and physical evolution of the collapse of prestellar cores until the formation of the first Larson core, in order to better characterize this early phase in the star formation process. We couple a state-of-the-art hydrodynamical model with full gas-grain chemistry, using different assumptions on the magnetic field strength and orientation. We extract the different components of each collapsing core (i.e., the central core, the outflow, the disk, the pseudodisk, and the envelope) to highlight their specific physical and chemical characteristics. Each component often presents a specific physical history, as well as a specific chemical evolution. From some species, the components can clearly be differentiated. The different core models can also be chemically differentiated. Our simulation suggests some chemical species as tracers of the different components of a collapsing prestellar dense core, and as tracers of the magnetic field characteristics of the core. From this result, we pinpoint promising key chemical species to be observed.

We have used the Herschel-HIFI instrument to observe both nuclear spin symmetries of amidogen (NH2) towards the high-mass star-forming regions W31C (G10.6-0.4), W49N (G43.2-0.1), W51 (G49.5-0.4) and G34.3+0.1. The aim is to investigate the ratio of nuclear spin types, the ortho-to-para ratio (OPR), of NH2. The excited NH2 transitions are used to construct radiative transfer models of the hot cores and surrounding envelopes in order to investigate the excitation and possible emission of the ground state rotational transitions of ortho-NH2 N_(K_a,K_c} J=1_(1,1) 3/2 - 0_(0,0) 1/2 and para-NH2 2_(1,2) 5/2 - 1_(0,1) 3/2$ used in the OPR calculations. Our best estimate of the average OPR in the envelopes lie above the high temperature limit of three for W49N, specifically 3.5 with formal errors of \pm0.1, but for W31C, W51, and G34.3+0.1 we find lower values of 2.5\pm0.1, 2.7\pm0.1, and 2.3\pm0.1, respectively. Such low values are strictly forbidden in thermodynamical equilibrium since the OPR is expected to increase above three at low temperatures. In the translucent interstellar gas towards W31C, where the excitation effects are low, we find similar values between 2.2\pm0.2 and 2.9\pm0.2. In contrast, we find an OPR of 3.4\pm0.1 in the dense and cold filament connected to W51, and also two lower limits of >4.2 and >5.0 in two other translucent gas components towards W31C and W49N. At low temperatures (T \lesssim 50 K) the OPR of H2 is <10^-1, far lower than the terrestrial laboratory normal value of three. In such a "para-enriched H2" gas, our astrochemical models can reproduce the variations of the observed OPR, both below and above the thermodynamical equilibrium value, by considering nuclear-spin gas-phase chemistry. The models suggest that values below three arise in regions with temperatures >20-25 K, depending on time, and values above three at lower temperatures.

We investigate the water deuteration ratio and ortho-to-para nuclear spin ratio of H2 (OPR(H2)) during the formation and early evolution of a molecular cloud, following the scenario that accretion flows sweep and accumulate HI gas to form molecular clouds. We follow the physical evolution of post-shock materials using a one-dimensional shock model, with post-processing gas-ice chemistry simulations. This approach allows us to study the evolution of the OPR(H2) and water deuteration ratio without an arbitrary assumption concerning the initial molecular abundances, including the initial OPR(H2). When the conversion of hydrogen into H2 is almost complete, the OPR(H2) is already much smaller than the statistical value of three due to the spin conversion in the gas phase. As the gas accumulates, the OPR(H2) decreases in a non-equilibrium manner. We find that water ice can be deuterium-poor at the end of its main formation stage in the cloud, compared to water vapor observed in the vicinity of low-mass protostars where water ice is likely sublimated. If this is the case, the enrichment of deuterium in water should mostly occur at somewhat later evolutionary stages of star formation, i.e., cold prestellar/protostellar cores. The main mechanism to suppress water ice deuteration in the cloud is the cycle of photodissociation and reformation of water ice, which efficiently removes deuterium from water ice chemistry. The removal efficiency depends on the main formation pathway of water ice. The OPR(H2) plays a minor role in water ice deuteration at the main formation stage of water ice.

We present chemical implications arising from spectral models fit to the Herschel/HIFI spectral survey toward the Orion Kleinmann-Low nebula (Orion KL). We focus our discussion on the eight complex organics detected within the HIFI survey utilizing a novel technique to identify those molecules emitting in the hottest gas. In particular, we find the complex nitrogen bearing species CH$_{3}$CN, C$_{2}$H$_{3}$CN, C$_{2}$H$_{5}$CN, and NH$_{2}$CHO systematically trace hotter gas than the oxygen bearing organics CH$_{3}$OH, C$_{2}$H$_{5}$OH, CH$_{3}$OCH$_{3}$, and CH$_{3}$OCHO, which do not contain nitrogen. If these complex species form predominantly on grain surfaces, this may indicate N-bearing organics are more difficult to remove from grain surfaces than O-bearing species. Another possibility is that hot (T$_{\rm kin}$$\sim$300 K) gas phase chemistry naturally produces higher complex cyanide abundances while suppressing the formation of O-bearing complex organics. We compare our derived rotation temperatures and molecular abundances to chemical models, which include gas-phase and grain surface pathways. Abundances for a majority of the detected complex organics can be reproduced over timescales $\gtrsim$ 10$^{5}$ years, with several species being under predicted by less than 3$\sigma$. Derived rotation temperatures for most organics, furthermore, agree reasonably well with the predicted temperatures at peak abundance. We also find that sulfur bearing molecules which also contain oxygen (i.e. SO, SO$_{2}$, and OCS) tend to probe the hottest gas toward Orion KL indicating the formation pathways for these species are most efficient at high temperatures.

We report additional detections of the chloronium molecular ion, H$_2$Cl$^+$, toward four bright submillimeter continuum sources: G29.96, W49N, W51, and W3(OH). With the use of the HIFI instrument on the Herschel Space Observatory, we observed the $2_{12}-1_{01}$ transition of ortho-H$_2^{35}$Cl$^+$ at 781.627 GHz in absorption toward all four sources. Much of the detected absorption arises in diffuse foreground clouds that are unassociated with the background continuum sources and in which our best estimates of the $N({\rm H_2Cl^+})/N({\rm H})$ ratio lie in the range $(0.9 - 4.8) \times 10^{-9}$. These chloronium abundances relative to atomic hydrogen can exceed the predictions of current astrochemical models by up to a factor of 5. Toward W49N, we have also detected the $2_{12}-1_{01}$ transition of ortho-H$_2^{37}$Cl$^+$ at 780.053 GHz and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ at 485.418 GHz. These observations imply $\rm H_2^{35}Cl^+/H_2^{37}Cl^+$ column density ratios that are consistent with the solar system $^{35}$Cl/$^{37}$Cl isotopic ratio of 3.1, and chloronium ortho-to-para ratios consistent with 3, the ratio of spin statistical weights.

Chemical models used to study the chemical composition of the gas and the ices in the interstellar medium are based on a network of chemical reactions and associated rate coefficients. These reactions and rate coefficients are partially compiled from data in the literature, when available. We present in this paper kida.uva.2014, a new updated version of the kida.uva public gas-phase network first released in 2012. In addition to a description of the many specific updates, we illustrate changes in the predicted abundances of molecules for cold dense cloud conditions as compared with the results of the previous version of our network, kida.uva.2011.

We have observed five sulphur-bearing molecules in foreground diffuse molecular clouds lying along the sight-lines to five bright continuum sources. We have used the GREAT instrument on SOFIA to observe the 1383 GHz $^2\Pi_{3/2} J=5/2-3/2$ transitions of SH towards the star-forming regions W31C, G29.96-0.02, G34.3+0.1, W49N and W51, detecting foreground absorption towards all five sources; and the EMIR receivers on the IRAM 30m telescope at Pico Veleta to detect the H$_2$S 1(10)-1(01), CS J=2-1 and SO 3(2)-2(1) transitions. In nine foreground absorption components detected towards these sources, the inferred column densities of the four detected molecules showed relatively constant ratios, with N(SH)/N(H$_2$S) in the range 1.1 - 3.0, N(CS)/N(H$_2$S) in the range 0.32 - 0.61, and N(SO)/N(H$_2$S) in the range 0.08 - 0.30. The observed SH/H$_2$ ratios - in the range (0.5-2.6) $\times 10^{-8}$ - indicate that SH (and other sulphur-bearing molecules) account for << 1% of the gas-phase sulphur nuclei. The observed abundances of sulphur-bearing molecules, however, greatly exceed those predicted by standard models of cold diffuse molecular clouds, providing further evidence for the enhancement of endothermic reaction rates by elevated temperatures or ion-neutral drift. We have considered the observed abundance ratios in the context of shock and turbulent dissipation region (TDR) models. Using the TDR model, we find that the turbulent energy available at large scale in the diffuse ISM is sufficient to explain the observed column densities of SH and CS. Standard shock and TDR models, however, fail to reproduce the column densities of H$_2$S and SO by a factor of about 10; more elaborate shock models - in which account is taken of the velocity drift, relative to H$_2$, of SH molecules produced by the dissociative recombination of H$_3$S$^+$ - reduce this discrepancy to a factor ~ 3.

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 present 0.8-mm band molecular images and spectra obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) toward one of the nearest galaxies with an active galactic nucleus (AGN), NGC 1068. Distributions of CO isotopic species ($^{13}$CO and C$^{18}$O) $\it{J}$ = 3--2, CN $\it{N}$ = 3--2 and CS $\it{J}$ = 7--6 are observed toward the circumnuclear disk (CND) and a part of the starburst ring with an angular resolution of $\sim$1.$^{\prime\prime}$3 $\times$ 1.$^{\prime\prime}$2. The physical properties of these molecules and shock-related molecules such as HNCO, CH$_{3}$CN, SO, and CH$_{3}$OH detected in the 3-mm band were estimated using rotation diagrams under the assumption of local thermodynamic equilibrium. The rotational temperatures of the CO isotopic species and the shock-related molecules in the CND are, respectively, 14--22 K and upper limits of 20--40 K. Although the column densities of the CO isotopic species in the CND are only from one-fifth to one-third of that in the starburst ring, those of the shock-related molecules are enhanced by a factor of 3--10 in the CND. We also discuss the chemistry of each species, and compare the fractional abundances in the CND and starburst ring with those of Galactic sources such as cold cores, hot cores, and shocked molecular clouds in order to study the overall characteristics. We find that the abundances of shock-related molecules are more similar to abundances in hot cores and/or shocked clouds than to cold cores. The CND hosts relatively complex molecules, which are often associated with shocked molecular clouds or hot cores. Because a high X-ray flux can dissociate these molecules, they must also reside in regions shielded from X-rays.

Understanding water deuterium fractionation is important for constraining the mechanisms of water formation in interstellar clouds. Observations of HDO and H$_2^{18}$O transitions were carried out towards the high-mass star-forming region G34.26+0.15 with the HIFI instrument onboard the Herschel Space Observatory, as well as with ground-based single-dish telescopes. Ten HDO lines and three H$_2^{18}$O lines covering a broad range of upper energy levels (22-204 K) were detected. We used a non-LTE 1D analysis to determine the HDO/H$_2$O ratio as a function of radius in the envelope. Models with different water abundance distributions were considered in order to reproduce the observed line profiles. The HDO/H$_2$O ratio is found to be lower in the hot core ($\sim$3.5 $\times$ 10$^{-4}$ - 7.5 $\times$ 10$^{-4}$) than in the colder envelope ($\sim$1.0 $\times$ 10$^{-3}$ - 2.2 $\times$ 10$^{-3}$). This is the first time that a radial variation of the HDO/H$_2$O ratio has been found to occur in a high-mass source. The chemical evolution of this source was modeled as a function of its radius and the observations are relatively well reproduced. The comparison between the chemical model and the observations leads to an age of $\sim$10$^5$ years after the infrared dark cloud stage.

The understanding of interstellar nitrogen chemistry has improved significantly with recent results from the Herschel Space Observatory. To set even better constraints, we report here on deep searches for the NH+ ground state rotational transition J=1.5-0.5 of the ^2Pi_1/2 lower spin ladder, with fine-structure transitions at 1013 and 1019 GHz, and the para-NH2- 1_1,1-0_0,0 rotational transition at 934 GHz towards Sgr B2(M) and G10.6-0.4 using Herschel-HIFI. No clear detections of NH+ are made and the derived upper limits are <2*10^-12 and <7*10^-13 in Sgr B2(M) and G10.6-0.4, respectively. The searches are complicated by the fact that the 1013 GHz transition lies only -2.5 km/s from a CH2NH line, seen in absorption in Sgr B2(M), and that the hyperfine structure components in the 1019 GHz transition are spread over 134 km/s. Searches for the so far undetected NH2- anion turned out to be unfruitful towards G10.6-0.4, while the para-NH2- 1_1,1-0_0,0 transition was tentatively detected towards Sgr B2(M) at a velocity of 19 km/s. Assuming that the absorption occurs at the nominal source velocity of +64 km/s, the rest frequency would be 933.996 GHz, offset by 141 MHz from our estimated value. Using this feature as an upper limit, we found N(p-NH2-)<4*10^11 cm^-2. The upper limits for both species in the diffuse line-of-sight gas are less than 0.1 to 2 % of the values found for NH, NH2, and NH3 towards both sources. Chemical modelling predicts an NH+ abundance a few times lower than our present upper limits in diffuse gas and under typical Sgr B2(M) envelope conditions. The NH2- abundance is predicted to be several orders of magnitudes lower than our observed limits, hence not supporting our tentative detection. Thus, while NH2- may be very difficult to detect in interstellar space, it could, be possible to detect NH+ in regions where the ionisation rates of H2 and N are greatly enhanced.

(Abridged) The birth environment of the Sun will have influenced the conditions in the pre-solar nebula, including the attainable chemical complexity, important for prebiotic chemistry. The formation and distribution of complex organic molecules (COMs) in a disk around a T Tauri star is investigated for two scenarios: (i) an isolated disk, and (ii) a disk irradiated externally by a nearby massive star. The chemistry is calculated along the accretion flow from the outer disk inwards using a comprehensive network. Two simulations are performed, one beginning with complex ices and one with simple ices only. For the isolated disk, COMs are transported without major alteration into the inner disk where they thermally desorb into the gas reaching an abundance representative of the initial assumed ice abundance. For simple ices, COMs efficiently form on grain surfaces under the conditions in the outer disk. Gas-phase COMs are released into the molecular layer via photodesorption. For the irradiated disk, complex ices are also transported inwards; however, they undergo thermal processing caused by the warmer conditions in the irradiated disk which tends to reduce their abundance along the accretion flow. For simple ices, grain-surface chemistry cannot synthesise COMs in the outer disk because the necessary grain-surface radicals, which tend to be particularly volatile, are not sufficiently abundant on the grain surfaces. Gas-phase COMs are formed in the inner region of the irradiated disk via gas-phase chemistry induced by the desorption of strongly bound molecules such as methanol; hence, the abundances are not representative of the initial molecular abundances injected into the outer disk. These results suggest that the composition of comets formed in isolated disks may differ from those formed in externally irradiated disks with the latter composed of more simple ices.

Sensitive observations with ALMA allow astronomers to observe the detailed distributions of molecules with relatively weak intensity in nearby galaxies. In particular, we report distributions of several molecular transitions including shock and dust related species ($^{13}$CO $J$ = 1--0, C$^{18}$O $J$ = 1--0, $^{13}$CN $N$ = 1--0, CS $J$ = 2--1, SO $J_N$ = 3$_2$--2$_1$, HNCO $J_{Ka,Kc}$ = 5$_{0,5}$--4$_{0,4}$, HC$_3$N $J$ = 11--10, 12--11, CH$_3$OH $J_K$ = 2$_K$--1$_K$, and CH$_3$CN $J_K$ = 6$_K$--5$_K$) in the nearby Seyfert 2 galaxy NGC 1068 observed with the ALMA early science program. The central $\sim$1 arcmin ($\sim$4.3 kpc) of this galaxy was observed in the 100 GHz region covering $\sim$96--100 GHz and $\sim$108--111 GHz with an angular resolution of $\sim4"\times2"$ (290 pc$\times$140 pc) to study the effects of an active galactic nucleus and its surrounding starburst ring on molecular abundances. Here, we present images and report a classification of molecular distributions into three main categories: (1) Molecules concentrated in the circumnuclear disk (CND) (SO $J_N$ = 3$_2$--2$_1$, HC$_3$N $J$ = 11--10, 12--11, and CH$_3$CN $J_K$ = 6$_K$--5$_K$), (2) Molecules distributed both in the CND and the starburst ring (CS $J$ = 2--1 and CH$_3$OH $J_K$ = 2$_K$--1$_K$), (3) Molecules distributed mainly in the starburst ring ($^{13}$CO $J$ = 1--0 and C$^{18}$O $J$ = 1--0). Since most of the molecules such as HC$_3$N observed in the CND are easily dissociated by UV photons and X-rays, our results indicate that these molecules must be effectively shielded. In the starburst ring, the relative intensity of methanol at each clumpy region is not consistent with those of $^{13}$CO, C$^{18}$O, and CS. This difference is probably caused by the unique formation and destruction mechanisms of CH$_3$OH.

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.

HNC and HCN, typically used as dense gas tracers in molecular clouds, are a pair of isomers that have great potential as a temperature probe because of temperature dependent, isomer-specific formation and destruction pathways. Previous observations of the HNC/HCN abundance ratio show that the ratio decreases with increasing temperature, something that standard astrochemical models cannot reproduce. We have undertaken a detailed parameter study on which environmental characteristics and chemical reactions affect the HNC/HCN ratio and can thus contribute to the observed dependence. Using existing gas and gas-grain models updated with new reactions and reaction barriers, we find that in static models the H + HNC gas-phase reaction regulates the HNC/HCN ratio under all conditions, except for very early times. We quantitively constrain the combinations of H abundance and H + HNC reaction barrier that can explain the observed HNC/HCN temperature dependence and discuss the implications in light of new quantum chemical calculations. In warm-up models, gas-grain chemistry contributes significantly to the predicted HNC/HCN ratio and understanding the dynamics of star formation is therefore key to model the HNC/HCN system.

[Abridged] Ethylene oxide and its isomer acetaldehyde are important complex organic molecules because of their potential role in the formation of amino acids. Despite the fact that acetaldehyde is ubiquitous in the interstellar medium, ethylene oxide has not yet been detected in cold sources. We aim to understand the chemistry of the formation and loss of ethylene oxide in hot and cold interstellar objects (i) by including in a revised gas-grain network some recent experimental results on grain surfaces and (ii) by comparison with the chemical behaviour of its isomer, acetaldehyde. We test the code for the case of a hot core. The model allows us to predict the gaseous and solid ethylene oxide abundances during a cooling-down phase prior to star formation and during the subsequent warm-up phase. We can therefore predict at what temperatures ethylene oxide forms on grain surfaces and at what temperature it starts to desorb into the gas phase. The model reproduces the observed gaseous abundances of ethylene oxide and acetaldehyde towards high-mass star-forming regions. In addition, our results show that ethylene oxide may be present in outer and cooler regions of hot cores where its isomer has already been detected. Despite their different chemical structures, the chemistry of ethylene oxide is coupled to that of acetaldehyde, suggesting that acetaldehyde may be used as a tracer for ethylene oxide towards cold cores.

(Abridged) Protoplanetary disks are vital objects in star and planet formation, possessing all the material which may form a planetary system orbiting the new star. We investigate the synthesis of complex organic molecules (COMs) in disks to constrain the achievable chemical complexity and predict species and transitions which may be observable with ALMA. We have coupled a 2D model of a protoplanetary disk around a T Tauri star with a gas-grain chemical network including COMs. We compare compare synthesised line intensities and calculated column densities with observations and determine those COMs which may be observable in future. COMs are efficiently formed in the disk midplane via grain-surface chemical reactions, reaching peak grain-surface fractional abundances 1e-6 - 1e-4 that of the H nuclei number density. COMs formed on grain surfaces are returned to the gas phase via non-thermal desorption; however, gas-phase species reach lower fractional abundances than their grain-surface equivalents, 1e-12 - 1e-7. Including the irradiation of grain mantle material helps build further complexity in the ice through the replenishment of grain-surface radicals which take part in further grain-surface reactions. There is reasonable agreement with several line transitions of H2CO observed towards several T Tauri star-disk systems. The synthesised line intensities for CH3OH are consistent with upper limits determined towards all sources. Our models suggest CH3OH should be readily observable in nearby protoplanetary disks with ALMA; however, detection of more complex species may prove challenging. Our grain-surface abundances are consistent with those derived from cometary comae observations providing additional evidence for the hypothesis that comets (and other planetesimals) formed via the coagulation of icy grains in the Sun's natal disk.

Water is observed throughout the universe, from diffuse interstellar clouds to protoplanetary disks around young stars, and from comets in our own solar system and exoplanetary atmospheres to galaxies at high redshifts. This review summarizes the spectroscopy and excitation of water in interstellar space as well as the basic chemical processes that form and destroy water under interstellar conditions. Three major routes to water formation are identified: low temperature ion-molecule chemistry, high-temperature neutral-neutral chemistry and gas-ice chemistry. The rate coefficients of several important processes entering the networks are discussed in detail; several of them have been determined only in the last decade through laboratory experiments and theoretical calculations. Astronomical examples of each of the different chemical routes are presented using data from powerful new telescopes, in particular the Herschel Space Observatory. Basic chemical physics studies remain critically important to analyze astronomical data.

It is currently assumed that infrared dark clouds (IRDCs) represent the earliest evolutionary stages of high-mass stars ($>$ 8 M$_{\odot}$). Submillimeter and millimeter-wave studies performed over the past 15 years show that IRDCs possess a broad variety of properties, and hence a wide range of problems and questions that can be tackled. In this paper, we report an investigation of the molecular composition and chemical processes in two groups of IRDCs. Using the Mopra, APEX, and IRAM radio telescopes over the last four years, we have collected molecular line data for CO, H$_2$CO, HNCO, CH$_3$CCH, CH$_3$OH, CH$_3$CHO, CH$_3$OCHO, and CH$_3$OCH$_3$. For all of these species we estimated molecular abundances. We then undertook chemical modeling studies, concentrating on the source IRDC028.34+0.06, and compared observed and modeled abundances. This comparison showed that to reproduce observed abundances of complex organic molecules (COMs), a 0-D gas-grain model with constant physical conditions is not sufficient. We achieved greater success with the use of a warm-up model, in which warm-up from 10 K to 30 K occurs following a cold phase.

We discuss models that astrochemists have developed to study the chemical composition of the interstellar medium. These models aim at computing the evolution of the chemical composition of a mixture of gas and dust under as- trophysical conditions. These conditions, as well as the geometry and the physical dynamics, have to be adapted to the objects being studied because different classes of objects have very different characteristics (temperatures, densities, UV radia- tion fields, geometry, history etc); e.g., proto-planetary disks do not have the same characteristics as protostellar envelopes. Chemical models are being improved continually thanks to comparisons with observations but also thanks to laboratory and theoretical work in which the individual processes are studied.

According to traditional gas-phase chemical models, O2 should be abundant in molecular clouds, but until recently, attempts to detect interstellar O2 line emission with ground- and space-based observatories have failed. Following the multi-line detections of O2 with low abundances in the Orion and rho Oph A molecular clouds with Herschel, it is important to investigate other environments, and we here quantify the O2 abundance near a solar-mass protostar. Observations of O2, at 487 GHz toward a deeply embedded low-mass Class 0 protostar, NGC 1333-IRAS 4A, are presented, using the HIFI instrument on the Herschel Space Observatory. Complementary data of the chemically related NO and CO molecules are obtained as well. The high spectral resolution data are analysed using radiative transfer models to infer column densities and abundances, and are tested directly against full gas-grain chemical models. The deep HIFI spectrum fails to show O2 at the velocity of the dense protostellar envelope, implying one of the lowest abundance upper limits of O2/H2 at <6x10^-9 (3 sigma). However, a tentative (4.5 sigma) detection of O2 is seen at the velocity of the surrounding NGC 1333 molecular cloud, shifted by 1 km/s relative to the protostar. For the protostellar envelope, pure gas-phase models and gas-grain chemical models require a long pre-collapse phase (~0.7-1x10^6 years), during which atomic and molecular oxygen are frozen out onto dust grains and fully converted to H2O, to avoid overproduction of O2 in the dense envelope. The same model also reproduces the limits on the chemically related NO molecule. The tentative detection of O2 in the surrounding cloud is consistent with a low-density PDR model with small changes in reaction rates. The low O2 abundance in the collapsing envelope around a low-mass protostar suggests that the gas and ice entering protoplanetary disks is very poor in O2.

Observations of deuterated species are useful in probing the temperature, ionization level, evolutionary stage, chemistry, and thermal history of astrophysical environments. The analysis of data from ALMA and other new telescopes requires an elaborate model of deuterium fractionation. This paper presents a publicly available chemical network with multi-deuterated species and an extended, up-to-date set of gas-phase and surface reactions. To test this network, we simulate deuterium fractionation in diverse interstellar sources. Two cases of initial abundances are considered: i) atomic except for H2 and HD, and ii) molecular from a prestellar core. We reproduce the observed D/H ratios of many deuterated molecules, and sort the species according to their sensitivity to temperature gradients and initial abundances. We find that many multiply-deuterated species produced at 10 K retain enhanced D/H ratios at temperatures $\la 100$ K. We study how recent updates to reaction rates affect calculated D/H ratios, and perform a detailed sensitivity analysis of the uncertainties of the gas-phase reaction rates in the network. We find that uncertainties are generally lower in dark cloud environments than in warm IRDCs and that uncertainties increase with the size of the molecule and number of D-atoms. A set of the most problematic reactions is presented. We list potentially observable deuterated species predicted to be abundant in low- and high-mass star-formation regions.

The recent discovery of terrestrial-type organic species such as methyl formate and dimethyl ether in the cold interstellar gas has proved that the formation of organic matter in the Galaxy begins at a much earlier stage of star formation than was thought before. This discovery represents a challenge for astrochemical modelers. The abundances of these molecules cannot be explained by the previously developed "warm-up" scenario, in which organic molecules are formed via diffusive chemistry on surfaces of interstellar grains starting at 30 K, and then released to the gas at higher temperatures during later stages of star formation. In this article, we investigate an alternative scenario in which complex organic species are formed via a sequence of gas-phase reactions between precursor species formed on grain surfaces and then ejected into the gas via efficient reactive desorption, a process in which non-thermal desorption occurs as a result of conversion of the exothermicity of chemical reactions into the ejection of products from the surface. The proposed scenario leads to reasonable if somewhat mixed results at temperatures as low as 10 K and may be considered as a step towards the explanation of abundances of terrestrial-like organic species observed during the earliest stages of star formation.

We use a high-temperature chemical network to derive the molecular abundances in axisymmetric accretion disk models around active galactic nuclei (AGNs) within 100 pc using simple radial and vertical density and temperature distributions motivated by more detailed physical models. We explore the effects of X-ray irradiation and cosmic ray ionization on the spatial distribution of the molecular abundances of CO, CN, CS, HCN, HCO+, HC3N, C2H, and c-C3H2 using a variety of plausible disk structures. These simple models have molecular regions with a layer of X-ray dominated regions, a midplane without the strong influence of X-rays, and a high-temperature region in the inner portion with moderate X-ray flux where families of polyynes (C$_{\rm n}$H$_{2}$) and cyanopolyynes can be enhanced. For the high midplane density disks we explore, we find that cosmic rays produced by supernovae do not significantly affect the regions unless the star formation efficiency significantly exceeds that of the Milky Way. We highlight molecular abundance observations and ratios that may distinguish among theoretical models of the density distribution in AGN disks. Finally, we assess the importance of the shock crossing time and the accretion time relative to the formation time for various chemical species. Vertical column densities are tabulated for a number of molecular species at both the characteristic shock crossing time and steady state. Although we do not attempt to fit any particular system or set of observations, we discuss our models and results in the context of the nearby AGN NGC 1068.

The observed gas-phase molecular inventory of hot cores is believed to be significantly impacted by the products of chemistry in interstellar ices. In this study, we report the construction of a full macroscopic Monte Carlo model of both the gas-phase chemistry and the chemistry occurring in the icy mantles of interstellar grains. Our model treats icy grain mantles in a layer-by-layer manner, which incorporates laboratory data on ice desorption correctly. The ice treatment includes a distinction between a reactive ice surface and an inert bulk. The treatment also distinguishes between zeroth and first order desorption, and includes the entrapment of volatile species in more refractory ice mantles. We apply the model to the investigation of the chemistry in hot cores, in which a thick ice mantle built up during the previous cold phase of protostellar evolution undergoes surface reactions and is eventually evaporated. For the first time, the impact of a detailed multilayer approach to grain mantle formation on the warm-up chemistry is explored. The use of a multilayer ice structure has a mixed impact on the abundances of organic species formed during the warm-up phase. For example, the abundance of gaseous HCOOCH3 is lower in the multilayer model than in previous grain models that do not distinguish between layers (so-called "two phase" models). Other gaseous organic species formed in the warm-up phase are affected slightly. Finally, we find that the entrapment of volatile species in water ice can explain the two-jump behavior of H2CO previously found in observations of protostars.

We investigate the molecular evolution and D/H abundance ratios that develop as star formation proceeds from a dense-cloud core to a protostellar core, by solving a gas-grain reaction network applied to a 1-D radiative hydrodynamic model with infalling fluid parcels. Spatial distributions of gas and ice-mantle species are calculated at the first-core stage, and at times after the birth of a protostar. Gas-phase methanol and methane are more abundant than CO at radii $r\lesssim 100$ AU in the first-core stage, but gradually decrease with time, while abundances of larger organic species increase. The warm-up phase, when complex organic molecules are efficiently formed, is longer-lived for those fluid parcels in-falling at later stages. The formation of unsaturated carbon chains (warm carbon-chain chemistry) is also more effective in later stages; C$^+$, which reacts with CH$_4$ to form carbon chains, increases in abundance as the envelope density decreases. The large organic molecules and carbon chains are strongly deuterated, mainly due to high D/H ratios in the parent molecules, determined in the cold phase. We also extend our model to simulate simply the chemistry in circumstellar disks, by suspending the 1-D infall of a fluid parcel at constant disk radii. The species CH$_3$OCH$_3$ and HCOOCH$_3$ increase in abundance in $10^4-10^5$ yr at the fixed warm temperature; both also have high D/H ratios.

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.

We have used the Herschel-HIFI instrument to observe interstellar nitrogen hydrides along the sight-lines towards W49N and G10.6-0.4 in order to elucidate the production pathways leading to nitrogen-bearing species in diffuse gas. All detections show absorption by foreground material over a wide range of velocities, as well as absorption associated directly with the hot-core source itself. As in the previously published observations towards G10.6-0.4, the NH, NH2 and NH3 spectra towards W49N show strikingly similar and non-saturated absorption features. We decompose the absorption of the foreground material towards W49N into different velocity components in order to investigate whether the relative abundances vary among the velocity components, and, in addition, we re-analyse the absorption lines towards G10.6-0.4 in the same manner. Abundances, with respect to molecular hydrogen, in each velocity component are estimated using CH. The analysis points to a co-existence of the nitrogen hydrides in diffuse or translucent interstellar gas with a high molecular fraction. Towards both sources, we find that NH is always at least as abundant as both o-NH2 and o-NH3, in sharp contrast to previous results for dark clouds. We find relatively constant N(NH)/N(o-NH3) and N(o-NH2)/N(o-NH3) ratios with mean values of 3.2 and 1.9 towards W49N, and 5.4 and 2.2 towards G10.6-0.4, respectively. The mean abundance of o-NH3 is ~2x10^-9 towards both sources. The nitrogen hydrides also show linear correlations with CN and HNC towards both sources, and looser correlations with CH. The upper limits on the NH+ abundance indicate column densities < 2 - 14 % of N(NH). Surprisingly low values of the ammonia ortho-to-para ratio are found in both sources, ~0.5 - 0.7 +- 0.1. This result cannot be explained by current models as we had expected to find a value of unity or higher.

(Abridged) We have observed velocity resolved spectra of four ro-vibrational far-infrared transitions of C3 between the vibrational ground state and the low-energy nu2 bending mode at frequencies between 1654--1897 GHz using HIFI on board Herschel, in DR21(OH), a high mass star forming region. Several transitions of CCH and c-C3H2 have also been observed with HIFI and the IRAM 30m telescope. A gas and grain warm-up model was used to identify the primary C3 forming reactions in DR21(OH). We have detected C3 in absorption in four far-infrared transitions, P(4), P(10), Q(2) and Q(4). The continuum sources MM1 and MM2 in DR21(OH) though spatially unresolved, are sufficiently separated in velocity to be identified in the C3 spectra. All C3 transitions are detected from the embedded source MM2 and the surrounding envelope, whereas only Q(4) & P(4) are detected toward the hot core MM1. The abundance of C3 in the envelope and MM2 is \sim6x10^{-10} and \sim3x10^{-9} respectively. For CCH and c-C3H2 we only detect emission from the envelope and MM1. The observed CCH, C3, and c-C3H2 abundances are most consistent with a chemical model with n(H2)\sim5x10^{6} cm^-3 post-warm-up dust temperature, T_max =30 K and a time of \sim0.7-3 Myr. Post warm-up gas phase chemistry of CH4 released from the grain at t\sim 0.2 Myr and lasting for 1 Myr can explain the observed C3 abundance in the envelope of DR21(OH) and no mechanism involving photodestruction of PAH molecules is required. The chemistry in the envelope is similar to the warm carbon chain chemistry (WCCC) found in lukewarm corinos. The observed lower C3 abundance in MM1 as compared to MM2 and the envelope could be indicative of destruction of C3 in the more evolved MM1. The timescale for the chemistry derived for the envelope is consistent with the dynamical timescale of 2 Myr derived for DR21(OH) in other studies.

We simulate the chemistry of infrared dark clouds (IRDCs) with a model in which the physical conditions are homogeneous and time-independent. The chemistry is solved as a function of time with three networks: one purely gas-phase, one that includes accretion and desorption, and one, the complete gas-grain network, that includes surface chemistry in addition. We compare our results with observed molecular abundances for two representative IRDCs -- IRDC013.90-1 and IRDC321.73-1 -- using the molecular species N$_2$H$^+$, HC$_3$N, HNC, HCO$^+$, HCN, C$_2$H, NH$_3$ and CS. IRDC013.90-1 is a cold IRDC, with a temperature below 20 K, while IRDC321.73-1 is somewhat warmer, in the range 20 - 30 K. We find that the complete gas-grain model fits the data very well, but that the goodness-of-fit is not sharply peaked at a particular temperature. Surface processes are important for the explanation of the high gas-phase abundance of N$_2$H$^+$ in IRDC321.73-1. The general success of the 0-D model in reproducing single-dish observations of our limited sample of 8 species shows that it is probably sufficient for an explanation of this type of data. To build and justify more complicated models, including spatial temperature and density structure, contraction, and heating, we require high-resolution interferometric data.

We report the first detection of interstellar mercapto radicals, obtained along the sight-line to the submillimeter continuum source W49N. We have used the GREAT instrument on SOFIA to observe the 1383 GHz Doublet Pi 3/2 J = 5/2 - 3/2 lambda doublet in the upper sideband of the L1 receiver. The resultant spectrum reveals SH absorption in material local to W49N, as well as in foreground gas, unassociated with W49N, that is located along the sight-line. For the foreground material at velocities in the range 37 - 44 km/s with respect to the local standard of rest, we infer a total SH column density ~ 2.6 E+12 cm-2, corresponding to an abundance of ~ 7 E-9 relative to H2, and yielding an SH/H2S abundance ratio ~ 0.13. The observed SH/H2S abundance ratio is much smaller than that predicted by standard models for the production of SH and H2S in turbulent dissipation regions and shocks, and suggests that the endothermic neutral-neutral reaction SH + H2 -> H2S + H must be enhanced along with the ion-neutral reactions believed to produce CH+ and SH+ in diffuse molecular clouds.

We present a novel chemical database for gas-phase astrochemistry. Named the KInetic Database for Astrochemistry (KIDA), this database consists of gas-phase reactions with rate coefficients and uncertainties that will be vetted to the greatest extent possible. Submissions of measured and calculated rate coefficients are welcome, and will be studied by experts before inclusion into the database. Besides providing kinetic information for the interstellar medium, KIDA is planned to contain such data for planetary atmospheres and for circumstellar envelopes. Each year, a subset of the reactions in the database (kida.uva) will be provided as a network for the simulation of the chemistry of dense interstellar clouds with temperatures between 10 K and 300 K. We also provide a code, named Nahoon, to study the time-dependent gas-phase chemistry of 0D and 1D interstellar sources.

An emerging theme in modern astrophysics is the connection between astronomical observations and the underlying physical phenomena that drive our cosmos. Both the mechanisms responsible for the observed astrophysical phenomena and the tools used to probe such phenomena - the radiation and particle spectra we observe - have their roots in atomic, molecular, condensed matter, plasma, nuclear and particle physics. Chemistry is implicitly included in both molecular and condensed matter physics. This connection is the theme of the present report, which provides a broad, though non-exhaustive, overview of progress in our understanding of the cosmos resulting from recent theoretical and experimental advances in what is commonly called laboratory astrophysics. This work, carried out by a diverse community of laboratory astrophysicists, is increasingly important as astrophysics transitions into an era of precise measurement and high fidelity modeling.

In previous warm-up chemical models of the low-mass star-forming region L1527, we investigated the evolution of carbon-chain unsaturated hydrocarbon species when the envelope temperature is slightly elevated to $T\approx 30$ K. These models demonstrated that enhanced abundances of such species can be explained by gas-phase ion-molecule chemistry following the partial sublimation of methane from grain surfaces. We also concluded that the abundances of hydrocarbon radicals such as the C$_{\rm n}$H family should be further enhanced as the temperatures increase to higher values, but this conclusion stood in contrast with the lack of unambiguous detection of these species toward hot core and corino sources. Meanwhile, observational surveys have identified C$_2$H, C$_4$H, CH$_3$CCH, and CH$_3$OH toward hot corinos (especially IRAS 16293-2422) as well as towards L1527, with lower abundances for the carbon chain radicals and higher abundances for the other two species toward the hot corinos. In addition, the {\it Herschel Space Telescope} has detected the bare linear chain C$_3$ in 50 K material surrounding young high-mass stellar objects. To understand these new results, we revisit previous warm-up models with an augmented gas-grain network that incorporated reactions from a gas-phase network that was constructed for use with increased temperature up to 800 K. Some of the newly adopted reactions between carbon-chain species and abundant H$_2$ possess chemical activation energy barriers. The revised model results now better reproduce the observed abundances of unsaturated carbon chains under hot-corino (100 K) conditions and make predictions for the abundances of bare carbon chains in the 50 K regions observed by Herschel HIFI.

We calculate cosmic ray transport with a collisional Boltzmann Transport Equation, including E-M forces. We apply the Crank-Nicholson Method to solve this equation. At each time step, the spatial distribution of cosmic rays is applied to the ZEUS 2D magnetohydrodynamics model, which is then utilized to calculate the resulting E-M field. Finally, the field is applied to the initial equation. This sequence is repeated over many time steps until a steady state is reached. We consider results from t = 0 until steady state for an isotropic low energy cosmic ray flux, and also for an enhanced cosmic ray flux impinging on one side of a molecular cloud. This cosmic ray flux is used to determine an ionization rate of interstellar hydrogen by cosmic rays, zeta. Astrochemical implications are briefly mentioned.

Context: Observations of small carbon-bearing molecules such as CCH, C4H, c-C3H2, and HCO in the Horsehead Nebula have shown these species to have higher abundances towards the edge of the source than towards the center. Aims: Given the determination of a wide range of values for zeta (s-1), the total ionization rate of hydrogen atoms, and the proposal of a column-dependent zeta(N_H), where N_H is the total column of hydrogen nuclei, we desire to determine if the effects of zeta(N_H) in a single object with spatial variation can be observable. We chose the Horsehead Nebula because of its geometry and high density. Method: We model the Horsehead Nebula as a near edge-on photon dominated region (PDR), using several choices for zeta, both constant and as a function of column. The column-dependent zeta functions are determined by a Monte Carlo model of cosmic ray penetration, using a steep power-law spectrum and accounting for ionization and magnetic field effects. We consider a case with low-metal elemental abundances as well as a sulfur-rich case. Results: We show that use of a column-dependent zeta(N_H) of 5(-15) s-1 at the surface and 7.5(-16) s-1 at Av = 10 on balance improves agreement between measured and theoretical molecular abundances, compared with constant values of zeta.

We have examined methanol emission from Orion KL with of the {\em Herschel}/HIFI instrument, and detected two methanol bands centered at 524 GHz and 1061 GHz. The 524 GHz methanol band (observed in HIFI band 1a) is dominated by the isolated $\Delta$J$=$0, K$=-4\rightarrow$-3, v$_t$$=$0 Q branch, and includes 25 E-type and 2 A-type transitions. The 1061 GHz methanol band (observed in HIFI band 4b) is dominated by the $\Delta$J$=$0, K$=7\rightarrow$6, v$_t$$=$0 Q branch transitions which are mostly blended. We have used the isolated E-type v$_t$$=$0 methanol transitions to explore the physical conditions in the molecular gas. With HIFI's high velocity resolution, the methanol emission contributed by different spatial components along the line of sight toward Orion KL (hot core, low velocity flow, and compact ridge) can be distinguished and studied separately. The isolated transitions detected in these bands cover a broad energy range (upper state energy ranging from 80 K to 900 K), which provides a unique probe of the thermal structure in each spatial component. The observations further show that the compact ridge is externally heated. These observations demonstrate the power of methanol lines as probes of the physical conditions in warm regions in close proximity to young stars.

We survey the current situation regarding chemical modelling of the synthesis of molecules in the interstellar medium. The present state of knowledge concerning the rate coefficients and their uncertainties for the major gas-phase processes -- ion-neutral reactions, neutral-neutral reactions, radiative association, and dissociative recombination -- is reviewed. Emphasis is placed on those reactions that have been identified, by sensitivity analyses, as 'crucial' in determining the predicted abundances of the species observed in the interstellar medium. These sensitivity analyses have been carried out for gas-phase models of three representative, molecule-rich, astronomical sources: the cold dense molecular clouds TMC-1 and L134N, and the expanding circumstellar envelope IRC +10216. Our review has led to the proposal of new values and uncertainties for the rate coefficients of many of the key reactions. The impact of these new data on the predicted abundances in TMC-1 and L134N is reported. Interstellar dust particles also influence the observed abundances of molecules in the interstellar medium. Their role is included in gas-grain, as distinct from gas-phase only, models. We review the methods for incorporating both accretion onto, and reactions on, the surfaces of grains in such models, as well as describing some recent experimental efforts to simulate and examine relevant processes in the laboratory. These efforts include experiments on the surface-catalysed recombination of hydrogen atoms, on chemical processing on and in the ices that are known to exist on the surface of interstellar grains, and on desorption processes, which may enable species formed on grains to return to the gas-phase.

We report observations of the reactive molecular ions OH$^+$, H$_2$O$^+$, and H$_3$O$^+$ towards Orion KL with Herschel/HIFI. All three $N=1-0$ fine-structure transitions of OH$^+$ at 909, 971, and 1033GHz and both fine-structure components of the doublet {\it ortho}-H$_2$O$^+$ $1_{11}-0_{00}$ transition at 1115 and 1139GHz were detected; an upper limit was obtained for H$_3$O$^+$. OH$^+$ and H$_2$O$^+$ are observed purely in absorption, showing a narrow component at the source velocity of 9 kms$^{-1}$, and a broad blueshifted absorption similar to that reported recently for HF and {\it para}-H$_{2}^{18}$O, and attributed to the low velocity outflow of Orion KL. We estimate column densities of OH$^+$ and H$_2$O$^+$ for the 9 km s$^{-1}$ component of $9 \pm 3 \times 10^{12}$cm$^{-2}$ and $7 \pm 2 \times 10^{12}$cm$^{-2}$, and those in the outflow of $1.9 \pm 0.7 \times 10^{13}$cm$^{-2}$ and $1.0 \pm 0.3 \times 10^{13}$cm$^{-2}$. Upper limits of $2.4\times 10^{12}$cm$^{-2}$ and $8.7\times 10^{12}$cm$^{-2}$ were derived for the column densities of {\it ortho} and {\it para}-H$_3$O$^+$ from transitions near 985 and 1657GHz. The column densities of the three ions are up to an order of magnitude lower than those obtained from recent observations of W31C and W49N. The comparatively low column densities may be explained by a higher gas density despite the assumption of a very high ionization rate.

In contrast with numerous studies on the physical and chemical structure of low- and high-mass protostars, much less is known about their intermediate-mass counterparts, a class of objects that could help to elucidate the mechanisms of star formation on both ends of the mass range. We present the first results from a rich HIFI spectral dataset on an intermediate-mass protostar, OMC2-FIR4, obtained in the CHESS (Chemical HErschel SurveyS of star forming regions) key programme. The more than 100 methanol lines detected between 554 and 961 GHz cover a range in upper level energy of 40 to 540 K. Our physical interpretation focusses on the hot core, but likely the cold envelope and shocked regions also play a role in reality, because an analysis of the line profiles suggests the presence of multiple emission components. An upper limit of 10^-6 is placed on the methanol abundance in the hot core, using a population diagram, large-scale source model and other considerations. This value is consistent with abundances previously seen in low-mass hot cores. Furthermore, the highest energy lines at the highest frequencies display asymmetric profiles, which may arise from infall around the hot core.

We report the detection of more than 48 velocity-resolved ground rotational state transitions of H2(16)O, H2(18)O, and H2(17)O - most for the first time - in both emission and absorption toward Orion KL using Herschel/HIFI. We show that a simple fit, constrained to match the known emission and absorption components along the line of sight, is in excellent agreement with the spectral profiles of all the water lines. Using the measured H2(18)O line fluxes, which are less affected by line opacity than their H2(16)O counterparts, and an escape probability method, the column densities of H2(18)O associated with each emission component are derived. We infer total water abundances of 7.4E-5, 1.0E-5, and 1.6E-5 for the plateau, hot core, and extended warm gas, respectively. In the case of the plateau, this value is consistent with previous measures of the Orion-KL water abundance as well as those of other molecular outflows. In the case of the hot core and extended warm gas, these values are somewhat higher than water abundances derived for other quiescent clouds, suggesting that these regions are likely experiencing enhanced water-ice sublimation from (and reduced freeze-out onto) grain surfaces due to the warmer dust in these sources.

Observations of HDO are an important complement for studies of water, because they give strong constraints on the formation processes -- grain surfaces versus energetic process in the gas phase, e.g. in shocks. The HIFI observations of multiple transitions of HDO in Sgr~B2(M) presented here allow the determination of the HDO abundance throughout the envelope, which has not been possible before with ground-based observations only. The abundance structure has been modeled with the spherical Monte Carlo radiative transfer code RATRAN, which also takes radiative pumping by continuum emission from dust into account. The modeling reveals that the abundance of HDO rises steeply with temperature from a low abundance ($2.5\times 10^{-11}$) in the outer envelope at temperatures below 100~K through a medium abundance ($1.5\times 10^{-9}$) in the inner envelope/outer core, at temperatures between 100 and 200~K, and finally a high abundance ($3.5\times 10^{-9}$) at temperatures above 200~K in the hot core.

To investigate the accretion and feedback processes in massive star formation, we analyze the shapes of emission lines from hot molecular cores, whose asymmetries trace infall and expansion motions. The high-mass star forming region SgrB2(M) was observed with Herschel/HIFI (HEXOS key project) in various lines of HCN and its isotopologues, complemented by APEX data. The observations are compared to spherically symmetric, centrally heated models with density power-law gradient and different velocity fields (infall or infall+expansion), using the radiative transfer code RATRAN. The HCN line profiles are asymmetric, with the emission peak shifting from blue to red with increasing J and decreasing line opacity (HCN to H$^{13}$CN). This is most evident in the HCN 12--11 line at 1062 GHz. These line shapes are reproduced by a model whose velocity field changes from infall in the outer part to expansion in the inner part. The qualitative reproduction of the HCN lines suggests that infall dominates in the colder, outer regions, but expansion dominates in the warmer, inner regions. We are thus witnessing the onset of feedback in massive star formation, starting to reverse the infall and finally disrupting the whole molecular cloud. To obtain our result, the THz lines uniquely covered by HIFI were critically important.

In contrast to extensively studied dense star-forming cores, little is known about diffuse gas surrounding star-forming regions. We study molecular gas in the high-mass star-forming region NGC6334I, which contains diffuse, quiescent components that are inconspicuous in widely used molecular tracers such as CO. We present Herschel/HIFI observations of CH toward NGC6334I observed as part of the CHESS key program. HIFI resolves the hyperfine components of its J=3/2-1/2 transition, observed in both emission and absorption. The CH emission appears close to the systemic velocity of NGC6334I, while its measured linewidth of 3 km/s is smaller than previously observed in dense gas tracers such as NH3 and SiO. The CH abundance in the hot core is 7 10^-11, two to three orders of magnitude lower than in diffuse clouds. While other studies find distinct outflows in, e.g., CO and H2O toward NGC6334I, we do not detect outflow signatures in CH. To explain the absorption signatures, at least two absorbing components are needed at -3.0 and +6.5 km/s with N(CH)=7 10^13 and 3 10^13 cm^-2. Two additional absorbing clouds are found at +8.0 and 0.0 km/s, both with N(CH)=2 10^13 cm^-2. Turbulent linewidths for the four absorption components vary between 1.5 and 5.0 km/s in FWHM. We constrain physical properties of our CH clouds by matching our CH absorbers with other absorption signatures. In the hot core, molecules such as H2O and CO trace gas that is heated and dynamically influenced by outflow activity, whereas CH traces more quiescent material. The four CH absorbers have column densities and turbulent properties consistent with diffuse clouds: two are located near NGC6334, and two are unrelated foreground clouds. Local density and dynamical effects influence the chemical composition of physical components of NGC6334, causing some components to be seen in CH but not in other tracers, and vice versa.

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

We present the first high spectral resolution observations of Orion KL in the frequency ranges 1573.4 - 1702.8 GHz (band 6b) and 1788.4 - 1906.8 GHz (band 7b) obtained using the HIFI instrument on board the Herschel Space Observatory. We characterize the main emission lines found in the spectrum, which primarily arise from a range of components associated with Orion KL including the hot core, but also see widespread emission from components associated with molecular outflows traced by H2O, SO2, and OH. We find that the density of observed emission lines is significantly diminished in these bands compared to lower frequency Herschel/HIFI bands.

Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star forming regions, and in particular the Class 0 protostar IRAS 16293-2422. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe the fundamental 1,1,1 - 0,0,0 transition of the ortho-D2O molecule, inaccessible from the ground, and to determine the ortho-to-para D2O ratio. We have detected the fundamental transition of the ortho-D2O molecule at 607.35 GHz towards IRAS 16293-2422. The line is seen in absorption with a line opacity of 0.62 +/- 0.11 (1 sigma). From the previous ground-based observations of the fundamental 1,1,0 - 1,0,1 transition of para-D2O seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 +/- 0.05 (1 sigma). We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 sigma level of uncertainty, to be compared with the thermal equilibrium value of 2:1.

We present observations of twelve rotational transitions of H2O-16, H2O-18, and H2O-17 toward the massive star-forming region NGC 6334 I, carried out with Herschel/HIFI as part of the guaranteed time key program Chemical HErschel Surveys of Star forming regions (CHESS). We analyze these observations to obtain insights into physical processes in this region. We identify three main gas components (hot core, cold foreground, and outflow) in NGC 6334 I and derive the physical conditions in these components. The hot core, identified by the emission in highly excited lines, shows a high excitation temperature of 200 K, whereas water in the foreground component is predominantly in the ortho- and para- ground states. The abundance of water varies between 4 10^-5 (outflow) and 10^-8 (cold foreground gas). This variation is most likely due to the freeze-out of water molecules onto dust grains. The H2O-18/H2O-17 abundance ratio is 3.2, which is consistent with the O-18/O-17 ratio determined from CO isotopologues. The ortho/para ratio in water appears to be relatively low 1.6(1) in the cold, quiescent gas, but close to the equilibrium value of three in the warmer outflow material (2.5(0.8)).

The HIFI instrument on board the Herschel Space Observatory has been used to observe interstellar nitrogen hydrides along the sight-line towards G10.6-0.4 in order to improve our understanding of the interstellar chemistry of nitrogen. We report observations of absorption in NH N=1-0, J=2-1 and ortho-NH2 1_1,1-0_0,0. We also observed ortho-NH3 1_0-0_0, and 2_0-1_0, para-NH3 2_1-1_1, and searched unsuccessfully for NH+. All detections show emission and absorption associated directly with the hot-core source itself as well as absorption by foreground material over a wide range of velocities. All spectra show similar, non-saturated, absorption features, which we attribute to diffuse molecular gas. Total column densities over the velocity range 11-54 km/s are estimated. The similar profiles suggest fairly uniform abundances relative to hydrogen, approximately 6*10^-9, 3*10^-9, and 3*10^-9 for NH, NH2, and NH3, respectively. These abundances are discussed with reference to models of gas-phase and surface chemistry.

We report a detection of the fundamental rotational transition of hydrogen fluoride in absorption towards Orion KL using Herschel/HIFI. After the removal of contaminating features associated with common molecules ("weeds"), the HF spectrum shows a P-Cygni profile, with weak redshifted emission and strong blue-shifted absorption, associated with the low-velocity molecular outflow. We derive an estimate of 2.9 x 10^13 cm^-2 for the HF column density responsible for the broad absorption component. Using our best estimate of the H2 column density within the low-velocity molecular outflow, we obtain a lower limit of ~1.6 x 10^-10 for the HF abundance relative to hydrogen nuclei, corresponding to 0.6% of the solar abundance of fluorine. This value is close to that inferred from previous ISO observations of HF J=2--1 absorption towards Sgr B2, but is in sharp contrast to the lower limit of 6 x 10^-9 derived by Neufeld et al. (2010) for cold, foreground clouds on the line of sight towards G10.6-0.4.

We discuss the detection of absorption by interstellar hydrogen fluoride (HF) along the sight line to the submillimeter continuum sources W49N and W51. We have used Herschel's HIFI instrument in dual beam switch mode to observe the 1232.4762 GHz J = 1 - 0 HF transition in the upper sideband of the band 5a receiver. We detected foreground absorption by HF toward both sources over a wide range of velocities. Optically thin absorption components were detected on both sight lines, allowing us to measure - as opposed to obtain a lower limit on - the column density of HF for the first time. As in previous observations of HF toward the source G10.6-0.4, the derived HF column density is typically comparable to that of water vapor, even though the elemental abundance of oxygen is greater than that of fluorine by four orders of magnitude. We used the rather uncertain N(CH)-N(H2) relationship derived previously toward diffuse molecular clouds to infer the molecular hydrogen column density in the clouds exhibiting HF absorption. Within the uncertainties, we find that the abundance of HF with respect to H2 is consistent with the theoretical prediction that HF is the main reservoir of gas-phase fluorine for these clouds. Thus, hydrogen fluoride has the potential to become an excellent tracer of molecular hydrogen, and provides a sensitive probe of clouds of small H2 column density. Indeed, the observations of hydrogen fluoride reported here reveal the presence of a low column density diffuse molecular cloud along the W51 sight line, at an LSR velocity of ~ 24kms-1, that had not been identified in molecular absorption line studies prior to the launch of Herschel.

We present initial results from the Herschel GT key program: Herschel observations of EXtra-Ordinary Sources (HEXOS) and outline the promise and potential of spectral surveys with Herschel/HIFI. The HIFI instrument offers unprecedented sensitivity, as well as continuous spectral coverage across the gaps imposed by the atmosphere, opening up a largely unexplored wavelength regime to high-resolution spectroscopy. We show the spectrum of Orion KL between 480 and 560 GHz and from 1.06 to 1.115 THz. From these data, we confirm that HIFI separately measures the dust continuum and spectrally resolves emission lines in Orion KL. Based on this capability we demonstrate that the line contribution to the broad-band continuum in this molecule-rich source is ~20-40% below 1 THz and declines to a few percent at higher frequencies. We also tentatively identify multiple transitions of HD18O in the spectra. The first detection of this rare isotopologue in the interstellar medium suggests that HDO emission is optically thick in the Orion hot core with HDO/H2O ~ 0.02. We discuss the implications of this detection for the water D/H ratio in hot cores.

We have observed CH absorption lines ($J=3/2, N=1 \leftarrow J=1/2, N=1$) against the continuum source Sgr~B2(M) using the \textit{Herschel}/HIFI instrument. With the high spectral resolution and wide velocity coverage provided by HIFI, 31 CH absorption features with different radial velocities and line widths are detected and identified. The narrower line width and lower column density clouds show `spiral arm' cloud characteristics, while the absorption component with the broadest line width and highest column density corresponds to the gas from the Sgr~B2 envelope. The observations show that each `spiral arm' harbors multiple velocity components, indicating that the clouds are not uniform and that they have internal structure. This line-of-sight through almost the entire Galaxy offers unique possibilities to study the basic chemistry of simple molecules in diffuse clouds, as a variety of different cloud classes are sampled simultaneously. We find that the linear relationship between CH and H$_2$ column densities found at lower $A_V$ by UV observations does not continue into the range of higher visual extinction. There, the curve flattens, which probably means that CH is depleted in the denser cores of these clouds.

We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S). The H$_2$Cl$^+$ column densities are compared to those of the chemically-related species HCl. The derived HCl/H$_2$Cl$^+$ column density ratios, $\sim$1--10, are within the range predicted by models of diffuse and dense Photon Dominated Regions (PDRs). However, the observed H$_2$Cl$^+$ column densities, in excess of $10^{13}$~cm$^{-2}$, are significantly higher than the model predictions. Our observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints for astrochemical models.

We present Herschel/HIFI observations of the fundamental rotational transitions of ortho- and para-H$_2^{16}$O and H$_2^{18}$O in absorption towards Sagittarius~B2(M) and W31C. The ortho/para ratio in water in the foreground clouds on the line of sight towards these bright continuum sources is generally consistent with the statistical high-temperature ratio of 3, within the observational uncertainties. However, somewhat unexpectedly, we derive a low ortho/para ratio of $2.35 \pm 0.35$, corresponding to a spin temperature of $\sim$27~K, towards Sagittarius~B2(M) at velocities of the expanding molecular ring. Water molecules in this region appear to have formed with, or relaxed to, an ortho/para ratio close to the value corresponding to the local temperature of the gas and dust.

We report the detection of absorption by interstellar hydroxyl cations and water cations, along the sight-line to the bright continuum source W49N. We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 972 GHz N = 1 - 0 transition of OH+ and the 1115 GHz 1(11) - 0(00) transition of ortho-H2O+. The resultant spectra show absorption by ortho-H2O+, and strong absorption by OH+, in foreground material at velocities in the range 0 to 70 km/s with respect to the local standard of rest. The inferred OH+/H2O+ abundance ratio ranges from ~ 3 to ~ 15, implying that the observed OH+ arises in clouds of small molecular fraction, in the 2 - 8% range. This conclusion is confirmed by the distribution of OH+ and H2O+ in Doppler velocity space, which is similar to that of atomic hydrogen, as observed by means of 21 cm absorption measurements, and dissimilar from that typical of other molecular tracers. The observed OH+/H abundance ratio of a few E-8 suggests a cosmic ray ionization rate for atomic hydrogen of (0.6 - 2.4) E-16 s-1, in good agreement with estimates inferred previously for diffuse clouds in the Galactic disk from observations of interstellar H3+ and other species.

H2O+ has been observed in its ortho- and para- states toward the massive star forming core Sgr B2(M), located close to the Galactic center. The observations show absorption in all spiral arm clouds between the Sun and Sgr B2. The average o/p ratio of H2O+ in most velocity intervals is 4.8, which corresponds to a nuclear spin temperature of 21 K. The relationship of this spin temperature to the formation temperature and current physical temperature of the gas hosting H2O+ is discussed, but no firm conclusion is reached. In the velocity interval 0 to 60 km/s, an ortho/para ratio of below unity is found, but if this is due to an artifact of contamination by other species or real is not clear.

The outflow driven by the low-mass class 0 protostar L1157 is the prototype of the so-called chemically active outflows. The bright bowshock B1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and chemical processes are strongly interdependent. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key program "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the excitation conditions in the shock region. The CO 5-4 and H2O lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 1b spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window help establish the origin of the emission detected in the main-beam of HIFI, and the physical conditions in the shock.}{Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission arises from two regions with distinct physical conditions: an extended, warm (100K), dense (3e5 cm-3) component at low-velocity, which dominates the water line flux in Band~1; a secondary component in a small region of B1 (a few arcsec) associated with high-velocity, hot (> 400 K) gas of moderate density ((1.0-3.0)e4 cm-3), which appears to dominate the flux of the water line at 179mu observed with PACS. The water abundance is enhanced by two orders of magnitude between the low- and the high-velocity component, from 8e-7 up to 8e-5. The properties of the high-velocity component agree well with the predictions of steady-state C-shock models.

We report the detection of absorption lines by the reactive ions OH+, H2O+ and H3O+ along the line of sight to the submillimeter continuum source G10.6$-$0.4 (W31C). We used the Herschel HIFI instrument in dual beam switch mode to observe the ground state rotational transitions of OH+ at 971 GHz, H2O+ at 1115 and 607 GHz, and H3O+ at 984 GHz. The resultant spectra show deep absorption over a broad velocity range that originates in the interstellar matter along the line of sight to G10.6$-$0.4 as well as in the molecular gas directly associated with that source. The OH+ spectrum reaches saturation over most velocities corresponding to the foreground gas, while the opacity of the H2O+ lines remains lower than 1 in the same velocity range, and the H3O+ line shows only weak absorption. For LSR velocities between 7 and 50 kms$^{-1}$ we estimate total column densities of $N$(OH+) $> 2.5 \times 10^{14}$ cm$^{-2}$, $N$(H2O+) $\sim 6 \times 10^{13}$ cm$^{-2}$ and $N$(H3O+) $\sim 4.0 \times 10^{13}$ cm$^{-2}$. These detections confirm the role of O$^+$ and OH$^+$ in initiating the oxygen chemistry in diffuse molecular gas and strengthen our understanding of the gas phase production of water. The high ratio of the OH+ by the H2O+ column density implies that these species predominantly trace low-density gas with a small fraction of hydrogen in molecular form.

We report the detection of strong absorption by interstellar hydrogen fluoride along the sight-line to the submillimeter continuum source G10.6-0.4 (W31C). We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 1232.4763 GHz J=1-0 HF transition in the upper sideband of the Band 5a receiver. The resultant spectrum shows weak HF emission from G10.6-0.4 at LSR velocities in the range -10 to -3 km/s, accompanied by strong absorption by foreground material at LSR velocities in the range 15 to 50 km/s. The spectrum is similar to that of the 1113.3430 GHz 1(11)-0(00) transition of para-water, although at some frequencies the HF (hydrogen fluoride) optical depth clearly exceeds that of para-H2O. The optically-thick HF absorption that we have observed places a conservative lower limit of 1.6E+14 cm-2 on the HF column density along the sight-line to G10.6-0.4. Our lower limit on the HF abundance, 6E-9 relative to hydrogen nuclei, implies that hydrogen fluoride accounts for between ~ 30 and 100% of the fluorine nuclei in the gas phase along this sight-line. This observation corroborates theoretical predictions that - because the unique thermochemistry of fluorine permits the exothermic reaction of F atoms with molecular hydrogen - HF will be the dominant reservoir of interstellar fluorine under a wide range of conditions.

We identify a prominent absorption feature at 1115 GHz, detected in first HIFI spectra towards high-mass star-forming regions, and interpret its astrophysical origin. The characteristic hyperfine pattern of the H2O+ ground-state rotational transition, and the lack of other known low-energy transitions in this frequency range, identifies the feature as H2O+ absorption against the dust continuum background and allows us to derive the velocity profile of the absorbing gas. By comparing this velocity profile with velocity profiles of other tracers in the DR21 star-forming region, we constrain the frequency of the transition and the conditions for its formation. In DR21, the velocity distribution of H2O+ matches that of the [CII] line at 158\mu\m and of OH cm-wave absorption, both stemming from the hot and dense clump surfaces facing the HII-region and dynamically affected by the blister outflow. Diffuse foreground gas dominates the absorption towards Sgr B2. The integrated intensity of the absorption line allows us to derive lower limits to the H2O+ column density of 7.2e12 cm^-2 in NGC 6334, 2.3e13 cm^-2 in DR21, and 1.1e15 cm^-2 in Sgr B2.

Because of new telescopes that will dramatically improve our knowledge of the interstellar medium, chemical models will have to be used to simulate the chemistry of many regions with diverse properties. To make these models more robust, it is important to understand their sensitivity to a variety of parameters. In this article, we report a study of the sensitivity of a chemical model of a cold dense core, with homogeneous and time-independent physical conditions, to variations in the following parameters: initial chemical inventory, gas temperature and density, cosmic-ray ionization rate, chemical reaction rate coefficients, and elemental abundances. From the results of the parameter variations, we can quantify the sensitivity of the model to each parameter as a function of time. Our results can be used in principle with observations to constrain some parameters for different cold clouds. We also attempted to use the Monte Carlo approach with all parameters varied collectively. Within the parameter ranges studied, the most critical parameters turn out to be the reaction rate coefficients at times up to 4e5 yr and elemental abundances at later times. At typical times of best agreement with observation, models are sensitive to both of these parameters. The models are less sensitive to other parameters such as the gas density and temperature. The improvement of models will require that the uncertainties in rate coefficients of important reactions be reduced. As the chemistry becomes better understood and more robust, it should be possible to use model sensitivities concerning other parameters, such as the elemental abundances and the cosmic ray ionization rate, to yield detailed information on cloud properties and history. Nevertheless, at the current stage, we cannot determine the best values of all the parameters simultaneously based on purely observational constraints.

Context: Chemical models of dense cloud cores often utilize the so-called pseudo-time-dependent approximation, in which the physical conditions are held fixed and uniform as the chemistry occurs. In this approximation, the initial abundances chosen, which are totally atomic in nature except for molecular hydrogen, are artificial. A more detailed approach to the chemistry of dense cold cores should include the physical evolution during their early stages of formation. Aims: Our major goal is to investigate the initial synthesis of molecular ices and gas-phase molecules as cold molecular gas begins to form behind a shock in the diffuse interstellar medium. The abundances calculated as the conditions evolve can then be utilized as reasonable initial conditions for a theory of the chemistry of dense cores. Methods: Hydrodynamic shock-wave simulations of the early stages of cold core formation are used to determine the time-dependent physical conditions for a gas-grain chemical network. We follow the cold post-shock molecular evolution of ices and gas-phase molecules for a range of visual extinction up to AV ~ 3, which increases with time. At higher extinction, self-gravity becomes important. Results: As the newly condensed gas enters its cool post-shock phase, a large amount of CO is produced in the gas. As the CO forms, water ice is produced on grains, while accretion of CO produces CO ice. The production of CO2 ice from CO occurs via several surface mechanisms, while the production of CH4 ice is slowed by gas-phase conversion of C into CO.

Methanol and its precursor formaldehyde are among the most studied organic molecules in the interstellar medium and are abundant in the gaseous and solid phases. We recently developed a model to simulate CO hydrogenation via H atoms on interstellar ice surfaces, the most important interstellar route to H2CO and CH3OH, under laboratory conditions. We extend this model to simulate the formation of both organic species under interstellar conditions, including freeze-out from the gas and hydrogenation on surfaces. Our aim is to compare calculated abundance ratios with observed values and with the results of prior models. Simulations under different conditions, including density and temperature, have been performed. We find that H2CO and CH3OH form efficiently in cold dense cores or the cold outer envelopes of young stellar objects. The grain mantle is found to have a layered structure with CH3OH on top. The species CO and H2CO are found to exist predominantly in the lower layers of ice mantles where they are not available for hydrogenation at late times.

The WGLA of the AAS (http://www.aas.org/labastro/) promotes collaboration and exchange of knowledge between astronomy and planetary sciences and the laboratory sciences (physics, chemistry, and biology). Laboratory data needs of ongoing and next generation planetary science missions are carefully evaluated and recommended in this white paper submitted by the WGLA to Planetary Decadal Survey.

We have investigated the role of molecular anion chemistry in pseudo-time dependent chemical models of dark clouds. With oxygen-rich elemental abundances, the addition of anions results in a slight improvement in the overall agreement between model results and observations of molecular abundances in TMC-1 (CP). More importantly, with the inclusion of anions, we see an enhanced production efficiency of unsaturated carbon-chain neutral molecules, especially in the longer members of the families CnH, CnH2, and HCnN. The use of carbon-rich elemental abundances in models of TMC-1 (CP) with anion chemistry worsens the agreement with observations obtained in the absence of anions.

Laboratory astrophysics and complementary theoretical calculations are the foundations of astronomy and astrophysics and will remain so into the foreseeable future. The mission enabling impact of laboratory astrophysics ranges from the scientific conception stage for airborne and space-based observatories, all the way through to the scientific return of these missions. It is our understanding of the under-lying physical processes and the measurements of critical physical parameters that allows us to address fundamental questions in astronomy and astrophysics. In this regard, laboratory astrophysics is much like detector and instrument development at NASA. These efforts are necessary for the success of astronomical research being funded by NASA. Without concomitant efforts in all three directions (observational facilities, detector/instrument development, and laboratory astrophysics) the future progress of astronomy and astrophysics is imperiled. In addition, new developments in experimental technologies have allowed laboratory studies to take on a new role as some questions which previously could only be studied theoretically can now be addressed directly in the lab. With this in mind we, the members of the AAS Working Group on Laboratory Astrophysics (WGLA), have prepared this White Paper on the laboratory astrophysics infrastructure needed to maximize the scientific return from NASA's space and Earth sciences program.

Laboratory astrophysics and complementary theoretical calculations are the foundations of astronomy and astrophysics and will remain so into the foreseeable future. The impact of laboratory astrophysics ranges from the scientific conception stage for ground-based, airborne, and space-based observatories, all the way through to the scientific return of these projects and missions. It is our understanding of the under-lying physical processes and the measurements of critical physical parameters that allows us to address fundamental questions in astronomy and astrophysics. In this regard, laboratory astrophysics is much like detector and instrument development at NASA, NSF, and DOE. These efforts are necessary for the success of astronomical research being funded by the agencies. Without concomitant efforts in all three directions (observational facilities, detector/instrument development, and laboratory astrophysics) the future progress of astronomy and astrophysics is imperiled. In addition, new developments in experimental technologies have allowed laboratory studies to take on a new role as some questions which previously could only be studied theoretically can now be addressed directly in the lab. With this in mind we, the members of the AAS Working Group on Laboratory Astrophysics, have prepared this State of the Profession Position Paper on the laboratory astrophysics infrastructure needed to ensure the advancement of astronomy and astrophysics in the next decade.

As the panel on Planetary Systems and Star Formation (PSF) is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of solar system bodies (other than the Sun) and extrasolar planets, debris disks, exobiology, the formation of individual stars, protostellar and protoplanetary disks, molecular clouds and the cold ISM, dust, and astrochemistry. Central to the progress in these areas are the corresponding advances in laboratory astro- physics which are required for fully realizing the PSF scientific opportunities in the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics and chemistry which produce the observed spectra and describe the astrophysical processes. We discuss four areas of laboratory astrophysics relevant to the PSF panel: atomic, molecular, solid matter, and plasma physics. Section 2 describes some of the new opportunities and compelling themes which will be enabled by advances in laboratory astrophysics. Section 3 provides the scientific context for these opportunities. Section 4 discusses some experimental and theoretical advances in laboratory astrophysics required to realize the PSF scientific opportunities of the next decade. As requested in the Call for White Papers, we present in Section 5 four central questions and one area with unusual discovery potential. We give a short postlude in Section 6.

As the Galaxies across Cosmic Time (GCT) panel is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of the formation, evolution, and global properties of galaxies and galaxy clusters, as well as active galactic nuclei and QSOs, mergers, star formation rate, gas accretion, and supermassive black holes. Central to the progress in these areas are the corresponding advances in laboratory astrophysics that are required for fully realizing the GCT scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics that produce the observed astrophysical processes. The 5 areas of laboratory astrophysics that we have identified as relevant to the CFP panel are atomic, molecular, solid matter, plasma, nuclear, and particle physics. In this white paper, we describe in Section 2 some of the new scientific opportunities and compelling scientific themes that will be enabled by advances in laboratory astrophysics. In Section 3, we provide the scientific context for these opportunities. Section 4 briefly discusses some of the experimental and theoretical advances in laboratory astrophysics required to realize the GCT scientific opportunities of the next decade. As requested in the Call for White Papers, Section 5 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 6.

As the Galactic Neighborhood (GAN) panel is fully aware, the next decade will see major advances in our understanding of this area of research. To quote from their charge, these advances will occur in studies of the galactic neighborhood, including the structure and properties of the Milky Way and nearby galaxies, and their stellar populations and evolution, as well as interstellar media and star clusters. Central to the progress in these areas are the corresponding advances in laboratory astrophysics that are required for fully realizing the GAN scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics and chemistry that produces the observed astrophysical processes. The 5 areas of laboratory astrophysics that we have identified as relevant to the GAN panel are atomic, molecular, solid matter, plasma, and nuclear physics. In this white paper, we describe in Section 2 some of the new scientific opportunities and compelling scientific themes that will be enabled by advances in laboratory astrophysics. In Section 3, we provide the scientific context for these opportunities. Section 4 briefly discusses some of the experimental and theoretical advances in laboratory astrophysics required to realize the GAN scientific opportunities of the next decade. As requested in the Call for White Papers, Section 5 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 6.

As the Cosmology and Fundamental Physics (CFP) panel is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of the early universe, the microwave background, the reionization and galaxy formation up to virialization of protogalaxies, large scale structure, the intergalactic medium, the determination of cosmological parameters, dark matter, dark energy, tests of gravity, astronomically determined physical constants, and high energy physics using astronomical messengers. Central to the progress in these areas are the corresponding advances in laboratory astrophysics which are required for fully realizing the CFP scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics which produce the observed astrophysical processes. The 5 areas of laboratory astrophysics which we have identified as relevant to the CFP panel are atomic, molecular, plasma, nuclear, and particle physics. Here, Section 2 describes some of the new scientific opportunities and compelling scientific themes which will be enabled by advances in laboratory astrophysics. In Section 3, we provide the scientific context for these opportunities. Section 4 briefly discusses some of the experimental and theoretical advances in laboratory astrophysics required to realize the CFP scientific opportunities of the next decade. As requested in the Call for White Papers, Section 5 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 6.

As the Stars and Stellar Evolution (SSE) panel is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of the Sun as a star, stellar astrophysics, the structure and evolution of single and multiple stars, compact objects, SNe, gamma-ray bursts, solar neutrinos, and extreme physics on stellar scales. Central to the progress in these areas are the corresponding advances in laboratory astrophysics, required to fully realize the SSE scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics that produces the observed astrophysical processes. The 6 areas of laboratory astrophysics, which we have identified as relevant to the CFP panel, are atomic, molecular, solid matter, plasma, nuclear physics, and particle physics. In this white paper, we describe in Section 2 the scientific context and some of the new scientific opportunities and compelling scientific themes which will be enabled by advances in laboratory astrophysics. In Section 3, we discuss some of the experimental and theoretical advances in laboratory astrophysics required to realize the SSE scientific opportunities of the next decade. As requested in the Call for White Papers, Section 4 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 5.

Chemical networks used for models of interstellar clouds contain many reactions, some of them with poorly determined rate coefficients and/or products. In this work, we report a method for improving the predictions of molecular abundances using sensitivity methods and ab initio calculations. Based on the chemical network osu.2003, we used two different sensitivity methods to determine the most important reactions as a function of time for models of dense cold clouds. Of these reactions, we concentrated on those between C and C3 and between C and C5, both for their effect on specific important species such as CO and for their general effect on large numbers of species. We then used ab initio and kinetic methods to determine an improved rate coefficient for the former reaction and a new set of products, plus a slightly changed rate coefficient for the latter. Putting our new results in a pseudo-time-dependent model of cold dense clouds, we found that the abundances of many species are altered at early times, based on large changes in the abundances of CO and atomic C. We compared the effect of these new rate coefficients/products on the comparison with observed abundances and found that they shift the best agreement from 3e4 yr to (1-3)e5 yr.

Aims. We attempt to understand the presence of gas phase CO below its freezing temperature in circumstellar disks. We study two promising mechanisms to explain this phenomenon: turbulent mixing and photodesorption. Methods. We compute the chemical evolution of circumstellar disks including grain surface reactions with and without turbulent mixing and CO photodesorption. Results. We show that photodesorption significantly enhances the gas phase CO abundance, by extracting CO from the grains when the visual extinction remains below about 5 magnitudes. However the resulting dependence of column density on radial distance is not consistent with observations so far. We propose that this inconsistency could be the result of grain growth. On the other hand, the influence of turbulent mixing is not found to be straightforward. The efficiency of turbulent mixing depends upon a variety of parameters, including the disk structure. For the set of parameters we chose, turbulent mixing is not found to have any significant influence on the CO column density.

We present a study of the chemistry of a dense photon-dominated region (PDR) using a time-dependent chemical model. Our major interest is to study the spatial distribution of complex molecules such as hydrocarbons and cyanopolyynes in the cool dense material bordering regions where star formation has taken place. Our standard model uses a homogeneous cloud of density 2x10e4 cm-3 and temperature T=40 K, which is irradiated by a far-ultraviolet radiation field of intermediate intensity, given by X=100. We find that over a range of times unsaturated hydrocarbons (e.g., C2H, C4H, C3H2) have relatively high fractional abundances in the more external layers of the PDR, whereas their abundances in the innermost layers are several orders of magnitudes lower. On the other hand, molecules that are typical of late-time chemistry are usually more abundant in the inner parts of the PDR. We also present results for models with different density, temperature, intensity of the radiation field and initial fractional abundances. Our results are compared with both high- and moderate-angular resolution observations of the Horsehead nebula. Our standard model is partially successful in reproducing the observations. Additional models run with different physical parameters are able to reproduce the abundance of many of the observed molecules, but we do not find a single model that fits all the observations at the same time. We discuss the suitability of a time-dependent model of a dense PDR such as ours as an estimator of the age of a PDR, provided that enough observational data exist.

The low-mass protostellar region L1527 is unusual because it contains observable abundances of unsaturated carbon-chain molecules including CnH radicals, H2Cn carbenes, cyanopolyynes, and the negative ions C4H- and C6H-, all of which are more associated with cold cores than with protostellar regions. Sakai et al. suggested that these molecules are formed in L1527 from the chemical precursor methane, which evaporates from the grains during the heat-up of the region. With the gas-phase osu.03.2008 network extended to include negative ions of the families Cn-, and CnH-, as well as the newly detected C3N-, we modeled the chemistry that occurs following methane evaporation at T~ 25-30 K. We are able to reproduce most of the observed molecular abundances in L1527 at a time of ~5000 yr. At later times, the overall abundance of anions become greater than that of electrons, which has an impact on many organic species and ions. The anion-to-neutral ratio in our calculation is in good agreement with observation for C6H- but exceeds the observed ratio by more than three orders of magnitude for C4H-. In order to explain this difference, further investigation is needed on the rate coefficients for electron attachment and other reactions regarding anions.

Sakai et al. have observed long-chain unsaturated hydrocarbons and cyanopolyynes in the low-mass star-forming region L1527, and have attributed this result to a gas-phase ion-molecule chemistry, termed ``Warm Carbon Chain Chemistry'', which occurs during and after the evaporation of methane from warming grains. The source L1527 is an envelope surrounding a Class 0/I protostar with regions that possess a slightly elevated temperature of ~30 K. The molecules detected by Sakai et al. are typically associated only with dark molecular clouds, and not with the more evolved hot corino phase. In order to determine if L1527 is chemically distinct from a dark cloud, we compute models including various degrees of heating. The results indicate that the composition of L1527 is somewhat more likely to be due to ``Warm Carbon Chain Chemistry'' than to be a remnant of a colder phase. If so, the molecular products provide a signature of a previously uncharacterized early phase of low mass star formation, which can be characterized as a ``lukewarm'' corino. We also include predictions for other molecular species that might be observed toward candidate lukewarm corino sources. Although our calculations show that unsaturated hydrocarbons and cyanopolyynes can be produced in the gas phase as the grains warm up to 30 K, they also show that such species do not disappear rapidly from the gas as the temperature reaches 200 K, implying that such species might be detected in hot corinos and hot cores.

Gas-phase processes were long thought to be the key formation mechanisms for complex organic molecules in star-forming regions. However, recent experimental and theoretical evidence has cast doubt on the efficiency of such processes. Grain-surface chemistry is frequently invoked as a solution, but until now there have been no quantitative models taking into account both the high degree of chemical complexity and the evolving physical conditions of star-forming regions. Here, we introduce a new gas-grain chemical network, wherein a wide array of complex species may be formed by reactions involving radicals. The radicals we consider (H, OH, CO, HCO, CH3, CH3O, CH2OH, NH and NH2) are produced primarily by cosmic ray-induced photodissociation of the granular ices formed during the colder, earlier stages of evolution. The gradual warm-up of the hot core is crucial to the formation of complex molecules, allowing the more strongly-bound radicals to become mobile on grain surfaces. This type of chemistry is capable of reproducing the high degree of complexity seen in Sgr B2(N), and can explain the observed abundances and temperatures of a variety of previously detected complex organic molecules, including structural isomers. Many other complex species are predicted by this model, and several of these species may be detectable in hot cores. Differences in the chemistry of high- and low-mass star-formation are also addressed; greater chemical complexity is expected where evolution timescales are longer.

Virtually all detailed gas-phase models of the chemistry of dense interstellar clouds exclude polycyclic aromatic hydrocarbons (PAH's). This omission is unfortunate because from the few studies that have been done on the subject, it is known that the inclusion of PAH's can affect the gas-phase chemistry strongly. We have added PAH's to our network to determine the role they play in the chemistry of cold dense cores. In the models presented here, we include radiative attachment to form PAH-, mutual neutralization between PAH anions and small positively-charged ions, and photodetachment. We also test the sensitivity of our results to changes in the size and abundance of the PAH's. Our results confirm that the inclusion of PAH's changes many of the calculated abundances of smaller species considerably. In TMC-1, the general agreement with observations is significantly improved contrary to L134N. This may indicate a difference in PAH properties between the two regions. With the inclusion of PAH's in dense cloud chemistry, high-metal elemental abundances give a satisfactory agreement with observations. As a result, we do not need to decrease the observed elemental abundances of all metals and we do not need to vary the elemental C/O ratio in order to produce large abundances of carbon species in TMC-1 (CP).

We investigate molecular evolution in a star-forming core that is initially a hydrostatic starless core and collapses to form a low-mass protostar. The results of a one-dimensional radiation-hydrodynamics calculation are adopted as a physical model of the core. We first derive radii at which CO and large organic species sublimate. CO sublimation in the central region starts shortly before the formation of the first hydrostatic core. When the protostar is born, the CO sublimation radius extends to 100 AU, and the region inside $\lesssim 10$ AU is hotter than 100 K, at which some large organic species evaporate. We calculate the temporal variation of physical parameters in infalling shells, in which the molecular evolution is solved using an updated gas-grain chemical model to derive the spatial distribution of molecules in a protostellar core. The shells pass through the warm region of $10 -100$ K in several $\times$ $10^4$ yr, and fall into the central star $\sim 100$ yr after they enter the region where $T \gtrsim 100$ K. We find that large organic species are formed mainly via grain-surface reactions at temperatures of $20 -40$ K and then desorbed into the gas-phase at their sublimation temperatures. Carbon-chain species can be formed by a combination of gas-phase reactions and grain-surface reactions following the sublimation of CH$_4$. Our model also predicts that CO$_2$ is more abundant in isolated cores, while gas-phase large organic species are more abundant in cores embedded in ambient clouds.

We study how uncertainties in the rate coefficients of chemical reactions in the RATE06 database affect abundances and column densities of key molecules in protoplanetary disks. We randomly varied the gas-phase reaction rates within their uncertainty limits and calculated the time-dependent abundances and column densities using a gas-grain chemical model and a flaring steady-state disk model. We find that key species can be separated into two distinct groups according to the sensitivity of their column densities to the rate uncertainties. The first group includes CO, C$^+$, H$_3^+$, H$_2$O, NH$_3$, N$_2$H$^+$, and HCNH$^+$. For these species, the column densities are not very sensitive to the rate uncertainties but the abundances in specific regions are. The second group includes CS, CO$_2$, HCO$^+$, H$_2$CO, C$_2$H, CN, HCN, HNC and other, more complex species, for which high abundances and abundance uncertainties co-exist in the same disk region, leading to larger scatters in the column densities. However, even for complex and heavy molecules, the dispersion in their column densities is not more than a factor of ~4. We perform a sensitivity analysis of the computed abundances to rate uncertainties and identify those reactions with the most problematic rate coefficients. We conclude that the rate coefficients of about a hundred of chemical reactions need to be determined more accurately in order to greatly improve the reliability of modern astrochemical models. This improvement should be an ultimate goal for future laboratory studies and theoretical investigations.

Although still poorly understood, the chemistry that occurs on the surfaces of interstellar dust particles profoundly affects the growth of molecules in the interstellar medium. An important set of surface reactions produces icy mantles of many monolayers in cold and dense regions. The monolayers are dominated by water ice, but also contain CO, CO_{2}, and occasionally methanol as well as minor constituents. In this paper, the rate of production of water-ice dominated mantles is calculated for different physical conditions of interstellar clouds and for the first time images of the morphology of interstellar ices are presented. For this purpose, the continuous-time random-walk Monte Carlo simulation technique has been used. The visual extinction, density, and gas and grain temperatures are varied. It is shown that our stochastic approach can reproduce the important observation that ice mantles only grow in the denser regions.

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.

The recent detection of the hydrocarbon anion C6H- in the interstellar medium has led us to investigate the synthesis of hydrocarbon anions in a variety of interstellar and circumstellar environments. We find that the anion/neutral abundance ratio can be quite large, on the order of at least a few percent, once the neutral has more than five carbon atoms. Detailed modeling shows that the column densities of C6H- observed in IRC+10216 and TMC-1 can be reproduced. Our calculations also predict that other hydrocarbon anions, such as C4H- and C8H-, are viable candidates for detection in IRC+10216, TMC-1 and photon-dominated regions such as the Horsehead Nebula.

Aims: The gas-phase abundance of methanol in dark quiescent cores in the interstellar medium cannot be explained by gas-phase chemistry. In fact, the only possible synthesis of this species appears to be production on the surfaces of dust grains followed by desorption into the gas. Yet, evaporation is inefficient for heavy molecules such as methanol at the typical temperature of 10 K. It is necessary then to consider non-thermal mechanisms for desorption. But, if such mechanisms are considered for the production of methanol, they must be considered for all surface species. Methods: Our gas-grain network of reactions has been altered by the inclusion of a non-thermal desorption mechanism in which the exothermicity of surface addition reactions is utilized to break the bond between the product species and the surface. Our estimated rate for this process derives from a simple version of classical unimolecular rate theory with a variable parameter only loosely contrained by theoretical work. Results: Our results show that the chemistry of dark clouds is altered slightly at times up to 10^6 yr, mainly by the enhancement in the gas-phase abundances of hydrogen-rich species such as methanol that are formed on grain surfaces. At later times, however, there is a rather strong change. Instead of the continuing accretion of most gas-phase species onto dust particles, a steady-state is reached for both gas-phase and grain-surface species, with significant abundances for the former. Nevertheless, most of the carbon is contained in an undetermined assortment of heavy surface hydrocarbons. Conclusions: The desorption mechanism discussed here will be better constrained by observational data on pre-stellar cores, where a significant accretion of species such as CO has already occurred.

Aim: To determine whether or not gas-phase chemical models with homogeneous and time-independent physical conditions explain the many observed molecular abundances in astrophysical sources, it is crucial to estimate the uncertainties in the calculated abundances and compare them with the observed abundances and their uncertainties. Non linear amplification of the error and bifurcation may limit the applicability of chemical models. Here we study such effects on dense cloud chemistry. Method: Using a previously studied approach to uncertainties based on the representation of rate coefficient errors as log normal distributions, we attempted to apply our approach using as input a variety of different elemental abundances from those studied previously. In this approach, all rate coefficients are varied randomly within their log normal (Gaussian) distribution, and the time-dependent chemistry calculated anew many times so as to obtain good statistics for the uncertainties in the calculated abundances. Results: Starting with so-called ``high-metal'' elemental abundances, we found bimodal rather than Gaussian like distributions for the abundances of many species and traced these strange distributions to an extreme sensitivity of the system to changes in the ratio of the cosmic ray ionization rate zeta\_He for He and that for molecular hydrogen zeta\_H2. The sensitivity can be so extreme as to cause a region of bistability, which was subsequently found to be more extensive for another choice of elemental abundances. To the best of our knowledge, the bistable solutions found in this way are the same as found previously by other authors, but it is best to think of the ratio zeta\_He/zeta\_H2 as a control parameter perpendicular to the ''standard'' control parameter zeta/n\_H.

Aims: The production of saturated organic molecules in hot cores and corinos is not well understood. The standard approach is to assume that, as temperatures heat up during star formation, methanol and other species evaporate from grain surfaces and undergo a warm gas-phase chemistry at 100 K or greater to produce species such as methyl formate, dimethyl ether, and others. But a series of laboratory results shows that protonated ions, typical precursors to final products in ion-molecule schemes, tend to fragment upon dissociative recombination with electrons rather than just ejecting a hydrogen atom. Moreover, the specific proposed reaction to produce protonated methyl formate is now known not to occur at all. Methods: We utilize a gas-grain chemical network to probe the chemistry of the relatively ignored stage of hot core evolution during which the protostar switches on and the temperature of the surrounding gas and dust rises from 10 K to over 100 K. During this stage, surface chemistry involving heavy radicals becomes more important as surface hydrogen atoms tend to evaporate rather than react. Results: Our results show that complex species such as methyl formate, formic acid, and dimethyl ether can be produced in large abundance during the protostellar switch-on phase, but that both grain-surface and gas-phase processes help to produce most species. The longer the timescale for protostellar switch-on, the more important the surface processes.

Low-mass protostars form from condensations inside molecular clouds when gravity overwhelms thermal and magnetic supporting forces. The first phases of the formation of a solar-type star are characterized by dramatic changes not only in the physical structure but also in the chemical composition. Since PPIV (e.g., Langer et al.), exciting new developments have occurred in our understanding of the processes driving this chemical evolution. These developments include two new discoveries : 1) extremely enhanced molecular deuteration, which is caused by the freeze-out of heavy-element-bearing molecules onto grain mantles during the Prestellar Core and Class 0 source phases; and 2) hot corinos, which are warm and dense regions at the center of Class 0 source envelopes and which are characterized by a multitude of complex organic molecules. In this chapter we will review these two new topics, and will show how they contribute to our understanding of the first phases of solar-type stars.

The gas-phase chemistry of dark clouds has been studied with a treatment of uncertainties caused both by errors in individual rate coefficients and uncertainties in physical conditions. Moreover, a sensitivity analysis has been employed to attempt to determine which reactions are most important in the chemistry of individual species. The degree of overlap between calculated errors in abundances and estimated observational errors has been used as an initial criterion for the goodness of the model and the determination of a best 'chemical' age of the source. For the well-studied sources L134N and TMC-1CP, best agreement is achieved at so-called "early times" ~10$^{5}$ yr, in agreement with previous calculations but here put on a firmer statistical foundation. A more detailed criterion for agreement, which takes into account the degree of disagreement, is also proposed. Poorly understood but critical classes of reactions are delineated, especially reactions between ions and polar neutrals. Such reactions will have to be understood better before the chemistry can be made more secure. Nevertheless, the level of agreement is low enough to indicate that a static picture of physical conditions without consideration of interactions with grain surfaces is inappropriate for a complete understanding of the chemistry.

Continuous-time, random-walk Monte Carlo simulations of H2 formation on grains have been performed for surfaces that are stochastically heated by photons. We have assumed diffuse cloud conditions and used a variety of grains of varying roughness and size based on olivine. The simulations were performed at different optical depths. We confirmed that small grains (r <= 0.02 micron) have low modal temperatures with strong fluctuations, which have a large effect on the efficiency of the formation of molecular hydrogen. The grain size distribution highly favours small grains and therefore H2 formation on these particles makes a large contribution to the overall formation rate for all but the roughest surfaces. We find that at A_V=0 only the roughest surfaces can produce the required amount of molecular hydrogen, but by A_V=1, smoother surfaces are possible alternatives. Use of a larger value for the evaporation energy of atomic hydrogen, but one still consistent with experiment, allows smoother surfaces to produce more H2.

We report a detailed attempt to model the ortho-to-para abundance ratio of c-C3H2 so as to reproduce observed values in the cores of the well-known source TMC-1. According to observations, the ortho-to-para ratios vary, within large uncertainties, from a low of near unity to a high of approximately three depending on the core. We used the osu.2003 network of gas-phase chemical reactions augmented by reactions that specifically consider the formation, depletion, and interconversion of the ortho and para forms of the c-C3H2 and its precursor ion c-C3H3+. We investigated the sensitivity of the calculated ortho-to-para ratio for c-C3H2 to a large number of factors. For the less evolved cores C, CP, and D, we had no difficulty reproducing the observed ortho-to-para ratios of 1-2.

It is not common to consider the role of uncertainties in the rate coefficients used in interstellar gas-phase chemical models. In this paper, we report a new method to determine both the uncertainties in calculated molecular abundances and their sensitivities to underlying uncertainties in the kinetic data utilized. The method is used in hot core models to determine if previous analyses of the age and the applicable cosmic-ray ionization rate are valid. We conclude that for young hot cores ($\le 10^4$ yr), the modeling uncertainties related to rate coefficients are reasonable so that comparisons with observations make sense. On the contrary, the modeling of older hot cores is characterized by strong uncertainties for some of the important species. In both cases, it is crucial to take into account these uncertainties to draw conclusions from the comparison of observations with chemical models.

The gravitational collapse of a spherical cloud core is investigated by numerical calculations. The initial conditions of the core lie close to the critical Bonnor-Ebert sphere with a central density of \sim 10^4 cm^{-3} in one model (alpha=1.1), while gravity overwhelms pressure in the other (alpha=4.0), where alpha is the internal gravity-to-pressure ratio. The alpha=1.1 model shows reasonable agreement with the observed velocity field in prestellar cores. Molecular distributions in cores are calculated by solving a chemical reaction network that includes both gas-phase and grain-surface reactions. When the central density of the core reaches 10^5 cm^{-3}, carbon-bearing species are significantly depleted in the central region of the alpha=1.1 model, while the depletion is only marginal in the other model. The two different approaches encompass the observed variations of molecular distributions in different prestellar cores, suggesting that molecular distributions can be probes of contraction or accumulation time scales of cores. The central enhancement of the NH3/N2H+ ratio, which is observed in some prestellar cores, can be reproduced under certain conditions by adopting recently measured branching fractions for N2H+ recombination. Various molecular species, such as CH3OH and CO2, are produced by grain-surface reactions. The ice composition depends sensitively on the assumed temperature. Multi-deuterated species are included in our most recent gas-grain chemical network. The deuterated isotopomers of H3+ are useful as probes of the central regions of evolved cores, in which gas-phase species with heavy elements are strongly depleted. At 10 K, our model can reproduce the observed abundance ratio of ND3/NH3, but underestimates the isotopic ratios of deuterated to normal methanol.

We report a theoretical study of sulphur chemistry, as applied to hot cores, where S-bearing molecular ratios have been previously proposed and used as chemical clocks. As in previous models, we follow the S-bearing molecular composition after the injection of grain mantle components into the gas phase. For this study, we developed a time-dependent chemical model with up-to-date reaction rate coefficients. We ran several cases, using different realistic chemical compositions for the grain mantles and for the gas prior to mantle evaporation. The modeling shows that S-bearing molecular ratios depend very critically on the gas temperature and density, the abundance of atomic oxygen, and, most importantly, on the form of sulphur injected in the gas phase, which is very poorly known. Consequently, ratios of S-bearing molecules cannot be easily used as chemical clocks. However, detailed observations and careful modeling of both physical and chemical structure can give hints on the source age and constrain the mantle composition (i.e. the form of sulphur in cold molecular clouds) and, thus, help to solve the mystery of the sulphur depletion. We analyse in detail the cases of Orion and IRAS16293-2422. The comparison of the available observations with our model suggests that the majority of sulphur released from the mantles is mainly in, or soon converted into, atomic form.

Motivated by recent observations of deuterated molecules in the interstellar medium, we examine the production of HD and D$_2$ molecules on dust grain surfaces. A mechanism for the enhanced production of these deuterated molecules is studied. This mechanism applies on grain surfaces on which D atoms stick more strongly than H atoms, under conditions of low flux and within a suitable range of temperatures. It is shown that under these conditions the production rates of HD and D$_2$ are greatly enhanced (vs. the H$_2$ production rate) compared with the expected rates based on the adsorption of gas-phase atomic abundances of D and H. The enhancement in the formation rate of HD is comparable with the enhancement due to gas-phase ion-molecule reactions in diffuse clouds.

In this paper, we compare the results of the modeling of a protoplanetary disk chemical evolution obtained with the UMIST95 and ``New Standard Model'' (NSM) chemical databases. Assuming the same initial conditions, it is found that the substitution of one chemical network by another causes almost no difference for the disk ionization degree. In contrast, the NSM and UMIST95 abundances of CO can differ by a factor of a hundred at some regions of the disk surface. However, relevant CO vertical column densities differ much less, at most by a factor of a few. In addition, we synthesize the single-dish CO(J=3-2) line by means of the 2D line radiative transfer for both considered chemical networks. It is shown that the intensity of this line in the case of the UMIST95 abundances is lower compared to the NSM case by about 15%.

We report the first detection of triply-deuterated methanol, with 12 observed transitions, towards the low-mass protostar IRAS 16293-2422, as well as multifrequency observations of 13CH3OH, used to derive the column density of the main isotopomer CH3OH. The derived fractionation ratio [CD3OH]/[CH3OH] averaged on a 10'' beam is 1.4%. Together with previous CH2DOH and CHD2OH observations, the present CD3OH observations are consistent with a formation of methanol on grain surfaces, if the atomic D/H ratio is 0.1 to 0.3 in the accreting gas. Such a high atomic ratio can be reached in the frame of gas-phase chemical models including all deuterated isotopomers of H3+.

We have studied deuterium fractionation on interstellar grains with the use of an exact method known as the direct master equation approach. We consider conditions pertinent to dense clouds at late times when the hydrogen is mostly in molecular form and a large portion of the gas-phase carbon has already been converted to carbon monoxide. Hydrogen, oxygen and deuterium atoms, as well as CO molecules, are allowed to accrete on to dust particles and react there to produce various stable molecules. The surface abundances, as well as the abundance ratios between deuterated and normal isotopomers, are compared with those calculated with the Monte Carlo approach. We find that the agreement between the Monte Carlo and the direct master equation methods can be made as close as desired. Compared with previous examples of the use of the direct master equation approach, our present method is much more efficient. It should now be possible to run large-scale gas-grain models in which the diffusive dust chemistry is handled `exactly'.

We have used the master equation approach to study a moderately complex network of diffusive reactions occurring on the surfaces of interstellar dust particles. This network is meant to apply to dense clouds in which a large portion of the gas-phase carbon has already been converted to carbon monoxide. Hydrogen atoms, oxygen atoms, and CO molecules are allowed to accrete onto dust particles and their chemistry is followed. The stable molecules produced are oxygen, hydrogen, water, carbon dioxide (CO2), formaldehyde (H2CO), and methanol (CH3OH). The surface abundances calculated via the master equation approach are in good agreement with those obtained via a Monte Carlo method but can differ considerably from those obtained with standard rate equations.

The molecular evolution that occurs in collapsing prestellar cores is investigated. To model the dynamics, we adopt the Larson-Penston (L-P) solution and analogues with slower rates of collapse. For the chemistry, we utilize the new standard model (NSM) with the addition of deuterium fractionation and grain-surface reactions treated via the modified rate approach. The use of surface reactions distinguishes the present work from our previous model. We find that these reactions efficiently produce H2O, H2CO, CH3OH, N2, and NH3 ices. In addition, the surface chemistry influences the gas-phase abundances in a variety of ways. The current reaction network along with the L-P solution allows us to reproduce satisfactorily most of the molecular column densities and their radial distributions observed in L1544. The agreement tends to worsen with models that include strongly delayed collapse rates. Inferred radial distributions in terms of fractional abundances are somewhat harder to reproduce. In addition to our standard chemical model, we have also run a model with the UMIST gas-phase chemical network. The abundances of gas-phase S-bearing molecules such as CS and CCS are significantly affected by uncertainties in the gas-phase chemical network. In all of our models, the column density of N2H+ monotonically increases as the central density of the core increases during collapse from 3 10^4 cm-3 to 3 10^7 cm-3. Thus, the abundance of this ion can be a probe of evolutionary stage. Molecular D/H ratios in assorted cores are best reproduced in the L-P picture with the conventional rate coefficients for fractionation reactions. If we adopt the newly measured and calculated rate coefficients, the D/H ratios, especially N2D+/N2H+, become significantly lower than the observed values.

We report the first detection of doubly-deuterated methanol (CHD2OH), as well as firm detections of the two singly-deuterated isotopomers of methanol (CH2DOH and CH3OD), towards the solar-type protostar IRAS16293-2422. From the present multifrequency observations, we derive the following abundance ratios: [CHD2OH]/[CH3OH] = 0.2 +/- 0.1, [CH2DOH]/[CH3OH] = 0.9 +/- 0.3, [CH3OD]/[CH3OH] = 0.04 +/- 0.02. The total abundance of the deuterated forms of methanol is greater than that of its normal hydrogenated counterpart in the circumstellar material of IRAS16293-2422, a circumstance not previously encountered. Formaldehyde, which is thought to be the chemical precursor of methanol, possesses a much lower fraction of deuterated isotopomers (~ 20%) with respect to the main isotopic form in IRAS16293-2422. The observed fractionation of methanol and formaldehyde provides a severe challenge to both gas-phase and grain-surface models of deuteration. Two examples of the latter model are roughly in agreement with our observations of CHD2OH and CH2DOH if the accreting gas has a large (0.2-0.3) atomic D/H ratio. However, no gas-phase model predicts such a high atomic D/H ratio, and hence some key ingredient seems to be missing.

The formation of singly and doubly deuterated isotopomers of formaldehyde and of singly, doubly, and multiply deuterated isotopomers of methanol on interstellar grain surfaces has been studied with a semi-empirical modified rate approach and a Monte Carlo method in the temperature range 10-20 K. Agreement between the results of the two methods is satisfactory for all major and many minor species throughout this range. If gas-phase fractionation can produce a high abundance of atomic deuterium, which then accretes onto grain surfaces, diffusive surface chemistry can produce large abundances of deuterated species, especially at low temperatures and high gas densities. Warming temperatures will then permit these surface species to evaporate into the gas, where they will remain abundant for a considerable period. We calculate that the doubly deuterated molecules CHD2OH and CH2DOD are particularly abundant and should be searched for in the gas phase of protostellar sources. For example, at 10 K and high density, these species can achieve up to 10-20% of the abundance of methanol.

We have investigated the two-dimensional (R,Z) distribution of deuterated molecular species in circumstellar disks around young stellar objects. The abundance ratios between singly deuterated and normal molecules (``D/H ratios'') in disks evolve in a similar way as in molecular clouds. Fractionation is caused by rapid exchange reactions that are exothermic because of energy differences between deuterated and normal species. In the midplane region, where molecules are heavily depleted onto grain surfaces, the D/H ratios of gaseous molecules are higher than at larger heights. The D/H ratios for the vertical column densities of NH3, H2O, and HCO+ are sensitive to the temperature, and decrease significantly with decreasing radial distance for R < 300 AU. The analogous D/H ratios for CH4 and H2CO, on the other hand, are not very sensitive to the temperature in the range (T=10-50 K) we are concerned with, and do not decrease with decreasing R at R > 50 AU. The D/H column-density ratios also depend on disk mass. In a disk with a larger mass, the ratios of deuterated species to normal species are higher, because of heavier depletion of molecules onto grains. In the second part of the paper, we report molecular column densities for disks embedded in ambient cloud gas. Our results suggest that CN and HCO+ can be tracers of gaseous disks, especially if the central object is a strong X-ray source. Our results also suggest that the radial distributions of CN, C2H, HCN, and H2CO may vary among disks depending on the X-ray luminosity of the central star.

We have investigated molecular distributions in protoplanetary disks, adopting a disk model with a temperature gradient in the vertical direction. The model produces sufficiently high abundances of gaseous CO and HCO+ to account for line observations of T Tauri stars using a sticking probability of unity and without assuming any non-thermal desorption. In regions of radius R > 10 AU, with which we are concerned, the temperature increases with increasing height from the midplane. In a warm intermediate layer, there are significant amounts of gaseous molecules owing to thermal desorption and efficient shielding of ultraviolet radiation by the flared disk. The column densities of HCN, CN, CS, H2CO, HNC and HCO+ obtained from our model are in good agreement with the observations of DM Tau, but are smaller than those of LkCa15. Molecular line profiles from our disk models are calculated using a 2-dimensional non-local-thermal-equilibrium (NLTE) molecular-line radiative transfer code for a direct comparison with observations. Deuterated species are included in our chemical model. The molecular D/H ratios in the model are in reasonable agreement with those observed in protoplanetary disks.

We have investigated the evolution and distribution of molecules in collapsing prestellar cores via numerical chemical models, adopting the Larson-Penston solution and its delayed analogues to study collapse. Molecular abundances and distributions in a collapsing core are determined by the balance among the dynamical, chemical and adsorption time scales. When the central density n_H of a prestellar core with the Larson-Penston flow rises to 3 10^6 cm^{-3}, the CCS and CO column densities are calculated to show central holes of radius 7000 AU and 4000 AU, respectively, while the column density of N2H+ is centrally peaked. These predictions are consistent with observations of L1544. If the dynamical time scale of the core is larger than that of the Larson-Penston solution owing to magnetic fields, rotation, or turbulence, the column densities of CO and CCS are smaller, and their holes are larger than in the Larson-Penston core with the same central gas density. On the other hand, N2H+ and NH3 are more abundant in the more slowly collapsing core. Therefore, molecular distributions can probe the collapse time scale of prestellar cores. Deuterium fractionation has also been studied via numerical calculations. The deuterium fraction in molecules increases as a core evolves and molecular depletion onto grains proceeds. When the central density of the core is n_H=3 10^6 cm^{-3}, the ratio DCO+/HCO+ at the center is in the range 0.06-0.27, depending on the collapse time scale and adsorption energy; this range is in reasonable agreement with the observed value in L1544.