High-energy interactions of charged black holes in full general relativity. II. Near-extremal merger remnants and universality with the irreducible mass
Abstract
A viable model for the dense matter equation of state above the nuclear saturation density includes a hadron-to-quark phase transition at densities relevant to compact objects. In this case, stable hybrid hadron-quark stars can arise. An even more interesting scenario is one where the hadron-to-quark phase transition results in the emergence of a third branch of stable compact objects (in addition to white dwarfs and neutron stars). Inherent to the presence of a third family of compact stars is the existence of twin stars—hybrid stars with the same mass as the corresponding neutron stars but with smaller radii. Interestingly, the neutron star–twin star scenario is consistent with GW170817. If twin stars exist in nature, it raises a question about the mechanism that leads to their formation. Here, we explore gravitational collapse as a pathway to the formation of low-mass twin stars. We perform fully general relativistic simulations of the collapse of a stellar iron core, modeled as a cold degenerate gas, to investigate whether the end product is a neutron star or a twin star. Our simulations show that, even with unrealistically large perturbations in the initial conditions, the core bounces well below the hadron-to-quark phase transition density, if the initial total rest mass is in the twin star range. Following cooling, these configurations produce neutron stars. We find that twin stars can potentially form due to mass loss, e.g., through winds, from a slightly more massive hybrid star that was initially produced in the collapse of a more massive core, or if the maximum neutron star mass is below the Chandrasekhar mass limit. The challenge in producing twin stars in gravitational collapse, in conjunction with the fine-tuning required because of their narrow mass range, suggests the rarity of twin stars in nature.
Copyright and License
© 2025 American Physical Society.
Acknowledgement
We thank Vikram Manikantan for his feedback on figures displayed in this work, and Maria Mutz for useful discussions and feedback on the manuscript. This work was in part supported by NSF Grant No. PHY-2145421 and NASA Grant No. 80NSSC24K0771. This work was supported by Advanced Cyberinfrastructure Coordination Ecosystem: Services & Support (ACCESS) allocation TG-PHY190020 and Frontera allocation PHY23009. ACCESS is funded by NSF Awards No. 2138259, No. 2138286, No. 2138307, No. 2137603, and No. 2138296 under the Office of Advanced Cyberinfrastructure. The simulations were performed on stampede2 and frontera, funded by NSF Awards No. 1540931 and No. 1818253, respectively, at the Texas Advanced Computing Center (TACC).
Files
Name | Size | Download all |
---|---|---|
md5:f40fa1290d2b97c0c3c275d99aad737d
|
808.6 kB | Preview Download |
Additional details
- National Science Foundation
- PHY-2145421
- National Science Foundation
- 2138259
- National Science Foundation
- 2138286
- National Science Foundation
- 2138307
- National Science Foundation
- 2137603
- National Science Foundation
- 2138296
- National Science Foundation
- 1540931
- National Science Foundation
- 1818253
- National Aeronautics and Space Administration
- 80NSSC24K0771
- Frontera Energy
- PHY23009
- Advanced Cyberinfrastructure Coordination Ecosystem: Services & Support
- TG-PHY190020
- Texas Advanced Computing Center
- Available
-
2025-05-12Published online
- Caltech groups
- Division of Geological and Planetary Sciences (GPS)
- Publication Status
- Published