[论文解读] Unified Halo-Independent Formalism Derived From Convex Hulls

Graciela B. Gelmini, Samuel J. Witte+1 more|arXiv (Cornell University)|Jul 21, 2017

Blind Source Separation Techniques被引用 2

一句话总结

本文提出了一种基于凸包理论的统一晕独立形式化方法，通过将速度分布表示为狄拉克δ函数的和，以最大化暗物质探测中的似然函数。该方法表明，仅在特定条件下最佳拟合晕函数是唯一的，否则需构建退化带；并显示假设各向同性会得到与非各向同性拟合不同的分段线性晕函数。

ABSTRACT

Using the Fenchel-Eggleston theorem for convex hulls (an extension of the Caratheodory theorem), we prove that any likelihood can be maximized by either a dark matter 1- speed distribution F(v) in Earth's frame or 2- Galactic velocity distribution fgal(), consisting of a sum of delta functions. The former case applies only to time-averaged rate measurements and the maximum number of delta functions is (−1), where is the total number of data entries. The second case applies to any harmonic expansion coefficient of the time-dependent rate and the maximum number of terms is . Using time-averaged rates, the aforementioned form of F(v) results in a piecewise constant unmodulated halo function 0BF(vmin) (which is an integral of the speed distribution) with at most (-1) downward steps. The authors had previously proven this result for likelihoods comprised of at least one extended likelihood, and found the best-fit halo function to be unique. This uniqueness, however, cannot be guaranteed in the more general analysis applied to arbitrary likelihoods. Thus we introduce a method for determining whether there exists a unique best-fit halo function, and provide a procedure for constructing either a pointwise confidence band, if the best-fit halo function is unique, or a degeneracy band, if it is not. Using measurements of modulation amplitudes, the aforementioned form of fgal(), which is a sum of Galactic streams, yields a periodic time-dependent halo function BF(vmin, t) which at any fixed time is a piecewise constant function of vmin with at most downward steps. In this case, we explain how to construct pointwise confidence and degeneracy bands from the time-averaged halo function. Finally, we show that requiring an isotropic Galactic velocity distribution leads to a Galactic speed distribution F(u) that is once again a sum of delta functions, and produces a time-dependent BF(vmin, t) function (and a time-averaged 0BF(vmin)) that is piecewise linear, differing significantly from best-fit halo functions obtained without the assumption of isotropy.

研究动机与目标

• 开发一种统一的形式化方法，以在不假设特定晕模型的前提下最大化暗物质直接探测中的似然函数。
• 确定在似然最大化过程中，最佳拟合晕函数何时唯一，何时出现退化。
• 基于时间平均和时间依赖的速率测量，构建晕函数的逐点置信带或退化带。
• 探讨假设银河系速度分布各向同性如何改变最终得到的晕函数形状。

提出的方法

• 应用Fenchel-Eggleston定理（其为Carathéodory定理的推广），将似然最大化所需的速度分布表示为有限个狄拉克δ函数的和。
• 推导出：对于时间平均速率，未调制的晕函数BF(vmin)为分段常数，最多有(n−1)个下降台阶，其中n为数据点数。
• 对于时间依赖速率的谐波展开系数，银河系速度分布fgal()被表示为最多n个狄拉克δ函数的和。
• 通过分析最佳拟合晕函数的唯一性，利用凸包技术评估退化性，构建置信带或退化带。
• 证明：若假设银河系速度分布各向同性，则晕函数F(u)为狄拉克δ函数的和，从而得到分段线性的时间依赖型最佳拟合BF(vmin, t)。
• 利用时间平均晕函数，推断调制幅度测量中时间依赖晕函数的置信带或退化带。

实验结果

研究问题

• RQ1 在晕独立似然最大化中，最佳拟合晕函数在何种条件下是唯一的？
• RQ2 当最佳拟合晕函数出现退化时，如何构建置信带？
• RQ3 当将银河系速度分布建模为狄拉克δ函数之和时，时间依赖晕函数BF(vmin, t)的函数形式为何？
• RQ4 假设银河系速度分布各向同性，与非各向同性拟合相比，对最佳拟合晕函数形状有何影响？
• RQ5 能否利用时间平均晕函数推断调制幅度测量中时间依赖晕函数的性质？

主要发现

• 最佳拟合晕函数仅在特定条件下唯一；否则必须使用退化带而非置信带。
• 对于时间平均速率测量，未调制的晕函数BF(vmin)为分段常数，最多有(n−1)个下降台阶，其中n为数据条目数。
• 当对时间依赖速率的谐波展开系数建模时，银河系速度分布fgal()被表示为最多n个狄拉克δ函数的和。
• 从流状fgal()导出的时间依赖晕函数BF(vmin, t)在任一固定时间下，对vmin为分段常数，最多有n个下降台阶。
• 假设各向同性会导致晕函数F(u)为狄拉克δ函数之和，从而产生分段线性的时间依赖型BF(vmin, t)，与非各向同性最佳拟合函数显著不同。
• 可基于所提出的正式化方法，利用时间平均晕函数构建时间依赖晕函数的逐点置信带或退化带。