Impact of ice growth on the physical and chemical properties of dense cloud cores

We investigated the effect of time-dependent ice growth on dust grains on the opacity and hence on the dust temperature in a collapsing molecular cloud core, with the aim of quantifying the effect of the dust temperature variations on ice abundances as well as the evolution of the collapse. To perform the simulations, we employed a one-dimensional collapse model that self-consistently and time-dependently combines hydrodynamics with chemical and radiative transfer simulations. The dust opacity was updated on the fly based on the ice growth as a function of the location in the core. The results of the fully dynamical model were compared against simulations run with different values of fixed ice thickness. We found that the ice thickness increases quickly and reaches a saturation value (as a result of a balance between adsorption and desorption) of approximately 90 monolayers in the central core (volume density ~104 cm−3), and several tens of monolayers at a volume density of ~103 cm−3, after only a few 105 yr of evolution. The results thus exclude the adoption of thin (approximately ten monolayer) ices in molecular cloud simulations except at very short timescales. However, the differences in abundances and the dust temperature between the fully dynamic simulation and those with a fixed dust opacity are small; abundances change between the solutions generally within a factor of two only. The assumptions on the dust opacity do have an effect on the collapse dynamics through the influence of the photoelectric effect on the gas temperature, and the simulations take a different time to reach a common central density. This effect is, however, small as well. In conclusion, carrying out chemical simulations using a dust temperature corresponding to a fixed opacity seems to be a good approximation. Still, although at least in the present case its effect on the overall results is limited – as long as the grains are monodisperse – ice growth should be considered to obtain the most accurate representation of the collapse dynamics. We have found in a previous work that considering a grain size distribution leads to a complicated ice composition that depends on the grain size nonlinearly. With this in mind, we will carry out a follow-up study where the influence of the grain size on the present simulation setup is investigated.