search

hamers / mse

Multiple Stellar Evolution -- A Population Synthesis Code for Multiple-Star Systems
GNU General Public License v3.0
9 stars 4 forks source link

Possibly wrong high mass black holes #46

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version b036225f5951c1be691051a3ddbf22758d995943 (with 0 supernova kicks)

These systems end up with high mass black holes, which probably (not sure!) might be wrong, especially system 42187. Examples - 5068 (31 MSun), 42187 (44 MSun), 63284 (31 MSun). I re-ran them with the lates version and still saw the same results.

index_system = 5068
configuration = "[2,2]"
N_bodies = 4
masses = [11.480709507969538, 3.6457095442634584, 7.160267781685773, 34.51804481639648]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [1.9077232939054682, 0.29143687635844623, 39.58910363532721]
eccentricities = [0.9244502029730212, 0.3037430079738928, 0.10758301005433966]
inclinations = [0.5408156367311555, 2.0549600465725244, 1.2114510361461213]
arguments_of_pericentre = [5.0265337306641475, 4.832210273222782, 0.664336346748185]
longitudes_of_ascending_node = [4.193160664920405, 5.666932986590698, 1.2585027759919258]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

index_system = 42187
configuration = "[2,2]"
N_bodies = 4
masses = [17.955099856415227, 25.43861143816699, 70.44729020382495, 70.54221254128075]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [2.2720248981610838, 3.4103402126941993, 30.21357327009001]
eccentricities = [0.7710879979782606, 0.7846935884744111, 0.35153389358305764]
inclinations = [1.5892470192609987, 1.5349852150676497, 1.9012925213102971]
arguments_of_pericentre = [2.6606954235213, 4.887883871569553, 3.689148814395094]
longitudes_of_ascending_node = [5.710763958698144, 2.0541324114926, 5.296973911249262]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

index_system = 63284
configuration = "[2,2]"
N_bodies = 4
masses = [33.10094308475692, 5.529951813777178, 65.87393448664203, 39.860439925592296]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [5.019474555774855, 0.4384528664467132, 72.94295687449474]
eccentricities = [0.8062429673991279, 0.5007754835232592, 0.4880394392603675]
inclinations = [1.356931270359145, 2.0678048806045957, 1.3406356779928115]
arguments_of_pericentre = [3.247602606101761, 0.6490004544830226, 0.7808334231283096]
longitudes_of_ascending_node = [6.244496779767114, 4.638393053125136, 1.4325320669797317]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]