Discovering unequal mass IMBHs [4.5 yr]

FoM

Science investigation and Observational requirement: SI 3.2.4b from Science requirement document.

Detection of MBHB with source-frame mass in the range \(10^4-10^5\) M\(_\odot\) and formation redshift in the range $z=1 $ with SNR\(>20\)

Measuring ['chirp mass'] with a given relative precision in \%.

- GREAT 0.0 [20 10 5] : Ms = 20000.0 ; z = 1.0: chirp mass

- GREAT 0.0 [20 10 5] : Ms = 200000.0 ; z = 1.0: chirp mass

- OK 148.68 [20 100 150] : Ms = 20000.0 ; z = 1: SNR

- GREAT 834.33 [20 100 150] : Ms = 200000.0 ; z = 1: SNR

- OK 0.01 [0.3 0.05 0.005] : Ms = 20000.0 ; z = 1: mass ratio rel. error

- GREAT 0.0 [0.3 0.05 0.005] : Ms = 200000.0 ; z = 1: mass ratio rel. error

- GREAT 0.01 [0.2 0.05 0.01] : Ms = 20000.0 ; z = 1: spin 2 error

- GREAT 0.01 [0.2 0.05 0.01] : Ms = 200000.0 ; z = 1: spin 2 error

- GREAT 0.0 [0.2 0.05 0.01] : Ms = 20000.0 ; z = 1: chirp mass

- GREAT 0.0 [0.2 0.05 0.01] : Ms = 200000.0 ; z = 1: chirp mass

Summary

Ms z snr tdet_s dMc/Mc dq/q dDL/DL dz/z dM dm1s/m1s dm2s/m2s dm1/m1 dm2/m2 dchieff dMf/Mf dchi1 dchi2 diota
20000.000000 1.000000 148.676383 2646988.702886 0.000015 0.005755 0.030317 0.024486 9.798005 0.012781 0.012809 0.000291 0.000213 0.000360 0.000005 0.000771 0.009097 0.020577
200000.000000 1.000000 834.325955 3583299.748314 0.000013 0.001634 0.010383 0.008381 307.635971 0.004437 0.004436 0.000909 0.000635 0.000692 0.000010 0.001269 0.005177 0.003367

MBHB with waveform IMRPhenomHM of

  • mass = [20000.0, 200000.0] solar mass
  • mass ratio in [10, 10] (log uniform),
  • redshift = [1.0]
  • t0 in [0.10 - 0.95] (uniform)
  • chi in [0.00 - 0.00] (uniform))

Conducting monte-carlo over 1000 trials, randomizing spins, sky position and orientation

Assuming SNR threshold for detection: 20.0

For derived parameters (dz/z, dMf/Mf), random samples are drawn from the inverse Fisher, derived quantities are computed and we extract their standard deviation and mean.

For the following waveform parameters (dMc/Mc, dDL/DL, dq/q), we use a log10 parameterization.

We report observer frame errors.

Detailed results

Ms=20000.0, z=1.0

Ms=200000.0, z=1.0

waveform_params

'''

  • minf : 2e-05
  • maxf : 1
  • fend : None
  • tmin : 0.0
  • tmax : None
  • phiref : 0.0
  • fref_for_phiref : 0.0
  • toffset : 0.0
  • TDI : TDI2AET
  • acc : 0.0001
  • order_fresnel_stencil : 0
  • approximant : IMRPhenomHM
  • responseapprox : full
  • frozenLISA : False
  • TDIrescaled : False
  • LISAconst : {'OrbitOmega': 1.9909865927683783e-07, 'OrbitPhi0': 0.0, 'Orbitt0': 0.0, 'OrbitR': 149597870700.0, 'OrbitL': 2500000000.0}
  • LISAnoise : [[2.00000000e-05 3.74868279e-45 3.74868279e-45 1.45384261e-55] [2.10000000e-05 3.51079457e-45 3.51079457e-45 1.76691006e-55] [2.20000000e-05 3.31246021e-45 3.31246021e-45 2.12793518e-55] ... [9.99998000e-01 1.35664204e-36 1.35664204e-36 2.82767611e-36] [9.99999000e-01 1.35666152e-36 1.35666152e-36 2.82788487e-36] [1.00000000e+00 1.35668095e-36 1.35668095e-36 2.82809354e-36]]
  • timetomerger_max : 4.5
  • modes : [(2, 2), (2, 1), (3, 3), (4, 4)] '''

Input parameters


- log : log/SO2_parameter_estimation.imbh2.log
- config : config.yml
- subconfig : 
- duration : 4.5 yr
- duty_cycle : 0.82
- armlength : 2500000000.0 m
- noise : SciRDv1
- show : False
- fast : 0
- tdi2 : 1
- freqs : [2e-05, 1]
- force_duration : 4.5_yr
- report : SciRDv1_redbook/report_SO2a_parameter_estimation.imbh2.md
- snr_thresh : 20
- fom_range : [20, 10, 5]
- fom_parameter : ['dMc/Mc']
- fom_title : Discovering unequal mass IMBHs
- fom_or : ['SI 3.2.4b']
- fom_objective : Detectability/PE for MBHB
- fom_config : SO2/parameter_estimation.imbh2.yml
- output : SciRDv1_redbook/data_SO2a_parameter_estimation.imbh2.pkl
- mass_range : [20000.0, 200000.0]
- z_range : [1.0]
- spin_range : [0, 0]
- mass_ratio_range : [10, 10]
- fom_obj : Detection of MBHB with source-frame mass in the range $10^4-10^5$ M$_\odot$ and formation redshift in the range $z=1 $ with SNR$>20$
- foms : [{'system': 0, 'value': 'snr', 'extra': 150, 'good': 100, 'bad': 20, 'title': 'Ms = 20000.0 ; z = 1: SNR', 'inverse': False, 'descr': None}, {'system': 1, 'value': 'snr', 'extra': 150, 'good': 100, 'bad': 20, 'title': 'Ms = 200000.0 ; z = 1: SNR', 'inverse': False, 'descr': None}, {'system': 0, 'value': 'dq/q', 'extra': 0.005, 'good': 0.05, 'bad': 0.3, 'title': 'Ms = 20000.0 ; z = 1: mass ratio rel. error', 'descr': None}, {'system': 1, 'value': 'dq/q', 'extra': 0.005, 'good': 0.05, 'bad': 0.3, 'title': 'Ms = 200000.0 ; z = 1: mass ratio rel. error', 'descr': None}, {'system': 0, 'value': 'dchi2', 'extra': 0.01, 'good': 0.05, 'bad': 0.2, 'title': 'Ms = 20000.0 ; z = 1: spin 2 error', 'descr': None}, {'system': 1, 'value': 'dchi2', 'extra': 0.01, 'good': 0.05, 'bad': 0.2, 'title': 'Ms = 200000.0 ; z = 1: spin 2 error', 'descr': None}, {'system': 0, 'value': 'dMc/Mc', 'extra': 0.01, 'good': 0.05, 'bad': 0.2, 'title': 'Ms = 20000.0 ; z = 1: chirp mass', 'descr': None}, {'system': 1, 'value': 'dMc/Mc', 'extra': 0.01, 'good': 0.05, 'bad': 0.2, 'title': 'Ms = 200000.0 ; z = 1: chirp mass', 'descr': None}]

Download

PE data imbh2_m_20000.0_z_1.0

Full data set imbh2_m_20000.0_z_1.0

PE data imbh2_m_200000.0_z_1.0

Full data set imbh2_m_200000.0_z_1.0