IR

Off[General::"spell1"]
Off[General::"spell"]

SetDirectory[NotebookDirectory[]];

Import data

hseis=Import["hseisTypeA.dat"];

hsusp=Import["hsusp20KSap.dat"];

hmirror=Import["hmirror20KSap.dat"];

QN=Import["QN780W23kgDCDRSE.dat"];

hgas=Import["gas1e-7.dat"];

Parameters

Gravity=6.673*10^(-11);(* Gravitational constant (m^3/Kg/s^2) *)
solarmass=1.989*10^(30);(* mass of Sun (Kg) *)
Mpc=365.2422*24*60*60*c*3.26*10^6;(* 1 Mpc (m) *)

mass1=1.4*solarmass;(* mass of a compact object (Kg) *)
mass2=1.4*solarmass;(* mass of other compace object (Kg) *)
mtotal=mass1+mass2;(* total moss of binary system (Kg) *)
mreduce=mass1*mass2/mtotal;(* reduced mass of binary system (Kg) *)
fmin=10;(* low cut off frequency of integral (Hz) *)
fmax=c^3/(6^(3/2)*Pi*Gravity*mtotal);(* high cut off frequency of integral (Hz) *)
SNR=8;(* threshold of Signal to Noise Ratio *)

fmax

hseis[[1000]][[1]]

hseis[[3196]][[1]]

Integration

integ=(∑_(ff=1000)^3196)(1/(QN[[ff]][[1]]^(7/3)*(QN[[ff]][[2]]^2+hsusp[[ff]][[2]]^2+hmirror[[ff]][[2]]^2+hseis[[ff]][[2]]^2+hgas[[ff]][[2]]^2))*(QN[[ff+1]][[1]]-QN[[ff]][[1]]));

Inspiral Range

Sky average (new); 0.442478 is given by Kanda-san.

(1/SNR)*2*0.442478*(Gravity/c^3)^(5/6)*c*(5/96*mreduce)^(1/2)*(mtotal/Pi^2)^(1/3)*Sqrt[4*integ]/Mpc