Surface Of Comet Research Articles

Electron‐Induced Radiolysis of Water Ice and the Buildup of Oxygen

AbstractIrradiation by energetic ions, electrons, and UV photons induces sputtering and chemical processes (radiolysis) in the surfaces of icy moons, comets, and icy grains. Laboratory experiments, both of ideal surfaces and of more complex and realistic analog samples, are crucial to understand the interaction of surfaces of icy moons and comets with their space environment. This study shows the first results of mass spectrometry measurements from porous water ice regolith samples irradiated with electrons as a representative analogy to water‐ice rich surfaces in the solar system. Previous studies have shown that most electron‐induced O radiolysis products leave the ice as and and that can be trapped under certain conditions in the irradiated ice. Our new laboratory experiments confirm these findings. Moreover, they quantify residence times and saturation levels of in originally pure water ice. O may also be released from the water ice by irradiation, but the quantification of the released O is more difficult and the total amount is sensitive to the electron flux and energy.

ALMA and GMRT Studies of Dust Continuum Emission and Spectral Lines Toward Oort Cloud Comet C/2022 E3 (ZTF)

The atomic and molecular compounds of cometary ices serve as valuable knowledge into the chemical and physical properties of the outer solar nebula, where comets are formed. From the cometary atmospheres, the atoms and gas-phase molecules arise mainly in three ways: (i) the outgassing from the nucleus, (ii) the photochemical process, and (iii) the sublimation of icy grains from the nucleus. In this paper, we present the radio and millimeter wavelength observation results of Oort cloud non-periodic comet C/2022 E3 (ZTF) using the Giant Metrewave Radio Telescope (GMRT) band L and the Atacama Large Millimeter/submillimeter Array (ALMA) band 6. We do not detect continuum emissions and an emission line of atomic hydrogen (H I) at rest frequency 1420 MHz from this comet using the GMRT. Based on ALMA observations, we detect the dust continuum emission and rotational emission lines of methanol (CH3OH) from comet C/2022 E3 (ZTF). From the dust continuum emission, the dust production (Afρ) activity of comet ZTF is 2280 ± 50 cm. Based on LTE spectral modeling, the column density and excitation temperature of CH3OH toward C/2022 E3 (ZTF) are (4.50 ± 0.25) × 1014 cm−2 and 70 ± 3 K respectively. The integrated emission maps show that CH3OH was emitted from the coma region of the comet. The production rate of CH3OH toward C/2022 E3 (ZTF) is (7.32 ± 0.64) × 1026 molecules s−1. The fractional abundance of CH3OH with respect to H2O in the coma of the comet is 1.52%. We also compare our derived abundance of CH3OH with the existence modeled value, and we see the observed and modeled values are quite similar. We claim that CH3OH is formed via the subsequential hydrogenation of formaldehyde (H2CO) on the grain surface of comet C/2022 E3 (ZTF).

Grain polydispersity and non-sphericity effects on gas flow through granular beds using measurements and modelling

ABSTRACT The quality of cometary surface activity simulations and erosion models of icy moons depends on a good knowledge of the surface layer permeability to gas flow. Therefore, we study various models of the Knudsen diffusion coefficient and the viscous permeability, which are used to describe the flow of rarefied gases through porous materials. Usually, these models are expressed for monodisperse packed beds. In this work, we describe a generalization to polydisperse packed beds and compare them with experimental results. In addition, we analyse non-spherical packings to test how well the recently developed models are applicable. For this purpose, the gas flow parameters of these samples are measured in a dedicated measurement set-up. Special attention had to be paid to biases in measuring the porosity and the pressure drop in the sample, which are discussed in detail. Our measurements confirm that the Knudsen diffusion coefficient is inversely proportional to the specific surface area of the grains and that the viscous permeability is inversely proportional to the specific surface area squared. Further, we were able to identify a relation between the gas flow parameters, represented by a parameter β, which seems to be an indicator of the mean orientation of the grains. The findings give further evidence of the importance of the grain size distribution and the grain shape for rarefied gas flow. In particular, the results show that the standard polydisperse model is not sufficient when a considerable part of the pore space consists of traps or other rarely percolated parts.

Ejection and dynamics of aggregates in the coma of comet 67P/Churyumov-Gerasimenko

The process of gas-driven ejection of refractory materials from cometary surfaces continues to pose a challenging question in cometary science. The activity modeling of comet 67P/Churyumov-Gerasimenko, based on data from the Rosetta mission, has significantly enhanced our comprehension of cometary activity. But thermophysical models have difficulties in simultaneously explaining the production rates of various gas species and dust. It has been suggested that different gas species might be responsible for the ejection of refractory material in distinct size ranges. This work focuses on investigating the abundance and the ejection mechanisms of large aggregates (gtrsim 1 cm) from the comet nucleus. We aim to determine their properties and map the distribution of their source regions across the comet surface. This can place constraints on activity models for comets. We examined 189 images acquired at five epochs by the OSIRIS/NAC instrument on board the Rosetta spacecraft. Our goal was to identify bright tracks produced by individual aggregates as they traversed the camera field of view. In parallel, we generated synthetic images based on the output of dynamical simulations involving various types of aggregates. By comparing these synthetic images with the observations, we determined the characteristics of the simulated aggregates that most closely resemble the observations. We have identified over 30000 tracks present in the OSIRIS images, derived constraints on the characteristics of the aggregates, and mapped their origins on the nucleus surface. The aggregates have an average radius of $ cm and a bulk density consistent with that of the comet's nucleus. Due to their size, gas drag exerts only a minor influence on their dynamical behavior, so an initial velocity is needed to bring them into the camera field of view. The source regions of these aggregates are predominantly located near the boundaries of distinct terrains on the surface.

New measurements of the hemispherical albedo of cometary analogues: Influence on the amplitude of the change of cometary brightness during its outburst

The paper presents new experimental results of measurements of the hemispherical albedo of a potential comet analogue and an analysis of its influence on the change in the cometary brightness. Measurements of the spectral properties of the tested analogues were performed using a Cary 5000 spectrometer in the wavelength range from 200 to 2500 nm. Two samples were used as a dust analogue: a mixture of sand and ash (sample A) and a mixture of sand and charcoal (sample B), respectively. The obtained hemispherical albedo values for the tested samples ranged from 1.20% to 32.25%. In particular, for sample B, the obtained hemispherical albedo values are comparable to the actual reflectance for comet 67P. Based on the obtained spectral properties of the tested samples, a thermodynamic model of cometary activity was developed and used, among others, to calculate the power absorbed by the tested analogue and the sublimation flux occurring through the porous structure of the fine-grained dust present on the cometary surface. The obtained results were used to calculate the change in the cometary brightness during its outburst, the upper limit of which was equal to 5.36 magnitude.

Localised sublimation under cometary dust mantles

Observations of cometary surfaces have indicated that the amount of exposed surficial water ice is insufficient to supply the measured water production rates within the inner solar system. The measured water production rates are also a small fraction of that which would arise from freely subliming pure water ice bodies of similar surface area. A plausible mechanism, given the now well-established, high porosity of comets is that water ice is subliming below the visible surface with the produced gas escaping through the pores in a non-subliming hot dust matrix. If this process leads to an additional heating of the evolving gas, this may provide an elegant explanation for the temperature of the gas required to match measured terminal velocities. We investigate this process using a series of numerical models based on the Direct Simulation Monte Carlo (DSMC) method for rarified gas dynamics. We show that inhomogeneities in the source are not evident in the flow field once the thickness of the porous layer reaches a specific size and that lateral flow through the porous layer (parallel to the surface) can broaden the apparent source size when the source is small and isolated. The temperatures of the gas flow at the visible surface can be significantly increased and will respond to the temperature distribution within the non-volatile porous layer. The influence of the structural properties of the porous layer (how the porous layer is mechanically constructed) is shown to have only limited influence on the final gas parameters at the surface.

Aligned fractures on asteroid Ryugu as an indicator of thermal fracturing

Context. Asteroid and comet surfaces are exposed to a complex environment that includes low gravity, high temperature gradients, and a bombardment of micrometeorites and cosmic rays. Surface material exposed to this environment evolves in a specific way depending on various factors such as the bodies’ size, heliocentric distance, and composition. Fractures in boulders, as seen on asteroid Ryugu, can help to determine and constrain the dominant processes eroding small-body surface materials. It is also possible to estimate fracture growth timescales based on the abundance and length of fractures in boulders. Aims. We analyse the number, orientation, and length of fractures on asteroid Ryugu to establish the relation between the fractures and the processes that may have formed them. We also compare our results to similar investigations conducted on other small bodies and estimate the timescale of fracture growth. Methods. 198 high-resolution Hayabusa2 images of asteroid Ryugu suitable for our fracture analysis were selected and map-projected. Within these images, fractures in boulders were manually mapped using the QGIS software. The fracture coordinates were extracted and the fractures’ orientation and length were computed for 1521 identified fractures. Results. Fractures in boulders on asteroid Ryugu are found to be preferentially north-south aligned, suggesting a formation through thermal erosion. Modeling the fracture length indicates a fracture growth timescale of 30 000 to 40 000 yr, slightly younger than ages found previously for asteroid Bennu. The errors in these ages, due to uncertainties about the thermophysical parameters used in this model, are substantial (−33 000 yr +250 000 yr). However, even with these large errors, the model suggests that thermal fracturing is a geologically fast process. These times are not too dissimilar to those quoted in the literature for Ryugu and Bennu, since similar thermophysical material parameters for Ryugu and Bennu seem likely.

Cometary outbursts in the Oort cloud

Comet nuclei in the outer Solar system are constantly irradiated by cosmic rays at low temperatures. Accumulated high concentrations of radicals can undergo fast recombination with significant heating of cometary surface layers. We present the model of comet activity at large heliocentric distances caused by the recombination of radicals. We found that the considered mechanism can cause activity of comets in distant regions of the Solar system, even at the Oort cloud distances. Outbursts in distant comet reservoirs can be a new source of dust and ice particles contributing to the recently discovered anomalous diffuse light in the cosmic extragalactic background optic light and the unexpectedly high flux of dust particles detected by the New Horizons dust counter at the edge of the Kuiper belt. The orbits of small-radii comets in the Oort cloud are highly influenced by cometary outbursts. This effect may account for the observed decrease in the number of small-radius long-period comets.

Cold electrons at a weakly outgassing comet

ABSTRACT Throughout the Rosetta mission, cold electrons (<1 eV) were measured in the coma of comet 67P/Churyumov–Gerasimenko. Cometary electrons are produced at ∼10 eV through photoionization or through electron-impact ionization collisions. The cold electron population is formed by cooling the warm population through inelastic electron–neutral collisions. Assuming radial electron outflow, electrons are collisional with the neutral gas coma below the electron exobase, which only formed above the comet surface in near-perihelion high-outgassing conditions (Q > 3 × 1027 s−1). However, the cold population was identified at low outgassing (Q < 1026 s−1), when the inner coma was not expected to be collisional. We examine cooling of electrons at a weakly outgassing comet, using a 3D collisional model of electrons at a comet. Electron paths are extended by trapping in an ambipolar electric field and by gyration around magnetic field lines. This increases the probability of electrons undergoing inelastic collisions with the coma and becoming cold. We demonstrate that a cold electron population can be formed and sustained, under weak outgassing conditions (Q = 1026 s−1), once 3D electron dynamics are accounted for. Cold electrons are produced in the inner coma through electron–neutral collisions and transported tailwards by an E × B drift We quantify the efficiency of trapping in driving electron cooling, with trajectories typically 100 times longer than expected from ballistic radial outflow. Based on collisional simulations, we define an estimate for a region where a cold electron population can form, bounded by an electron cooling exobase. This estimate agrees well with cold electron measurements from the Rosetta Plasma Consortium.

Upper limit of the solar wind protons backscattering efficiency from Comet 67P/Churyumov-Gerasimenko

Context. Solar wind ions backscattering is a fundamental plasma-surface interaction process that may occur on all celestial bodies exposed to the solar wind and lacking a significant atmosphere or magnetosphere. Yet, observations have been limited to the regolith-covered Moon and Phobos, one of the Martian moons. Aims. We aim to expand our knowledge of the process to include comets by investigating the backscattering of solar wind protons from the surface of comet 67P/Churyumov-Gerasimenko. Methods. We used one of the ion spectrometers on board ESA’s Rosetta spacecraft to search for evidence of backscattered solar wind protons from the cometary surface. The signal of interest was expected to be very weak and several statistical treatments of the data were essential to eliminate any influence from background noise and instrumental effects. Due to limited knowledge of the signal location within the observed parameter space, we conducted a statistical analysis to identify the most probable conditions for detecting the signal. Results. No significant solar wind backscattered protons were ever observed by the instrument. The statement applies to the large spectrum of observation conditions. An upper limit of the backscattered proton flux is given, as well as an upper limit of the backscattering efficiency of 9 × 10−4. Conclusions. The surface of comet 67P/Churyumov-Gerasimenko distinguishes itself as a notably weak reflector of solar wind protons, with its backscattering efficiency, at most, as large as the lowest observed backscattering efficiency from the lunar regolith.

Preservation of pristine materials under impact craters formed on comet nuclei

Comet nuclei are the most primitive small bodies in the solar system, and would contain highly volatiles and amorphous water ice. However, comet nuclei may lose their primordial properties after high-velocity impacts with other small bodies. Such high-velocity impacts might also cause the characteristic circular structures and pinnacles on the comet-nucleus surface. However, the true origin of these geological structure remains uncertain. We here sought to clarify the shape of impact craters and the effects of post-shock heating on the comet surfaces. We conducted high-velocity impact experiments at velocities of 3–5.8 km s−1 on porous-ice targets with 40%, 50% and 60% porosity, and measured the post-shock temperature around the impact crater. The hard, dense refrozen layer formed on the crater wall resembled the pinnacles and sharp pit rims observed on comet nuclei, suggesting that such structures on comets may originate from impacts. In addition, the crater size scaling law for porous ice targets was successfully obtained in the strength-dominated regime. From the temperature profiles obtained, the energy partition coefficient of the projectile kinetic energy to the post-shock heat, γ, was estimated as follows: γ=0.08, 0.27 and 0.63 for 40%, 50% and 60% porosity at 4.2 km s−1 impact velocity, and γ=0.58 and 0.22 for 50% porosity at 3 km s−1 and 5.8 km s−1 impact velocities. These results, together with our crater size scaling law for porous ice, were used to estimate the temperature distribution induced by post-shock heat on the comet-nucleus surface. Volatile depletion and crystallization of amorphous water ice can occur at distances less than twice the crater radius from the crater center—that is, the post-shock heat effects are confined to the area around the crater. This might suggest that the comet nucleus surface remains relatively primitive after impacts.

Spectrophotometric properties of CoPhyLab’s dust mixtures

ABSTRACT Objective: In the framework of the Cometary Physics Laboratory (CoPhyLab) and its sublimation experiments of cometary surface analogues under simulated space conditions, we characterize the properties of intimate mixtures of juniper charcoal and SiO2 chosen as a dust analogue (Lethuillier et al. 2022). We present the details of these investigations for the spectrophotometric properties of the samples. Methods: We measured these properties using a hyperspectral imager and a radio-goniometer. From the samples’ spectra, we evaluated reflectance ratios and spectral slopes. From the measured phase curves, we inverted a photometric model for all samples. Complementary characterizations were obtained using a pycnometer, a scanning electron microscope and an organic elemental analyser. Results: We report the first values for the apparent porosity, elemental composition, and VIS-NIR spectrophotometric properties for juniper charcoal, as well as for intimate mixtures of this charcoal with the SiO2. We find that the juniper charcoal drives the spectrophotometric properties of the intimate mixtures and that its strong absorbance is consistent with its elemental composition. We find that SiO2 particles form large and compact agglomerates in every mixture imaged with the electron microscope, and its spectrophotometric properties are affected by such features and their particle-size distribution. We compare our results to the current literature on comets and other small Solar system bodies and find that most of the characterized properties of the dust analogue are comparable to some extent with the spacecraft-visited cometary nuclei, as well as to Centaurs, Trojans, and the bluest TNOs.

The influence of hemispherical albedo on the dynamics of cometary particles

ABSTRACT This paper presents a new approach to determining the dimensions of particles that can be lifted to a coma by the sublimation flux of water ice. Due to the diverse colour of the cometary surfaces, in the presented calculations we take into account the influence of the reflectance coefficient (hemispherical albedo) on the size of particles lifted from the surface of the nucleus. A Cary 5000 spectrometer was used to measure hemispherical albedo. Spectroscopic measurements used dust particles that consisted of pure silicate or were a mixture of pure silicate and organic compounds such as charcoal, soot, and ash. In addition, these admixtures acted as a means of darkening the tested sample. Depending on the colour of the sample and the wavelength, the measured value of hemispherical albedo oscillated in the range from 4.93 per cent to 48.65 per cent. A thermodynamic model of cometary matter emission was developed based on the results of hemispherical albedo measurements. The performed numerical simulations show that the temperature decreases with the increase of the hemispherical albedo, which translates into the sublimation flux and the dynamics of the dust-ice particles present on the surface of the cometary nucleus. Taking into account the extreme values of the hemispherical albedo, the differences in the size of particles carried away from the surface of the nucleus are in the order of centimeters.

Pits on Jupiter-family Comets and the Age of Cometary Surfaces

Large and deep depressions, also known as pits, are observed at the surface of all Jupiter-family comets (JFCs) imaged by spacecraft missions. They offer the opportunity to glimpse the subsurface characteristics of comet nuclei and study the complex interplay between surface structures and cometary activity. This work investigates the evolution of pits at the surface of 81P/Wild 2, 9P/Tempel 1, and 103P/Hartley 2, in continuation of the work by Benseguane et al. on 67P/Churyumov–Gerasimenko. Pits are selected across the surface of each nucleus, and high-resolution shape models are used to compute the energy they receive. A thermal evolution model is applied to constrain how cometary activity sustained under current illumination conditions could modify them. Similar to what was found for 67P, we show that erosion resulting from water-driven activity is primarily controlled by seasonal patterns that are unique to each comet as a consequence of their shape and rotational properties. However, progressive erosion sustained after multiple perihelion passages is not able to carve any of the observed pits. Instead, cometary activity tends to erase sharp morphological features; they become wider and shallower over time. Our results reinforce the evolutionary sequence evidenced from independent measurables to transform “young” cometary surfaces, with sharp surface topography prone to outbursts, into “old” cometary surfaces. Finally, we suggest that the mechanism at the origin of the pits on JFCs should be able to carve these structures in a region of the solar system where water ice does not sublimate; the Centaur phase thus appears critical to understand JFC surface properties.

Non-gravitational Mechanism of Comets’ Ejection from the Oort Cloud Due to Cometary Outbursts

Cometary nuclei located in the Oort cloud accumulate high concentration of radicals in surface layers under cosmic ray irradiation at low temperatures. Recombination of radicals induced by an increase in the surface temperature of a comet by a close passing star, O/B stars, or nearby supernovae leads to the heating of the ice layer with the releasing of volatiles from the amorphous ice. When high gas pressure builds up beneath the cometary surface, dust and gas are ejected. The resulting jet of gas and dust can change the comet’s orbit in the Oort cloud. The studied non-gravitational mechanism can effectively expel comets with a radius of ≤1 km from the Oort cloud into the inner part of the Solar system. The total effect of cometary outbursts on the stability of cometary orbits during the evolution of Solar system can result in a decrease in the number of long-period small-radius comets.

Main Results from the ISSI International Team “Characterization of 67P Cometary Activity”

The ESA/Rosetta mission accompanied the Jupiter Family Comet 67P/Churyumov-Gerasimenko and provided a huge amount of data which are providing important results about cometary activity mechanisms. We summarize the results obtained within the ISSI International Team Characterization of 67P cometary activity, which studied dust and gas ejection in different stages of the comet’s orbit, by means of a data fusion between instruments onboard the Rosetta orbiter, i.e., the OSIRIS camera, the VIRTIS imaging spectrometer, the GIADA dust detector, the MIDAS atomic force microscope, the COSIMA dust mass spectrometer, and the ROSINA gas mass spectrometer, supported by numerical models and experimental work. The team reconstructed the motion of the dust particles ejected from the comet surface, finding a correlation between dust ejection and solar illumination as well as larger occurrence of fluffy (pristine) particles in less processed and more pebble-rich terrains. Dust activity is larger in ice-rich terrains, indicating that water sublimation is the dominant activity process during the perihelion phase. The comparison of dust fluxes of different particle size suggests a link between dust morphology and ejection speed, generation of micrometric dust from fragmentation of millimetric dust, and homogeneity of physical properties of compact dust particles across the 67P surface. The comparison of fluxes of refractory and ice particles suggests the occurrence of a small amount of ice in fluffy particles, which is released when they are fragmented. A new model of cometary activity has been finally developed, according to which the comet nucleus includes Water-Ice-Enriched Blocks (WEBs), that, when exposed by CO2 activity, are the main sources of water sublimation and dust ejection.

Chemical evolution of electron-bombarded crystalline water ices at different temperatures using the procoda code

Abstract Water ices are a common component of cold space environments, including molecular and protostellar clouds, and the frozen surfaces of moons, planets, and comets. When exposed to ionizing and/or thermal processing, they become a nursery for new molecular species and are also responsible for their desorption to the gas-phase. Crystalline water ice, produced by the deposition of gaseous water at warm (80–150 K) surfaces or by the heating of cold amorphous water ice (up to ∼150 K), is also regularly detected by astronomical observations. Here, we employed the procoda code to map the chemical evolution of 5 keV electron-bombarded crystalline water-ices at different temperatures (12, 40, 60 and 90 K). The chemical network considered a total of 61 coupled reactions involving nine different chemical species within the ice. Among the results, we observe that the average calculated effective rate constants for radiation-induced dissociation decrease as the ice´s temperature increases. The abundance of molecular species in the ice at chemical equilibrium and its desorption to gas-phase depend on both the temperature of the ice. H2O molecules are the dominant desorbed species, with a desorption yield of about 1 molecule per 100 electrons, which seems to be enhanced for warmer crystalline ices. The obtained results can be employed in astrochemical models to simulate the chemical evolution of interstellar and planetary environments. These findings have implications for astrochemistry and astrobiology, providing insight into crucial chemical processes and helping us understand the chemistry in cold regions in space.

Cometary surface dust layers built out of millimetre-scale aggregates: dependence of modelled cometary gas production on the layer transport properties

ABSTRACT The standard approach to obtaining knowledge about the properties of the surface layer of a comet from observations of gas production consists of two stages. First, various thermophysical models are used to calculate gas production for a few sets of parameters. Second, a comparison of observations and theoretical predictions is performed. This approach is complicated because the values of many model characteristics are known only approximately. Therefore, it is necessary to investigate the sensitivity of the simulated outgassing to variations in the properties of the surface layer. This problem was recently considered by us for aggregates up to tens of microns in size. For millimetre-size aggregates, a qualitative extension of the method used to model the structural characteristics of the layer is required. It is also necessary to study the role of radiative thermal conductivity, which may play an important role for such large particles. We investigated layers constructed from large aggregates and having various thicknesses and porosity and evaluated the effective sublimation of water ice at different heliocentric distances. For radiative conductivity, approximate commonly used models and the complicated model based on the dense-medium radiative transfer theory were compared. It was shown that for millimetre-size aggregates careful consideration of the radiative thermal conductivity is required since this mechanism of energy transfer may change the resulting gas productivity by several times. We demonstrate that our model is more realistic for an evolved comet than simple models parameterizing the properties of the cometary surface layer, yet maintains comparable computational complexity.

Modelling of the processes of dunes formation on the surface of comet 67P/Churyumov–Gerasimenko

ABSTRACTThe paper discusses selected mechanisms that can significantly contribute to the formation of dune structures on the surface of comets. In these considerations, mechanisms related to the dynamics of dust on the surface of the nucleus were taken into account. In the first step, the maximum dimensions of particles that can be entrained into the cometary atmosphere as a result of local sublimation of water ice have been determined. Then, three processes related to the movement of dust on the surface of the comet nucleus were considered, such as dust migration, a landslide that can turn into a local avalanche, and fluidization of unstable dust. The thermodynamic process responsible for initiating these mechanisms is ice sublimation. These mechanisms cause the dust to move by liquefaction of the given layer. Due to the complex and regular structure of a cometary nucleus, moving dust can suddenly stop flowing. For individual mechanisms responsible for the formation of dune structures, the maximum angular width for a given particle, particle mobility coefficient, and static fluidization pressure were determined. It should be noted that in those parts of the comet nucleus that are illuminated by the Sun, these mechanisms can occur almost simultaneously, leading to the formation of dune structures by piling up particles.

Mapping the evolution of pure CO2 ices irradiated by ions, UV, and electrons using the upgraded PROCODA code (employing an effective rate constant ordering by thermochemistry data)

Chemical reactions and desorption processes are being triggered by incoming ionizing radiation over astrophysical ices in cold space environments. The quantification of these processes is crucial to achieve a detailed understanding of the underlying chemistry occurring within the ice. With this goal, we have upgraded the PROCODA code (Pilling et al. 2022a) which solves a system of coupled differential equations and describes the evolution of the molecular abundances under processing by radiation, now including an effective rate constant (ERCs) ordering by employing thermochemistry data taken from literature. This methodology helps to identify the most important reactions within the reaction network and therefore decreases the degeneracy of the solutions and enhancing the accuracy of the calculations. Here, we described the chemical evolution of four irradiated pure CO2 considering 11 different chemical species, 100 reaction routes and 11 radiation-induced desorption processes. The best-fit models provide the effective rate constants, several desorption parameters, as well as, the characterization of the chemical equilibrium (CE) phase. A comparison with previous code version was given and indicates that the ordering of rate constants by thermochemistry data is more important when more energy is deposited in the ice. The current work present more realistic values for the effective rate constants and a better characterization of the CE phase, such data can be used to refine astrochemical models to better describe cold space environments in the presence of incoming ionizing radiation field such molecular clouds and protoplanetary regions and the surface of comets and frozen moons and planets.