We present deep CCS and HC$_7$N observations of the L1495-B218 filaments in the Taurus molecular cloud obtained using the K-band focal plane array on the 100m Green Bank Telescope. We observed the L1495-B218 filaments in CCS $J_N$ = 2$_1$$-$1$_0$ and HC$_7$N $J$ = 21$-$20 with a spectral resolution of 0.038 km s$^{-1}$ and an angular resolution of 31$''$. We observed strong CCS emission in both evolved and young regions and weak emission in two evolved regions. HC$_7$N emission is observed only in L1495A-N and L1521D. We find that CCS and HC$_7$N intensity peaks do not coincide with NH$_3$ or dust continuum intensity peaks. We also find that the fractional abundance of CCS does not show a clear correlation with the dynamical evolutionary stage of dense cores. Our findings and chemical modeling indicate that the fractional abundances of CCS and HC$_7$N are sensitive to the initial gas-phase C/O ratio, and they are good tracers of young condensed gas only when the initial C/O is close to solar value. Kinematic analysis using multiple lines including NH$_3$, HC$_7$N, CCS, CO, HCN, \& HCO$^+$ suggests that there may be three different star formation modes in the L1495-B218 filaments. At the hub of the filaments, L1495A/B7N has formed a stellar cluster with large-scale inward flows (fast mode), while L1521D, a core embedded in a filament, is slowly contracting due to its self-gravity (slow mode). There is also one isolated core that appears to be marginally stable and may undergo quasi-static evolution (isolated mode).

Based on the interstellar CO/H2 ratio, carbon monoxide-based censuses of protoplanetary disks in Lupus, sigma Orionis, and Chamaeleon I found no disks more massive than the minimum-mass solar nebula, which is inconsistent with the existence of exoplanets more massive than Jupiter. Observations and models are converging on the idea that ionization-driven chemistry depletes carbon monoxide in T-Tauri disks. Yet the extent of CO depletion depends on the incident flux of ionizing radiation, and some T-Tauri stars may have winds strong enough to shield their disks from cosmic rays. There is also a range of X-ray luminosities possible for a given stellar mass. Here we use a suite of chemical models, each with a different incident X-ray or cosmic-ray flux, to assess whether CO depletion is a typical outcome for T-Tauri disks. We find that CO dissociation in the outer disk is a robust result for realistic ionization rates, with abundance reductions between 70% and 99.99% over 2 Myr of evolution. Furthermore, after the initial dissociation epoch, the inner disk shows some recovery of the CO abundance from CO2 dissociation. In highly ionized disks, CO recovery in the inner disk combined with depletion in the outer disk creates a centrally peaked CO abundance distribution. The emitting area in rare CO isotopologues may be an indirect ionization indicator: in a cluster of disks with similar ages, those with the most compact CO isotopologue emission see the highest ionization rates.

We investigate the abiotic production of oxygen and its photochemical byproduct ozone through water vapor photolysis in moist atmospheres of temperate terrestrial exoplanets. The amount of water vapor available for photolysis in the middle atmosphere of a planet can be limited by an atmospheric cold-trap, the formation of which largely depends on the amount of non-condensable gases. We study this effect using a photochemical model coupled to a 1D radiative-convective equilibrium model in atmospheres with N2, CO2 and H2O as the main constituents. We find that in atmospheres with a low N2 inventory, water vapor mixing ratios in the middle atmosphere can be over two orders of magnitude higher compared to atmospheres with an Earth-like N2 inventory. Without a strong surface sink, the non-condensable oxygen can build up rapidly, drying out the upper atmosphere. With a moderate surface sink, the planet can approach a steady state with significant oxygen mixing ratios in which oxygen production is balanced by surface uptake. We use a radiative transfer model to study the spectroscopic fingerprint of these atmospheres in transit observations. Spectral signatures of abiotic oxygen and ozone can be of comparable magnitude as in spectra of Earth seen as an exoplanet. Middle atmospheric water vapor is unlikely to be a usable indicator of the abiotic origin of oxygen because of the influence of oxygen on the water vapor distribution. This suggests that atmospheric oxygen and ozone cannot be used as binary bioindicators and their interpretation will likely require atmospheric and planetary models.

The protonated form of CO2, HOCO+, is assumed to be an indirect tracer of CO2 in the millimeter/submillimeter regime since CO2 lacks a permanent dipole moment. Here, we report the detection of two rotational emission lines (4 0,4-3 0,3) and (5 0,5-4 0,4) of HOCO+ in IRAS 16293-2422. For our observations, we have used EMIR heterodyne 3 mm receiver of the IRAM 30m telescope. The observed abundance of HOCO+ is compared with the simulations using the 3-phase NAUTILUS chemical model. Implications of the measured abundances of HOCO+ to study the chemistry of CO2 ices using JWST-MIRI and NIRSpec are discussed as well.

Turbulence is the leading candidate for angular momentum transport in protoplanetary disks and therefore influences disk lifetimes and planet formation timescales. However, the turbulent properties of protoplanetary disks are poorly constrained observationally. Recent studies have found turbulent speeds smaller than what fully-developed MRI would produce (Flaherty et al. 2015, 2017). However, existing studies assumed a constant CO/H2 ratio of 0.0001 in locations where CO is not frozen-out or photo-dissociated. Our previous studies of evolving disk chemistry indicate that CO is depleted by incorporation into complex organic molecules well inside the freeze-out radius of CO. We consider the effects of this chemical depletion on measurements of turbulence. Simon et al. (2015) suggested that the ratio of the peak line flux to the flux at line center of the CO J=3-2 transition is a reasonable diagnostic of turbulence, so we focus on that metric, while adding some analysis of the more complex effects on spatial distribution. We simulate the emission lines of CO based on chemical evolution models presented in Yu et al. (2016), and find that the peak-to-trough ratio changes as a function of time as CO is destroyed. Specifically, a CO-depleted disk with high turbulent velocity mimics the peak-to-trough ratios of a non-CO-depleted disk with lower turbulent velocity. We suggest that disk observers and modelers take into account the possibility of CO depletion when using line peak-to-trough ratios to constrain the degree of turbulence in disks. Assuming that CO/H2 = 0.0001 at all disk radii can lead to underestimates of turbulent speeds in the disk by at least 0.2 km/s.

Using forward models for representative exoplanet atmospheres and a radiometric instrument model, we have generated synthetic observational data to explore how well the major C- and O-bearing chemical species (CO, CO2, CH4, and H2O), important for determining atmospheric opacity and radiation balance, can be constrained by transit measurements as a function of spectral wavelength coverage. This work features simulations for a notional transit spectroscopy mission and compares two cases for instrument spectral coverage (wavelength coverage from 0.5-2.5 ${\mu}$m and 0.5-5 ${\mu}$m. The simulation is conducted on a grid with a range of stellar magnitudes and incorporates a full retrieval of atmospheric model parameters. We consider a range of planets from sub-Neptunes to hot Jupiters and include both low and high mean molecular weight atmospheres. We find that including the 2.5-5 ${\mu}$m wavelength range provides a significant improvement in the degree of constraint on the retrieved molecular abundances: up to ~3 orders of magnitude for a low mean molecular weight atmosphere (${\mu}$=2.3) and up to a factor of ~6 for a high mean molecular weight atmosphere (${\mu}$=28). These decreased uncertainties imply that broad spectral coverage between the visible and the mid-infrared is an important tool for understanding the chemistry and composition of exoplanet atmospheres. This analysis suggests that the JWST/NIRSpec 0.6-5 ${\mu}$m prism spectroscopy mode, or similar wavelength coverage with possible future missions, will be an important resource for exoplanet atmospheric characterization.

Gas in protostellar disks provides the raw material for giant planet formation and controls the dynamics of the planetesimal-building dust grains. Accurate gas mass measurements help map the observed properties of planet-forming disks onto the formation environments of known exoplanets. Rare isotopologues of carbon monoxide (CO) have been used as gas mass tracers for disks in the Lupus star-forming region, with an assumed interstellar CO/H$_2$ abundance ratio. Unfortunately, observations of T-Tauri disks show that CO abundance is not interstellar---a finding reproduced by models that show CO abundance decreasing both with distance from the star and as a function of time. Here we present radiative transfer simulations that assess the accuracy of CO-based disk mass measurements. We find that the combination of CO chemical depletion in the outer disk and optically thick emission from the inner disk leads observers to underestimate gas mass by more than an order of magnitude if they use the standard assumptions of interstellar CO/H$_2$ ratio and optically thin emission. Furthermore, CO abundance changes on million-year timescales, introducing an age/mass degeneracy into observations. To reach factor of a few accuracy for CO-based disk mass measurements, we suggest that observers and modelers adopt the following strategies: (1) select the low-$J$ transitions; (2) observe multiple CO isotopologues and use either intensity ratios or normalized line profiles to diagnose CO chemical depletion; and (3) use spatially resolved observations to measure the CO abundance distribution.

We present results of an investigation into the formation of nitrogen-bearing molecules in the atmosphere of Titan. We extend a previous model (Li et al. 2015, 2016) to cover the region below the tropopause, so the new model treats the atmosphere from Titan's surface to an altitude of 1500 km. We consider the effects of condensation and sublimation using a continuous, numerically stable method. This is coupled with parameterized treatments of the sedimentation of the aerosols and their condensates, and the formation of haze particles. These processes affect the abundances of heavier species such as the nitrogen-bearing molecules, but have less effect on the abundances of lighter molecules. Removal of molecules to form aerosols also plays a role in determining the mixing ratios, in particular of HNC, HC3N and HCN. We find good agreement with the recently detected mixing ratios of C2H5CN, with condensation playing an important role in determining the abundance of this molecule below 500 km. Of particular interest is the chemistry of acrylonitrile (C2H3CN) which has been suggested by Stevenson et al. (2015) as a molecule that could form biological membranes in an oxygen-deficient environment. With the inclusion of haze formation we find good agreement of our model predictions of acrylonitrile with the available observations.

The gas near the midplanes of planet-forming protostellar disks remains largely unprobed by observations due to the high optical depth of commonly observed molecules such as CO and H$_2$O. However, rotational emission lines from rare molecules may have optical depths near unity in the vertical direction, so that the lines are strong enough to be detected, yet remain transparent enough to trace the disk midplane. Here we present a chemical model of an evolving T-Tauri disk and predict the optical depths of rotational transitions of $^{12}$C$^{16}$O, $^{13}$C$^{16}$O, $^{12}$C$^{17}$O and $^{12}$C$^{18}$O. The MRI-active disk is primarily heated by the central star due to the formation of the dead zone. CO does not freeze out in our modeled region within $70$AU around a sunlike star. However, the abundance of CO decreases because of the formation of complex organic molecules (COM), producing an effect that can be misinterpreted as the "snow line". These results are robust to variations in our assumptions about the evolution of the gas to dust ratio. The optical depths of low-order rotational lines of C$^{17}$O are around unity, making it possible to see into the disk midplane using C$^{17}$O. Combining observations with modeled C$^{17}$O$/$H$_2$ ratios, like those we provide, can yield estimates of protoplanetary disks' gas masses.

We present a new method for assessing the intrinsic 3D shape of prestellar cores from molecular column densities. We have employed hydrodynamic simulations of contracting, isothermal cores considering three intrinsic geometries: spherical, cylindrical/filamentary and disk-like. We have coupled our hydrodynamic simulations with non-equilibrium chemistry. We find that a) when cores are observed very elongated (i.e. for aspect ratios $\le$ 0.15) the intrinsic 3D geometry can be probed by their 2D molecular emission maps, since these exhibit significant qualitative morphological differences between cylindrical and disk-like cores. Specifically, if a disk-like core is observed as a filamentary object in dust emission, then it will be observed as two parallel filaments in $\rm{N_2H^{+}}$; b) for cores with higher aspect ratios (i.e. 0.15 $\sim$ 0.9) we define a metric $\Delta$ that quantifies whether a molecular column density profile is centrally peaked, depressed or flat. We have identified one molecule ($\rm{CN}$) for which $\Delta$ as a function of the aspect ratio probes the 3D geometry of the core; and c) for cores with almost circular projections (i.e. for aspect ratios $\sim$ 1), we have identified three molecules ($\rm{OH}$, $\rm{CO}$ and $\rm{H_2CO}$) that can be used to probe the intrinsic 3D shape by close inspection of their molecular column density radial profiles. We alter the temperature and the cosmic-ray ionization rate and demonstrate that our method is robust against the choice of parameters.

Conditions in the protosolar nebula have left their mark in the composition of cometary volatiles, thought to be some of the most pristine material in the solar system. Cometary compositions represent the end point of processing that began in the parent molecular cloud core and continued through the collapse of that core to form the protosun and the solar nebula, and finally during the evolution of the solar nebula itself as the cometary bodies were accreting. Disentangling the effects of the various epochs on the final composition of a comet is complicated. But comets are not the only source of information about the solar nebula. Protostellar disks around young stars similar to the protosun provide a way of investigating the evolution of disks similar to the solar nebula while they are in the process of evolving to form their own solar systems. In this way we can learn about the physical and chemical conditions under which comets formed, and about the types of dynamical processing that shaped the solar system we see today. This paper summarizes some recent contributions to our understanding of both cometary volatiles and the composition, structure and evolution of protostellar disks.

The ratio of mass and magnetic flux determines the relative importance of magnetic and gravitational forces in the evolution of molecular clouds and their cores. Its measurement is thus central in discriminating between different theories of core formation and evolution. Here we discuss the effect of chemical depletion on measurements of the mass-to-flux ratio using the same molecule (OH) both for Zeeman measurements of the magnetic field and the determination of the mass of the region. The uncertainties entering through the OH abundance in determining separately the magnetic field and the mass of a region have been recognized in the literature. It has been proposed however that, when comparing two regions of the same cloud, the abundance will in both cases be the same. We show that this assumption is invalid. We demonstrate that when comparing regions with different densities, the effect of OH depletion in measuring changes of the mass-to-flux ratio between different parts of the same cloud can even reverse the direction of the underlying trends (for example, the mass-to-flux ratio may appear to decrease as we move to higher density regions). The systematic errors enter primarily through the inadequate estimation of the mass of the region.

ABRIDGED: Context: HI and CO large scale surveys of the Milky Way trace the diffuse atomic clouds and the dense shielded regions of molecular hydrogen clouds. However, until recently, we have not had spectrally resolved C+ surveys to characterize the photon dominated interstellar medium, including, the H2 gas without C, the CO-dark H2, in a large sample of clouds. Aims: To use a sparse Galactic plane survey of the 1.9 THz [C II] spectral line from the Herschel Open Time Key Programme, Galactic Observations of Terahertz C+ (GOT C+), to characterize the H2 gas without CO in a statistically significant sample of clouds. Methods: We identify individual clouds in the inner Galaxy by fitting [CII] and CO isotopologue spectra along each line of sight. We combine these with HI spectra, along with excitation models and cloud models of C+, to determine the column densities and fractional mass of CO-dark H2 clouds. Results: We identify 1804 narrow velocity [CII] interstellar cloud components in different categories. About 840 are diffuse molecular clouds with no CO, 510 are transition clouds containing [CII] and 12CO, but no 13CO, and the remainder are dense molecular clouds containing 13CO emission. The CO-dark H2 clouds are concentrated between Galactic radii 3.5 to 7.5 kpc and the column density of the CO-dark H2 layer varies significantly from cloud-to-cloud with an average 9X10^(20) cm-2. These clouds contain a significant fraction of CO-dark H2 mass, varying from ~75% for diffuse molecular clouds to ~20% for dense molecular clouds. Conclusions: We find a significant fraction of the warm molecular ISM gas is invisible in HI and CO, but is detected in [CII]. The fraction of CO-dark H2 is greatest in the diffuse clouds and decreases with increasing total column density, and is lowest in the massive clouds.

Comparison of linewidths of spectral line profiles of ions and neutral molecules have been recently used to estimate the strength of the magnetic field in turbulent star-forming regions. However, the ion (HCO+) and neutral (HCN) species used in such studies may not be necessarily co-evolving at every scale and density and may thus not trace the same regions. Here, we use coupled chemical/dynamical models of evolving prestellar molecular cloud cores including non-equilibrium chemistry, with and without magnetic fields, to study the spatial distribution of HCO+ and HCN, which have been used in observations of spectral linewidth differences to date. In addition, we seek new ion-neutral pairs that are good candidates for such observations because they have similar evolution and are approximately co-spatial in our models. We identify three such good candidate pairs: HCO+/NO, HCO+/CO, and NO+/NO.

We study the effect that non-equilibrium chemistry in dynamical models of collapsing molecular cloud cores has on measurements of the magnetic field in these cores, the degree of ionization, and the mean molecular weight of ions. We find that OH and CN, usually used in Zeeman observations of the line-of-sight magnetic field, have an abundance that decreases toward the center of the core much faster than the density increases. As a result, Zeeman observations tend to sample the outer layers of the core and consistently underestimate the core magnetic field. The degree of ionization follows a complicated dependence on the number density at central densities up to 10^5 cm^{-3} for magnetic models and 10^6 cm^{-3} in non-magnetic models. At higher central densities the scaling approaches a power-law with a slope of -0.6 and a normalization which depends on the cosmic-ray ionization rate {\zeta} and the temperature T as ({\zeta}T)^1/2. The mean molecular weight of ions is systematically lower than the usually assumed value of 20 - 30, and, at high densities, approaches a value of 3 due to the asymptotic dominance of the H3+ ion. This significantly lower value implies that ambipolar diffusion operates faster.

We combine dynamical and non-equilibrium chemical modeling of evolving prestellar molecular cloud cores, and explore the evolution of molecular abundances in the contracting core. We model both magnetic cores, with varying degrees of initial magnetic support, and non-magnetic cores, with varying collapse delay times. We explore, through a parameter study, the competing effects of various model parameters in the evolving molecular abundances, including the elemental C/O ratio, the temperature, and the cosmic-ray ionization rate. We find that different models show their largest quantitative differences at the center of the core, whereas the outer layers, which evolve slower, have abundances which are severely degenerate among different dynamical models. There is a large range of possible abundance values for different models at a fixed evolutionary stage (central density), which demonstrates the large potential of chemical differentiation in prestellar cores. However, degeneracies among different models, compounded with uncertainties induced by other model parameters, make it difficult to discriminate among dynamical models. To address these difficulties, we identify abundance ratios between particular molecules, the measurement of which would have maximal potential for discrimination among the different models examined here. In particular, we find that the ratios between NH3 and CO; NH2 and CO; NH3 and HCO+ are sensitive to the evolutionary timescale, and that the ratio between HCN and OH is sensitive to the C/O ratio. Finally, we demonstrate that measurements of the central deviation (central depletion or enhancement) of abundances of certain molecules are good indicators of the dynamics of the core.

The influence of turbulent mixing on the chemistry of the interstellar medium has so far received little attention. Previous studies of this effect have suggested that it might play an important role in mixing the various phases of the interstellar medium. In this paper we examine the potential effects of turbulent diffusion on the deuterium chemistry within molecular clouds. We find that such mixing acts to reduce the efficiency of deuteration in these clouds by increasing the ionization fraction and reducing freeze-out of heavy molecules. This leads to lower abundances for many deuterated species. We also examine the influence of turbulent mixing on the transition from atomic hydrogen to H2 and from atomic deuterium to HD near the cloud edge. We find that including turbulent diffusion in our models serves to push these transitions deeper into the cloud and helps maintain a higher atomic fraction throughout the cloud envelope. Based on these findings, we propose a new process to account for the significant scatter in the observed atomic D/H ratio for galactic sightlines extending beyond the Local Bubble. Although several mechanisms have been put forward to explain this scatter, they are unable to fully account for the range in D/H values. We suggest a scenario in which turbulent mixing of atomic and molecular gas at the edges of molecular clouds causes the observed atomic D/H ratio to vary by a factor of ~2.

Aims: Little information exists on the circumstellar molecular abundance stratifications of many molecules. The aim is to study the circumstellar chemical abundance pattern of 11 molecules and isotopologs ($^{12}$CO, $^{13}$CO, SiS, $^{28}$SiO, $^{29}$SiO, $^{30}$SiO, HCN, CN, CS, SO, SO$_2$) in the oxygen-rich evolved star IK~Tau. Methods: We have performed an in-depth analysis of a large number of molecular emission lines excited in the circumstellar envelope around IK~Tau. The analysis is done based on a non-local thermodynamic equilibrium (non-LTE) radiative transfer analysis, which calculates the temperature and velocity structure in a self-consistent way. The chemical abundance pattern is coupled to theoretical outer wind model predictions including photodestruction and cosmic ray ionization. Not only the integrated line intensities, but also the line shapes, are used as diagnostic tool to study the envelope structure. Results: The deduced wind acceleration is much slower than predicted from classical theories. SiO and SiS are depleted in the envelope, possibly due to the adsorption onto dust grains. For HCN and CS a clear difference with respect to inner wind non-equilibrium predictions is found, either indicating uncertainties in the inner wind theoretical modeling or the possibility that HCN and CS (or the radical CN) participate in the dust formation. The low signal-to-noise profiles of SO and CN prohibit an accurate abundance determination; the modeling of high-excitation SO$_2$ lines is cumbersome, possibly related to line misidentifications or problems with the collisional rates. The SiO isotopic ratios ($^{29}$SiO/$^{28}$SiO and $^{30}$SiO/$^{28}$SiO) point toward an enhancement in $^{28}$SiO compared to results of classical stellar evolution codes. Predictions for H$_2$O lines in the spectral range of the Herschel/HIFI mission are performed. [abbreviated]

We present the results of models of the chemistry, including deuterium, in the inner regions of protostellar disks. We find good agreement with recent gas phase observations of several (non--deuterated) species. We also compare our results with observations of comets and find that in the absence of other processing e.g. in the accretion shock at the surface of the disk, or by mixing in the disk, the calculated D/H ratios in ices are higher than measured and reflect the D/H ratio set in the molecular cloud phase. Our models give quite different abundances and molecular distributions to other inner disk models because of the differences in physical conditions in the model disk. This emphasizes how changes in the assumptions about the density and temperature distribution can radically affect the results of chemical models.

We present mid-infrared spectral maps of the NGC 1333 star forming region, obtained with the the Infrared Spectrometer on board the Spitzer Space Telescope. Eight pure H2 rotational lines, from S (0) to S (7), are detected and mapped. The H2 emission appears to be associated with the warm gas shocked by the multiple outflows present in the region. A comparison between the observed intensities and the predictions of detailed shock models indicates that the emission arises in both slow (12 - 24 km/s) and fast (36 - 53 km/s) C-type shocks with an initial ortho-to-para ratio of ~ 1. The present H2 ortho-to-para ratio exhibits a large degree of spatial variations. In the post-shocked gas, it is usually about 2, i.e. close to the equilibrium value (~ 3). However, around at least two outflows, we observe a region with a much lower (~ 0.5) ortho-to-para ratio. This region probably corresponds to gas which has been heated-up recently by the passage of a shock front, but whose ortho-to-para has not reached equilibrium yet. This, together with the low initial ortho-to-para ratio needed to reproduce the observed emission, provide strong evidence that H2 is mostly in para form in cold molecular clouds. The H2 lines are found to contribute to 25 - 50% of the total outflow luminosity, and thus can be used to ascertain the importance of star formation feedback on the natal cloud. From these lines, we determine the outflow mass loss rate and, indirectly, the stellar infall rate, the outflow momentum and the kinetic energy injected into the cloud over the embedded phase. The latter is found to exceed the binding energy of individual cores, suggesting that outflows could be the main mechanism for core disruption.

To date, there is no core accretion simulation that can successfully account for the formation of Uranus or Neptune within the observed 2-3 Myr lifetimes of protoplanetary disks. Since solid accretion rate is directly proportional to the available planetesimal surface density, one way to speed up planet formation is to take a full accounting of all the planetesimal-forming solids present in the solar nebula. By combining a viscously evolving protostellar disk with a kinetic model of ice formation, we calculate the solid surface density in the solar nebula as a function of heliocentric distance and time. We find three effects that strongly favor giant planet formation: (1) a decretion flow that brings mass from the inner solar nebula to the giant planet-forming region, (2) recent lab results (Collings et al. 2004) showing that the ammonia and water ice lines should coincide, and (3) the presence of a substantial amount of methane ice in the trans-Saturnian region. Our results show higher solid surface densities than assumed in the core accretion models of Pollack et al. (1996) by a factor of 3 to 4 throughout the trans-Saturnian region. We also discuss the location of ice lines and their movement through the solar nebula, and provide new constraints on the possible initial disk configurations from gravitational stability arguments.

We investigate the gas-phase and grain-surface chemistry in the inner 30 AU of a typical protoplanetary disk using a new model which calculates the gas temperature by solving the gas heating and cooling balance and which has an improved treatment of the UV radiation field. We discuss inner-disk chemistry in general, obtaining excellent agreement with recent observations which have probed the material in the inner regions of protoplanetary disks. We also apply our model to study the isotopic fractionation of carbon. Results show that the fractionation ratio, 12C/13C, of the system varies with radius and height in the disk. Different behaviour is seen in the fractionation of different species. We compare our results with 12C/13C ratios in the Solar System comets, and find a stark contrast, indicative of reprocessing.

We present two new in situ core accretion simulations of Saturn with planet formation timescales of 3.37 Myr (model S0) and 3.48 Myr (model S1), consistent with observed protostellar disk lifetimes. In model S0, we assume rapid grain settling reduces opacity due to grains from full interstellar values (Podolak 2003). In model S1, we do not invoke grain settling, instead assigning full interstellar opacities to grains in the envelope. Surprisingly, the two models produce nearly identical formation timescales and core/atmosphere mass ratios. We therefore observe a new manifestation of core accretion theory: at large heliocentric distances, the solid core growth rate (limited by Keplerian orbital velocity) controls the planet formation timescale. We argue that this paradigm should apply to Uranus and Neptune as well.

We present new models of the deuterium chemistry in protoplanetary disks, including, for the first time, multiply deuterated species. We use these models to explore whether observations in combination with models can give us clues as to which desorption processes occur in disks. We find, in common with other authors, that photodesorption can allow strongly bound molecules such as HDO to exist in the gas phase in a layer above the midplane. Models including this process give the best agreement with the observations. In the midplane, cosmic ray heating can desorb weakly bound molecules such as CO and N$_2$. We find the observations suggest that N$_2$ is gaseous in this region, but that CO must be retained on the grains to account for the observed DCO$^+$/HCO$^+$. This could be achieved by CO having a higher binding energy than N$_2$ (as may be the case when these molecules are accreted onto water ice) or by a smaller cosmic ray desorption rate for CO than assumed here, as suggested by recent theoretical work. For gaseous molecules the calculated deuteration can be greatly changed by chemical processing in the disk from the input molecular cloud values. On the grains singly deuterated species tend to retain the D/H ratio set in the molecular cloud, whereas multiply deuterated species are more affected by the disk chemistry. Consequently the D/H ratios observed in comets may be partly set in the parent cloud and partly in the disk, depending on the molecule.

Benzene (c-C6H6) formation in the inner 3 AU of a protostellar disk can be efficient, resulting in high abundances of benzene in the midplane region. The formation mechanism is different to that found in interstellar clouds and in protoplanetary nebulae, and proceeds mainly through the reaction between allene (C3H4) and its ion. This has implications for PAH formation, in that some fraction of PAHs seen in the solar system could be native rather than inherited from the interstellar medium.

The dynamics and chemistry of protostellar disks are likely to be intricately linked, with dynamical processes altering the chemical composition, and chemistry, in turn, controlling the ionization structure and hence the ability of the magneto-rotational instability to drive the disk turbulence. Here we present the results from the first chemical models of the outer regions (R > 100 AU) of protoplanetary disks to consider the effects of turbulence driven diffusive mixing in the vertical direction. We show that vertical diffusion can greatly affect the column densities of many species, increasing them by factors of up to two orders of magnitude. Previous disk models have shown that disks can be divided into three chemically distinct layers, with the bulk of the observed molecular emission coming from a region between an atomic/ionic layer on the surface of the disk and the midplane regoin where the bulk of molecules are frozen onto grains. Diffusion retains this three layer structure, but increases the depth of the molecular layer by bringing atoms and atomic ions form by photodissociation in the surface layers into the shielded molecular layer where molecules can reform. For other species, notably NH3 and N2H+, the column densities are relatively unaffected by diffusion. These species peak in abundance near the midplane where most other molecules are heavily depleted, rather than in the molecular layer above. Diffusion only affects the abundances of those molecules with peak abundances in the molecular layer. We find that diffusion does not affect the ionization fraction of the disk. We compare the calculated column densities to observations of DM Tau, LkCa 15 and TW Hya and find good agreement for many molecules with a diffusion coefficient of 1e18 cm^2 s^-1.

The effects of turbulence on the mixing of gases and dust in the outer Solar nebula are examined using 3-D MHD calculations in the shearing-box approximation with vertical stratification. The turbulence is driven by the magneto-rotational instability. The magnetic and hydrodynamic stresses in the turbulence correspond to an accretion time at the midplane about equal to the lifetimes of T Tauri disks, while accretion in the surface layers is thirty times faster. The mixing resulting from the turbulence is also fastest in the surface layers. The mixing rate is similar to the rate of radial exchange of orbital angular momentum, so that the Schmidt number is near unity. The vertical spreading of a trace species is well-matched by solutions of a damped wave equation when the flow is horizontally-averaged. The damped wave description can be used to inexpensively treat mixing in 1-D chemical models. However, even in calculations reaching a statistical steady state, the concentration at any given time varies substantially over horizontal planes, due to fluctuations in the rate and direction of the transport. In addition to mixing species that are formed under widely varying conditions, the turbulence intermittently forces the nebula away from local chemical equilibrium. The different transport rates in the surface layers and interior may affect estimates of the grain evolution and molecular abundances during the formation of the Solar system.

The determination of the fractional electron abundance, Xe, in protostellar cores relies on observations of molecules, such as DCO+, H13CO+ and CO, and on chemical models to interpret their abundance. Studies of protostellar cores have revealed significant variations of Xe from core to core within a range 10^-8<Xe<10^-6. The physical origin of these large variations in Xe is not well understood, unless unlikely variations in the cosmic ray ionization rate or ad hoc values of metal depletion are assumed. In this work we explore other potential causes of these variations in Xe, including core age, extinction and density. We compute numerically the intensity of the radiation field within a density distribution generated by supersonic turbulence. Taking into account the lines of sight in all directions, the effective visual extinction in dense regions is found to be always much lower than the extinction derived from the column density along a fixed line of sight. Dense cores with volume and column densities comparable to observed protostellar cores have relatively low mass-averaged visual extinction, 2mag <= A_V <= 5mag, such that photo-ionization can sometimes be as important as cosmic ray ionization. Chemical models, including gas-grain chemistry and time dependent gas depletion and desorption, are computed for values of visual extinction in the range 2mag <= A_V <= 6mag, and for a hydrogen gas density of 10^4cm^-3$, typical of protostellar cores. The models presented here can reproduce some of the observed variations of ion abundance from core to core as the combined effect of visual extinction and age variations. The range of electron abundances predicted by the models is relatively insenstive to density over 10^4 to 10^6 cm^{-3}.

The gas-phase non-equilibrium chemistry of the inner wind of the oxygen-rich Mira variable IK Tau (NML Tau) is investigated using a physio-chemical model describing the periodically shocked gas close to the stellar photosphere. We predict the formation in large amounts of a few carbon species in the inner envelope, in particular CO$_2$, which has been recently detected with ISO in the spectra of several oxygen-rich semi-regular and Mira stars. The theoretical abundances are also in excellent agreement with values derived from millimeter and sub-millimeter observations, pointing to the fact that some carbon species in oxygen stars do form from shock chemistry in the inner layers and travel the envelope as ``parent'' species.