Hadron-quark phase transitions along proto-neutron-star evolution

The paper investigates the construction of hot hybrid equations of state (EOSs) to model hadron-quark phase transitions in proto-neutron stars, focusing on different evolutionary stages of a proto-hybrid star (a proto-neutron star with a quark matter core). The authors employ two distinct Skyrme effective interactions (SkI4 and Ska) for hadronic matter and the vector MIT bag model for quark matter. They consider three evolutionary stages: (1) the first stage immediately after formation, characterized by trapped neutrinos, an entropy per baryon S b = 1 and a lepton fraction Y l = 0.4; (2) the second stage after neutrino diffusion, characterized by S b = 2 and no neutrinos; and (3) the final cold stage at zero temperature.

For the first stage, the authors model the EOS by considering global conservation of the lepton fraction, which results in an extended mixed phase. For the second and final stages, they use both the Maxwell construction (for local electric-charge conservation, LCN) and the Gibbs construction (for global electric-charge conservation, GCN). A key methodological point is that, for the LCN case in the second stage, the authors allow the two phases to have different entropies per baryon to avoid a temperature discontinuity at the phase coexistence point, which would otherwise violate thermal equilibrium. They verify the thermodynamic consistency of this approach using the Clausius-Clapeyron equation, finding that the entropy difference between the hadronic and quark phases is of order unity (e.g., S b q = 3.04 for SkI4 and S b q = 2.81 for Ska, compared to S b h = 2).

The resulting hybrid EOSs are then used to solve the Tolman-Oppenheimer-Volkov (TOV) equations to describe hybrid star configurations. The authors find that, in the LCN case, the EOS for the first evolutionary stage predicts a larger maximum mass than the EOSs for the subsequent stages, while in the GCN case all three EOSs predict similar maximum masses. They also observe the emergence of thermal twin star solutions (stars with identical mass but different radii) in the LCN case for the second stage.

Finally, the authors study the structural evolution of proto-hybrid stars using constant rest mass sequences. They find that, regardless of whether electric charge is globally or locally conserved, the earlier stages of a hybrid star's life may play a crucial role in determining its maximum possible gravitational mass in later stages. Specifically, they show that there is no astrophysical path leading to the maximum-mass configurations predicted by the cold EOSs, because a star formed with a rest mass equal to the maximum rest mass of the cold case would become dynamically unstable during the intermediate (second) stage of evolution. This implies that even if a cold hybrid EOS predicts equilibrium configurations with masses that satisfy astronomical constraints, those configurations may never be realized in nature if any preceding configuration along the evolutionary sequence is dynamically unstable.

Improvements to AI Systems:

1. Thermodynamic Consistency Enforcement in Multi-Phase EOS Models

• Improvement: Integrate the paper’s method of allowing phase-dependent entropy (e.g., S b q not equal to S b h) into AI-driven equation-of-state (EOS) generators. Train the AI to automatically detect and resolve temperature discontinuities at phase boundaries by applying the Clausius-Clapeyron equation as a constraint.

• Capability: The AI can now produce physically consistent hybrid EOSs for finite-temperature astrophysical simulations, avoiding unphysical jumps in thermal equilibrium without manual intervention.

2. Evolutionary-Stage-Aware Stability Prediction

• Improvement: Enhance AI models that predict neutron star stability by incorporating multi-stage evolutionary sequences (trapped-neutrino, neutrino-free, cold) rather than relying solely on cold EOSs. Use the paper’s finding that intermediate stages can render maximum-mass cold configurations dynamically unstable.

• Capability: The AI can now flag “false stable” configurations by simulating the full thermal and compositional trajectory, thereby preventing overestimation of maximum gravitational mass in stellar evolution codes.

3. Twin-Star Solution Discovery via Constrained Optimization

• Improvement: Train AI to search for thermal twin star solutions (same mass, different radii) by explicitly varying entropy per baryon and lepton fraction in the LCN (local charge neutrality) case, using the paper’s observation of such solutions in the second stage.

• Capability: The AI can automatically generate candidate twin-star configurations for gravitational-wave and electromagnetic follow-up, aiding in the identification of phase-transition signatures.

4. Charge-Conservation-Mode-Aware EOS Interpolation

• Improvement: Build AI models that learn to switch between Maxwell (LCN) and Gibbs (GCN) constructions based on the evolutionary stage and neutrino trapping status, as demonstrated in the paper. This includes handling global vs. local charge conservation in the mixed phase.

• Capability: The AI can produce seamless EOS transitions for merger simulations, where both charge-conservation modes may be relevant at different timescales, improving the realism of binary neutron star merger models.

5. Rest-Mass Sequence Instability Early-Warning System

• Improvement: Implement an AI-based early-warning module that, given a cold EOS, computes constant rest-mass sequences across all evolutionary stages and identifies whether any intermediate stage violates dynamical stability. This directly uses the paper’s conclusion that no astrophysical path may lead to cold maximum-mass configurations.

• Capability: The AI can now pre-screen any proposed cold hybrid EOS for astrophysical viability, saving computational time in large-scale surveys of nuclear matter models.

6. Automated Entropy-Difference Calibration

• Improvement: Use the paper’s numerical results (e.g., S b q = 3.04 for SkI4) to train a regression model that predicts the entropy difference between hadronic and quark phases as a function of the Skyrme interaction parameters and bag constant.

• Capability: The AI can rapidly estimate entropy jumps for new EOS parameter sets without running full TOV or phase-coexistence calculations, accelerating model exploration.

7. Multi-Stage TOV Solver Integration

• Improvement: Enhance AI-driven TOV solvers to accept stage-dependent EOSs (with different lepton fractions, entropies, and charge-conservation rules) and output mass-radius relations for each stage, including the maximum mass and twin-star branches.

• Capability: The AI can produce complete evolutionary tracks of proto-hybrid stars in a single run, enabling direct comparison with observational constraints from supernova remnants and gravitational-wave events.

The new era of multi-messenger astronomy requires the accurate and self-consistent derivation of the nuclear equation of state at high temperature. In the present work, we focused on the calculation of hot hybrid equations of state, studying the different evolution stages of a proto-neutron star with a quark matter core (proto-hybrid star). For the hadronic matter we used two distinct Skyrme effective interactions, while for quark matter the well-known vector MIT bag model was employed. To model the era of trapped neutrinos in the system we considered the global conservation of lepton fraction which resulted in an equation of state with an extended mixed phase. For periods following the neutrino diffusion phase of a proto-neutron star, the equations of state were modelled using both the Maxwell and the Gibbs construction depending on the assumption for either local or global electric-charge conservation. With the use of the derived hybrid models, we solved the Tolman-Oppenheimer-Volkov equations to describe the corresponding hybrid star configurations. Finally, we investigated how the structure of proto-hybrid stars evolves, using constant rest mass sequences. We found that regardless of whether electric-charge is globally or locally conserved, the earlier stages of a hybrid star's life may play a crucial role on the determination of its maximum possible gravitational mass in later stages.