2025-11-22T11:49:15.930592

Dark Matter Subhalos and Higher Order Catastrophes in Gravitational Wave Lensing

Vujeva, Ezquiaga, Gilman et al.

Gravitational lensing is an invaluable probe of the nature of dark matter, and the structures it forms. Lensed gravitational waves in particular allow for unparalleled sensitivity to small scale structures within the lenses, due to the precise time resolution in combination with the continuous monitoring of the entire sky. In this work, we show two distinct ways of using strongly lensed gravitational waves to identify the presence of dark matter subhalos: \emph{i)} through higher order caustics generating high relative magnification ($Î¼_r > 2$), short time delay image pairs that break the caustic universality relations of single dark matter halos, which occur for $\sim 1-10$ percent of strongly lensed events in our cold dark matter models, and \emph{ii)} through the presence of more than three highly magnified images, which occur for $\sim 0.01-1$ percent of the same simulated events. We find that these results are highly sensitive to the concentrations of subhalos in our simulations, and more mildly to their number densities. The presence of low-mass subhalos increases the probability of observing wave-optics lensing in lensed gravitational waves, which is studied by solving the diffraction integral with the stationary phase approximation, as well as numerically. We also report distinct quantitative and qualitative differences in the distributions of relative magnifications and time delays for subhalo populations with increased number densities or concentrations. With the upcoming detection of strongly lensed events by ground- and space- based detectors, comparisons against these simulated distributions will provide insight into the nature of dark matter.

academic

Dark Matter Subhalos and Higher Order Catastrophes in Gravitational Wave Lensing

基本信息

论文ID: 2510.14953
标题: Dark Matter Subhalos and Higher Order Catastrophes in Gravitational Wave Lensing
作者: Luka Vujeva, Jose María Ezquiaga, Daniel Gilman, Srashti Goyal, Miguel Zumalacárregui
分类: astro-ph.CO, astro-ph.GA, gr-qc
发表时间: October 17, 2025
论文链接: https://arxiv.org/abs/2510.14953

摘要

本文研究了引力透镜效应对引力波的影响，特别是暗物质子晕如何通过强引力透镜现象影响引力波观测。研究表明，可以通过两种方式识别暗物质子晕的存在：1）通过高阶焦散产生的高相对放大率（μr > 2）、短时延的像对，这些像对打破了单一暗物质晕的焦散普适性关系；2）通过产生超过三个高度放大的像。研究发现这些结果对子晕的浓度高度敏感，对数密度的敏感性相对较低。

研究背景与动机

问题定义

传统的引力透镜研究主要关注单一、光滑的暗物质晕（如SIS、NFW或SIE轮廓），但实际的星系尺度透镜中包含大量暗物质子晕。这些子晕的存在会显著改变引力波透镜的观测特征，但此前缺乏系统性的研究。

研究重要性

暗物质探测：引力波透镜为探测小尺度暗物质结构提供了独特的工具，能够探测到质量低至10⁷M☉的暗物质子晕
精确时间分辨率：引力波探测器具有毫秒级的时间精度，远超电磁波透镜的天级时延精度
暗物质模型区分：不同暗物质模型（冷、温、模糊暗物质等）在小尺度结构上存在显著差异

现有方法局限性

大多数研究假设光滑的单一透镜
缺乏对子晕影响的系统性建模
电磁波观测受空间分辨率限制，难以区分高阶焦散效应

核心贡献

首次系统研究：首次全面研究了暗物质子晕对强引力波透镜的影响
识别标准：提出了两个明确的子晕识别标准：μr > 2的短时延像对和超过三个高度放大的像
高阶焦散理论：将燕尾和蝴蝶焦散理论应用于引力波透镜分析
波动光学效应：研究了子晕对引力波衍射和干涉现象的影响
参数敏感性分析：量化了观测概率对子晕浓度和数密度的依赖关系

方法详解

任务定义

研究暗物质子晕如何影响强引力波透镜的观测特征，特别是相对放大率（μr）和时延（ΔT）的分布关系。

模型架构

复合透镜模型

主晕模型：使用椭圆奇异等温球面（eSIS）描述星系尺度主晕
- 透镜势：ψ(x⃗) = √(x₁² + x₂²/q²)
- 爱因斯坦半径：~9.4 kpc（透镜红移zL = 0.5）
- 椭率：q = 0.9
子晕分布：使用pyHalo代码生成子晕性质
- 质量函数：d²N/dlogmdA = Σsub(m/m₀)^(-α)
- NFW密度轮廓：ρ(r) = ρs/(r/rs)(1 + r/rs)²
- 浓度-质量关系：采用Diemer-Joyce关系

波动光学处理

衍射积分：F(ω) = (ω/2πi)∫d²x⃗e^(iωT(x⃗,y⃗))
几何光学近似：在高频极限下使用驻相近似
时域放大因子：通过GLoW代码计算I(τ)和G(τ)

技术创新点

高阶焦散分析

燕尾焦散（A₄）：由扰动折叠焦散产生
- 生成函数：V = (1/5)x⁵ + (1/3)ax³ + (1/2)bx² + cx
- 最多产生4个实像
蝴蝶焦散（A₅）：由扰动尖点焦散产生
- 生成函数：V = (1/6)x⁶ + (1/4)ax⁴ + (1/3)bx³ + (1/2)cx² + dx
- 最多产生5个实像

普适性关系的破缺

传统单一透镜的焦散遵循严格的普适性关系：

折叠焦散：μr = 1（恒定）
尖点焦散：μr ≤ 2（上限）

子晕的存在破坏了这些关系，使得μr > 2成为可能。

实验设置

模拟参数

子晕质量范围：10⁷ - 10⁹M☉
浓度模型：
- 基准：Diemer-Joyce关系（cDJ）
- 高浓度：c = 3×cDJ, 6×cDJ
数密度参数：Σsub = 0.025, 0.050, 0.075

选择标准

放大率阈值：|μ| > 5（仅考虑临界曲线附近的子晕）
源面采样：最小放大率μsr > 25的高分辨率区域

计算工具

pyHalo：子晕性质生成
GLoW：引力波透镜计算（包括波动光学效应）

实验结果

主要统计结果

模型	P(μr > 2)	P(Nim > 5)
基准模型	9.0×10⁻³	2.8×10⁻⁴
c = 3×cDJ	8.0×10⁻²	6.4×10⁻³
c = 6×cDJ	8.7×10⁻²	1.6×10⁻²
Σsub = 0.050	1.0×10⁻²	1.7×10⁻⁴
Σsub = 0.075	4.0×10⁻²	3.1×10⁻⁴