Dust spectral index and temperature based on GNILC

[zonca]

See comparison between GNILC and PySM 2 here: https://nbviewer.jupyter.org/gist/zonca/c6aa3e83f0151666e33a28cc84c0ba20

Suggestion by Hans Kristian (referring to the PySM 2 templates) was:

1) threshold beta_dust between 1.4 and 1.7, and smooth to ~2 degrees FWHM, to suppress the worst CO, point source and CIB effects 2) Smooth T_dust also to ~2 degrees, to reduce the worst Planck noise sensitivity

However no idea how to proceed. @brandonshensley @NicolettaK?

[NicolettaK]

Hi @zonca. I have the new templates for thermal dust beta and temperature that we are currently using for LiteBIRD. They are obtained by smoothing the Commander maps at 2°. No need to reduce the beta range as it is already ok after smoothing. Shall I share them with you?

[zonca]

no @NicolettaK thanks, I think it would be better to use the GNILC products for templates / spectral indices / dust temperature

[giuspugl]

I've finalized a procedure to inject smaller angular scales to the spectral parameters.

• we use the Planck GNILC maps of Bd and Td.
• fit a power law in the range $\ell \in [100,400]$ which corresponds to the range we fit the Intensity GNILC 2015 map
• for the whole analysis we evaluate fullsky power spectra with ANAFAST
• Modulate by total intensity map both Bd and Td (?)

[zonca]

@giuspugl the notebook for spectral index and T dust is at https://github.com/galsci/pysm/blob/f3b79810d22805cf5b0904eb13f65df6c6b6dbed/docs/preprocess-templates/gnilc_dust_spectral_index_T_dust.ipynb Do you have an executed version of that notebook for reference? (in a gist) I checked in this Gist (https://gist.github.com/giuspugl/d36a0e529a3b4a89c19bce5b58c9494a) but the cells related to spectral index and T_dust are not executed.

Instead of recomputing the modulation factor in that notebook, can I load the modulation map for T that I have created in the notebook which creates the templates? Or do they differ somehow?

[zonca]

I am referring to the modulation alms saved in cell 44 of https://nbviewer.org/gist/zonca/b0ec6cc14fa4c22ded497cc445a7ea46#Define-Modulation-maps

[giuspugl]

Sorry about that ! Here you can find an updated notebook https://gist.github.com/giuspugl/d36a0e529a3b4a89c19bce5b58c9494a i made sure it is compiled also for the spectral parameters too .

Instead of recomputing the modulation factor in that notebook, can I load the modulation map for T that I have created in the notebook which creates the templates? Or do they differ somehow?

Yes! the modulation for Td and Beta_dust is the same as the intensity one.

[zonca]

What do you think about the reusing the same modulation map?

On Thu, Jan 13, 2022, 09:24 Giuseppe Puglisi @.***> wrote:

Sorry about that ! Here you can find an updated notebook https://gist.github.com/giuspugl/d36a0e529a3b4a89c19bce5b58c9494a i made sure it is compiled also for the spectral parameters too .

Instead of recomputing the modulation factor in that notebook, can I load the modulation map for T that I have created in the notebook which creates the templates? Or do they differ somehow? Yes! the modulation for Td and Beta_dust is the same as the intensity one.

— Reply to this email directly, view it on GitHub https://github.com/galsci/pysm/issues/69#issuecomment-1012348349, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAC5Q4STE4BNPT7PPUH7FHDUV4DF5ANCNFSM4TVVVJUQ . You are receiving this because you were mentioned.Message ID: @.***>

[brandonshensley]

What do you think about the reusing the same modulation map?

Reusing the modulation map seems fine to me as long as the fluctuations themselves are different. Physically it is sensible for the amplitude of fluctuations to be related to the total intensity.

[zonca]

Implemented in https://github.com/galsci/pysm/pull/104