pith. sign in

arxiv: 1906.10630 · v1 · pith:EYNKDTDKnew · submitted 2019-06-25 · 🌌 astro-ph.HE

Timing Properties of Shocked Accretion Flows around Neutron Stars -- II. Viscous Disks and Boundary Layers

Ayan Bhattacharjee , Sandip K. Chakrabarti This is my paper

Pith reviewed 2026-05-25 15:53 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords accretion flowsneutron starsboundary layersviscositySPH simulationsCENBOLoutflowsHMXBs
0
0 comments X

The pith

Viscous sub-Keplerian flows around neutron stars form multiple boundary layers whose structure changes with viscosity strength, always including a centrifugal barrier layer.

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

The paper simulates viscous accretion flows onto weakly magnetic neutron stars using smoothed particle hydrodynamics and includes radiative cooling. It reports that low viscosity produces only a Normal Boundary Layer to match rotation, moderate viscosity adds a Radiative Keplerian Disk layer between the incoming flow and that boundary layer, and higher viscosity leaves only those two layers while a centrifugal-pressure-dominated layer always appears farther out. The work establishes the first self-consistent transition from sub-Keplerian inflow to these layered structures for neutron stars and links the layers to matter ejection. The results apply directly to wind-fed systems such as high-mass X-ray binaries.

Core claim

With the introduction of a small viscosity in a sub-Keplerian flow, only a single Normal Boundary Layer forms to adjust the rotational velocity component. With the increase of viscosity, the region extends radially and beyond some critical value, a RAdiative KEplerian Disk forms between the sub-Keplerian flow and the NBOL. When viscosity is increased further only NBOL and RAKED remain. In all such cases, the CENtrifugal pressure dominated BOundary Layer is formed away from the star. This is the first self-consistent study where such a transition from sub-Keplerian flows has been reported for neutron stars, and the results identify the connection between accretion and ejection of matter.

What carries the argument

Formation of multiple boundary layers (NBOL, RAKED, CENBOL) in the presence of viscosity and cooling, tracked through Smoothed Particle Hydrodynamics simulations of sub-Keplerian inflow.

If this is right

  • A single Normal Boundary Layer appears at low viscosity to adjust rotation.
  • A Radiative Keplerian Disk layer appears between the sub-Keplerian flow and the Normal Boundary Layer once viscosity exceeds a threshold.
  • The centrifugal-pressure-dominated boundary layer forms away from the star in every viscosity regime.
  • The layered structure produces a direct connection between accretion and ejection of matter.
  • The same sequence governs disk, boundary-layer, and outflow formation in wind-dominated neutron-star systems.

Where Pith is reading between the lines

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

  • The viscosity-dependent layering may set the characteristic frequencies observed in the power spectra of neutron-star X-ray binaries.
  • If real microphysical viscosity behaves similarly, the same transitions could appear in other compact-object accretion flows when cooling is efficient.
  • The model implies that changes in effective viscosity could switch the system between different outflow states without altering the central object.

Load-bearing premise

The simulations assume a specific prescription for viscosity and radiative cooling that produces the reported layer transitions.

What would settle it

X-ray timing observations of a wind-accreting neutron star system such as Cir X-1 that either show or fail to show the predicted sequence of radial layer transitions as accretion rate or effective viscosity changes.

Figures

Figures reproduced from arXiv: 1906.10630 by Ayan Bhattacharjee, Sandip K. Chakrabarti.

Figure 1
Figure 1. Figure 1: (a) Velocity vector vrrˆ+ vzzˆ (arrow heads) with Mach number in the colour bar, for the flow configuration C1 at time t = 0.2800s and (b) corresponding contours of constant Mach number range showing the outer and inner shocks. (c) Velocity vector vrrˆ + vzzˆ (arrow heads) with Mach number in the colour bar, for the flow configuration C2 at time t = 0.1400s and (d) corresponding contours of constant Mach n… view at source ↗
read the original abstract

We use Smoothed Particle Hydrodynamics to study viscous accretion flows around a weakly magnetic neutron star. We show the formation of multiple ``boundary" layers in presence of both cooling and viscosity. We find that with the introduction of a small viscosity in a sub-Keplerian flow, much like the wind accretion in HMXBs such as Cir X-1, only a single Normal Boundary Layer (NBOL) forms to adjust the rotational velocity component. With the increase of viscosity, the region extends radially and beyond some critical value, a RAdiative KEplerian Disk/layer (RAKED) forms between the sub-Keplerian flow and the NBOL. When viscosity is increased further only NBOL and RAKED remain. In all such cases, the CENtrifugal pressure dominated BOundary Layer (CENBOL) is formed, away from the star, as in the case of black holes. This is the first self-consistent study where such a transition from sub-Keplerian flows has been reported for neutron stars. We also identify the connection between accretion and ejection of matter, following the Two-Component Advective Flow for black holes, for neutron stars. The results are crucial in the understanding of the formation of disks, boundary layers and outflows in wind dominated neutron star systems.

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 uses Smoothed Particle Hydrodynamics simulations to examine viscous accretion flows with radiative cooling onto weakly magnetic neutron stars, initialized in sub-Keplerian rotation. It reports that increasing the viscosity parameter produces a sequence of boundary-layer structures: a single Normal Boundary Layer (NBOL) at low viscosity, followed by the additional formation of a Radiative Keplerian Disk/layer (RAKED) at intermediate values, with only NBOL and RAKED remaining at higher viscosity; a Centrifugal pressure dominated Boundary Layer (CENBOL) forms away from the star in all cases. The work claims this constitutes the first self-consistent demonstration of such transitions for neutron stars and identifies links between accretion and ejection analogous to the Two-Component Advective Flow model for black holes.

Significance. If the reported structures prove robust under the adopted viscosity and cooling prescriptions, the results would provide a useful descriptive mapping of how sub-Keplerian flows adjust to form multiple boundary layers around neutron stars. This extends boundary-layer phenomenology previously explored for black holes to neutron-star systems and supplies a framework for interpreting disk formation and outflows in wind-fed high-mass X-ray binaries. The approach relies on forward integration of the hydrodynamic equations with added terms rather than any fitted or circular construction.

major comments (2)
  1. [Abstract] Abstract: the central claims rest on the emergence of specific layer sequences (NBOL alone, then NBOL+RAKED, then NBOL+RAKED) as viscosity is increased, yet the manuscript supplies no information on numerical resolution, convergence tests, error bars on layer locations, or sensitivity to the viscosity and cooling parameters. Without these, the quantitative reliability of the reported transitions cannot be assessed.
  2. [Results] Simulation description (throughout results): the transitions are stated to depend on viscosity strength, but no tests are presented showing that the layer locations or the critical viscosity values remain stable under changes in particle number, artificial viscosity coefficients, or cooling rate. This is load-bearing because the sequence itself is the primary result.
minor comments (1)
  1. [Abstract] Abstract: the acronym RAKED is defined with mixed capitalization ('RAdiative KEplerian Disk/layer'); a consistent typographic convention would improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments. The concerns about numerical resolution, convergence, and parameter sensitivity are valid and will be addressed by expanding the methods and results sections in the revised manuscript.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the central claims rest on the emergence of specific layer sequences (NBOL alone, then NBOL+RAKED, then NBOL+RAKED) as viscosity is increased, yet the manuscript supplies no information on numerical resolution, convergence tests, error bars on layer locations, or sensitivity to the viscosity and cooling parameters. Without these, the quantitative reliability of the reported transitions cannot be assessed.

    Authors: We agree that the manuscript does not currently supply explicit details on numerical resolution, convergence tests, error bars on layer locations, or sensitivity analyses. In the revised version we will add a dedicated paragraph in the numerical methods section specifying the SPH particle number employed, the artificial viscosity coefficients, and any resolution or parameter-sensitivity checks that were performed. Where feasible we will also report approximate uncertainties on the radial locations of the identified layers and discuss how the cooling prescription affects the transitions. revision: yes

  2. Referee: [Results] Simulation description (throughout results): the transitions are stated to depend on viscosity strength, but no tests are presented showing that the layer locations or the critical viscosity values remain stable under changes in particle number, artificial viscosity coefficients, or cooling rate. This is load-bearing because the sequence itself is the primary result.

    Authors: We concur that demonstrating stability of the layer sequence and critical viscosity thresholds under changes in particle number, artificial viscosity, and cooling rate is necessary. We will perform and document additional runs that vary these quantities and include the outcomes in the revised results section to show that the reported progression (NBOL alone, NBOL+RAKED, NBOL+RAKED with CENBOL always present) persists. This will directly address the robustness concern. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results from forward SPH integration

full rationale

The paper reports outcomes of Smoothed Particle Hydrodynamics simulations of sub-Keplerian flows with added viscosity and cooling terms. The reported boundary-layer transitions (NBOL, RAKED, CENBOL) are direct numerical results under the chosen prescriptions, not quantities fitted from the same data or defined in terms of themselves. No equation or step reduces the layer formation to a self-citation chain, ansatz smuggled via citation, or renaming of a known result. The central claim is a descriptive statement about what the simulations exhibit, which remains independent of the input parameters. Minor self-citations to prior TCAF work by the same group are present but not load-bearing for the simulation outcomes themselves.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The model rests on standard hydrodynamic equations plus domain assumptions about viscosity and cooling; no new particles or forces are introduced, but the specific functional forms of viscosity and cooling act as adjustable elements whose values determine the reported transitions.

free parameters (2)
  • viscosity strength parameter
    Viscosity is varied across regimes to produce the sequence of layer formations; its functional form and magnitude are not fixed by first principles in the abstract.
  • cooling rate parameter
    Cooling is included alongside viscosity and is required for the radiative Keplerian layer to appear.
axioms (2)
  • standard math Standard Navier-Stokes equations govern the accretion flow under the SPH discretization
    SPH is used to integrate the hydrodynamic equations with added viscosity and cooling terms.
  • domain assumption The neutron star is weakly magnetic so that magnetic effects can be neglected
    Stated in the title and abstract as the regime under study.

pith-pipeline@v0.9.0 · 5779 in / 1379 out tokens · 29104 ms · 2026-05-25T15:53:30.034874+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.

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    We use Smoothed Particle Hydrodynamics to study viscous accretion flows around a weakly magnetic neutron star. We show the formation of multiple “boundary” layers in presence of both cooling and viscosity... with the introduction of a small viscosity... only a single Normal Boundary Layer (NBOL) forms... With the increase of viscosity... a RAdiative KEplerian Disk/layer (RAKED) forms...

  • IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    We increased the viscosity parameter from 0.075 (C1) to 0.15 (C2) to 0.3 (C3) and kept injected λinj = 1.7

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

46 extracted references · 46 canonical work pages · 2 internal anchors

  1. [1]

    - [1] #1 = = ^ ^ ^ .\!\!^ d .\!\!^ h .\!\!^ m .\!\!^ s .\!\!^ @mss

    thebibliography [1] 20pt to REFERENCES 6pt =0pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key on reference command Each re...

    work page doi:10.5281/zenodo.831784 2017

  2. [2]

    K., & Molteni, D., 2002, Journal of Astrophysics and Astronomy, 23, 155

    Acharya, K., Chakrabarti, S. K., & Molteni, D., 2002, Journal of Astrophysics and Astronomy, 23, 155

    work page 2002

  3. [3]

    K.\ 2019, arXiv e-prints , arXiv:1905.01538

    Banerjee, A., Bhattacharjee, A., & Chakrabarti, S. K.\ 2019, arXiv e-prints , arXiv:1905.01538

    work page arXiv 2019

  4. [4]

    Banerjee, I., Bhattacharjee, A., Banerjee, A., et al.\ 2019, arXiv e-prints , arXiv:1904.11644

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  5. [5]

    F., Paltani, S., & Skinner, G

    Barret, D., Klu \'z niak, W., Olive, J. F., Paltani, S., & Skinner, G. K., 2005, MNRAS, 357, 1288

    work page 2005

  6. [6]

    C., 2005, MNRAS, 361, 855

    Barret, D., Olive, J.-F., & Miller, M. C., 2005, MNRAS, 361, 855

    work page 2005

  7. [7]

    (Eds.), ASSP, 53, 93, Springer (Heidelberg)

    Bhattacharjee, A., 2018, in Exploring the Universe: From Near Space to Extra-galactic Mukhopadhyay, B & Sasmal, S. (Eds.), ASSP, 53, 93, Springer (Heidelberg)

    work page 2018

  8. [8]

    K., 2017, MNRAS, 466, 1372

    Bhattacharjee, A., Banerjee, I., Banerjee, A., Debnath, D., & Chakrabarti, S. K., 2017, MNRAS, 466, 1372

    work page 2017

  9. [9]

    K., 2017, MNRAS, 472, 1361

    Bhattacharjee, A., & Chakrabarti, S. K., 2017, MNRAS, 472, 1361

    work page 2017

  10. [10]

    K.\ 2019, ApJ, 873, 119

    Bhattacharjee, A., & Chakrabarti, S. K.\ 2019, ApJ, 873, 119

    work page 2019

  11. [11]

    Bhattacharyya, S., & Chakrabarty, D.\ 2017, , 835, 4

    work page 2017

  12. [12]

    Belloni, T., Psaltis, D., & van der Klis, M., 2002, ApJ, 572, 392

    work page 2002

  13. [13]

    Bildsten, L.\ 1998, American Institute of Physics Conference Series, 431, 299

    work page 1998

  14. [14]

    Boutloukos, S., van der Klis, M., Altamirano, D., et al., 2006, ApJ, 653, 1435

    work page 2006

  15. [15]

    K., 1989, MNRAS, 240, 7

    Chakrabarti, S. K., 1989, MNRAS, 240, 7

    work page 1989

  16. [16]

    K., 1990, MNRAS, 243, 610

    Chakrabarti, S. K., 1990, MNRAS, 243, 610

    work page 1990

  17. [17]

    Chakrabarti, 1995, in Ann. NY Acad. Sci., 759, Seventeenth Texas Symposium on Relativistic Astrophysics and Cosmology, ed. H. Bohringer, G.E. Morfil & J. Truemper, 546

    work page 1995

  18. [18]

    K., 1996, ApJ, 464, 664

    Chakrabarti, S. K., 1996, ApJ, 464, 664

    work page 1996

  19. [19]

    K., 1997, ApJ, 484, 313

    Chakrabarti, S. K., 1997, ApJ, 484, 313

    work page 1997

  20. [20]

    Study of Accretion processes Around Black Holes becomes Science: Tell Tale Observational Signatures of Two Component Advective Flows

    Chakrabarti, S. K., 2017, Proceedings of 14th Marcel Grossman meeting on General relativity, M. Bianchi, R. Ruffini, R.Jantzen (Eds.), 369-384 (World Scientific:Singapore), (arXiv:1604.05955)

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  21. [21]

    K., Acharyya, K., & Molteni, D., 2004, AAP, 421, 1

    Chakrabarti, S. K., Acharyya, K., & Molteni, D., 2004, AAP, 421, 1

    work page 2004

  22. [22]

    K., & Molteni, D.\ 1993, , 417, 671

    Chakrabarti, S. K., & Molteni, D.\ 1993, , 417, 671

    work page 1993

  23. [23]

    K., & Molteni, D., 1995, MNRAS, 272, 80

    Chakrabarti, S. K., & Molteni, D., 1995, MNRAS, 272, 80

    work page 1995

  24. [24]

    K., & Sahu, S

    Chakrabarti, S. K., & Sahu, S. A., 1997, AAP, 323, 382

    work page 1997

  25. [25]

    K.\ 2016, , 462, 3502

    Deb, A., Giri, K., & Chakrabarti, S. K.\ 2016, , 462, 3502

    work page 2016

  26. [26]

    K., & Mondal, S.\ 2014, , 440, L121

    Debnath, D., Chakrabarti, S. K., & Mondal, S.\ 2014, , 440, L121

    work page 2014

  27. [27]

    Ertan, \"U .\ 2018, , 479, L12

    work page 2018

  28. [28]

    K., 2013, MNRAS, 430, 2836

    Giri, K., & Chakrabarti, S. K., 2013, MNRAS, 430, 2836

    work page 2013

  29. [29]

    A.\ 2017, , 471, 4827

    G \"u gercino g lu, E., & Alpar, M. A.\ 2017, , 471, 4827

    work page 2017

  30. [30]

    K., 1998, MNRAS, 299, 799

    Lanzafame, G., Molteni, D., & Chakrabarti, S. K., 1998, MNRAS, 299, 799

    work page 1998

  31. [31]

    W., 2002, ApJ, 580, 423

    Mauche, C. W., 2002, ApJ, 580, 423

    work page 2002

  32. [32]

    M \'e ndez, M., 2006, MNRAS, 371, 1925

    work page 2006

  33. [33]

    K., 2001, ApJl, 563, L57

    Molteni, D., Acharya, K., Kuznetsov, O., Bisikalo, D., & Chakrabarti, S. K., 2001, ApJl, 563, L57

    work page 2001

  34. [34]

    K., 1994, ApJ, 425, 161

    Molteni, D., Lanzafame, G., & Chakrabarti, S. K., 1994, ApJ, 425, 161

    work page 1994

  35. [35]

    K., 1996, ApJ, 470, 460

    Molteni, D., Ryu, D., & Chakrabarti, S. K., 1996, ApJ, 470, 460

    work page 1996

  36. [36]

    K., 1996, ApJ, 457, 805

    Molteni, D., Sponholz, H., & Chakrabarti, S. K., 1996, ApJ, 457, 805

    work page 1996

  37. [37]

    J., 1992, ARAA, 30, 543

    Monaghan, J. J., 1992, ARAA, 30, 543

    work page 1992

  38. [38]

    J., 1980, AAP, 88, 23

    Paczy \'n sky, B., & Wiita, P. J., 1980, AAP, 88, 23

    work page 1980

  39. [39]

    Revnivtsev, M., & Mereghetti, S.\ 2015, , 191, 293

    work page 2015

  40. [40]

    K., & Molteni, D., 1997, ApJ, 474, 378

    Ryu, D., Chakrabarti, S. K., & Molteni, D., 1997, ApJ, 474, 378

    work page 1997

  41. [41]

    Sanna, A., Riggio, A., Burderi, L., et al.\ 2017, , 469, 2

    work page 2017

  42. [42]

    Shang, J.-R., Debnath, D., Chatterjee, D., et al.\ 2019, , 875, 4

    work page 2019

  43. [43]

    I., & Sunyaev, R

    Shakura, N. I., & Sunyaev, R. A., 1973, AAP, 24, 337

    work page 1973

  44. [44]

    @esa ( ) , n @biblabelnum##1 ##1

    \@ifclassloaded aguplus natbib The aguplus class already includes natbib coding, so you should not add it explicitly Type <Return> for now, but then later remove the command natbib from the document \@ifclassloaded nlinproc natbib The nlinproc class already includes natbib coding, so you should not add it explicitly Type <Return> for now, but then later r...

    work page

  45. [45]

    @stdbsttrue NAT@ctr \@lbibitem[ NAT@ctr ] \@lbibitem[#1]#2 \@extra@b@citeb \@ifundefined br@#2\@extra@b@citeb \@namedef br@#2 \@nameuse br@#2\@extra@b@citeb \@ifundefined b@#2\@extra@b@citeb @num @parse #2 [ @natanchorstart #2\@extra@b@citeb \@biblabel @num @natanchorend] @ifcmd#1()()\@nil #2 @lbibitem\@undefined @lbibitem\@lbibitem \@lbibitem[#1]#2 @lbib...

    work page

  46. [46]

    @open @close @open @close and [1] URL: #1 \@ifundefined chapter * \@mkboth \@ifundefined NAT@sectionbib * \@mkboth * \@mkboth\@gobbletwo \@ifclassloaded amsart * \@ifclassloaded amsbook * \@ifundefined bib@heading @heading NAT@ctr thebibliography [1] @ \@biblabel NAT@ctr \@bibsetup #1 NAT@ctr 0 @openbib .11em \@plus.33em \@minus.07em 4000 4000 `\.=1000 \@...

    work page