[论文解读] Waveform Reconstruction of Core-Collapse Supernovae Gravitational-Waves with Ensemble Empirical Mode Decomposition

Yong Yuan, Xi-Long Fan+3 more|arXiv (Cornell University)|Sep 12, 2023

Gamma-ray bursts and supernovaePhysics and Astronomy被引用 3

一句话总结

本论文提出使用集合经验模态分解（EEMD）从核心坍缩超新星（CCSN）重建引力波（GW）波形，以应对其复杂时频结构带来的挑战。通过分解磁旋转和中微子驱动机制产生的模拟GW信号，研究显示前六阶固有模态函数（IMFs）之和可在匹配分数阈值为0.75时实现波形重建，当归一化振幅为5×10⁻²¹时，重建的误报概率（FAPR）最低可达4×10⁻³，实现约37 kpc的重建距离。

ABSTRACT

The gravitational waves (GW) from core-collapse supernovae (CCSN) have been proposed as a probe to investigate physical properties inside of the supernova. However, how to search and extract the GW signals from core-collapse supernovae remains an open question due to its complicated time-frequency structure. In this paper, we apply the Ensemble Empirical Mode Decomposition (EEMD) method to decompose and reconstruct simulated GW data generated by magnetorotational mechanism and neutrino-driven mechanism within the advanced LIGO, using the match score as the criterion for assessing the quality of the reconstruction. The results indicate that by decomposing the data, the sum of the first six intrinsic mode functions (IMFs) can be used as the reconstructed waveform. To determine the probability that our reconstructed waveform corresponds to a real GW waveform, we calculate the false alarm probability of reconstruction (FAPR). By setting the threshold of the match score to be 0.75, we obtain FAPR of GW sources at a distance of 5 kpc and 10 kpc to be $6 imes10^{-3}$ and $1 imes10^{-2}$ respectively. If we normalize the maximum amplitude of the GW signal to $5 imes10^{-21}$, the FAPR at this threshold is $4 imes10^{-3}$. Furthermore, in our study, the reconstruction distance is not equivalent to the detection distance. When the strain of GW reaches $7 imes 10^{-21}$, and the match score threshold is set at 0.75, we can reconstruct GW waveform up to approximately 36 kpc.

研究动机与目标

• 为解决核心坍缩超新星（CCSN）引力波信号重建的挑战，其具有复杂、非平稳的时频结构。
• 评估集合经验模态分解（EEMD）在从模拟数据中重建CCSN GW波形方面的性能。
• 利用重建的误报概率（FAPR）量化重建波形的可靠性。
• 确定在真实信噪比条件下，EEMD可成功重建CCSN GW信号的最大距离。
• 比较不同CCSN机制和信号振幅（包括归一化振幅）下的重建性能。

提出的方法

• 对由磁旋转和中微子驱动机制生成的核心坍缩超新星模拟引力波信号应用EEMD进行分解。
• 通过累加EEMD识别出的前六阶固有模态函数（IMFs）来重构原始波形。
• 将重构波形与注入波形之间的匹配分数作为重建质量的主要度量指标。
• 计算重建的误报概率（FAPR）以评估重构信号的统计显著性。
• 将GW信号的最大振幅归一化为5×10⁻²¹，以消除FAPR评估中振幅依赖性偏差。
• 改变源距离（5 kpc、10 kpc，直至37 kpc），以确定在匹配分数阈值为0.75时的有效重建范围。

实验结果

研究问题

• RQ1 EEMD能否有效重建具有复杂时频特征的核心坍缩超新星引力波信号？
• RQ2 使用EEMD进行准确波形重建时，最优的固有模态函数（IMFs）数量是多少？
• RQ3 在不同距离和振幅下，EEMD重构的CCSN信号的重建误报概率（FAPR）是多少？
• RQ4 与先进LIGO的探测距离相比，重建距离如何？
• RQ5 EEMD能否在不同CCSN机制和信号振幅（包括归一化振幅）下可靠地重构波形？

主要发现

• EEMD分解所得前六阶固有模态函数（IMFs）之和可提供最精确的核心坍缩超新星引力波波形重建。
• 在匹配分数阈值为0.75时，模拟信号在5 kpc处的FAPR为1×10⁻²，10 kpc处为3×10⁻²。
• 在将最大振幅归一化为5×10⁻²¹后，FAPR降低至4×10⁻³，表明该重建方法具有高可靠性。
• 当应变达到7×10⁻²¹且匹配分数阈值为0.75时，基于EEMD的重建方法可成功实现约37 kpc距离的GW波形重建。
• 重建距离（37 kpc）小于探测距离，表明重建能力受限于信号质量，而不仅限于可检测性。
• 由于振幅较低和探测器响应特性，约一半的模拟波形在5 kpc和10 kpc处无法成功重建。