A Salpeter IMF and an NFW halo: Disentangling the dark and stellar mass of an elliptical galaxy through precise lens modelling of a double-source-plane system

기술 요약: 잭팟 렌즈(Jackpot Lens)에서의 암흑 물질과 항성 질량 분리

문제 정의
거대 은하의 암흑 물질(DM) 분포를 결정하는 것은 도전적인 과제이다. 암흑 물질은 빛을 내지 않기 때문에 중력적 영향으로부터 추론해야 하기 때문이다. 강한 중력 렌즈 현상은 총 질량에 대한 직접적인 탐사 도구를 제공하지만, 서로 다른 질량 프로파일이 동일한 렌싱 관측치를 생성할 수 있는 질량-시트 변환(mass-sheet transformation, MST) 퇴화 문제로 인해 어려움을 겪는다. 또한, 초기 질량 함수(IMF) 및 잠재적인 반경 방향 구배에 따라 달라지는 항성 질량-광도 비($M/L$)의 불확실성으로 인해 항성과 암흑 물질 성분을 분리하는 작업은 복잡하다. 이중 소스 평면 렌스(DSPL) 시스템인 "잭팟" 렌스(J0946+1006)에 대한 기존 연구들은 내부 암흑 물질 기울기($\gamma_{in}$)에 대해 상충하는 결과를 내놓았다. 일부는 매우 가파른 첨점($\gamma \approx 1.7$)을 시사한 반면, 다른 연구들은 나바로-프렌크-화이트(NFW) 예측과 일치하는 프로파일을 찾아냈다. 이러한 불일치는 상수 항성 $M/L$ 비, 구형 대칭성, 또는 실제의 복잡한 복합 질량 분포를 포착하지 못할 수 있는 단일 멱법칙(single power-law) 프로파일의 사용과 같은 단순화된 모델링 가정에서 기인하는 경우가 많다.

방법론
저자들은 항성과 암흑 물질 성분을 명시적으로 분리하는 복합 질량 모델을 사용하여 J0946+1006에 대한 강한 렌싱 분석을 제시하며, DSPL의 독특한 기하학적 구조를 활용하여 MST 퇴화를 해소한다.

• 데이터: 본 분석은 Hubble 우주 망원경(HST) ACS F814W 이미징과 별의 속도 분산을 측정하기 위한 Multi Unit Spectroscopic Explorer (MUSE) 적분 필드 분광법을 활용한다. 이 시스템은 전경 렌스($z=0.222$)와 두 개의 배경 소스($z=0.609$ (S1) 및 $z=2.045$ (S2))로 구성된다.
• 빛 모델링: 렌스 광도는 15개의 성분을 가진 다중 가우시안 확장(Multi-Gaussian Expansion, MGE)을 사용하여 모델링된다. 소스 광도는 매개변수적(타원형 가우시안) 방식과 비매개변수적(Matérn 공분산을 가진 픽셀화된 가우시안 프로세스) 방식을 모두 사용하여 모델링된다.
• 질량 모델링:
  • 항성 성분: 항성 질량은 전역 질량-광도 비($M/L$)를 통해 빛 모델과 연결되며, 자유로운 반경 방향 구배($\nabla M/L$)를 허용한다.
  • 암흑 물질 성분: 기본 모델은 타원형 일반화 NFW(gNFW) 프로파일을 사용하며, 계산 속도를 높이기 위해 3D MGE(30개의 가우시안)로 근사화하였다. 내부 기울기($\gamma_{in}$)와 척도 반경($r_s$)은 자유 매개변수이다.
  • 대안 모델: 비교를 위해 암흑 물질 성분에 타원형 멱법칙(EPL)을 사용하는 모델(Star+EPL)도 피팅한다.
• 추론 전략: 본 연구는 GPU 가속 및 자동 미분을 위해 JAX 기반으로 구축된 Herculens 코드를 사용한다. 워크플로우는 사후 확률 모드를 찾기 위한 확률적 변분 추론(SVI)을 거쳐, 전체 사후 분포를 매핑하기 위한 깁스 샘플링 내의 해밀토니안 몬테카를로(NUTS)를 포함한다.
• 퇴화 해소: DSPL 기하학은 MST를 제약하는 데 사용되는데, 첫 번째 소스 평면에 적합한 질량 시트는 두 번째 소스의 이미징을 재현하는 데 실패하기 때문이다. 또한, 저자들은 선 시선 속도 분산($\sigma_{los} = 249.1 \pm 2.1$ km s$^{-1}$) 데이터와 $\Lambda$CDM 시뮬레이션으로부터 얻은 농도-질량($c-M$) 사전 확률을 통합하여 사후 확률의 가중치를 재설정한다.

주요 기여

1. DSPL 제약을 이용한 복합 모델링: 본 논문은 두 성분(항성 + gNFW)의 유연한 질량 모델을 DSPL 시스템에 적용하며, 두 번째 소스 평면을 활용하여 운동학적 가정에만 의존하지 않고 MST를 부분적으로 해소한다.
2. 유연한 항성 $M/L$: 상수 항성 $M/L$을 가정했던 이전 분석들과 달리, 본 연구는 반경 방향 구배와 자유로운 전역 정규화를 허용하여 이러한 변화에 대한 IMF 추론의 견고성을 테스트한다.
3. 내부 기울기 재평가: 본 연구는 매우 가파른 내부 암흑 물질 첨점($\gamma \approx 1.7$)에 대한 기존의 주장을, 그러한 가파름이 물리적 실체가 아니라 멱법칙 모델링 가정의 산물일 수 있음을 보여줌으로써 반박한다.

결과

• 항성 질량 및 IMF: 데이터는 대략적으로 일정한 항성 질량-광도 비($\nabla M/L \approx 0$)를 선호한다. 추론된 항성 질량은 $M_\star = 4.4^{+0.25}_{-0.39} \times 10^{11} M_\odot$이다. IMF 불일치 매개변수는 $\alpha_{SPS} = 0.96^{+0.16}_{-0.14}$ (운동학 및 $c-M$ 사전 확률 포함)이며, 이는 Salpeter IMF와는 일치하고 Chabrier IMF와는 약 $3\sigma$ 수준에서 긴장 관계에 있음을 나타낸다.
• 암흑 물질 프로파일: 암흑 물질 헤일로의 내부 기울기는 $\gamma_{in} = 1.04^{+0.10}_{-0.14}$로 제한되며, 이는 표준 NFW 예측($\gamma=1$)과 일치한다. 이는 $\gamma \approx 1.7$의 가파른 결과를 냈던 이전 연구들과 대조된다. 저자들은 이전의 가파른 결과가 EPL 모델의 한계, 즉 gNFW 프로파일이 허용하는 낮은 기울기의 매개 변수 공간에 접근할 수 없기 때문이라고 설명한다.
• 헤일로 질량 및 SHMR: 총 헤일로 질량은 $M_{200} = 1.11^{+0.37}_{-0.32} \times 10^{13} M_\odot$이며, 이에 따른 항성-헤일로 질량비는 $\log_{10}(M_{200}/M_\star) = 1.41^{+0.13}_{-0.14}$이다. 이 값은 유사한 적색편이에서의 일반적인 문헌 관계보다 약 0.1 dex 높으며, 이 차이는 해당 관계의 고유 산포와 맞먹는 수준이다.
• 모델 비교: 암흑 물질을 모델링하기 위해 EPL 프로파일을 강제했을 때(Star+EPL), 모델은 가파른 기울기($\gamma \approx 1.61$)를 회복하여 이전의 결과를 재현한다. 그러나 gNFW 모델은 이러한 가파름이 호(arc)가 조사하는 좁은 반경 범위에 대한 근사적 산물임을 보여준다.

의의 및 주장
본 논문은 DSPL 기하학이 항성 운동학의 기능과 맞먹는 질량 프로파일 제약력을 제공하며, 운동학적 데이터를 얻기 어려운 고적색편이 렌스에 대한 실질적인 대안이 될 수 있다고 주장한다. 결과는 충분한 유연성을 가진 모델을 사용하고 여러 소스 평면을 통해 MST를 제약할 경우, 거대 조기형 은하의 내부 암흑 물질 프로파일이 NFW 예측과 일치함을 시사한다. 저자들은 잭팟 렌스가 본질적으로 가파른 암흑 물질 첨점을 필요로 하지 않으며, 이전 연구에서의 가파름은 모델링 선택의 결과였을 가능성이 높다고 결론짓는다. 그들은 이 접근 방식이 Euclid, 중국 우주 정거장 망원경(CSST), 그리고 Nancy Grace Roman 우주 망원경과 같은 차세대 대규모 조사에서 수행될 대규모 DSPL 분석의 템플릿 역할을 할 것이라고 상정한다. 이러한 조사는 인구 통계적 수준의 연구를 통해 고유한 암흑 물질 헤일로 특성을 통계적으로 회복하고 선택 편향을 해결할 수 있게 할 것이다.