Issue 
A&A
Volume 626, June 2019



Article Number  A18  
Number of page(s)  14  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201834414  
Published online  05 June 2019 
SuperEddington accretion discs with advection and outflows around magnetized neutron stars
^{1}
Tuorla Observatory, Department of Physics and Astronomy, 20014 University of Turku, Finland
email: anna.chashkina@utu.fi, juri.poutanen@utu.fi
^{2}
Sternberg Astronomical Institute, Lomonosov Moscow State University, Universitetsky pr. 13, 119992 Moscow, Russia
^{3}
Space Research Institute of the Russian Academy of Sciences, Profsoyuznaya str. 84/32, 117997 Moscow, Russia
^{4}
Nordita, KTH Royal Institute of Technology and Stockholm University, Roslagstullsbacken 23, 10691 Stockholm, Sweden
Received:
10
October
2018
Accepted:
24
March
2019
We present a model for a superEddington accretion disc around a magnetized neutron star taking into account advection of heat and the mass loss by the wind. The model is semianalytical and predicts radial profiles of all the basic physical characteristics of the accretion disc. The magnetospheric radius is found as an eigenvalue of the problem. When the inner disc is in radiationpressuredominated regime but does not reach its local Eddington limit, advection is mild, and the radius of the magnetosphere depends weakly on the accretion rate. Once it approaches the local Eddington limit the disc becomes advectiondominated, and the scaling for the magnetospheric radius with the mass accretion rate is similar to the classical Alfvén relation. Allowing for the mass loss in a wind leads to an increase in the magnetospheric radius. Our model can be applied to a wide variety of magnetized neutron stars accreting close to or above their Eddington limits: ultraluminous Xray pulsars, Be/Xray binaries in outbursts, and other systems. In the context of our model we discuss the observational properties of NGC 5907 X1, the brightest ultraluminous pulsar currently known, and NGC 300 ULX1, which is apparently a Be/Xray binary experiencing a very bright superEddington outburst.
Key words: accretion, accretion disks / stars: neutron / stars: magnetic field / Xrays: binaries
© ESO 2019
1. Introduction
Mass transfer rates in binary systems may vary in very broad limits, from very low to amounts vastly exceeding the Eddington limit of the accretor. Highly superEddington accretion rates are not surprising in binary systems containing a neutron star (NS) and a massive star filling its Roche lobe. In this case, the mass transfer runs on a relatively short thermal timescale of a massive star, which is hundreds of thousands of years. The mass transfer does not quench because the high mass ratio makes a binary system more likely to remain in contact or even to tighten, thus increasing the mass transfer rate (see e.g. Pavlovskii & Ivanova 2015; Pavlovskii et al. 2017).
The strong magnetic field makes it possible to transport the accreted matter deep into the gravitational well to the NS surface. Presumably, the matter falls to the NS surface in a thin curtain (see Fig. 1a in Basko & Sunyaev 1976). The high ratio of the radiating surface area to the volume, and the reduced scattering cross section in the strong magnetic field are invoked to explain the observed excess above the Eddington limit (see e.g. Mushtukov et al. 2015; Kawashima et al. 2016).
It has recently been realized that some of the ultraluminous Xray sources (ULX) in nearby galaxies are actually magnetized NSs accreting well above their Eddington limits. Bachetti et al. (2014) using NuSTAR data discovered coherent pulsations in the ULX X2 in the galaxy M82. Later, two similar objects, NGC 7793 P13 (Israel et al. 2017a; Fürst et al. 2016) and NGC 5907 X1 (Israel et al. 2017b), were found. NGC 5907 X1 is exceptional in its luminosity, exceeding 10^{41} erg s^{−1} during some of the observations. These are the prototypical members of the ULX pulsar (ULXP) family. More recently, other ULXPs have been found. Pulsations with a period of 20–30 s were discovered in NGC 300 ULX1, the supernova impostor SN2010da (Villar et al. 2016), with a very strong spinup of ṗ = −1.75 × 10^{−7} s s^{−1} (Bachetti et al. 2018). Also, M51 ULX7 was identified as a ULXP (G. Israel, priv. comm.); and finally, the first ULXP in the Milky Way, Swift J0243.6+6124, was discovered (Kennea et al. 2017; Tsygankov et al. 2018; WilsonHodge et al. 2018). Though only a few persistent superEddington objects have been robustly identified as NSs to date, it is quite natural to expect more supercritically accreting NSs among the more general class of ULXs (see Kaaret et al. 2017).
Studies of superEddington accretion discs started simultaneously with the standard disc theory by Shakura & Sunyaev (1973). They suggested that, for the accretion rate Ṁ_{0} exceeding a certain critical value Ṁ_{cr}, a wind emanates inside the spherization radius R_{sph}, driven by the radiation pressure, and removes just the right amount of matter to keep the disc at the local Eddington limit. The luminosity of the supercritical disc in this case exceeds the Eddington luminosity L_{Edd} by a logarithmic factor 1 + ln(Ṁ_{0}/Ṁ_{cr}), and the accretion rate decreases with radius as Ṁ(R) = Ṁ_{0} R/R_{sph}. We refer to such a scenario as “classical mass loss”. Implications of this scenario to the discs around magnetized NSs were considered by Lipunov (1982), and recently by Grebenev (2017) and Mushtukov et al. (2019). Remarkably, because the properties of the inner disc in this model are independent of the outer boundary conditions, the size of the magnetosphere and the luminosity at high Ṁ_{0} converge to universal values dependent on the NS magnetic moment μ alone.
Later it was realized that it is important to take into account other effects related to deviations from the thin disc approximation. The most significant departure from the standard disc model is the heat advection, which violates the locality assumption in the energy balance of the disc. The role of advection at low accretion rates was considered, for example, by Ichimaru (1977), Narayan & Yi (1994, 1995), and Abramowicz et al. (1995, 1996), and at high accretion rates by Begelman & Meier (1982), Abramowicz et al. (1988), and Beloborodov (1998). Numerical simulations of supercritical accretion onto a black hole were performed by Eggum et al. (1988), Ohsuga et al. (2005), Okuda et al. (2005), and more recently by Ohsuga & Mineshige (2011), McKinney et al. (2014), Sa̧dowski et al. (2014), Sa̧dowski & Narayan (2016), and Ogawa et al. (2017). Semianalytical models for supercritical discs including heat advection and outflows were constructed by Lipunova (1999), Kitabatake et al. (2002), Fukue (2004), and Poutanen et al. (2007).
While most of this theoretical work on supercritical accretion was devoted to discs around black holes, the discovery of ULXPs has drawn attention to magnetized NSs accreting at high rates. Numerical magnetohydrodynamic (MHD) simulations of subcritical disc accretion onto magnetized neutron stars have been carried out in Romanova et al. (2002, 2003), Long et al. (2005), who studied extensively the properties of the magnetospheric MHD flow and the processes at the inner boundary of the disc, and the formation of funnelshaped flows in particular. They consider magnetically driven winds, which are particularly favoured in the propeller regime (see review by Romanova & Owocki 2015). In our work we focus instead on radiatively driven outflows, which are expected to be launched by supercritical discs. It is important to note that even in highly supercritical ULXPs like M82 X2 and NGC 7793 P13, the accretion disc outside the magnetosphere may remain in a subcritical regime. Chashkina et al. (2017; hereafter CAP) considered accretion onto a magnetized NS in a regime that can be characterized as intermediate, with an accretion rate in the range (10^{−9}−10^{−6}) M_{⊙} yr^{−1} for μ ∼ 10^{31}−10^{32} G cm^{3} inferred for ULXPs (see Tsygankov et al. 2016). The disc outside the magnetosphere in this case remains geometrically thin, nearly Keplerian, and does not lose considerable amounts of matter. This allowed us to use certain results of the standard disc theory such as local radiation energy balance and simple conservation laws for mass and angular momentum. This approach works well if the formally calculated spherization radius is smaller than the Alfvén radius. However, superEddington NSs with lower magnetic moments or those accreting at higher rates should possess a supercritical advective disc, whose radiation pressure becomes sufficient to unbind part of the accreting material. Recently, simulations of supercritical accretion onto a nonmagnetized NS have been performed by Abarca et al. (2018) in an attempt to explain nonpulsating ULXs.
In the present work, we develop a model of a supercritical accretion disc around a magnetized NS, applicable to a broader class of objects, including NSs with pulsarscale magnetic fields, luminosities hundreds of times higher than the Eddington limit, and in particular extremely bright sources like NGC 5907 X1. The structure of the paper is as follows. In Sect. 2 we describe our model and present the main equations describing disc accretion accounting for the effects of advection and the mass and momentum loss in a wind. Section 3 is devoted to the results; we present the dependences of the magnetospheric radius on the mass accretion rate and NS magnetic moment. We also consider the effects of the pulsar spin and the irradiation of the disc by the central source on the disc structure. We discuss our results and apply them to particular ULXPs in Sect. 4. We conclude in Sect. 5.
2. Model
2.1. Basic equations
We consider a NS of mass M, radius R_{NS}, and magnetic moment μ. We assume it is an aligned rotator with an unperturbed dipolar magnetic field inside the magnetospheric radius, which is equal to the inner radius of the disc R_{in}. Outside this radius, the magnetic field lines are opened by the ideally conducting accretion disc, as proposed by Parfrey et al. (2016). A simple sketch of the adopted structure of the disc is shown in Fig. 1. All the interactions between the disc and the magnetosphere are assumed to occur in a narrow strip at the edge of the disc and are reduced to two boundary conditions introduced in CAP and later in this section. One of the conditions, the pressure balance, as shown in 3D MHD simulations by Kulkarni & Romanova (2013) for lower accretion rates, corresponds to the radius from which the matter starts to flow along a funnel stream.
Fig. 1. Structure of an accretion disc around a ULXP. For very high mass accretion rates, the inner parts of the disc (inside the spherization radius R_{sph} shown here) enter the superEddington accretion regime. Inside R_{sph} the thin disc model is not applicable. Mass loss in a wind is shown by blue arrows. The accretion column, where most of the energy is released, is shown by yellow cones, and the red wavy lines refer to the radiation of the column that may affect the inner disc pressure balance. The red vertical line marks the effective boundary between the disc and the magnetosphere at R_{in}. 
In our model, we include three effects important for an accretion disc at a near or supercritical accretion rate: (i) advection, (ii) mass loss in the wind, and (iii) angular momentum loss in the wind. Advection is related to the increasing photon diffusion timescales at high mass accretion rates. It alters the energy equation. The matter loss in the wind makes the disc accretion rate dependent on radius, and also affects the angular momentum and the energy conservation equations because the wind carries some angular momentum and energy. Poutanen et al. (2007) assumed that the specific angular momentum of the wind is equal to that of the matter in the disc. However, if the wind is magnetized, it can efficiently remove the angular momentum from the disc, as in the centrifugally driven wind model of Blandford & Payne (1982).
We follow the basic framework of CAP and retain some of the important assumptions of the model. First of all, the boundary conditions at the inner disc edge remain the same. We take the torque balance in the form
and assume pressure balance, which yields for the viscous stress tensor the following relation:
Here, Ṁ_{in} is the mass accretion rate at the inner edge of the disc; Ω_{NS} is the angular velocity of the NS; Ω_{in}, H_{in}, and are the angular velocity, the halfthickness, and the vertically integrated rϕcomponent of the viscous stress tensor at the inner boundary of the disc R_{in}; the dimensionless constant k_{t} parametrizes the efficiency of the angular momentum removal by the magnetic and viscous torques; and α is dimensionless viscosity parameter as defined by Shakura & Sunyaev (1973). The luminosity L = ηṀ_{in} c^{2} is released close to the NS in the accretion column with efficiency η as measured by an imaginary observer in the equatorial plane at the inner edge of the disc. Because the accretion column emits anisotropically, the efficiency η may differ from the angleaveraged efficiency of an isotropic source and from the efficiency measured by an observer at infinity. The radial structure of the disc is described by the angular frequency Ω, thickness H, and the viscous stress W_{rϕ}. The angular frequency and the viscous stress should conform to the boundary conditions (1) and (2). Equation (1) is similar to the inner boundary condition in Spruit & Taam (1993) and Rappaport et al. (2004) in the case of a slowly rotating NS. The additional term in the RHS is radiation drag that can be important at high luminosities L > L_{Edd}. We note that Eq. (1) is valid only for the case of the disc rotating faster than the magnetosphere. For the propeller case (Ω_{in} < Ω_{NS}), the boundary condition would be different. Moreover, one could consider a radial distribution of the magnetic torque applied to the disc, as was done, for example, by Kluzniak & Rappaport (2007). On the other hand, Matt & Pudritz (2005b) argue that if Ω_{in} > Ω_{NS}, the magnetic interaction can be limited to a narrow ring.
We take the radial component of the momentum equation in the form
where is the surface density and is the vertically integrated pressure. As in CAP, we ignore the dynamic term v_{R} dv_{R}/dR, which is suppressed by a factor of α^{2}(H/R)^{2} ≪ 1 with respect to the pressure gradient term.
In all the equations for the radial structure we use vertically integrated quantities. In the equations for the vertical structure we assume a fixed vertical effective polytropic index n, defined through the equation for the density ρ = ρ_{c}[1 − (z/H)^{2}]^{n} and for the pressure P = P_{c}[1 − (z/H)^{2}]^{n + 1} (Paczynski & Jaroszynski 1978). This allows us to compute analytically the connection factors between the midplane quantities such as ρ_{c} and P_{c}, and the corresponding vertically integrated values Σ and Π, as described in Appendix A. For the vertical temperature profile used in the advection term calculations we use the expression T = T_{c}[1 − (z/H)^{2}]^{(n + 1)/4}, which is valid for the radiationpressuredominated case as T ∝ p^{1/4}; here T_{c} is the central temperature on the disc.
The angular momentum conservation equation is modified by an additional term corresponding to the angular momentum outflow in the wind,
where ψ ≥ 1 allows us to scale up the net angular momentum lost in the wind. If ψ = 1, the net angular momentum in the wind is equal to that in the disc (Lipunova 1999; Poutanen et al. 2007). Larger ψ can appear in magnetocentrifugal winds, where the subAlfvénic part of the transonic flow rotates approximately rigidly, thus increasing the angular momentum loss by a factor of ∼(R_{wA}/R)^{2}, where R_{wA} is the cylindrical radius of the Alfvén surface, and R is the launch point of the wind streamline (see also Matt & Pudritz 2005a who used this approach for magnetospheric outflows). In the centrifugally driven wind model of Blandford & Payne (1982), for instance, R_{wA}/R may be as high as 5. Simulations of the centrifugally driven winds (e.g. Ustyugova et al. 1999) suggest even greater lever lengths at low accretion rates, while at high accretion rates when the wind is radiatively driven the effect of magnetic stresses is smaller (see e.g. Proga 2000).
The mass accretion rate derivative dṀ/dR needs to be calculated using some additional equations describing the physics of the wind launching. We assume that some fraction ϵ_{w} ≤ 1 of the energy leaving the disc with radiation is spent to accelerate the outflow (Lipunova 1999; Poutanen et al. 2007),
where Ω_{K} is Keplerian angular velocity and T_{eff} is effective temperature. The model of optically thick, energydriven wind developed by Poutanen et al. (2007) gives similar results. Physically, there could be wind losses everywhere in the disc, but we take into account only the continual radiationdriven wind which would operate within the spherization radius R_{sph} defined by the condition on the minimal relative thickness H/R > (H/R)_{cr}. Because the details of wind launching are uncertain, we assume (H/R)_{cr} = 1. We note that for a disc around a black hole, the wind mass loss is determined only by the mass of the black hole and the outer accretion rate. For a disc around a magnetized neutron star, the mass lost in the wind also depends on the magnetic moment of the star and its spin period.
We use the αviscosity prescription for the vertically integrated viscous stress:
The midplane pressure can be expressed as a sum of the radiation and gas pressures,
where m̃ is the mean particle mass and (Eq. (A.5)). Here we have taken into account the αviscosity prescription and the relations between the central and vertically integrated quantities implied by the adopted vertical structure (Eqs. (A.3) and (A.4)). We do not use Eq. (7) directly, but as a supplement to the energy equation (see Sects. 2.2 and 2.3). Treatment of advection is a new part of the model that requires a separate consideration.
2.2. Advection
In a general case with advection, some fraction of energy is radially transported. The energy flux carried by radiation diffusion in the vertical direction is no longer equal to the local energy release,
where the total energy released in the disc at radius R per unit area is
and the radiation flux from both sides of the disc is (see Eq. (A.10))
Here, κ ≃ 0.34 cm^{2} g^{−1} is the Thomson scattering opacity. The advected flux Q_{adv} can be viewed (as done in e.g. Lipunova 1999) as the flux of heat carried with the flow, and thus can be expressed through the specific entropy per particle s
where the radial velocity (taken by its absolute value) is v_{R} = Ṁ/(2πRΣ). Taking into account Eq. (7) we can expand the expression above as (see Appendix B)
where the coefficients 𝒮, 𝒫, 𝒬, and ℛ are given by Eqs. (B.8)–(B.11). The energy balance can be written using Eqs. (8), (9), and (12) as
This relation can be rewritten for the derivative of the angular velocity as
where the coefficients C_{Ω, Σ, w, T, free} are given by expressions (B.12)–(B.16).
2.3. Solving the disc equations
Unlike in our previous work (CAP), here there are five variables for which we need to solve differential equations: the angular frequency Ω, the tangential stress W_{rϕ}, the surface density Σ, the mass accretion rate Ṁ, and the midplane temperature T_{c}. There are also seven free global parameters: the viscosity parameter α, the wind efficiency parameter ϵ_{w}, the wind magnetization parameter ψ, the polytropic index n, the spin period p and magnetic moment of a NS μ, and the accretion efficiency η. We calculate the disc structure from its outer edge to the inner boundary.
We solve the five equations listed below:

The derivative of vertically integrated viscous stress that can be obtained from the radial Euler Eq. (3) and αviscosity prescription (6):

Whenever the outflow condition (H/R > (H/R)_{cr}) is satisfied, the mass accretion rate changes with radius according to Eq. (5):

The differential equation for the central temperature can be obtained by taking derivative of the pressure Eq. (7):
Here, we used the thickness of the disc H following from the hydrostatic equilibrium (see Eq. (A.2)) together with the gradient of W_{rϕ} substituted from Eq. (14). We also use the ratio of the gas pressure to the total pressure in the equatorial plane expressed using the adopted vertical structure (see Eq. (A.11)):

To determine the radial dependence of the angular velocity, we substitute Eqs. (15) and (16) into angular momentum conservation Eq. (4):
The second term on the righthand side is switched off for a thin disc (H/R < (H/R)_{cr}) when there are no outflows.

The differential equation for the surface density was obtained from the advection Eq. (14) substituting all other derivatives from Eqs. (15)–(19):
The disc structure equations are further converted to a more compact dimensionless form in Appendix C.
Our solution should satisfy a number of boundary conditions. The accretion rate at the outer boundary of the integration region should be equal to Ṁ_{0}, which is one of the global parameters of the simulation. Also two variables, Ω and W_{rϕ}, should satisfy boundary conditions (1) and (2). In order to achieve that we use two adjustable parameters. The first one is the tangential stress at the outer boundary . The second one is the magnetospheric radius R_{in}, which we parametrize through the ratio
where
is the Alfvén radius (see e.g. Elsner & Lamb 1977 and Sect. 6.3 in Frank et al. 2002). As a result of the iteration procedure we get the five sought for functions and the value for ξ. Because the accretion rate at the inner boundary of the disc is physically more relevant, we also use a differently normalized version of the relative magnetospheric radius:
3. Results
3.1. Global parameters
The relative magnetospheric size ξ is found as one of the two eigenvalues of the problem, the other being . Other global parameters obtained in the model include several quantities at the magnetospheric boundary: the relative thickness of the disc (H/R)_{in}, the fraction of mass reaching the magnetosphere Ṁ_{in}/Ṁ_{0}, and the advection fraction (Q_{adv}/Q^{+})_{in}. We also track the maximum thickness of the disc (H/R)_{max}.
All the simulations were made for the following values of the parameters: α = 0.1, k_{t} = 0.5, ϵ_{w} = 0.5, ψ = 1, n = 1, and η = 0 if not stated otherwise. We first consider a slowly rotating accretor with the spin period equal to ten equilibrium periods p = 10p_{eq},^{1} where the equilibrium period (Lipunov 1992; Illarionov & Sunyaev 1975) is defined as the spin period for which the Alfvén radius R_{A} equals the corotation radius ,
where μ_{30} = μ/10^{30} G cm^{3} is the dimensionless NS magnetic moment and m = M/1.4M_{⊙} is the normalized NS mass. We also use the dimensionless mass accretion rate ṁ_{0} = Ṁ_{0}/Ṁ_{Edd} normalized by the Eddington value Ṁ_{Edd} = 4πGM/cκ. In Fig. 2 we show the contours of ξ_{eff} in the ṁ_{0}−μ plane, covering two orders of magnitude in magnetic moment and five orders of magnitude in the mass accretion rate.
Fig. 2. Relative magnetospheric radius ξ_{eff} given by Eq. (23) is colourcoded and shown by contours on the ṁ_{0}− μ plane. The spin period is p = 10 p_{eq}. The luminosity scale is calculated assuming efficiency η = 0.1, but irradiation effects on the disc structure were ignored. 
The magnetospheric radius and the thickness of the disc are tightly related. There is good agreement with the results of CAP, as we can see from Fig. 3. The relative magnetospheric size ξ_{eff} remains a monotonic function of (H/R)_{in} and behaves in approximate accordance with Eq. (57) in CAP (within 5%). This longperiod approximation works fine far from the equilibrium period, for instance, in outbursts.
Fig. 3. Relative magnetospheric radii ξ (dashed blue curve) and ξ_{eff} (green solid) as functions of the disc thickness at the magnetospheric boundary. The dotted red line corresponds to the longperiod asymptotic given by Eq. (57) in CAP. Magnetic moment of the NS was set to μ_{30} = 100. 
Not only mass but also angular momentum is lost in the wind. If there is no angular momentum flow due to stresses in the wind (ψ = 1), the fraction of angular momentum expelled from the disc depends only on the distribution of the mass loss over the radial coordinate. The total angular momentum flux through an annulus in the disc is composed of the angular momentum carried by the matter in the disc and of the viscous torque acting on the annulus,
where ω = Ω/Ω_{K}. Unless some of this angular momentum is removed in the wind, l_{in} = l(R_{in}) should be equal to l_{0} = l(R_{out}). In Fig. 4 we show the ratio of the angular momentum fluxes l_{in}/l_{0} as a function of the ratio of mass accretion rates Ṁ_{in}/Ṁ_{0}. Both ratios start at unity for conservative thin disc accretion and then decrease as Ṁ_{0} increases. The slope of the curve in Fig. 4 shows the evolution of the mean net angular momentum in the wind. With increasing mass accretion rate, the outflow involves larger radii. Characteristic radii losing most of the angular momentum is approximately equal to the spherization radius that leads to the scaling (see Fig. 4):
Fig. 4. Fraction of the angular momentum flux retained in the superEddington disc l_{in}/l_{0} as a function of the fraction of mass reaching the magnetosphere. The dotted black line corresponds to equal ratios (l_{in}/l_{0} = Ṁ_{in}/Ṁ_{0}, as expected if the net angular momentum is constant); the solid green line corresponds to the expected scaling of a growing spherization radius (Eq. (26)). 
The amount of angular momentum lost in a centrifugal wind is enhanced approximately proportionally to ψ.
3.2. Dependence on the spin period
The transition to the propeller regime, when the accretion flow can no longer spin up the NS, can be traced using the fastness parameter ω_{s} defined as
When ω_{s} = 1, the inner rim of the disc rotates exactly with the same frequency as the magnetosphere, making our first boundary condition marginally satisfiable.
The main effect of increasing the fastness parameter (or decreasing the spin period of the accretor) on the properties of our solution is the increasing ratio ξ of the size of the magnetosphere to the Alfvén radius. The factor ξ_{eff} increases by about 40% between the slowly rotating NS case and the propeller limit (see Fig. 5).
Fig. 5. Parameter ξ_{eff} as a function of fastness parameter for magnetic dipole moments μ_{30} = 1 and μ_{30} = 10 and mass accretion rates ṁ_{0} = 10^{4} and ṁ_{0} = 5 × 10^{3}. 
Near the corotation, when R_{co} = R_{in}, the disc thickness approaches zero and the boundary condition for the viscous stresses (2) reduces to the zerotorque condition used in standard disc theory. This allows us to compare our results directly to some of the results obtained using the codes designed for black hole accretion by Poutanen et al. (2007), among others. Their model takes into account the advection effects and the outflows from the disc. We compare the spherization radius defined as the maximum distance from the NS where the condition for outflows H > R is fulfilled.
In Fig. 6, we show the spherization radius normalized by the dimensionless mass accretion rate (green dots) and compare our results to the zerotorque case. For ϵ_{w} = 0.5 and ṁ_{0} = 10^{4}, Eq. (21) in Poutanen et al. (2007) predicts R_{sph} ≃ 0.575 ṁ_{0}(GM/c^{2}), in reasonable agreement with our results close to equilibrium.
Fig. 6. Ratio of the spherization radius r_{sph} (in units of gravitational radius R_{g}) to the dimensionless mass accretion rate ṁ_{0} = 10^{4} as a function of the fastness parameter (green dots). The red horizontal line gives the result from Eq. (21) in Poutanen et al. (2007). 
3.3. Effects of irradiation
All the previous results were calculated without irradiation effects by setting η = 0 (or L = 0 in Eqs. (1) and (2)). The real efficiency affecting the pressure balance is probably of the same order as the integrated accretion efficiency, though strong anisotropy of the radiation from the column is not excluded. Confirming the result of CAP, we find that radiation from the accretion column is an important factor affecting the structure of the disc, in particular, the radius of the magnetosphere. The magnetospheric radius increases by up to 30% for ṁ_{0} ∼ 10^{4}, as is shown in Fig. 7. The effect grows rapidly with the mass accretion rate and with the magnetic moment. The latter is a consequence of a rapidly growing ratio P_{rad}/P_{mag} with the disc inner radius in the dipole approximation (see e.g. CAP, Eq. (63)).
Fig. 7. Colourcoded contours of relative correction to the magnetospheric radius caused by irradiation ξ(η = 0.1)/ξ(η = 0)−1 shown in the ṁ_{0} − μ plane. 
3.4. Dependence on other parameters
There are adjustable parameters that influence the structure of the disc. In Table 1 we show the global properties of the discs for different parameters. The fiducial model here is the model with α = 0.1, ϵ_{w} = 0.5, n = 1, η = 0, and ψ = 1. All the models are calculated for a NS with magnetic field μ = 10^{30} G cm^{3} and accretion rate ṁ_{0} = 3000. Other models differ from the fiducial model by a single parameter.
Disc properties.
As in CAP, we find the viscosity parameter to be an important factor that alters the structure of the disc. For the parameters under consideration, changing α from 0.1 to 0.5 results in a twofold decrease in disc thickness and ξ, and quenches wind formation. Changing the vertical disc structure by increasing the effective polytropic index n also makes the disc thinner and less likely to form outflows.
3.5. Effects of advection and wind
Under the assumptions we use, including the adopted vertical structure and the minimum disc thickness for wind launching, advection starts to play a role rather early, when the entire disc is still subcritical. As a consequence, the disc thickness stabilizes at H ≃ R (see Fig. 8). A slimming effect of advection was noted earlier by Abramowicz et al. (1988), Beloborodov (1998), Lipunova (1999), and Lasota et al. (2016), among others. The local Eddington limit (H = R) is reached at the critical mass accretion rate of
i.e. 5–6 times higher than in CAP (see their Eq. (66)).
Advection starts to dominate in the energy balance already below this limit, at ṁ_{0} ∼ 10^{3} (see Fig. 9), making the inner disc a huge reservoir of heat. For ṁ_{0} above the limit given by Eq. (28), most of the gravitational energy released in the disc is stored as heat. To illustrate this, we calculated the cumulative luminosity of the disc integrated from some radius R to the outer radius . Similarly, by integrating Q_{rad} and Q_{adv}, we can define the cumulative radiative and advection powers. The total cumulative power is shown by a blue dashed line in Fig. 10; at the inner radius R_{in} it is in full agreement with theoretical prediction, L_{theor} = GMṀ/2R_{in}. We note that the ratio of the disc luminosity to the luminosity of the accretion column can be as low as ∼R_{NS}/R_{in}.
Fig. 8. Relative disc thickness as a function of the radius for a model with μ = 10^{30} G cm^{3}, ṁ_{0} = 3 × 10^{3}, and p = 0.67 s. Our results are shown by the red dotted curve, whereas the solid green and blue curves correspond to the asymptotics for the zones B and A of the standard disc, respectively. 
Fig. 9. Advection effects for the model shown in Fig. 8. The advection flux fraction Q_{adv}/Q^{+} and the radiated fraction Q_{rad}/Q^{+} are shown by the red solid and green dashed curves, respectively. 
Fig. 10. Cumulative luminosities as functions of radius. The total power generated in the disc is shown by the dashed blue curve. The solid red and dotted black curves show the radiative and advected luminosities, respectively. The parameters of the model are μ_{30} = 1 and ṁ_{0} = 3000. 
In Fig. 11, we plot the fraction of the initial mass accretion rate remaining in the disc for different ṁ_{0} as a function of radius. At very high accretion rates, the supercritical wind blows away a considerable part of the accretion material and operates at all radii within R_{sph}. At some intermediate rates, 2000 ≲ ṁ_{0} ≲ 3000, there is a prominent subcritical region near R_{in} where there is no wind. This is caused by the nonmonotonic dependence of the disc height on radius. We note that the amount of the blownaway material depends on the condition used to switch on the wind. The wind launching relies on complex physical processes and can be described only very approximately in 1D models.
Fig. 11. Fraction of the mass accretion rate reaching radius R for a NS with magnetic moment μ_{30} = 1. The lines from top to bottom correspond to different ṁ_{0} in the interval 2000–3500 with steps of 500. 
3.6. Magnetospheric radius for different accretion regimes
One of the most important outputs of our model is the radius of the magnetosphere. In many standard models (Ghosh et al. 1977; Koenigl 1991; Wang 1996; Kluzniak & Rappaport 2007), the magnetospheric radius is supposed to scale with the Alfvén radius, corresponding to the constant ξ in our notation. It is worth noticing that MHD calculations of magnetospheric flows with relatively weak accretion rates (Bessolaz et al. 2008; Kulkarni & Romanova 2013) demonstrate similar values of ξ ∼ 0.4 to those in Fig. 2. Chashkina et al. (2017) have demonstrated that the radius of the magnetosphere interacting with a thin radiationpressuredominated disc depends weakly on the mass accretion rate. Accounting for the effects of advection and wind losses makes the picture more complicated.
The dependence of the magnetospheric radius (in units of gravitational radius R_{g} = GM/c^{2}) on the mass accretion rate is shown in Fig. 12 for a wide range of accretion rates. The inner disc regions of most Xray pulsars are in the gaspressuredominated regime. As the accretion rate increases, radiation pressure becomes important. For pulsarscale magnetic fields, μ ∼ 10^{30} G cm^{3}, this happens at luminosities of a few times L_{Edd}, which are quite reachable, for instance, in Be/Xray binaries during strong outbursts like the recent superEddington outburst of SMC X3 (Townsend et al. 2017; Tsygankov et al. 2017; Weng et al. 2017). As was shown in CAP (see Eq. (61)), the magnetospheric radius becomes almost independent of the accretion rate, if the radiation pressure dominates at the inner edge of a subcritical disc:
Fig. 12. Magnetospheric radius in units of R_{g} as a function of the accretion rate for a NS with magnetic moment μ_{30} = 1. Parts of the black solid curve with different slopes correspond to the different regimes of accretion near the magnetospheric boundary. Two standard solutions () are plotted with the grey dotted lines: ξ = 0.5 and ξ = 1 (spherically symmetric case). 
Thus, provided there is a direct measurement of the magnetospheric radius, for example, from quasiperiodic oscillations, we can directly estimate the magnetic field of a NS, with a weak dependence on the viscosity parameter α.
The inner disc radius is defined mainly by the balance of pressures. The pressure inside the disc is related to its thickness. Thus, the dependence of H/R on ṁ_{0} is crucial for the behaviour of r_{in}(ṁ_{0}). As the accretion rate increases, advection starts to play an important role. The relative thickness of the disc is no longer proportional to ṁ_{0}, and the magnetospheric radius again depends on ṁ_{0}. The interplay between wind losses and advection makes the radius dependence on mass accretion rate shallower than the ξ = const approximation historically proposed for spherical accretion, but stronger than r_{in} = const.
With the increasing magnetic field, the magnetospheric radius increases, and all the boundaries between different regimes shift to higher accretion rates, as shown in Fig. 13. In addition, the length of the plateau corresponding to the thin radiationpressuredominated inner disc gradually decreases with the magnetic field, becoming effectively zero at μ ∼ 10^{32} G cm^{3}. Weakly magnetized objects, on the other hand, should have a prominent region of constant magnetosphere size. The plateau starts when radiation pressure begins to dominate over gas pressure at the radius of the magnetosphere. Position of the boundary between gas and radiationpressuredominated regions of the standard disc scales with mass accretion rate as R_{ab} ∝ ṁ^{16/21} (Shakura & Sunyaev 1973). Because R_{in} ∝ μ^{4/9} in the radiationpressuredominated regime (see Eq. (29)), the left boundary of the plateau depends on the magnetic moment as ṁ_{in,left} ∝ μ^{7/12}. The right boundary of the plateau is determined by advection effects. Advection becomes important when Q_{adv} ∼ Q_{rad}, which implies H_{in} ∼ R_{in}. Hence the radius at which the inner disc becomes advective scales linearly with the mass accretion rate, and the accretion rate at the right end of the plateau is ṁ_{in,right} ∝ μ^{4/9}. Thus, the length of the plateau slowly decreases with magnetic moment as ṁ_{in,right}/ṁ_{in,left} ∝ μ^{−5/36}.
It is interesting to compare our results with the classical prescriptions (Ghosh et al. 1977; Wang 1996). Figure 12 demonstrates that the classical dependences are much steeper, having the slope of δ ≡ d log r_{in}/d log ṁ = −2/7. In the case of μ = 10^{30} G cm^{3}, the slope is δ ≈ −0.21 in the gaspressure dominated case, δ ≈ −0.07 in the radiationpressuredominated regime, and δ ≈ −0.16 when advection dominates. The evolution of the local slope δ is plotted in Fig. 14. At low accretion rates, corresponding to the gaspressuredominated thin disc, we get δ ∼ −0.23 in accordance with the results of CAP (see their Sect. 5.1). The maximum value of δ depends on how prominent the thin radiationpressuredominated part of the disc is for a given magnetic field, changing from nearly zero for small μ to about −0.1 for magnetarscale fields. The highest mass accretion rates tend to reproduce much steeper dependences, approaching δ ≃ −2/7 ≃ −0.29. We note that, for a comparatively narrow range of parameters (relevant for T Tauri star magnetospheres), Kulkarni & Romanova (2013) have obtained δ ≃ −0.2 from their 3D MHD simulations.
Fig. 14. Slope of the dependence of the magnetospheric radius on the accretion rate for the models shown in Fig. 13. 
4. Application to ULXPs
4.1. NGC 5907 X1 as a ULXP with a supercritical accretion disc
The ULXP NGC 5907 X1 has a huge period derivative, even after averaging in time: its period has changed from 1.43 to 1.13 s over the 10 years of observations (Israel et al. 2017b). It is also remarkable that its luminosity exceeds 10^{41} erg s^{−1} during some of the observations. The maximum detected period derivative (by absolute value) was ṗ = −5 × 10^{−9} s s^{−1}, about an order of magnitude higher than the average value for this object. The high value of p̄ was put forward as an argument for this object being in a pure spinup state, with the unconstrained braking term in the angular momentum equation being negligible.
The rate of change of the total angular momentum of a NS can be written as
where I ≃ (1−2) × 10^{45} g cm^{2} is the moment of inertia of the NS, and and K_{sd} are spinup and spindown torques. Ignoring the unknown spindown contribution provides the constraint
resulting in an inequality for accretion rate
where I_{45} = I/10^{45} g cm^{2}, and ṗ_{−12} = ṗ/(10^{−12} s s^{−1}). Alternatively, having an independent estimate for the accretion rate, we then can set a lower limit for the magnetic field of this object. For a maximum accretion rate of ∼6000, corresponding to isotropic Xray emission, the normalized magnetic moment μ_{30} cannot be less than ∼0.06.
In the opposite case, when a NS is close to equilibrium, its magnetospheric radius is about the size of the corotation radius. Thus, we can put an upper limit on the magnetic field, which suggests that the NS is still in the accretion regime. In terms of mass accretion rate, this condition may be written as
Figure 15 shows these two limits as applied to NGC 5907 X1. The spinup line shows the lower limit for Ṁ_{0} set by inequality (32) using period derivative ṗ = −5 × 10^{−9} s s^{−1}. The area below this line is forbidden unless some additional spinup process is present. The propeller line shows the boundary of the region where the magnetospheric radius is equal to the corotation radius. The region below this line is prohibited because the accretion disc rotates more slowly than the magnetosphere of the NS, and no stable accretion is possible. The maximum bolometric luminosity of this object is about 2 × 10^{41} erg s^{−1} (Israel et al. 2017b), which leads to a lower limit on the mass accretion rate ṁ_{0} > ṁ_{in} ∼ 6500 assuming efficiency η ∼ 0.15. This value of efficiency does not take into account any beaming effects that in principle can alter the observed efficiency value. This figure provides evidence for a truly high mass accretion rate in NGC 5907 X1, as significantly low ṁ_{0} ≲ 10^{3} are forbidden for any magnetic moments. Therefore, any beaming exceeding a factor of 5–6 is unlikely because with increasing beaming, the red line in Fig. 15 will move down and there will be no allowed region left. Another argument against strong beaming is the observed high pulse fractions and nearly sinusoidal profiles, inconsistent with collimation by a wind, and indicating that the emission site is seen during a large fraction of the spin period.
Fig. 15. Restrictions for the position of NGC 5907 X1 in the ṁ_{0} − μ plane. The red horizontal line corresponds to the mass accretion rate of ṁ_{0} = 6500 (accretion efficiency η ≃ 0.15). The black solid line corresponds to p_{eq} = p = 1.137 s. The dotted spinup line shows the lower limit for ṁ_{in} set by inequality (32) using period derivation ṗ = 5 × 10^{−9} s s^{−1}. We use here M = 1.4 M_{⊙}, R = 10 km, I = 1.5 × 10^{45} g cm^{2}, and p = 1.137 s. 
We can set an upper limit for the magnetic field of NGC 5907 X1 as μ ≤ 7.5 × 10^{31} G cm^{3} and μ ≤ 5.45 × 10^{31} G cm^{3} if we take into account irradiation from the column (see Sect. 3.3). There is evidence for bimodal distribution in luminosities of NGC 5907 X1 (see Fig. S2 in Israel et al. 2017b) that can be interpreted as a manifestation of the propeller effect (similarly to M82 X2 in Tsygankov et al. 2016). This would mean that the source is close to the propeller line. Then instead of an upper limit on the magnetic moment, we get its accurate estimate. Beaming does not play a major role as magnetic field weakly depends on ṁ_{0} (as , see Eq. (33)). The disc is expected to be supercritical (i.e. having winds) in its inner parts if μ_{30} ≲ 14 (see Eq. (28)). At a pulsar magnetic field μ_{30} ∼ 1, the expected outflow rate from the disc is about 20% of the inflow rate, or 5 × 10^{−6} M_{⊙} yr^{−1}.
4.2. NGC 300 ULX1: a ULXP with a strong spinup
The source NGC 300 ULX1 was identified as a ULXP only recently by Carpano et al. (2018). It is characterized by a moderate peak luminosity of L ≃ 4.7 × 10^{39} erg s^{−1}, but the observed spinup rate of this source is exceptionally high: its spin period has changed from 126 to 19 s in about 4 years and is still decreasing (see Vasilopoulos et al. 2018). Its highest spinup rate is estimated as ṗ ≃ −5.5 × 10^{−7} s s^{−1}, which is the highest ṗ ever observed from an accreting NS. The long period makes NGC 300 ULX1 a promising candidate for a NS in a pure spinup state. This can be seen in Fig. 16, where the observed properties of the object are in good agreement with a pulsarscale μ ∼ 10^{30} G cm^{3} of a NS that is currently being rapidly spun up. As in the case of NGC 5907 X1 considered above, pure spinup gives us a lower limit on the magnetic field while the propeller limit gives us an upper limit. Thus, the magnetic field of NGC 300 ULX1 is in the range μ = (1.5−85) × 10^{30} G cm^{3}, which is not affected much by the effects of irradiation. The range is very wide, but sets an upper limit for the beaming factor of 2.5 only, consistent with the constraints set by Binder et al. (2018) based on the number of He IIionizing quanta. The lower boundary here corresponds to pure spinup and implies an equilibrium period of about p_{eq} ∼ 0.2 s. To reach such a period, the object needs to accrete at the observed rate for about t ∼ p_{s}/ṗ ∼ 1 yr.
Fig. 16. Restrictions for the position of NGC 300 ULX1 in the ṁ_{0} − μ plane. The red horizontal line corresponds to the luminosity of 4.7 × 10^{39} erg s^{−1} assuming accretion efficiency η ≃ 0.15. The black solid line corresponds to p_{eq} = p = 20 s. The dotted spinup line shows the lower limit for ṁ_{0} set by inequality (32) using period derivation ṗ = 5.5 × 10^{−7} s s^{−1} (Carpano et al. 2018). We use here M = 1.4 M_{⊙}, R = 10 km, I = 1.5 × 10^{45} g cm^{2}, and p = 20 s. 
At the moment of the first observations of NGC 300 ULX1 in 2010, the NS should have had a period even longer than the 45 s measured in 2016. The question then arises why the NS had such a long period to start with. This period may have originated from a long episode of a very low rate of wind accretion if we assume that this NS was born rotating much faster. To estimate the minimum timescale of the spindown to such a long period, we use Eq. (30), neglect any spinup torque acting on the NS, and parametrize the spindown torque as (Lipunov 1992). The time required to spin down to the observed period is about years. This means that, in order to explain a recently observed long spin period, the system with a normal pulsarscale magnetic field should have spent tens of thousands of years accreting at a very low rate, constantly in propeller regime. The very high luminosity that we observe now is an exceptionally rare event in this system. This is consistent with the extraordinary properties of the optical B[e] transient SN2010da; this object is associated with the disappearance of a huge amount of hot dust observed before the outburst, very high peak optical luminosity, and a bright B[e] supergiant observed after the event (Villar et al. 2016). All this fits well with a catastrophic event in a previously faint Be/Xray binary.
5. Conclusions
Our model provides a simple, physically motivated description of accretion onto a NS with a magnetosphere, where the interaction between the accretion disc and the magnetosphere of the NS is reduced to two boundary conditions. This allows us to reconstruct the structure of the disc and obtain a relative magnetospheric radius ξ, which is important for describing the spin evolution of magnetized NSs and for interpreting observational data on Xray pulsars.
Depending on the dipole magnetic moment of a NS and the mass accretion rate, the inner parts of the accretion disc may appear in different regimes. For classical Xray pulsars (μ ∼ 10^{30} G cm^{3}, ṁ ≲ 10), the accretion disc remains thin and gaspressuredominated, which implies a nearly classical scaling R_{in} ∝ ṁ − 0.22. As the mass accretion rate increases, a large portion of the disc can exist in a radiationpressuredominated regime. Unlike in the radiationpressuredominated disc without advection where the inner radius is independent of accretion rate (CAP), the present model with advection yields R_{in} ∝ ṁ^{−(0.05−0.1)}. For ULXPs, where the mass accretion rates reach ṁ_{0}ṁ_{0}∼ṁ_{0}10^{3}−10^{4}, the effects of advection and mass loss by the wind make the scalings similar to a spherical accretion case, resulting in a trend of about R_{in} ∝ ṁ^{−(0.2−0.3)}. However, if the strength of the magnetic field is one to two orders of magnitude higher than the usual pulsar values, the magnetospheric radius is larger, making the inner disc again geometrically thin and leading to a nearly flat dependence of the magnetospheric radius on ṁ.
Predictions for the magnetospheric radius can be tested through timing analysis of the stochastic component of the variability of Xray sources, where quasiperiodic oscillations and breaks in the powerdensity spectra likely trace the dynamical timescales at the inner rim of the disc. Another test is spin period dynamics. We can compare the observed ṗ with theoretical predictions, as was done here for NGC 5907 X1 and NGC 300 ULX1. The constraints we obtain from the observations of these two objects allow for a rather broad range of magnetic moments. However, the observational data set upper limits for beaming (not more than a factor of 5–6), confirming that ULXPs are intrinsically very luminous objects rather than Xray sources whose luminosity is amplified by an order of magnitude or more as a result of anisotropy of their emission.
When p ≫ p_{eq}, the exact value of the spin period barely affects ξ_{eff}; see Sect. 3.2.
Acknowledgments
This research was supported by the grant 14.W03.31.0021 of the Ministry of Science and Higher Education of the Russian Federation. We also acknowledge the support from the Program of development of M.V. Lomonosov Moscow State University (Leading Scientific School “Physics of stars, relativistic objects and galaxies”) (GL, PA) and from the Väisälä foundation (AC). The authors thank Konstantin Malanchev for the fruitful discussion. AC and PA thank the lecturers and participants of the Astrosoma summer school (http://astrosoma.org) who provided an excellent environment for the development of this paper.
References
 Abarca, D., Kluźniak, W., & Sa̧dowski, A. 2018, MNRAS, 479, 3936 [NASA ADS] [CrossRef] [Google Scholar]
 Abramowicz, M. A., Czerny, B., Lasota, J. P., & Szuszkiewicz, E. 1988, ApJ, 332, 646 [NASA ADS] [CrossRef] [Google Scholar]
 Abramowicz, M. A., Chen, X., Kato, S., Lasota, J.P., & Regev, O. 1995, ApJ, 438, L37 [NASA ADS] [CrossRef] [Google Scholar]
 Abramowicz, M. A., Chen, X. M., Granath, M., & Lasota, J. P. 1996, ApJ, 471, 762 [NASA ADS] [CrossRef] [Google Scholar]
 Bachetti, M., Harrison, F. A., Walton, D. J., et al. 2014, Nature, 514, 202 [NASA ADS] [CrossRef] [Google Scholar]
 Bachetti, M., Grefenstette, B. W., & Walton, D. J. 2018, ATel, 11282 [Google Scholar]
 Basko, M. M., & Sunyaev, R. A. 1976, MNRAS, 175, 395 [NASA ADS] [CrossRef] [Google Scholar]
 Begelman, M. C., & Meier, D. L. 1982, ApJ, 253, 873 [NASA ADS] [CrossRef] [Google Scholar]
 Beloborodov, A. M. 1998, MNRAS, 297, 739 [NASA ADS] [CrossRef] [Google Scholar]
 Bessolaz, N., Zanni, C., Ferreira, J., Keppens, R., & Bouvier, J. 2008, A&A, 478, 155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Binder, B., Levesque, E. M., & DornWallenstein, T. 2018, ApJ, 863, 141 [NASA ADS] [CrossRef] [Google Scholar]
 Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883 [NASA ADS] [CrossRef] [Google Scholar]
 Carpano, S., Haberl, F., Maitra, C., & Vasilopoulos, G. 2018, MNRAS, 476, L45 [NASA ADS] [CrossRef] [Google Scholar]
 Chashkina, A., Abolmasov, P., & Poutanen, J. 2017, MNRAS, 470, 2799 [NASA ADS] [CrossRef] [Google Scholar]
 Eggum, G. E., Coroniti, F. V., & Katz, J. I. 1988, ApJ, 330, 142 [NASA ADS] [CrossRef] [Google Scholar]
 Elsner, R. F., & Lamb, F. K. 1977, ApJ, 215, 897 [NASA ADS] [CrossRef] [Google Scholar]
 Frank, J., King, A., & Raine, D. J. 2002, Accretion Power in Astrophysics, 3rd edn. (Cambridge, UK: Cambridge University Press), 398 [Google Scholar]
 Fukue, J. 2004, PASJ, 56, 569 [NASA ADS] [CrossRef] [Google Scholar]
 Fürst, F., Walton, D. J., Harrison, F. A., et al. 2016, ApJ, 831, L14 [NASA ADS] [CrossRef] [Google Scholar]
 Ghosh, P., Pethick, C. J., & Lamb, F. K. 1977, ApJ, 217, 578 [NASA ADS] [CrossRef] [Google Scholar]
 Grebenev, S. A. 2017, Astron. Lett., 43, 464 [NASA ADS] [CrossRef] [Google Scholar]
 Ichimaru, S. 1977, ApJ, 214, 840 [NASA ADS] [CrossRef] [Google Scholar]
 Illarionov, A. F., & Sunyaev, R. A. 1975, A&A, 39, 185 [NASA ADS] [Google Scholar]
 Israel, G. L., Papitto, A., Esposito, P., et al. 2017a, MNRAS, 466, L48 [NASA ADS] [CrossRef] [Google Scholar]
 Israel, G. L., Belfiore, A., Stella, L., et al. 2017b, Science, 355, 817 [NASA ADS] [CrossRef] [Google Scholar]
 Kaaret, P., Feng, H., & Roberts, T. P. 2017, ARA&A, 55, 303 [Google Scholar]
 Kato, S., Fukue, J., & Mineshige, S. 2008, BlackHole Accretion Disks  Towards a New Paradigm (Kyoto: Kyoto University Press) [Google Scholar]
 Kawashima, T., Mineshige, S., Ohsuga, K., & Ogawa, T. 2016, PASJ, 68, 83 [NASA ADS] [CrossRef] [Google Scholar]
 Kennea, J. A., Lien, A. Y., Krimm, H. A., Cenko, S. B., & Siegel, M. H. 2017, ATel, 10809 [Google Scholar]
 Kitabatake, E., Fukue, J., & Matsumoto, K. 2002, PASJ, 54, 235 [NASA ADS] [CrossRef] [Google Scholar]
 Kluzniak, W., & Rappaport, S. 2007, ApJ, 671, 1990 [NASA ADS] [CrossRef] [Google Scholar]
 Koenigl, A. 1991, ApJ, 370, L39 [NASA ADS] [CrossRef] [Google Scholar]
 Kulkarni, A. K., & Romanova, M. M. 2013, MNRAS, 433, 3048 [NASA ADS] [CrossRef] [Google Scholar]
 Lasota, J.P., Vieira, R. S. S., Sadowski, A., Narayan, R., & Abramowicz, M. A. 2016, A&A, 587, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lipunov, V. M. 1982, Sov. Ast., 26, 54 [NASA ADS] [Google Scholar]
 Lipunov, V. M. 1992, Astrophysics of Neutron Stars (Berlin: SpringerVerlag), 108 [Google Scholar]
 Lipunova, G. V. 1999, Astron. Lett., 25, 508 [NASA ADS] [Google Scholar]
 Long, M., Romanova, M. M., & Lovelace, R. V. E. 2005, ApJ, 634, 1214 [NASA ADS] [CrossRef] [Google Scholar]
 Matt, S., & Pudritz, R. E. 2005a, ApJ, 632, L135 [NASA ADS] [CrossRef] [Google Scholar]
 Matt, S., & Pudritz, R. E. 2005b, MNRAS, 356, 167 [NASA ADS] [CrossRef] [Google Scholar]
 McKinney, J. C., Tchekhovskoy, A., Sadowski, A., & Narayan, R. 2014, MNRAS, 441, 3177 [NASA ADS] [CrossRef] [Google Scholar]
 Mushtukov, A. A., Suleimanov, V. F., Tsygankov, S. S., & Poutanen, J. 2015, MNRAS, 454, 2539 [Google Scholar]
 Mushtukov, A. A., Ingram, A., Middleton, M., Nagirner, D. I., & van der Klis, M. 2019, MNRAS, 484, 687 [NASA ADS] [CrossRef] [Google Scholar]
 Narayan, R., & Yi, I. 1994, ApJ, 428, L13 [NASA ADS] [CrossRef] [Google Scholar]
 Narayan, R., & Yi, I. 1995, ApJ, 452, 710 [NASA ADS] [CrossRef] [Google Scholar]
 Ogawa, T., Mineshige, S., Kawashima, T., Ohsuga, K., & Hashizume, K. 2017, PASJ, 69, 33 [NASA ADS] [CrossRef] [Google Scholar]
 Ohsuga, K., & Mineshige, S. 2011, ApJ, 736, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Ohsuga, K., Mori, M., Nakamoto, T., & Mineshige, S. 2005, ApJ, 628, 368 [NASA ADS] [CrossRef] [Google Scholar]
 Okuda, T., Teresi, V., Toscano, E., & Molteni, D. 2005, MNRAS, 357, 295 [NASA ADS] [CrossRef] [Google Scholar]
 Paczynski, B., & Jaroszynski, M. 1978, Acta Astron., 28, 111 [NASA ADS] [Google Scholar]
 Parfrey, K., Spitkovsky, A., & Beloborodov, A. M. 2016, ApJ, 822, 33 [NASA ADS] [CrossRef] [Google Scholar]
 Pavlovskii, K., & Ivanova, N. 2015, MNRAS, 449, 4415 [NASA ADS] [CrossRef] [Google Scholar]
 Pavlovskii, K., Ivanova, N., Belczynski, K., & Van, K. X. 2017, MNRAS, 465, 2092 [NASA ADS] [CrossRef] [Google Scholar]
 Poutanen, J., Lipunova, G., Fabrika, S., Butkevich, A. G., & Abolmasov, P. 2007, MNRAS, 377, 1187 [NASA ADS] [CrossRef] [Google Scholar]
 Proga, D. 2000, ApJ, 538, 684 [NASA ADS] [CrossRef] [Google Scholar]
 Rappaport, S. A., Fregeau, J. M., & Spruit, H. 2004, ApJ, 606, 436 [NASA ADS] [CrossRef] [Google Scholar]
 Romanova, M. M., & Owocki, S. P. 2015, Space Sci. Rev., 191, 339 [NASA ADS] [CrossRef] [Google Scholar]
 Romanova, M. M., Ustyugova, G. V., Koldoba, A. V., & Lovelace, R. V. E. 2002, ApJ, 578, 420 [NASA ADS] [CrossRef] [Google Scholar]
 Romanova, M. M., Ustyugova, G. V., Koldoba, A. V., Wick, J. V., & Lovelace, R. V. E. 2003, ApJ, 595, 1009 [NASA ADS] [CrossRef] [Google Scholar]
 Sa̧dowski, A., & Narayan, R. 2016, MNRAS, 456, 3929 [NASA ADS] [CrossRef] [Google Scholar]
 Sa̧dowski, A., Narayan, R., McKinney, J. C., & Tchekhovskoy, A. 2014, MNRAS, 439, 503 [NASA ADS] [CrossRef] [Google Scholar]
 Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 [NASA ADS] [Google Scholar]
 Spruit, H. C., & Taam, R. E. 1993, ApJ, 402, 593 [CrossRef] [Google Scholar]
 Townsend, L. J., Kennea, J. A., Coe, M. J., et al. 2017, MNRAS, 471, 3878 [NASA ADS] [CrossRef] [Google Scholar]
 Tsygankov, S. S., Mushtukov, A. A., Suleimanov, V. F., & Poutanen, J. 2016, MNRAS, 457, 1101 [NASA ADS] [CrossRef] [Google Scholar]
 Tsygankov, S. S., Doroshenko, V., Lutovinov, A. A., Mushtukov, A. A., & Poutanen, J. 2017, A&A, 605, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tsygankov, S. S., Doroshenko, V., Mushtukov, A. A., Lutovinov, A. A., & Poutanen, J. 2018, MNRAS, 479, L134 [NASA ADS] [CrossRef] [Google Scholar]
 Ustyugova, G. V., Koldoba, A. V., Romanova, M. M., Chechetkin, V. M., & Lovelace, R. V. E. 1999, ApJ, 516, 221 [NASA ADS] [CrossRef] [Google Scholar]
 Vasilopoulos, G., Haberl, F., Carpano, S., & Maitra, C. 2018, A&A, 620, L12 [NASA ADS] [EDP Sciences] [Google Scholar]
 Villar, V. A., Berger, E., Chornock, R., et al. 2016, ApJ, 830, 11 [NASA ADS] [CrossRef] [Google Scholar]
 Wang, Y.M. 1996, ApJ, 465, L111 [NASA ADS] [CrossRef] [Google Scholar]
 Weng, S.S., Ge, M.Y., Zhao, H.H., et al. 2017, ApJ, 843, 69 [NASA ADS] [CrossRef] [Google Scholar]
 WilsonHodge, C. A., Malacaria, C., Jenke, P. A., et al. 2018, ApJ, 863, 9 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Disc vertical structure
For the vertical structure of the disc, we assume the form ρ = ρ_{c}(1 − x^{2})^{n}, where x = z/H and n can be viewed as an effective vertical polytropic index. The vertical hydrostatic relation
implies a similar law for the vertical pressure profile P = P_{c}(1−x^{2})^{n + 1}, and a scaling relation for the disc thickness
where the surface density and vertically integrated pressure are related to the midplane quantities as
Here
A similar model for the vertical structure was considered by Kato et al. (2008). In the radiationpressuredominated regime T ∝ P^{1/4}, and therefore T ∝ (1 − x^{2})^{(n + 1)/4}. This approximation is not valid near the disc surface because T_{eff} ≠ 0. The vertical radiative energy flux is determined by the vertical radiation diffusion equation,
where ϵ = aT^{4} is the radiation energy density and D = c/(3κρ) is the diffusion coefficient. We denote the total energy release per unit surface area by Q^{+}, and the total energy leaving the two sides of the disc by
The diffusion approximation allows us to connect the effective temperature to the temperature gradient inside the disc as
Taking into account Eqs. (A.3) and (A.4), this expression implies
and hence
The midplane gas pressure is P_{g} = k T_{c}ρ_{c}/m̃. The gastototal pressure ratio then equals
where m̃ is the mean particle mass (for completely ionized gas of solar metallicity ṁ ≈ 0.6m_{p}).
Appendix B: Derivation of advection equations
The advective flux is given by expression (11), which contains the specific (per particle) dimensionless entropy
The radial derivative of the entropy, under our assumptions about the vertical structure, can be written as
where
and
with β being the gastototal pressure ratio introduced in Eq. (A.11). Vertical integration of Eq. (11) allows us to write the advective flux as
where the dimensionless constants are
To shorten the notations in Eq. (13), we also use the following combinations:
Appendix C: Dimensionless notations and equations
Here we list the dimensionless parameters and combinations used throughout the paper. Our notations here are identical to those in CAP. We normalize the NS mass as
The radius and the disc thickness r = R/R_{g} and h = H/R_{g} are measured in units of the gravitational radius R_{g}. The angular frequency is normalized by the local Keplerian frequency as
The characteristic magnetic moments of NSs lie in the range 10^{28}−10^{32} G cm^{3}, hence we normalize μ as
The mass accretion rate is normalized by the Eddington value as
where
It is convenient to express the surface density in the units of the inverse opacity κ^{−1}. This quantity also has the physical meaning of the disc vertical optical depth
The dimensionless version of the vertically integrated tangential stress can be constructed as
For temperatures we use the normalization
where
The inner radius of the disc can be normalized either by the gravitational or by the Alfvén radius
where the dimensionless Alfvén radius is
and
We also introduce the natural time unit
that can be viewed as a Keplerian rotation period at R_{g}, and the dimensionless factor
The physical meaning of χ is the square of the dimensionless speed of sound (c_{s}/c)^{2} corresponding to the characteristic temperature T_{*}.
Here we give all the equations in dimensionless form, as they were used to calculate the disc structure. The angular velocity at the inner boundary can be found from Eq. (1):
where p is the neutron star period in seconds. From Eq. (2) we can find the stress tensor at the boundary of the disc
The ratio of the gas pressure to the total pressure is, according to Eq. (A.11),
The differential Eqs. (15)–(20) in the dimensionless form are
Taking into account the sign of dΩ/dR < 0, we rewrite Eq. (20) as
where the dimensionless versions of coefficients (B.12)–(B.15) are
All Tables
All Figures
Fig. 1. Structure of an accretion disc around a ULXP. For very high mass accretion rates, the inner parts of the disc (inside the spherization radius R_{sph} shown here) enter the superEddington accretion regime. Inside R_{sph} the thin disc model is not applicable. Mass loss in a wind is shown by blue arrows. The accretion column, where most of the energy is released, is shown by yellow cones, and the red wavy lines refer to the radiation of the column that may affect the inner disc pressure balance. The red vertical line marks the effective boundary between the disc and the magnetosphere at R_{in}. 

In the text 
Fig. 2. Relative magnetospheric radius ξ_{eff} given by Eq. (23) is colourcoded and shown by contours on the ṁ_{0}− μ plane. The spin period is p = 10 p_{eq}. The luminosity scale is calculated assuming efficiency η = 0.1, but irradiation effects on the disc structure were ignored. 

In the text 
Fig. 3. Relative magnetospheric radii ξ (dashed blue curve) and ξ_{eff} (green solid) as functions of the disc thickness at the magnetospheric boundary. The dotted red line corresponds to the longperiod asymptotic given by Eq. (57) in CAP. Magnetic moment of the NS was set to μ_{30} = 100. 

In the text 
Fig. 4. Fraction of the angular momentum flux retained in the superEddington disc l_{in}/l_{0} as a function of the fraction of mass reaching the magnetosphere. The dotted black line corresponds to equal ratios (l_{in}/l_{0} = Ṁ_{in}/Ṁ_{0}, as expected if the net angular momentum is constant); the solid green line corresponds to the expected scaling of a growing spherization radius (Eq. (26)). 

In the text 
Fig. 5. Parameter ξ_{eff} as a function of fastness parameter for magnetic dipole moments μ_{30} = 1 and μ_{30} = 10 and mass accretion rates ṁ_{0} = 10^{4} and ṁ_{0} = 5 × 10^{3}. 

In the text 
Fig. 6. Ratio of the spherization radius r_{sph} (in units of gravitational radius R_{g}) to the dimensionless mass accretion rate ṁ_{0} = 10^{4} as a function of the fastness parameter (green dots). The red horizontal line gives the result from Eq. (21) in Poutanen et al. (2007). 

In the text 
Fig. 7. Colourcoded contours of relative correction to the magnetospheric radius caused by irradiation ξ(η = 0.1)/ξ(η = 0)−1 shown in the ṁ_{0} − μ plane. 

In the text 
Fig. 8. Relative disc thickness as a function of the radius for a model with μ = 10^{30} G cm^{3}, ṁ_{0} = 3 × 10^{3}, and p = 0.67 s. Our results are shown by the red dotted curve, whereas the solid green and blue curves correspond to the asymptotics for the zones B and A of the standard disc, respectively. 

In the text 
Fig. 9. Advection effects for the model shown in Fig. 8. The advection flux fraction Q_{adv}/Q^{+} and the radiated fraction Q_{rad}/Q^{+} are shown by the red solid and green dashed curves, respectively. 

In the text 
Fig. 10. Cumulative luminosities as functions of radius. The total power generated in the disc is shown by the dashed blue curve. The solid red and dotted black curves show the radiative and advected luminosities, respectively. The parameters of the model are μ_{30} = 1 and ṁ_{0} = 3000. 

In the text 
Fig. 11. Fraction of the mass accretion rate reaching radius R for a NS with magnetic moment μ_{30} = 1. The lines from top to bottom correspond to different ṁ_{0} in the interval 2000–3500 with steps of 500. 

In the text 
Fig. 12. Magnetospheric radius in units of R_{g} as a function of the accretion rate for a NS with magnetic moment μ_{30} = 1. Parts of the black solid curve with different slopes correspond to the different regimes of accretion near the magnetospheric boundary. Two standard solutions () are plotted with the grey dotted lines: ξ = 0.5 and ξ = 1 (spherically symmetric case). 

In the text 
Fig. 13. Same as Fig. 12, but for different magnetic moments μ = 10^{30}−10^{32} G cm^{3}. 

In the text 
Fig. 14. Slope of the dependence of the magnetospheric radius on the accretion rate for the models shown in Fig. 13. 

In the text 
Fig. 15. Restrictions for the position of NGC 5907 X1 in the ṁ_{0} − μ plane. The red horizontal line corresponds to the mass accretion rate of ṁ_{0} = 6500 (accretion efficiency η ≃ 0.15). The black solid line corresponds to p_{eq} = p = 1.137 s. The dotted spinup line shows the lower limit for ṁ_{in} set by inequality (32) using period derivation ṗ = 5 × 10^{−9} s s^{−1}. We use here M = 1.4 M_{⊙}, R = 10 km, I = 1.5 × 10^{45} g cm^{2}, and p = 1.137 s. 

In the text 
Fig. 16. Restrictions for the position of NGC 300 ULX1 in the ṁ_{0} − μ plane. The red horizontal line corresponds to the luminosity of 4.7 × 10^{39} erg s^{−1} assuming accretion efficiency η ≃ 0.15. The black solid line corresponds to p_{eq} = p = 20 s. The dotted spinup line shows the lower limit for ṁ_{0} set by inequality (32) using period derivation ṗ = 5.5 × 10^{−7} s s^{−1} (Carpano et al. 2018). We use here M = 1.4 M_{⊙}, R = 10 km, I = 1.5 × 10^{45} g cm^{2}, and p = 20 s. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.