使用 CAMB 和 Healpy 生成 E 和 B 模式地图答案

【问题标题】：Producing E and B mode maps using CAMB and Healpy使用 CAMB 和 Healpy 生成 E 和 B 模式地图
【发布时间】：2016-03-12 20:04:17
【问题描述】：

我一直在尝试使用 CAMB 和 healpy 制作 E 和 B 模式地图，但是到目前为止我制作的东西看起来并不多。

我正在使用 CAMB 生成模拟地图，链接如下： http://lambda.gsfc.nasa.gov/toolbox/tb_camb_form.cfm, 并且我保留了所有默认参数，而不是选中张量 cl 和不同 r 的框。

我一直在使用以下代码来制作地图：

# read in the power spectra produced by camb, ordering is: [TT, EE, BB, TE]
cl = hp.read_cl('scalar_and_tensor_power_spectra_produced_by_CAMB.fits')

# set the random seed for phase information
numpy.random.seed(some number) #Im using 1 for now

#Produce E mode and B mode maps
EEmap = hp.synfast(cl[1], NSIDE) 
BBmap = hp.synfast(cl[2], NSIDE)

快速谷歌搜索显示与我得到的结果截然不同的结果，您可以在以下链接中看到： http://tinypic.com/view.php?pic=ifvrpt&s=9#.VmY7Jx8So8o http://tinypic.com/view.php?pic=a5ngx&s=9#.VmY7KR8So8o

我希望在 E 和 B 中看到某些图案，因为它们是渐变和卷曲成分，这在我在网上找到的图像中很明显，但在我的产品中却没有。

所以我有三个问题：

我用于生成地图的逻辑是否有错误？
是否需要额外的过滤才能看到预期的 E 和 B 模式？
- 类似于以下链接： http://cosmologist.info/polar/EBsupport.html
在 healpy 以外的 healpix 版本中，mollview 中有一个参数称为 POLARIZATION (http://healpix.jpl.nasa.gov/html/idlnode33.htm#idl:mollview:polarization) 其中设置 POLARIZATION = 3 允许您将偏振信息重叠绘制为无头矢量。这就是我认为他们在这里所做的事情： http://bicepkeck.org/media/b_over_b_rect_BICEP2.png &这里的第一张图片： https://inspirehep.net/record/787860/plots, 有没有办法在healpy中做到这一点？由于 POLARIZATION 参数似乎不包含在此版本中。

【问题讨论】：

标签： healpy

【解决方案1】：

请参阅synalm http://healpy.readthedocs.org/en/1.x/generated/healpy.sphtfunc.synalm.html#healpy.sphtfunc.synalm 的帮助，如果您使用hp.synfast(cl) 来创建温度图，其中cl 的温度。您应该为每个组件创建一个 Cl 列表，即 TT、EE、BB、TE，然后将所有 Cl 设为零且只有一个非零。
您可以使用假设为 1 度的高斯光束平滑输出映射，以模拟真实仪器如何看待该信号。
不支持在healpy中绘制vectos

【讨论】：

我认为这可能源于对 E 和 B 模式如何工作的知识不足。如果 synfast(cl) 生成的 IQU 映射中只有 B 的贡献（因为其他所有内容都设置为零），那么说 B = sqrt(Q^2 + U^2) 是否正确？或类似的东西？
不，B在另一个基数上，就是当你进行谐波变换时，它总是用于功率谱域。 sqrt(Q^2 + U^2)是P，即总极化幅度。
当我做我认为你说的事情时：synfast([zeros, zeros, BB cl's, zeros])，我得到一个 IQU 图（我认为），其强度分量为零且 Q 非零和 U 组件。我该如何解释这个？这种方法不是制作E图和B图的方法吗？或者我只是对什么是 E 和 B 的“地图”感到困惑？
正确，你得到一个 IQU 图，Q 和 U 中的极化信号是你的 B 模式的来源。 E 和 B 在角功率谱域中，Q 和 U 在图中。查看background.uchicago.edu/~whu/polar/webversion/polar.html 或一本关于宇宙学的书籍（Liddle 或 Partridge 或 Durer）
我认为这对我来说现在（相对而言）都有意义。感谢您的帮助。