pith. sign in

arxiv: 1906.10234 · v1 · pith:Q3OGVUO4new · submitted 2019-06-24 · 🌌 astro-ph.GA · astro-ph.HE

The full evolution of supernova remnants in low and high density ambient media

Santiago Jimenez , Guillermo Tenorio-Tagle , Sergiy Silich This is my paper

Pith reviewed 2026-05-25 16:49 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.HE
keywords supernova remnantsradiative coolingthin-shell approximationshock dynamicsinterstellar mediumhigh-density mediaSedov-Taylor phase
0
0 comments X

The pith

In dense media above 5 times 10 to the fifth particles per cubic centimeter, radiative cooling stops the reverse shock before it reaches the explosion center.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper models the complete evolution of supernova remnants from the ejecta-dominated phase through the momentum-dominated stage in homogeneous ambient media of varying densities. It employs a thin-shell approximation that includes the ejected gas configuration and radiative cooling to track shock dynamics. In high-density cases the cooling causes the thermal pressure to drop rapidly, which prevents the Sedov-Taylor phase and stops the reverse shock short of the center. A reader would care because this changes how much mechanical and thermal feedback supernovae deliver to dense interstellar environments such as star-forming regions.

Core claim

For ambient densities n0 greater than 5 times 10 to the 5 per cubic centimeter, strong radiative cooling causes the thermal pressure in the shocked gas to fall rapidly enough that the reverse shock never reaches the center of the explosion, thereby inhibiting the Sedov-Taylor stage and limiting the overall feedback that the remnant can provide.

What carries the argument

The thin-shell approximation that incorporates the ejected-gas configuration and radiative cooling to evolve the remnant from ejecta-dominated to momentum-dominated stages across a wide density range.

If this is right

  • Supernova remnants supply significantly less feedback to high-density interstellar environments than they do at lower densities.
  • The Sedov-Taylor expansion phase is skipped entirely once radiative cooling becomes dominant.
  • The remnant reaches the momentum-dominated stage with a reduced total energy and momentum budget.
  • The internal structure of the remnant is altered because the reverse shock fails to cross the center.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This density threshold may mark a transition where supernova feedback becomes inefficient at dispersing dense molecular clouds.
  • Galactic evolution models that assume a universal Sedov-Taylor phase may overestimate energy injection in the densest star-forming regions.
  • Comparisons between observed supernova remnant sizes and velocities in varying density environments could directly test the predicted cutoff.

Load-bearing premise

The thin-shell approximation remains valid and accurately captures the dynamics when strong radiative cooling is present in high-density homogeneous media.

What would settle it

A high-resolution hydrodynamic simulation or direct observation of a young supernova remnant in an ambient medium denser than 5 times 10 to the 5 cm^{-3} that shows whether the reverse shock reaches the geometric center at late times.

Figures

Figures reproduced from arXiv: 1906.10234 by Guillermo Tenorio-Tagle, Santiago Jimenez, Sergiy Silich.

Figure 1
Figure 1. Figure 1: The structure of a SNR. The left panel shows the initial condition at t = t0. The right panel presents the structure of the SNR at a later time t > t0. See the text for a discussion on the labels of this scheme. stage, a very thin, cold and dense shell is formed at the outer edge of the SNR. In order to preserve pressure, the density increases in response to the sudden fall of the post-shock temperature. T… view at source ↗
Figure 2
Figure 2. Figure 2: The evolution of the shocks radii for a SNR with E0 = 1051 erg, n0 = 1 cm−3 , Mej = 3M⊙ and index n = 2. The top panel presents the case of the leading shock RLS and the bottom panel the reverse shock RRS. The solid lines show our results and the dashed and dotted lines are the analytic and numerical radii obtained by TM99, respectively. The starred variables at the axis are dimensionless variables as defi… view at source ↗
Figure 3
Figure 3. Figure 3: Evolution of the total, thermal, kinetic and radiated energies (solid, dotted, dashed, and dash-dotted lines, respectively) for a SNR evolving in the ambient medium with density n0 = 1 cm−3 . The three evolutionary stages are separated by the dashed and solid vertical lines. The dashed vertical line marks the moment when the reverse shock reaches the center of the explosion and the solid vertical line the … view at source ↗
Figure 4
Figure 4. Figure 4: The rate of energy loss (left panel) from the outer (dotted line) and inner shells (solid line) for the case of n0 = 1 cm−3 . Right panel: the leading (dotted line) and reverse (solid line) shocks velocities as function of time. Note that V˜RS is the reverse shock velocity in the frame of the unshocked ejecta (see equation 8) and VLS is the leading shock velocity in the rest frame. The vertical lines indic… view at source ↗
Figure 5
Figure 5. Figure 5: The kinetic (dashed line) and thermal (solid line) en￾ergies of a SNR evolving in an ambient medium with density n0 = 107 cm−3 . after the explosion, which is the time when the leading shock becomes radiative. This leads to the rapid remnant evolution as the thermal energy dramatically decreases in a short time interval. The fact that the remnant energies are decaying for most of the time covered by our ca… view at source ↗
Figure 6
Figure 6. Figure 6: The shocked ambient gas and the shocked ejecta densities and the cooling rates as a function of time. The left and right upper panels present the shocked ejecta ns2 and shocked ambient gas ns1 densities for different ambient gas densities (shown in the legend in units of cm−3 ). The left and right lower panels present the corresponding rates of energy losses Q2 and Q1. . 4 SNR EVOLUTION IN DIFFERENT AMBIEN… view at source ↗
Figure 7
Figure 7. Figure 7: The thermal energy of SNRs evolving into an ISM with different densities (listed in the legend in units of cm−3 ). the reverse shock at larger distances from the center of the explosion. The density of the shocked ambient gas ns1 (right upper panel) is initially close to the adiabatic strong shock limit 4n0. However, it increases orders of magnitude upon strong radiative cooling at the transition time tsf … view at source ↗
Figure 8
Figure 8. Figure 8: Top panel: Fraction of the thermalized ejecta Ms2/Mej as a function of time for different values of the ambient gas density (shown in the legend). The open squares indicate the time when the leading shock becomes radiative (i.e., when t = tsf ) for each case. For low density models this occur after full thermalization of the ejecta (Ms2/Mej = 1), in contrast with the large density models when Ms2/Mej < 1 a… view at source ↗
Figure 9
Figure 9. Figure 9: Evolution of the reverse (solid line) and leading (dotted line) shock radii for the densities n0 = 1, 104 , 5 × 105 and 107 cm−3 , respectively. The radiative cooling impact the reverse shock dynam￾ics. For low density cases, the reverse shock promptly reaches the center of the explosion (see upper panels of [PITH_FULL_IMAGE:figures/full_fig_p009_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: The total momentum deposited by the shocked gas for different values of the ambient gas density. The horizontal line shows the initial momentum of the ejecta pej . 11 1 1 11    1   1     1 11 1 1 11    1    1 11 1 [PITH_FULL_IMAGE:figures/full_fig_p010_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: The evolution of the reverse shock position for models with different chemical compositions. Panels a, b,c, and d present the reverse shock position as a function of time for ambient gas densities 104 cm−3 , 5 × 105 cm−3 , 3.5 × 106 cm−3 and 107 cm−3 , respectively. MNRAS 000, 1–14 (2019) [PITH_FULL_IMAGE:figures/full_fig_p010_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: The time evolution of the leading shock velocity VLS for the cases presented in the legend. The right y-axis shows the post-shock temperature TLS. The solid, dotted, dashed and dash-dotted vertical lines present the thin-shell formation time tsf for the cases (105 cm−3 , M1), (105 cm−3 , M4), (1 cm−3 , M1) and (1 cm−3 , M4), respectively. istration number 613136) and by the Sistema Nacional de Investigado… view at source ↗
read the original abstract

Supernova explosions and their remnants (SNRs) drive important feedback mechanisms that impact considerably the galaxies that host them. Then, the knowledge of the SNRs evolution is of paramount importance in the understanding of the structure of the interstellar medium (ISM) and the formation and evolution of galaxies. Here we study the evolution of SNRs in homogeneous ambient media from the initial, ejecta-dominated phase, to the final, momentum-dominated stage. The numerical model is based on the Thin-Shell approximation and takes into account the configuration of the ejected gas and radiative cooling. It accurately reproduces well known analytic and numerical results and allows one to study the SNR evolution in ambient media with a wide range of densities $n_{0}$. It is shown that in the high density cases, strong radiative cooling alters noticeably the shock dynamics and inhibits the Sedov-Taylor stage, thus limiting significantly the feedback that SNRs provide to such environments. For $n_{0}>5 \times 10^{5}$ cm$^{-3}$, the reverse shock does not reach the center of the explosion due to the rapid fall of the thermal pressure in the shocked gas caused by strong radiative cooling.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript develops a thin-shell numerical model for supernova remnant evolution from the ejecta-dominated phase through the momentum-dominated stage in homogeneous media. The model incorporates ejecta configuration and radiative cooling, reproduces standard analytic and numerical benchmarks at moderate densities, and concludes that for n0 > 5×10^5 cm^{-3} strong cooling causes rapid thermal-pressure collapse, inhibiting the Sedov-Taylor phase and preventing the reverse shock from reaching the explosion center.

Significance. If the thin-shell closure holds under strong radiative losses, the work supplies an efficient tool for mapping SNR feedback across the full density range encountered in star-forming regions. The approach correctly identifies that cooling can truncate the energy-conserving stage and thereby reduce momentum injection, a result with direct implications for ISM structure and galactic evolution models.

major comments (2)
  1. [Model description] Model description (thin-shell equations): the headline result that the reverse shock fails to reach the center for n0>5×10^5 cm^{-3} follows directly from integrating the thin-shell momentum and energy equations once radiative losses are inserted; however, no test is presented showing that the thin-shell geometry and contact-discontinuity tracking remain accurate once the cooling length becomes comparable to the shell thickness.
  2. [Validation section] Validation section: the statement that the code 'accurately reproduces well known analytic and numerical results' is supported only for lower-density regimes; no 1D or 2D hydrodynamical comparison is reported at n0>10^5 cm^{-3}, the regime in which the new dynamical claim is made.
minor comments (1)
  1. [Methods] Notation for the cooling function and the precise definition of the thin-shell thickness should be stated explicitly in the methods to allow independent reproduction.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive report and the opportunity to address the points raised. We respond to each major comment below and indicate where revisions will be made to the manuscript.

read point-by-point responses

  1. Referee: [Model description] Model description (thin-shell equations): the headline result that the reverse shock fails to reach the center for n0>5×10^5 cm^{-3} follows directly from integrating the thin-shell momentum and energy equations once radiative losses are inserted; however, no test is presented showing that the thin-shell geometry and contact-discontinuity tracking remain accurate once the cooling length becomes comparable to the shell thickness.

    Authors: We agree that the headline result follows directly from the thin-shell equations with radiative losses included. The thin-shell approximation is a standard and widely used closure for SNR evolution studies; our implementation tracks the contact discontinuity and incorporates cooling self-consistently. However, we acknowledge that an explicit test of geometric fidelity when the cooling length approaches the shell thickness is not provided. Such a test would require full 1D/2D hydrodynamical runs at the relevant densities, which lies outside the scope of the present thin-shell study. We will revise the model-description section to include an explicit discussion of this assumption and its potential limitations in the high-cooling regime. revision: partial

  2. Referee: [Validation section] Validation section: the statement that the code 'accurately reproduces well known analytic and numerical results' is supported only for lower-density regimes; no 1D or 2D hydrodynamical comparison is reported at n0>10^5 cm^{-3}, the regime in which the new dynamical claim is made.

    Authors: The analytic benchmarks we reproduce (Sedov-Taylor, ejecta-dominated phase) are density-independent in their scaling relations, and our code matches them at the moderate densities where those solutions apply. The new dynamical claim at n0 > 5×10^5 cm^{-3} arises precisely from the inclusion of strong radiative losses, which are absent from the standard benchmarks. We will revise the validation section to (i) state the density range of the comparisons explicitly and (ii) clarify that the high-density behavior constitutes a model prediction under the thin-shell closure rather than a direct numerical validation. revision: partial

Circularity Check

0 steps flagged

No significant circularity; standard thin-shell model with independent validation claims

full rationale

The paper builds its results on the established Thin-Shell approximation plus radiative cooling, explicitly stating that the model reproduces known analytic and numerical results. No load-bearing steps reduce by construction to fitted parameters, self-definitions, or self-citation chains. The central claim about reverse-shock behavior at high n0 follows from the model's integrated equations without the derivation being equivalent to its inputs. This is a normal non-finding for a numerical implementation of a standard framework.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The model rests on the thin-shell approximation and inclusion of radiative cooling in homogeneous media; no free parameters or invented entities are mentioned in the abstract.

axioms (1)
  • domain assumption The thin-shell approximation accurately models SNR evolution including radiative cooling effects across the density range studied.
    Stated as the basis of the numerical model in the abstract.

pith-pipeline@v0.9.0 · 5741 in / 1207 out tokens · 31150 ms · 2026-05-25T16:49:54.556150+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

62 extracted references · 62 canonical work pages · 1 internal anchor

  1. [1]

    write newline

    " write newline "" before.all 'output.state := FUNCTION fin.entry write newline FUNCTION new.block output.state before.all = 'skip after.block 'output.state := if FUNCTION new.sentence output.state after.block = 'skip output.state before.all = 'skip after.sentence 'output.state := if if FUNCTION not #0 #1 if FUNCTION and 'skip pop #0 if FUNCTION or pop #1...

    work page

  2. [2]

    V., Leitner S

    Agertz O., Kravtsov A. V., Leitner S. N., Gnedin N. Y., 2013, @doi [ ] 10.1088/0004-637X/770/1/25 , http://cdsads.u-strasbg.fr/abs/2013ApJ...770...25A 770, 25

    work page doi:10.1088/0004-637x/770/1/25 2013

  3. [3]

    Bianchi S., Schneider R., 2007, @doi [ ] 10.1111/j.1365-2966.2007.11829.x , http://adsabs.harvard.edu/abs/2007MNRAS.378..973B 378, 973

    work page doi:10.1111/j.1365-2966.2007.11829.x 2007

  4. [4]

    S., Silich S

    Bisnovatyi-Kogan G. S., Silich S. A., 1995, @doi [Reviews of Modern Physics] 10.1103/RevModPhys.67.661 , http://adsabs.harvard.edu/abs/1995RvMP...67..661B 67, 661

    work page doi:10.1103/revmodphys.67.661 1995

  5. [5]

    M., Wright E

    Blondin J. M., Wright E. B., Borkowski K. J., Reynolds S. P., 1998, @doi [ ] 10.1086/305708 , http://adsabs.harvard.edu/abs/1998ApJ...500..342B 500, 342

    work page doi:10.1086/305708 1998

  6. [6]

    J., Lyerly W

    Borkowski K. J., Lyerly W. J., Reynolds S. P., 2001, @doi [ ] 10.1086/319011 , http://adsabs.harvard.edu/abs/2001ApJ...548..820B 548, 820

    work page doi:10.1086/319011 2001

  7. [7]

    A., 1977, @doi [Annual Review of Astronomy and Astrophysics] 10.1146/annurev.aa.15.090177.001135 , 15, 175

    Chevalier R. A., 1977, @doi [ ] 10.1146/annurev.aa.15.090177.001135 , http://adsabs.harvard.edu/abs/1977ARA

    work page doi:10.1146/annurev.aa.15.090177.001135 1977

  8. [8]

    A., 1982a, @doi [ ] 10.1086/160126 , http://adsabs.harvard.edu/abs/1982ApJ...258..790C 258, 790

    Chevalier R. A., 1982a, @doi [ ] 10.1086/160126 , http://adsabs.harvard.edu/abs/1982ApJ...258..790C 258, 790

    work page doi:10.1086/160126

  9. [9]

    A., 1982b, @doi [ ] 10.1086/160167 , http://adsabs.harvard.edu/abs/1982ApJ...259..302C 259, 302

    Chevalier R. A., 1982b, @doi [ ] 10.1086/160167 , http://adsabs.harvard.edu/abs/1982ApJ...259..302C 259, 302

    work page doi:10.1086/160167

  10. [10]

    F., McKee C

    Cioffi D. F., McKee C. F., Bertschinger E., 1988, @doi [ ] 10.1086/166834 , http://adsabs.harvard.edu/abs/1988ApJ...334..252C 334, 252

    work page doi:10.1086/166834 1988

  11. [11]

    T., McKee C

    Draine B. T., McKee C. F., 1993, @doi [ ] 10.1146/annurev.aa.31.090193.002105 , http://adsabs.harvard.edu/abs/1993ARA

    work page doi:10.1146/annurev.aa.31.090193.002105 1993

  12. [12]

    T., Salpeter E

    Draine B. T., Salpeter E. E., 1979, @doi [ ] 10.1086/157165 , http://cdsads.u-strasbg.fr/abs/1979ApJ...231...77D 231, 77

    work page doi:10.1086/157165 1979

  13. [13]

    G., 2017, @doi [ ] 10.3847/1538-4357/836/1/80 , http://adsabs.harvard.edu/abs/2017ApJ...836...80E 836, 80

    Elmegreen B. G., 2017, @doi [ ] 10.3847/1538-4357/836/1/80 , http://adsabs.harvard.edu/abs/2017ApJ...836...80E 836, 80

    work page doi:10.3847/1538-4357/836/1/80 2017

  14. [14]

    C., Stewart G

    Fabian A. C., Stewart G. C., Brinkmann W., 1983, in Danziger J., Gorenstein P., eds, IAU Symposium Vol. 101, Supernova Remnants and their X-ray Emission. pp 119--124

    work page 1983

  15. [15]

    Falle S. A. E. G., 1975, @doi [ ] 10.1093/mnras/172.1.55 , http://cdsads.u-strasbg.fr/abs/1975MNRAS.172...55F 172, 55

    work page doi:10.1093/mnras/172.1.55 1975

  16. [16]

    Ferreira S. E. S., de Jager O. C., 2008, @doi [ ] 10.1051/0004-6361:20077824 , http://cdsads.u-strasbg.fr/abs/2008A

    work page doi:10.1051/0004-6361:20077824 2008

  17. [17]

    Franco J., Tenorio-Tagle G., Bodenheimer P., Rozyczka M., 1991, @doi [ ] 10.1086/132884 , http://adsabs.harvard.edu/abs/1991PASP..103..803F 103, 803

    work page doi:10.1086/132884 1991

  18. [18]

    Gaffet B., 1978, @doi [ ] 10.1086/156506 , http://adsabs.harvard.edu/abs/1978ApJ...225..442G 225, 442

    work page doi:10.1086/156506 1978

  19. [19]

    F., 1973, @doi [ ] 10.1093/mnras/161.1.47 , http://adsabs.harvard.edu/abs/1973MNRAS.161...47G 161, 47

    Gull S. F., 1973, @doi [ ] 10.1093/mnras/161.1.47 , http://adsabs.harvard.edu/abs/1973MNRAS.161...47G 161, 47

    work page doi:10.1093/mnras/161.1.47 1973

  20. [20]

    F., 1975, @doi [ ] 10.1093/mnras/171.2.263 , http://adsabs.harvard.edu/abs/1975MNRAS.171..263G 171, 263

    Gull S. F., 1975, @doi [ ] 10.1093/mnras/171.2.263 , http://adsabs.harvard.edu/abs/1975MNRAS.171..263G 171, 263

    work page doi:10.1093/mnras/171.2.263 1975

  21. [21]

    Haid S., Walch S., Naab T., Seifried D., Mackey J., Gatto A., 2016, @doi [ ] 10.1093/mnras/stw1082 , http://adsabs.harvard.edu/abs/2016MNRAS.460.2962H 460, 2962

    work page doi:10.1093/mnras/stw1082 2016

  22. [22]

    Hamilton A. J. S., Sarazin C. L., 1984, @doi [ ] 10.1086/162145 , http://adsabs.harvard.edu/abs/1984ApJ...281..682H 281, 682

    work page doi:10.1086/162145 1984

  23. [23]

    C., Krause O., Rieke G

    Hines D. C., Krause O., Rieke G. H., Fan X., Blaylock M., Neugebauer G., 2006, @doi [ ] 10.1086/504109 , http://adsabs.harvard.edu/abs/2006ApJ...641L..85H 641, L85

    work page doi:10.1086/504109 2006

  24. [24]

    P., Rakowski C

    Hughes J. P., Rakowski C. E., Burrows D. N., Slane P. O., 2000, @doi [ ] 10.1086/312438 , http://cdsads.u-strasbg.fr/abs/2000ApJ...528L.109H 528, L109

    work page doi:10.1086/312438 2000

  25. [25]

    M., 2012, @doi [ ] 10.1088/0004-637X/746/2/130 , http://adsabs.harvard.edu/abs/2012ApJ...746..130H 746, 130

    Hwang U., Laming J. M., 2012, @doi [ ] 10.1088/0004-637X/746/2/130 , http://adsabs.harvard.edu/abs/2012ApJ...746..130H 746, 130

    work page doi:10.1088/0004-637x/746/2/130 2012

  26. [26]

    C., 2015, @doi [ ] 10.1088/0004-637X/802/2/99 , http://adsabs.harvard.edu/abs/2015ApJ...802...99K 802, 99

    Kim C.-G., Ostriker E. C., 2015, @doi [ ] 10.1088/0004-637X/802/2/99 , http://adsabs.harvard.edu/abs/2015ApJ...802...99K 802, 99

    work page doi:10.1088/0004-637x/802/2/99 2015

  27. [27]

    M., Hwang U., 2003, @doi [ ] 10.1086/378268 , http://adsabs.harvard.edu/abs/2003ApJ...597..347L 597, 347

    Laming J. M., Hwang U., 2003, @doi [ ] 10.1086/378268 , http://adsabs.harvard.edu/abs/2003ApJ...597..347L 597, 347

    work page doi:10.1086/378268 2003

  28. [28]

    Lecture Notes 1997 Swiss Society for Astrophysics and Astronomy

    LeVeque R., Steiner O., Gautschy A., Mihalas D., Dorfi E., M \"u ller E., 2006, Computational Methods for Astrophysical Fluid Flow: Saas-Fee Advanced Course 27. Lecture Notes 1997 Swiss Society for Astrophysics and Astronomy. Saas-Fee Advanced Course, Springer Berlin Heidelberg, https://books.google.com.mx/books?id=xOYFCAAAQBAJ

    work page 2006

  29. [29]

    P., Cen R., Bryan G

    Li M., Ostriker J. P., Cen R., Bryan G. L., Naab T., 2015, @doi [ ] 10.1088/0004-637X/814/1/4 , http://adsabs.harvard.edu/abs/2015ApJ...814....4L 814, 4

    work page doi:10.1088/0004-637x/814/1/4 2015

  30. [30]

    Dust Production in the High-Redshift Universe

    Marchenko S. V., 2006, in Lamers H. J. G. L. M., Langer N., Nugis T., Annuk K., eds, Astronomical Society of the Pacific Conference Series Vol. 353, Stellar Evolution at Low Metallicity: Mass Loss, Explosions, Cosmology. p. 299 ( @eprint astro-ph/0511147 )

    work page internal anchor Pith review Pith/arXiv arXiv 2006

  31. [31]

    Mart \' nez-Gonz \'a lez S., Tenorio-Tagle G., Silich S., 2016, @doi [ ] 10.3847/0004-637X/816/1/39 , http://cdsads.u-strasbg.fr/abs/2016ApJ...816...39M 816, 39

    work page doi:10.3847/0004-637x/816/1/39 2016

  32. [32]

    Martizzi D., Faucher-Gigu \`e re C.-A., Quataert E., 2015, @doi [ ] 10.1093/mnras/stv562 , http://cdsads.u-strasbg.fr/abs/2015MNRAS.450..504M 450, 504

    work page doi:10.1093/mnras/stv562 2015

  33. [33]

    F., Ostriker J

    McKee C. F., Ostriker J. P., 1977, @doi [ ] 10.1086/155667 , http://adsabs.harvard.edu/abs/1977ApJ...218..148M 218, 148

    work page doi:10.1086/155667 1977

  34. [34]

    F., Truelove J

    McKee C. F., Truelove J. K., 1995, @doi [ ] 10.1016/0370-1573(94)00106-D , http://cdsads.u-strasbg.fr/abs/1995PhR...256..157M 256, 157

    work page doi:10.1016/0370-1573(94)00106-d 1995

  35. [35]

    C., Marinacci F., 2017, @doi [ ] 10.1093/mnras/stx467 , http://cdsads.u-strasbg.fr/abs/2017MNRAS.468.1505M 468, 1505

    McKinnon R., Torrey P., Vogelsberger M., Hayward C. C., Marinacci F., 2017, @doi [ ] 10.1093/mnras/stx467 , http://cdsads.u-strasbg.fr/abs/2017MNRAS.468.1505M 468, 1505

    work page doi:10.1093/mnras/stx467 2017

  36. [36]

    R., Dwek E., Slavin J

    Micelotta E. R., Dwek E., Slavin J. D., 2016, @doi [ ] 10.1051/0004-6361/201527350 , http://adsabs.harvard.edu/abs/2016A

    work page doi:10.1051/0004-6361/201527350 2016

  37. [37]

    J., 2015, @doi [ ] 10.1051/0004-6361/201525644 , http://cdsads.u-strasbg.fr/abs/2015A

    Micha owski M. J., 2015, @doi [ ] 10.1051/0004-6361/201525644 , http://cdsads.u-strasbg.fr/abs/2015A

    work page doi:10.1051/0004-6361/201525644 2015

  38. [38]

    Dover Books on Physics, Dover Publications, https://books.google.com.mx/books?id=GVK8AQAAQBAJ

    Mihalas D., Mihalas B., 2013, Foundations of Radiation Hydrodynamics. Dover Books on Physics, Dover Publications, https://books.google.com.mx/books?id=GVK8AQAAQBAJ

    work page 2013

  39. [39]

    L., Dunne L., Eales S

    Morgan H. L., Dunne L., Eales S. A., Ivison R. J., Edmunds M. G., 2003, @doi [ ] 10.1086/379639 , http://adsabs.harvard.edu/abs/2003ApJ...597L..33M 597, L33

    work page doi:10.1086/379639 2003

  40. [40]

    Nozawa T., Kozasa T., Tominaga N., Maeda K., Umeda H., Nomoto K., Krause O., 2010, @doi [ ] 10.1088/0004-637X/713/1/356 , http://adsabs.harvard.edu/abs/2010ApJ...713..356N 713, 356

    work page doi:10.1088/0004-637x/713/1/356 2010

  41. [41]

    P., McKee C

    Ostriker J. P., McKee C. F., 1988, @doi [Reviews of Modern Physics] 10.1103/RevModPhys.60.1 , http://adsabs.harvard.edu/abs/1988RvMP...60....1O 60, 1

    work page doi:10.1103/revmodphys.60.1 1988

  42. [42]

    Petruk O., 2006, ArXiv Astrophysics e-prints, http://adsabs.harvard.edu/abs/2006astro.ph..4057P

    work page 2006

  43. [43]

    Cambridge University Press, https://books.google.com.mx/books?id=1aAOdzK3FegC

    Press W., 2007, Numerical Recipes 3rd Edition: The Art of Scientific Computing. Cambridge University Press, https://books.google.com.mx/books?id=1aAOdzK3FegC

    work page 2007

  44. [44]

    C., Cox D

    Raymond J. C., Cox D. P., Smith B. W., 1976, @doi [ ] 10.1086/154170 , http://adsabs.harvard.edu/abs/1976ApJ...204..290R 204, 290

    work page doi:10.1086/154170 1976

  45. [45]

    Rowlands K., Gomez H. L., Dunne L., Arag \'o n-Salamanca A., Dye S., Maddox S., da Cunha E., van der Werf P., 2014, @doi [ ] 10.1093/mnras/stu605 , http://cdsads.u-strasbg.fr/abs/2014MNRAS.441.1040R 441, 1040

    work page doi:10.1093/mnras/stu605 2014

  46. [46]

    M., Kosenko D., Kaastra J

    Schure K. M., Kosenko D., Kaastra J. S., Keppens R., Vink J., 2009, @doi [ ] 10.1051/0004-6361/200912495 , http://adsabs.harvard.edu/abs/2009A

    work page doi:10.1051/0004-6361/200912495 2009

  47. [47]

    I., 1946, Dokl

    Sedov L. I., 1946, Dokl. Akad. Nauk SSSR, http://adsabs.harvard.edu/abs/1999A

    work page 1946

  48. [48]

    A., 1992, @doi [ ] 10.1007/BF00646764 , http://adsabs.harvard.edu/abs/1992Ap

    Silich S. A., 1992, @doi [ ] 10.1007/BF00646764 , http://adsabs.harvard.edu/abs/1992Ap

    work page doi:10.1007/bf00646764 1992

  49. [49]

    Silich S., Tenorio-Tagle G., 2018, @doi [ ] 10.1093/mnras/sty1383 , http://adsabs.harvard.edu/abs/2018MNRAS.478.5112S 478, 5112

    work page doi:10.1093/mnras/sty1383 2018

  50. [50]

    S., Seward F

    Slane P., Chen Y., Schulz N. S., Seward F. D., Hughes J. P., Gaensler B. M., 2000, @doi [ ] 10.1086/312589 , http://adsabs.harvard.edu/abs/2000ApJ...533L..29S 533, L29

    work page doi:10.1086/312589 2000

  51. [51]

    A., 2017, @doi [ ] 10.1093/mnras/stw2978 , http://adsabs.harvard.edu/abs/2017MNRAS.465.3793T 465, 3793

    Tang X., Chevalier R. A., 2017, @doi [ ] 10.1093/mnras/stw2978 , http://adsabs.harvard.edu/abs/2017MNRAS.465.3793T 465, 3793

    work page doi:10.1093/mnras/stw2978 2017

  52. [52]

    R., 2009, @doi [ ] 10.1111/j.1365-2966.2009.15167.x , http://adsabs.harvard.edu/abs/2009MNRAS.398..607V 398, 607

    Tang S., Wang Q. D., 2009, @doi [ ] 10.1111/j.1365-2966.2009.15099.x , http://cdsads.u-strasbg.fr/abs/2009MNRAS.397.2106T 397, 2106

    work page doi:10.1111/j.1365-2966.2009.15099.x 2009

  53. [53]

    Tenorio-Tagle G., Bodenheimer P., Franco J., Rozyczka M., 1990, , http://adsabs.harvard.edu/abs/1990MNRAS.244..563T 244, 563

    work page 1990

  54. [54]

    Tenorio-Tagle G., Rozyczka M., Franco J., Bodenheimer P., 1991, @doi [ ] 10.1093/mnras/251.2.318 , http://adsabs.harvard.edu/abs/1991MNRAS.251..318T 251, 318

    work page doi:10.1093/mnras/251.2.318 1991

  55. [55]

    Tenorio-Tagle G., Silich S., Mart \' nez-Gonz \'a lez S., Mu \ n oz-Tu \ n \'o n C., Palou s J., W \"u nsch R., 2013, @doi [ ] 10.1088/0004-637X/778/2/159 , http://cdsads.u-strasbg.fr/abs/2013ApJ...778..159T 778, 159

    work page doi:10.1088/0004-637x/778/2/159 2013

  56. [56]

    Tenorio-Tagle G., Mu \ n oz-Tu \ n \'o n C., Silich S., Cassisi S., 2015, @doi [ ] 10.1088/2041-8205/814/1/L8 , http://adsabs.harvard.edu/abs/2015ApJ...814L...8T 814, L8

    work page doi:10.1088/2041-8205/814/1/l8 2015

  57. [57]

    Terlevich R., Tenorio-Tagle G., Franco J., Melnick J., 1992, @doi [ ] 10.1093/mnras/255.4.713 , http://adsabs.harvard.edu/abs/1992MNRAS.255..713T 255, 713

    work page doi:10.1093/mnras/255.4.713 1992

  58. [58]

    Thornton K., Gaudlitz M., Janka H.-T., Steinmetz M., 1998, @doi [ ] 10.1086/305704 , http://adsabs.harvard.edu/abs/1998ApJ...500...95T 500, 95

    work page doi:10.1086/305704 1998

  59. [59]

    Todini P., Ferrara A., 2001, @doi [ ] 10.1046/j.1365-8711.2001.04486.x , http://adsabs.harvard.edu/abs/2001MNRAS.325..726T 325, 726

    work page doi:10.1046/j.1365-8711.2001.04486.x 2001

  60. [60]

    K., McKee C

    Truelove J. K., McKee C. F., 1999, @doi [ ] 10.1086/313176 , http://adsabs.harvard.edu/abs/1999ApJS..120..299T 120, 299

    work page doi:10.1086/313176 1999

  61. [61]

    Walch S., Naab T., 2015, @doi [ ] 10.1093/mnras/stv1155 , http://cdsads.u-strasbg.fr/abs/2015MNRAS.451.2757W 451, 2757

    work page doi:10.1093/mnras/stv1155 2015

  62. [62]

    Wiersma R. P. C., Schaye J., Smith B. D., 2009, @doi [ ] 10.1111/j.1365-2966.2008.14191.x , http://adsabs.harvard.edu/abs/2009MNRAS.393...99W 393, 99

    work page doi:10.1111/j.1365-2966.2008.14191.x 2009