A numerical approach to model chemistry of complex organic molecules in a protoplanetary disk

Abstract Multiphase astrochemical modeling presents a numerical challenge especially for the simulation of objects with the wide range of physical parameters such as protoplanetary disks. We demonstrate an implementation of the analytical Jacobian for the numerical integration of the system of differential rate equations that govern chemical evolution in star-forming regions. The analytical Jacobian allowed us to greatly improve the stability of the code in protoplanetary disk conditions. We utilize the MONACO code to study the evolution of abundances of chemical species in protoplanetary disks. The chemical model includes 670 species and 6,015 reactions in the gas phase and on interstellar grains. The specific feature of the utilized chemical model is the inclusion of low-temperature chemical processes leading to the formation of complex organic molecules (COMs), included previously in the models of chemistry of COMs in prestellar clouds. To test the impact of analytical Jacobian on the stability of numerical simulations of chemical evolution in protoplanetary disks, we calculated the chemical composition of the disk using a two-phase model and four variants of the chemical reaction network, three values of the surface diffusion rates, and two types of the initial chemical composition. We also show a preliminary implementation of the analytical Jacobian to a three-phase model.