pith. machine review for the scientific record. sign in

arxiv: 2604.19236 · v2 · submitted 2026-04-21 · 🌌 astro-ph.HE · astro-ph.SR

Recognition: unknown

Numerical Studies of Accretion Flows onto a Neutron Star Engulfed in a Massive Star

Daiyu Sakurai , Ryuichiro Akaho , Shoichi Yamada
Authors on Pith no claims yet

Pith reviewed 2026-05-10 02:37 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.SR
keywords common envelope evolutionneutron star accretionbow shocksdrag forcebinary evolutionhydrodynamical simulationsorbital decaymassive star envelope
0
0 comments X

The pith

Simulations show nested bow shocks around an engulfed neutron star enhance drag by 10-100 times and can reverse its direction.

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

The paper runs general-relativistic hydrodynamical simulations of a neutron star moving through the envelope of a massive star using a multi-layer domain decomposition to bridge the enormous scale gap from the star's surface to the accretion radius. It discovers that nested bow shocks with alternating orientations appear in most cases examined. This flow pattern produces a drag force much larger than the standard Bondi-Hoyle-Lyttleton estimate and allows the net force to point forward under some conditions rather than always backward. Accurate drag forces matter because they control the rate of orbital shrinkage during common envelope evolution and thus the properties of the resulting compact binaries.

Core claim

In a survey of ten models covering diverse envelope conditions, the accretion flow develops nested bow shocks whose orientations alternate. This arrangement differs from the usual single-shock picture and increases the drag force exerted on the neutron star by one to two orders of magnitude relative to the Bondi-Hoyle-Lyttleton formula. The net force can also reverse direction depending on the envelope density and velocity, contrary to the expectation that drag always slows the companion.

What carries the argument

nested bow shocks with alternating orientations in the accretion flow, which alter momentum transfer and produce enhanced or reversed net force on the neutron star

If this is right

  • Drag prescriptions in common envelope evolution models must be increased by factors of 10 to 100.
  • Orbital decay during the common envelope phase occurs much faster than in standard calculations.
  • The net force reversal means some neutron stars may not spiral inward under certain envelope conditions.
  • Predictions for merger timescales and final separations of neutron star binaries require revision.
  • The compactness of post-common-envelope systems can differ from earlier estimates.

Where Pith is reading between the lines

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

  • Force reversal in some envelopes could eject the neutron star outward, producing wider orbits than expected.
  • Stronger overall drag shortens the common envelope phase and may reduce the number of tight binaries that survive.
  • Including radiation or magnetic fields in future runs might break the nested shock pattern and lower the drag enhancement.
  • Comparing the revised drag rates to observed populations of post-common-envelope X-ray binaries offers a direct test.

Load-bearing premise

The multi-layer domain-decomposition approach resolves the accretion flow accurately over the 10,000 to 10 million scale range without introducing artifacts, and the modeled envelope conditions capture the essential physics without effects such as magnetic fields or radiation.

What would settle it

A higher-resolution simulation or one including magnetic fields that fails to produce alternating nested bow shocks and instead yields drag forces matching the Bondi-Hoyle-Lyttleton value would falsify the central claim.

Figures

Figures reproduced from arXiv: 2604.19236 by Daiyu Sakurai, Ryuichiro Akaho, Shoichi Yamada.

Figure 1
Figure 1. Figure 1: Schematic illustration of an NS engulfed in an RSG envelope during the spiral-in phase. The accretion radius 𝑅𝑎 gives the effective cross section of the NS, characterizing the gravitationally-focused flow near the NS. Since 𝑅𝑎 ≪ 𝑎 (orbital separation), the flow can be modelled locally by using the physical conditions of the RSG envelope as the boundary condition. which gives the spatial scale over which th… view at source ↗
Figure 2
Figure 2. Figure 2: presents the radial profiles of density 𝜌, Mach number M, dimensionless parameter 𝜖𝜌, and helium mass fraction for the RSG models with zero-age main-sequence masses of 25, 30, 35, and 40 𝑀⊙, evaluated at the onset of core helium burning. In what follows, we focus on the model with 𝑀ZAMS = 40 𝑀⊙, whose total mass at that stage is 𝑀RSG = 32.8 𝑀⊙. The density changes by ∼10 orders of magnitude. The NS moves a… view at source ↗
Figure 3
Figure 3. Figure 3: Radial profiles of the gravitational binding energy, 𝐸bind (𝑟 ), and the orbital energy budget, 𝛼CE |Δ𝐸orb (𝑟 )| , both employed in the 𝛼-formalism, for different choices of 𝛼th and 𝛼CE. accretion velocity so that the accretion radius and hence the com￾putational domain can be reduced (Cruz-Osorio & Rezzolla 2020; Kaaz et al. 2023). In addition, the common use of Newtonian gravity and idealized equations o… view at source ↗
Figure 4
Figure 4. Figure 4: Schematic illustration of the multi-layer domain-decomposition strategy (𝑛layer = 3) used to resolve the multiscale accretion flow. The com￾putational domain is divided into concentric spherical shells. The adjacent layers are partially overlapped (red regions). The hydrodynamical equations are solved in each layer individually and consecutively with the boundary conditions given by the solutions in the ad… view at source ↗
Figure 5
Figure 5. Figure 5: Quasi-steady structure for model test-Heos. Top row: Mach number and streamlines (left half) and density and velocity vectors (right half). Bottom row: adiabatic index (left half) and the time-variability diagnostic defined in Section 4.1 (right half). Panels from left to right zoom into the inner region by roughly a factor of ten, so that the shock fronts and time variability can be seen clearly. Numbered… view at source ↗
Figure 6
Figure 6. Figure 6: Quasi-steady structures in models test-gamma5/3 (top row) and test-gamma4/3 (bottom row). For each row, the left half of each panel displays the Mach number with streamlines, and the right half shows the density with velocity vectors as in [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Representative realistic models corresponding to different locations in the RSG envelope. Top, middle, and bottom rows correspond to 𝜌∞ = 1, 10−5 , and 10−8 g cm−3 . As in [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 9
Figure 9. Figure 9: Time evolution of the mass accretion rate 𝑀¤ (𝑅NS ) at the NS surface for model 𝜌1e0 with 𝜌∞ = 1 g cm−3 . The blue curve shows the instantaneous value; the red dashed line indicates the time-averaged value; the orange shaded band corresponds to the ±1𝜎 range [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Same as [PITH_FULL_IMAGE:figures/full_fig_p012_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Time evolution of the drag force at the NS surface for model 𝜌1e0 with 𝜌∞ = 1 g cm−3 . Colours denote different contributions (27 in Sec. 3.3): 𝐹mom,𝑧 (blue), 𝐹pre,𝑧 (red), 𝐹dyn,𝑧 (magenta), and the net force 𝐹net (cyan). The cyan dashed line and shaded band show the time-averaged value and ±1𝜎 dispersion of 𝐹net, respectively. MNRAS 000, 1–20 (2026) [PITH_FULL_IMAGE:figures/full_fig_p012_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Same as [PITH_FULL_IMAGE:figures/full_fig_p013_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Illustrative spiral-in trajectories of an NS inside an RSG envelope, starting from the same initial condition at 𝑎 = 0.95 RRSG ≃ 2300 𝑅⊙. The black curve is obtained by a numerical integration of Eqs. (30)–(33) with the canonical BHL formula for 40 yr; the red curve uses the time-averaged drag force inferred in this work and corresponds to 4 yr. Both trajectories start from the same initial condition. The… view at source ↗
Figure 14
Figure 14. Figure 14: shows the resulting detection reach for 𝜏 = 1 day. Among the three models considered here, only the innermost case, 𝜌1e0, corresponding to log(𝑎/𝑅⊙) = 0.06, reaches 10 kpc for HK, with a [PITH_FULL_IMAGE:figures/full_fig_p014_14.png] view at source ↗
read the original abstract

Massive stars commonly form binaries that can evolve into compact systems via common envelope evolution (CEE), a critical but poorly understood phase -- especially when the companion is a neutron star. Understanding the drag force exerted on a neutron star during CEE is a key to the quantitative evaluation of orbital decay, merger timescale, and compactness of the resultant binary. In this paper, we conduct general-relativistic hydrodynamical simulations under a novel strategy of multi-layer domain-decomposition to treat the vast disparity of $10^4$--$10^7$ between the neutron star radius and the accretion radius. Our 10-model survey spans diverse physical conditions that the neutron star encounters in the envelope of a massive star. We find that nested bow shocks with alternating orientations commonly form. This configuration is qualitatively different from those in the conventional picture and results in an enhancement of the drag force by one to two orders of magnitude from what the Bondi--Hoyle--Lyttleton formula predicts. Moreover, the direction of the net force can reverse depending on the envelope conditions, contrary to the standard picture in which the drag always decelerates the companion. These results will serve as a basis for improvements of the drag force prescription in CEE modeling, and have implications for binary evolution theory.

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

3 major / 2 minor

Summary. The paper conducts general-relativistic hydrodynamical simulations of accretion onto a neutron star engulfed in a massive star's envelope during common-envelope evolution. It employs a novel multi-layer domain-decomposition strategy to span the 10^4--10^7 dynamic range between the neutron-star radius and the accretion radius. A survey of 10 models under varied envelope conditions reports the formation of nested bow shocks with alternating orientations, yielding drag forces enhanced by 1--2 orders of magnitude relative to the Bondi--Hoyle--Lyttleton prediction and allowing net-force reversal depending on envelope conditions.

Significance. If the numerical results prove robust, the work would supply a revised, physically motivated drag prescription for CEE modeling with neutron-star companions, directly affecting orbital-decay rates, merger timescales, and the compactness of post-CEE binaries. The reported qualitative departure from the conventional BHL picture, including force reversal, is a potentially important shift. The scale-bridging method itself is a technical strength that could be adopted more broadly if validated.

major comments (3)
  1. [Numerical Method] Numerical method (multi-layer decomposition): The central claims rest on the multi-layer domain-decomposition technique, yet the manuscript supplies no explicit interface-flux continuity tests, global conservation diagnostics, or side-by-side comparisons against single-domain reference runs at overlapping scales. Without these, it remains possible that alternating bow shocks and the reported drag enhancement are numerical artifacts induced at layer boundaries across the 10^4--10^7 range.
  2. [Results] Results section (10-model survey): The abstract and survey description state that nested alternating bow shocks form commonly and boost drag by 1--2 orders of magnitude, but no resolution studies, convergence metrics, or error bars on the time-averaged drag force are provided. This absence prevents assessment of whether the quantitative enhancement and force-reversal behavior are converged and physical.
  3. [Results] Comparison to analytic limits: The manuscript contrasts its findings with the Bondi--Hoyle--Lyttleton formula but does not report direct numerical tests against the BHL solution in the appropriate supersonic, uniform-medium limit at the same resolution and domain setup. Such a benchmark is required to establish that the reported deviations are not method-specific.
minor comments (2)
  1. [Abstract] The abstract and introduction would benefit from a concise statement of the numerical resolution employed in each layer and the Courant-limited time-stepping strategy.
  2. [Figures] Figure captions describing shock structures should explicitly label the radial extents of each domain layer to allow readers to judge interface placement relative to the bow shocks.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their thorough review and constructive comments on our manuscript. We agree that further validation of the numerical method and additional quantitative assessments of the results will strengthen the paper. We address each major comment below and have prepared revisions accordingly.

read point-by-point responses

  1. Referee: Numerical method (multi-layer decomposition): The central claims rest on the multi-layer domain-decomposition technique, yet the manuscript supplies no explicit interface-flux continuity tests, global conservation diagnostics, or side-by-side comparisons against single-domain reference runs at overlapping scales. Without these, it remains possible that alternating bow shocks and the reported drag enhancement are numerical artifacts induced at layer boundaries across the 10^4--10^7 range.

    Authors: We acknowledge that explicit validation of the multi-layer domain-decomposition is important for supporting the central claims. Although the underlying GR hydrodynamics code has been validated in earlier studies, we agree that interface-specific tests should be shown here. In the revised manuscript we will add interface-flux continuity tests across layer boundaries, global conservation diagnostics (mass and energy), and side-by-side comparisons against single-domain reference runs at overlapping scales. Full single-domain runs across the entire 10^4--10^7 dynamic range remain computationally prohibitive, which motivated the multi-layer strategy, but the requested overlapping-scale benchmarks will be included to demonstrate continuity and rule out boundary artifacts. revision: yes

  2. Referee: Results section (10-model survey): The abstract and survey description state that nested alternating bow shocks form commonly and boost drag by 1--2 orders of magnitude, but no resolution studies, convergence metrics, or error bars on the time-averaged drag force are provided. This absence prevents assessment of whether the quantitative enhancement and force-reversal behavior are converged and physical.

    Authors: We agree that resolution studies, convergence metrics, and error estimates are required to substantiate the quantitative drag enhancements and force-reversal results. In the revised manuscript we will include resolution studies for representative models from the survey, showing how the time-averaged drag force behaves with increasing resolution. We will also report error bars or uncertainty ranges based on temporal variations during the averaging window and, where feasible, multiple realizations. revision: yes

  3. Referee: Comparison to analytic limits: The manuscript contrasts its findings with the Bondi--Hoyle--Lyttleton formula but does not report direct numerical tests against the BHL solution in the appropriate supersonic, uniform-medium limit at the same resolution and domain setup. Such a benchmark is required to establish that the reported deviations are not method-specific.

    Authors: We recognize the value of a direct numerical benchmark against the BHL analytic solution in the supersonic uniform-medium limit. In the revised manuscript we will add a dedicated benchmark subsection that performs such a simulation using the same resolution and domain configuration as the envelope models, and we will compare the measured drag force directly to the BHL prediction to confirm that the method recovers the expected result in that limit. revision: yes

Circularity Check

0 steps flagged

No circularity: results are direct numerical outputs

full rationale

The paper reports outcomes from general-relativistic hydrodynamical simulations employing a multi-layer domain-decomposition technique to handle scale disparities. No analytical derivation chain, parameter fitting, or predictive formulas are presented that could reduce to the inputs by construction. Claims about nested bow shocks, drag enhancement, and force reversal follow directly from the simulation survey across 10 models. No self-citations, ansatzes, or renamings of known results serve as load-bearing steps. The method is described as novel but is not justified via circular self-reference; any concerns about numerical artifacts fall under correctness or validation rather than circularity in the derivation.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The claims depend on the numerical method handling the scale disparity correctly and the physical assumptions about the envelope being sufficient.

axioms (2)
  • standard math The equations of general relativistic hydrodynamics accurately describe the accretion flow.
    Basis for the simulations described.
  • domain assumption The envelope of the massive star can be modeled with simplified density and velocity profiles spanning diverse conditions.
    The 10 models use various physical conditions but specific profiles are assumed.

pith-pipeline@v0.9.0 · 5534 in / 1421 out tokens · 55093 ms · 2026-05-10T02:37:03.456225+00:00 · methodology

discussion (0)

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

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Hyperaccreting Neutron Stars inside Massive Envelopes: The Implausibility of Thorne-\.Zytkow Objects

    astro-ph.HE 2026-04 unverdicted novelty 7.0

    Hypercritical accretion onto neutron stars embedded in massive envelopes leads to rapid collapse into black holes rather than stable Thorne-Zytkow objects.

Reference graph

Works this paper leans on

74 extracted references · 61 canonical work pages · cited by 1 Pith paper · 1 internal anchor

  1. [1]

    1939", month =

    The collision of two stars. , year = "1939", month = "Apr", volume =

    1939

  2. [2]

    , year = 1944, month = jan, volume =

    On the mechanism of accretion by stars. , year = 1944, month = jan, volume =. doi:10.1093/mnras/104.5.273 , adsurl =

    work page doi:10.1093/mnras/104.5.273 1944

  3. [3]

    1952", month =

    On Spherically Symmetrical Accretion. , year = "1952", month = "May", volume =

    1952

  4. [4]

    Structure and Evolution of Close Binary Systems , year = 1976, editor =

    Common Envelope Binaries. Structure and Evolution of Close Binary Systems , year = 1976, editor =

    1976

  5. [5]

    , keywords =

    Double white dwarfs as progenitors of R Coronae Borealis stars and type I supernovae. , year = 1984, month = feb, volume =. doi:10.1086/161701 , adsurl =

    work page doi:10.1086/161701 1984

  6. [6]

    , year = 1985, month = nov, volume =

    Hydrodynamic calculations of axisymmetric accretion flow. , year = 1985, month = nov, volume =. doi:10.1093/mnras/217.2.367 , adsurl =

    work page doi:10.1093/mnras/217.2.367 1985

  7. [7]

    , keywords =

    Common Envelope Evolution and Double Cores of Planetary Nebulae. , year = 1990, month = jul, volume =. doi:10.1086/168974 , adsurl =

    work page doi:10.1086/168974 1990

  8. [8]

    C., & Chevalier, R

    Steady Spherical Hypercritical Accretion onto Neutron Stars. , year = 1991, month = jul, volume =. doi:10.1086/170272 , adsurl =

    work page doi:10.1086/170272 1991

  9. [9]

    , keywords =

    Neutron Star Accretion in a Stellar Envelope. , keywords =. doi:10.1086/186905 , adsurl =

    work page doi:10.1086/186905

  10. [10]

    Three-dimensional hydrodynamic Bondi-Hoyle accretion. V. Specific heat ratio 1.01, nearly isothermal flow. , keywords =

  11. [11]

    , keywords =

    Dynamical Friction in a Gaseous Medium. , keywords =. doi:10.1086/306858 , archivePrefix =. astro-ph/9810324 , adsurl =

    work page doi:10.1086/306858

  12. [12]

    , year = 2000, month = aug, volume =

    On the energy equation and efficiency parameter of the common envelope evolution. , year = 2000, month = aug, volume =

    2000

  13. [13]

    , keywords =

    Common Envelope Evolution of Massive Binary Stars. , keywords =. doi:10.1146/annurev.astro.38.1.113 , adsurl =

    work page doi:10.1146/annurev.astro.38.1.113

  14. [14]

    X., & Swesty, F

    The Accuracy, Consistency, and Speed of an Electron-Positron Equation of State Based on Table Interpolation of the Helmholtz Free Energy. , keywords =. doi:10.1086/313304 , adsurl =

    work page doi:10.1086/313304

  15. [15]

    , keywords =

    A review of Bondi-Hoyle-Lyttleton accretion. , keywords =. doi:10.1016/j.newar.2004.06.001 , archivePrefix =. astro-ph/0406166 , primaryClass =

    work page doi:10.1016/j.newar.2004.06.001 2004

  16. [16]

    , keywords =

    A fresh look at the unstable simulations of Bondi-Hoyle-Lyttleton accretion. , keywords =. doi:10.1051/0004-6361:20042201 , archivePrefix =. astro-ph/0502168 , primaryClass =

    work page doi:10.1051/0004-6361:20042201

  17. [17]

    , year = 2007, month = aug, volume =

    Dynamical Friction of a Circular-Orbit Perturber in a Gaseous Medium. , year = 2007, month = aug, volume =. doi:10.1086/519302 , archivePrefix =. 0705.0084 , adsurl =

    work page doi:10.1086/519302 2007

  18. [18]

    , keywords =

    The Interaction of Stellar Objects within a Common Envelope. , keywords =. doi:10.1086/526343 , archivePrefix =. 0710.3631 , primaryClass =

    work page doi:10.1086/526343

  19. [19]

    , keywords =

    Simulating the Common Envelope Phase of a Red Giant Using Smoothed-particle Hydrodynamics and Uniform-grid Codes. , keywords =. doi:10.1088/0004-637X/744/1/52 , archivePrefix =. 1107.5072 , primaryClass =

    work page doi:10.1088/0004-637x/744/1/52

  20. [20]

    , keywords =

    An AMR Study of the Common-envelope Phase of Binary Evolution. , keywords =. doi:10.1088/0004-637X/746/1/74 , archivePrefix =. 1107.3889 , primaryClass =

    work page doi:10.1088/0004-637x/746/1/74

  21. [21]

    Double Compact Objects. I. The Significance of the Common Envelope on Merger Rates. , keywords =. doi:10.1088/0004-637X/759/1/52 , archivePrefix =. 1202.4901 , primaryClass =

    work page doi:10.1088/0004-637x/759/1/52

  22. [22]

    E., de Koter, A., et al

    Binary Interaction Dominates the Evolution of Massive Stars. Science , keywords =. doi:10.1126/science.1223344 , archivePrefix =. 1207.6397 , primaryClass =

    work page doi:10.1126/science.1223344

  23. [23]

    2013, A&A Rv, 21, 59, doi: 10.1007/s00159-013-0059-2

    Common envelope evolution: where we stand and how we can move forward. , keywords =. doi:10.1007/s00159-013-0059-2 , archivePrefix =. 1209.4302 , primaryClass =

    work page doi:10.1007/s00159-013-0059-2

  24. [24]

    A., & Yungelson, L

    The Evolution of Compact Binary Star Systems. Living Reviews in Relativity , keywords =. doi:10.12942/lrr-2014-3 , archivePrefix =. 1403.4754 , primaryClass =

    work page doi:10.12942/lrr-2014-3 2014

  25. [25]

    , keywords =

    Asymmetric Accretion Flows within a Common Envelope. , keywords =. doi:10.1088/0004-637X/803/1/41 , archivePrefix =. 1410.3823 , primaryClass =

    work page doi:10.1088/0004-637x/803/1/41

  26. [26]

    , keywords =

    Recombination Energy in Double White Dwarf Formation. , year = 2015, month = jun, volume =. doi:10.1093/mnrasl/slv043 , archivePrefix =. 1503.02750 , adsurl =

    work page doi:10.1093/mnrasl/slv043 2015

  27. [27]

    T., R¨ opke, F

    Hydrodynamic Moving-mesh Simulations of the Common Envelope Phase in Binary Stellar Systems. , year = 2016, month = jan, volume =. doi:10.3847/2041-8205/816/1/L9 , archivePrefix =. 1512.04529 , adsurl =

    work page doi:10.3847/2041-8205/816/1/l9 2016

  28. [28]

    , archivePrefix = "arXiv", eprint =

    The first gravitational-wave source from the isolated evolution of two stars in the 40-100 solar mass range. Nature , year = 2016, month = jun, volume =. doi:10.1038/nature18322 , archivePrefix =. 1602.04531 , adsurl =

    work page doi:10.1038/nature18322 2016

  29. [29]

    doi:10.1142/9692 , adsurl =

    Numerical Relativity. doi:10.1142/9692 , adsurl =

    work page doi:10.1142/9692

  30. [30]

    , keywords =

    Common Envelope Wind Tunnel: Coefficients of Drag and Accretion in a Simplified Context for Studying Flows around Objects Embedded within Stellar Envelopes. , keywords =. doi:10.3847/1538-4357/aa6117 , archivePrefix =. 1704.02372 , primaryClass =

    work page doi:10.3847/1538-4357/aa6117

  31. [31]

    The Lives and Death-Throes of Massive Stars , year = 2017, editor =

    The multiplicity of massive stars: a 2016 view. The Lives and Death-Throes of Massive Stars , year = 2017, editor =. doi:10.1017/S1743921317003209 , archivePrefix =. 1703.01608 , primaryClass =

    work page doi:10.1017/s1743921317003209 2016

  32. [32]

    , keywords =

    The effect of binding energy and resolution in simulations of the common envelope binary interaction. , keywords =. doi:10.1093/mnras/sty794 , archivePrefix =. 1706.09786 , primaryClass =

    work page doi:10.1093/mnras/sty794

  33. [33]

    , year = 2018, month = dec, volume =

    On the formation history of Galactic double neutron stars. , year = 2018, month = dec, volume =. doi:10.1093/mnras/sty2463 , archivePrefix =. 1805.07974 , adsurl =

    work page doi:10.1093/mnras/sty2463 2018

  34. [34]

    J., Ramirez-Ruiz, E., et al

    The Complete Evolution of a Neutron-Star Binary through a Common Envelope Phase Using 1D Hydrodynamic Simulations. , keywords =. doi:10.3847/2041-8213/ab40d1 , archivePrefix =. 1907.12573 , primaryClass =

    work page doi:10.3847/2041-8213/ab40d1 2041

  35. [35]

    2019, ApJS, 243, 10, doi: 10.3847/1538-4365/ab2241

    Modules for Experiments in Stellar Astrophysics (MESA): Pulsating Variable Stars, Rotating Stars, and Applied Numerical Methods. , keywords =. doi:10.3847/1538-4365/ab2241 , archivePrefix =. 1903.01426 , primaryClass =

    work page doi:10.3847/1538-4365/ab2241 1903

  36. [36]

    , keywords =

    How drag force evolves in global common envelope simulations. , keywords =. doi:10.1093/mnras/stz2813 , archivePrefix =. 1908.06195 , primaryClass =

    work page doi:10.1093/mnras/stz2813 1908

  37. [37]

    , keywords =

    Common-envelope evolution with an asymptotic giant branch star. , year = 2020, month = dec, volume =. doi:10.1051/0004-6361/202038992 , archivePrefix =. 2007.11000 , adsurl =

    work page doi:10.1051/0004-6361/202038992 2020

  38. [38]

    , keywords =

    Common-envelope Dynamics of a Stellar-mass Black Hole: General Relativistic Simulations. , keywords =. doi:10.3847/1538-4357/ab89aa , archivePrefix =. 2004.13782 , primaryClass =

    work page doi:10.3847/1538-4357/ab89aa 2004

  39. [39]

    , keywords =

    Common Envelope Wind Tunnel: The Effects of Binary Mass Ratio and Implications for the Accretion-driven Growth of LIGO Binary Black Holes. , keywords =. doi:10.3847/1538-4357/ab9ac6 , archivePrefix =. 1910.13333 , primaryClass =

    work page doi:10.3847/1538-4357/ab9ac6 1910

  40. [40]

    2020, Common Envelope Evolution, doi: 10.1088/2514-3433/abb6f0

    Common Envelope Evolution. doi:10.1088/2514-3433/abb6f0 , adsurl =

    work page doi:10.1088/2514-3433/abb6f0

  41. [41]

    , keywords =

    Impact of massive binary star and cosmic evolution on gravitational wave observations I: black hole-neutron star mergers. , keywords =. doi:10.1093/mnras/stab2716 , archivePrefix =. 2103.02608 , primaryClass =

    work page doi:10.1093/mnras/stab2716

  42. [42]

    , year = 2022, month = apr, volume =

    Bipolar planetary nebulae from common-envelope evolution of binary stars. , year = 2022, month = apr, volume =. doi:10.1051/0004-6361/202142478 , archivePrefix =. 2110.13177 , adsurl =

    work page doi:10.1051/0004-6361/202142478 2022

  43. [43]

    M., Schneider, F

    From 3D hydrodynamic simulations of common-envelope interaction to gravitational-wave mergers. , keywords =. doi:10.1051/0004-6361/202142731 , archivePrefix =. 2111.12112 , primaryClass =

    work page doi:10.1051/0004-6361/202142731

  44. [44]

    2022, ApJL, 937, L42, doi: 10.3847/2041-8213/ac9519

    A Two-stage Formalism for Common-envelope Phases of Massive Stars. , year = 2022, month = oct, volume =. doi:10.3847/2041-8213/ac9519 , archivePrefix =. 2207.07140 , adsurl =

    work page doi:10.3847/2041-8213/ac9519 2022

  45. [45]

    , keywords =

    Common envelopes in massive stars: towards the role of radiation pressure and recombination energy in ejecting red supergiant envelopes. , keywords =. doi:10.1093/mnras/stac049 , archivePrefix =. 2111.00923 , primaryClass =

    work page doi:10.1093/mnras/stac049

  46. [46]

    , keywords =

    Three-dimensional simulations of the jet feedback mechanism in common envelope jets supernovae. , keywords =. doi:10.1093/mnras/stac1341 , archivePrefix =. 2112.01459 , primaryClass =

    work page doi:10.1093/mnras/stac1341

  47. [47]

    2019, MNRAS, 488, 5615, doi: 10.1093/mnras/stz2013

    Companion-launched jets and their effect on the dynamics of common envelope interaction simulations. , keywords =. doi:10.1093/mnras/stz2013 , archivePrefix =. 1902.03931 , primaryClass =

    work page doi:10.1093/mnras/stz2013 1902

  48. [48]

    M., Tomida, K., White, C

    The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers. , keywords =. doi:10.3847/1538-4365/ab929b , archivePrefix =. 2005.06651 , primaryClass =

    work page doi:10.3847/1538-4365/ab929b 2005

  49. [49]

    S., Berger, E., Stevenson, S., et al

    Impact of massive binary star and cosmic evolution on gravitational wave observations - II. Double compact object rates and properties. , keywords =. doi:10.1093/mnras/stac1677 , archivePrefix =. 2112.05763 , primaryClass =

    work page doi:10.1093/mnras/stac1677

  50. [50]

    2022, PhR, 955, 1, doi: 10.1016/j.physrep.2022.01.003

    Merging stellar-mass binary black holes. , keywords =. doi:10.1016/j.physrep.2022.01.003 , archivePrefix =. 1806.05820 , primaryClass =

    work page doi:10.1016/j.physrep.2022.01.003 2022

  51. [51]

    , keywords =

    Protoneutron Star Convection Simulated with a New General Relativistic Boltzmann Neutrino Radiation Hydrodynamics Code. , keywords =. doi:10.3847/1538-4357/acad76 , archivePrefix =. 2206.01673 , primaryClass =

    work page doi:10.3847/1538-4357/acad76

  52. [52]

    Living Reviews in Computational Astrophysics , keywords =

    Simulations of common-envelope evolution in binary stellar systems: physical models and numerical techniques. Living Reviews in Computational Astrophysics , keywords =. doi:10.1007/s41115-023-00017-x , archivePrefix =. 2212.07308 , primaryClass =

    work page doi:10.1007/s41115-023-00017-x

  53. [53]

    D., Acernese, F., et al

    GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo during the Second Part of the Third Observing Run. Physical Review X , keywords =. doi:10.1103/PhysRevX.13.041039 , archivePrefix =. 2111.03606 , primaryClass =

    work page doi:10.1103/physrevx.13.041039

  54. [54]

    2024, ARA&A, 62, 21, doi: 10.1146/annurev-astro-052722-105936

    The Evolution of Massive Binary Stars. , keywords =. doi:10.1146/annurev-astro-052722-105936 , archivePrefix =. 2311.01865 , primaryClass =

    work page doi:10.1146/annurev-astro-052722-105936

  55. [55]

    2024, ApJ, 977, 200, doi: 10.3847/1538-4357/ad9561

    Multizone Modeling of Black Hole Accretion and Feedback in 3D GRMHD: Bridging Vast Spatial and Temporal Scales. , keywords =. doi:10.3847/1538-4357/ad9561 , archivePrefix =. 2405.13887 , primaryClass =

    work page doi:10.3847/1538-4357/ad9561

  56. [56]

    Nature Astronomy , keywords =

    Observations of high-order multiplicity in a high-mass stellar protocluster. Nature Astronomy , keywords =. doi:10.1038/s41550-023-02181-9 , archivePrefix =. 2401.06545 , primaryClass =

    work page doi:10.1038/s41550-023-02181-9

  57. [57]

    , keywords =

    Constraining common envelope evolution in binary neutron star formation with combined galactic and gravitational-wave observations. , keywords =. doi:10.1093/mnrasl/slaf105 , archivePrefix =. 2508.14397 , primaryClass =

    work page doi:10.1093/mnrasl/slaf105

  58. [58]

    GWTC-4.0: Population Properties of Merging Compact Binaries

    GWTC-4.0: Population Properties of Merging Compact Binaries. arXiv e-prints , keywords =. doi:10.48550/arXiv.2508.18083 , archivePrefix =. 2508.18083 , primaryClass =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2508.18083

  59. [59]

    Impact of rotation, accretion, and stratification

    Local simulations of common-envelope dynamical inspiral. Impact of rotation, accretion, and stratification. arXiv e-prints , keywords =. doi:10.48550/arXiv.2601.04188 , archivePrefix =. 2601.04188 , primaryClass =

    work page doi:10.48550/arxiv.2601.04188

  60. [60]

    Z ytkow Objects through Neutrinos

    Identifying Thorne-\. Z ytkow Objects through Neutrinos. , year = 2025, month = apr, volume =

    2025

  61. [61]

    , year = 2026, month = jan, volume =

    r -process Nucleosynthesis from Hyperaccreting Neutron Stars in Common Envelopes. , year = 2026, month = jan, volume =. doi:10.3847/1538-4357/ae2311 , adsurl =

    work page doi:10.3847/1538-4357/ae2311 2026

  62. [62]

    , year = 2021, month = oct, volume =

    Prospects of Gravitational Wave Detections from Common Envelope Evolution with LISA. , year = 2021, month = oct, volume =. doi:10.3847/1538-4357/ac1110 , adsurl =

    work page doi:10.3847/1538-4357/ac1110 2021

  63. [63]

    , year = 2020, month = apr, volume =

    Gravitational waves from inspirals of compact objects in binary common-envelope evolution. , year = 2020, month = apr, volume =

    2020

  64. [64]

    Neutrino Energy Loss in Stellar Interiors. VII. Pair, Photo-, Plasma, Bremsstrahlung, and Recombination Neutrino Processes. , year = 1996, month = feb, volume =. doi:10.1086/192264 , adsurl =

    work page doi:10.1086/192264 1996

  65. [65]

    Diffuse supernova neutrino background search at Super-Kamiokande. Phys. Rev. D , year = 2021, month = dec, volume =

    2021

  66. [66]

    JCAP , year = 2022, month = oct, volume =

    Prospects for Detecting the Diffuse Supernova Neutrino Background with JUNO. JCAP , year = 2022, month = oct, volume =

    2022

  67. [67]

    , title =

    Lions, P.-L. , title =. First International Symposium on Domain Decomposition Methods for Partial Differential Equations , editor =

  68. [68]

    and Valli, A

    Quarteroni, A. and Valli, A. , title =. 1999 , doi =

    1999

  69. [69]

    Schwarz, H. A. , title =. Vierteljahrsschrift der Naturforschenden Gesellschaft in Z

  70. [70]

    and Widlund, O

    Toselli, A. and Widlund, O. B. , title =. 2005 , doi =

    2005

  71. [71]

    , keywords =

    Detectability of Late-time Supernova Neutrinos with Fallback Accretion onto Protoneutron Star. , keywords =. doi:10.3847/1538-4357/ad118c , archivePrefix =. 2304.11150 , primaryClass =

    work page doi:10.3847/1538-4357/ad118c

  72. [72]

    doi:10.3847/1538-4357/acc7a1

    Jet Formation in 3D GRMHD Simulations of Bondi-Hoyle-Lyttleton Accretion. , keywords =. doi:10.3847/1538-4357/acc7a1 , archivePrefix =. 2201.11753 , primaryClass =

    work page doi:10.3847/1538-4357/acc7a1

  73. [73]

    , keywords =

    Three-dimensional radiation hydrodynamics simulations of wandering intermediate-mass black holes considering the anisotropic radiation and dust sublimation. , keywords =. doi:10.1093/mnras/stae195 , archivePrefix =. 2402.00369 , primaryClass =

    work page doi:10.1093/mnras/stae195

  74. [74]

    , keywords =

    Gaseous dynamical friction under radiative feedback: do intermediate-mass black holes speed up or down?. , keywords =. doi:10.1093/mnras/staa1338 , archivePrefix =. 2002.08017 , primaryClass =

    work page doi:10.1093/mnras/staa1338 2002