내부 구조가 거대 외계행성 금속성에 미치는 영향

초록

본 연구는 44개의 가스 거대 외계행성을 대상으로 세 가지 내부 구조 가설(핵+대류층, 희석 핵, 완전 혼합)을 적용해 진화 모델을 계산하고, 각 구조가 추정되는 전체 금속 함량(Z) 및 질량‑금속성 관계에 미치는 영향을 조사하였다. 모든 구조에서 무거운 원소 질량(MZ)은 행성 질량(M)과 양의 상관관계를 보였으며, 전체 금속성(Z)은 질량이 증가할수록 감소하는 경향을 나타냈다. 희석 핵 모델은 핵‑대류층 모델과 거의 동일한 금속성을 제공했으며, 대기 금속성 증가와 비정상적 대류(비단열) 모델은 추정 금속성을 최대 35 %까지 상승시켰다. 저질량(0.2 MJ 이하) 행성들이 금속성‑질량 반비례 관계를 주도하고, 1 MJ 이상 행성에서도 예상보다 높은 금속성(Z≈0.1–0.3)이 관측되었다. 향후 PLATO와 Ariel 등 미션을 통한 정밀 질량·반지름·대기 조성 측정이 내부 구조와 형성 모델을 제약하는 데 핵심이 될 것이다.

상세 분석

이 논문은 가스 거대 외계행성의 내부 구조가 행성 전체 금속 함량(Z) 추정에 미치는 영향을 정량적으로 평가한다. 연구팀은 44개의 행성을 선정했으며, 질량 0.12–5.98 MJ, 반지름·연령·평형 온도 등 관측값의 불확실성을 5 % 이하(반지름)·10 % 이하(질량)로 제한하였다. 내부 구조 가설은 크게 세 가지로 나뉜다. 첫 번째는 전통적인 핵‑대류층(CE) 모델로, 무거운 원소는 별도의 핵에 집중되고 대류층은 순수 H‑He라고 가정한다. 두 번째는 완전 혼합(FM) 모델로, 무거운 원소가 행성 전체에 균일하게 섞여 있다. 세 번째는 희석 핵(DC) 모델로, 핵‑대류층 경계가 뚜렷하지 않고 무거운 원소가 질량 좌표에 따라 점진적으로 감소하는 구배를 가진다. 또한 대기 금속성(Z_atm)을 별 금속성(Z★) 혹은 3 × Z★로 설정한 변형 모델(DCA, DCA3)도 고려하였다.

모델링에는 CEPAM 코드를 사용해 1‑D 구조·열진화 방정식을 풀었으며, 물·암석은 SESAME, H‑He는 최신 CMS19+HG23 방정식을 적용했다. 대기에서는 비회색 복사 전이와 두 밴드(가시·IR) 평균 불투명도를 사용했으며, 무거운 원소가 증가하면 불투명도가 상승해 냉각 속도가 늦어진다. 모든 구조에 대해 기본적으로는 단열(P‑T) 프로파일을 가정했지만, 비단열(비정상 대류) 경우도 일부 실험했다.

금속성 추정은 역전파 방식으로 수행되었다. 관측된 질량·반지름·연령을 입력으로, 금속성 Z를 변화시켜 모델 반지름이 관측값과 일치하도록 루트‑파인딩(Brent)으로 해결한다. 각 행성에 대해 5 000개의 무작위 샘플을 생성해 Z의 사후 분포를 얻고, 평균·표준편차를 산출했다.

결과적으로, FM 모델은 대기·대류층에 무거운 원소가 섞여 평균 분자량과 불투명도가 증가하므로, 동일한 반지름을 맞추기 위해 더 큰 MZ가 필요했다. 반면 CE와 DC 모델은 비슷한 Z 값을 보였으며, DC 모델은 핵‑대류층 경계가 흐려져 무거운 원소가 외부로 퍼져 밀도가 상승하지만, 전체 금속성은 CE와 거의 동일하거나 약간 낮았다. 대기 금속성을 Z★ 혹은 3 × Z★로 높이면 불투명도 상승 효과가 커져 추정 Z가 최대 35 %까지 증가했다. 특히 저질량 행성(0.2 MJ 이하)은 높은 금속성을 보이며 전체 Z‑M 반비례 관계를 주도했으며, 1 MJ 이상 행성에서도 Z≈0.1–0.3이라는 예상보다 높은 금속성이 관측되었다.

이러한 결과는 행성 형성 이론에 중요한 제약을 제공한다. 코어‑축적 모델은 저질량 행성에서 높은 금속성을 기대하지만, 고질량 행성에서도 높은 금속성이 존재한다는 점은 가스 흡수 단계에서의 고체 물질 함유량이 크거나, 행성 내부 혼합이 활발히 일어났을 가능성을 시사한다. 향후 PLATO와 Ariel이 제공할 정밀 질량·반지름·대기 조성 데이터는 비단열 효과와 대기 불투명도 파라미터를 더 정확히 측정하게 해, 희석 핵 모델의 존재 여부와 그 폭을 구체화할 수 있을 것이다.