POSYDON-code / POSYDON

POSYDON is a next-generation single and binary-star population synthesis code incorporating full stellar structure and evolution modeling with the use of MESA.
BSD 3-Clause "New" or "Revised" License
30 stars 19 forks source link

Undefined binaries in flow: ('NS', 'H-rich_Core_H_burning', 'disrupted', 'CC2') #458

Open mkruckow opened 6 days ago

mkruckow commented 6 days ago

When doing some tests for v2.0pre1, I have see several cases of

They happend both in IF and NN runs. It is questionable how a 'H-rich_Core_H_burning' is send to CC.

Here the full error of one example case:

Traceback (most recent call last):
  File "/home/users/k/kruckow/POSYDON/posydon/popsyn/binarypopulation.py", line 325, in _safe_evolve
    binary.evolve()
  File "/home/users/k/kruckow/POSYDON/posydon/binary_evol/binarystar.py", line 219, in evolve
    self.run_step()
  File "/home/users/k/kruckow/POSYDON/posydon/binary_evol/binarystar.py", line 238, in run_step
    raise ValueError("Undefined next step given stars/binary states {}.".format(total_state))
ValueError: Undefined next step given stars/binary states ('NS', 'H-rich_Core_H_burning', 'disrupted', 'CC2').

Failed Binary Initial Conditions:
S1 mass: 7.302102710931483
S2 mass: 0.7879377046773028
S1 state: H-rich_Core_H_burning
S2 state: H-rich_Core_H_burning
orbital period: 1719.5155572498547
eccentricity: 0.0
binary state: detached
binary event: detached
S1 natal kick array: [500.3825406264961 2.9365934718580298 1.2255496882454837 5.953462883339887]
S2 natal kick array: [None None None None]
maxbriel commented 2 days ago

@mkruckow Is this related to our v2 issue with the single stars and the EEPs? (PR#243) There we had to evolve the v2 single stars past the age of the Universe to get all their EEPs to show up.

Till what point do we evolve the v1 version of single stars?

mkruckow commented 2 days ago

Good point, this might in fact be the issue which we solved for v2 by evolving the stars further as it was done before (incl. v1). As we can't solve it on the grid data level for v1, shell we try to catch the error and/or mention it in the v1 documentation?