Multidimensional distributions¶
COSMIC implements multidimensionally distributed initial binaries according to Moe & Di Stefano 2017. The python code used in COSMIC to create this sample was written by Mads Sorenson, and is based on the IDL codes written to accompany Moe & Di Stefano 2017.
The multidimensional initial binary data is sampled in COSMIC as follows:
In [1]: from cosmic.sample.initialbinarytable import InitialBinaryTable
In [2]: from cosmic.sample.sampler import multidim
Tip
The multidimensional sampler arguments can be found here get_multidim_sampler() or by using the multidimensional sampler help call (help(multidim.get_multidim_sampler))
The random seed is used to reproduce your initial sample, since there are several stochastic processes involved in the muldimensional sample. As in the independent sampler, the final_kstar1 and final_kstar2 inputs are lists containing the kstar types that the evolved population should contain.
The multidimensional sample is generated as follows:
In [3]: InitialBinaries, mass_singles, mass_binaries, n_singles, n_binaries = InitialBinaryTable.sampler('multidim', final_kstar1=[11], final_kstar2=[11], rand_seed=2, nproc=1, SF_start=13700.0, SF_duration=0.0, met=0.02, size=10)
In [4]: print(InitialBinaries)
kstar_1 kstar_2 mass_1 mass_2 ... bhspin_1 bhspin_2 tphys binfrac
0 1.0 1.0 7.400447 3.295040 ... 0.0 0.0 0.0 0.824843
1 1.0 1.0 1.067083 1.032837 ... 0.0 0.0 0.0 0.408716
2 1.0 1.0 1.041342 0.809598 ... 0.0 0.0 0.0 0.408716
3 1.0 1.0 12.530773 9.774547 ... 0.0 0.0 0.0 0.872301
4 1.0 1.0 0.993994 0.987391 ... 0.0 0.0 0.0 0.407540
5 1.0 1.0 2.410236 1.072114 ... 0.0 0.0 0.0 0.521047
6 1.0 1.0 3.411454 1.891869 ... 0.0 0.0 0.0 0.593455
7 1.0 1.0 2.432400 2.010995 ... 0.0 0.0 0.0 0.521047
8 1.0 1.0 2.398650 0.975678 ... 0.0 0.0 0.0 0.521047
9 1.0 1.0 1.290047 0.986051 ... 0.0 0.0 0.0 0.421378
10 0.0 15.0 0.183103 0.000000 ... 0.0 0.0 0.0 1.000000
11 1.0 15.0 0.943782 0.000000 ... 0.0 0.0 0.0 1.000000
12 1.0 15.0 0.723793 0.000000 ... 0.0 0.0 0.0 1.000000
13 0.0 15.0 0.115933 0.000000 ... 0.0 0.0 0.0 1.000000
14 0.0 15.0 0.129371 0.000000 ... 0.0 0.0 0.0 1.000000
15 0.0 15.0 0.319392 0.000000 ... 0.0 0.0 0.0 1.000000
16 0.0 15.0 0.301592 0.000000 ... 0.0 0.0 0.0 1.000000
17 0.0 15.0 0.170359 0.000000 ... 0.0 0.0 0.0 1.000000
18 1.0 15.0 0.801332 0.000000 ... 0.0 0.0 0.0 1.000000
19 0.0 15.0 0.175066 0.000000 ... 0.0 0.0 0.0 1.000000
20 1.0 15.0 1.427620 0.000000 ... 0.0 0.0 0.0 1.000000
21 0.0 15.0 0.474281 0.000000 ... 0.0 0.0 0.0 1.000000
22 0.0 15.0 0.414853 0.000000 ... 0.0 0.0 0.0 1.000000
23 0.0 15.0 0.192693 0.000000 ... 0.0 0.0 0.0 1.000000
24 1.0 15.0 0.925249 0.000000 ... 0.0 0.0 0.0 1.000000
25 0.0 15.0 0.103579 0.000000 ... 0.0 0.0 0.0 1.000000
26 0.0 15.0 0.155512 0.000000 ... 0.0 0.0 0.0 1.000000
27 1.0 15.0 1.510590 0.000000 ... 0.0 0.0 0.0 1.000000
28 1.0 15.0 0.838802 0.000000 ... 0.0 0.0 0.0 1.000000
29 0.0 15.0 0.404639 0.000000 ... 0.0 0.0 0.0 1.000000
30 0.0 15.0 0.159762 0.000000 ... 0.0 0.0 0.0 1.000000
31 0.0 15.0 0.388535 0.000000 ... 0.0 0.0 0.0 1.000000
32 1.0 15.0 6.080843 0.000000 ... 0.0 0.0 0.0 1.000000
33 1.0 15.0 1.362545 0.000000 ... 0.0 0.0 0.0 1.000000
34 0.0 15.0 0.211245 0.000000 ... 0.0 0.0 0.0 1.000000
35 0.0 15.0 0.275150 0.000000 ... 0.0 0.0 0.0 1.000000
36 1.0 15.0 1.976093 0.000000 ... 0.0 0.0 0.0 1.000000
37 0.0 15.0 0.561326 0.000000 ... 0.0 0.0 0.0 1.000000
38 0.0 15.0 0.210801 0.000000 ... 0.0 0.0 0.0 1.000000
39 1.0 15.0 0.981359 0.000000 ... 0.0 0.0 0.0 1.000000
40 1.0 15.0 3.327484 0.000000 ... 0.0 0.0 0.0 1.000000
41 0.0 15.0 0.149577 0.000000 ... 0.0 0.0 0.0 1.000000
42 0.0 15.0 0.598543 0.000000 ... 0.0 0.0 0.0 1.000000
43 0.0 15.0 0.138489 0.000000 ... 0.0 0.0 0.0 1.000000
44 1.0 15.0 1.128125 0.000000 ... 0.0 0.0 0.0 1.000000
45 0.0 15.0 0.083211 0.000000 ... 0.0 0.0 0.0 1.000000
46 0.0 15.0 0.122140 0.000000 ... 0.0 0.0 0.0 1.000000
47 0.0 15.0 0.360367 0.000000 ... 0.0 0.0 0.0 1.000000
48 0.0 15.0 0.542877 0.000000 ... 0.0 0.0 0.0 1.000000
49 0.0 15.0 0.288410 0.000000 ... 0.0 0.0 0.0 1.000000
50 0.0 15.0 0.316816 0.000000 ... 0.0 0.0 0.0 1.000000
51 0.0 15.0 0.176489 0.000000 ... 0.0 0.0 0.0 1.000000
52 0.0 15.0 0.670359 0.000000 ... 0.0 0.0 0.0 1.000000
53 1.0 15.0 0.759643 0.000000 ... 0.0 0.0 0.0 1.000000
[54 rows x 38 columns]
Note
NOTE that in the multidimensional case, the binary fraction is a parameter in the sample. This results in the size of the initial binary data matching the size provided to the sampler. As in the independent sampling case, we keep track of the total sampled mass of singles and binaries as well as the total number of single and binary stars to scale the simulated population to astrophysical populations.
Note
NOTE that you can also keep singles for the multidim sampelr as well. As with the independent sampler, the singles will be added onto the end of the InitialBinaryTable where kstar_1 will host the singles, kstar_2 will be filled with 15s only, and all orbital properties (e.g. porb or ecc) will be indicated with -1.