The physical evolution of Young Stellar Objects (YSOs) is accompanied by an enrichment of the molecular complexity, mainly triggered by the heating and energetic processing of the astrophysical ices. In this paper, a study of how the ice column density varies across the protostellar evolution has been performed. Tabulated data of H$_2$O, CO$_2$, CH$_3$OH, HCOOH observed by ground- and space-based telescopes toward 27 early-stage YSOs were taken from the literature. The observational data shows that ice column density and spectral index ($\alpha$), used to classify the evolutionary stage, are well correlated. A 2D continuum radiative transfer simulation containing bare and grains covered by ices at different levels of cosmic-ray processing were used to calculate the Spectral Energy Distributions (SEDs) in different angle inclinations between face-on and edge-on configuration. The H$_2$O:CO$_2$ ice mixture was used to address the H$_2$O and CO$_2$ column density variation whereas the CH$_3$OH and HCOOH are a byproduct of the virgin ice after the energetic processing. The simulated spectra were used to calculate the ice column densities of YSOs in an evolutionary sequence. As a result, the models show that the ice column density variation of HCOOH with $\alpha$ can be justified by the envelope dissipation and ice energetic processing. On the other hand, the ice column densities are mostly overestimated in the cases of H$_2$O, CO$_2$ and CH$_3$OH, even though the physical and cosmic-ray processing effects are taken into account.

The survival of astrophysical ices in the star-forming regions depends on the suitability of temperature, density and radiation conditions. In this paper, the role of Interstellar Radiation Field (ISRF) on ices in Elias 29 is addressed. This object is the most luminous protostar in $\rho$ Oph E molecular cloud, and is surrounded by many Young Stellar Objects distant only a few arcminutes. In addition, other two brightness BV stars (S1 and HD 147889) enhance the external irradiation in Elias 29. This study was carried out by using the Monte Carlo radiative transfer code RADMC-3D assuming an internal and external irradiation. As result, we found that HD 147889 dominates the ISRF, instead of the closest protostars, and contributes to enhance the external irradiation in 44 times the standard value. Furthermore, remarkable effects are observed in FIR spectrum, as well as in near-IR image. Additionally, the snowline positions of volatile compounds, such as CO, O$_2$, N$_2$ and CH$_4$ are redefined to a toroidal-shaped morphology in the envelope, with low FUV flux (10$^{-7}$ ergs cm$^{-2}$ s$^{-1}$). In such scenario, the formation of complex molecules as result of hydrogenation or oxygenation of volatile species are expected to be severely affected.

In dense and cold regions of the interstellar medium (ISM), molecules may be adsorbed onto dust grains to form the ice mantles. Once formed, they can be processed by ionizing radiation coming from stellar or interstellar medium leading to formation of several new molecules in the ice. Among the different kind of ionizing radiation, cosmic rays play an important role in the solid-phase chemistry because of the large amount of energy deposited in the ices. The physicochemical changes induced by the energetic processing of astrophysical ices are recorded in a intrinsic parameter of the matter called complex refractive index (CRI). In this paper, we present for the first time a catalogue containing 39 complex refractive indices (n, k) in the infrared from 2.0 - 16.6 micrometer for 13 different water-containing ices processed in laboratory by cosmic ray analogs. The calculation was done by using the NKABS (acronym of determination of N and K from ABSorbance data) code, which employs the Lambert-Beer and Kramers-Kronig equations to calculate the values of n and k. The results are also available at the website: http://www1.univap.br/gaa/nkabs-database/data.htm. As test case, a H2O:NH3:CO2:CH4 ice was employed in a radiative transfer simulation of a prototoplanetary disk to show that these data are indispensable to reproduce the spectrum of YSOs containing ices.

We present the study of the effects of high energy cosmic rays (CRs) over the astrophysical ices, observed toward the embedded class I protostar Elias 29, by using computational modeling and laboratory data. Its spectrum was observed with {\it Infrared Space Observatory - ISO}, covering 2.3 - 190 $\mu$m. The modeling employed the three-dimensional Monte Carlo radiative transfer code RADMC-3D (Dullemond et al. 2012) and laboratory data of bombarded ice grains by CRs analogs, and unprocessed ices (not bombarded). We are assuming that Elias 29 has a self-irradiated disk with inclination $i =$ 60$^{\circ}$, surrounded by an envelope with bipolar cavity. The results show that absorption features toward Elias 29, are better reproduced by assuming a combination between unprocessed astrophysical ices at low temperature (H$_2$O, CO, CO$_2$) and bombarded ices (H$_2$O:CO$_2$) by high energy CRs. Evidences of the ice processing around Elias 29 can be observed by the good fitting around 5.5-8.0 $\mu$m, by polar and apolar ice segregation in 15.15-15.25 $\mu$m, and by presence of the CH$_4$ and HCOOH ices. Given that non-nitrogen compounds were employed in this work, we assume that absorption around 5.5-8.0 $\mu$m should not be associated with NH$_4^+$ ion (Shutte & Khanna2003), but more probably with aliphatic ethers (e.g. R1-OCH$_2$-R2), CH$_3$CHO and related species. The results obtained in this paper are important, because they show that the environment around protostars is better modeled considering processed samples and, consequently, demonstrates the chemical evolution of the astrophysical ices.