Projects of the Sherpas team

The sherpas group is involved in theoretical and observational studies of accretion and ejection around young stars, compact objects as well as high energy processes : relativistic jets, high energy emission, gamma ray bursts. Our main activities are: theoretical studies and numerical simulations of magnetised plasmas, studies of acceleration and emission processes of high and very high energy particles (X rays, gamma rays, cosmic rays), and participation in instrumental collaborations: satellites observing in X-ray (XMM, NuSTAR, Chandra, IXPE) and gamma-ray (Fermi), Cherenkov telescopes (H.E.S.S, CTA). Our activities fall within the themes selected by the ASTRONET network, designed to direct the 5-25 years prospective of research in astrophysics of the European Union: do we understand the extremes of the Universe?

Observation of high-energy emission

Motivations

Observations of the high-energy sky reveal many compact objects, Galactic and extragalactic black holes, neutron stars and white dwarfs, and allow us to study the radiation processes associated with these objects and their extreme environment. Observational diagnostics (morphology, variability, flux and spectrum) reveal the physical phenomena at work in these objects and are essential to test the theoretical predictions developed in our team.

(Left) Combined radio/optical/X image of the Centaurus A Active Galaxy Nucleus. (Right) Artist’s views of the inner regions of an active galactic nucleus (bottom) and of an X-ray binary (top). These two systems show the same emission regions (disk, corona, wind, jet) while the difference in the black hole mass and size is generally greater than six or nine orders of magnitude.

Objectives

Although compact objects have been studied for several decades, many questions remain unanswered:

  • What is the geometry of the central regions close to the compact object? What are their dynamics?
  • What are the radiative processes at the origin of the observed high-energy emission? What is the origin of its variability?
  • What can we deduce from the signatures of winds and jets on the ejection processes ? on the accretion-ejection links ?
  • What are the observational differences/ similarities between active galactic nuclei and X-ray binaries?
  • What can we learn from the populations of detected X-ray and Gamma-ray sources?

Methodology

We access the physical conditions close to the compact object through multi-wavelength observations, covering the whole electromagnetic spectrum from radio to gamma rays. Our group is more specifically involved, as PI and co-I, in high energy observations with high energy satellites like XMM-Newton, Chandra, NuSTAR and IXPE.

Contacts

Pierre-Olivier Petrucci, Maïca Clavel, Guillaume Dubus, Gilles Henri, Dilruwan Dehiwalage Don, Vittoria Gianolli, Maxime Parra

Kinetic processes and particle acceleration

Overview

There is overwhelming observational evidence that compact objects (white dwarfs, neutron stars, black holes) are involved in some of the most extreme particle acceleration phenomena in the Universe. Particle spectra escaping from these cosmic accelerators are broad steep power laws extending to tremendous energies, many orders of magnitude above the particles rest mass energies (i.e., they are ultra-relativistic). There are currently two main acceleration mechanisms under scrutiny by the community: diffusive shock acceleration and magnetic reconnection.

Global 3D PIC simulation of an inclined pulsar magnetosphere (Cerutti et al. 2016)
Full general relativistic PIC simulation of a rotating (Kerr) black hole magnetosphere (Crinquand et al. 2021)

Particle acceleration in relativistic anisotropic magnetized shocks (Cerutti & Giacinti 2020)

Objectives

The work of the team in this field tries to answer the following fundamental questions:

  • How are particles accelerated within relativistic magnetic reconnection sites?
  • How and where are particle accelerated within neutron star and black hole magnetospheres?
  • How are particle accelerated at relativistic, magnetized collisionless shocks?
  • What are the observable signatures?

Astrophysical objects of interest

  • Pulsars, winds and pulsar wind nebulae
  • Stellar mass and supermassive black holes (X-ray binaries, AGN, Sgr A*)
  • Relativistic jets, blazars
  • Gamma-ray bursts
  • Nova outbursts

Methods

Studying particle acceleration requires a complete description of plasma phenomena at the kinetic scale, i.e., below the particle Larmor radius scale. The particle-in-cell (PIC) technique is a numerical method well-suited for this purpose: it solves the motion of a large collection of charged particles along with the evolution of the electromagnetic fields and their coupling between them. The team uses and develops the Zeltron code, a state-of-the-art PIC code optimized for studying astrophysical relativistic plasmas.

The ERC-SPAWN project

The ERC-SPAWN project (Simulating Particle Acceleration WithiN black hole magnetospheres, 2020-2025) aims at producing an ab-initio model of black hole magnetospheres. This work is essential to interpret current and future horizon-scale observations by Gravity@VLTI and the Event Horizon Telescope of nearby supermassive black holes.

Contacts

Benoît Cerutti, Guillaume Dubus, John Mehlhaff, Enzo Figueiredo, Valentina Richard Romei, Adrien Soudais

Modelling accretion-ejection phenomena

Motivations

The team’s work on accretion-ejection phenomena seeks to model and interpret the emission of astrophysical objects with accretion disks. We now know that accretion phenomena are almost systematically associated with ejection, in the form of jets or winds. A specificity of the Sherpas team is to focus on the interaction between the large-scale magnetic field and the accretion disc to explain this link, using both phenomenological models, observations and high-performance numerical simulations.

Objectives

The work in this area seeks to answer the following questions:

  • By what mechanism do these disks accrete (hydro/MHD turbulence, ejection...)?
  • What are the links between ejection and accretion phenomena?
  • How does the magnetic field evolve in these objects and what would be the associated signatures?
  • What are the links between these accretion-ejection processes and the radiative processes occurring in the environment of compact objects (X-binaries, AGNs, cataclysmic variables, FuOr)?
Model of a disc forming in a compact binary system (dwarf novae/X-binary). The secondary star is on the right and is losing material that spirals around the primary star in the centre of the simulation.

Methods

For the modelling part, we use both semi-analytical solutions (self-similar solutions), phenomenological models of disks, and direct numerical simulations in the magnetohydrodynamic approximation (IDEFIX, PLUTO and SNOOPY codes).

Modelling of a turbulent disc with the PLUTO code (seen from the edge) showing both accretion and matter ejection, PLUTO code (Jacquemin-Ide et al. 2021).



Semi-analytical models allow an easier comparison with observations. This is for example the case of the self-similar Jet Emitting Disk (JED) solutions developed in the sherpas team for years. They can be compared to X-ray observations of compact objects such as AGNs or X-ray binaries.

Fitting of the X-ray spectra of the X-binary GX 339-4 during its 2010-2011 outburst with the ejecting disk model. The black markers are the observed values (Clavel et al. 2016) while the coloured lines display the results of the fits: green, blue, yellow and red for the hard, hard-intermediate, soft-intermediate and soft states, respectively. Left: Light curves. Right: Disk Fraction Luminosity Diagram (see Marcel et al. 2019).

Contacts

Guillaume Dubus, Jonathan Ferreira, Geoffroy Lesur, Pierre-Olivier Petrucci, Nicolas Scepi, Marc Van den Bossche, Nathan Zimniak

Dynamics of protoplanetary disks

Motivations

Protoplanetary disks are cold, low ionisation environments. They can nevertheless be described as a non-ideal plasma, subject to a coupling to the magnetic field depending on the location inside the disc. It is now thought that this coupling is necessary to explain the accretion rates in protoplanetary disks but also the molecular winds that are almost systematically observed.

Objectives

The objective of this research theme is to develop non-ideal plasma models of these disks from first principles of physics, in order to deduce observational signatures of the dynamics of these objects. In particular, we aim to:

  • Understand how these objects accrete, despite their very low ionisation fraction
  • Understand how the winds interact with the dust in the disc (sedimentation, entrainment and transport)
  • Study the interaction between the plasma and the forming planets still buried in the disc, in particular the migration of the planets, depending on their mass.
Flow and density map of the outflow in the vicinity of a planet with a mass of Jupiter in formation (on the right side of the figure) in a magnetised disc. Model by G. Wafflard-Fernandez, code IDEFIX.
Cross-section of the cavity of a magnetised transition disc, code PLUTO

Methods

To study these phenomena, we use numerical methods adapted to the solution of the compressible magnetohydrodynamics (MHD) equations. In practice, we use the Idefix code, which is a new generation Godunov finite volume code, adapted to massively parallel accelerated super-computers (e.g. using GPUs) such as Jean Zay at Idris or Adastra at CINES.

The ERC-MHDiscs project

The ERC MHDiscs project aims to perform ab-initio simulations of protoplanetary disks taking into account the interaction with the surrounding magnetic field and the low degree of ionisation of these objects. In particular, we aim to obtain observable quantitative predictions that would indicate the presence of planets buried in these discs, or magnetohydrodynamic phenomena that could be precursors to planetary formation.

Contacts

Geoffroy Lesur, Nicolas Scepi, Gaylor Wafflard-Fernandez, Jonah Mauxion

Astrocladistics and astrostatistics

Motivations

Big data are an inescapable reality in astrophysics today, and the adapted tools are starting to spread and develop: Statistics, Big Data, Data Mining, Machine Learning, Deep Learning, AI (Artificial Intelligence), these terms are now part of the landscape of astrophysicists. Still in its infancy a decade ago, astrostatistics now has its own dedicated journal that Didier Fraix-Burnet created in 2019.

A particularly important area is the real-time classification of new observations because the rate and volume of data acquisition by large modern telescopes and detectors no longer allow for visual or manual analysis. For this purpose, astronomers have mainly tried to make the algorithms reproduce our knowledge (supervised classification), mainly through neural networks or Deep Learning. This allows a priori to classify new observations very quickly in the already known classes. On the condition of having a very large learning base.

Objectives

Understand the power of unsupervised classification, i.e. the ability of the algorithm to build its learning base by itself. The objective of unsupervised classification is to search for structures in a data space with n objects and p variables. These structures characterize either similarities (groups or classes) or relationships such as evolutionary paths. The former are rather studied by statistical methods (machine learning, data mining), the latter are based on graph theory and are the prerogative of phylogenetic methods (path minimization).

Methods

Astrocladistics
Initiated in 2001 by Didier Fraix-Burnet, astrocladistics concerns the use of phylogenetic methods in astrophysics, cladistics, also called Maximum Parsimony, being one of them but certainly the most general. Initially thought for the case of galaxy diversification, astrocladistics has been successfully applied on all kinds of astrophysical entities, such as stars, stellar clusters (globular and open), and even Jupiter and Saturn satellites! To read more about astrocladistics, visit D. Fraix-Burnet’s dedicated web-blog.

Cladogram of 14 Dwarf galaxies of the Local Group obtained with 24 characters (observables and derived quantities). Bootstrap values (above) and Decay indices (below) are indicated for each node. The outgroup (SagDig) has been chosen because it contains the lowest amount of metallic material, suggesting that it is made up of more primordial material (see Fraix-Burnet et al. 2006). Author D. Fraix-Burnet. Licence CC-BY. Image credits.
Physical interpretation of the phylogenetic classification established from a thousand galaxies (Fraix-Burnet et al. 2012).

Astrostatistics
The unsupervised approach that we implement is based on the technique of mixing Gaussian distributions in a latent subspace. In other words, we look for groupings of various Gaussian distributions not in the space of the variables, but in a subspace where these groupings are clearer. The algorithm used is called Fisher-EM and was developed by Charles Bouveyron and is available under R.
We applied it to a sample of 702 248 optical spectra of galaxies (from the SDSS), each spectrum containing 1437 monochromatic fluxes. We find an optimum of 86 classes, of which the 37 largest cover 99% of the sample (Fraix-Burnet et al. 2021). A first quick analysis shows the physical relevance of these classes in classically used diagrams.

Average spectra and dispersion for the 86 classes
Distribution in the color-magnitude diagram of the 86 classes

Julien Dubois has confirmed the ability of this algorithm to discover physically relevant classes from simulated spectra, and he is working on proposing a fine interpretation of each of the SDSS classes in order to create a true objective atlas of galaxy spectra, and on studying the evolution of this classification using higher redshift samples.

Contacts

Didier Fraix-Burnet, Hugo Chambon