However, what you really want must be specified different to get the full range of possible bright stars, especially if you are using a non-delta burst star formation history:
final_kstar1 = [14] final_kstar2 = [2,3,4,5,6,7,8,9,10,11,12,13,14] initBin_n, m_singles_n, m_binaries_n, n_singles_n, n_binaries_n = InitialBinaryTable.sampler('independent', final_kstar1, final_kstar2, binfrac_model=0.5, primary_model='kroupa01', ecc_model='sana12', porb_model='sana12', SF_start=13700.0, SF_duration=0.0, met=0.02, size=10000)
this produces:
I think the solution is to only place mass cuts for kstars in [10,11,12,13,14] and nothing else.
katiebreivik
commented
4 years ago
If you, e.g., want to get BH + Giant star companions, you would specify an initial population as:
final_kstar1 = [14] final_kstar2 = [2,3,4,5,6,] initBin_n, m_singles_n, m_binaries_n, n_singles_n, n_binaries_n = InitialBinaryTable.sampler('independent', final_kstar1, final_kstar2, binfrac_model=0.5, primary_model='kroupa01', ecc_model='sana12', porb_model='sana12', SF_start=13700.0, SF_duration=0.0, met=0.02, size=10000)this produces:
However, what you really want must be specified different to get the full range of possible bright stars, especially if you are using a non-delta burst star formation history:
final_kstar1 = [14] final_kstar2 = [2,3,4,5,6,7,8,9,10,11,12,13,14] initBin_n, m_singles_n, m_binaries_n, n_singles_n, n_binaries_n = InitialBinaryTable.sampler('independent', final_kstar1, final_kstar2, binfrac_model=0.5, primary_model='kroupa01', ecc_model='sana12', porb_model='sana12', SF_start=13700.0, SF_duration=0.0, met=0.02, size=10000)this produces:
I think the solution is to only place mass cuts for kstars in [10,11,12,13,14] and nothing else.