Determination of the mass loss rate and the terminal velocity of stellar winds. I. Genetic algorithm for automatic line profile fitting



Document title: Determination of the mass loss rate and the terminal velocity of stellar winds. I. Genetic algorithm for automatic line profile fitting
Journal: Revista mexicana de astronomía y astrofísica
Database: PERIÓDICA
System number: 000377901
ISSN: 0185-1101
Authors: 1
1
Institutions: 1Universidad Nacional Autónoma de México, Instituto de Astronomía, México, Distrito Federal. México
Year:
Season: Abr
Volumen: 41
Number: 1
Pages: 121-129
Country: México
Language: Inglés
Document type: Artículo
Approach: Experimental
Spanish abstract Se presenta un nuevo método de ajuste automático de perfiles P Cyg observados en espectros UV de vientos estelares. La función fuente de las líneas se calcula usando la aproximación de Sobolev y el flujo emitido se obtiene a través de la solución exacta de la ecuación de transporte (similar al método SEI descrito por Lamer et al. 1987). La calidad del ajuste se evalúa usando el estimador de verosimilitud. La maximización de la verosimilitud se hace a través de un algoritmo genético. Las ventajas de nuestro método comparado con otros métodos similares son su estabilidad y su poca sensibilidad de las condiciones iniciales. Además, el método garantiza la localización del máximo global de la superficie de verosimilitud, lo que no se puede hacer con los algoritmos clásicos. Presentamos la implementación del método, las pruebas de su funcionalidad usando datos sintéticos y reales y la estimación de los intervalos de confianza de los resultados
English abstract A new method for automatic fitting of P Cyg line profiles in UV spectra of stellar winds is presented. The line source function is calculated using Sobolev's approximation and the emergent flux is obtained by exact integration of the equation of the radiation transport (similar to the SEI method described by Lamers et al. (1987)). The quality of the fit is evaluated using the likelihood estimator. The maximization of the likelihood is done by a genetic algorithm. The advantages of our method with respect to other similar approaches are its robustness and its insensibility to the initial guess. In addition, the algorithm guarantees the localization of the global maximum of the likelihood hypersurface, which is not the case for classical minimization algorithms. Here we present an implementation of the genetic algorithm for line profile fitting, its tests on both synthetic and real data and an estimation of the confidence limits of the results
Disciplines: Física y astronomía
Keyword: Astronomía,
Estrellas,
Pérdida de masa,
Viento estelar,
Velocidad
Keyword: Physics and astronomy,
Astronomy,
Stars,
Mass loss,
Stellar wind,
Velocity
Full text: Texto completo (Ver HTML)