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.17 [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.15 [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.504012 0.546657 1396.379712 0.490958 0.510632 0.097520 0.087238 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.416777 3845.441877 0.366849 0.409885 0.115054 0.105756 0.029716 0.002919 0.180990 0.256332 0.262478
100000.000000 12.000000 131.875959 24496.871836 0.000374 0.025316 0.102252 0.091440 4765.788970 0.085224 0.088805 0.013267 0.012024 0.005070 0.000300 0.023299 0.033676 0.040693
5000000.000000 12.000000 65.171609 5995.355414 0.019861 0.038163 0.759741 0.900386 1179811.156712 0.738202 0.731565 0.017127 0.039468 0.051369 0.001989 0.147183 0.221789 0.044944

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 4.71195758e-45 4.71195758e-45 6.70580303e-55] [2.10000000e-05 4.71334803e-45 4.71334803e-45 8.15094275e-55] [2.20000000e-05 4.71480631e-45 4.71480631e-45 9.81798275e-55] ... [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