Fitting infrared ice spectra with genetic modelling algorithms. Presenting the ENIIGMA fitting tool

Abstract

Context. A variety of laboratory ice spectra simulating different chemical environments, ice morphology as well as thermal and energetic processing are demanded to provide an accurate interpretation of the infrared spectra of protostars. To answer which combination of laboratory data best fit the observations, an automated statistically-based computational approach becomes necessary. Aims. To introduce a new approach, based on evolutionary algorithms, to search for molecules in ice mantles via spectral decomposition of infrared observational data with laboratory ice spectra. Methods. A publicly available and open-source fitting tool, called ENIIGMA (dEcompositioN of Infrared Ice features using Genetic Modelling Algorithms), is introduced. The tool has dedicated Python functions to carry out continuum determination of the protostellar spectra, silicate extraction, spectral decomposition and statistical analysis to calculate confidence intervals and quantify degeneracy. As an assessment of the code, several tests were conducted with known ice samples and constructed mixtures. A complete analysis of the Elias 29 spectrum was performed as well. Results. The ENIIGMA fitting tool can identify the correct ice samples and their fractions in all checks with known samples tested in this paper. Concerning the Elias 29 spectrum, the broad spectral range between 2.5-20 $\mu$m was successfully decomposed after continuum determination and silicate extraction. This analysis allowed the identification of different molecules in the ice mantle, including a tentative detection of CH${}_3$CH${}_2$OH. Conclusions. The ENIIGMA is a toolbox for spectroscopy analysis of infrared spectra that is well-timed with the launch of the James Webb Space Telescope. Additionally, it allows for exploring the different chemical environments and irradiation fields in order to correctly interpret astronomical observations.

23 pages, 19 figures, 3 tables. Accepted for publication in A&A

References in corpus (15)

• Array Programming with NumPy
• The evolution of amorphous hydrocarbons in the ISM: dust modelling from a new vantage point
• The c2d Spitzer spectroscopic survey of ices around low-mass young stellar objects II: CO2
• Atom addition reactions in interstellar ice analogues
• Desorption of CO and O2 interstellar ice analogs
• A non-energetic mechanism for glycine formation in the interstellar medium
• ALMA observations of the kinematics and chemistry of disc formation
• Physical and chemical fingerprint of protostellar disc formation
• Infrared spectroscopy of HCOOH in interstellar ice analogues
• The Composition of Interstellar Grains Toward Zeta Ophiuchi: Constraining the Elemental Budget Near the Diffuse-Dense Cloud Transition
• Properties of protostars in the Elephant Trunk globule IC 1396A
• Fitting infrared ice spectra with genetic modelling algorithms. Presenting the ENIIGMA fitting tool
• Two dimensional ice mapping of molecular cores
• Computational modeling of the class I low-mass protostar Elias 29 applying optical constants of ices processed by high energy cosmic ray analogs
• The use of genetic algorithm to model protoplanetary discs

Cited by in corpus (1)

• Fitting infrared ice spectra with genetic modelling algorithms. Presenting the ENIIGMA fitting tool