Accretion Column Research Articles

Probing the energy and luminosity-dependent spectro-timing properties of RX J0440.9+4431 with AstroSat

ABSTRACT The Be/X-ray binary pulsar RX J0440.9+4431 went through a giant outburst in December 2022 with a peak flux of $\sim$2.3 Crab in 15–50 keV. We studied the broad-band timing and spectral properties of RX J0440.9+4431 using four AstroSat observations, where the source transited between subcritical and supercritical accretion regimes. Pulsations were detected significantly above 100 keV. The pulse profiles were found to be highly luminosity- and energy-dependent. A significant evolution in the pulse profile shape near the peak of the outburst indicates a possible change in the accretion mode and beaming patterns of RX J0440.9+4431. The rms pulsed fraction was luminosity- and energy-dependent, with a concave-like feature around 20–30 keV. The depth of this feature varied with luminosity, indicating changes in the accretion column height and proportion of reflected photons. The broad-band continuum spectra were best fitted with a two-component Comptonization model with a blackbody component or a two-blackbody component model with a thermal Comptonization component. A quasi-periodic oscillation (QPO) at 60 mHz was detected at a luminosity of $2.6 \times 10^{37}$ erg s$^{-1}$, which evolved into 42 mHz at $1.5 \times 10^{37}$ erg s$^{-1}$. The QPO rms were found to be energy dependent with an overall increasing trend with energy. For the first time, we found the QPO frequency varying with photon energy in an X-ray pulsar, which poses a challenge in explaining the QPO with current models such as the Keplarian and beat frequency model. Hence, more physically motivated models are required to understand the physical mechanism behind the mHz QPOs.

A quantitative explanation of the cyclotron-line variation in accreting magnetic neutron stars of super-critical luminosity

Context. Magnetic neutron stars (NSs) often exhibit a cyclotron resonant scattering feature (CRSF) in their X-ray spectra. Cyclotron lines are believed to be generated in the radiative shock in the accretion column. High-luminosity NSs show a smooth anti-correlation between the cyclotron-line centroid (ECRSF) and X-ray luminosity (LX). Aims. It has been pointed out that the observed ECRSF − LX smooth anti-correlation in high-luminosity NSs is in tension with the theoretically predicted one, if the radiative shock is the site of cyclotron-line formation. The shock height increases approximately linearly with luminosity, while the dipolar magnetic field drops as the cubic power of distance, thereby implying that the cyclotron-line energy ought to decrease significantly when the luminosity increases by, say, an order of magnitude, which is contrary to observations. Since there is no other candidate site for the cyclotron-line formation, we re-examine the predicted rate of change of the cyclotron-line energy with luminosity at the radiative shock, taking a closer look at the Physics involved. Methods. We developed a purely analytical model describing the overall dependence of the observed cyclotron energy centroid on the shock front’s height, including both the relativistic boosting and the gravitational redshift effects in our considerations. The relativistic boosting effect is due to the mildly relativistic motion of the accreting plasma upstream with respect to the shock’s reference frame. Reults. We find that the cyclotron-line energy varies with (a) the shock height due to the dipolar magnetic field, (b) the Doppler boosting between the shock and bulk-motion frames, and (c) the gravitational redshift. We show that the relativistic effects noticeably weaken the predicted ECRSF − LX anti-correlation. We use our model to fit the data of the X-ray source V0332+53 that exhibits a weak negative correlation and demonstrate that the model fits the data impressively well, thereby alleviating the tension between observations and theory. Conclusions. The reported weak anti-correlation between cyclotron-line centroid and X-ray luminosity in the supercritical accretion regime may be explained by the combination of the variation of the magnetic-field strength along the accretion column, the effect of Doppler boosting, and the gravitational redshift. As a result of these effects, the actual magnetic field on the surface of the neutron star may be a factor of ∼2 larger than the naively inferred value from the observed CRSF.

Spectral evolution of RX J0440.9+4431 during the 2022–2023 giant outburst observed with Insight-HXMT

In 2022–2023, the Be/X-ray binary X-ray pulsar RX J0440.9+4431 underwent a Type II giant outburst, reaching a peak luminosity of LX ∼ × 1037 erg s−1. In this work, we utilized Insight-HXMT data to analyze the spectral evolution of RX J0440.9+4431 during the giant outburst. By analyzing the variation in the X-ray spectrum during the outburst using standard phenomenological models, we find that as the luminosity approaches the critical luminosity, the spectrum becomes flatter, with the photon enhancement predominantly concentrated around ∼2 keV and 20–40 keV. The same behavior has also been noted in Type II outbursts from other sources. While the phenomenological models provide good fits to the spectrum, it is sometimes difficult to gain insight into details of the fundamental accretion physics using this approach. Hence, we also analyzed spectra obtained during high and low phases of the outburst using a new, recently developed physics-based theoretical model that allows us to study the variations in the physical parameters during the outburst, such as the temperature, density, and magnetic field strength. Application of the theoretical model reveals that the observed spectrum is dominated by Comptonized bremsstrahlung emission emitted from the column walls in both the high and low states. We show that the spectral flattening observed at high luminosities results from a decrease in the electron temperature, combined with a compactification of the emission zone, which reduces the efficiency of bulk Comptonization. We also demonstrate that when the source is at maximum luminosity, the spectrum tends to harden around the peak of the pulse profile, and we discuss possible theoretical explanations for this behavior. We argue that the totality of the behavior in this source can be explained if the accretion column is in a quasi-critical state at the time of maximum luminosity during the outburst.

Gas Dynamics in the Herbig Ae Star HD 163296 in 2024

We observed the Herbig Ae star HD 163296 over 9 nights in 2024. On 4 of the nights, the He i line at 1.083 μm exhibited a profile rarely seen in the star. Rather than the typical P Cygni profile indicating a net gas outflow, the profile consisted of a deep absorption line with weak emission peaks on either side. These seem to be consistent with either an outflowing bipolar jet, a circular ring of gas orbiting close to the star, and/or accretion columns connecting the inner disk to the star. Given the history of the ejection of Herbi-Haro object in the star’s bipolar jets, these observations should coincide with the very beginning of the ejection event. Coordinated observations with multiple observing techniques over the next decade are warranted.

A dusty magnetospheric stream explaining the light curves of the dipper objects: Finding a new inclination threshold to produce dippers

Context. The so-called “dippers” are young stellar objects that exhibit dimming episodes in their optical light curves. The common interpretation for the occurrence of these dips is that dusty regions periodically or quasi-periodically cross the line of sight toward the object. Aims. We develop a model where we assume that these regions are located at the intersection of the magnetospheric stream with the disk. The stream is fed by gas and dust coming from the disk. As the material follows the magnetic field lines above the disk plane, it forms an opaque screen that partially blocks the stellar emission. The amount of extinction caused by the material crossing the line of sight depends on the abundance and location of the dust along the stream, which depends on the degree of dust evaporation due to the heating by the star. Methods. We run hydrodynamical simulations of dusty accretion streams to produce synthetic dipper light curves for a sample of low-mass young stars still accreting from their disk according to evolutionary models. We compare the distribution of the light curve amplitudes between the synthetic sample and observed samples of dippers from various star-forming regions. Results. Dust evaporation along the accretion column drives the distribution of photometric amplitudes. Our results suggest that most of the observed dippers correspond to systems seen at high inclination. However, dust survival within accretion columns may also produce dippers at lower inclination, down to about 45°. We find that the dust temperature arising from stellar irradiation should be increased by a factor 1.6 to find consistency between the fraction of dippers our model predicts in star-forming regions and the observed fraction of 20–30%. Conclusions. Transient dust survival in accretion columns appear as an alternative (or complementary) mechanism to inner disk warp occultation in order to account for low-inclination dippers in star-forming regions.

Super-Eddington Magnetized Neutron Star Accretion Flows: A Self-similar Analysis

The properties of super-Eddington accretion disks exhibit substantial distinctions from the sub-Eddington ones. In this paper, we investigate the accretion process of a magnetized neutron star (NS) surrounded by a super-Eddington disk. By constructing self-similar solutions for the disk structure, we study in detail an interaction between the NS magnetosphere and the inner region of the disk, revealing that this interaction takes place within a thin boundary layer. The magnetosphere truncation radius is found to be approximately proportional to the Alfvén radius, with a coefficient ranging between 0.34–0.71, influenced by the advection and twisting of a magnetic field, NS rotation, and radiation emitted from an NS accretion column. Under super-Eddington accretion, the NS can readily spin up to become a rapid rotator. The proposed model can be employed to explore the accretion and evolution of NSs in diverse astrophysical contexts, such as ultraluminous X-ray binaries or active galactic nucleus disks.

Neutron star accretion events in AGN discs: mutimessenger implications

ABSTRACT This paper investigates the accretion of neutron stars (NSs) in active galactic nucleus (AGN) accretion discs. We classify potential accretion modes of NSs in AGN discs, proposing a hierarchical model of NS accretion: accretion flow from the Bondi sphere to accretion columns. The accretion of NSs in AGN discs differs from that of BHs, especially within the scale of the NS’s magnetosphere due to its hard surface and magnetic field. As the accretion flow approaches the magnetosphere, the magnetic fields guide the accretion flow to form accretion columns, primarily dominated by neutrinos. While neutrinos generated from single NS accretion may not have observable effects, considering the all-sky background, they contribute to the neutrino background in the sub-MeV energy range comparable to that of supernova explosions. NS accretion may also lead to the generation of mass quadrupole moments, consequently generating gravitational waves (GWs). The GWs, which exhibit characteristic effects like periodic modulations and echoes, could be observed by third-generation GW detectors. The emission of neutrinos and GWs carries away energy and angular momentum brought by accretion, reducing the feedback effect on the AGN disc. This results in an exceptionally high NS accretion rate, leading to a collapse time-scale shorter than the migration-merge time-scale, making it less likely that binary NS mergers originate from AGN discs.

Winds of change: The nuclear and galaxy-scale outflows and the X-ray variability of 2MASS 0918+2117

Context. In this work, we test feedback propagation models on the test case of 2MASS 0918+2117 (2M0918), a z = 0.149 X-ray variable AGN that shows tentative evidence for nuclear ultra-fast outflows (UFOs) in a 2005 XMM-Newton observation. We also investigate whether UFOs can be related to the observed X-ray variability. Aims. We observed 2M0918 with XMM-Newton and NuSTAR in 2020 to confirm the presence of and characterize the UFOs. We performed a kinematic analysis of the publicly available 2005 SDSS optical spectrum to reveal and measure the properties of galaxy-scale ionized outflows. Furthermore, we constructed 20-year-long light curves of observed flux, line-of-sight column density, and intrinsic accretion rate from the spectra of the first four SRG/eROSITA all-sky surveys and archival observations from Chandra and XMM-Newton. Methods. We detect UFOs with v ∼ 0.16c and galaxy-scale ionized outflows with velocities of ∼700 km s−1. We also find that the drastic X-ray variability (factors > 10) can be explained in terms of variable obscuration and variable intrinsic luminosity. Results. Comparing the energetics of the two outflow phases, 2M0918 is consistent with momentum-driven wind propagation. 2M0918 expands the sample of AGN with both UFOs and ionized gas winds from 5 to 6 and brings the sample of AGN hosting multiscale outflows to 19, contributing to a clearer picture of feedback physics. From the variations in accretion rate, column density, and ionization level of the obscuring medium, we propose a scenario that connects obscurers, an accretion enhancement, and the emergence of UFOs.

Return to the Forgotten Ultraluminous X-Ray Source: A Broadband NICER+NuSTAR Study of NGC 4190 ULX-1

We observed the nearby and relatively understudied ultraluminous X-ray source (ULX) NGC 4190 ULX-1 jointly with Neutron Star Interior Composition Explorer (NICER) and NuSTAR to investigate its broadband spectrum, timing properties, and spectral variation over time. We found NGC 4190 ULX-1 to have a hard spectrum characterized by two thermal components (with temperatures ∼0.25 and ∼1.6 keV) and a high-energy excess typical of the ULX population although the spectrum turns over at an unusually low energy. While no pulsations were detected (with pulsed fraction 3σ upper limits of 16% for NICER and 35% for NuSTAR), the source shows significant stochastic variability, and the covariance spectrum indicates the presence of a high-energy cutoff power-law component, potentially indicative of an accretion column. Additionally, when fitting archival XMM-Newton data with a similar model, we find that the luminosity–temperature evolution of the hot thermal component follows the behavior of a super-Eddington slim disk though the expected spectral broadening for such a disk is not seen, suggesting that the inner accretion disk may be truncated by a magnetic field. Therefore, despite the lack of detected pulsations, there is tantalizing evidence for NGC 4190 ULX-1 being a candidate neutron star accretor although further broadband observations will be required to confirm this behavior.

Cyclotron line formation in the radiative shock of an accreting magnetized neutron star

Context. Magnetic neutron stars (NSs) often exhibit a cyclotron resonant scattering feature (CRSF) in their X-ray spectra. Accretion onto their magnetic poles is responsible for the emergence of X-rays, but the site of the CRSF formation is still a puzzle. A promising candidate for high-luminosity sources has always been the radiative shock in the accretion column. Nevertheless, no quantitative calculations of spectral formation at the radiative shock have been performed so far. Aims. It is well accepted that, in the accretion column of a high-luminosity, accreting magnetic NS, a radiative shock is formed. Here, we aim to explore the scenario where the shock is the site of the cyclotron-line formation. We studied spectral formation at the radiative shock and the emergent spectral shape across a wide range of the parameter space and determined which parameters hold an important role in shaping a prominent CRSF. Methods. We developed a Monte Carlo (MC) code based on the forced first collision numerical scheme to conduct radiation transfer simulations at the radiative shock. The seed photons were due to bremsstrahlung and were emitted in the post-shock region. We properly treated bulk-motion Comptonization in the pre-shock region, thermal Comptonization in the post-shock region, and resonant Compton scattering in both regions. We adopted a fully relativistic scheme for the interaction between radiation and electrons, employing an appropriate polarization-averaged differential cross-section. As a result, we calculated the angle- and energy-dependent emergent X-ray spectrum from the radiative shock, focusing on both the CRSF and the X-ray continuum, under diverse conditions. The accretion column was characterized by cylindrical symmetry, and the radiative shock was treated as a mathematical discontinuity. Results. We find that a power law, hard X-ray continuum, and a CRSF are naturally produced by the first-order Fermi energization as the photons criss-cross the shock. The depth and the width of the CRSF depend mainly on the transverse optical depth and the post-shock temperature. We show that the cyclotron-line energy centroid is shifted by ∼(20 − 30)% to lower energies compared to the classical cyclotron energy; this is due to the Doppler boosting between the shock reference frame and the bulk-motion frame. We demonstrate that a “bump” feature arises in the right wing of the CRSF due to the up-scattering of photons by the accreting plasma and extends to higher energies for larger optical depths and post-shock temperatures. Conclusions. We conclude that resonant Compton scattering of photons by electrons in a radiative shock is efficient in producing a power-law X-ray continuum with a high-energy cutoff accompanied by a prominent CRSF. The implications of the Doppler effect on the centroid of the emergent absorption feature must be considered if an accurate determination of the magnetic field strength is desired.

Flickering pulsations in bright X-ray pulsars: the evidence of gravitationally lensed and eclipsed accretion column

ABSTRACT It is expected that the extreme mass accretion rate onto strongly magnetized neutron stars results in the appearance of accretion columns above the stellar surface. For a distant observer, rotation of a star results in periodic variations of X-ray flux. Because the mass accretion rate fluctuates around the average value, the pulse profiles are not stable and demonstrate fluctuations as well. In the case of bright X-ray pulsars, however, pulse fluctuations are not solely attributed to variations in the mass accretion rate. They are also influenced by the variable height of the columns, which is dependent on the mass accretion rate. This study delves into the process of pulse profile formation in bright X-ray pulsars, taking into account stochastic fluctuations in the mass accretion rate, the corresponding variations in accretion column geometry, and gravitational bending. Our analysis reveals that potential eclipses of accretion columns by a neutron star during their spin period should manifest specific features in pulse profile variability. Applying a novel pulse profile analysis technique, we successfully detect these features in the bright X-ray transient V 0332+53 at luminosities $\gtrsim 2\times 10^{38}\, {\rm erg\ \rm s^{-1}}$. This detection serves as compelling evidence for the eclipse of an accretion column by a neutron star. Detection of the eclipse places constraints on the relation between neutron star mass, radius, and accretion column height. Specifically, we can establish an upper limit on the accretion column height, which is crucial for refining theoretical models of extreme accretion.

Pulsed Iron Line Emission from the First Galactic Ultraluminous X-Ray Pulsar Swift J0243.6+6124

We report the phase-resolved spectral results of the first Galactic pulsating ultraluminous X-ray source (PULX) Swift J0243.6+6124, modeling its 2017–2018 outburst peak using data collected by the Hard X-ray Modulation Telescope (Insight-HXMT). The broad energy coverage of Insight-HXMT allows us to obtain a more accurate spectral continuum to reduce the coupling of broad iron line profiles with other components. We use three different continuum spectrum models but obtain similar iron line results. For the first time, we detect the pulse characteristics of the broad iron line in a PULX. The variation in the width and intensity of this iron line with σ ∼ 1.2–1.5 keV has a phase offset of about 0.25 from the pulse phase. We suggest that the uneven irradiation of the thick inner disk by the accretion column produces the modulated variation of the broad iron line. In addition, the nonpulsed narrow line is suggested to come from the outer disk region.

A survey for variable young stars with small telescopes – IX. Evolution of spot properties on YSOs in IC 5070

ABSTRACT We present spot properties on 32 periodic young stellar objects in IC 5070. Long term, ∼5 yr, light curves in the V, R, and I-bands are obtained through the HOYS (Hunting Outbursting Young Stars) citizen science project. These are dissected into 6 months long slices, with 3 months oversampling, to measure 234 sets of amplitudes in all filters. We fit 180 of these with reliable spot solutions. Two thirds of spot solutions are cold spots, the lowest is 2150 K below the stellar temperature. One third are warm spots that are above the stellar temperature by less than ∼2000 K. Cold and warm spots have maximum surface coverage values of 40 per cent, although only 16 per cent of warm spots are above 20 per cent surface coverage as opposed to 60 per cent of the cold spots. Warm spots are most likely caused by a combination of plages and low-density accretion columns, most common on objects without inner disc excess emission in K − W2. Five small hot spot solutions have <3 per cent coverage and are 3000–5000 K above the stellar temperature. These are attributed to accretion, and four of them occur on the same object. The majority of our objects are likely to be accreting. However, we observe very few accretion hot spots as either the accretion is not stable on our time-scale or the photometry is dominated by other features. We do not identify cyclical spot behaviour on the targets. We additionally identify and discuss a number of objects that have interesting amplitudes, phase changes, or spot properties.

The enigmatic dance of the HD 189733A system: A quest for accretion

Context. Several studies suggest that the emission properties of a star can be significantly affected by its interaction with a nearby planet through magnetic fields or interaction between their respective winds. However, the actual observability of these effects remains a subject of debate. An illustrative example is the HD 189733A system: certain characteristics of its emissions have been interpreted as indicative of ongoing interactions between the star and its associated planet. Other studies attribute these characteristics to the coronal activity of the star. Aims. In this study we aimed to investigate whether the observed stellar X-ray flare events, which appear to be in phase with the planetary period in the HD 189733A system, could be attributed to the accretion of the planetary wind onto the stellar surface or if they resulted from an interaction between the planetary and stellar winds. Methods. We developed a 3D magnetohydrodynamic model with the PLUTO code that describes the system HD 189733A , including the central host star and its hot Jupiter along with their respective winds. The effects of gravity and the magnetic fields of both the star and the planet are taken into account. Results. Our analysis reveals that, in the cases examined in this study, the accretion scenario is only viable when the stellar magnetic field strength is at 5 G and the planetary magnetic field strength is at 1 G. In this scenario, the Rayleigh-Taylor instabilities lead to the formation of an accretion column that connects the star to the planet. Once formed the accretion column remains stable for the duration of the simulation. The accretion column produces an accretion rate of the order of 1012 g s−1 and shows an average density of about 107 cm−3. In the other case explored, the accretion column does not form because the Rayleigh-Taylor instability is suppressed by the stronger magnetic field intensities assumed for both the star and the planet. We synthesized the emission resulting from the shocked planetary wind and found that the total X-ray emission ranges from 5 × 1023 to 1024 erg s−1. Conclusions. In the case of accretion, the emission originating from the hotspot cannot be distinguished from the coronal activity. Also, the interaction between the planetary and stellar winds cannot be responsible for the X-ray emission, as the total emission produced is about four orders of magnitude lower than the average X-ray luminosity of the star.

Confirmation of the presence of a CRSF in the NICER spectrum of X 1822-371

Aims. X 1822-371 is an eclipsing binary system with a period close to 5.57 h and an orbital period derivative Ṗorb of 1.42(3) × 10−10 s s−1. The extremely high value of its Ṗorb is compatible with a super-Eddington mass transfer rate from the companion star and, consequently, an intrinsic luminosity at the Eddington limit. The source is also an X-ray pulsar, it shows a spin frequency of 1.69 Hz and is in a spin-up phase with a spin frequency derivative of 7.4 × 10−12 Hz s−1. Assuming a luminosity at the Eddington limit, a neutron star magnetic field strength of B = 8 × 1010 G is estimated. However, a direct measure of B could be obtained observing a CRSF in the energy spectrum. Analysis of XMM-Newton data suggested the presence of a cyclotron line at 0.73 keV, with an estimated magnetic field strength of B = (8.8 ± 0.3)×1010 G. Methods. Here we analyze the 0.3–50 keV broadband spectrum of X 1822-371 combining a 0.3–10 keV NICER spectrum and a 4.5–50 keV NuSTAR spectrum to investigate the presence of a cyclotron absorption line and the complex continuum emission spectrum. Results. The NICER spectrum confirms the presence of a cyclotron line at 0.66 keV. The continuum emission is modeled with a Comptonized component, a thermal component associated with the presence of an accretion disk truncated at the magnetospheric radius of 105 km and a reflection component from the disk blurred by relativistic effects. Conclusions. We confirm the presence of a cyclotron line at 0.66 keV inferring a NS magnetic field of B = (7.9 ± 0.5)×1010 G and suggesting that the Comptonized component originates in the accretion columns.

Broad-band noise and quasi-periodic oscillation characteristics of the X-ray pulsar RX J0440.9+4431

ABSTRACT We present a comprehensive timing analysis on the Be/X-ray binary pulsar RX J0440.9+4431 using observations from NICER and Insight-HXMT during the 2022–2023 outburst. The power density spectrum (PDS) of RX J0440.9+4431 exhibits typical aperiodic variability in X-ray flux across a wide frequency range. During a supercritical accretion state, we detect quasi-periodic oscillations (QPOs) at 0.2–0.5 Hz in the light curves of five pulses for RX J0440.9+4431. The observed QPOs manifest during flares, while the flares appear at the peaks of the pulse profiles on a time-scale of seconds and are primarily caused by an increase in hard photons. These flares can be explained by increased material ingestion in the accretion column at a fixed phase, primarily generating hard photons. Alternatively, an increase in accretion rate, independent of phase, may result in highly beamed hard photons within the accretion column, causing the flares. We argue the origin of QPOs to instabilities within the accretion flow. Additionally, we find that the break frequencies in the noise power spectra align well with $\propto L_{\mathrm{x}}^{3 / 7}$ across three orders of magnitude in the luminosity, which points to a relatively strong magnetic field in RX J0440.9+4431, estimated to be ~1013 G.

Insights into the phase-dependent cyclotron line feature in XTE J1946+274: an AstroSat and Insight-HXMT view

ABSTRACT XTE J1946+274 is a Be/X-ray binary with a 15.8 s spin period and 172 d orbital period. Using RXTE/PCA data of the 1998 outburst, a cyclotron line around 37 keV was reported. The presence of this line, its dependence on the pulse phase, and its variation with luminosity have been of some debate since. In this work, we present the re-analysis of two AstroSat observations: one made during the rising phase of the 2018 outburst and the other during the declining phase of the 2021 outburst. We also present a new analysis of the Insight-HXMT observations of the source at the peak of the 2018 outburst. We find the source to be spinning up over the course of the outburst and spinning down between the two outbursts. We report the presence of a higher cyclotron line energy using the 2018 AstroSat observation (∼45 keV) and 2018 Insight-HXMT observation (∼50 keV) and a line at ∼40 keV during the declining phase of the 2021 outburst using data from AstroSat. We also investigate the pulse phase dependence of the cyclotron line parameters and find that the line is significantly detected in all the phases of both AstroSat observations, along with showing variation with the pulse phase. This differs from the previous results reported using BeppoSAX and NuSTAR. We explain this behaviour of the cyclotron line to be due to photon spawning and different accretion column radii at the two poles of this neutron star.

A global 3D simulation of magnetospheric accretion – I. Magnetically disrupted discs and surface accretion

ABSTRACT We present a 3D ideal MHD simulation of magnetospheric accretion on to a non-rotating star. The accretion process unfolds with intricate 3D structures driven by various mechanisms. First, the disc develops filaments at the magnetospheric truncation radius (RT) due to magnetic interchange instability. These filaments penetrate deep into the magnetosphere, form multiple accretion columns, and eventually impact the star at ∼30o from the poles at nearly the free-fall speed. Over 50 per cent (90 per cent) of accretion occurs on just 5 per cent (20 per cent) of the stellar surface. Secondly, the disc region outside RT develops large-scale magnetically dominated bubbles, again due to magnetic interchange instability. These bubbles orbit at a sub-Keplerian speed, persisting for a few orbits while leading to asymmetric mass ejection. The disc outflow is overall weak because of mostly closed field lines. Thirdly, magnetically supported surface accretion regions appear above the disc, resembling a magnetized disc threaded by net vertical fields, a departure from traditional magnetospheric accretion models. Stellar fields are efficiently transported into the disc region due to above instabilities, contrasting with the ‘X-wind’ model. The accretion rate on to the star remains relatively steady with a 23 per cent standard deviation. The periodogram reveals variability occurring at around 0.2 times the Keplerian frequency at RT, linked to the large-scale magnetic bubbles. The ratio of the spin-up torque to $\dot{M}(GM_*R_T)^{1/2}$ is around 0.8. Finally, after scaling the simulation, we investigate planet migration in the inner protoplanetary disc. The disc driven migration is slow in the MHD turbulent disc beyond RT, while aerodynamic drag plays a significant role in migration within RT.

Trends in torques acting on the star during a star-disk magnetospheric interaction

Aims. We assess the modification of angular momentum transport in various configurations of star-disk accreting systems based on numerical simulations with different parameters. In particular, we quantify the torques exerted on a star by the various components of the flow and field in our simulations of a star-disk magnetospheric interaction. Methods. In a suite of resistive and viscous numerical simulations, we obtained results using different stellar rotation rates, dipole magnetic field strengths, and resistivities. We probed a part of the parameter space with slowly rotating central objects, up to 20% of the Keplerian rotation rate at the equator. Different components of the flow in star-disk magnetospheric interaction were considered in the study: a magnetospheric wind (i.e., the “stellar wind”) ejected outwards from the stellar vicinity, matter infalling onto the star through the accretion column, and a magnetospheric ejection launched from the magnetosphere. We also took account of trends in the total torque in the system and in each component individually. Results. We find that for all the stellar magnetic field strengths, B⋆, the anchoring radius of the stellar magnetic field in the disk is extended with increasing disk resistivity. The torque exerted on the star is independent of the stellar rotation rate, Ω⋆, in all the cases without magnetospheric ejections. In cases where such ejections are present, there is a weak dependence of the anchoring radius on the stellar rotation rate, with both the total torque in the system and torque on the star from the ejection and infall from the disk onto the star proportional to Ω⋆B3. The torque from a magnetospheric ejection is proportional to Ω⋆4. Without the magnetospheric ejection, the spin-up of the star switches to spin-down in cases involving a larger stellar field and faster stellar rotation. The critical value for this switch is about 10% of the Keplerian rotation rate.

Apparent luminosity and pulsed fraction affected by gravitational lensing of accretion columns in bright X-ray pulsars

ABSTRACT The luminosity of X-ray pulsars is their key parameter determining the geometry and physical conditions of the accretion flow both on the spatial scales of a binary system and on much smaller scales of emitting regions located close to the stellar surface. Traditionally, the luminosity of X-ray pulsars is estimated out of the X-ray energy flux averaged over the pulsed period and the estimated distance to the source. Due to the anisotropy of X-ray emission, the luminosity estimated on the base of the observed pulse profile can differ from the actual one. Super-cw2ritical X-ray pulsars with accretion columns are of particular interest because the X-ray flux from columns is a matter of strong gravitational lensing by a neutron star. Using toy model of an accretion column, we simulate beam patterns in super-critical X-ray pulsars, construct theoretical pulse profiles for different geometries and mutual orientations of pulsars and distant observers and show that despite strong light bending, the typical deviation of the apparent luminosity from the actual one is $\sim 20~{{\ \rm per \, cent}}$ only, and in $\sim 90~{{\ \rm per \, cent}}$ of cases, the apparent luminosity 0.8L ≲ Lapp ≲ 1.25L. However, the shape of the pulse profiles is strongly affected by the geometry of the emitting region. We show that the appearance and growth of accretion columns tend to be accompanied by an increase of observed pulsed fraction, which is in agreement with the recent observations of bright X-ray transients.