지난주에는 SDSS 자료로 엑셀로 HR도를 그려봤는데 오늘은 SIMBAD라는 곳에서 자료를 다운받아 TOPCAT라는 프로그램으로 HR도를 그려봤다.
과정은 너무 길지만 잊지 않도록 중요한 부분만 정리해둬야죠.
- SIMBA Dbasicquery에서 천체검색하기
2. children을 눌러 데이터를 추출한다
3. 그러자 이렇게 많은 데이터가 나오는데 1147개의 천체에 대한 자료가 나왔다. M45와 M67은 9000개 이상의 데이터가 추출돼 있어 시간이 걸린다.그런데 prob.member라는 별이 이 성단에 속해 있을 가능성도 동시에 표출되는데 0~100%까지 모두 추출되기 때문에 코딩에 따라 확률이 낮은 별은 없애야 하는데 지금은 대략적인 그림을 그리기 때문에 패스. 맨 아래로 내려가 Votable 형식으로 데이터 다운로드를 한다.
4.TOPCAT에서 Openanewtable-System Browser에서 다운로드 받은 파일을 열면 프로그램에 포함된다.
5. 위 메뉴에서 Displaycolumnmetadata에 들어가면 성단에 대한 모든 데이터가 포함되어 있는데 필요한 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에 관측했는지 알 수 없기 때문에 쉽게 판단할 수 없다.