클러스터 가장자리에서 떠도는 λ 보티스 별

초록

λ 보티스 별은 가벼운 원소는 태양과 비슷하지만 철족 원소가 결핍된 A형 별이다. 기존에는 은하간 매질(ISM)에서 금속이 적은 가스를 흡수해 이런 화학적 특성이 생긴다고 보았지만, 이 별들이 방사능이 강한 대규모 별이 존재하는 열린 클러스터 안에서도 발견되면서 그 가설에 의문이 제기되었다. 본 연구는 N‑body 시뮬레이션을 이용해 클러스터 내 A형 별들의 궤도와 방사능 환경을 추적하고, 별이 조력 반경(Jacobi radius) 밖으로 탈출했을 때 ISM을 흡수할 수 있는지, 그리고 그 후 다시 클러스터 중심부로 돌아와도 방사능이 충분히 억제되지 않는지를 검증한다. 결과는 일부 A형 별이 클러스터 외부에서 충분한 B‑H‑L(보디‑호일‑리틀턴) 흡수율을 얻고, 이후 넓은 이진계(>100 au)를 형성할 수 있음을 보여준다. 따라서 이진계에서 구성원 별의 화학적 차이가 관측되더라도 ISM 흡수 시나리오를 배제할 수 없으며, 클러스터 내·외부에서 λ 보티스 별이 존재할 수 있음을 시사한다.

상세 분석

이 논문은 λ 보티스(λ Boötis) 별의 화학적 특성이 은하간 매질(ISM)에서 금속이 거의 없는 가스를 장기간(≈1 Myr) 흡수함으로써 형성된다는 기존 가설을, 별이 밀집된 열린 클러스터 환경에서도 유지될 수 있는지를 정량적으로 검증한다. 먼저 저자는 150, 500, 1000, 3000개의 별을 포함하는 네 가지 규모의 별 형성 영역을 설정하고, 마스체르거(2013) IMF를 이용해 질량을 샘플링한다. A형 별은 1.8–2.5 M⊙ 범위로 정의하고, 초기 위치와 속도는 프랙탈 차원 D = 1.6을 갖는 박스 프랙탈 방법으로 생성해 실제 젊은 별 형성 구역에서 관측되는 공간·운동학적 서브구조를 재현한다. 각 시뮬레이션은 α_vir = 0.3(서브볼트)으로 설정해 1 Myr 이내에 중심 집중형 클러스터가 형성되도록 한다.

동역학적 진화는 4차 Hermite 적분기(kira, Starlab)를 이용해 100 Myr까지 진행되며, 별 질량 손실은 포함하지 않지만, 5 M⊙ 초과 질량을 가진 별들의 FUV·EUV 방사능을 사후 처리한다. 핵심 가정은 A형 별이 조력 반경(r_J) 바깥으로 탈출하면 주변 ISM가 ‘프리시스트(pristine)’ 상태라 가정하고, 그때 보디‑호일‑리틀턴(BHL) 공식으로 흡수율을 계산한다. 흡수 반경 r_acc = √(2.5 G M_★/v_rel²)와 ISM 밀도 ρ = 10 cm⁻³를 사용해 \dot{M}_acc = π r_acc² ρ v_rel을 구한다. 40 Myr 이후에는 Type Ia 초신성으로 인한 철 풍부 가스가 퍼지므로, 그 시점 이후의 흡수는 금속 결핍 가스를 가정하지 않는다.

디스크 손실은 FUV 플럭스를 Habing 단위(G₀)로 변환한 뒤, FRIED 모델을 적용해 \dot{M}FUV를 추정한다. 초기 디스크 질량은 별 질량의 10%이며, 반경은 r_disc = 200 au × (M★/M_⊙)^0.3 로 설정한다. 20 M⊙ 초과 별에 대해서는 EUV에 의한 추가 손실을 Johnstone(1998)식으로 보정한다. 디스크가 완전히 소실되면 이후 ISM 흡수가 불가능하다고 가정한다.

시뮬레이션 결과, 전체 A형 별 중 약 5–10%가 조력 반경을 넘어 장기간(수 Myr) 체류하면서 \dot{M}_acc이 10⁻¹⁴–10⁻⁹ M⊙ yr⁻¹ 사이에 도달한다는 것이 확인되었다. 이때 해당 별들은 평균 FUV 플럭스가 G₀ ≈ 10–100 수준에 머물러, 디스크 소실까지는 이르지 않아 충분한 물질을 흡수할 수 있다. 또한, 일부 별은 외부에서 다시 클러스터 중심부로 재진입하지만, 재진입 시점의 방사능은 여전히 G₀ < 10³ 수준으로, λ 보티스 화학 패턴을 억제하기엔 부족하다.

이와 더불어, 시뮬레이션은 외부에서 흡수한 A형 별이 다른 A형 별과의 중력 상호작용을 통해 >100 au 거리의 넓은 이진계를 형성할 확률이 ≈2%임을 보여준다. 이는 관측적으로 “동일 클러스터 내에서 화학적 차이를 보이는 이진 구성원”이 존재하더라도 ISM 흡수 시나리오를 배제할 수 없다는 중요한 함의를 제공한다.

결론적으로, λ 보티스 별이 클러스터 내부·외부에서 존재할 수 있는 동역학적 경로와, 방사능 억제 효과가 충분히 약한 상황을 정량적으로 제시함으로써, 기존의 ‘클러스터 내에서는 λ 보티스 별이 존재할 수 없다’는 주장에 반증을 제시한다. 다만, 실제 ISM 밀도와 방사능 변동, 디스크 물리(예: 자기장, 복사압) 등을 포함한 고해상도 방사선-유체 시뮬레이션이 필요하다는 점을 남긴다.