astropaint.profiles package

Submodules

astropaint.profiles.Battaglia16 module

library containing [projected] gas profiles

E. Schaan: Tau profile, from Battaglia 2016 watch typos in paper. This code is correct.

astropaint.profiles.Battaglia16.kSZ_T(R, R_200c, M_200c, v_r, redshift, *, T_cmb=2.7251)

kinetic Sunyaev Zeldovich effect #TODO: add reference

astropaint.profiles.Battaglia16.ne_2D(R, R_200c, M_200c, redshift)
2D physical electron number density profile in 1/(Mpc/h)^3

assuming fully ionized gas and primordial He abundance

Parameters

  • R (arraylike, float) – 2D projected distance from the center of the halo

  • R_200c (float [Mpc/h]) – comoving radius of the halo

  • M_200c (float [Msun/h]) – mass of the halo

  • redshift – redshift of the halo

astropaint.profiles.Battaglia16.ne_3D(r, R_200c, M_200c, redshift)

3D physical electron number density profile in 1/(Mpc/h)^3 assuming fully ionized gas and primordial He abundance

Parameters

  • r (arraylike, float) – 3D distance from the center of the halo

  • R_200c (float [Mpc/h]) – comoving radius of the halo

  • M_200c (float [Msun/h]) – mass of the halo

  • redshift – redshift of the halo

astropaint.profiles.Battaglia16.rho_gas_2D(R, R_200c, M_200c, redshift)

2D physical gas density profile in (Msun/h) / (Mpc/h)^3

Parameters

  • R (arraylike, float) – 2D projected distance from the center of the halo

  • R_200c (float [Mpc/h]) – comoving radius of the halo

  • M_200c (float [Msun/h]) – mass of the halo

  • redshift – redshift of the halo

astropaint.profiles.Battaglia16.rho_gas_2D_interp(R, R_200c, M_200c, redshift)

Interpolated version of rho_gas_2D 2D physical gas density profile in (Msun/h) / (Mpc/h)^3

Parameters

  • R (arraylike, float) – 2D projected distance from the center of the halo

  • R_200c (float [Mpc/h]) – comoving radius of the halo

  • M_200c (float [Msun/h]) – mass of the halo

  • redshift – redshift of the halo

astropaint.profiles.Battaglia16.rho_gas_3D(r, R_200c, M_200c, redshift)

3D physical gas density profile in (Msun/h) / (Mpc/h)^3

Parameters

  • r (arraylike, float) – 3D distance from the center of the halo

  • R_200c (float [Mpc/h]) – comoving radius of the halo

  • M_200c (float [Msun/h]) – mass of the halo

  • redshift – redshift of the halo

astropaint.profiles.Battaglia16.tau_2D(R, R_200c, M_200c, redshift)
Thompson scattering optical depth 2D projected profile in 1/(Mpc/h) comoving

ie you get the 2D tau profile by projecting this profile along the physical (not comoving) radial coordinate assuming fully ionized gas and primordial He abundance

Parameters

  • R (arraylike, float) – 2D projected distance from the center of the halo

  • R_200c (float [Mpc/h]) – comoving radius of the halo

  • M_200c (float [Msun/h]) – mass of the halo

  • redshift – redshift of the halo

astropaint.profiles.Battaglia16.tau_2D_interp(R, R_200c, M_200c, redshift)

Interpolated version of tau_2D: Thompson scattering optical depth 2D projected profile in 1/(Mpc/h) comoving ie you get the 2D tau profile by projecting this profile along the physical (not comoving) radial coordinate assuming fully ionized gas and primordial He abundance

Parameters

  • R (arraylike, float) – 2D projected distance from the center of the halo

  • R_200c (float [Mpc/h]) – comoving radius of the halo

  • M_200c (float [Msun/h]) – mass of the halo

  • redshift – redshift of the halo

astropaint.profiles.Battaglia16.tau_3D(r, R_200c, M_200c, redshift)

Thompson scattering optical depth 3D profile in 1/(Mpc/h) comoving ie you get the 2D tau profile by projecting this profile along the physical (not comoving) radial coordinate assuming fully ionized gas and primordial He abundance

Parameters

  • r (arraylike, float) – 3D distance from the center of the halo

  • R_200c (float [Mpc/h]) – comoving radius of the halo

  • M_200c (float [Msun/h]) – mass of the halo

  • redshift – redshift of the halo

astropaint.profiles.NFW module

library containing [projected] halo profiles

astropaint.profiles.NFW.BG(R_vec, c_200c, R_200c, M_200c, theta, phi, v_th, v_ph, *, T_cmb=2.7251)

Birkinshaw-Gull effect aka moving lens aka Rees-Sciama (moving gravitational potential)

astropaint.profiles.NFW.deflection_angle(R, c_200c, R_200c, M_200c, *, suppress=True, suppression_R=8)

calculate the deflection angle of a halo with NFW profile Using Eq 6 in Baxter et al 2015 (1412.7521)

Parameters

  • R – distance from the center of halo [Mpc]

  • c_200c – halo concentration parameter

  • R_200c – halo 200c radius in [Mpc]

  • M_200c – halo 200c mass of halo in M_sun

Returns

Return type

the deflection angle at distance R from the center of halo

astropaint.profiles.NFW.kSZ_T(R, rho_s, R_s, v_r, *, T_cmb=2.7251)

kinetic Sunyaev Zeldovich effect #TODO: add reference

astropaint.profiles.NFW.rho_2D(R, rho_s, R_s)

3D NFW profile intgrated along the line of sight

Returns

surface mass density

Return type

[M_sun/Mpc^2]

astropaint.profiles.NFW.rho_2D_bartlemann(R, rho_s, R_s)

projected NFW mass profile Eq. 7 in Bartlemann 1996: https://arxiv.org/abs/astro-ph/9602053

Returns

surface mass density

Return type

[M_sun/Mpc^2]

astropaint.profiles.NFW.rho_2D_interp(R, rho_s, R_s)

3D NFW profile intgrated along a sampled number of line of sights and then interpolated

Returns

surface mass density

Return type

[M_sun/Mpc^2]

astropaint.profiles.NFW.rho_3D(r, rho_s, r_s)

Calculate the 3D NFW density profile #TODO: add reference Eq.

Parameters

  • r – distance from the center

  • rho_s – density at radius r_s

  • r_s – characterisic radius R_200c/c_200c

Returns

Return type

rho = 4 * rho_s * r_s ** 3 / r / (r + r_s) ** 2

astropaint.profiles.NFW.tau_2D(R, rho_s, R_s)

projected NFW tau profile Eq. 7 in Battaglia 2016 :

Returns

tau

Return type

[NA]

astropaint.profiles.spherical module

library containing simple [projected] spherical halo profiles

astropaint.profiles.spherical.constant_density_2D(R, constant)

return a constant value at every input r :param R: distance from the center :type R: [Mpc] :param constant: multiplicative constant

Returns

Return type

constant

astropaint.profiles.spherical.linear_density_2D(R, intercept, slope)

return a R*constant at every input R :param R: distance from the center :type R: [Mpc] :param intercept: intercept of the line :param slope: slope of the line

Returns

Return type

intercept

astropaint.profiles.spherical.solid_sphere_2D(R, M_200c, R_200c)

projected mass density of uniform sphere :param r: distance from the center :type r: [Mpc] :param M_200c: total mass of the sphere :type M_200c: [M_sun] :param R_200c: total radius (edge) of the sphere :type R_200c: [Mpc]

Returns

Return type

Sigma = M_200c /2/pi * sqrt(R_tot**2 - r**2)/R_tot**3

Module contents