Assembly bias and the redshift evolution of intrinsic alignments for LRGs

기술 요약: LRG에 대한 조립 편향(Assembly bias) 및 적색편이 진화

문제 제기
은하의 고유 정렬(Intrinsic alignments, IA)은 약한 중력 렌즈 효과 조사에서 주요한 천체물리학적 계통 오차 요인을 구성하며, 적절히 모델링되지 않을 경우 우주론적 파라미터 추론에 편향을 일으킬 수 있다. NLA(Non-Linear Alignment) 모델과 같은 현재의 모델링 노력은 흔히 정렬 진폭이 오직 헤일로 질량에만 의존한다고 가정한다. 그러나 Euclid 및 LSST와 같이 통계적 정밀도가 높은 차세대 조사의 등장과 함께, 관측적 불확실성에 의해 이전에 숨겨져 왔던 이차적 의존성들이 유의미해질 수 있다. 특히, 헤일로 질량 이외의 형성 이력에 따른 헤일로 특성의 의존성인 조립 편향(assembly bias)의 영향과 정렬 신호의 적색편이 진화에 대한 조사가 필요하며, 이를 위해 대규모 유체역학 시뮬레이션을 이용한 조사가 요구된다.

방법론
저자들은 중력적 피드백 변동(fiducial feedback variation)과 바리온 물리학을 포함하는 FLAMINGO 코스몰로지컬 유체역학 시뮬레이션(Schaye et al. 2023)의 2.8 Gpc 박스 크기 실행 버전을 활용한다. 본 연구는 $M_h > 10^{12} M_\odot$인 헤일로와 300개 이상의 항성 입자($\sim 2.22 \times 10^{11} M_\odot$)를 가진 은하로 정의되는 LRG(Luminous Red Galaxy) 유사 중심 은하 샘플에 초점을 맞춘다.

주요 방법론적 단계는 다음과 같다:

1. 형성 적색편이 정의: 형성 적색편이($z_f$)는 헤일로가 현재 질량의 절반($0.5 M_{200\text{mean}}$)을 가졌던 시점의 적색편이로 정의된다. 이 정의는 헤일로 농축도($c_{200\text{crit}}$)와 강한 상관관계를 보임으로써 검증되었다.
2. 상관 함수 추정: 일반화된 Landy-Szalay 추정자를 사용하여 위치-위치($gg$) 및 위치-형태($g+$) 상관 함수를 계산하고, IATheory 코드를 통해 NLA 프레임워크 내에서 모델링한다.
3. 조립 편향 분석: 고정된 질량에서 정렬 진폭($A_1$)이 $z_f$에 의존하는지를 격리하기 위해 샘플을 헤일로 질량과 형성 적색편이별로 빈(bin)을 나눈다.
4. 적색편이 진화 모델링: 다양한 시뮬레이션 스냅샷($z=0$ ~ $z=1.5$)에 걸쳐 정렬 진폭과 헤일로 질량 사이의 멱함수 관계($A_1 \propto M^{\beta_M}$)를 피팅한다.
5. 은하 추적: 단순한 질량 성장으로부터 적색편이 진화를 분리하기 위해, HBT-HERONS 서브헤일로 파인더의 Track_id를 사용하여 중심 은하들을 여러 스냅샷에 걸쳐 추적한다. 샘플은 초기 적색편이($z_i$)에서 선택되며, 병합되거나 항성 입자 임계값 미만으로 떨어지는 객체는 제거된 채 $z=0$까지 추적된다.

주요 결과

• 조립 편향 의존성: 정렬 진폭은 헤일로 질량뿐만 아니라 형성 적색편이에도 의존한다. 고정된 헤일로 질량에서, 더 일찍 형성된(높은 $z_f$) 헤일로와 그 중심 은하들은 체계적으로 더 높은 정렬 진폭을 보인다. 이 효과는 $A_1$–질량 관계에서 수직 이동으로 나타나지만, 멱함수의 기울기는 거의 변하지 않는다.
• 파라미터의 적색편이 진화: 질량-정렬 멱함수의 진폭($\alpha_M$)은 적색편이가 증가함에 따라 꾸준히 증가한다. 기울기($\beta_M$)는 $z \sim 1$까지 감소하다가 이후 평탄해진다. 저자들은 $\alpha_M(z)$와 $\beta_M(z)$에 대한 경험적인 이차식 피팅(식 13 및 14)을 제공한다.
• 질량-적색편이 모델: 주어진 적색편이와 평균 헤일로 질량을 가진 LRG 샘플의 정렬 진폭을 예측하기 위해 결합된 질량-적색편이 모델(식 12–14)이 제안된다.
• 추적 분석: 높은 초기 적색편이($z_i$)로부터 은하를 추적할 때, 정렬 진폭은 질량 축적과 일치하게 적색편이가 감소함에 따라 증가한다. 그러나 질량-적색편이 모델 예측값으로 정규화한 후에도 체계적인 오프셋이 남아 있다: 높은 $z_i$에서 추적된 은하들은 낮은 $z_i$에서 추적된 은하들보다 더 높은 정렬 진폭을 유지한다. 이는 조립 편향이 주어진 질량과 적색편이에서 정렬 진폭에 환원 불가능한 산포(scatter)를 도입함을 확인시켜 준다.
• 멱함수 편차: 멱함수 피트는 저질량($< 10^{13} M_\odot$) 및 고질량($> 3 \times 10^{14} M_\odot$) 영역에서 벗어난다. 저질량 영역에서는 일부 중력 전용 결과와 대조적으로 감소가 관찰되는 반면, 고질량 영역에서는 지속적인 헤일로 형성에 기인한 감소가 관찰된다.

의의 및 주장
본 논문은 조립 편향이 고정된 헤일로 질량에서의 진폭에 근본적인 퍼짐(spread)을 도입하기 때문에, 약한 중력 렌즈 조사에서의 IA 모델링이 정렬 진폭에 대해 임의로 작은 사전 범위(prior range)에 의존해서는 안 된다고 주장한다. 기존 모델(NLA 등)의 함수 형태가 조립 편향을 설명하기 위해 수정될 필요는 없을지라도, 이 이차적 의존성을 수용하기 위해 진폭 파라미터에 대한 사전 제약(prior constraints)은 넓혀져야 한다.

또한, 본 결과는 중력 전용 시뮬레이션을 이용한 HOD(Halo Occupation Distribution) 기반 정렬 모델의 구현을 복잡하게 만든다. HOD 모델은 일반적으로 질량만을 기준으로 특성을 할당하므로, 명시적으로 포함되지 않는 한 형성 시기 의존성을 포착하지 못한다. 본 연구는 IA의 적색편이 진화 모델을 위한 시뮬레이션 기반 가이드라인을 제공하며, LRG 유사 샘플에 대한 질량-적색편이 의존성 정렬 진폭에 대한 경험적 피팅을 제공한다.

한계점
저자들은 FLAMINGO 시뮬레이션의 해상도가 은하 디스크의 진화를 해결하기에는 불충분하다(디스크 분율이 $z \sim 1.5$까지 $1\%$ 미만으로 유지됨)고 언급한다. 따라서 본 결과는 순수 LRG 샘플(타원 은하)에 엄격히 적용되며, 디스크 은하로의 타원 은하 진화 또는 디스크 은하의 정렬 신호는 고려하지 않는다. 질량-적색편이 모델은 $10^{13}$에서 $10^{14.5} M_\odot$ 사이의 헤일로 질량과 $z \sim 1.5$까지의 적색편이에 대해 유효하다.