1. “$m<1\,\mathrm{eV}$ 就是 wave-like”需要怎样限定？

• 真正的局域判据首先是相空间占据数，而不是单独的质量截断。令每个内部自由度的典型相干胞占据数为 $\mathcal N\sim n_{\rm DM}\lambda_{\rm coh}^3/g_s$。对 $\rho_{\rm DM}=0.4\,{\rm GeV\,cm^{-3}}$、Standard Halo Model 的 $v_0=220\,{\rm km\,s^{-1}}$，Cheong 等得到
  $$\mathcal N\simeq {1.22\times10^{29}\over g_s}\left({1\,\mu{\rm eV}\over m}\right)^4,$$
  所以 $m=1\,{\rm eV}$ 时仍有 $\mathcal N\simeq1.2\times10^5/g_s$，而 $\mathcal N\sim1$ 的波—粒边界约在 $m\simeq18.7\,{\rm eV}/g_s^{1/4}$。因此“远低于 eV”是很安全的高占据区，但 $1\,\mathrm{eV}$ 不是基本界线；边界还依赖局域密度、速度分布、简并度和探测器分辨率。[1]
• 还要比较 de Broglie／相干尺度与探测器尺度。固定速度粒子的 $\lambda_{\rm dB}=2\pi/(mv)$，而相干长度由速度分布导致的去相位决定，两者只在速度宽度与平均速度同阶时才相差 $O(1)$。相干时间由线宽 $\Delta\omega\sim m v^2$ 决定，参数上 $\tau_{\rm coh}\sim1/(m v^2)$。[7] 对上述 halo 模型，严格按一阶相关函数定义可得 $\tau_{\rm coh}\simeq2.80\,{\rm ms}(1\,\mu{\rm eV}/m)$、相干长度 $d_{\rm coh}\simeq674\,{\rm m}(1\,\mu{\rm eV}/m)$，相应品质因子 $Q_{\rm DM}\simeq6.8\times10^5$。[1]
• 高占据仍不是“纯 coherent state”的判据。若某组模的密度算符有正的 $P(\alpha)$，正常排序的相关函数可以由 classical stochastic amplitudes 重现；但 number state 的 $P$ 分布是奇异且非正的，即使 $N\gg1$ 也不因此成为 Glauber coherent state。标准搜索把暗物质写成带随机 Rayleigh 振幅和均匀随机相位的平面波叠加，是一个额外的统计假设，而不是从 $m<1\,\mathrm{eV}$ 自动推出的定理。[1][7]

2. 必须分开的四个“态”层次

1. 宇宙产生机制给出的全局量子态。misalignment axion/ALP 的均匀零模在开始振荡后常被称作宏观 coherent field；Nelson 与 Scholtz 明确说可把标量零模“看作”宏观粒子的 coherent state，但这是一种对经典背景的量子解释，并不证明今天银河晕每个局域模仍处于单一纯 coherent state。[4] 对暴胀期的轻标量扰动，线性量子演化会产生强 two-mode squeezing；大 squeezing 使场振幅表现为近经典随机变量，但 squeezed state 与 Glauber coherent state 不是同一量子态。[3]
2. 一个银河晕或一个相干 patch 的局域态。引力聚集、模混合、环境退相干和只保留实验可达自由度都会改变原初态。Cheong 等给出的 Gaussian density operator
   $$\rho_k=\int {d^2\alpha\over\pi N_k}\,e^{-|\alpha|^2/N_k}\,|\alpha\rangle\langle\alpha|$$
   是局域 virialized halo 的合理候选模型，但作者明确强调尚无从宇宙产生到局域密度矩阵的完整证明。它是许多 coherent states 的混合态，单模上具有 chaotic/thermal 型的数统计，并非一个未知相位的单一纯 coherent state。[1]
3. 对未知相位、速度分布和 halo realization 平均后的 reduced/statistical state。即使每次 realization 是 $|\alpha e^{i\theta}\rangle$，若实验无相位参考并对 $\theta$ 平均，就得到
   $$\rho_{\rm ph}=\int_0^{2\pi}{d\theta\over2\pi}|\alpha e^{i\theta}\rangle\langle\alpha e^{i\theta}|=e^{-N}\sum_{n=0}^{\infty}{N^n\over n!}|n\rangle\langle n|.$$
   这消除了 number-basis 的非对角元，却仍保持 Poisson 计数和 $g^{(2)}(0)=1$。它与上一项 Gaussian $P$ mixture 不同：后者对 $|\alpha|$ 也作 Rayleigh/Gaussian 平均，单模理想极限给出 bunching $g^{(2)}(0)=2$。这是本任务根据密度矩阵直接展开得到的区分；Cheong 等的 Eq. (7) 与 Sec. V 给出后一种 Gaussian mixture 及其 $g^{(2)}$。[1]
4. classical stochastic field approximation。它是计算某组低阶相关函数的有效表示，不必等同于任何特定纯态。随机波模型写成速度模、随机相位和 Rayleigh 振幅的和，可直接给出线宽、功率谱以及跨传感器协方差；但仅拟合这些量通常只能识别统计相关函数，不能反演密度算符的全部非对角元或判定纯度。[1][7]

3. 不同候选粒子的文献依据与 caveat

• misalignment axion/ALP：经典均匀初值和宏观占据支持“condensate/coherent field”的语言，且局域实验常用随机速度模叠加；然而从原初零模到 virialized halo 的局域纯态没有已完成的推导。最保守的说法是：misalignment 给出高占据经典背景的初始条件，不自动选定今天的局域 Glauber 纯态。[4][1]
• 暴胀产生的轻标量：在近自由、线性演化下，跨 $\mathbf k$ 与 $-\mathbf k$ 的模式可形成强 squeezed state；大 squeezing 会压低某一正交分量并给出固定的 classical phase-space trajectory。重返视界、辐射/物质环境和非线性成晕可能导致退相干或模混合，所以不能把 primordial squeezing 直接宣称为银河实验室里可观测的 nonclassical squeezing。[3]
• hidden photon：暴胀涨落机制主要产生纵向模；Graham、Mardon 与 Rajendran 从量子真空涨落/Bogoliubov 演化得到高占据 Gaussian 场，并在重返视界后用 coherent classical-field modes 传播其功率谱。论文同时指出，最简单的 canonical massive vector 的均匀 misalignment 机制并不像标量那样有效，除非引入曲率耦合或其他 UV 机制。[5] 把其原初 Bogoliubov Gaussian 态称为 squeezed、同时把今天局域 halo 另作随机场处理，是本任务据此作出的分层推断；把两者进一步等同于局域 pure coherent state 则没有文献证明。
• number/Fock-like state：它不是通常 halo 形成模型的默认预测，但作为反例非常重要：它说明 $N\gg1$ 与 classical wave-like observables 并不唯一指定量子态。只要实验只看两点函数，就可能与具有相同 $N_k$ 的其他态退化；需要数涨落或更高阶相关才能排除。[1]

4. 单个选定模式的量子统计：哪些相同，哪些不同？

以下表格是本任务的标准量子光学推导。取 $N=\langle a^\dagger a\rangle$、$g^{(2)}(0)=\langle a^\dagger a^\dagger aa\rangle/N^2$，并假设理想单模、无损耗；实际 halo 是多模的，有限带宽、时间分箱、损耗和噪声会稀释差异。Cheong 等在 Secs. I、V、VII 给出相同的 coherent/Gaussian/Fock 诊断框架。[1]

• 两点函数/功率谱的退化：若不同态被调到相同 $N_k$，且实验对随机相位平均，许多线性探测器只看到 $\langle\phi(x)\phi(y)\rangle$ 或等价 PSD；此时 pure coherent realization、phase-randomized mixture、Gaussian random field，甚至大 $N$ 的 Fock-like 模型可能给出相同或近似相同的二阶信号。差异转移到一点函数的相位稳定性、数涨落、$g^{(2)}$ 及更高阶连通相关。[1]
• 大 $N$ 不让所有差别消失：Fock 的 $g^{(2)}=1-1/N$ 的确趋近 coherent 的 1，但其 ${\rm Var}\,n=0$ 与 coherent 的 $N$ 仍不同；thermal 的 $g^{(2)}=2$ 也不随 $N$ 消失。实际可见性由探测效率、模数、积分窗口和加性量子噪声决定，而非仅由 $N$ 决定。[1]

5. 实验实际上“测到”什么？

• 线性场幅/相位：相位敏感读出在一个相干时间内可测窄带正交分量；长时间未知相位平均后，线性均值通常归零，实验转而累积功率。Marsh 指出典型 haloscope 输出可写成转换功率 $P=V\Gamma_E n$，其线宽 $\Delta\omega/\omega\sim v^2\sim10^{-6}$；长时间 frequency-domain 搜索往往对多种底层态模型不敏感。[2]
• 谐振腔功率或吸收/单量子激发率：轴子或暗光子场对腔模/量子比特提供近共振 drive，实验测输出功率、额外 photon counts、跃迁概率或事件率。这些量通常是耦合和场的二次量，直接约束谱密度，而不是暗物质密度矩阵的层析。[2][9][11]
• Fukuda 相关 transmon hidden-photon 方案：Chen、Fukuda 等把暗光子写成经典振荡场，动能混合产生 $E_X\propto\epsilon\sqrt{2\rho_{\rm DM}}$ 的 AC electric drive；它在共振时驱动 transmon，基本协议读出 $p_{g\to e}\simeq\sin^2(\eta t)$。因此被测量的是量子比特激发概率/重复事件数；在该协议中未知 drive 相位不会变成暗物质量子态层析。[9]
• Fukuda 的 QEC-like wave-DM noise mitigation：该工作假设一个在 de Broglie 波长内对多传感器共享、但未知相位的 classical oscillatory drive。理想的 $Y$ 读出在线性阶对 $\epsilon\tau\cos\phi$ 敏感；相位平均使它归零，标准无相位参考策略改测二阶的激发数。QEC-like 控制把 bit-flip/excitation noise 导出信号子空间，并在不同 $N$ 区域恢复相应灵敏度标度；这是噪声抑制和多传感器增益，不是对暗物质纯度或 nonclassicality 的测量。[10]
• 原子钟、磁共振和干涉仪：线性 portal 可使 $\alpha$、质量比或能级周期调制；例如 clock 读出 transition-frequency shift，CASPEr 型方案读出由 $\nabla\phi$ 产生的有效磁场和磁化。它们通常访问场的一点时间序列及其两点相关/PSD。[6][2]
• 二次耦合与“intensity interferometry”：二次 portal 对 $\phi^2$ 或 $(\nabla\phi)^2$ 敏感，含 DC 与 $2m$ 分量；传感器网络可交叉相关这类 stochastic intensity 以压低局域噪声。Masia-Roig 等的 Eq. (13) 是二次信号的一阶时空 coherence，并类比 Hanbury Brown–Twiss；它不是自动等于暗物质湮灭算符定义的量子光学 $g^{(2)}$，所以观测到该相关本身并不证明非经典量子态。[8]
• 跨传感器相关/相位：两站保留 Fourier 实部与虚部及同步相位，可测非对角协方差 $c_{ij}(\omega),s_{ij}(\omega)$，从而恢复单站功率谱丢失的速度方向信息；最敏感的 baseline 约为 de Broglie 尺度。该方案检验空间 coherence 和 $f(\mathbf v)$，但若分析仍采用 Gaussian random-field likelihood，它本身不判定 underlying density operator 是 pure 还是 mixed。[7]

6. 现有搜索能区分 pure coherent、phase-randomized mixture 或 Fock-like state 吗？

• 一般不能。标准功率、跃迁率或长时间 averaged PSD 主要固定 $N_k$ 与两点函数。pure coherent $|\alpha\rangle$ 和固定幅度的 phase-averaged coherent mixture 有完全相同的 Poisson number distribution；没有外部相位参考时，它们的区别藏在 number basis 的非对角元。Marsh 对 coherent/incoherent 模型也发现，频域平均下功率与线宽可以一致，可能的差别应在相干时间内的 time-domain statistics 中寻找。[2]
• Fock-like state 原则上可区分，现实中极难。理想的 sub-Poisson counting 或 $g^{(2)}(0)<1$ 是 nonclassical witness；但大 $N$ 时 Fock 的偏离只有 $1/N$，而暗物质—探测器耦合极弱、转换损耗和探测器自身量子噪声会主导。Cheong 等把测量 higher-order coherence 视为可能路径，同时明确指出所需噪声控制非常严苛。[1]
• photon counting 要区分“计数谁”。Agrawal 等已在 transmon-cavity 系统中制备腔模 $|n\rangle$，用 number-resolved QND 读出寻找 $|n\rangle\to|n+1\rangle$，并利用 $(n+1)$ 受激增强搜索暗光子。这里 Fock state 是人工制备的电磁腔态，counted photon 是转换后的腔光子；该实验提高对 classical DM drive 的灵敏度，但没有证明暗物质输入态是 Fock state。[11]
• 原则上可行的区分测量：(i) 有稳定外部相位参考的 homodyne/heterodyne，在单个相干 patch 内测 number-basis coherence 或 quadrature squeezing；(ii) 校准完整 transduction channel 后，对转换 quanta 做 number-resolved counting、HBT $g^{(2)}$ 与更高阶 factorial moments；(iii) 寻找 $g^{(2)}<1$、sub-Poisson、quadrature squeezing 或非正 $P$ 的 nonclassicality witness；(iv) 以 de Broglie 尺度分布、严格同步的多站相位读出检查空间 coherence。前三项才直接指向量子态差异；第四项主要重建 classical phase-space correlations，除非进一步测到不能由任意正 $P$ 模型解释的不等式违反。这里关于完整实验组合的表述是本任务的研究建议，而非任何单篇论文已经完成的方案。

7. 三层结论不能混为一谈

以下是本任务为避免过度宣称而采用的操作性分层。

1. 探测到 classical oscillatory signal：看到与 halo 线宽、调制、空间相关和耦合选择定则一致的窄带 drive。这是绝大多数现行搜索的直接目标。
2. 确认信号来自量子化的暗物质场：单个量子探测器的 click 或单个转换 photon 还不够，因为经典 drive 同样能使量子传感器离散跃迁。必须排除普通 classical stochastic source，并验证与量子化交换/统计有关、且不能由正 $P$ classical field 重现的观测。
3. 测量暗物质量子态或 nonclassicality：需要在已知 detector channel、损耗和噪声下重建多个矩或密度矩阵信息，并观察相位 coherence、squeezing、sub-Poisson 或其他 nonclassical witness。现有主流 wave-DM searches 尚未达到这一层；Cheong 等的量子密度矩阵框架是在说明应怎样提出这个问题，而不是报告已经完成的暗物质量子态层析。[1]

8. 最具体的下一步

若要把这个问题变成可计算方案，建议先选定一种耦合与一个可测的输出模，写出开放量子系统 channel $\rho_{\rm DM}\mapsto\rho_{\rm out}$，再比较五个输入模型在相同 $N_k$ 和相同 halo PSD 下给出的输出 $g^{(2)}$、factorial cumulants 与 quadrature covariance。首要 benchmark 是证明某个统计量在 pure coherent 与 phase-averaged coherent 间仍有差异；若没有外部相位参考，这两者在所有 phase-invariant number measurements 上应完全退化。随后再把多模、有限相干时间、效率、dark counts 和传感器 backaction 纳入 Fisher-information 或 likelihood 分析。这一段是本任务的推导与研究建议。

References