Radiation-dominated polar emitting region of an accreting X-ray pulsar -- I. Polarization- and spectrum-dependent structure, and the emergent continuum

1. Problem Statement

The paper addresses the complex physics of accretion columns in X-ray pulsars, specifically in the radiation-dominated regime (luminosities $L \gtrsim 10^{36}–10^{37}$ erg s$^{-1}$). In these systems, matter from a companion star falls onto the magnetic poles of a neutron star, forming a column where radiation pressure exceeds gas pressure.

Key challenges in modeling this regime include:

• Non-linear coupling: The dynamics of the flow, the thermal structure of the plasma, and the radiative transfer are tightly coupled.
• Polarization dependence: Strong magnetic fields make the plasma birefringent, supporting two normal modes (ordinary and extraordinary) with different scattering cross-sections.
• Shock physics: The formation of a radiation-dominated shock wave and the subsequent "settling mound" requires accurate treatment of bulk motion effects (Comptonization) and induced scattering.
• Observational discrepancies: Existing models often fail to explain the correlation between cyclotron resonant scattering features (CRSF) and X-ray luminosity, or the specific shape of high-energy spectral tails.

2. Methodology

The author develops a self-consistent, three-dimensional (geometrically 2D) numerical model of an axisymmetric accretion column.

• Governing Equations:
  • Momentum Equation: Describes the flow deceleration driven by the radiation pressure gradient, neglecting gas pressure and viscosity.
  • Radiative Transfer Equations: Solved in the Fokker-Planck approximation for two polarization modes (ordinary and extraordinary). The equations include:
    • Spatial diffusion and advection.
    • Thermal Comptonization.
    • First-order and Second-order Bulk Comptonization: The latter accounts for the heating of electrons and photon energy gain due to the bulk velocity of the plasma ($m_e v^2/3$ term).
    • Induced Compton Scattering: Crucial for establishing electron temperature in high-density regions.
    • Mode Switching: Transitions between polarization modes.
    • Free-free processes and Cyclotron emission (via electron-proton collisions).
• Geometry: The model uses spherical polar coordinates $(r, \theta)$ appropriate for a dipole magnetic field ($B \propto r^{-3}$). It considers both "filled" funnels and "hollow" columns.
• Electron Temperature: Calculated locally using the Compton equilibrium approximation, modified to include bulk heating effects. A minimum temperature floor ($T^* \approx 1$ keV) is set in pre-shock regions where Compton equilibrium breaks down.
• Numerical Scheme:
  • Time-relaxation method using an explicit finite-difference scheme.
  • Parallelized computation (OpenMP) on high-performance CPUs.
  • Boundary conditions include free-fall velocity at the top, blackbody radiation at the base, and free escape at the sidewalls.

3. Key Contributions

• Self-Consistent Polarization Modeling: Unlike previous 1D or grey-approximation models, this work solves the radiative transfer for both polarization modes simultaneously with the hydrodynamic structure, allowing for frequency- and polarization-dependent scattering cross-sections.
• Inclusion of Induced Compton Scattering: The paper demonstrates that induced scattering is not negligible; it plays a critical role in establishing the electron temperature in the settling mound, ensuring the plasma reaches local thermodynamic equilibrium (LTE).
• Second-Order Bulk Comptonization: The model explicitly incorporates the $m_e v^2/3$ term in the Doppler shift, showing its necessity for generating high-energy spectral tails and heating the shock wave.
• Realistic Shock Structure: The model resolves the shock wave and the settling mound self-consistently, avoiding the assumption of a pre-defined velocity profile or temperature distribution.

4. Key Results

• Column Structure:
  • The height of the accretion column does not scale linearly with the accretion rate ($\dot{M}$) as predicted by grey approximations. Instead, due to frequency-dependent cross-sections, the column height saturates at high accretion rates ($\dot{M}_{17} \gtrsim 3$).
  • The shock wave is steeper and the settling mound is lower in a dipole field compared to cylindrical channel models.
• Thermal Structure:
  • Shock Heating: The electron temperature in the shock wave reaches 10–20 keV (or higher), significantly exceeding the local Compton temperature of the post-shock region. This heating is driven by bulk motion effects (first and second-order Comptonization).
  • Post-Shock Zone: Deep in the settling mound, the electron temperature approaches the Compton temperature ($T_C \sim T_0$), which is close to the blackbody temperature of the base but modified by radiation density.
  • Role of Induced Scattering: Without induced scattering, the model fails to produce self-consistent steady-state solutions for moderate-to-high accretion rates, as the electron temperature would be underestimated.
• Emergent Spectrum and Polarization:
  • Spectral Shape: The total emergent spectrum exhibits a high-energy power-law tail. The extraordinary mode dominates near the spectral peak, while the ordinary mode dominates at very high energies (though it tends toward saturation).
  • Polarization Fraction (PF): The calculated PF ranges from 0.05 to 0.5, depending on energy and accretion rate. At low accretion rates, the ordinary mode spectrum can be close to the extraordinary mode, reducing the net polarization.
  • High-Energy Tails: The high-energy tails are formed primarily by second-order bulk Comptonization and thermal Comptonization in the shock.
• Comparison with Observations:
  • The model successfully reproduces spectra with high-energy humps (20–30 keV) observed in sources like 1A 0535+262 and V 0332+53.
  • It suggests that the negative correlation between CRSF energy and luminosity (observed in some sources) might be explained by distortions in the magnetic field strength within the column rather than just changes in column height.

5. Significance

This work represents a significant advancement in the theoretical understanding of accreting X-ray pulsars:

1. Physical Consistency: It provides a fully self-consistent solution where the flow dynamics, temperature, and radiation field are mutually dependent, moving beyond the "fixed profile" approximations of earlier studies.
2. Mechanism Identification: It confirms that bulk Comptonization (specifically second-order effects) and induced scattering are the dominant mechanisms shaping the thermal structure and high-energy spectrum, rather than simple thermal equilibrium.
3. Observational Interpretation: The results offer a framework to interpret complex observational data, such as the shape of the X-ray continuum and the behavior of cyclotron lines, suggesting that simple 1D models are insufficient for high-luminosity pulsars.
4. Future Directions: The paper highlights the need for future MHD simulations to account for magnetic field distortions and relativistic radiative transfer to fully resolve discrepancies in CRSF behavior and spectral fitting.

In summary, Gornostaev presents a sophisticated numerical framework that reveals the intricate interplay between radiation pressure, magnetic fields, and bulk plasma motion in X-ray pulsars, providing a more accurate physical basis for interpreting high-energy astrophysical observations.