RT info:eu-repo/semantics/masterThesis
T1 Relaciones de temperatura para el ion de Cl2+ en regiones HII y la determinación de abundancias químicas
A1 Orte García, Maialen
AB H II regions and planetary nebulae are created when hot stars ionise nearby gas clouds.H II regions are ionised with high-energy UV radiation from young and massive stars, andthey are primarily composed of hydrogen, but we can also find elements such as helium,carbon, nitrogen and oxygen, to name some of the most common ones.Although they may produce thousands of stars, generally less than the 10 % of theavailable gas is converted into stars, while the rest of it is dispersed by phenomena likestellar winds, radiation pressure and supernova explosions.Elements heavier than H and He are formed from stellar nucleosynthesis processes, andare sent back to interstellar medium through different mechanisms such as gas ejections,all of them based on the same idea: the outgassing of stars during the last stages of theirlives.These stars in their initial and final phases emit photons that can ionize the surrounding gas, increasing the amount of free electrons that are in the ionized nebulae. Theseelectrons interact creating what we call thermal equilibrium and photoionization equilibrium.Considering these processes, we can determine the chemical composition of the gas,which depends on three parameters: the electron density (ne), temperature (Te) and thechemical abundance of the ion that emits the radiation. So, if we can determine the firsttwo parameters, ne and Te, we can obtain the ionic abundances, once we have the intensityof the appropiate emission lines.The main objective of this study is Chlorine (Cl), which belongs to group VIIA of theperiodic table, known as halogens. Spectral lines of Cl are quite unusual in stellar spectra,so the determinations of its abundance in stars are very few. Because of this, most ofthe information that is available nowadays about Cl abundance comes from analyzingemission lines in ionized nebulae.The calculation of Cl abundance is based on measuring the intensity of the [Cl III]line doublet, at 5517 and 5537 ˚A, which are the brightest lines of Cl in the optical range.Moreover, in common conditions of an H II region, we expect that only Cl+ and Cl3+ havea significant contribution to total Cl abundance, as well as Cl2+, which is the dominantone.Nevertheless, the most common situation is that where we can only get [Cl III] linesto obtain total Cl abundance, and this is when we use the ionization correction factors(ICF), which correct the contribution to the total abundance of an element given by thoseions we can not observe.But there is still one more problem to get through with Cl: the correct determinationof the electron temperature representative of the nebular volume where Cl2+ is present.When it cannot be obtained through [Cl III] lines, we have to use valid diagnoses for otherions, such as O2+ and N+.First of all, we consider photoionization models, which through a numerical code, solve all the physical processes inside ionized nebulae, and calculate the physical conditions(Te, ne) and emitted spectra. In this case, we use the model grid from Vale Asari et al.(2016), available at the Mexican Million Models Database (3MdB) (Morisset, 2015). Moreover, we also use the constraints proposed by Amayo (2019) to this same grid, based onobservational arguments, which, in the end, make the model grid a realistic one.One of the aims of this study, is to compare the results with real observations, in orderto conclude if the temperature we choose is valid or not for representing Te(Cl2+). For it, we have an observational sample of 65 spectra of 33 extragalactic and 15 Galactic H IIregions, obtained from different publications of our research group. For these spectra, wetake the physical conditions calculated in M´endez-Delgado et al. (2022a).Through the photoionization models, we explore the behavior of Te as a function of theionization degree, and, as we expected, we see that Te(N+) better represents zones witha lower ionization degree, and vice versa for Te(O2+). So we think about a temperaturethat can represent the most common conditions of ionized nebulae in the volume whereCl2+ is present, as a combination of Te(O2+) and Te(N+). We define the representativetemperature, which depends on P instead of the ionization degree, since P is directlyobtained through observational data, given by the ratio of line intensities.We represent Te(Cl2+) as a function of these three temperatures, Te(N+), Te(O2+) andTrep, and we see that Trep is the one that better fits the results. On the other hand, byanalyzing the exact values of the standard deviation that results when we use Te(N+) orTe(O2+), we see that it is lower for Te(N+) for P ≤ 0.3, and lower for Te(O2+) for P > 0.3.This leads us to define a new temperature, named reference temperature.After representing Te(Cl2+) as a function of all the temperatures, we see that both,Trep and TP , are more appropiate for representing Te(Cl2+), since they both fit betterthan Te(N+) or Te(O2+).Once we establish the appropiate temperatures, we process the data of the 65 spectraand start determining the ionic abundances. For it, we use PyNeb (Luridiana et al., 2015),a numerical code made for calculus in ionized nebulae, which gives us the ionic abundanceof an ion through the getIonAbundance() function.In order to obtain associated uncertainties, we also calculate these abundances througha Monte Carlo simulation of 500 experiments, so that we can use this value for the restof the study.In our data, all the 65 spectra have intensities for [Cl III] emission lines, so we calculatethe ionic abundance of Cl2+ for all of them. Moreover, there are 14 spectra in whichthe intensities for both [Cl II] and [Cl IV] emission lines are also given, so, for these, wecan obtain the ionic abundances for Cl+ and Cl3+. This allows us to calculate the totalabundance of Cl by just adding the ionic abundances of Cl+, Cl2+ and Cl3+, in additionto the method using the ICF.Moreover, we will also calculate the Cl/O ratio, which is particularly interesting, dueto some studies reveal that its value must be constant. The reason for this is that Cl isproduced in the same zones where the O burning happens, so both elements have a similarnucleosynthetic origin (Esteban et al., 2015). For a reference, we use the value proposedby Lodders (2019) for the Sun, log(Cl/O)⊙ = −3.50 ± 0.09.We start calculating the total abundance of Cl as the sum of the ionic abundances,for each proposed temperature, and the 14 spectra for which we have calculated the threeneeded abundances. In fact, we add one more spectrum to this case: the Herbig-Haroobject HH204, from the Orion nebula, due to its ionization degree (parameterized byO2+/O) being nearly zero. For it, we expect Cl/H=Cl+/H++Cl2+/H+.Once we obtain the total abundance of Cl, we calculate the Cl/O ratio for these15 spectra, for each proposed temperature. For analyzing all these results, we representthem graphically, as follows: 12 + log(Cl/H) as a function of 12 + log(O/H), log(Cl/O) asa function of O2+/O, and log(Cl/O) as a function of 12 + log(O/H). Although the resultswe obtain for this method seem to be correct, we cannot draw a strong conclusion, dueto the low number of points available for this method.Next, we calculate the total abundance through the ICF. In this case, we use the ICF proposed by Amayo et al. (2021), who also gives mathematical expressions to estimatethe uncertainties associated to this factor. Once we get Cl/O, we continue by calculatingthe total abundance of Cl. Then, we repeat the three representations described before,and analyze our results. We have to take into account that, with this method, the numberof points is considerably higher, since we have used the 65 spectra available, excludingM42-HH204.In the results given in the ICF method, we see that the bias of the total abundance ofCl as a function of O is quite similar to the one obtained with the sum of ionic abundances,for all the 4 proposed temperatures. On the other hand, for the Cl/O ratio, we see that,even though the weighted average of the represented points is nearly -3.5, the ratio doesnot remain constant. We must also consider that the uncertainties given for this methodare considerably higher, especially for those points with a lower ionization degree.Finally, we end our analysis by representing the total abundance of Cl obtained throughthe ICF, as a function of the one calculated as the sum of ionic abundances. This allows usto see how the ICF underestimates the total abundance of Cl, since nearly all the points,for the 4 used temperatures, are located under the line representing a 1:1 relationship.As we found at the beginning of this study, the temperatures that are more appropiaterepresenting Te(Cl2+) are Trep and TP . However, the uncertainties given for TP are a bithigher than those obtained for Trep, especially for spectra with a low ionization degree.Trep is a linear combination of Te(N+) and Te(O2+), so the uncertainties are calculated using the error propagation, where the uncertainty associated to P takes part and, therefore, makes a difference.
YR 2023
FD 2023
LK http://riull.ull.es/xmlui/handle/915/32443
UL http://riull.ull.es/xmlui/handle/915/32443
LA es
NO Máster en Astrofísica
DS Repositorio institucional de la Universidad de La Laguna
RD 01-jun-2024