Gravitational-Lensing-Dark-Matter-Mapping

以星光的彎曲稱量不可見之物

暗物質不發光、不吸收、不散射電磁波——它與可見世界的唯一互動方式是引力。愛因斯坦的廣義相對論告訴我們，質量彎曲時空，而光線沿著彎曲時空的測地線傳播。因此，當來自遙遠背景星系的光線穿過一個前景大質量結構（如星系團）時，光線的路徑會被彎曲，背景星系的影像在我們眼中變得被輕微拉伸並沿切線方向對齊——這就是弱引力透鏡（Weak Gravitational Lensing）效應。

弱透鏡是當代宇宙學中最有力的暗物質測繪工具，因為它不依賴任何關於暗物質粒子本性的假設——它只測量總質量的投影密度，無論這些質量是由重子物質還是暗物質構成。透過對數百萬個背景星系的形狀進行統計分析，天文學家可以重建前景宇宙結構的二維質量分佈圖，其精確度僅受觀測星系數量的泊松噪聲限制。

剪切場與質量重構

弱透鏡效應的觀測量是剪切（shear）γ——背景星系橢圓率的切向分量。但單一星系的內在橢圓率（intrinsic ellipticity）遠大於弱透鏡剪切訊號（σ_intrinsic ≈ 0.3 vs γ ≈ 0.01–0.1），因此必須對大量背景星系進行統計平均才能提取剪切場。這就是所謂的「形狀噪聲」（shape noise）——弱透鏡測量的根本統計限制。

從剪切場γ(θ)到質量面密度κ(θ)（會聚，convergence）的重構是一個二維反演問題，其數學基礎是Kaiser & Squires (1993)的反演公式。在傅立葉空間中，κ與γ之間存在簡單的代數關係——κ̃(ℓ)正比於γ̃(ℓ)與一個核函數的卷積。這一反演的挑戰在於質量-片簡併（mass-sheet degeneracy）——κ的零模無法從剪切場中唯一確定。

從星系團到宇宙網

弱透鏡的應用跨越了從星系團尺度（~10¹² M☉）到宇宙大尺度結構（~10¹⁵ M☉）的廣闊質量範圍。在星系團尺度，弱透鏡可以用來校準質量-可觀測量關係（如質量-X射線溫度關係），從而將星系團巡天轉化為精確宇宙學約束。在宇宙學尺度，宇宙剪切（cosmic shear）——大尺度結構對背景星系的集體弱透鏡效應——是測量σ₈（物質密度擾動振幅）與Ω_m（物質密度參數）的最靈敏探針之一。

Euclid衛星、Vera C. Rubin天文台（LSST）與羅曼太空望遠鏡（Roman Space Telescope）將在未來十年內將弱透鏡星系的觀測數量從數百萬提升至數十億，使宇宙剪切測量的統計精度達到亞百分比級，足以區分不同的暗能量模型與修改引力理論。

弱透鏡剪切場模擬

以下Python程式演示了簡單的NFW暗物質暈弱透鏡剪切場計算。

import numpy as np

class NFWLensingModel:
    """NFW (Navarro-Frenk-White)暗物質暈的弱透鏡模型"""

    def __init__(self, M200=1e14, c=5.0, z_lens=0.3, z_src=1.0):
        self.M200 = M200      # 暈質量 (M☉)
        self.c = c            # 集中度參數
        self.z_lens = z_lens
        self.z_src = z_src

    def sigma_crit(self):
        """計算臨界面密度 (簡化宇宙學)"""
        D_l = self.z_lens * 3000 / 0.67  # Mpc (近似)
        D_s = self.z_src * 3000 / 0.67
        D_ls = D_s - D_l
        return (1.66e-27 * D_s) / (D_l * D_ls)  # kg/m²

    def shear_profile(self, r):
        """計算NFW暈在投影半徑r處的剪切訊號 (簡化)"""
        r_s = self.M200 / self.c  # 特徵半徑 (kpc)
        x = r / r_s
        # NFW的表面密度核
        def f(x):
            if x < 1:
                return 1 - (2/np.sqrt(1-x**2)) * np.arctanh(np.sqrt(1-x**2))
            elif x == 1: return 1/3
            else: return 1 - (2/np.sqrt(x**2-1)) * np.arctan(np.sqrt(x**2-1))
        delta_sigma = 1e12 * np.array([f(xi) for xi in x])
        return delta_sigma / self.sigma_crit()

model = NFWLensingModel(M200=1e14, c=5.0)
r = np.logspace(0, 3, 50)
shear = model.shear_profile(r)
print(f"峰值剪切: {np.max(shear):.4f} at r={r[np.argmax(shear)]:.0f} kpc")

超越標準模型的探針

弱透鏡不僅可以測繪暗物質的空間分佈，還可以檢驗暗物質的粒子本質。如果暗物質具有非零的自相互作用截面（self-interacting dark matter），星系團核心的暗物質密度剖面將偏離標準的NFW尖峰（cusp）分佈，呈現較為平坦的核心（core）。弱透鏡對星系團質量分佈的高解析度測量是區分這些暗物質模型的最直接手段之一。