Issue 
A&A
Volume 544, August 2012



Article Number  A123  
Number of page(s)  13  
Section  Astrophysical processes  
DOI  https://doi.org/10.1051/00046361/201219065  
Published online  10 August 2012 
Spectral formation in accreting Xray pulsars: bimodal variation of the cyclotron energy with luminosity
^{1} School of Physics, Astronomy, and Computational Sciences, MS 5C3, George Mason University, 4400 University Drive, Fairfax, VA, USA
email: pbecker@gmu.edu
^{2} Institut für Astronomie und Astrophysik, Abt. Astronomie, Universität Tübingen, Sand 1, 72076 Tübingen, Germany
^{3} LeibnizInstitut für Astrophysik Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany
^{4} Department of Electronics and Computer Science, Nagano National College of Technology, 716 Tokuma, 3818550 Nagano, Japan
^{5} ISDC Data Center for Astrophysics, Université de Genève, Chemin d’Ecogia 16, 1290 Versoix, Switzerland
^{6} AIM (UMR 7158 CEA/DSM – CNRS – Université Paris Diderot) Irfu/Service d’Astrophysique, 91191 GifsurYvette, France
^{7} European Space Agency, European Space Astronomy Centre, PO Box 78, 28691 Villanueva de la Cañada, Madrid, Spain
^{8} Space Science Division, Naval Research Laboratory, Washington, DC, USA
^{9} Dr. Karl RemeisObservatory and Erlangen Centre for Astroparticle Physics, Sternwartstr. 7, 96049 Bamberg, Germany
Received: 17 February 2012
Accepted: 22 May 2012
Context. Accretionpowered Xray pulsars exhibit significant variability of the cyclotron resonance scattering feature (CRSF) centroid energy on pulsetopulse timescales, and also on much longer timescales. Two types of spectral variability are observed. For sources in group 1, the CRSF energy is negatively correlated with the variable source luminosity, and for sources in group 2, the opposite behavior is observed. The physical basis for this bimodal behavior is currently not well understood.
Aims. We explore the hypothesis that the accretion dynamics in the group 1 sources is dominated by radiation pressure near the stellar surface, and that Coulomb interactions decelerate the gas to rest in the group 2 sources.
Methods. We derive a new expression for the critical luminosity, L_{crit}, such that radiation pressure decelerates the matter to rest in sources with Xray luminosity L_{X} > L_{crit}. The formula for L_{crit} is based on a simple physical model for the structure of the accretion column in luminous Xray pulsars that takes into account radiative deceleration, the energy dependence of the cyclotron cross section, the thermodynamics of the accreting gas, the dipole structure of the pulsar magnetosphere, and the diffusive escape of radiation through the column walls. We show that for typical neutron star parameters, , where B_{12} is the surface magnetic field strength in units of 10^{12} G.
Results. The formula for the critical luminosity is evaluated for five sources, using the maximum value of the CRSF centroid energy to estimate the surface magnetic field strength B_{12}. The results confirm that the group 1 sources are supercritical (L_{X} > L_{crit}) and the group 2 sources are subcritical (L_{X} < L_{crit}), although the situation is less clear for those highly variable sources that cross over the line L_{X} = L_{crit}. We also explain the variation of the CRSF energy with luminosity as a consequence of the variation of the characteristic emission height. The sign of this dependence is opposite in the supercritical and subcritical cases, hence creating the observed bimodal behavior.
Conclusions. We have developed a new model for the critical luminosity in accretionpowered Xray pulsars that explains the bimodal dependence of the CRSF centroid energy on the Xray luminosity L_{X}. Our model provides a physical basis for the observed variation of the CRSF energy as a function of L_{X} for both the group 1 (supercritical) and the group 2 (subcritical) sources as a result of the variation of the emission height in the column.
Key words: stars: neutron / pulsars: general / radiative transfer / accretion, accretion disks
© ESO, 2012
1. Introduction
Xray binary pulsars (XRBPs) were first observed by Giacconi et al. (1971) and Tananbaum et al. (1972), and now include many of the brightest sources in the Xray sky. In XRBPs, the main sequence companion star transfers matter to the neutron star via Roche lobe overflow, or via a strong stellar wind (Frank et al. 2002). The gas forms an accretion disk around the neutron star, and the material spirals inward until the pressure of the star’s dipole magnetic field becomes comparable to the ram pressure of the matter in the disk. This occurs at the Alfvén radius, located several thousand kilometers out in the accretion disk. The fullyionized accreting plasma is entrained by the magnetic field at the Alfvén radius, and from there the matter is guided through the magnetosphere, forming accretion columns at one or both of the magnetic poles of the star. As the star spins, the inclination angle between the star’s magnetic axis and the axis of the accretion disk changes, and therefore the Alfvén radius varies with the spin period of the star.
Fig. 1 Schematic illustration of the geometry of the accretion column and the variation of the characteristic emission height and emission beam pattern with increasing luminosity L_{X}: a) subcritical, L_{X} < L_{coul} < L_{crit}, pencil beam; b) subcritical, L_{X} ≲ L_{coul} < L_{crit}, pencil beam; c) subcritical, L_{coul} < L_{X} ≲ L_{crit}, intermediate beam pattern; d) supercritical, L_{X} ≳ L_{crit}, fan beam. 
The observed Xray emission is powered by the conversion of gravitational potential energy into kinetic energy, which is then transferred to the radiation field via electron scattering, and ultimately escapes through the walls of the column. The structure of the accretion column is maintained by the strong magnetic field, with a surface strength B_{∗} ≳ 10^{12} G, which results in a magnetic pressure far exceeding that of either the gas or the radiation field. The high incident speed of the freelyfalling plasma, ~ 0.6 c, creates very high temperatures, T ~ 10^{8} K. However, the observed Xray pulsar spectra are highly nonthermal, indicating that the accreting gas is unable to equilibrate during the accretion timescale. In this situation, bulk and thermal Comptonization play key roles in establishing the shape of the observed spectra (Becker & Wolff 2007).
The Xray spectra of many XRBPs contain cyclotron resonant scattering features (CRSFs) appearing as absorption lines. The features are caused by resonant scattering of photons off plasma electrons whose energy is quantized according to their Landau level (see e.g. Trümper et al. 1978; Isenberg et al. 1998; ArayaGóchez & Harding 2000). The CRSFs, when detected, provide a direct measurement of the magnetic field strength at the characteristic altitude of the Xray emission. The energy of the fundamental line and the spacing between the harmonics are approximately proportional to the Bfield strength.
Many XRBPs display Xray spectra that vary significantly with luminosity on timescales much longer than the pulsation period. In particular, variations in the energy of CRSFs as a function of luminosity on timescales of days to months have been detected in V 0332+53 (Mowlavi et al. 2006; Tsygankov et al. 2010), 4U 0115+63 (Mihara et al. 2004; Tsygankov et al. 2007), and Her X1 (Staubert et al. 2007; Vasco et al. 2011). In addition to the longerterm variability, there is also mounting evidence for pulsetopulse variability, in which the spectral hardness, the centroid energy of the CRSF, and the luminosity vary in a correlated way (Klochkov et al. 2011). This shorttimescale variability is likely related to the nonstationarity of the accretion flow, perhaps indicating that the entrainment of matter from the disk onto the magnetic field lines results in filaments and blobs of accreting gas which are then channeled onto the star in a nonuniform way.
The data from both longterm and shortterm (pulsetopulse) observations point to the existence of two types/modes of spectral variability (see discussion by Klochkov et al. 2011). For sources in group 1 (e.g., V 0332+53), the centroid energy of the CRSF is negatively correlated with luminosity. For sources in group 2 (e.g., Her X1), the opposite behavior is observed. The type of spectral variability is likely driven by the mode of accretion, which in turn is determined by the luminosity (see discussion in Staubert et al. 2007). Staubert et al. (2007) and Klochkov et al. (2011) have proposed that for sources in group 1, the deceleration of the flow to rest at the stellar surface is accomplished by the pressure of the radiation field, and in the group 2 sources the deceleration occurs via Coulomb interactions. In this interpretation, a given source falls in one group or the other depending on the value of its Xray luminosity, L_{X}, relative to the critical luminosity, L_{crit}. The hypothesis is that the group 1 sources are supercritical, with Xray luminosity L_{X} > L_{crit}, and the group 2 sources are subcritical (L_{X} < L_{crit}).
The theory predicts that sources in their supercritical state should display a negative correlation between the luminosity and the cyclotron energy, while sources in the subcritical state should display the reverse behavior. Geometrically, the variation of the CRSF energy with luminosity is connected with variation of the characteristic emission height, which is the altitude in the accretion column where the cyclotron absorption feature is imprinted on the observed spectrum. The variation of the emission height as a function of luminosity in the subcritical and supercritical cases is indicated schematically in Fig. 1.
The general picture described above provides a qualitative basis for the interpretation of the observed correlated variation of the CRSF centroid energy with Xray luminosity in some XRBPs. However, in order to obtain a quantitative understanding of these observations, one must develop a more detailed physical model for the critical luminosity, and for the dependence of the CRSF energy on the luminosity in the subcritical and supercritical regimes. The first goal of this paper is to derive a new expression for the critical luminosity, taking into account the magnetospheric connection between the radius of the accretion column and the Alfvén radius in the disk, and the energy and angle dependence of the cyclotron scattering cross section. The second goal is to examine the dependence of the CRSF centroid energy on the luminosity in the subcritical and supercritical sources.
The remainder of the paper is organized as follows. In Sect. 2, we obtain a fundamental expression for the critical luminosity that depends on the stellar mass, radius, and surface magnetic field strength. In Sect. 3 we develop simple physical models for the variation of the characteristic emission height as a function of the luminosity for subcritical and supercritical sources. In Sect. 4, we evaluate the critical luminosity based on measurements of the CRSF centroid energy for several XRBPs. We use our subcritical and supercritical models for the variation of the emission height to predict the variation of the CRSF energy as a function of luminosity. The predicted spectral variability is compared with the observational data for each source. We discuss our results and draw conclusions in Sect. 5.
2. Critical luminosity
2.1. Eddington luminosity and radiative deceleration
First we recall the definition of the standard Eddington luminosity, L_{Edd}, for spherically symmetric accretion onto a central mass M_{∗}. If the accreting gas is fullyionized hydrogen, we obtain (1)where σ_{T} is the Thomson cross section, m_{p} is the proton mass, c is the speed of light, and G is the gravitational constant. When the Xray luminosity L_{X} = L_{Edd}, the rate at which momentum is transferred to the gas via Compton scattering balances the gravitational force. Hence if L_{X} > L_{Edd}, the net force is in the outward direction and the gas decelerates as it falls toward the central mass.
We need to make two adjustments to Eq. (1) in order to compute the effective Eddington limit, , appropriate for treating Xray pulsar accretion columns. The first adjustment is to replace the Thomson cross section σ_{T} with σ_{}, which represents the mean scattering cross section for photons propagating parallel to the magnetic field axis. The second adjustment is to take the geometry of the accretion flow into account by reducing the luminosity by the ratio of the column crosssectional area divided by the surface area of the star. Employing these corrections yields for the effective Eddington limit (2)where R_{∗} is the stellar radius and r_{0} denotes the radius of the accretion column, which we assume to have an approximately cylindrical geometry.
Based on Eq. (2), Basko & Sunyaev (1976) concluded that for Xray luminosities L_{X} ≳ 10^{36} erg s^{1}, the incident, freelyfalling gas is decelerated by a vertical flux of radiation that is locally superEddington. The scattering of the incident radiation removes kinetic energy from the electrons (and from the protons via Coulomb coupling), thereby decelerating the gas. Although the radiation flux seen by the gas inside the column is propagating upward, the Xrays that ultimately carry away the kinetic energy actually escape through the walls of the column, rather than the top, unless the luminosity L_{X} ≲ 10^{36} erg s^{1}.
The deceleration of the gas begins when the freelyfalling material encounters a radiationdominated shock whose height above the star increases with increasing luminosity, reaching an altitude of several kilometers for L_{X} ~ 10^{37−38} erg s^{1}. As the gas passes through the shock, the accretion velocity is reduced by a factor of ~7. Even though the radiation flux inside the column is superEddington, the height of the radiationdominated shock is stable. This reflects the fact that the shock is a wave structure, and is not composed of a fixed population of particles. Matter moves through the shock and decelerates, but the shock height remains fixed, unless the luminosity changes. Unlike a classical gasmediated shock, the radiationdominated shock is not discontinuous, and instead has a thickness that is a few times larger than the meanfree path for radiation scattering.
Below the radiationdominated shock, the matter is further decelerated in the hydrostatic “sinking regime”, in which the remaining momentum is transferred to the radiation field and radiated away through the column walls (Basko & Sunyaev 1976). The specific mechanism accomplishing the final deceleration to rest at the stellar surface in the sinking region depends on the luminosity of the accretion flow (see Fig. 1). At the highest luminosities, L_{X} ~ 10^{37 − 38} erg s^{1}, the radiation field accomplishes the deceleration all the way down to the stellar surface (Basko & Sunyaev 1976). At intermediate luminosities L_{X} ~ 10^{36 − 37} erg s^{1}, the final phase of deceleration may occur via Coulomb breaking in a plasma cloud just above the stellar surface (Nelson et al. 1993). It is expected that at very low luminosities, L_{X} ≲ 10^{34 − 35} erg s^{1}, there is no radiationdominated shock at all, and the material passes through a conventional gasmediated shock at altitude z = h_{g} before striking the stellar surface (Langer & Rappaport 1982).
The angular pattern of the emitted radiation also depends on the luminosity (see Fig. 1). In highluminosity sources (L_{X} ~ 10^{37 − 38} erg s^{1}), the emitted radiation primarily escapes through the column walls in the sinking region, forming a “fan beam” (Davidson 1973). For lowluminosity sources (L_{X} ≲ 10^{35} erg s^{1}), the emission escapes from the top of the column, forming a “pencil beam” (Burnard et al. 1991; Nelson et al. 1993). In the intermediate range, L_{X} ≲ 10^{35 − 37} erg s^{1}, the emission pattern may be a hybrid combination of these two types (Blum & Kraus 2000).
Focusing on the highluminosity case for now, we can estimate the luminosity required to decelerate the gas to rest at the stellar surface by considering the physical processes occurring in the sinking region below the radiationdominated shock. The accreting matter approaches the top of the shock with the incident freefall velocity, which we approximate using the value at the stellar surface, (3)Advection is dominant over diffusion in the shock, and therefore very little radiation energy escapes through the walls of the accretion column in the vicinity of the shock (Burnard et al. 1991). Hence the shock jump conditions are well approximated by the standard RankineHugoniot relations for a gas with adiabatic index γ = 4/3 (Basko & Sunyaev 1976). In this case, the matter leaves the shock with the postshock velocity (4)where we have assumed that the shock is strong, which is reasonable in the luminous sources (Becker 1998). If the altitude of the radiationdominated shock above the stellar surface is H, and the gas decelerates at a constant rate a from the postshock velocity v_{ff}/7 to rest at the stellar surface in the dynamical time t_{dyn}, then we can write the simple kinematical relations (5)Upon elimination of t_{dyn}, we obtain for the required upward acceleration (6)Since the effective gravity is reduced by the pressure of the radiation field, the net acceleration can also be related to the luminosity L_{X} via (7)Setting Eqs. (6) and (7) equal and solving for L_{X} yields (8)Substituting for using Eq. (2), we obtain for the critical luminosity (9)Our goal is to express the parameters r_{0}, σ_{}, and H appearing on the righthand side of Eq. (9) in terms of observable quantities.
2.2. Radiationdominated shock height
The altitude, H, of the radiationdominated shock can be estimated by considering the relationship between the dynamical timescale for deceleration, t_{dyn}, and the photon escape timescale, t_{diff}, which is the mean time it takes the photons to diffuse through the walls of the accretion column. In order for the gas to come to rest at the stellar surface, these two timescales must be comparable in the sinking region below the shock, which is a general property of accretion flows onto white dwarf stars and neutron stars (e.g., Imamura et al. 1987). Combining Eqs. (4) and (5), we obtain for the dynamical time (10)The escape timescale for the photons to diffuse through the column walls is estimated by writing (11)where n_{e} is the electron number density, is the photon diffusion velocity perpendicular to the column axis, and τ_{ ⊥ } is the perpendicular optical thickness. The Thomson cross section σ_{T} is appropriate for photons propagating perpendicular to the column axis (Wang & Frank 1981). The electron number density n_{e} appearing in Eq. (11) can be eliminated using the mass conservation relation, (12)where Ṁ is the accretion rate and v is the inflow velocity, defined to be positive. Combining relations, we can express the escape time through the walls as (13)The deceleration in the sinking region begins on the downstream side of the shock, and therefore we set v = v_{ps} in Eq. (13) and equate t_{dyn} and t_{esc} to obtain (14)which is essentially the same result obtained by Burnard et al. (1991). Expressing the accretion rate in terms of the luminosity using the relation (15)yields the equivalent expression (16)This confirms that the shock is located a few kilometers above the stellar surface in the luminous sources with L_{X} ~ 10^{37−38} erg s^{1} (Basko & Sunyaev 1976). It follows that R_{∗}/(49H) ≪ 1 for sources close to or above the critical luminosity, and therefore Eq. (9) reduces to (17)in agreement with Burnard et al. (1991). Note that in this limit, the critical luminosity simply reduces to the effective Eddington value given by Eq. (2).
2.3. Connection between column radius and Alfvén radius
In this section, we wish to relate the critical luminosity L_{crit} in Eq. (17) to the magnetic field strength at the stellar surface, B_{∗}, by utilizing the connection between the radius of the accretion column, r_{0}, and the Alfvén radius in the disk, R_{A}. The inclination angle between the axis of the accretion disk and the star’s magnetic axis varies with a period equal to the pulsar’s spin period, and this causes an associated variation of the Alfvén radius. However, for our purposes here, an adequate approximation is obtained by using Eq. (13) from Lamb et al. (1973), which yields (18)where the constant Λ = 1 for spherical accretion and Λ < 1 for disk accretion. A variety of uncertainties are folded into Λ, such as the spinaveraging of R_{A} and the possible role of plasma shielding and other magnetospheric effects. Based on Eq. (2) from Harding et al. (1984), Λ can be approximated in the disk application using (19)where α < 1 denotes the ShakuraSunyaev parameter (Shakura & Sunyaev 1973). Although it is difficult to estimate α with any certainty, we generally expect to find α ~ 0.01 − 0.1. We therefore set Λ = 0.1 in our numerical applications.
The Alfvén radius in the disk is connected with the outer surface of the accretion column through the dipole shape of the pulsar magnetosphere. The equation for the shape of the critical field line as a function of the polar angle θ is given by the standard dipole formula (20)Setting the radius R equal to the stellar radius R_{∗} yields for the critical angle at the outer edge of the accretion column (21)Using the smallangle relation θ_{c} ≈ sinθ_{c}, we obtain for the column radius (22)By substituting for the Alfvén radius in Eq. (22) using Eq. (18), we find that the expression for the column radius r_{0} can be rewritten in cgs units as (23)Using Eq. (23) to substitute for r_{0} in Eq. (17), and setting L_{X} = L_{crit}, we obtain for the critical luminosity the new expression (24)where we have also introduced the Thomson cross section σ_{T} as a convenient scaling for the parallel scattering cross section σ_{}. The next step is to evaluate the cross section ratio σ_{}/σ_{T} in terms of observable source parameters.
2.4. Electron scattering cross section for parallel propagation
In typical Xray pulsars, most of the observed radiation is emitted at energies below the cyclotron energy, E_{cyc}. Hence the crosssection ratio σ_{}/σ_{T} can be roughly approximated using (e.g., Arons et al. 1987). (25)where is a measure of the mean energy of the photons propagating parallel to the magnetic field, and the cyclotron energy E_{cyc} is given by (26)The mean photon energy in Eq. (25) can be estimated observationally by integrating the spectrum for a given source. However, in luminous Xray pulsars, most of the observed radiation escapes through the walls of the accretion column, perpendicular to the magnetic field, and therefore the observed spectrum may not be representative of the distribution of photons propagating along the column axis. As an alternative, we can estimate based on the thermal structure of the accreting gas. Specifically, we assume that (27)where T_{eff} is the effective temperature of the radiation in the postshock region, k is Boltzmann’s constant, and the constant w depends on the shape of the spectrum inside the column. We expect that w ~ 1−3, with the lower value corresponding to bremsstrahlung and the upper value to a Planck spectrum. Detailed models suggest that the spectrum inside the column is dominated by bremsstrahlung emission (Becker & Wolff 2007), and therefore we will set w = 1 in the numerical results presented later.
The effective temperature is related to the postshock radiation pressure, P_{r}, via (28)The value of P_{r} can be estimated using the momentum balance relation (29)where ρ_{ff} and v_{ff} (Eq. (3)) denote the upstream mass density and velocity, respectively, just above the shock. Eliminating P_{r} between Eqs. (28) and (29), and substituting for Ṁ and r_{0} using Eqs. (15) and (23), we find that (30)Combining Eqs. (25) − (27), and (30), we obtain for the required cross section ratio the result (31)Using this result to substitute for σ_{   }/σ_{T} in Eq. (24) yields the final expression for the critical luminosity as a function of the surface magnetic field strength, (32)For typical neutron star parameters, with M_{∗} = 1.4 M_{⊙}, R_{∗} = 10 km, Λ = 0.1, and w = 1, we obtain , where B_{12} is the surface magnetic field strength in units of 10^{12} G. In Sect. 4 we plot the critical luminosity and compare it with the variability data for several XRBPs.
3. Variation of emission height
The new expression for the critical luminosity given by Eq. (32) allows us to separate accretionpowered Xray pulsars into subcritical and supercritical categories. Our hypothesis is that in the subcritical sources with variable luminosity L_{X}, the cyclotron energy E_{cyc} will exhibit a positive correlation with L_{X}, and in the supercritical sources the reverse behavior will be observed. The observed CRSF is imprinted on the spectrum at the altitude where most of the emitted radiation escapes from the accretion column. In order to quantify the expected behaviors in the subcritical and supercritical regimes, we must therefore explore the expected variation of the emission height as a function of luminosity for both types of sources. The geometry of the super and subcritical sources is depicted schematically in Fig. 1.
3.1. Supercritical sources
In the supercritical sources (luminosity L_{X} ≳ L_{crit}), radiation pressure dominates the flow dynamics all the way to the stellar surface. Inside the radiationdominated shock, the infalling matter begins to decelerate as it first encounters the “cushion” of radiation hovering at the shock altitude (Davidson 1973). At this altitude, there is a local balance between downward advection and upward diffusion of radiation, and therefore the photon distribution is roughly static. Most of the kinetic energy of the accretion flow is carried away by the scattered radiation, which is likely to be beamed in the downward direction due to specialrelativistic aberration (e.g., Ferrigno et al. 2009). Below the shock altitude, the photons are trapped by advection, although they eventually manage to escape by diffusing through the walls of the column. Hence the observed radiation does not escape from the shock altitude H, but rather from a lower altitude.
Our goal here is to estimate the typical altitude, denoted by h_{s}, at which the photons diffuse through the column walls to form the observed Xray spectrum in the supercritical case. We assume that the observed CRSF is imprinted at this altitude, because at higher altitudes the photons have not had enough time to diffuse through the column and escape through the walls. Conversely, at lower altitudes, the increasing density of the gas in the column inhibits the escape of radiation. We therefore expect that the CRSF energy will reflect the cyclotron energy at the altitude z = h_{s} where the final deceleration phase begins.
We can estimate the emission height h_{s} in the supercritical sources by ensuring that all of the kinetic energy is radiated through the walls by the escaping photons in the altitude range 0 < z < h_{s} (Basko & Sunyaev 1976). Working in the frame comoving with the mean vertical velocity of the radiation in the accretion column, we note that the fraction of the radiation escaping through the walls in the comoving time interval dt′ is equal to dt′/t_{esc}, where t_{esc} is the escape time given by Eq. (13). The requirement that all of the radiation escapes by the time the matter reaches the stellar surface is therefore expressed by the integral condition (33)Using Eq. (13) to substitute for t_{esc} yields (34)where the effective velocity for the photon transport is defined by (35)The flow is expected to be almost perfectly “trapped” in the region below z′ = h_{s}, meaning that advection and diffusion are nearly balanced, leaving very little net flux of radiation (Becker 1998). This implies that the effective velocity v_{eff} is much smaller than the flow velocity v. We define the parameter ξ as the ratio of these two velocities, (36)We demonstrate in Appendix A that the value of ξ can be estimated using (37)where ℳ_{∞} denotes the incident (upstream) Mach number of the flow with respect to the radiation sound speed. In the sinking region below the shock, the effective velocity approaches zero as the gas settles onto the stellar surface. The relationship between the upstream Mach number ℳ_{∞} and the Xray luminosity L_{X} is plotted in Fig. 12 from Becker (1998). For the parameter range of interest here, it is sufficient to adopt a constant value for ξ in the range ξ ~ 10^{2} − 10^{3}. The low value for the effective velocity tends to make the emission region more compact in the supercritical sources.
Combining relations, we find that (38)and therefore the altitude of the emission region is given by (39)where the final result follows from Eq. (15). We can also express h_{s} in cgs units using (40)Note that the emission height in the supercritical sources varies in proportion to the luminosity L_{X} (see Fig. 1).
Based on Eqs. (16) and (40), we conclude that (41)and therefore the characteristic height of emission in the supercritical sources is located far below the altitude of the radiationdominated shock. Eq. (40) indicates that the height of the emission region h_{s} scales in proportion to the luminosity L_{X} in the supercritical sources, which is consistent with the observed behavior in the group 1 sources (Klochkov et al. 2011).
3.2. Subcritical sources
In the subcritical sources (luminosity L_{X} ≲ L_{crit}), the matter still passes through a radiationdominated shock, which accomplishes the initial deceleration, but the pressure of the radiation is insufficient to bring the matter to rest at the stellar surface (Basko & Sunyaev 1976). In this case, the final stopping occurs via direct Coulomb interactions close to the base of the accretion column (Burnard et al. 1991). Our goal in this section is to estimate the characteristic emission height for the subcritical sources, denoted by h_{c}, which is the altitude at which Coulomb interactions begin to decelerate the plasma to rest. The emission is expected to be concentrated in this region because this is essentially the first opportunity that the radiation inside the column has to diffuse through the walls. At higher altitudes, the radiation is swept along by advection, and there is not enough time for the photons to escape. At lower altitudes, the radiation is trapped in the column due to the increasing density, and therefore the emission through the walls tapers off as z → 0. Hence we expect that in the subcritical case, the CRSF energy is imprinted at the altitude z = h_{c} where the strong Coulomb deceleration begins.
The Thomson optical depth, τ_{∗}, required to stop the flow via Coulomb interactions can be estimated in the typical pulsar magnetic field regime using Eq. (3.34) from Nelson et al. (1993) to write (42)where the maximum excited Landau level, n_{max}, is given by (43)A summary of the derivation leading up to Eq. (42) is provided in Appendix B. Adopting typical Xray pulsar parameters, we find that τ_{∗} ~ 20, which is the value utilized in our numerical examples.
We can use Eq. (42) to estimate the emission height in the subcritical sources, h_{c}, as follows. The Thomson depth τ as a function of the altitude z measured from the stellar surface is computed using (44)where ρ = n_{e}m_{p} is the mass density, given by (see Eq. (12)) (45)Assuming that the gas decelerates uniformly in the Coulomb stopping region (starting at altitude h_{c}) from the postshock velocity v_{ps}, we find that the required deceleration is given by (cf. Eq. (6)). The velocity profile v(z) associated with the constant deceleration a is computed using the standard kinematical relations (46)where the negative sign appears because we have defined v and a to be positive. Setting v_{ps} = v_{ff}/7 (see Eq. (4)), we obtain for the deceleration . Upon integration of Eq. (46), we therefore find that the velocity profile in the Coulomb stopping region is given by (47)where v_{ff} is evaluated using Eq. (3).
Using Eq. (47) to substitute for v(z) in Eq. (45) and carrying out the integration in Eq. (44), we obtain for the optical depth profile (48)Finally, setting z = h_{c} and τ = τ_{∗}, we find that the Thomson optical depth required for Coulomb stopping is given by (49)which can be rearranged to obtain for the subcritical emission height (50)Substituting for Ṁ and r_{0} using Eqs. (15) and (23), respectively, yields the equivalent cgs expression (51)This result indicates that the emission height in the subcritical sources decreases with increasing luminosity, which is consistent with the behavior observed in the group 2 sources (Staubert et al. 2007). As indicated in Fig. 1, in the subcritical case, a decrease in the luminosity causes the beam pattern to transition from a pure fan configuration to a hybrid pattern that includes a pencil component. At very low luminosities, the shock essentially sits on the stellar surface, and the emission occurs via the pencil component only (Burnard et al. 1991; Nelson et al. 1993).
Fig. 2 Variability of the cyclotron line energy with luminosity for different sources. The blue dashed line represents the critical luminosity, plotted by setting L_{X} = L_{crit} and E_{cyc} = E_{∗}, where L_{crit} is evaluated using Eq. (55). The red dashed line represents the Coulomb stopping luminosity, plotted by setting L_{X} = L_{coul} and E_{cyc} = E_{∗}, where L_{coul} is evaluated using Eq. (59). Left: observations on a pulsetopulse timescale. Right: observations on longer timescales. 
Source sample characteristics.
3.3. Transition from Coulomb stopping to gas shock
Equation (51) gives the height of the emission zone in the subcritical case, based on the assumption that the final deceleration to rest at the stellar surface occurs via Coulomb interactions. We can estimate the minimum luminosity (or accretion rate) such that Coulomb interactions are capable of stopping the flow by requiring that h_{c} < H, where H is the altitude of the radiationdominated shock. If this condition is violated, then the incident velocity of the gas entering the Coulomb deceleration region becomes essentially the full freefall velocity, v_{ff}, rather than the reduced postshock velocity, v_{ps} = v_{ff}/7. The associated drop in the density causes the Thomson depth τ to fall below the value τ_{∗} ~ 20 required for the gas to be effectively stopped via Coulomb interactions. It is not completely clear what happens in this case, but we expect that the final phase of deceleration will occur via passage through a gasmediated shock near the stellar surface (Langer & Rappaport 1982).
By combining Eqs. (14), (15) and (50), we can show that the condition h_{c} < H implies that (52)Substituting for r_{0} using Eq. (23) and rearranging, we obtain (53)where (54)For luminosities L_{X} ≲ L_{coul}, we expect that the characteristic emission height settles down onto the stellar surface. For very low luminosities, L_{X} ≲ 10^{34−35} erg s^{1}, the radiationdominated shock and the Coulomb atmosphere both dissipate, and the matter strikes the stellar surface after passing through a gasmediated shock (Langer & Rappaport 1982), as indicated in Fig. 1.
Fig. 3 Same as Fig. 2 except with the xaxis rescaled as L_{X}/L_{crit} and the yaxis rescaled as E_{cyc}/E_{∗}. The values of E_{∗} and L_{crit} used for each source are listed in Table 2. The solid curves associated with each group of data represent the theoretical predictions for the variation of the CRSF energy E_{cyc} as a function of L_{X} for each source, computed using Eq. (58), with the emission height h evaluated using Eq. (40) for the supercritical sources and Eq. (51) for the subcritical sources. 
Theoretical parameters for each source, based on analysis of the pulsetopulse variability and the longerterm variability.
4. Applications
Our final result for the critical luminosity as a function of the surface magnetic field strength B_{∗} is given by Eq. (32), which can be rewritten as (55)where (56)denotes the surface value of the CRSF energy. This relation is indicated by the dashed blue line in Fig. 2. Sources to the right of this line are radiating supercritically, and consequently radiation pressure accomplishes the deceleration all the way down to the stellar surface. For sources to the left of this line, the final deceleration occurs via Coulomb interactions.
It is now interesting to compute the critical luminosity for a number of sources and to compare it with the observed variability of E_{cyc} as a function of L_{X} for subcritical and supercritical cases. We have selected for this study the sources for which the behavior of the cyclotron line energy with luminosity has been studied sufficiently well, including both the variations on long timescales (days to months) and short timescales (“pulsetopulse”). Our source sample is listed in Table 1, where we also summarize the reported E_{cyc} − L_{X} behavior (positive/negative correlation or no dependence) for each source. The two panels in Fig. 2 depict the variability of E_{cyc} with luminosity both on “pulsetopulse” (left) and longer (right) timescales for all of the sources considered here. The corresponding references for the data are given in Table 1. The data for the longerterm variability are based on studies of the phaseaveraged emission conducted either during outbursts of the transient sources (Tsygankov et al. 2007, 2010), or following the longterm variation of the emission from the persistent sources (Staubert et al. 2007). For the pulsetopulse variability we refer to the results of Klochkov et al. (2011) who have shown that for a set of pulsars the cyclotron line energy varies with the amplitude of individual pulses. This amplitude most probably reflects the instantaneous mass accretion rate.
The luminosities for A 0535+26, 4 U0115+63 and Her X1 were calculated by integrating the flux of each source over a nearly identical energy range ΔE (see Table 1). For V 0332+53 (Tsygankov et al. 2010) and GX 3041 (Yamamoto et al. 2011) published values for L_{X} and E_{cyc} were taken, where the luminosity has been calculated over a slightly larger energy range [3 − 100] keV. The luminosities were computed using the source distances listed in Table 1. It should be noted that uncertainties in the distances will create additional uncertainties in our computed luminosities, which have not been considered here. However, the typical effective uncertainties of ~ 10% in the calculated luminosities do not strongly affect our results.
One can see that the sign of the correlation between the cyclotron line energy and the luminosity, when detected, is consistent between the longterm and pulsetopulse studies. This suggests that both are reflecting the same underlying physics. It is clear that the variation of E_{cyc} with L_{X} for sources on each side of the critical line L_{X} = L_{crit} is qualitatively consistent with the theoretical predictions, i.e., the correlation between E_{cyc} and L_{X} is negative for the supercritical sources, and positive for the subcritical ones, reflecting the positive variation of the emission height h = h_{s} (Eq. (40)) with increasing L_{X} for the former sources and the negative variation of the emission height h = h_{c} (Eq. (51)) for the latter sources. In the next section, we explore the variation of E_{cyc} as a function of L_{X} in more detail.
4.1. Variation of E_{cyc} with L_{X}
The observed value of E_{cyc} is connected with the local field strength B at the emission altitude h via Eq. (26), where B has the dipole dependence (57)The corresponding variation of E_{cyc} as a function of the emission height h is therefore given by (58)where E_{∗} is the surface value for the cyclotron energy, computed using Eq. (56). Note that we have neglected the variation of the gravitational redshift factor, which is reasonable given the level of approximation employed here (Staubert et al. 2007). In applying Eq. (58) to the supercritical and subcritical cases, we set h = h_{s} (Eq. (40)) and h = h_{c} (Eq. (51)), respectively. Hence Eq. (58) can be used to develop theoretical predictions for the variation of E_{cyc} as a function of L_{X} for supercritical and subcritical sources.
In Fig. 3 we replot the Fig. 2 data based on a rescaling of the vertical and horizontal axes using the parameters E_{∗} and L_{crit}, respectively, which are related to each other via Eq. (55). The value of L_{crit} used for each source is listed in Table 2, along with the corresponding value for the surface cyclotron energy E_{∗} obtained by comparing the theoretical variation of E_{cyc} with the observed variation for each source. In computing L_{crit}, we assume for all sources the same canonical neutron star mass and radius values M_{∗} = 1.4 M_{⊙} and R_{∗} = 10 km, and we set Λ = 0.1 and w = 1 based on the theoretical considerations discussed above. The vertical dashed line marked L_{X} = L_{crit} separates the sources into their sub and supercritical luminosity states. It should be noted that, as M_{∗} and R_{∗} are also input parameters for L_{crit}, the exact positioning of the source data on the xaxis is driven also by the assumed canonical mass and radius values, which might in fact differ between the individual sources.
Figure 3 also includes curves representing the expected theoretical variation of the CRSF energy E_{cyc} as a function of the luminosity L_{X}, computed using Eq. (58) with either h = h_{s} (Eq. (40)) for the supercritical sources or h = h_{c} (Eq. (51)) for the subcritical ones. We again adopt the canonical values M_{∗} = 1.4 M_{⊙}, R_{∗} = 10 km, Λ = 0.1, and w = 1, and we set the Coulomb stopping optical depth using τ_{∗} = 20. The values of E_{∗} and ξ are varied for each source so as to improve the agreement with the data (see Table 2). The parameter ξ is only relevant for the supercritical sources. The values of ξ reported in Table 2 for V 0332+53 and 4U 0115+63 are in the range ξ ~ 10^{2}−10^{3}, as expected for marginally trapped accretion columns (Becker 1998). Based on the results depicted in Fig. 3, we conclude that the agreement between the observed variation of E_{cyc} and that predicted by the theoretical models developed here is reasonably close for both the supercritical and subcritical sources.
A special case is 4U 0115+63. This source was previously observed to show an anticorrelation (e.g. Tsygankov et al. 2007 and references therein). The right panels of our Figs. 2 and 3 (observations on longer timescales) include the results for E_{cyc} and L_{X} obtained by Tsygankov et al. (2007). However, recent studies have shown that the presence of the anticorrelation depends on the choice of the continuum model (Müller et al. 2011).
4.2. Comparison of L_{X} with L_{coul}
It is also interesting to compare the observed luminosities with the minimum value, L_{coul}, required for Coulomb stopping to decelerate the flow to rest at the stellar surface, given by Eq. (54). By combining Eqs. (54) and (56), we find that L_{coul} is related to E_{∗} via (59)This relation is indicated by the dashed red line in Fig. 2. For sources to the left of this line, we expect that the effect of Coulomb interactions is reduced, and the final stopping occurs via passage through a discontinuous, gasmediated shock (Langer & Rappaport 1982). Hence we anticipate that the emission region approaches the stellar surface as L_{X} is reduced below L_{coul}. This interpretation is consistent with the observational data plotted in Figs. 2 and 3, which indicate that E_{cyc} approaches a constant value in the limit L_{X}/L_{coul} ≪ 1.
Following our hypothesis, sources in their supercritical state (L_{X}/L_{crit} ≳ 1) should display a negative correlation between the luminosity and the cyclotron energy while sources in the subcritical state (L_{X}/L_{crit} ≲ 1) should display the reverse behavior. V 0332+53 and Her X1 in their super and subcritical luminosity states nicely fit into that hypothesis on both long and very short timescales. In particular, we note that the model parameters listed in Table 2 for Her X1 are the same for both the longterm and pulsetopulse data. Hence the model developed here for the variation of E_{cyc} as a function of L_{X} provides a robust connection with the data across the entire range of luminosity variation for this source, which implies that the underlying physical mechanism of variation is the same for the longterm and pulsetopulse variations.
The subcritical source A 0535+26 shows no significant trend on long timescales (Fig. 2, right), perhaps due to the fact that the luminosity is always close to or below the Coulomb stopping limit, L_{coul}, which suggests that we should expect little variation of E_{cyc} with L_{X}. On pulsetopulse timescales, A 0535+26 shows some suggestion of a positive correlation (Fig. 2, left), as expected for a subcritical source. The positive correlation suggested by the shorttimescale data may reflect the fact that the luminosity is somewhat higher than L_{coul}, which places it in the subcritical regime according to our theory. However, we note that our model for the variation of E_{cyc} with L_{X} does not work well for A 0535+26, unless we choose an unreasonably large value for E_{∗} relative to the observational data. We believe this reflects the inapplicability of our model in very lowluminosity sources with L_{X} ≲ L_{coul}. For GX 3041, only an indication of a positive E_{cyc} − L_{X} correlation, consistent with its subcritical state, can be seen, as also reported by Yamamoto et al. (2011). No pulsetopulse spectra are yet available for GX 3041.
5. Conclusions
We have examined the hypothesis that observed bimodal variability of the CRSF energy E_{cyc} with luminosity L_{X} in accretionpowered Xray pulsars reflects the dominant mode of accretion, as proposed by Staubert et al. (2007) and Klochkov et al. (2011). In particular, we have derived an expression for the critical luminosity L_{crit} such that the dynamics in the supercritical sources is determined by the radiation pressure, and the dynamics in the subcritical sources is determined by a combination of radiation pressure and either Coulomb interactions or gas pressure. The detailed formula for L_{crit} is given by Eq. (32), but essentially we find that for typical neutron star parameters, , where B_{12} is the surface magnetic field strength in units of 10^{12} G.
The formula for the critical luminosity was evaluated for 5 sources, based on the maximum value for the CRSF centroid energy for each source, E_{∗}, which is treated as a variable parameter in our approach. The results obtained for E_{∗} are close to the maximum observed values for the CRSF energy. The results depicted in Fig. 2 confirm that L_{X} > L_{crit} in the group 1 sources and L_{X} < L_{crit} in the group 2 sources. The situation is less clear for highly variable sources with luminosity L_{X} that crosses over the line L_{X} = L_{crit}, such as V 0332+53 and 4U 0115+63. These two sources display a negative correlation between E_{cyc} and L_{X} in the supercritical regime, as expected, but the trend does not reverse as predicted by our model when L_{X} < L_{crit}. This suggests that these sources may actually always remain supercritical, despite the fact that they cross the vertical line in Fig. 3. This behavior can be accommodated within our model by slightly changing the parameters Λ and w in Eq. (55). Or, alternatively, the behavior of these sources could indicate that their mass and radius values deviate from the canonical values assumed here.
We have developed simple physical models describing the quantitative variation of E_{cyc} with L_{X} in the supercritical and subcritical sources, given by Eq. (58), with the emission height h set using h = h_{s} (Eq. (40)) for the supercritical sources (group 1), and h = h_{c} (Eq. (51)) for the subcritical sources (group 2). In Fig. 3 the formulas we derived for E_{cyc} as a function of L_{X} were compared with the data for the supercritical and subcritical sources. The agreement between the theoretical predictions and the data suggests that our fundamental model for the physical processes operating in these systems is essentially correct. It is important to note that our formulas for the emission heights h_{c} and h_{s} are not equipped to handle the transcritical case with L_{X} ~ L_{crit}, and therefore further work is required in order to treat sources such as V 0332+53 and 4U 0115+63.
Although the observational picture is still not complete, especially for sources with highly variable luminosities, nonetheless we believe that the emerging bimodal paradigm for the variability of the CRSF energy with luminosity in XRBPs supports the hypothesis that we are seeing direct evidence for two different accretion regimes, depending on whether the luminosity is above or below the corresponding value of L_{crit} for the given surface magnetic field strength B_{∗}. The agreement between the theoretical predictions and the observational data is rather surprising given the level of approximation employed here. We therefore believe that these ideas can provide a useful framework for future detailed modeling of XRBPs with variable luminosities, as well as motivation for further observations.
Acknowledgments
The authors gratefully acknowledge generous support provided by the ISSI in Bern, Switzerland, during the course of this work. M.T.W. acknowledges support from the US Office of Naval Research. I.C. acknowledges financial support from the French Space Agency CNES through CNRS. The authors are also grateful for assistance and useful comments from S. Müller, B. West, K. Wolfram, and A. Bodaghee. We also thank S. Tsygankov for providing observational data in digital form. Finally, we are grateful to the referee, Lev Titarchuk, for a careful reading of the manuscript and several insightful suggestions for improvement.
References
 ArayaGóchez, R. A., & Harding, A. K. 2000, ApJ, 544, 1067 [NASA ADS] [CrossRef] [Google Scholar]
 Arons, J., Klein, R. I., & Lea, S. M. 1987, ApJ, 312, 666 [NASA ADS] [CrossRef] [Google Scholar]
 Basko, M. M., & Sunyaev, R. A. 1976, MNRAS, 175, 395 [NASA ADS] [CrossRef] [Google Scholar]
 Becker, P. A. 1998, ApJ, 498, 790 [NASA ADS] [CrossRef] [Google Scholar]
 Becker, P. A., & Wolff, M. T. 2007, ApJ, 654, 435 [NASA ADS] [CrossRef] [Google Scholar]
 Blum, S., & Kraus, U. 2000, ApJ, 529, 968 [NASA ADS] [CrossRef] [Google Scholar]
 Burnard, D., Arons, J., & Klein, R. 1991, ApJ, 367, 575 [CrossRef] [Google Scholar]
 Caballero, I., Kretschmar, P., Santangelo, A., et al. 2007, A&A, 465, L21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Davidson, K. 1973, Nature Phys. Sci., 246, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Ferrigno, C., Becker, P. A., Segreto, A., Mineo, T., & Santangelo, A. 2009, A&A, 498, 825 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Frank, J., King, A. R., & Raine, D. J. 2002, Accretion Power in Astrophysics (Cambridge: Cambridge Univ. Press) [Google Scholar]
 Giacconi, R., Gursky, H., Kellogg, E., et al. 1971, ApJ, 167, L67 [NASA ADS] [CrossRef] [Google Scholar]
 Harding, A. K., Kirk, J. G., Galloway, D. J., Meszaros, P., et al. 1984, ApJ, 278, 369 [NASA ADS] [CrossRef] [Google Scholar]
 Imamura, J. N., Durisen, R. H., Lamb, D. Q., & Weast, G. J. 1987, ApJ, 313, 298 [NASA ADS] [CrossRef] [Google Scholar]
 Isenberg, M., Lamb, D. Q., & Wang, J. C. 1998, ApJ, 493, 154 [NASA ADS] [CrossRef] [Google Scholar]
 Klochkov, D., Staubert, R., Santangelo, A., et al. 2011, A&A, 532, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mihara, T., Makishima, K., & Nagase, F. 2004, ApJ, 610, 390 [NASA ADS] [CrossRef] [Google Scholar]
 Mowlavi, N., Kreykenbohm, I., Shaw, S. E., et al. 2006, A&A, 451, 187 [Google Scholar]
 Müller, S., Obst, M., Kreykenbohm, I., et al. 2011, The Xray Universe 2011, Presentations of the Conference held in Berlin, Germany [Google Scholar]
 Negueruela, I., & Okazaki, A. T. 2001, A&A, 369, 108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Negueruela, I., Roche, P., Fabregat, J., & Coe, M. J. 1999, MNRAS, 307, 695 [NASA ADS] [CrossRef] [Google Scholar]
 Nelson, R., Salpeter, E., & Wassermann, I. 1993, ApJ, 418, 874 [NASA ADS] [CrossRef] [Google Scholar]
 Lamb, F. K., Pethick, C. J., & Pines, D. 1973, ApJ, 184, 271 [NASA ADS] [CrossRef] [Google Scholar]
 Langer, S. H., & Rappaport, S. 1982, ApJ, 257, 733 [NASA ADS] [CrossRef] [Google Scholar]
 Parkes, G. E., Murdin, P. G., & Mason, K. O. 1980, MNRAS, 190, 537 [NASA ADS] [CrossRef] [Google Scholar]
 Reynolds, A. P., Quaintrell, H., Still, M. D., et al. 1997, MNRAS, 288, 43 [NASA ADS] [CrossRef] [Google Scholar]
 Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 [NASA ADS] [Google Scholar]
 Staubert, R., Shakura, N. I., Postnov, K., et al. 2007, A&A, 465, L25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Steele, I. A., Negueruela, I., Coe, M. J., & Roche, P. 1998, MNRAS, 297, L5 [NASA ADS] [CrossRef] [Google Scholar]
 Tananbaum, H., Gursky, H., Kellogg, E. M., et al. 1972, ApJ, 174, L143 [NASA ADS] [CrossRef] [Google Scholar]
 Trümper, J., Pietsch, W., Reppin, C., et al. 1978, ApJ, 219, L105 [NASA ADS] [CrossRef] [Google Scholar]
 Tsygankov, S. S., Lutovinov, A. A., Churazov, E. M., & Sunyaev, R. A. 2007, Astron. Lett., 33, 368 [NASA ADS] [CrossRef] [Google Scholar]
 Tsygankov, S. S., Lutovinov, A. A., & Serber, A. V. 2010, MNRAS, 401, 1628 [NASA ADS] [CrossRef] [Google Scholar]
 Vasco, D., Klochkov, D., & Staubert, R. 2011, A&A, 532, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Wang, Y.M., & Frank, J. 1981, A&A, 93, 25 [Google Scholar]
 Yamamoto, T., Sugizaki, M., Mihara, T., et al. 2011, PASJ, 63, 751 [Google Scholar]
Appendix A: Effective velocity
The characteristic emission height in the supercritical sources treated in Sect. 3.1 is estimated by comparing the vertical transport time for the radiation with the mean escape time for the photons to diffuse out through the walls of the accretion column. This requires a determination of the “effective velocity”, v_{eff}, which is defined as the net photon transport velocity in the vertical direction, taking into account the competing effects of advection and diffusion. The former process tends to drag photons downward toward the stellar surface, and the latter process tends to transport photons in the opposite direction, upward through the accretion column. The relationship between the effective velocity v_{eff} and the flow velocity v is expressed by the dimensionless parameter ξ, defined by (A.1)In “trapped” regions of the flow, vertical advection and diffusion are nearly balanced, and consequently v_{eff} ≪ v and ξ ≪ 1 (Becker 1998). Trapping tends to occur in the lower, hydrostatic region of the accretion column in the supercritical sources. In this situation, the photons tend to “hover” in a small altitude range until they escape through the walls of the accretion column. Hence this effect reduces the size of the emission region in the supercritical sources treated in Sect. 3.1.
The gas enters the top of the accretion column moving supersonically, but it must come to rest at the stellar surface, and therefore the flow passes through a sonic point somewhere in the column. The sonic point is located in the middle of the radiationdominated shock, where the flow begins to decelerate from the incident freefall velocity v_{ff} (Eq. (3)). Hence the sonic point represents the top of the hydrostatic sinking region, where the radiation tends to escape, and we will therefore estimate the value of ξ using conditions there.
In order to determine the flow velocity at the sonic point, it is useful to consider the conservation of mass and momentum in the hydrostatic region of the column. We have (A.2)where J and I denote the fluxes of mass and momentum, respectively. These two fluxes are conserved in the roughly cylindrical, hydrostatic portion of the accretion column.
We can use Eqs. (A.2) to obtain a relationship between the flow velocity v and the radiation Mach number, ℳ, defined by (A.3)where a denotes the radiation sound speed, given by (A.4)The result obtained is (A.5)where v_{c} denotes the flow velocity at the radiation sonic point, where ℳ = 1.
The value of ξ can be estimated by examining the vertical propagating of the photons in a radiationdominated accretion column described by the exact dynamical solution obtained by Basko & Sunyaev (1976) and Becker (1998). This solution assumes a cylindrical geometry in the hydrostatic lower region of the accretion column. The total radiation energy flux in the vertical direction is given by (A.6)where P_{r} is the radiation pressure, and the first and second terms on the righthand side represent advection and diffusion, respectively. We define v_{eff} by writing (A.7)so that v_{eff} represents the “effective” bulk velocity that would yield the correct energy flux.
By combining Eqs. (A.6) and (A.7), we obtain (A.8)or, equivalently, (A.9)where we have eliminated the electron number density using the relation J = n_{e}m_{p}v. The pressure P_{r} can be expressed in terms of the flow velocity v by using Eqs. (A.2) to write (A.10)Using this relation to substitute for the pressure P_{r} in Eq. (A.9) yields (A.11)The exact solution for the flow velocity profile v(z) in a cylindrical accretion column is given by (Basko & Sunyaev 1976; Becker 1998) (A.12)where (A.13)denotes the value of the dimensionless total energy flux at the sonic point, and (A.14)is the altitude of the sonic point above the stellar surface.
Combining Eqs. (A.11) − (A.14), we obtain after some algebra (A.15)In particular, at the sonic point (z = z_{st}), we have v = v_{c}, and therefore our result for ξ reduces to (A.16)The dimensionless energy flux is related to the incident Mach number of the flow, ℳ_{∞}, via (Becker 1998) (A.17)Combining Eqs. (A.16) and (A.17) yields for the value of ξ at the radiation sonic point (A.18)We use this relation in Sect. 3.1, where we estimate the height of the emission region in the supercritical sources.
Appendix B: Coulomb stopping depth
In the subcritical sources treated in Sect. 3.2, radiation pressure is insufficient to decelerate the flow to rest at the stellar surface. In this regime, the final deceleration likely occurs via Coulomb interactions between the infalling plasma and the mound of dense gas that has built up just above the stellar surface. We can estimate the vertical extent of this region, and therefore obtain an approximation of the characteristic emission altitude in the subcritical sources, by computing the Thomson optical depth, τ, measured from the stellar surface, and setting it equal to the value required to stop the flow, denoted by τ_{∗}. Nelson et al. (1993) carried out a detailed calculation of τ_{∗}, and the final result is presented in their equation Eq. (3.34). In this section we summarize the derivation.
In a magnetized pulsar accretion column, the Coulomb stopping of the gas occurs via coupling between infalling protons and stationary electrons in the mound. The corresponding rate of change of the proton kinetic energy, (B.1)is given by Eq. (3.31) from Nelson et al. (1993), which states that (B.2)where lnΛ_{c} is the Coulomb logarithm and v is the velocity of the protons. Note that the righthand side of Eq. (B.2) is positive in our sign convention since the value of E_{p} decreases with decreasing altitude.
Using Eq. (B.1) to substitute for v yields the equivalent form (B.3)We can transform from the altitude dz to the Thomson depth dτ using (B.4)which yields (B.5)Treating the Coulomb logarithm as a constant and integrating with respect to τ, we obtain the solution (B.6)where the stopping depth, τ_{∗}, is defined by (B.7)and the incident proton kinetic energy, E_{0}, is equal to the freefall value, (B.8)Substituting for the Thomson cross section, σ_{T}, in Eq. (B.7) using (B.9)yields the equivalent result (B.10)In magnetized pulsar accretion columns, with discrete Landau levels, the Coulomb logarithm is given by Eq. (3.32) from Nelson et al. (1993), which states that (B.11)where the maximum excited Landau level, n_{max}, is given by (B.12)Combining Eqs. (B.10) and (B.11) and substituting for v_{ff} using Eq. (3) gives the final result, (B.13)in agreement with Eq. (3.34) from Nelson et al. (1993). For typical Xray pulsar parameters, we obtain τ_{∗} ~ 20, and this is the value utilized in computing the characteristic emission height in the subcritical sources in Sect. 3.2.
All Tables
Theoretical parameters for each source, based on analysis of the pulsetopulse variability and the longerterm variability.
All Figures
Fig. 1 Schematic illustration of the geometry of the accretion column and the variation of the characteristic emission height and emission beam pattern with increasing luminosity L_{X}: a) subcritical, L_{X} < L_{coul} < L_{crit}, pencil beam; b) subcritical, L_{X} ≲ L_{coul} < L_{crit}, pencil beam; c) subcritical, L_{coul} < L_{X} ≲ L_{crit}, intermediate beam pattern; d) supercritical, L_{X} ≳ L_{crit}, fan beam. 

In the text 
Fig. 2 Variability of the cyclotron line energy with luminosity for different sources. The blue dashed line represents the critical luminosity, plotted by setting L_{X} = L_{crit} and E_{cyc} = E_{∗}, where L_{crit} is evaluated using Eq. (55). The red dashed line represents the Coulomb stopping luminosity, plotted by setting L_{X} = L_{coul} and E_{cyc} = E_{∗}, where L_{coul} is evaluated using Eq. (59). Left: observations on a pulsetopulse timescale. Right: observations on longer timescales. 

In the text 
Fig. 3 Same as Fig. 2 except with the xaxis rescaled as L_{X}/L_{crit} and the yaxis rescaled as E_{cyc}/E_{∗}. The values of E_{∗} and L_{crit} used for each source are listed in Table 2. The solid curves associated with each group of data represent the theoretical predictions for the variation of the CRSF energy E_{cyc} as a function of L_{X} for each source, computed using Eq. (58), with the emission height h evaluated using Eq. (40) for the supercritical sources and Eq. (51) for the subcritical sources. 

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.