Contain Mathematica script to calculate amplification factor
pyCBC_function.py
basic functions to calculate waveforms and lensingLISA.py
LISA sensitivities from here
Codes used in the paper arXiv:2104.07055
MSD_2.py
general code to compute lensingMSD_matrix_2.py
code to plot matrix of SNRsSNR_regimes.py
code to compute SNR calculation and plotSNR_regimes_2.py
edited version of "SNR_regimes.py"
multi_lenses.py
code to compute lensing of sigle signal by multiple lenses
SIS_NFW.py
SNR.py
calculate SNRlens_models.py
lens_profiles.py
compute and plot lens mass profilesparam_lens_model.py
contain data about lens modelsphase.py
compute phase calculations and plotsplots.py
compute different kind of plots
amps files
..._FP means for phase and with the normalization factor phi_m
..._norm means long array and with the normalization factor phi_m
..._FP_nn means for phase and WITHOUT normalization
*... * means long array and WITHOUT normalization
m1 = m2 = 5*10^7
zS = 1
dt = 1
N = 2157032
N_mt = 3157032
zL = 0.5
theta_E_SIS = 2.62394 * 10^-7
lens | sigma [km/s] | r @ y=1 : theta_E*DL [m] |
---|---|---|
SIS | 66.256 | 9.62299*10^18 |
NIS | 66.10 | 9.62299*10^18 |
lm | rho_s [kg/m^3] | rs [m] |
NFW | 7.954495462000312*10^-22 | 9.065935*10^20 |
NFW_2 | 7.954495*10^-22 | 3.4313*10^20 |
Hernquist | 2.87*10^-21 | 3.5948*10^20 |
GNFW 0 | 1.10*10^-21 | 3.3531*10^20 |
GNFW 2 | 4.938784307737533*10^-22 | 1.517823*10^20 |
BUR | 7.323121259691678*10^-22 | 3.43436*10^20 |
- ySIS = 1.00000 (red)
- yNIS = 1.00000
- yNFW = 0.01061 (blue)
- yNF2 = 0.02804 (* df_clr[9] , orangered (zs=0.55) *)
- yHER = 0.02677
- yGN0 = 0.02870
- yGN2 = 0.06340 (green)
- yBUR = 0.02802
z_L = 0.25
- ySIS = 1.00000 (red)
- yNFW = 0.010830 (blue) (0.00680)
- yNF2 = 0.028615 (* df_clr[9] , orangered (zs=0.55) *) (0.01796)
- yGN2 = 0.064674 (green) (0.04059)
z_L = 0.15
lens | sigma [km/s] | r @ y=1 : theta_E*DL [m] |
---|---|---|
SIS | 69.750 | 8.685542809103445*10^18 |
lm | rho_s [kg/m^3] | rs [m] |
NFW | 5.340197838131539*10^-22 | 1.474791032921173*10^21 |
NFW_2 | 5.340197838131539*10^-22 | 4.4068076149185295*10^20 |
GNFW 2 | 3.315620133481915*10^-22 | 1.930856447560086*10^20 |
- ySIS = 1.00000 (red)
- yNFW = 0.005889 (blue)
- yNF2 = 0.019709 (* df_clr[9] , orangered (zs=0.55) *)
- yGN2 = 0.044983 (green)
- ySIS = 3.14100 (red)
- yNFW = 0.03340 (blue)
- yNF2 = 0.08809 (* df_clr[9] , orangered (zs=0.55) *)
- yGN2 = 0.19914 (green)
- ySIS = 0.10000
- yNIS = 0.10000
- yNFW = 0.00106
- yNF2 = 0.00280
- yHER = 0.00268
- yGN0 = 0.00287
- yGN2 = 0.00634
- yBUR = 0.00280
- ySIS - 2.00000
- yNIS - 2.00000
- yNFW - 0.02123
- yNF2 - 0.05609
- yHER - 0.05354
- yGN0 - 0.05740
- yGN2 - 0.12680
- yBUR - 0.05604
zL = 0.5
theta_E_SIS = 1.86 * 10^-7
lens | sigma [km/s] | r @ y=1 : theta_E*DL [m] |
---|---|---|
SIS | 55.71 | 6.80448*10^18 |
NIS | 55.58 | 6.80448*10^18 |
(lm) | rho_s [kg/m^3] | rs [m] |
NFW | 7.954495462000312*10^-22 | 8.39304*10^20 |
NFW_2 | 7.954495*10^-22 | 2.88202*10^20 |
Hernquist | 2.87*10^-21 | 3.32636*10^20 |
GNFW 0 | 1.10*10^-21 | 2.82031*10^20 |
GNFW 2 | 4.938784307737533*10^-22 | 1.2713*10^20 |
BUR | 7.323121259691678*10^-22 | 2.8819*10^20 |
- ySIS - 1.00000
- yNIS - 1.00000
- yNFW - 0.00811
- yNF2 - 0.02361
- yHER - 0.02046
- yGN0 - 0.02413
- yGN2 - 0.05352
- yBUR - 0.02356
- ySIS = 0.10000
- yNIS = 0.10000
- yNFW = 0.00081
- yNF2 = 0.00236
- yHER = 0.00205
- yGN0 = 0.00241
- yGN2 = 0.00535
- yBUR = 0.00236