BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Guillaume Dubus

Guillaume Dubus appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-ph.GAastro-ph.SRastro-ph.COastro-phastro-ph.IMphysics.plasm-ph

Catalog footprint

What is connected

29works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2026arXiv

Magnetic winds in resistive compact binary discs

Dwarf Novae and low-mass X-ray binaries are eruptive binary systems comprised of a Roche-lobe overflowing solar-type star and an accreting compact object. Their recurrence time can be explained by a low-accreting phase, the quiescence, during which the angular momentum transport parameter is inferred to be $α\approx 0.01$ by the Disc Instability Model. Non-magnetics mechanisms, such as spiral wave transport, only achieve angular momentum transport an order of magnitude too low, at best, because these discs are so thin in quiescence. During this phase, the Magneto-rotational Instability is known to be suppressed by the increased resistivity of the weakly ionised plasma. Studying these thin magnetised discs is a numerical challenge because of the wide range of scales to be resolved. Thanks to the new GPU-accelerated code Idefix, we produce global 3D MHD simulations of a very thin disc $(H/R = 0.01)$ for the first time. We explore the possibility that an MHD wind arises and increases accretion in low magnetic Reynolds number $(\mathrm{Rm}\approx100)$ and realistic plasma parameter ($β\approx 1000$) regimes. We observe that the MRI is only quenched in the resistive disc bulk but survives in the disc atmosphere. This drives strong accretion and wind launching. We quantify the efficiency of the resulting wind and measure its global effect on the disc. We explore the effect of the initial disc magnetisation and compare the accretion/ejection regime with and without resistivity.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Intra-pulse variability induced by plasmoid formation in pulsar magnetospheres

Pulsars show irregularities in their pulsed radio emission that originate from propagation effects and the intrinsic activity of the source. In this work, we investigate the role played by magnetic reconnection and the formation of plasmoids in the pulsar wind current sheet as a possible source of intrinsic pulse-to-pulse variability in the incoherent, high-energy emission pattern. We used a two-dimensional particle-in-cell simulation of an orthogonal pulsar magnetosphere restricted to the plane perpendicular to the star spin axis. We evolved the solution for several tens of pulsar periods to gather a statistically significant sample of synthetic pulse profiles. The formation of plasmoids leads to strong pulse-to-pulse variability in the form of multiple short, bright subpulses, which appear only on the leading edge of each main pulse. These secondary peaks of emission are dominated by the dozen plasmoids that can grow up to macroscopic scales. They emerge from the high end of the hierarchical merging process occurring along the wind current layer. The flux of the subpulses is correlated with their width in phase. Although the full-scale separation is not realistic, we argue that the simulation correctly captures the demographics and the properties of the largest plasmoids, and therefore of the brightest subpulses. The prediction of subpulses at specific pulse phases provides a new observational test of the magnetic reconnection scenario as the origin of the pulsed incoherent emission. High-time-resolution observations of the Crab pulsar in the optical range may be the most promising source to target for this purpose.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetic field transport in compact binaries

Dwarf novae (DNe) and low mass X-ray binaries (LMXBs) show eruptions that are thought to be due to a thermal-viscous instability in their accretion disk. These eruptions provide constraints on angular momentum transport mechanisms. We explore the idea that angular momentum transport could be controlled by the dynamical evolution of the large scale magnetic field. We study the impact of different prescriptions for the magnetic field evolution on the dynamics of the disk. This is a first step in confronting the theory of magnetic field transport with observations. We develop a version of the disk instability model that evolves the density, the temperature and the large scale vertical magnetic flux together. We take into account the accretion driven by turbulence or by a magnetized outflow. To evolve the magnetic flux, we use a toy model with physically motivated prescriptions depending mainly on the local magnetization. We find that allowing magnetic flux to be advected inwards provides the best agreement with DNe lightcurves. This leads to a hybrid configuration with an inner magnetized disk, driven by angular momentum losses to an MHD outflow, sharply transiting to an outer weakly-magnetized turbulent disk, where the eruptions are triggered. The dynamical impact is equivalent to truncating a viscous disk so that it does not extend down to the compact object, with the truncation radius dependent on the magnetic flux and evolving as $\dot{M}^{-2/3}$. Models of DNe and LMXBs lightcurves typically require the outer, viscous disk to be truncated in order to match observations. There is no generic explanation for this truncation. We propose that it is a natural outcome of the presence of large-scale magnetic fields in both DNe and LMXBs, the magnetic flux accumulating towards the center to produce a magnetized disk with a fast accretion timescale.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multi-dimensional simulations of ergospheric pair discharges around black holes

Black holes are known to launch powerful relativistic jets and emit highly variable gamma radiation. How these jets are loaded with plasma remains poorly understood. Spark gaps are thought to drive particle acceleration and pair creation in the black-hole magnetosphere. In this paper, we perform 2D axisymmetric general-relativistic particle-in-cell simulations of a monopole black-hole magnetosphere with a realistic treatment of inverse Compton scattering and pair production. We find that the magnetosphere can self-consistently fill itself with plasma and activate the Blandford-Znajek mechanism. A highly time-dependent spark gap opens near the inner light surface which injects pair plasma into the magnetosphere. These results may account for the high-energy activity observed in active galactic nuclei and explain the origin of plasma at the base of the jet.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gamma-ray emission from binaries in context

More than a dozen binary systems are now established as sources of variable, high energy (HE, 0.1-100 GeV) gamma rays. Five are also established sources of very high energy (VHE, >100 GeV) gamma rays. The mechanisms behind gamma-ray emission in binaries are very diverse. My current understanding is that they divide up into four types of systems: gamma-ray binaries, powered by pulsar rotation; microquasars, powered by accretion onto a black hole or neutron star; novae, powered by thermonuclear runaway on a white dwarf; colliding wind binaries, powered by stellar winds from massive stars. Some of these types had long been suspected to emit gamma rays (microquasars), others have taken the community by surprise (novae). My purpose here is to provide a brief review of the current status of gamma-ray emission from binaries, in the context of related objects where similar mechanisms are at work (pulsar wind nebulae, active galactic nuclei, supernova remnants).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modelling the high-energy emission from gamma-ray binaries using numerical relativistic hydrodynamics

Detailed modeling of the high-energy emission from gamma-ray binaries has been propounded as a path to pulsar wind physics. Fulfilling this ambition requires a coherent model of the flow and its emission in the region where the pulsar wind interacts with the stellar wind of its companion. We developed a code that follows the evolution and emission of electrons in the shocked pulsar wind based on inputs from a relativistic hydrodynamical simulation. The code is used to model the well-documented spectral energy distribution and orbital modulations from LS 5039. The pulsar wind is fully confined by a bow shock and a back shock. The particles are distributed into a narrow Maxwellian, emitting mostly GeV photons, and a power law radiating very efficiently over a broad energy range from X-rays to TeV gamma rays. Most of the emission arises from the apex of the bow shock. Doppler boosting shapes the X-ray and VHE lightcurves, constraining the system inclination to $i\approx 35^{\rm o}$. There is a tension between the hard VHE spectrum and the level of X-ray to MeV emission, which requires differing magnetic field intensities that are hard to achieve with a constant magnetisation $σ$ and Lorentz factor $Γ_{p}$ of the pulsar wind. Our best compromise implies $σ\approx 1$ and $Γ_{p}\approx 5\times 10^{3}$, respectively higher and lower than the typical values in pulsar wind nebulae. The high value of $σ$ derived here, where the wind is confined close to the pulsar, supports the classical picture that has pulsar winds highly magnetised at launch. However, such magnetisations will require further investigations to be based on relativistic MHD simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Shining in the Dark: the Spectral Evolution of the First Black Holes

Massive Black Hole (MBH) seeds at redshift $z \gtrsim 10$ are now thought to be key ingredients to explain the presence of the super-massive ($10^{9-10} \, \mathrm{M_{\odot}}$) black holes in place $ < 1 \, \mathrm{Gyr}$ after the Big Bang. Once formed, massive seeds grow and emit copious amounts of radiation by accreting the left-over halo gas; their spectrum can then provide crucial information on their evolution. By combining radiation-hydrodynamic and spectral synthesis codes, we simulate the time-evolving spectrum emerging from the host halo of a MBH seed with initial mass $10^5 \, \mathrm{M_{\odot}}$, assuming both standard Eddington-limited accretion, or slim accretion disks, appropriate for super-Eddington flows. The emission occurs predominantly in the observed infrared-submm ($1-1000 \, \mathrm{μm}$) and X-ray ($0.1 - 100 \, \mathrm{keV}$) bands. Such signal should be easily detectable by JWST around $\sim 1 \, \mathrm{μm}$ up to $z \sim 25$, and by ATHENA (between $0.1$ and $10 \, \mathrm{keV}$, up to $z \sim 15$). Ultra-deep X-ray surveys like the Chandra Deep Field South could have already detected these systems up to $z \sim 15$. Based on this, we provide an upper limit for the $z \gtrsim 6$ MBH mass density of $ρ_{\bullet} \lesssim 2.5 \times 10^{2} \, \mathrm{M_{\odot} \, Mpc^{-3}}$ assuming standard Eddington-limited accretion. If accretion occurs in the slim disk mode the limits are much weaker, $ρ_{\bullet} \lesssim 7.6 \times 10^{3} \, \mathrm{M_{\odot} \, Mpc^{-3}}$ in the most constraining case.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The case for super-critical accretion onto massive black holes at high redshift

Short-lived intermittent phases of super-critical (super-Eddington) growth, coupled with star formation via positive feedback, may account for early growth of massive black holes (MBH) and coevolution with their host spheroids. We estimate the possible growth rates and duty cycles of these episodes, both assuming slim accretion disk solutions, and adopting the results of recent numerical simulations. The angular momentum of gas joining the accretion disk determines the length of the accretion episodes, and the final mass a MBH can reach. The latter can be related to the gas velocity dispersion, and in galaxies with low-angular momentum gas the MBH can get to a higher mass. When the host galaxy is able to sustain inflow rates at 1-100 msunyr, replenishing and circulation lead to a sequence of short (~1e4-1e7 years), heavily obscured accretion episodes that increase the growth rates, with respect to an Eddington-limited case, by several orders of magnitude. Our model predicts that the ratio of MBH accretion rate to star formation rate is 1e2 or higher, leading, at early epochs, to a ratio of MBH to stellar mass higher than the "canonical" value of ~1e-3, in agreement with current observations. Our model makes specific predictions that long-lived super-critical accretion occurs only in galaxies with copious low-angular momentum gas, and in this case the MBH is more massive at fixed velocity dispersion.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

X-ray Transients: Hyper- or Hypo-Luminous?

The disk instability picture gives a plausible explanation for the behavior of soft X-ray transient systems if self-irradiation of the disk is included. We show that there is a simple relation between the peak luminosity (at the start of an outburst) and the decay timescale. We use this relation to place constraints on systems assumed to undergo disk instabilities. The observable X-ray populations of elliptical galaxies must largely consist of long-lived transients, as deduced on different grounds by Piro and Bildsten (2002). The strongly-varying X-ray source HLX-1 in the galaxy ESO 243-49 can be modeled as a disk instability of a highly super-Eddington stellar-mass binary similar to SS433. A fit to the disk instability picture is not possible for an intermediate-mass black hole model for HLX-1. Other, recently identified, super-Eddington ULXs might be subject to disk instability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cosmic ray induced ionisation of a molecular cloud shocked by the W28 supernova remnant

Cosmic rays are an essential ingredient in the evolution of the interstellar medium, as they dominate the ionisation of the dense molecular gas, where stars and planets form. However, since they are efficiently scattered by the galactic magnetic fields, many questions remain open, such as where exactly they are accelerated, what is their original energy spectrum, and how they propagate into molecular clouds. In this work we present new observations and discuss in detail a method that allows us to measure the cosmic ray ionisation rate towards the molecular clouds close to the W28 supernova remnant. To perform these measurements, we use CO, HCO$^+$, and DCO$^+$ millimetre line observations and compare them with the predictions of radiative transfer and chemical models away from thermodynamical equilibrium. The CO observations allow us to constrain the density, temperature, and column density towards each observed position, while the DCO$^+$/HCO$^+$ abundance ratios provide us with constraints on the electron fraction and, consequently, on the cosmic ray ionisation rate. Towards positions located close to the supernova remnant, we find cosmic ray ionisation rates much larger ($\gtrsim100$) than those in standard galactic clouds. Conversely, towards one position situated at a larger distance, we derive a standard cosmic ray ionisation rate. Overall, these observations support the hypothesis that the $γ$ rays observed in the region have a hadronic origin. In addition, based on CR diffusion estimates, we find that the ionisation of the gas is likely due to $0.1 -1$ GeV cosmic rays. Finally, these observations are also in agreement with the global picture of cosmic ray diffusion, in which the low-energy tail of the cosmic ray population diffuses at smaller distances than the high-energy counterpart.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gamma-ray binaries and related systems

After initial claims and a long hiatus, it is now established that several binary stars emit high (0.1-100 GeV) and very high energy (>100 GeV) gamma rays. A new class has emerged called 'gamma-ray binaries', since most of their radiated power is emitted beyond 1 MeV. Accreting X-ray binaries, novae and a colliding wind binary (eta Car) have also been detected - 'related systems' that confirm the ubiquity of particle acceleration in astrophysical sources. Do these systems have anything in common ? What drives their high-energy emission ? How do the processes involved compare to those in other sources of gamma rays: pulsars, active galactic nuclei, supernova remnants ? I review the wealth of observational and theoretical work that have followed these detections, with an emphasis on gamma-ray binaries. I present the current evidence that gamma-ray binaries are driven by rotation-powered pulsars. Binaries are laboratories giving access to different vantage points or physical conditions on a regular timescale as the components revolve on their orbit. I explain the basic ingredients that models of gamma-ray binaries use, the challenges that they currently face, and how they can bring insights into the physics of pulsars. I discuss how gamma-ray emission from microquasars provides a window into the connection between accretion--ejection and acceleration, while eta Car and novae raise new questions on the physics of these objects - or on the theory of diffusive shock acceleration. Indeed, explaining the gamma-ray emission from binaries strains our theories of high-energy astrophysical processes, by testing them on scales and in environments that were generally not foreseen, and this is how these detections are most valuable.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Particle acceleration and non-thermal emission during the V407 Cygni nova outburst

On March 2010, the symbiotic binary V407 Cyg erupted as a result of a nova explosion. The event gave rise to a two-week long burst of >100MeV gamma-rays detected by Fermi/LAT, a unique observation testifying to particle acceleration in the system. The outburst can be considered a scaled-down supernova, with short dynamical time scale, and thus can constitute a test case for theories of the origin of galactic cosmic rays. We aim at determining the properties of the accelerated particles and identifying the origin of the high-energy radiation. We developed a model for diffusive shock acceleration and non-thermal emission in V407 Cyg, complemented by an evaluation of the thermal emission from the shocked plasma. We considered both leptonic and hadronic contributions to the non-thermal processes, and investigated the effect of many binary and nova parameters. The gamma-ray emission is mostly of leptonic origin and arises predominantly from inverse-Compton scattering of the nova light. Upscattering of the red giant photons is a minor contribution. Matching the light curve requires gas accumulation in the vicinity of the white dwarf, as a consequence of wind accretion, while the spectrum imposes particle scattering close to the Bohm limit in the upstream equipartition magnetic field. The nova accelerated protons (respectively electrons) with energies up to ~300GeV (respectively ~20GeV), for a total non-thermal energy 10^43 erg after two weeks, representing ~10% of the initial nova kinetic energy. The electron-to-proton ratio at injection is 6%. The V407 Cyg eruption can be understood from the same principles that are invoked for particle acceleration in supernova remnants, although without the need for strong magnetic field amplification. The population of novae in symbiotic systems is a negligible source of Galactic cosmic rays, and most likely not a class of TeV-emitters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulating gamma-ray binaries with a relativistic extension of RAMSES

Gamma-ray binaries are composed of a massive star and a rotation-powered pulsar with a highly relativistic wind. The collision between the winds from both objects creates a shock structure where particles are accelerated, resulting in the observed high energy emission. We study the impact of special relativity on the structure and stability of the colliding wind region and highlight the differences with colliding winds from massive stars. We focus on evolution with increasing values of the Lorentz factor of the pulsar wind, keeping in mind that current simulations are unable to reach the expected values of the pulsar wind Lorentz factors by orders of magnitude. We use high resolution numerical simulations with a relativistic extension to the hydrodynamics code RAMSES we have developed. Using 2D simulations, we focus on the region close to the binary, neglecting orbital motion. We use different values of the Lorentz factor of the pulsar wind, up to 16. We find analytic scaling relations between stellar wind collisions and gamma-ray binaries. They provide the position of the contact discontinuity. The position of the shocks strongly depends on the Lorentz factor, the relativistic wind is more collimated than expected based on non-relativistic simulations. Beyond a certain distance, the shocked flow is accelerated to its initial velocity and follows adiabatic expansion. We provide guidance for extrapolation towards more realistic values of the Lorentz factor of the pulsar wind. We extended the adaptive mesh refinement code RAMSES to relativistic hydrodynamics. This code is suited for the study of astrophysical objects such as pulsar wind nebulae, gamma-ray bursts or relativistic jets and will be part of the next public release of RAMSES. Using this code we performed simulations of gamma-ray binaries, highlighting the limits and possibilities of current hydrodynamic models of such systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Colliding wind binaries and gamma-ray binaries : relativistic version of the RAMSES code

Gamma-ray binaries are colliding wind binaries (CWB) composed of a massive star a non-accreting pulsar with a highly relativistic wind. Particle acceleration at the shocks results in emission going from extended radio emission to the gamma-ray band. The interaction region is expected to show common features with stellar CWB. Performing numerical simulations with the hydrodynamical code RAMSES, we focus on their structure and stability and find that the Kelvin-Helmholtz instability (KHI) can lead to important mixing between the winds and destroy the large scale spiral structure. To investigate the impact of the relativistic nature of the pulsar wind, we extend RAMSES to relativistic hydrodynamics (RHD). Preliminary simulations of the interaction between a pulsar wind and a stellar wind show important similarities with stellar colliding winds with small relativistic corrections.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Impact of orbital motion on the structure and stability of adiabatic shocks in colliding wind binaries

The collision of winds from massive stars in binaries results in the formation of a double-shock structure with observed signatures from radio to X-rays. We study the structure and stability of the colliding wind region as it turns into a spiral due to orbital motion. We focus on adiabatic winds, where mixing between the two winds is expected to be restricted to the Kelvin-Helmholtz instability (KHI). Mixing of the Wolf-Rayet wind with hydrogen-rich material is important for dust formation in pinwheel nebulae such as WR 104, where the spiral structure has been resolved in infrared. We use the hydrodynamical code RAMSES with an adaptive grid. A wide range of binary systems with different wind velocities and mass loss rates are studied with 2D simulations. A specific 3D simulation is performed to model WR 104. Orbital motion leads to the formation of two distinct spiral arms where the KHI develops differently. We find that the spiral structure is destroyed when there is a large velocity gradient between the winds, unless the collimated wind is much faster. We argue that the KHI plays a major role in maintaining or not the structure. We discuss the consequences for various colliding wind binaries. When stable, there is no straightforward relationship between the spatial step of the spiral, the wind velocities, and the orbital period. Our 3D simulation of WR 104 indicates that the colder, well-mixed trailing arm has more favourable conditions for dust formation than the leading arm. The single-arm infrared spiral follows more closely the mixing map than the density map, suggesting the dust-to-gas ratio may vary between the leading and trailing density spirals. However, the density is much lower than what dust formation models require. Including radiative cooling would lead to higher densities, and also to thin shell instabilities whose impact on the large structure remains unknown.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modelling colliding wind binaries with RAMSES, extension to special relativity

We present high resolution simulations with RAMSES of supersonic colliding stellar winds. The collision results in a double shock structure which is subject to different instabilities. The Kelvin-Helmholtz instability (KHI) introduces some mixing and variability. For isothermal winds, the Non-linear Thin Shell Instability violently affects the interaction region. Properly modelling these instabilities requires a high enough resolution and an adapted numerical method, especially when one of the winds strongly dominates the other one. At large scale, orbital motion is expected to turn the shocked zone into a spiral but we find that in some configurations the KHI may disrupt the spiral. A colliding wind structure is also expected in gamma-ray binaries composed of a massive star and a young pulsar which emits a highly relativistic wind. Numerical simulations are necessary to understand the geometry of such systems and should take into account the relativistic nature of the pulsar wind. We implemented a second order Godunov method to solve the equations of relativistic hydrodynamics with RAMSES, including the possibility for Adaptive Mesh Refinement. After a brief overview of our numerical implementation, we will present preliminary simulation that contrast the structure of the flow in gamma-ray and stellar binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Models of AM CVn star outbursts

Outbursting AM CVn stars exhibit outbursts similar to those observed in different types of dwarf novae. Their light-curves combine the characteristic features of SU UMa, ER UMa, Z Cam, and WZ Sge-type systems but also show a variety of properties never observed in dwarf novae. The compactness of AM CVn orbits and their unusual chemical composition make these systems valuable testbeds for outburst models. We aim for a better understanding of the role of helium in the accretion disc instability mechanism, testing the model for dwarf novae outbursts in the case of AM CVn stars, and aim to explain the outburst light-curves of these ultra-compact binaries. We calculated the properties of the hydrogen-free AM CVn stars using our previously developed numerical code adapted to the different chemical composition of these systems and supplemented with formulae accounting for mass transfer rate variations, additional sources of the disc heating, and the primary's magnetic field. We discovered how helium-dominated discs react to the thermal-viscous instability and were able to reproduce various features of the outburst cycles in the light-curves of AM CVn stars. The AM CVn outbursts can be explained by the suitably adapted dwarf-nova disc instability model but, as in the case of its application to hydrogen-dominated cataclysmic variables, one has to resort to additional mechanisms to account for the observed superoutbursts, dips, cycling states, and standstills. We show that the enhanced mass-transfer rate, due presumably to variable irradiation of the secondary, must not only be taken into account but is a determining factor that shapes AM CVn star outbursts. The cause of the variable secondary's irradiation has yet to be understood; the best candidate is the precession of a tilted/warped disc.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Pulsed high energy gamma-rays from thermal populations in the current sheets of pulsar winds

Context. More than one hundred GeV pulsars have been detected up to now by the LAT telescope on the Fermi gamma-ray observatory, showing peak energies around a few GeV. Current modelling proposes that the high energy emission comes from outer magnetospheric gaps, however radiation from the equatorial current sheet which separates the two magnetic hemispheres outside the light cylinder has also been investigated. Aims. In this paper we discuss the region right outside the light cylinder, or "near wind" zone. We investigate the possibility that synchrotron radiation emitted by thermal populations in the equatorial current sheet of the pulsar wind in this region can explain the lightcurves and spectra observed by Fermi/LAT. Methods. We use analytical estimates as well as detailed numerical computation to calculate the gamma-ray luminosities, lightcurves and spectra of gamma-ray pulsars. Results. Many of the characteristics of the gamma-ray pulsars observed by Fermi/LAT can be reproduced by our model, most notably the position of these objects in the P - Pdot diagram, and the range of gamma-ray luminosities. A testable result is a sub-exponential cutoff with an index b = 0.35. We also predict the existence of a population of pulsars with cutoff energies in the MeV range. These have systematically lower spindown luminosities than the Fermi/LAT detected pulsars. Conclusions. It is possible for relativistic populations of electrons and positrons in the current sheet of a pulsar's wind right outside the light cylinder to emit synchrotron radiation that peaks in the sub-GeV to GeV regime, with gamma-ray efficiencies similar to those observed for the Fermi/LAT pulsars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The gamma-ray emitting region of the jet in Cyg X-3

We study models of the gamma-ray emission of Cyg X-3 observed by Fermi. We calculate the average X-ray spectrum during the gamma-ray active periods. Then, we calculate spectra from Compton scattering of a photon beam into a given direction by isotropic relativistic electrons with a power-law distribution, both based on the Klein-Nishina cross section and in the Thomson limit. Applying the results to scattering of stellar blackbody radiation in the inner jet of Cyg X-3, we find that a low-energy break in the electron distribution at a Lorentz factor of ~ 300--1000 is required by the shape of the observed X-ray/gamma-ray spectrum in order to avoid overproducing the observed X-ray flux. The electrons giving rise to the observed \g-rays are efficiently cooled by Compton scattering, and the power-law index of the acceleration process is ~ 2.5--3. The bulk Lorentz factor of the jet and the kinetic power before the dissipation region depend on the fraction of the dissipation power supplied to the electrons; if it is ~ 1/2, the Lorentz factor is ~ 2.5, and the kinetic power is ~ 10^38 erg/s, which represents a firm lower limit on the jet power, and is comparable to the bolometric luminosity of Cyg X-3. Most of the power supplied to the electrons is radiated. The broad band spectrum constrains the synchrotron and self-Compton emission from the gamma-ray emitting electrons, which requires the magnetic field to be relatively weak, with the magnetic energy density < a few times 10^-3 of that in the electrons. The actual value of the magnetic field strength can be inferred from a future simultaneous measurement of the IR and gamma-ray fluxes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Absorption of high-energy gamma rays in Cygnus X-3

The microquasar Cygnus X-3 was detected at high energies by the gamma-ray space telescopes AGILE and Fermi. The gamma-ray emission is transient, modulated with the orbital period and seems related to major radio flares, i.e. to the relativistic jet. The GeV gamma-ray flux can be substantially attenuated by internal absorption with the ambient X-rays. In this study, we examine quantitatively the effect of pair production in Cygnus X-3 and put constraints on the location of the gamma-ray source. Cygnus X-3 exhibits complex temporal and spectral patterns in X-rays. During gamma-ray flares, the X-ray emission can be approximated by a bright disk black body component and a non-thermal tail extending in hard X-rays, possibly related to a corona above the disk. We calculate numerically the exact optical depth for gamma rays above a standard accretion disk. Emission and absorption in the corona are also investigated. GeV gamma rays are significantly absorbed by soft X-rays emitted from the inner parts of the accretion disk. The absorption pattern is complex and anisotropic. Isotropization of X-rays due to Thomson scattering in the companion star wind tends to increase the gamma-ray opacity. Gamma rays from the corona suffer from strong absorption by photons from the disk and cannot explain the observed high-energy emission, unless the corona is unrealistically extended. The lack of absorption feature in the GeV emission indicates that high-energy gamma rays should be located at a minimum distance ~10^8-10^10 cm from the compact object. The gamma-ray emission is unlikely to have a coronal origin.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

High resolution numerical simulations of unstable colliding stellar winds

We investigate the hydrodynamics of the interaction of two supersonic winds in binary systems. The collision of the winds creates two shocks separated by a contact discontinuity. The overall structure depends on the momentum flux ratio eta of the winds. We use the code RAMSES with adaptive mesh refinement to study the shock structure up to smaller values of eta, higher spatial resolution and greater spatial scales than have been previously achieved. 2D and 3D simulations, neglecting orbital motion, are compared to widely-used analytic results and their applicability is discussed. In the adiabatic limit, velocity shear at the contact discontinuity triggers the Kelvin-Helmholtz instability. We quantify the amplitude of the resulting fluctuations and find that they can be significant even with a modest initial shear. Using an isothermal equation of state leads to the development of thin shell instabilities. The initial evolution and growth rates enables us to formally identify the non-linear thin shell instability (NTSI) close to the binary axis. Some analogue of the transverse acceleration instability is present further away. The NTSI produces large amplitude fluctuations and dominates the long-term behaviour. We point out the computational cost of properly following these instabilities. Our study provides a basic framework to which the results of more complex simulations, including additional physical effects, can be compared.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Implication of the striped pulsar wind model for gamma-ray binaries

(abridged) Gamma-ray binaries are massive stars with compact object companions that are observed to emit most of their energy in the gamma-ray range. One of these binaries is known to contain a radio pulsar, PSR B1259-63. Synchrotron and inverse Compton emission from particles accelerated beyond the light cylinder in striped pulsar winds has been proposed to explain the X-ray to high energy (HE, $>$ 100 MeV) gamma-ray emission from isolated pulsars. This pulsar model extends naturally to binary environments, where seed photons for inverse Compton scattering are provided by the companion star. Here, we investigate the possibility of gamma-ray emission from PSR B1259-63 in the framework of the striped pulsar wind model. The orbital geometry of PSR B1259-63 is well constrained by observations and the double radio pulse suggests an almost orthogonal rotator so that the solid angle covered by the striped region is close to $4π$. We calculate the orbital and rotational phase-resolved spectral variability and light-curves to expect.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Near-infrared jet emission in the microquasar XTE J1550-564

Context: Microquasars are accreting Galactic sources that are also observed to launch relativistic jets. A key signature of the ejection is non-thermal radio emission. The level of this jet component at high frequencies is still poorly constrained. Aims: The X-ray binary and microquasar black hole candidate XTE J1550-564 exhibited a faint X-ray outburst in April 2003 during which it stayed in the X-ray low/hard state. We took optical and near-infrared (NIR) observations with the ESO/NTT telescope during this outburst to disentangle the various contributions to the spectral energy distribution (SED) and investigate the presence of a jet component. Methods: Photometric and spectroscopic observations allowed us to construct an SED and also to produce a high time-resolution lightcurve. Results: The SED shows an abrupt change of slope from the NIR domain to the optical. The NIR emission is attributed to non-thermal synchrotron emission from the compact, self-absorbed jet that is known to be present in the low/hard state. This is corroborated by the fast variability, colours, lack of prominent spectral features and evidence for intrinsic polarisation. The SED suggests the jet break from the optically thick to the thin regime occurs in the NIR. Conclusions: The simultaneous optical-NIR data allow an independent confirmation of jet emission in the NIR. The transition to optically thin synchrotron occurs at NIR frequencies or below, which leads to an estimated characteristic size greater than 2x10^8cm and magnetic field less than 5T for the jet base, assuming a homogeneous one-zone synchrotron model.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

High-energy radiation from the relativistic jet of Cygnus X-3

Cygnus X-3 is an accreting high-mass X-ray binary composed of a Wolf-Rayet star and an unknown compact object, possibly a black hole. The gamma-ray space telescope Fermi found definitive evidence that high-energy emission is produced in this system. We propose a scenario to explain the GeV gamma-ray emission in Cygnus X-3. In this model, energetic electron-positron pairs are accelerated at a specific location in the relativistic jet, possibly related to a recollimation shock, and upscatter the stellar photons to high energies. The comparison with Fermi observations shows that the jet should be inclined close to the line of sight and pairs should not be located within the system. Energetically speaking, a massive compact object is favored. We report also on our investigations of the gamma-ray absorption of GeV photons with the radiation emitted by a standard accretion disk in Cygnus X-3. This study shows that the gamma-ray source should not lie too close to the compact object.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modeling the three-dimensional pair cascade in binaries. Application to LS 5039

LS 5039 is a Galactic binary system emitting high and very-high energy gamma rays. The gamma-ray flux is modulated on the orbital period and the TeV lightcurve shaped by photon-photon annihilation. The observed very-high energy modulation can be reproduced with a simple leptonic model but fails to explain the flux detected by HESS at superior conjunction, where gamma rays are fully absorbed. The contribution from an electron-positron pair cascade could be strong and prevail over the primary flux at superior conjunction. The created pairs can be isotropized by the magnetic field, resulting in a three-dimensional cascade. The aim of this article is to investigate the gamma ray radiation from this pair cascade in LS 5039. This additional component could account for HESS observations at superior conjunction in the system. A semi-analytical and a Monte Carlo method for computing three-dimensional cascade radiation are presented and applied in the context of binaries. Three-dimensional cascade radiation contributes significantly at every orbital phase in the TeV lightcurve, and dominates close to superior conjunction. The amplitude of the gamma-ray modulation is correctly reproduced for an inclination of the orbit of about 40 degrees. Primary pairs should be injected close to the compact object location, otherwise the shape of the modulation is not explained. In addition, synchrotron emission from the cascade in X-rays constrains the ambient magnetic field to below 10 G. The radiation from a three-dimensional pair cascade can account for the TeV flux detected by HESS at superior conjunction in LS 5039, but the very-high energy spectrum at low fluxes remains difficult to explain in this model.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Reconciling observed GRB prompt spectra with synchrotron radiation ?

(abridged)Prompt GRB emission is often interpreted as synchrotron radiation from high-energy electrons accelerated in internal shocks. Fast synchrotron cooling predicts that the photon index below the spectral peak is alpha=-3/2. This differs significantly from the observed median value alpha \approx -1. We quantify the influence of inverse Compton and adiabatic cooling on alpha to understand whether these processes can reconcile the observations with a synchrotron origin. We use a time-dependent code that follows both the shock dynamics and electron energy losses. We investigate the dependence of alpha on the parameters of the model. Slopes between -3/2 and -1 are reached when electrons suffer IC losses in the Klein-Nishina regime. This does not necessarily imply a strong IC component in the Fermi/LAT range because scatterings are only moderately efficient. Steep slopes require that a large fraction (10-30%) of the dissipated energy is given to a small fraction (<~1%) of the electrons and that the magnetic energy density fraction remains low (<~ 0.1%). Values of alpha up to -2/3 can be obtained with relatively high radiative efficiencies (>50%) when adiabatic cooling is comparable with radiative cooling (marginally fast cooling). This requires collisions at small radii and/or with low magnetic fields. Amending the standard fast cooling scenario to account for IC cooling naturally leads to alpha up to -1. Marginally fast cooling may also account for alpha up to -2/3, although the conditions required are more difficult to reach. About 20% of GRBs show spectra with slopes alpha>-2/3. Other effects, not investigated here, such as a thermal component in the electron distribution or pair production by HE photons may further affect alpha. Still, the majority of observed GRB prompt spectra can be reconciled with a synchrotron origin, constraining the microphysics of mildly relativistic internal shocks.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

X-ray absorption and occultation in LS 5039

Gamma-ray binaries are systems containing a massive star and a compact object that have been detected up to TeV energies. The high energy emission could result from particle acceleration in the region where the stellar wind from the massive star interacts with the relativistic wind from a young pulsar. LS 5039 has the most compact orbit amongst gamma-ray binaries and its X-ray lightcurve shows a stable modulation synchronized with the orbital period. Photoelectric absorption of X-rays in the O star wind and occultation of the X-ray emitting region by the massive star can alter the X-ray lightcurve and spectrum along the orbit. Yet, the X-ray spectrum and lightcurve of LS 5039 do not show intrinsic absorption or X-ray eclipses. We study these effects in the framework of the pulsar wind scenario as a function of the binary inclination angle, the stellar wind mass-loss rate and the size of the X-ray emitter. An extended X-ray emission region >~ 3R_star appears necessary to reconcile the pulsar wind scenario with observations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The outbursts properties of AM CVn stars

We briefly summarize the observational properties of ultra-compact binaries called AM CVn stars. We analize their outbursts originating from the thermal-viscous instability in helium accretion disc. We present our preliminary results in applying the model of Dwarf Novae outbursts to helium discs. We can calculate models of outbursts of reasonable amplitude of 2 mag with a constant alfa parameter throughout the calculation. Setting the mass transfer rate close to its upper critical value produces model lightcurves that resemble short superoutbursts.

0 citations0 reviewsNo code linkedOpen access