Disentangling seed formation mechanisms [4.5 yr]

FoM

Science investigation and Observational requirement: SI 3.2.1 from Science requirement document.

Detection of GW inspiral with source-frame mass in the range \(10^3-10^5\) M\(_\odot\) and formation redshift in the range \(z = 10-15\) with SNR\(>10\)

- OK 15.43 [10.0 15.0 30.0] : SNR

System of mass 5e3 at z=12

- OK 23.48 [10.0 15.0 30.0] : SNR

System of mass 1e4 at z=12

- GREAT 131.88 [10.0 15.0 30.0] : SNR

System of mass 1e5 at z=12

- GREAT 65.06 [10.0 15.0 30.0] : SNR

System of mass 5e6 at z=12

- OK 0.5 [1.0 0.5 0.2] : Relative error on the distance luminosity

System of mass 5e3 at z=12

- GREAT 0.1 [0.75 0.3 0.1] : Relative error on the detector frame mass 1

System of mass 5e3 at z=12

- GREAT 0.09 [0.75 0.3 0.1] : Relative error on the detector frame mass 2

System of mass 5e3 at z=12

- OK 0.12 [0.75 0.3 0.1] : Relative error on the detector frame mass 1

System of mass 1e4 at z=12

- OK 0.11 [0.75 0.3 0.1] : Relative error on the detector frame mass 2

System of mass 1e4 at z=12

- GREAT 0.01 [0.75 0.3 0.1] : Relative error on the detector frame mass 1

System of mass 1e5 at z=12

- GREAT 0.01 [0.75 0.3 0.1] : Relative error on the detector frame mass 2

System of mass 1e5 at z=12

- GREAT 0.02 [0.75 0.3 0.1] : Relative error on the detector frame mass 1

System of mass 5e6 at z=12

- GREAT 0.04 [0.75 0.3 0.1] : Relative error on the detector frame mass 2

System of mass 5e6 at z=12

- OK 0.31 [0.7 0.4 0.3] : Asbolute error on spin 1

System of mass 5e3 at z=12

- GREAT 0.18 [0.7 0.4 0.3] : Asbolute error on spin 1

System of mass 1e4 at z=12

- GREAT 0.02 [0.7 0.4 0.3] : Asbolute error on spin 1

System of mass 1e5 at z=12

- GREAT 0.17 [0.7 0.4 0.3] : Asbolute error on spin 1

System of mass 5e6 at z=12

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
5000.000000 12.000000 15.434610 54742.791039 0.000213 0.421103 0.504011 0.536574 1396.376758 0.489936 0.513773 0.097520 0.087239 0.029493 0.002915 0.305421 0.564243 0.415041
10000.000000 12.000000 23.475441 45302.684104 0.000224 0.266935 0.417973 0.411384 3829.582222 0.363954 0.410412 0.115058 0.105768 0.029719 0.002909 0.181002 0.256330 0.262478
100000.000000 12.000000 131.875892 24496.871836 0.000386 0.025317 0.102627 0.090674 4784.060689 0.083288 0.089497 0.013256 0.011995 0.005092 0.000300 0.023358 0.033663 0.040747
5000000.000000 12.000000 65.062897 5706.570176 0.021115 0.038724 0.818693 1.016945 1243844.959322 0.783217 0.776696 0.017988 0.041805 0.058897 0.002078 0.165578 0.231216 0.046945

MBHB with waveform IMRPhenomHM of

  • mass = [5000.0, 10000.0, 100000.0, 5000000.0] solar mass
  • mass ratio in [1, 3] (log uniform),
  • redshift = [12]
  • 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: 10.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=5000.0, z=12

Ms=10000.0, z=12

Ms=100000.0, z=12

Ms=5000000.0, z=12

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 7.71748532e-41 7.71748532e-41 8.45230844e-54] [2.10000000e-05 5.48584269e-41 5.48584269e-41 7.52449171e-54] [2.20000000e-05 3.96262828e-41 3.96262828e-41 6.82098233e-54] ... [9.99998000e-01 2.73617388e-36 2.73617388e-36 5.73895037e-36] [9.99999000e-01 2.73619429e-36 2.73619429e-36 5.73933494e-36] [1.00000000e+00 2.73621460e-36 2.73621460e-36 5.73971931e-36]]
  • timetomerger_max : 4.5
  • modes : [(2, 2), (2, 1), (3, 3), (4, 4)] '''

Input parameters


- log : log/SO2_parameter_estimation.seeds.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.seeds.md
- snr_thresh : 10
- fom_range : []
- fom_parameter : []
- fom_title : Disentangling seed formation mechanisms
- fom_or : ['SI 3.2.1']
- fom_objective : Detectability/PE for MBHB
- fom_config : SO2/parameter_estimation.seeds.yml
- output : SciRDv1_redbook/data_SO2a_parameter_estimation.seeds.pkl
- mass_range : [5000.0, 10000.0, 100000.0, 5000000.0]
- z_range : [12]
- spin_range : [0, 0]
- mass_ratio_range : [1, 3]
- fom_obj : Detection of GW inspiral with source-frame mass in the range $10^3-10^5$ M$_\odot$ and formation redshift in the range $z = 10-15$ with SNR$>10$
- foms : [{'system': 0, 'value': 'snr', 'extra': 30.0, 'good': 15.0, 'bad': 10.0, 'title': 'SNR', 'descr': 'System of mass 5e3 at z=12', 'inverse': False}, {'system': 1, 'value': 'snr', 'extra': 30.0, 'good': 15.0, 'bad': 10.0, 'title': 'SNR', 'descr': 'System of mass 1e4 at z=12', 'inverse': False}, {'system': 2, 'value': 'snr', 'extra': 30.0, 'good': 15.0, 'bad': 10.0, 'title': 'SNR', 'descr': 'System of mass 1e5 at z=12', 'inverse': False}, {'system': 3, 'value': 'snr', 'extra': 30.0, 'good': 15.0, 'bad': 10.0, 'title': 'SNR', 'descr': 'System of mass 5e6 at z=12', 'inverse': False}, {'system': 0, 'value': 'dDL/DL', 'extra': 0.2, 'good': 0.5, 'bad': 1.0, 'title': 'Relative error on the distance luminosity', 'descr': 'System of mass 5e3 at z=12'}, {'system': 0, 'value': 'dm1/m1', 'extra': 0.1, 'good': 0.3, 'bad': 0.75, 'title': 'Relative error on the detector frame mass 1', 'descr': 'System of mass 5e3 at z=12'}, {'system': 0, 'value': 'dm2/m2', 'extra': 0.1, 'good': 0.3, 'bad': 0.75, 'title': 'Relative error on the detector frame mass 2', 'descr': 'System of mass 5e3 at z=12'}, {'system': 1, 'value': 'dm1/m1', 'extra': 0.1, 'good': 0.3, 'bad': 0.75, 'title': 'Relative error on the detector frame mass 1', 'descr': 'System of mass 1e4 at z=12'}, {'system': 1, 'value': 'dm2/m2', 'extra': 0.1, 'good': 0.3, 'bad': 0.75, 'title': 'Relative error on the detector frame mass 2', 'descr': 'System of mass 1e4 at z=12'}, {'system': 2, 'value': 'dm1/m1', 'extra': 0.1, 'good': 0.3, 'bad': 0.75, 'title': 'Relative error on the detector frame mass 1', 'descr': 'System of mass 1e5 at z=12'}, {'system': 2, 'value': 'dm2/m2', 'extra': 0.1, 'good': 0.3, 'bad': 0.75, 'title': 'Relative error on the detector frame mass 2', 'descr': 'System of mass 1e5 at z=12'}, {'system': 3, 'value': 'dm1/m1', 'extra': 0.1, 'good': 0.3, 'bad': 0.75, 'title': 'Relative error on the detector frame mass 1', 'descr': 'System of mass 5e6 at z=12'}, {'system': 3, 'value': 'dm2/m2', 'extra': 0.1, 'good': 0.3, 'bad': 0.75, 'title': 'Relative error on the detector frame mass 2', 'descr': 'System of mass 5e6 at z=12'}, {'system': 0, 'value': 'dchi1', 'extra': 0.3, 'good': 0.4, 'bad': 0.7, 'title': 'Asbolute error on spin 1', 'descr': 'System of mass 5e3 at z=12'}, {'system': 1, 'value': 'dchi1', 'extra': 0.3, 'good': 0.4, 'bad': 0.7, 'title': 'Asbolute error on spin 1', 'descr': 'System of mass 1e4 at z=12'}, {'system': 2, 'value': 'dchi1', 'extra': 0.3, 'good': 0.4, 'bad': 0.7, 'title': 'Asbolute error on spin 1', 'descr': 'System of mass 1e5 at z=12'}, {'system': 3, 'value': 'dchi1', 'extra': 0.3, 'good': 0.4, 'bad': 0.7, 'title': 'Asbolute error on spin 1', 'descr': 'System of mass 5e6 at z=12'}]

Download

PE data seeds_m_5000.0_z_12

Full data set seeds_m_5000.0_z_12

PE data seeds_m_10000.0_z_12

Full data set seeds_m_10000.0_z_12

PE data seeds_m_100000.0_z_12

Full data set seeds_m_100000.0_z_12

PE data seeds_m_5000000.0_z_12

Full data set seeds_m_5000000.0_z_12