The python codes contained in this folder have been used to reproduce the main the results shown in the paper Testing the Quasar Hubble Diagram with LISA Standard Sirens [1]. We implemented a Bayesian hypothesis test to understand how many LISA MBHBs Standard Siren observations are necessary to test the deviation from the CDM model found by Risaliti and Lusso with quasars observations [2].
In addition, we also implemented an MCMC to understand the constraining power of LISA standard sirens on matter density and Hubble parameter for a flat
CDM model.
This folder contains:
-
mcmc_flatLCDM.py produces the posterior samples matter density
and Hubble parameter h of a flat
CDM model given redshift, luminosity distance and luminosity distance uncertainty observations (for instance the dataset SSdatasets/median_dataset_15SS.dat)
execute with: python mcmc_flatLCDM.py SSdatasets/median_dataset_15SS.dat the posterior samples are stored in the folder SSdatasets with name median_dataset_15SS.dat.h5
-
evidence_generation.py produces the samples of the Bayes factor realizations for each fixed number of standard sirens
execute: python evidence_generation.py output: EvidenceData_BayesPlot.h5 where each column contains the samples of the Bayes factor realizations, and the column number corresponds to the number of standard sirens
-
LISA_vs_Quasars.nb Mathematica notebook used to construct catalogues of LISA MBHB Standard Sirens and to investigate the constraints on Hubble parameter and matter density
-
SSdatasets Folder containing the representative LISA MBHB Standard Siren observations and results obtained with LISA_vs_Quasars.nb
-
SS_catalogues.dat (Redshift, Luminosity Distance [Gpc], uncertainty on Luminosity Distance [Gpc]) observations of the catalogues used to construct the Bayes Factor distributions
Lorenzo Speri: [email protected]
Our codes used in Testing the Quasar Hubble Diagram with LISA Standard Sirens [1] to analyze quasar observations can be also provided, however their usage depends on the quasar data which are available from the corresponding author of [2] upon request. Please contact the author of [2] for the usage permission of their data and we will provide all the datasets and the codes to reproduce our analysis.
[1]: L. Speri, N. Tamanini, R. R. Caldwell, J. Gair, B. Wang. Testing the Quasar Hubble Diagram with LISA Standard Sirens https://journals.aps.org/prd/abstract/10.1103/PhysRevD.103.083526
[2]: G. Risaliti, E. Lusso. Cosmological constraints from the Hubble diagram of quasars at high redshifts. Nat Astron 3, 272–277 (2019). https://doi.org/10.1038/s41550-018-0657-z