Olli Sipilä

Astronomy & Astrophysics (A&A) is an international journal which publishes papers on all aspects of astronomy and astrophysics

Aims: We aim to investigate the abundances of light deuterium-bearing species such as HD, H2D+ and D2H+ in a gas-grain chemical model including an extensive description of deuterium and spin state chemistry, in physical conditions appropriate to the very centers of starless cores. Methods: We combine a gas-grain chemical model with radiative transfer calculations to simulate density and temperature structure in starless cores. The chemical model includes deuterated forms of species with up to 4 atoms and the spin states of the light species H2, H2+ and H3+ and their deuterated forms. Results: We find that HD eventually depletes from the gas phase because deuterium is efficiently incorporated to grain-surface HDO, resulting in inefficient HD production on grains. HD depletion has consequences not only on the abundances of e.g. H2D+ and D2H+, whose production depends on the abundance of HD, but also on the spin state abundance ratios of the various light species, when compared with the complete depletion model where heavy elements do not influence the chemistry. Conclusions: While the eventual HD depletion leads to the disappearance of light deuterium-bearing species from the gas phase in a relatively short timescale at high density, we find that at late stages of core evolution the abundances of H2D+ and D2H+ increase toward the core edge and the disributions become extended. The HD depletion timescale increases if less oxygen is initially present in the gas phase, owing to chemical interaction between the gas and the dust predecing the starless core phase. Our results are greatly affected if H2 is allowed to tunnel on grain surfaces, and therefore more experimental data not only on tunneling but also on the O + H2 surface reaction in particular is needed.

Astronomy & Astrophysics (A&A) is an international journal which publishes papers on all aspects of astronomy and astrophysics

Aims: We aim to simulate radial profiles of molecular abundances and the gas temperature in cold and heavily shielded starless cores by combining chemical and radiative transfer models. Methods: A determination of the dust temperature in a modified Bonnor-Ebert sphere is used to calculate initial radial molecular abundance profiles. The abundances of selected cooling molecules corresponding to two different core ages are then extracted to determine the gas temperature at two time steps. The calculation is repeated in an iterative process yielding molecular abundances consistent with the gas temperature. Line emission profiles for selected substances are calculated using simulated abundance profiles. Results: The gas temperature is a function of time; the gas heats up as the core gets older because the cooling molecules are depleted onto grain surfaces. The contributions of the various cooling molecules to the total cooling power change with time. Radial chemical abundance profiles are non-trivial: different species present varying degrees of depletion and in some cases inward-increasing abundances profiles, even at t > 10^5 years. Line emission simulations indicate that cores of different ages can present significantly different line emission profiles, depending on the tracer species considered. Conclusions: Chemical abundances and the associated line cooling power change as a function of time. Most chemical species are depleted onto grain surfaces at densities exceeding ~10^5 cm^-3. Notable exceptions are NH_3 and N2H^+; the latter is largely undepleted even at n_H~10^6 cm-3. On the other hand, chemical abundances are not significantly developed in regions of low gas density even at t~10^5 years, revealed by inward-increasing abundance gradients. The gas temperature can be significantly different from the dust temperature; this may have implications on core stability.

Astronomy & Astrophysics (A&A) is an international journal which publishes papers on all aspects of astronomy and astrophysics

Aims: We aim to derive a stability condition for non-isothermal Bonnor-Ebert spheres and compare the physical properties of critical non-isothermal and isothermal gas spheres. These configurations can serve as models for prestellar cores before gravitational collapse. Methods: A stability condition for non-isothermal spheres is derived by constructing an expression for the derivative of boundary pressure with respect to core volume. The temperature distribution is determined by means of radiative transfer calculations. Based on the stability analysis, we derive the physical parameters of critical cores for the mass range 0.1 - 5.0 M_sun. In addition, the properties of roughly Jupiter-mass cores are briefly examined. Results: At the low-mass end the critical non-isothermal sphere has lower central density and a slightly larger physical radius than the corresponding isothermal sphere (i.e. one with the same mass and average temperature). The temperature decrease towards the core centre becomes steeper towards smaller masses as the central density becomes higher. The slope depends on the adopted dust model. We find that the critical dimensionless radius increases above the isothermal value xi_0 = 6.45 for very low-mass cores (< 0.2 M_sun). However, in the mass-range studied here the changes are within 5% from the isothermal value. Conclusions: The density structures of non-isothermal and isothermal Bonnor-Ebert spheres for a given mass are fairly similar. However, the present models predict clear differences in the average temperatures for the same physical radius. Especially for low-mass cores, the temperature gradient probably has implications on the chemistry and the observed line emission. We also find that hydrostatic Jupiter-mass cores with radii less than 100 AU would have very high boundary pressures compared with typical pressures in the interstellar space.

Astronomy & Astrophysics (A&A) is an international journal which publishes papers on all aspects of astronomy and astrophysics

Context: The depletion of heavy elements in cold cores of interstellar molecular clouds can lead to a situation where deuterated forms of H_3^+ are the most useful spectroscopic probes of the physical conditions. Aims: The aim is to predict the observability of the rotational lines of H_2D^+ and D_2H^+ from prestellar cores. Methods: Recently derived rate coefficients for the H_3^+ + H_2 isotopic system were applied to the "complete depletion" reaction scheme to calculate abundance profiles in hydrostatic core models. The ground-state lines of H_2D^+(o) (372 GHz) and D_2H^+(p) (692 GHz) arising from these cores were simulated. The excitation of the rotational levels of these molecules was approximated by using the state-to-state coefficients for collisions with H_2. We also predicted line profiles from cores with a power-law density distribution advocated in some previous studies. Results: The new rate coefficients introduce some changes to the complete depletion model, but do not alter the general tendencies. One of the modifications with respect to the previous results is the increase of the D_3^+ abundance at the cost of other isotopologues. Furthermore, the present model predicts a lower H_2D^+ (o/p) ratio, and a slightly higher D_2H^+ (p/o) ratio in very cold, dense cores, as compared with previous modelling results. These nuclear spin ratios affect the detectability of the submm lines of H_2D^+(o) and D_2H^+(p). The previously detected H_2D^+ and D_2H^+ lines towards the core I16293E, and the H_2D^+ line observed towards Oph D can be reproduced using the present excitation model and the physical models suggested in the original papers.