BLAST 지도에서 픽셀 히스토그램을 이용한 서브밀리미터 은하 수 카운트 추정

초록

BLAST가 관측한 250·350·500 µm 파장의 혼잡 제한 지도에 대해 픽셀 히스토그램(P(D)) 최대우도 분석을 적용하였다. 소스 추출 대신 플럭스 밀도 노드에 전력법칙을 연결해 밝은 소스부터 혼잡 한계 이하까지 연속적인 수 카운트를 얻었으며, 250 µm에서는 –3.7, 350·500 µm에서는 –4.5의 급격한 기울기를 확인하고 0.015 Jy 이하에서는 완만해지는 전이를 보였다. 불확실성 및 상관관계 추정 방법도 제시해 모델 적합에 활용 가능하도록 하였다.

상세 분석

본 논문은 BLAST( Balloon‑borne Large Aperture Submillimeter Telescope) 가 수행한 약 10 deg² 규모의 GOODS‑S 필드 관측 데이터를 대상으로, 전통적인 소스 추출 방식이 갖는 혼잡(confusion) 한계와 선택 편향을 회피하기 위해 픽셀 강도 분포인 P(D) 히스토그램을 이용한 최대우도(MLE) 방법을 적용하였다. P(D) 분석은 각 픽셀에 기록된 플럭스 값이 실제 은하군의 플럭스 분포와 관측 잡음의 합성 결과라는 가정 하에, 가설적인 수 카운트 모델(플럭스‑노드와 그 사이를 연결하는 전력법칙)로부터 기대되는 픽셀 히스토그램을 계산하고, 실제 히스토그램과의 로그우도 차이를 최소화하는 파라미터를 찾는다.

핵심적인 기술적 요소는 다음과 같다. 첫째, 플럭스 노드를 자유롭게 배치하고 각 구간마다 독립적인 전력 지수(α)를 부여함으로써, 기존의 고정된 파워‑로우 모델보다 유연한 형태의 수 카운트를 재구성한다. 둘째, 관측 잡음은 지도 전체에 걸쳐 측정된 RMS 값을 이용해 가우시안 분포로 모델링하고, 지도 내 비균일한 노이즈 레벨을 반영하기 위해 각 픽셀별 가중치를 적용한다. 셋째, 우도 함수는 다중 파라미터 공간에서 뉴턴‑랩슨 혹은 마르코프 체인 몬테카를로(MCMC) 샘플링을 통해 최적화되며, 이 과정에서 파라미터 간 상관관계와 공분산 행렬을 동시에 추정한다.

결과적으로 250 µm 파장에서는 플럭스 범위 0.02–0.5 Jy에서 수 카운트가 dN/dS ∝ S⁻³·⁷ 정도의 급격한 감소를 보였으며, 350·500 µm에서는 dN/dS ∝ S⁻⁴·⁵ 로 더욱 가파른 경사를 나타냈다. 그러나 플럭스가 약 0.015 Jy 이하로 내려가면 전이가 일어나, 기울기가 완만해지는 ‘플랫’ 구간이 관측되었다. 이러한 전이는 은하의 진화 모델에서 저광도 은하군이 지배적인 역할을 함을 시사한다. 또한, 불확실성 추정에서는 포아송 잡음과 지도 내 비균일한 노이즈를 모두 고려한 풀리시안(Full‑Posterior) 분포를 사용했으며, 파라미터 간 강한 양의 상관관계가 존재함을 확인했다.

이 방법은 특히 Herschel와 같은 차세대 서브밀리미터 관측기의 대용량 지도에 적용할 경우, 소스 추출에 의존하지 않고 전체 플럭스 분포를 직접적으로 해석할 수 있어, 은하 형성·진화 모델의 정밀 검증에 유용하다.