별 형성의 관측적 제약 — 무엇이 측정값이고, 무엇이 모형을 통과한 2차 산물인가

별 형성 연구에서 가장 자주 인용되는 두 숫자가 있습니다. 은하의 별 형성률(SFR)과 분자기체 질량(M_H2)입니다. 케니컷-슈미트 관계도, 기체 소진 시간도, 자유낙하당 효율도 전부 이 두 값 위에 얹혀 있습니다.

그런데 망원경이 실제로 기록하는 것은 둘 중 어느 것도 아닙니다. 검출기에 남는 것은 특정 파장에서의 유량입니다. 별 형성률도 기체 질량도, 그 유량에 여러 겹의 변환을 적용한 뒤에 나오는 유도량입니다.

이 글은 별 형성의 정량적 결론들이 각각 어떤 변환에 의존하는지, 그리고 그 변환이 흔들릴 때 어떤 결론이 함께 흔들리는지를 정리합니다. 결론을 부정하려는 것이 아니라, 어떤 자릿수가 관측에서 오고 어떤 자릿수가 모형에서 오는지를 분리하려는 것입니다.

관측량과 유도량의 분리

먼저 층위를 나눠 두겠습니다.

표의 왼쪽 열이 논문 초록에 등장하는 양이고, 오른쪽 열이 그 값을 만든 재료입니다. 오른쪽이 바뀌면 왼쪽이 같이 움직입니다. 그리고 오른쪽 대부분은 다른 은하, 다른 금속함량, 다른 밀도 환경에서 그대로 성립한다고 검증된 적이 없습니다.

첫 번째 병목 — CO는 수소가 아니다

분자기체의 대부분은 H2입니다. 그런데 H2는 등핵 이원자 분자라 영구 쌍극자 모멘트가 없고, 저온 분자운의 온도(약 10 K)에서는 사실상 복사하지 않습니다. 가장 많은 물질이 보이지 않는 상태입니다.

그래서 CO를 대신 봅니다. CO의 J=1-0 전이는 저온에서도 여기되고 밝습니다. 문제는 CO 선속에서 H2 질량으로 넘어갈 때 필요한 X_CO 전환계수입니다.

은하계 원반에서 흔히 쓰는 표준값은 약 2×10²⁰ cm⁻² (K km s⁻¹)⁻¹ 이고, 리뷰가 제시하는 불확실도는 대략 ±30% 수준입니다(Bolatto, Wolfire & Leroy 2013). 여기까지는 비교적 안정적입니다.

문제는 이 값이 환경에 따라 계통적으로 이동한다는 점입니다. 금속함량이 낮으면 CO는 광해리로 파괴되는데 H2는 자기차폐로 살아남습니다. 그러면 CO가 없는 H2, 이른바 CO-dark 기체가 생깁니다. 이 경우 X_CO는 표준값보다 커집니다. 반대로 은하 중심부나 병합 은하처럼 온도가 높고 속도 분산이 큰 영역에서는 X_CO가 표준값보다 작아지는 경향이 보고됩니다.

여기서 중요한 점은 방향입니다. X_CO를 잘못 잡으면 기체 질량이 틀리고, 기체 질량이 틀리면 소진 시간과 자유낙하당 효율이 같은 비율로 틀립니다. 즉 “별 형성이 느리다”는 결론의 정량적 부분은 X_CO의 정확도를 넘을 수 없습니다.

두 번째 병목 — 보이는 별은 소수다

별 형성률 지표는 대부분 대질량성을 봅니다. Hα는 O·B형 별이 만든 이온화 광자를 재처리한 결과이고, 원자외선은 수천만 년 수명의 별을, 적외선은 그 별빛을 흡수해 재복사한 먼지를 봅니다.

그런데 질량으로 보면 별의 대부분은 저질량성입니다. 관측되는 것은 소수의 밝은 별이고, 인용되는 것은 전체 질량입니다. 이 간극을 메우는 것이 초기질량함수입니다.

거의 모든 SFR 보정계수는 특정 IMF를 전제로 계산됩니다. 그래서 같은 유량이라도 어떤 IMF를 쓰느냐에 따라 SFR이 달라집니다. 문헌에서 흔히 쓰는 세 가지는 이렇게 갈립니다.

• Salpeter (1955) — 단일 멱함수, dN/dm ∝ m^-2.35. 저질량 쪽까지 그대로 연장하면 저질량성을 과다 예측합니다.
• Kroupa (2001) — 구간별 멱함수. 0.5 M☉ 위에서 기울기 2.3, 0.08~0.5 M☉에서 1.3, 그 아래에서 0.3.
• Chabrier (2003) — 1 M☉ 아래를 로그정규 분포로, 위를 멱함수로.

Salpeter를 Kroupa나 Chabrier로 바꾸면 같은 광도에 대응하는 질량이 대략 1.5~1.7배 줄어듭니다. 이것은 관측이 바뀌어서가 아니라 가정이 바뀌어서 생기는 차이입니다. 서로 다른 IMF로 계산된 SFR을 나란히 놓고 비교하면 안 되는 이유가 여기 있습니다.

케니컷-슈미트 관계 — 지수 1.4의 위상

슈미트(1959)는 별 형성률 밀도가 기체 밀도의 거듭제곱을 따른다고 제안했고, 케니컷(1998)이 이를 은하 전역 면밀도로 관측 검증했습니다. 널리 인용되는 결과는 다음과 같습니다.

Σ_SFR ∝ Σ_gas^N, N = 1.4 ± 0.15

이 관계는 정상 원반은하와 폭발적 별생성 은하를 아우르는 넓은 범위에서 성립합니다. 여기서 두 가지를 구분해야 합니다.

관측적으로 견고한 부분은 “기체가 많은 곳에서 별이 더 많이 만들어진다”는 단조 증가 관계와 그 대략적인 기울기입니다. 이 부분은 추적자를 바꾸어도 살아남습니다.

모형 의존적인 부분은 지수의 정확한 값입니다. Σ_gas는 X_CO에, Σ_SFR은 IMF와 감광 보정에 각각 의존합니다. 두 축이 서로 다른 가정 위에 있으므로, 축 하나에 계통 기울기가 있으면 지수 자체가 이동합니다. 실제로 분해능·추적자·표본을 바꾼 후속 연구들에서 지수는 대략 1에서 2 사이로 흩어집니다.

그리고 결정적으로, 이것은 상관관계이지 인과 기구가 아닙니다. 지수 1.4가 곧 “별 형성이 자유낙하 시간에 의해 조절된다”는 뜻은 아닙니다. 같은 기울기를 만들어 내는 물리 모형이 여럿이며, 관계식 자체는 그중 어느 것도 배제하지 못합니다.

느린 별 형성 — 1퍼센트라는 숫자

별 형성의 가장 오래된 정량적 문제는 속도입니다.

은하계의 분자기체가 전부 자유낙하로 붕괴한다면, 예상되는 별 형성률은 실제 관측값(대략 1~2 M☉/yr 수준)보다 두 자릿수 이상 커야 합니다. 실제로는 그렇게 되지 않습니다.

이 간극을 정량화하는 양이 자유낙하당 효율 ε_ff 입니다. 구름이 한 번의 자유낙하 시간 동안 자기 질량의 몇 퍼센트를 별로 바꾸는가를 나타냅니다. Krumholz & Tan (2007)은 분자운, 적외선 은하, HCN으로 추적한 조밀 기체까지 밀도가 여러 자릿수에 걸친 계들을 함께 놓고, ε_ff가 대체로 몇 퍼센트 이하, 대략 1퍼센트 언저리에 머문다고 보고했습니다.

이 결과의 강한 점은 밀도가 4~5자릿수 바뀌는 동안에도 값이 크게 움직이지 않는다는 것입니다. 계통오차가 한 계에만 있었다면 이런 평탄함이 나오기 어렵습니다.

동시에 약한 점도 분명합니다. ε_ff는 M_gas와 SFR을 모두 필요로 하므로 앞의 두 병목을 그대로 상속합니다. 게다가 t_ff를 계산하려면 밀도를 정해야 하는데, 밀도는 구름의 경계를 어디로 잡느냐에 따라 달라집니다. 관측에서 “구름”은 자연이 준 경계가 아니라 알고리즘이 정의한 덩어리입니다.

IMF는 보편적인가 — 아직 닫히지 않은 논쟁

IMF가 보편적이라는 가정은 위 계산 대부분의 전제입니다. 그래서 이 가정이 흔들리면 파급이 넓습니다.

보편성을 지지하는 쪽의 근거는 은하계 내부입니다. 성단, 성협, 분자운 코어를 개별적으로 셌을 때 나오는 질량함수가 환경 차이에 비해 놀랄 만큼 비슷합니다. 리뷰(Offner et al. 2014)는 관측된 변화의 상당 부분이 표본 크기의 통계 요동과 계통효과로 설명 가능한 범위에 있다고 정리합니다. 성단에서 대질량성이 관측 한계에 걸리거나, 미분해 쌍성이 저질량성을 가리는 효과가 대표적입니다.

변화를 주장하는 쪽은 은하계 밖을 봅니다. van Dokkum & Conroy (2010)는 광도가 높은 타원은하의 적분광에서 저질량 왜성에 민감한 흡선 특징을 검출하고, 이 은하들이 저질량성이 상대적으로 많은(bottom-heavy) IMF를 가진다고 보고했습니다. 사실이라면 그 은하들의 항성질량은 표준 IMF로 얻은 값보다 상당히 커집니다.

이 논쟁이 아직 닫히지 않은 이유는 두 진영이 같은 것을 재고 있지 않기 때문입니다. 은하계 연구는 별을 하나씩 세지만, 타원은하 연구는 적분광의 흡선 지수를 재고 이를 항성 개체군 합성 모형으로 해석합니다. 후자에서는 IMF 효과가 금속함량·원소 존재비 패턴·나이와 축퇴됩니다. 즉 저질량성 초과 신호와 원소 존재비 이상을 관측만으로 완전히 분리하기 어렵습니다.

따라서 현재 상태는 “IMF가 변한다”도 “변하지 않는다”도 아닙니다. 은하계 규모에서는 보편성 가정이 잘 버티고, 극단적 환경으로 나가면 검증되지 않은 영역이라는 것이 정확한 서술입니다.

원시성 광도 문제 — 시간 평균이 숨기는 것

지금까지는 은하 규모였고, 개별 원시성 규모에서는 다른 종류의 제약이 걸립니다.

원시성이 정상적인 강착률로 꾸준히 물질을 받는다면 예상되는 광도가 나옵니다. 그런데 적외선 탐사에서 실제로 측정된 원시성 광도 분포는 예상보다 어두운 쪽으로 치우칩니다(Evans et al. 2009). 이것이 이른바 광도 문제입니다.

유력한 해석은 간헐적 강착입니다. 원시성이 대부분의 시간을 낮은 강착률로 보내다가, 짧은 폭발 구간에 물질의 상당 부분을 삼킨다는 것입니다. FU Ori형 폭발이 그 관측적 대응물로 논의됩니다(Audard et al. 2014).

이 해석이 중요한 이유는 방법론적입니다. 순간 광도를 시간 평균 강착률로 읽으면 안 된다는 뜻이기 때문입니다. 표본이 아무리 커도, 폭발 구간의 지속시간이 짧으면 관측 표본에는 거의 잡히지 않습니다. 관측된 분포가 실제 분포를 대표한다는 가정 자체가 여기서 깨집니다.

이 논리가 적용되지 않는 경우

지금까지의 서술을 “그러므로 모든 수치가 믿을 수 없다”로 읽으면 반대쪽으로 틀립니다. 위 제약들이 성립하지 않거나 약해지는 조건이 분명히 있습니다.

같은 계를 상대 비교할 때는 전환계수가 상쇄됩니다. 같은 은하 안에서 영역 A와 B의 기체 질량비를 볼 때, X_CO가 두 영역에서 같다고 볼 근거가 있다면 그 비율은 X_CO의 절대값과 무관합니다. 계통오차는 절대량에 걸리지 비율에 걸리지 않습니다.

IMF 가정이 상쇄되는 경우도 같습니다. 두 은하의 SFR 비를 같은 IMF로 계산했다면, IMF가 틀렸더라도 비는 대체로 살아남습니다. IMF가 결정적인 것은 절대 질량 규모입니다.

별을 하나씩 셀 수 있으면 IMF 가정이 필요 없습니다. 근거리 성단이나 성협처럼 개별 광도함수를 얻을 수 있는 계에서는 IMF가 가정이 아니라 측정 대상입니다. 이 경우 문제는 IMF가 아니라 완전도 한계와 쌍성 미분해입니다.

자유낙하당 효율은 밀도를 정의할 수 있을 때만 의미가 있습니다. 경계가 뚜렷한 조밀 코어에서는 t_ff가 비교적 잘 정의되지만, 계층적이고 프랙탈에 가까운 분자운 전체에 하나의 t_ff를 부여하는 것은 편의적 처리입니다. 이 경우 ε_ff는 물리량이라기보다 정의에 의존하는 지표로 읽어야 합니다.

그리고 케니컷-슈미트 관계는 예측 도구가 아니라 요약 통계입니다. 개별 영역의 별 형성률을 예측하는 데 쓰면 산포가 크게 나옵니다. 은하 전역으로 평균했을 때 성립하는 관계이며, 분해능을 높일수록 관계는 약해지고 산포는 커집니다.

다음 검증 경로

현재 열려 있는 관측적 진전 경로는 대체로 세 갈래입니다.

첫째, X_CO의 독립 교정. 먼지 연속파, 중성탄소 미세구조선, 감마선 등 CO에 의존하지 않는 추적자로 같은 영역의 기체 질량을 재고 CO 기반 값과 대조하는 접근입니다. 서로 다른 계통을 가진 방법이 같은 값에 수렴하면 그 자체가 강한 증거가 됩니다.

둘째, 분해된 규모에서의 통계. 근방 은하를 개별 분자운 규모로 분해한 목록이 축적되면서(예: PHANGS-ALMA 기반 분자운 목록, Rosolowsky et al. 2021), 구름 단위의 물성과 별 형성 활동을 직접 짝지을 수 있게 되었습니다. 은하 전역 평균에서 보이지 않던 산포의 구조가 여기서 드러납니다.

셋째, 적외선에서의 완전도 개선. 감광에 가려 세지 못하던 원시성과 저질량성을 더 깊이 검출하면, 광도 문제와 IMF 저질량 끝의 계통이 함께 좁혀집니다. 다만 이때에도 파장을 바꾸면 세는 대상이 바뀐다는 점은 남습니다. 서로 다른 파장의 표본을 합칠 때의 교차 교정이 새 계통오차원이 됩니다.

연구자용 Q&A

Q. SFR 값을 인용할 때 IMF를 명시해야 하는 이유는? 보정계수 자체가 IMF에 종속되기 때문입니다. Salpeter 기반 값과 Chabrier 기반 값은 대략 1.5~1.7배 차이가 나며, 이 차이는 관측 오차가 아니라 규약 차이입니다. 표에 IMF를 병기하지 않은 SFR은 다른 연구와 비교 불가능합니다.

Q. ε_ff가 여러 계에서 비슷한 것이 곧 보편 법칙의 증거인가? 직접적인 증거는 아닙니다. 밀도가 수 자릿수 걸친 범위에서 값이 평탄한 것은 시사적이지만, ε_ff는 M_gas·SFR·t_ff 세 유도량의 조합이므로 세 계통이 우연히 서로를 상쇄했을 가능성을 배제하려면 독립 교정이 필요합니다.

Q. 케니컷-슈미트 관계의 지수를 물리 모형의 판별 기준으로 쓸 수 있나? 권장되지 않습니다. 여러 물리 모형이 유사한 기울기를 재현하며, 축 양쪽의 계통 기울기가 지수를 직접 이동시킵니다. 지수보다는 산포의 구조와 분해능 의존성이 더 판별력이 있습니다.

Q. 타원은하의 bottom-heavy IMF 보고를 어떻게 취급해야 하나? 검출된 것은 저질량 왜성에 민감한 흡선 특징이고, IMF는 그 특징을 개체군 합성 모형으로 해석한 결과입니다. 금속함량·원소 존재비·나이와의 축퇴가 남아 있으므로, 현 단계에서는 확정된 IMF 변화가 아니라 모형 의존적 시사로 인용하는 것이 정확합니다.

Q. 분자운 수명 추정치의 주된 불확실도는? 구름의 경계 정의와 추적자 대응입니다. 서로 다른 추적자의 공간 분포에서 통계적으로 시계를 뽑는 방법은 두 추적자가 같은 물리 단계를 대표한다고 가정합니다. 이 가정은 환경마다 다르게 성립합니다.

Q. 원시성 광도 분포에서 시간 평균 강착률을 복원할 수 있나? 분포만으로는 불충분합니다. 폭발의 지속시간과 빈도를 독립적으로 제약해야 하며, 이는 폭발 중인 천체의 표본을 늘리거나 강착 이력의 화학적 흔적을 찾는 방식으로 접근되고 있습니다.

본 글은 공개된 논문과 리뷰에 근거한 해설이며, 인용된 수치는 원 논문의 값과 오차를 그대로 옮긴 것입니다. 논문 그림을 복제하지 않았고, 본문의 표와 설명은 편집팀이 작성했습니다. 연구에 인용하실 때는 원 논문을 직접 확인하시기 바랍니다. 최종 확인일: 2026-08-02.

작성: Inknebula 편집팀 — 측정값과 유도량을 나눠서 적습니다.

1차 출처

• Salpeter, E. E. (1955) ApJ 121, 161. doi:10.1086/145971
• Schmidt, M. (1959) ApJ 129, 243. doi:10.1086/146614
• Kennicutt, R. C. Jr. (1998) ApJ 498, 541–552. doi:10.1086/305588
• Kroupa, P. (2001) MNRAS 322, 231–246. doi:10.1046/j.1365-8711.2001.04022.x
• Chabrier, G. (2003) PASP 115, 763–795. doi:10.1086/376392
• Lada, C. J. & Lada, E. A. (2003) ARA&A 41, 57–115. doi:10.1146/annurev.astro.41.011802.094844
• Krumholz, M. R. & Tan, J. C. (2007) ApJ 654, 304–315. doi:10.1086/509101
• Evans, N. J. II et al. (2009) ApJS 181, 321–350. doi:10.1088/0067-0049/181/2/321
• van Dokkum, P. G. & Conroy, C. (2010) Nature 468, 940–942. doi:10.1038/nature09578
• Kennicutt, R. C. & Evans, N. J. (2012) ARA&A 50, 531–608. doi:10.1146/annurev-astro-081811-125610
• Bolatto, A. D., Wolfire, M. & Leroy, A. K. (2013) ARA&A 51, 207–268. doi:10.1146/annurev-astro-082812-140944
• Offner, S. S. R. et al. (2014) Protostars and Planets VI. doi:10.2458/azu_uapress_9780816531240-ch003
• Audard, M. et al. (2014) Protostars and Planets VI. doi:10.2458/azu_uapress_9780816531240-ch017
• Heyer, M. & Dame, T. M. (2015) ARA&A 53, 583–629. doi:10.1146/annurev-astro-082214-122324
• Rosolowsky, E. et al. (2021) MNRAS 502, 1218–1245. doi:10.1093/mnras/stab085

⚠️ 미인용 처리

• 은하계 별 형성률의 구체적 최신값과 오차 — 리뷰마다 표본·보정이 달라 본문에서는 자릿수(1~2 M☉/yr 수준)로만 서술했습니다.
• 자유낙하 예상 별 형성률과 실측의 정확한 비율 — 문헌마다 가정이 달라 “두 자릿수 이상”으로만 적었습니다.
• 분자운 수명의 구체적 추정치(수 Myr~수십 Myr) — 방법·표본 의존성이 커서 수치를 넣지 않고 불확실도 구조만 서술했습니다.
• X_CO 환경 의존성의 정량적 스케일링 관계 — 금속함량 의존 지수는 문헌 간 차이가 커서 방향(증가/감소)만 적었습니다.
• Chabrier IMF의 특성질량 수치 — 단일별 IMF와 계 IMF에서 값이 달라 본문에서는 함수 형태만 서술했습니다.