Continuous Wave Analysis of the NANOGrav 15-Year Dataset

This repository contains the data used in the analysis detailed within "<em>The NANOGrav 15-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries</em>" (DOI arXiv:2306.16222). It also provides instruction on how one could reproduce those results on their own. <code>jar</code> folder: contains pickled <code>enterprise</code> pulsar objects that contain each pulsar's TOA data and timing model information <code>v1p1_all_dict.json</code>: dictionary with measured white noise and red noise parameters for each pulsar - used to fix white noise parameters when running QuickCW <code>emp_dist_15yr_v1p1_bence_my_run_v4.3.pkl</code>: <code>enterprise_extension</code> <code>EmpiricalDistribution2D</code> objects for the amplitude and spectral index of each pulsar's red noise - used to make informed red noise proposals when running QuickCW <code>pulsar_distances_15yr.pkl</code>: dictionary containing 3 element list for each pulsar specifying its distance [kpc], statistical error of distance [kpc], and method of distance measurement [PX for parallax or DM for dispersion measure] <code>15yr_quickCW_detection.h5</code>: HDF5 file containing results from QuickCW detection runs (log uniform amplitude prior). Includes the following data: MCMC samples (<code>samples_cold</code>) log likelihood values corresponding to those samples (<code>log_likelihood</code> ) temperature ladder used for the parallel tempering (<code>Ts</code>) names of parameters corresponding to each column in the samples array (<code>par_names</code>) rate of different MCMC jumps being accepted (<code>acc_fraction</code>) diagonal fisher matrix used for some proposals (<code>fisher_diag</code>) <code>15yr_quickCW_UL.h5</code>: Same as <code>15yr_quickCW_detection.h5</code> but for upper limits based on runs with uniform amplitude prior. <code>15yr_cw_3d_limits_v4.npz</code>: npz file containing limits as a function of frequency and sky location. Specifically it has the following arrays: <code>F_edges</code>: array specifying the edges of the frequency bins <code>strain_limit_skies</code>: 2D array containing the strain upper limits as a function of frequency bin and sky pixel <code>strain_limit_sky_sigmas</code>: 2D array containing 1-sigma statistical errors on the strain upper limits <code>dist_limit_skies</code>: 2D array containing the luminosity distance limits in Mpc as a function of frequency bin and sky pixel <code>dist_limit_sky_sigmas</code>: 2D array containing 1-sigma statistical errors on the distance limits in Mpc本代码仓库包含了《NANOGrav 15年数据集：单个超大质量双黑洞并合产生的引力波的贝叶斯上限》（DOI arXiv:2306.16222）一文中详细分析所使用的全部数据，同时提供了复现该研究结果的操作指南。 `jar` 文件夹：内含经序列化（pickle）处理的`enterprise`脉冲星对象，其中存储了每颗脉冲星的到达时间（Time of Arrival, TOA）数据与计时模型信息。 `v1p1_all_dict.json`：存储每颗脉冲星实测白噪声与红噪声参数的字典，用于在运行QuickCW时固定白噪声参数。 `emp_dist_15yr_v1p1_bence_my_run_v4.3.pkl`：`enterprise_extension`库的`EmpiricalDistribution2D`二维经验分布对象，对应每颗脉冲星红噪声的振幅与谱指数，用于在运行QuickCW时生成合理的红噪声提议分布。 `pulsar_distances_15yr.pkl`：存储每颗脉冲星相关信息的字典，其中包含一个三元素列表，依次为脉冲星距离（单位：千秒差距[kpc]）、距离的统计误差（单位：千秒差距[kpc]）以及距离测量方法——`PX`代表视差测量，`DM`代表色散量测量。 `15yr_quickCW_detection.h5`：存储QuickCW探测运行结果的HDF5格式文件（采用对数均匀振幅先验），包含以下数据：马尔可夫链蒙特卡洛（Markov Chain Monte Carlo, MCMC）采样样本（`samples_cold`）、对应采样的对数似然值（`log_likelihood`）、并行回火所用的温度阶梯（`Ts`）、采样数组各列对应的参数名称（`par_names`）、不同MCMC跳转的接受率（`acc_fraction`）、部分提议所用的对角费舍尔矩阵（`fisher_diag`）。 `15yr_quickCW_UL.h5`：与`15yr_quickCW_detection.h5`结构一致，但对应采用均匀振幅先验的运行所得到的上限结果。 `15yr_cw_3d_limits_v4.npz`：存储随频率与天区位置变化的上限值的npz格式文件，具体包含以下数组：`F_edges`：频率分箱的边界数组；`strain_limit_skies`：二维数组，存储各频率分箱与天区像素对应的应变上限；`strain_limit_sky_sigmas`：二维数组，存储应变上限的1σ统计误差；`dist_limit_skies`：二维数组，存储各频率分箱与天区像素对应的光度距离上限（单位：百万秒差距[Mpc]）；`dist_limit_sky_sigmas`：二维数组，存储光度距离上限的1σ统计误差（单位：百万秒差距[Mpc]）