Hercules X-1 is one of the best studied highly magnetised neutron star X-ray binaries with a wealth of archival data. We present the discovery of an ionised wind in its X-ray spectrum when the source ...is in the high state. The wind detection is statistically significant in most of the XMM-Newton observations, with velocities ranging from 200 to 1000 km/s. Observed features in the iron K band can be explained by both wind absorption or by a forest of iron emission lines. However, we also detect nitrogen, oxygen and neon absorption lines at the same systematic velocity in the high-resolution RGS grating spectra. The wind must be launched from the accretion disc, and could be the progenitor of the UV absorption features observed at comparable velocities, but the latter likely originate at significantly larger distances from the compact object. We find strong correlations between the ionisation level of the outflowing material and the ionising luminosity as well as the super-orbital phase. If the luminosity is driving the correlation, the wind could be launched by a combination of Compton heating and radiation pressure. If instead the super-orbital phase is the driver for the variations, the observations are likely scanning the wind at different heights above the warped accretion disc. If this is the case, we can estimate the wind mass outflow rate, corrected for the limited launching solid angle, to be roughly 70% of the mass accretion rate.
The Galactic black hole X-ray binary MAXI J1820+070 had a bright outburst in 2018 when it became the second brightest X-ray source in the Sky. It was too bright for X-ray CCD instruments such as ...XMM-Newton and Chandra, but was well observed by photon-counting instruments such as NICER and NuSTAR. We report here on the discovery of an excess emission component during the soft state. It is best modelled with a blackbody spectrum in addition to the regular disk emission, modelled either as diskbb or kerrbb. Its temperature varies from about 0.9 to 1.1 keV which is about 30 to 80 per cent higher than the inner disc temperature of diskbb. Its flux varies between 4 and 12 percent of the disc flux. Simulations of magnetised accretion discs have predicted the possibility of excess emission associated with a non-zero torque at the Innermost Stable Circular Orbit (ISCO) about the black hole, which from other NuSTAR studies lies at about 5 gravitational radii or about 60 km (for a black hole mass is 8 M). In this case the emitting region at the ISCO has a width varying between 1.3 and 4.6 km and would encompass the start of the plunge region where matter begins to fall freely into the black hole.
We describe a Godunov-type magnetohydrodynamic (MHD) code based on the Miyoshi and Kusano (2005) solver which can be used to solve various astrophysical hydrodynamic and MHD problems. The energy ...equation is in the form of entropy conservation. The code has been implemented on several different coordinate systems: 2.5D axisymmetric cylindrical coordinates, 2D Cartesian coordinates, 2D plane polar coordinates, and fully 3D cylindrical coordinates. Viscosity and diffusivity are implemented in the code to control the accretion rate in the disk and the rate of penetration of the disk matter through the magnetic field lines. The code has been utilized for the numerical investigations of a number of different astrophysical problems, several examples of which are shown.
Recent radio emission, polarization, and Faraday rotation maps of the radio jet of the galaxy 3C 303 have shown that one knot of this jet has a {\it galactic}-scale electric current of \(\sim 3\times ...10^{18}\) Ampère flowing along the jet axis (Kronberg et al. 2011). We develop the theory of relativistic Poynting-flux jets which are modeled as a transmission line carrying a DC current \(I_0\), having a potential drop \(V_0\), and a definite impedance \({\cal Z}_0 =90(u_z/c)\Omega\), where \(u_z\) is the bulk velocity of the jet plasma. The electromagnetic energy flow in the jet is \({\cal Z}_0 I_0^2\). The observed current in 3C 303 can be used to calculate the electromagnetic energy flow in this magnetically dominated jet. Time-dependent but not necessarily small perturbations of a Poynting-flux jet - possibly triggered by a gas cloud penetrating the jet - are described by "telegrapher's equations," which predict the propagation speed of disturbances and the effective wave impedance \({\cal Z}\). The disturbance of a Poynting jet by the cloud gives rise to localized dissipation in the jet which may explain the enhanced synchrotron radiation in the knots of the 3C 303 jet, and in the apparently stationary knot HST-1 in the jet from the nucleus of the galaxy M87 (Biretta et al. 1999).
We have carried out and analyzed a set of axisymmetric MHD simulations of the evolution of a turbulent/diffusive accretion disc around an initially unmagnetized star. The disc is initially threaded ...by a weak magnetic field where the magnetic pressure is significantly less than the kinetic pressure in the disc. The viscosity and magnetic diffusivity are modelled by two "alpha" parameters, while the coronal region above the disc is treated using ideal MHD. The initial magnetic field is taken to consist of three poloidal field loops threading the disc. The motivation for this study is to understand the advection of disc matter and magnetic field by the turbulent/diffusive disc. At early times the innermost field loop twists and its field lines become open. The twisting of the opened field lines leads to the formation of both an inner collimated, magnetically-dominated jet, and at larger distances from the axis a matter dominated uncollimated wind. For later times, the strength of the magnetic field decreases owing to field reconnection and annihilation in the disc. For the early times, we have derived from the simulations both the matter accretion speed in the disc and the accretion speed of the magnetic field. We show that the derived matter accretion speed agrees approximately with the predictions of a model where the accretion speed is the sum of two terms, one due to the disc's viscosity (which gives a radial outflow of angular momentum in the disc), and a second due to the twisted magnetic field at the disc's surface (which gives a vertical outflow of angular momentum). For early times we find that the magnetic contribution is roughly twice the viscous contribution for the case where the alpha parameters are both equal to 0.1. At later times the magnetic contribution to the matter speed becomes small compared to the viscous contribution.
We study function spaces and extension results in relation with Dirichlet problems involving integrodifferential operators. For such problems, data are prescribed on the complement of a given domain ...Ω⊂Rd. We introduce a function space that serves as a trace space for nonlocal Dirichlet problems and study related extension results.
ABSTRACT
We study non-axisymmetric features of 3D line-driven winds in the Sobolev approximation, where the optical depth is calculated using the full velocity gradient. We find that non-axisymmetric ...density features, so-called clumps, form primarily at the base of the wind on super-Sobolev length-scales. The density of clumps differs by a factor of ∼3 from the azimuthal average, the magnitude of their velocity dispersion is comparable to the flow velocity, and they produce ∼20 per cent variations in the column density. Clumps may be observable because differences in density produce enhancements in emission and absorption profiles or through their velocity dispersion that enhances line broadening.
Abstract
We study mass outflows driven from accretion discs by radiation pressure due to spectral lines. To investigate non-axisymmetric effects, we use the athena++ code and develop a new module to ...account for radiation pressure driving. In 2D, our new simulations are consistent with previous 2D axisymmetric solutions by Proga et al., who used the zeus 2D code. Specifically, we find that the disc winds are time dependent, characterized by a dense stream confined to ∼45° relative to the disc mid-plane and bounded on the polar side by a less dense, fast stream. In 3D, we introduce a vertical, ϕ-dependent, subsonic velocity perturbation in the disc mid-plane. The perturbation does not change the overall character of the solution but global outflow properties such as the mass, momentum, and kinetic energy fluxes are altered by up to 100 per cent. Non-axisymmetric density structures develop and persist mainly at the base of the wind. They are relatively small, and their densities can be a few times higher than the azimuthal average. The structure of the non-axisymmetric and axisymmetric solutions differ also in other ways. Perhaps most importantly from the observational point of view are the differences in the so-called clumping factors, that serve as a proxy for emissivity due to two body processes. In particular, the spatially averaged clumping factor over the entire fast stream, while it is of a comparable value in both solutions, it varies about 10 times faster in the non-axisymmetric case.