README:

These data files contain the results of simulations of BBH mergers matching the first 3 GW events (BBH mergers) observed by aLIGO.

These were calculated using the fiducial model in the rapid binary population synthesis element of COMPAS (v0.1).

We also provide a couple of variations described in the paper which demonstrate the effect of uncertainties in binary evolution.

Any work making use of these simulations should cite:

Stevenson, S., Vigna-Gómez, A., Mandel, I., Barrett, J. W., Neijssel, C. J., Perkins, D., & de Mink, S. E. (2017). Formation of the first three gravitational-wave observations through isolated binary evolution. Nature Communications, 8, 14906, doi: 10.1038/ncomms14906 (2017).

URL:	https://www.nature.com/articles/ncomms14906
ADS: 	https://ui.adsabs.harvard.edu/abs/2017NatCo...814906S/abstract
arXiv:	https://arxiv.org/abs/1704.01352

For any queries, email:

compas@star.sr.bham.ac.uk 

Compact Object Mergers: Population Astrophysics and Statistics (COMPAS)
Website:	www.compas.science

File headers
---------------

Headers:

In each file:
The first line gives the units of quantities
The second line gives the 'key' name of that quantity

_initials.txt

ID              #-- Non unique identifier in COMPAS
SEED            #-- unique identifier in COMPAS
mass1           #-- ZAMS mass of initially more massive star in [Msol]
mass2           #-- ZAMS mass of initially less massive star in [Msol]
separation      #-- ZAMS semi major axis in [AU]
eccentricity    #-- ZAMS eccentricity
sigma_kick      #-- 1D SN kick velocity dispersion hyper parameter in [kms^-1]
LBV_multiplier  #-- LBV wind mass loss hyper parameter
CE_Alpha        #-- Common envelope alpha hyper parameter
Metallicity1    #-- Fractional metallicitiy 0.0 < Z < 1.0 of primary star. Zsol = 0.02
Metallicity2    #-- Fractional metallicitiy 0.0 < Z < 1.0 of secondary star. Zsol = 0.02


_mergers.txt

ID                              #-- Non unique identifier in COMPAS
SEED                            #-- Unique identifier in COMPAS
separationInitial               #-- ZAMS semi major axis in AU
eccentricityInitial             #-- ZAMS eccentricity
separationPrior2ndSN            #-- Semi major axis at the moment of the second supernova in [AU]
eccentricityPrior2ndSN          #-- eccentricity at the moment of the second supernova in [AU]
relativeVelocity2ndSN           #-- Relative orbital velocity at the moment of the second supernova in [m s^-1]
separationDCOFormation          #-- Semi major axis at double compact object formation in [AU]
eccentricityDCOFormation        #-- Eccentricity at double compact object formation in [AU]
Metallicity1                    #-- Fractional metallicity 0.0 < Z < 1.0 of primary star. Zsol = 0.02
Metallicity2                    #-- Fractional metallicity 0.0 < Z < 1.0 of secondary star. Zsol = 0.02
M1ZAMS                          #-- ZAMS mass of initially more massive star in [Msol]
totalMassDCOFormation1          #-- Total mass of star 1 at DCO formation in [Msol]
HeCoreMassDCOFormation1         #-- Helium core mass of star 1 at DCO formation in [Msol]
coreMassDCOFormation1           #-- core mass of star 1 at DCO formation in [Msol]
HeCoreMassCE1                   #-- Helium core mass of star 1 at common envelope in [Msol]
coreMassCE1                     #-- Core mass of star 1 at common envelope in [Msol]
drawnKick1                      #-- Kick velocity drawn at SN of star 1 in [km s^-1]
M1                              #-- Final compact object mass of initially more massive star in [Msol] 
stellarType1                    #-- Final stellar type of initially more massive star (according to Hurley et al 2000). NS = 13, BH = 14
M2ZAMS                          #-- ZAMS mass of initially less massive star in [Msol]
totalMassDCOFormation2          #-- Total mass of star 2 at DCO formation in [Msol]
HeCoreMassDCOFormation2         #-- Helium core mass of star 2 at DCO formation in [Msol]
coreMassDCOFormation2           #-- core mass of star 2 at DCO formation in [Msol]
HeCoreMassCE2                   #-- Helium core mass at common envelope in [Msol]
coreMassCE2                     #-- core mass at common envelope in [Msol]
drawnKick2                      #-- Kick velocity drawn at SN of star 2 in [km s^-1]
M2                              #-- Final compact object mass of initially less massive star in [Msol]
stellarType2                    #-- Final stellar type of initially less massive star (according to Hurley et al 2000). 13 = NS, 14 = BH
tc                              #-- Coalescence time due to gravitational wave emission from time of double compact object formation in [Myrs]
tform                           #-- Time to formation of double compact object in [Myrs]
flbv                            #-- multiplicative constant for Luminous Blue Variable (LBV) mass loss rates 
sigmaKick                       #-- 1D SN kick velocity dispersion hyper parameter in [km s^-1]
CEalpha                         #-- Common envelope alpha hyper parameter
kickDirectionPower              #-- hyper parameter describing SN kick direction distribution (not varied)
wolfRayetMultiplier             #-- multiplicative constant for Wolf Rayet (WR) mass loss rates (not varied)
RLOFSecondaryAfterCEE           #-- Boolean. Whether the secondary experiences RLOF after Common Envelope Event
PrimaryMTCase                   #-- Primary Mass Transfer (MT) case. 0 = no mass transfer, 1 = case A, 2 = case B, 3 = case C 
SecondaryMTCase                 #-- Secondary Mass Transfer (MT) case. 0 = no mass transfer, 1 = case A, 2 = case B, 3 = case C
preSNeOrbitalEnergy1            #-- PreSN  Orbital energy in [Msol AU^2 yr^-2]
postSNeOrbitalEnergy1           #-- PostSN Orbital energy in [Msol AU^2 yr^-2]
preSNeOrbitalEnergy2            #-- PreSN  Orbital energy in [Msol AU^2 yr^-2]
postSNeOrbitalEnergy2           #-- PostSN Orbital energy in [Msol AU^2 yr^-2]
CEflag                          #-- Boolean. True if the system experienced at least one Common Envelope Event
lambda1CE                       #-- Lambda parameter describing envelope binding energy for star 1 (not varied)
lambda2CE                       #-- Lambda parameter describing envelope binding energy for star 2 (not varied)
EccentricityPreCEE              #-- Eccentricity before Common Envelope Event
EccentricityPostCEE             #-- Eccentricity after Common Envelope Event 
SemiMajorAxisPreCEE             #-- Semi-major axis before Common Envelope Event
SemiMajorAxisPostCEE            #-- Semi-major axis after Common Envelope Event
RL1to2PreCEE                    #-- Roche lobe as seen by the primary (q=M1/M2) at the onset of the Common Envelope Event in [Rsol]
RL1to2PostCEE                   #-- Roche lobe as seen by the primary (q=M1/M2) after the envelope has been expelled in a Common Envelope Event in [Rsol]
RL2to1PreCEE                    #-- Roche lobe as seen by the secondary (q=M2/M1) at the onset of the Common Envelope Event in [Rsol]
RL2to1PostCEE                   #-- Roche lobe as seen by the secondary (q=M2/M1) after the envelope has been expelled in a Common Envelope Event in [Rsol]
optimisticCEFlag                #-- Boolean. True if any of the stars did not have a clear core-envelope separation (optimistic) during a Common Envelope Event

Reading files in Python
-------------------------

The files can be read in Python with 

import numpy as np
data = np.genfromtxt("dataFile.txt", names=True)

print(data.dtype.names) #-- Prints the list of headers (see above for description)
print(data['M1ZAMS']) #-- Print ZAMS masses of primary stars.
