지난 주에는 SDSS 자료로 엑셀로 HR도를 그려봤는데 오늘은 SIMBAD라는 곳에서 자료를 내려받아 TOPCAT라는 프로그램으로 HR도를 그려봤다.
과정은 매우 길지만 잊지 않도록 중요한 부분만 정리해 둬야겠어.
- SIMBAD basicquery에서 천체를 검색한다
2. children을 눌러 데이터를 추출하는
3. 그럼 이렇게 탁산의 데이터가 나오는데 1147개의 천체에 대한 자료가 나왔다. M45와 M67은 9000개 이상의 데이터가 추출되어 상당한 시간이 걸린다.그런데 prob.member와 별이 이 성단에 속해 있을 가능성도 함께 표출하지만 0100%까지 모두 추출되므로 코딩을 통해 확률이 낮은 별은 제거해야 하는데 지금은 대략적인 그림을 그리기 위해 패스. 맨 아래에 내려가서 Votable 형태로 데이터를 다운로드 한다.
4. TOPCAT에서 Open a new ta ble – System Browser로 다운로드 받은 파일을 열면 프로그램에 포함된다.
5. 위 메뉴에서 Display co lumn metadata에 들어가면 성단에 관련된 모든 데이터가 다 포함되어 있는데, 필요한 FLUX_B와 FLUX_V만 남기고 모두 체크를 지운다.
6. 지금 B-V 값을 계산해야 하는데 왼쪽 상단의 플러스 버튼을 눌러 그렇게 계산한 값을 테이블에 추가해 준다.
7.x, y축을 설정하고 그래프를 예쁘게 장식하면
이렇게 그래프를 그릴 수 있지만, 그 전주 엑셀에서 그린 그림보다는 좀 더 전문적인 느낌(?)이다.(그림 그리는 과정은 힘들지만 예쁘게 꾸미는 것은 언제든지 재미있다.)
같은 방법으로 M45와 M67의 HR도를 x, y축의 범위를 맞춰서 그려보면
이렇게 나와.M15는 거의 모든 별이 MS를 떠나 있는 늙은 별이며 M45는 V등급이 46등급인 별이 마침 MS를 떠나버리는 젊은 별, M67는 2214등급으로 전향점을 갖는 중간 별임을 확인할 수 있다. 실제 관측한 사진을 보면 M16은 붉은 빛을, M45는 푸른 빛을 내기 때문에 두 성단의 나이가 확연히 차이가 나는 것을 확인할 수 있다.
M15 위의 가지에서 두 갈래로 나뉘는 것처럼 보이는 것은 1)배경성단 때문인지, 2)구상성단의 multi-population이 관측된 것일 수 있으나 이 데이터가 얼마나 deep로 관측되었는지 모르기 때문에 쉽게 판단할 수 없다.