Modelling methanol and hydride formation in the JWST Ice Age era

• DOI: 10.48550/arxiv.2502.10123
• Publication Date: 2025-01-01
• Authors: Jiménez-Serra, Iz...; Megías, Andrés; Salaris, Joseph; Cuppen, Herma; Taillard, Angèle; Jin, Miwha; Wakelam, Valentine; Vasyunin, Anton I.; Caselli, Paola; Pendleton, Yvonne J.; Dartois, Emmanuel; Noble, Jennifer A.; Viti, Serena; Borshcheva, Katerina; Garrod, Robin T.; Lamberts, Thanja; Fraser, Helen; Melnick, Gary; McClure, Melissa; Rocha, Will; Drozdovskaya, Mar...; Lis, Dariusz C.
• Keywords: Astrophysics - Solar and Stellar Astrophysics; Astrophysics of Galaxies (astro-ph.GA); FOS: Physical sciences; Astrophysics - Astrophysics of Galaxies; Solar and Stellar Astrophysics (astro-ph.SR)

Accepted in A&A<br/>(Abridged) JWST observations have measured the ice composition toward two highly-extinguished field stars in the Chamaeleon I cloud. The observed extinction excess on the long-wavelength side of the H2O ice band at 3 micron has been attributed to a mixture of CH3OH with ammonia hydrates, which suggests that CH3OH ice could have formed in a water-rich environment with little CO depletion. Laboratory experiments and quantum chemical calculations suggest that CH3OH could form via the grain surface reactions CH3+OH and/or C+H2O in water-rich ices. However, no dedicated chemical modelling has been carried out thus far to test their efficiency and dependence on the astrochemical code employed. We model the ice chemistry in the Chamaeleon I cloud using a set of astrochemical codes (MAGICKAL, MONACO, Nautilus, UCLCHEM, and KMC simulations) to test the effects of the different code architectures and of the assumed ice chemistry. Our models show that the JWST ice observations are better reproduced for gas densities >1e5 cm-3 and collapse times >1e5 yr. CH3OH ice forms predominantly (>99%) via CO hydrogenation. The contribution of reactions CH3+OH and C+H2O, is negligible. The CO2 ice may form either via CO+OH or CO+O depending on the code. However, KMC simulations reveal that both mechanisms are efficient despite the low rate constant of the CO+O surface reaction. CH4 is largely underproduced for all codes except for UCLCHEM, for which a higher amount of atomic C is available during the initial translucent cloud phase. Large differences in the ice abundances are found at Tdust<12 K between diffusive and non-diffusive chemistry codes. This is due to the fact that non-diffusive chemistry takes over diffusive chemistry at such low Tdust. This could explain the rather constant ice chemical composition found in Chamaeleon I and other dense cores despite the different visual extinctions probed.<br/>