[논문 리뷰] Determination of the mass loss rate and the terminal velocity of stellar winds. I Genetic algorithm for automatic line profile fitting

Л. Георгиев, X. Hernández|ArXiv.org|2005. 01. 28.

Calibration and Measurement Techniques인용 수 3

한 줄 요약

이 논문은 우주선 풍의 UV 별풍 스펙트럼에서 P Cyg 선형형을 자동으로 피팅하기 위한 유전 알고리즘을 제시한다. Sobolev 근사와 정확한 복사 전달을 사용하여 질량 손실률과 최종 속도를 추정한다. 이 방법은 전반적인 가능도 최대값을 안정적으로 탐색하여 신호 대 잡음비가 낮고 해상도가 낮은 데이터에서도 신뢰할 수 있는 매개변수 추정이 가능하며, 몬테카를로 시뮬레이션을 통해 통계적으로 의미 있는 신뢰구간을 제공한다.

ABSTRACT

Terminal wind velocity and mass loss rate are the most fundamental parameters of stellar winds. Unfortunately, their determination calls for high resolution spectroscopy in a range of wavelengths spanning from the ultraviolet to the infrared. For weak and/or distant objects, this becomes unfeasible. Nevertheless, it is possible to obtain an accurate estimate of these parameters through a simplified study of the formation processes of resonant lines which show P Cyg profiles. In this case, the line profile is a complex function of 6 parameters, with the radiative transport treatment becoming relatively inexpensive. However, preforming a 6 dimensional parameter fit raises a number of problems if one seeks an objective and automatic procedure to yield the optimal values, from which wind velocities and mass loss rates can be estimated. Useing a Likelihood function to construct a well defined statistical estimator of the goodness of fit which corresponds to a given model, we turn to a genetic algorithm through which we find the global maximum of the 6 dimensional Likelihood hyper-surface. We here present the implementation of the method, its successful testing with synthetic line profiles, where the answer is known in advance, together with first results of its application to real data

연구 동기 및 목표

• UV 별풍 스펙트럼에서 P Cyg 선형형을 신뢰할 수 있고 자동으로 피팅할 수 있는 방법을 개발하는 것.
• 고차원 매개변수 공간에서 국소 최소값에 갇히기 쉬운 전통적 최소화 방법의 한계를 극복하는 것.
• 낮은 신호 대 잡음비와 낮은 해상도의 UV 스펙트럼으로부터 질량 손실률과 최종 속도를 정확하게 추정할 수 있도록 하는 것.
• 몬테카를로 시뮬레이션을 통해 피팅된 매개변수에 대한 객관적이고 통계적으로 의미 있는 신뢰구간을 제공하는 것.
• 초기 매개변수 추측에 대한 의존도를 줄여 다양한 관측 조건에서의 안정성을 확보하는 것.

제안 방법

• 이 방법은 구형 대칭 별풍에서 선의 소스 함수를 계산하기 위해 Sobolev 근사를 사용한다.
• 복사 전달은 통합을 통해 정확히 해결하며, SEI 방법(Lamers 등 1987)과 유사하게 비출 유량 프로파일을 계산한다.
• 적합도 평가를 위해 가능도 추정을 사용하며, 국소 최소값에 갇히지 않도록 유전 알고리즘을 사용해 최대화를 수행한다.
• 유전 알고리즘은 매개변수 격자 위에서 작동하며, 초기 범위가 넓기만 하면 되며, 초기 추측에 대한 민감도가 매우 낮다.
• 회복된 매개변수에 대해 몬테카를로 시뮬레이션을 적용하여 전체 신뢰구간과 공분산 행렬을 추정한다.
• 미래 연구에서는 유전 알고리즘으로 전역 최대값을 탐색한 후, 전통적 최소화 방법을 사용해 해를 정밀화하는 하이브리드 접근 방식을 적용할 수 있다.

실험 결과

연구 질문

• RQ1 유전 알고리즘이 고차원 매개변수 공간에서 국소 최소값에 갇히지 않고 UV 별풍 스펙트럼의 P Cyg 선형형을 안정적으로 피팅할 수 있는가?
• RQ2 낮은 신호 대 잡음비와 낮은 해상도의 UV 데이터로부터 질량 손실률과 최종 속도를 얼마나 정확하게 추정할 수 있는가?
• RQ3 이 방법의 성능이 초기 매개변수 추측에 얼마나 의존하는가?
• RQ4 피팅 과정에서 객관적으로 통계적으로 의미 있는 신뢰구간을 도출할 수 있는가?
• RQ5 이 방법은 간섭성 흡수선 또는 혼합 흡수선으로 인한 스펙트럼 오염을 얼마나 잘 처리하는가?

주요 결과

• 유전 알고리즘이 가능도 표면의 전역 최대값을 성공적으로 탐색하여 전통적 최소화 방법이 자주 갇히는 국소 최소값을 피하는 데 성공했다.
• 이 방법은 합성 스펙트럼과 실제 UV 스펙트럼 모두에서 안정적인 성능을 보였으며, 낮은 신호 대 잡음비와 스펙트럼 오염이 있는 경우에도 유사하게 잘 작동했다.
• LMC의 별 HD5980에서 P V 1117/1129 Å 선을, 상당한 간섭성 흡수선이 존재하는 상황에서도 신뢰성 있게 피팅하였으며, 오염된 스펙트럼 영역에 낮은 가중치를 할당함으로써 성능을 확보했다.
• 회복된 매개변수에 대한 몬테카를로 시뮬레이션을 통해 전체 공분산 행렬과 통계적으로 의미 있는 신뢰구간을 도출하였다.
• 이 방법은 질량 손실률과 이온화/ excitation 분수의 곱($\dot{M}q_{\text{ion}}q_{\text{ex}}$)을 추정할 수는 있으나, 저해상도 및 낮은 S/N 데이터만으로는 $\dot{M}$과 이 분수를 분리할 수 없다.
• 초기 매개변수 추측에 대해 민감도가 없으며, 초기 범위가 넓더라도 알고리즘이 전역 최대값 쪽으로 이동함을 확인하였다.