H-rich non-burning stars being directed to CO_HeMS grid
#218 opened on Dec 13, 2023
Repository metrics
- Stars
- (42 stars)
- PR merge metrics
- (PR metrics pending)
Description
H-rich_non_burning stars are being directed to the CO_HeMS step. These stars are then immediately re-directed from the CO_HeMS step, with mixed results.
Here is an example where the the binary continues to complete its evolution after the redirect step (the minority of cases):
And here is an example of the behavior for the majority of cases. The binary event is set as redirect_from_CO_HeMS two times consecutively, and its state is set to initial_RLOF in the detached step, which ends the evolution:
This issue is stemming from the fact that H-rich_non_burning are included in STAR_STATES_HE_RICH_EVOLVABLE in the flow, with the following comment:
# CE ejcetion happens istantanously, the star does not readjust before
# we infer the state, if core_definition_H_fraction=0.1 then surface_h1=0.1
# and the state is H-rich_non_burning which we stil want to evolve thorugh
# the step_CO_HeMS
Note: these stars are not designated in STAR_STATES_H_RICH, only STAR_STATES_HE_RICH_EVOLVABLE.
We have the following questions:
- Are these
H-rich_non_burningstars being classified correctly as such by the IF interpolator, i.e. should they actually be classified asH-rich_Shell_H_burningor similar? - If they are being classified correctly, should they be sent to
step_CO_HeMSat all?
Note: We should make sure the redirect_from_CO_HeMS event is not appearing twice consecutively in the evolutionary tracks regardless, as it is confusing for users to read/understand. @ezapartas mentioned this in issue #202.