분자 구름 속 먼지 응집과 파편화: 충돌이 입도 분포에 미치는 영향

초록

이 연구는 10 K, 밀도 10³–10⁷ cm⁻³ 범위의 정적 분자 구름에서 먼지 입자들의 충돌·응집·파편화를 상세히 모델링한다. 충돌 결과를 수치 시뮬레이션으로 만든 표를 이용해 코어그레이션 방정식을 풀어 입도 분포의 시간 변화를 추적한다. 초기에는 응집이 지배해 입도가 크게 성장하지만, 일정 시점에서 충돌 속도가 파편화 임계치를 넘으면 파편화가 우세해져 성장은 멈추고 작은 입자가 재생산된다. 최종적으로는 로그 크기 구간당 질량이 거의 일정한 평형 상태에 도달한다. 구름의 수명이 자유 낙하 시간 수준이면 변화가 미미하지만, 장기적인 지지 메커니즘이 존재하면 불투명도가 크게 감소한다.

상세 분석

본 논문은 두 단계로 구성된 코어그레이션 모델을 제시한다. 첫 번째 단계는 충돌 속도가 낮아 입자들이 주로 ‘스틱(sticking)’하여 부피와 내부 공극률이 증가하는 응집 단계이며, 두 번째 단계는 충돌 에너지가 파편화 임계값을 초과해 ‘샤터링(shattering)’, ‘브레이크(breakage)’, ‘에로전(erosion)’ 등 다양한 파편화 메커니즘이 작동하는 단계이다. 저자들은 개별 충돌 실험을 고해상도 디스크리트 엘레멘트 메서드(DDEM)와 분자 동역학(MD) 시뮬레이션으로 수행하고, 입자 재료(실리카, 수산화물 등)와 구조(프랙털 차원 D_f≈2.0, 2.5 등)에 따른 충돌 결과 매트릭스를 구축하였다. 이 매트릭스는 충돌 질량비, 상대 속도, 입자 내부 결합 강도 등을 입력으로 받아 ‘완전 응집’, ‘부분 파편화’, ‘전면 파편화’ 등 5가지 구역으로 구분된 결과를 제공한다.

코어그레이션 방정식은 질량 보존을 전제로 연속적인 질량 구간을 설정하고, 충돌 커널 K(m_i,m_j)=π(a_i+a_j)²Δv_ij·S_ij 로 정의한다. 여기서 Δv_ij는 난류와 열운동에 의해 결정되는 상대 속도이며, S_ij는 충돌 성공 확률(응집 혹은 파편화 비율)이다. 저자들은 온도 10 K에서 열운동에 의한 속도가 약 10 cm s⁻¹ 수준이며, 밀도 10⁵ cm⁻³ 이상에서는 난류에 의한 속도가 100 cm s⁻¹까지 상승한다는 점을 강조한다. 이러한 속도 스케일은 파편화 임계 에너지(≈10⁻⁸ erg)와 비교해볼 때, 밀도 10⁶ cm⁻³ 이상에서 파편화가 급격히 시작됨을 의미한다.

시뮬레이션 결과는 질량 분포가 초기에는 좁은 로그 정규형을 보이다가 평균 질량 ⟨m⟩이 시간에 따라 지수적으로 증가한다는 것을 보여준다. 그러나 ⟨m⟩이 특정 임계값을 초과하면 파편화가 지배적으로 전환되어 질량 스펙트럼이 ‘플랫(flat)’ 형태, 즉 로그 크기당 질량이 거의 일정한 상태에 수렴한다. 이 평형 상태는 ‘자기조절(self‑regulation)’ 메커니즘으로 해석될 수 있는데, 큰 입자는 파편화를 통해 작은 입자를 지속적으로 공급하고, 작은 입자는 다시 응집을 통해 성장한다는 순환이 이루어진다.

구름 수명에 대한 의존성도 중요한 결과로 제시된다. 자유 낙하 시간(~10⁵ yr) 수준에서는 응집이 충분히 진행되지 않아 입도 분포 변화가 미미하고, 광학적 깊이(κ)도 거의 변하지 않는다. 반면, 수백만 년에 이르는 장기 안정 구름에서는 평균 입자 크기가 수십 마이크론까지 성장하고, 전체 질량의 30–50 %가 1 µm 이상 입자에 집중된다. 이때 가시광·근적외선 영역의 평균 불투명도는 2배 이상 감소한다. 이러한 결과는 별 형성 초기 단계에서 코어 온도와 화학 반응률에 직접적인 영향을 미칠 수 있음을 시사한다.