It's written in the massive stars: The role of stellar physics in the formation of black holes

In the age of gravitational-wave (GW) sources and newly discovered local black holes (BH) and neutron stars (NS), understanding the fate of stars is a key question. Not every massive star is expected to successfully explode as a supernova and leave behind a NS; some stars form BHs. The remnant depends on explosion physics but also on the final core structure, often summarized by the compactness parameter or iron core mass, where high values have been linked to BH formation. Several groups have reported similar patterns in these parameters as a function of mass, characterized by a prominent compactness peak followed by another peak at higher masses, pointing to a common underlying physical mechanism. Here, we investigate its origin by computing single-star models from 17 to 50 solar masses with MESA. The first and second compactness increases originate from core carbon and neon burning, respectively, becoming neutrino dominated, which enhances the core contraction and ultimately increases the iron-core mass and compactness. An early core neon ignition during carbon burning, and an early silicon ignition during oxygen burning help counter the core contraction and decrease the final iron core mass and compactness. Shell mergers between C/Ne and O-burning shells further decrease the compactness and we show that they are due to an enhanced entropy production in these layers. We find that the final structure of massive stars is not random but already written in their cores at core helium exhaustion. The same mechanisms determine the final structure of any star in this core mass range, including binary products, though binary interactions systematical shift the range of expected BH formation. Finally, we discuss the role of stellar physics uncertainties and how to apply these findings to studies of GW sources. [Abridged]