화성 규모 충돌 결과 예측 모델: 합병·스키퍼·파괴 전이의 새로운 해석

초록

본 연구는 SPH 시뮬레이션을 통해 화성 크기 천체 간 충돌을 조사하고, 겹치는 질량이 흡수하는 에너지를 기반으로 한 새로운 분석 모델을 제시한다. 이 모델은 충돌 속도·각도·질량비에 따른 ‘합병’, ‘스키퍼(hit‑and‑run)’, ‘파괴(전멸)’ 전이를 정확히 예측하며, 각도 가중 평균식을 제공해 행성 형성 시뮬레이션에 바로 적용할 수 있다.

상세 분석

본 논문은 화성 규모(≈6×10²⁶ g) 두 천체의 충돌을 96가지 파라미터 조합으로 SPH(Smooth Particle Hydrodynamics) 시뮬레이션하였다. 핵심 물리 모델은 Tillotson 상태방정식을 사용해 철 핵(30 wt %)과 화강암 맨틀(70 wt %)을 구성했으며, 인공 점성(α=1.5, β=3.0)과 중력 트리 알고리즘을 통해 충돌 후 잔해의 질량과 운동을 정확히 추적했다. 충돌 파라미터는 질량비(M_tar/M_pro = 1, 3, 10, 30), 충돌 각도(θ = 30°, 45°, 60°), 그리고 충돌 속도 비(v_imp/v_esc = 1.1~16)이다.

시뮬레이션 결과는 충돌 에너지 E_imp와 두 체계 중력 에너지 E_2B의 비율(E_imp/E_2B)로 정규화했을 때, (M_lar−M_tar)/M_pro와 M_s/(M_tar+M_pro)의 변화를 통해 ‘합병’, ‘스키퍼’, ‘침식(erosive)’ 전이를 명확히 구분할 수 있음을 보여준다. 특히 낮은 각도(θ ≈ 30°)에서는 겹치는 부피가 커져 충격 파동에 의한 열에너지 전환이 크게 일어나, 동일한 v_imp에서도 파괴가 더 쉽게 일어난다. 반면 높은 각도에서는 겹치는 부피가 작아 에너지 손실이 적고, 동일한 v_imp에서도 스키퍼가 우세하다.

새로운 분석 모델은 겹치는 질량 M_t,o와 M_p,o를 기하학적 공식(부록 B)으로 계산하고, 이를 통해 겹치는 영역에 전달되는 에너지 E_over = γ′ E_imp을 정의한다. 여기서 γ′는 겹치는 질량 비율에 따라 결정되며, 충돌면에 수직인 성분을 고려해 E_over,⊥ = E_over cos²θ 로 보정한다. 잔여 운동 에너지 E_res = E_imp − E_over,⊥와 중력 결합 에너지 E_2B를 비교해 ‘합병’(E_res ≲ E_2B)과 ‘스키퍼’(E_res ≫ E_2B)를 구분한다.

질량 손실은 두 가지 메커니즘으로 나뉜다. ① 중력 주도 손실은 차원화된 중력 에너지 ˜E_g = E_g/E_over,⊥ 로 표현되며, F_g(x)=1−(1+x/0.15)^{−1} 형태의 함수로 질량 보존율을 근사한다. ② 기화·증발 손실은 특수 에너지 u_cv(주로 맨틀 물질)와 연계된 차원화된 기화 에너지 ˜E_v 로 기술한다. 고에너지 충돌에서는 ˜E_v가 지배해 급격한 질량 감소가 발생한다. 모델은 M_lar = (M_tar+M_pro) F_g(˜E_g,tot) (합병) 혹은 M_lar = M_tar F_g(˜E_g,tar)+M_s,gain (스키퍼) 로 구현되며, 시뮬레이션 데이터와 높은 일치도를 보인다.

또한 각도 가중 평균을 수행해 〈M_lar〉와 〈M_s〉에 대한 반경-각도 통합식을 도출하였다. 이는 행성 형성 N‑body 시뮬레이션에 직접 삽입 가능하도록 반경·속도·각도 분포를 고려한 반경-각도 평균 함수를 제공한다. 마지막으로, 동일한 에너지 기반 모델이 미세먼지 집합체 충돌 실험(먼지 응집)과도 일관된 결과를 보이며, 규모 확장성(scaling)이 뛰어남을 확인했다.

이러한 결과는 기존 Leinhardt & Stewart(2012) 혹은 Genda et al.(2017) 모델이 각도·기화 효과를 충분히 반영하지 못했던 한계를 보완하고, 행성 형성 단계에서 스키퍼 충돌이 차지하는 비중을 정량적으로 평가할 수 있게 한다.