Drifters on the edge of town: $λ$ Boötis stars in clusters

기술 요약: 도시 외곽의 표류자들: 클러스터 내의 λ Boötis 별들

문제 제述
λ Boötis (λ Boo) 별은 가벼운 원소(C, N, O, S)는 태양과 유사한 함량을 보이지만 철족 원소는 현저히 결핍된 화학적 특이성을 가진 A형 별이다. 이들에 대한 주요 이론적 설명은 이 별들이 성간 매질(ISM)로부터 금속 함량이 낮은 순수한 가스를 강착(accretion)한다는 가설을 세운다. 그러나 이 시나리오는 중요한 문제에 직면해 있다: 최근 성단 내부 또는 외곽에서 발견된 λ Boo 별들은, 그러한 환경 내 중간 질량 및 대질량 별(> 5 M⊙)들의 강력한 이온화 복사가 순수한 ISM 가스의 강착을 방해하거나 강착에 필요한 주변 성간 원반(circumstellar disc)을 제거할 것이라는 예상과 모순된다. 또한, λ Boo의 함량을 생성하는 데 필요한 강착률(10⁻¹⁴ ~ 10⁻⁹ M⊙ yr⁻¹)과 백색 왜다 관측을 통해 추론된 훨씬 낮은 강착률(< 10⁻¹⁷ M⊙ yr⁻¹) 사이에는 긴장 관계가 존재한다. 본 연구가 다루는 구체적인 문제는, 진화하는 클러스터 내 A형 별들의 역학적 이력이 클러스터의 고복사 환경으로부터 탈출하여 순수한 가스를 강착한 후 다시 돌아올 수 있는지 여부이며, 이를 통해 클러스터 환경 내 λ Boo 별의 존재와 ISM 강착 가설을 화해시킬 수 있는지이다.

연구 방법론
저자들은 100 Myr 동안 성단 형성 영역의 형성과 역학적 진화를 모델링하기 위해 순수 중력 N-body 시뮬레이션을 사용한다.

• 시뮬레이션 설정: 총 항성 수($N_*$)가 150, 500, 1000, 3000인 네 가지 세트의 시뮬레이션을 수행하였다. 초기 조건은 박스 프랙탈 방법(프랙탈 차원 $D=1.6$)을 통해 생성된 공간적 및 운동학적 하부 구조와 $\alpha_{vir} = 0.3$의 비리알 비율를 포함하였으며, 이는 초기 붕괴와 중심 집중형 클러스터 형성을 유도하였다. 항성 질량은 상한선 50 M⊙를 가진 Maschberger (2013) 초기 질량 함수(IMF)를 따랐다. 원시 쌍성은 제외되었다.
• 후처리 분석: 시뮬레이션은 다음 사항을 결정하기 위해 사후 분석되었다:
  1. 강착 가능성: A형 별들을 추적하여 이들이 야코비(Jacobi, 조석) 반지름($r_J$) 너머로 이동하는 시점을 식별하였다. $r_J$ 외부로 이동했을 때, 정지된 ISM 구름($\rho = 10$ cm⁻³)에 대한 항성의 속도를 기반으로 Bondi-Hoyle-Lyttleton (BHL) 강착률을 계산하였다.
  2. 복사 환경: 각 A-별에 입사되는 원자외선(FUV) 및 극자외선(EUV) 플럭스를 대질량 별(> 5 M⊙ 및 > 20 M⊙)과의 거리에 기초하여 계산하였다.
  3. 원반 생존: FRIED 그리드의 광증발(photoevaporation) 모델과 해석적 처방을 사용하여 A-별의 주변 성간 원반 질량 손실을 모델링하였다. 총 질량 손실이 초기 원반 질량(항성 질량의 10%로 가정)을 초과할 경우 원반이 "파괴"된 것으로 간주하였다. λ Boo 형성 후보가 되기 위해서는 원반을 유지하면서 ISM 물질을 강착해야 했다.
  4. 시간적 제약: Type Ia 초신위에 의한 ISM의 농축 가능성을 고려하여, 40 Myr 이후의 강착을 제한하는 제약을 두어 계산을 반복하였다.

주요 결과

• 역학적 이력: 모든 시뮬레이션 실현 과정에서, A-별의 대다수(90–98%)는 100 Myr의 진화 과정 중 어느 시점에 야코비 반지름 너머로 이동하였다. 이는 클러스터 핵으로부터의 방출(ejection) 또는 별들이 조석 반지름 밖으로 이동했다가 다시 내부 영역으로 돌아오는 "표류(drifting)"를 통해 발생하였다.
• 강착률: 야코비 반지름 너머로 이동한 A-별들은 $10^{-12}$에서 $10^{-5}$ M⊙ yr⁻¹ 범위의 BHL 강착률을 달성하였다. 높은 속도로 방출된 별들은 낮은 강착률($10^{-12}–10^{-10}$ M⊙ yr⁻¹)을 보였으나, 클러스터 중심 궤도를 따라 밖으로 표류했다가 다시 돌아온 별들은 더 높은 강착률($10^{-9}–10^{-5}$ M⊙ yr⁻¹)을 보였으며, 이는 λ Boo 함량을 설명하는 데 필요한 범위와 겹쳤다.
• 원반 생존 vs 복사: 주변 성간 원반의 생존은 클러스터의 총 항성 수(따라서 이온화원의 수)에 크게 의존하였다.
  • 저질량 영역($N_* = 150$)에서는 야코비 반지름 너머로 나간 A-별의 약 55%가 10 Myr 이상 원반을 유지하였다.
  • 고질량 영역($N_* = 3000$)에서는 강렬한 FUV/EUV 장으로 인해 A-별의 원반 유지율이 0%였으며, 이는 해당 환경에서 λ Boo 형성을 사실상 불가능하게 만든다.
  • 중간 규모 영역($N_* = 500, 1000$)에서는 일부 A-별(각각 8.8%와 4.8%)이 강착하는 동안 원반을 유지하였다.
• 쌍성계: 연구는 A-별을 포함하는 넓은 쌍성계(wide binary systems, >100 au)가 성단 형성 영역의 해체 과정 중에 포획을 통해 형성될 수 있음을 발견하였다. 결정적으로, 별은 순수한 ISM 물질을 강착하여 λ Boo 함량을 발달시킨 후 동반성을 포획하거나, 그 반대의 과정을 거칠 수 있다. 이러한 메커니즘은 구성 성분 간에 현저히 다른 화학적 함량을 가진 쌍성계를 형성할 수 있게 한다. 시뮬레이션된 이러한 쌍성의 분리 분포는 λ Boo 별을 포함하는 관측된 넓은 쌍성계(예: HD 87304 + CD-33 6615B)와 교차한다.

의의 및 주장
본 논문은 클러스터 내부 또는 근처에 λ Boo 별이 존재하는 것이 ISM 강착 시나리오를 본질적으로 부정하는 것은 아님을 입증한다고 주장한다. 저자들은 다음과 같이 논증한다:

1. 역학적 탈출은 충분하다: A-별이 순수한 가스를 강착하기 위해 영구적으로 방출될 필요는 없다. 이들은 조석 반지름 너머로 표류했다가 다시 클러스터 중심으로 돌아올 수 있다.
2. 복사 임계치: 강착 시나리오의 생존 가능성은 클러스터 질량과 밀접하게 연관되어 있다. 저질량 클러스터(최근 λ Boo 별이 발견된 Theia 139 및 HSC 1640과 유사한 환경)는 원반 생존을 허용할 만큼 약한 복사장을 갖는 반면, 고질량 클러스터는 그렇지 않다.
3. 쌍성 진단: 구성 성분 간의 상이한 함량을 보이는 쌍성계의 관측은 강착 가설을 무효화하지 않는다. 이러한 시스템은 강착 사건이 일어난 후에 역학적으로 형성될 수 있다.

저자들은 절대적인 강착률의 타당성에 대해서는 신중한 태도를 유지하며, 요구되는 강착률과 백색 왜다 한계 사이의 긴장이 여전히 해결되지 않았음을 언급하였다. 그들은 자신들의 역학 모델이 클러스터 내 λ Boo 별 형성의 가능성을 뒷받침하지만, 구체적인 강착률과 강착되는 원반 및 복사장의 물리적 상호작용을 검증하기 위해서는 추가적인 연구가 필요하다고 결론지었다.