Combined Orbital and Thermal Evolution of Oort Cloud Comets

Abstract: We present a fully integrated model of comet evolution that couples thermal and compositional processes with dynamical processes continuously, from formation to present-day activity. The combined code takes into account changes in orbital parameters that define the heliocentric distance as a function of time, which is fed into the thermal/compositional evolution code. The latter includes a set of volatile species, gas flow through the porous interior, crystallization of amorphous ice, sublimation and refreezing of volatiles in the pores. We follow the evolution of three models, with radii of 2, 10 and 50 km for 4.6 Gyr, through different dynamical epochs, starting in the vicinity of Neptune, moving to the Oort Cloud and after a long sojourn there, back inward to the planetary region. The initial composition includes a mixture of CO, CO2 ices, amorphous water ice with trapped CO and CO2, and dust.We find that the CO ice is completely depleted in the small object, but preserved in the larger ones from a depth of 500 m to the center, while the CO2 and the amorphous ice are entirely preserved. Of crucial importance is the change in CO abundance profiles during the cooling phase, as the objects migrate to the OC. Upon return from the Oort Cloud, the activity is driven by CO sublimation at large heliocentric distances (up to 50 au), by CO2 inward of 13 au and by gas released from crystallizing amorphous ice at about 7 au. We test the effect of radioactive heating by long-lived isotopes and find that it is negligible. Considering sub-solar temperatures and limited active areas, we show that CO2 production rates can exceed the detection limit as far out as 25 au.