We report on the development of MPI-AMRVAC version 2.0, which is an open-source framework for parallel, grid-adaptive simulations of hydrodynamic and magnetohydrodynamic (MHD) astrophysical ...applications. The framework now supports radial grid stretching in combination with adaptive mesh refinement (AMR). The advantages of this combined approach are demonstrated with one-dimensional, two-dimensional, and three-dimensional examples of spherically symmetric Bondi accretion, steady planar Bondi-Hoyle-Lyttleton flows, and wind accretion in supergiant X-ray binaries. Another improvement is support for the generic splitting of any background magnetic field. We present several tests relevant for solar physics applications to demonstrate the advantages of field splitting on accuracy and robustness in extremely low-plasma β environments: a static magnetic flux rope, a magnetic null-point, and magnetic reconnection in a current sheet with either uniform or anomalous resistivity. Our implementation for treating anisotropic thermal conduction in multi-dimensional MHD applications is also described, which generalizes the original slope-limited symmetric scheme from two to three dimensions. We perform ring diffusion tests that demonstrate its accuracy and robustness, and show that it prevents the unphysical thermal flux present in traditional schemes. The improved parallel scaling of the code is demonstrated with three-dimensional AMR simulations of solar coronal rain, which show satisfactory strong scaling up to 2000 cores. Other framework improvements are also reported: the modernization and reorganization into a library, the handling of automatic regression tests, the use of inline/online Doxygen documentation, and a new future-proof data format for input/output.
Context.
Low-luminosity supermassive and stellar-mass black holes (BHs) may be embedded in a collisionless and highly magnetized plasma. They show nonthermal flares indicative of particles being ...accelerated up to relativistic speeds by dissipative processes in the vicinity of the BH. During near-infrared flares from the supermassive BH Sagittarius A
*
(Sgr A
*
), the GRAVITY Collaboration detected circular motion and polarization evolution, which suggest the presence of transient synchrotron-emitting hot spots moving around the BH.
Aims.
We study 3D reconnecting current layers in the magnetosphere of spinning BHs to determine whether plasma-loaded flux ropes which are formed near the event horizon could reproduce the hot spot observations and help constrain the BH spin.
Methods.
We performed global 3D particle-in-cell simulations in Kerr spacetime of a pair plasma embedded in a strong and large-scale magnetic field originating in a perfectly conducting disk in prograde Keplerian rotation.
Results.
A cone-shaped current layer develops which surrounds the twisted open magnetic field lines threading the event horizon. Spinning magnetic field lines coupling the disk to the BH inflate and reconnect a few gravitational radii above the disk. This quasi-periodic cycle accelerates particles, which accumulate in a few macroscopic flux ropes rotating with the outermost coupling magnetic field line. Once flux ropes detach, they propagate in the current layer following what appears as a rapidly opening spiral when seen face-on. A single flux rope carries enough relativistic electrons and positrons to emit synchrotron radiation at levels suitable to reproduce the peak-luminosity of the flares of Sgr A
*
but it quickly fades away as it flows away.
Conclusions.
Our kinematic analysis of the flux ropes motion favors a BH spin of 0.65 to 0.8 for Sgr A
*
. The duration of the flares of Sgr A
*
can only be explained provided the underlying magnetic loop seeded in the disk mid-plane has a finite lifetime and azimuthal extension. In this scenario, the hot spot corresponds to a spinning arc along which multiple reconnection sites power the net emission as flux ropes episodically detach.
Abstract
Supergiant X-ray binaries (sgxb) host a compact object, often a neutron star (NS), orbiting an evolved O/B star. Mass transfer proceeds through the intense line-driven wind of the stellar ...donor, a fraction of which is captured by the gravitational field of the NS. The subsequent accretion process on to the NS is responsible for the abundant X-ray emission from sgxb. They also display peak-to-peak variability of the X-ray flux by a factor of a few 10–100, along with changes in the hardness ratios possibly due to varying absorption along the line of sight. We use recent radiation-hydrodynamic simulations of inhomogeneities (a.k.a. clumps) in the non-stationary wind of massive hot stars to evaluate their impact on the time-variable accretion process. For this, we run 3D hydrodynamic simulations of the wind in the vicinity of the accretor to investigate the formation of the bow shock and follow the inhomogeneous flow over several spatial orders of magnitude, down to the NS magnetosphere. In particular, we show that the impact of the wind clumps on the time variability of the intrinsic mass accretion rate is severely tempered by the crossing of the shock, compared to the purely ballistic Bondi–Hoyle–Lyttleton estimation. We also account for the variable absorption due to clumps passing by the line of sight and estimate the final effective variability of the column density and mass accretion rate for different orbital separations. Finally, we compare our results to the most recent analysis of the X-ray flux and the hardness ratio in Vela X-1.
Complex stellar winds from evolved stars
Stars less than eight times the mass of the Sun end their lives as planetary nebulae, structures of ionized gas thrown off by the star and heated by the ...exposed stellar core. Planetary nebulae are often bipolar in shape or contain complex morphological features such as rings or spirals. Decin
et al.
observed the stellar winds of 14 stars during their asymptotic giant branch (AGB) phase of stellar evolution, which immediately precedes the planetary nebula phase. They found morphologies in the AGB winds similar to planetary nebulae and demonstrated that they are produced by the influence of a binary companion on the AGB wind.
Science
, this issue p.
1497
Complex morphologies in the winds of asymptotic giant branch stars and planetary nebulae are produced by binary interactions.
Binary interactions dominate the evolution of massive stars, but their role is less clear for low- and intermediate-mass stars. The evolution of a spherical wind from an asymptotic giant branch (AGB) star into a nonspherical planetary nebula (PN) could be due to binary interactions. We observed a sample of AGB stars with the Atacama Large Millimeter/submillimeter Array (ALMA) and found that their winds exhibit distinct nonspherical geometries with morphological similarities to planetary nebulae (PNe). We infer that the same physics shapes both AGB winds and PNe; additionally, the morphology and AGB mass-loss rate are correlated. These characteristics can be explained by binary interaction. We propose an evolutionary scenario for AGB morphologies that is consistent with observed phenomena in AGB stars and PNe.
Bondi–Hoyle accretion configurations occur as soon as a gravitating body is immersed in an ambient medium with a supersonic relative velocity. From wind-accreting X-ray binaries to runaway neutron ...stars, such a regime has been witnessed many times and is believed to account for shock formation, the properties of which can be only marginally derived analytically. In this paper, we present the first results of the numerical characterization of the stationary flow structure of Bondi–Hoyle accretion on to a compact object, from the large-scale accretion radius down to the vicinity of the compact body. For different Mach numbers, we study the associated bow shock. It turns out that those simulations confirm the analytical prediction by Foglizzo & Ruffert concerning the topology of the inner sonic surface with an adiabatic index of 5/3. They also enable us to derive the related mass accretion rates, the position and the temperature of the bow shock, as function of the flow parameters, along with the transverse density and temperature profiles in the wake.
Context.
In the Milky Way, ∼18 Wolf-Rayet+O star (WR+O) binaries are known with estimates of their stellar and orbital parameters. Whereas black hole+O star (BH+O) binaries are thought to evolve from ...WR+O binaries, only one such system is known in the Milky Way. To resolve this disparity, it was suggested recently that upon core collapse, the WR stars receive large kicks such that most of the binaries are disrupted.
Aims.
We reassess this issue, with a particular emphasis on the uncertainty in predicting the X-ray emission from wind-accreting BHs in BH+O binaries, which is key to identifying such systems.
Methods.
BH+O systems are thought to be X-ray bright only when an accretion disk forms around the BHs. We followed the methodology of previous work and applied an improved analytic criterion for the formation of an accretion disk around wind accreting BHs. We then used stellar evolutionary models to predict the properties of the BH+O binaries which are expected to descend from the observed WR+O binaries if the WR stars would form BHs without a natal kick.
Results.
We find that disk formation sensitively depends on the O stars’ wind velocity, the amount of specific angular momentum carried by the wind, the efficiency of angular momentum accretion by the BH, and the spin of the BH. We show that whereas the assumption of a low wind velocity may lead to the prediction that most of the BH+O star binaries will have an extended X-ray bright period, this is not the case when typical wind velocities of O stars are considered. We find that a high spin of the BH can boost the duration of the X-ray active phase as well as the X-ray brightness during this phase. This produces a strong bias for detecting high mass BH binaries in X-rays with high BH spin parameters.
Conclusions.
We find that large BH formation kicks are not required to understand the sparsity of X-ray bright BH+O stars in the Milky Way. Probing for a population of X-ray silent BH+O systems with alternative methods can likely inform us about BH kicks and the necessary conditions for high energy emission from high mass BH binaries.
Context.
Sagittarius A*, the supermassive black hole at the center of our Galaxy, exhibits episodic near-infrared flares. The recent monitoring of three such events with the GRAVITY instrument has ...shown that some flares are associated with orbital motions in the close environment of the black hole. The GRAVITY data analysis indicates a super-Keplerian azimuthal velocity, while (sub-) Keplerian velocity is expected for the hot flow surrounding the black hole.
Aims.
We develop a semi-analytic model of the Sagittarius A* flares based on an ejected large plasmoid, inspired by recent particle-in-cell global simulations of black hole magnetospheres. We model the infrared astrometric and photometric signatures associated with this model.
Methods.
We considered a spherical macroscopic hot plasma region that we call a large plasmoid. This structure was ejected along a conical orbit in the vicinity of the black hole. This plasmoid was assumed to be formed by successive mergers of smaller plasmoids produced through magnetic reconnection that we did not model. Nonthermal electrons were injected into the plasmoid. We computed the evolution of the electron-distribution function under the influence of synchrotron cooling. We solved the radiative transfer problem associated with this scenario and transported the radiation along null geodesics of the Schwarzschild space time. We also took the quiescent radiation of the accretion flow into account, on top of which the flare evolves.
Results.
For the first time, we successfully account for the astrometric and flux variations of the GRAVITY data with a flare model that incorporates an explicit modeling of the emission mechanism. The prediction of our model and recent data agree well. In particular, the azimuthal velocity of the plasmoid is set by the magnetic field line to which it belongs, which is anchored in the inner parts of the accretion flow, hence the super-Keplerian motion. The astrometric track is also shifted with respect to the center of mass due to the quiescent radiation, in agreement with the difference measured with the GRAVITY data.
Conclusions.
These results support the hypothesis that magnetic reconnection in a black hole magnetosphere is a viable model for the infrared flares of Sagittarius A*.
Context. In high mass X-ray binaries, an accreting compact object orbits a high mass star, which loses mass through a dense and in homogeneous wind.
Aims. Using the compact object as an X-ray ...backlight, the time variability of the absorbing column density in the wind can be exploited in order to shed light on the micro-structure of the wind and obtain unbiased stellar mass-loss rates for high mass stars.
Methods. We developed a simplified representation of the stellar wind where all the matter is gathered in spherical “clumps” that are radially advected away from the star. This model enables us to explore the connections between the stochastic properties of the wind and the variability of the column density for a comprehensive set of parameters related to the orbit and to the wind micro-structure, such as the size of the clumps and their individual mass. In particular, we focus on the evolution with the orbital phase of the standard deviation of the column density and of the characteristic duration of enhanced absorption episodes. Using the porosity length, we derive analytical predictions and compare them to the standard deviations and coherence time scales that were obtained.
Results. We identified the favorable systems and orbital phases to determine the wind micro-structure. The coherence time scale of the column density is shown to be the self-crossing time of a single clump in front of the compact object. We thus provide a procedure to get accurate measurements of the size and of the mass of the clumps, purely based on the observable time variability of the column density.
Conclusions. The coherence time scale grants direct access to the size of the clumps, while their mass can be deduced separately from the amplitude of the variability. We further show how monitoring the variability at superior conjunctions can probe the onset of the clump-forming region above the stellar photosphere. If the high column density variations in some high mass X-ray binaries are due to unaccreted clumps which are passing by the line-of-sight, this would require high mass clumps to reproduce the observed peak-to-peak amplitude and coherence time scales. These clump properties are marginally compatible with the ones derived from radiative-hydrodynamics simulations. Alternatively, the following components could contribute to the variability of the column density: larger orbital scale structures produced by a mechanism that has yet to be identified or a dense environment in the immediate vicinity of the accretor, such as an accretion disk, an outflow, or a spherical shell surrounding the magnetosphere of the accreting neutron star.
We present the results of a search through the photometric database of Kepler eclipsing binaries looking for evidence of hierarchical triple-star systems. The presence of a third star orbiting the ...binary can be inferred from eclipse timing variations. We apply a simple algorithm in an automated determination of the eclipse times for all 2157 binaries. The "calculated" eclipse times, based on a constant period model, are subtracted from those observed. The resulting O - C (observed minus calculated times) curves are then visually inspected for periodicities in order to find triple-star candidates. After eliminating false positives due to the beat frequency between the ~1/2 hr Kepler cadence and the binary period, 39 candidate triple systems were identified. The periodic O - C curves for these candidates were then fit for contributions from both the classical Roemer delay and so-called physical delay, in an attempt to extract a number of the system parameters of the triple. We discuss the limitations of the information that can be inferred from these O - C curves without further supplemental input, e.g., ground-based spectroscopy. Based on the limited range of orbital periods for the triple-star systems to which this search is sensitive, we can extrapolate to estimate that at least 20% of all close binaries have tertiary companions.
Bright and eclipsing, the high-mass X-ray binary Vela X-1 offers a unique opportunity to study accretion onto a neutron star from clumpy winds of O/B stars and to disentangle the complex accretion ...geometry of these systems. In Chandra-HETGS spectroscopy at orbital phase ~0.25, when our line of sight towards the source does not pass through the large-scale accretion structure such as the accretion wake, we observe changes in overall spectral shape on timescales of a few kiloseconds. This spectral variability is, at least in part, caused by changes in overall absorption and we show that such strongly variable absorption cannot be caused by unperturbed clumpy winds of O/B stars. We detect line features from high and low ionization species of silicon, magnesium, and neon whose strengths and presence depend on the overall level of absorption. These features imply a co-existence of cool and hot gas phases in the system, which we interpret as a highly variable, structured accretion flow close to the compact object such as has been recently seen in simulations of wind accretion in high-mass X-ray binaries.