Gyrokinetic modelling of the nonlinear interaction between energetic particle-driven instabilities and microturbulence in tokamak plasmas

Tokamak plasmas are strongly nonlinear systems far from thermodynamic equilibrium, in which instabilities of very different spatial scales coexist, ranging from large-scale macroscopic oscillations to microturbulence. The presence of energetic ions produced by fusion reactions or by auxiliary heating further enhances these instabilities through wave–particle resonances. Microturbulence is responsible for heat and particle transport in the thermal plasma, while instabilities driven by energetic particles can induce their radial transport and, consequently, their losses. Both phenomena degrade the performance of present tokamak plasmas, and possibly also those of burning plasmas such as ITER.
Recent results, however, show that these instabilities, which have long been studied separately, can interact nonlinearly, and that this interaction may lead to an unexpected improvement of plasma confinement.
The objective of this project is to investigate these multiscale interactions using the gyrokinetic code GTC, which is able to simultaneously simulate turbulence and energetic-particle-driven instabilities. This work aims to improve the understanding of the nonlinear mechanisms governing plasma confinement and to identify optimal regimes for future fusion plasmas.

Impact of irradiation parameters on the alpha’ phase formation in oxide dispersion strengthened steels

Ferritic-martensitic oxide dispersion strengthened steels (ODS steels) are materials of great interest in the nuclear industry. Predominantly composed of iron and chromium, these materials can become brittle due to the precipitation of a chromium-rich phase, called a', under irradiation. This phase, known to be sensitive to irradiation conditions, provides an ideal topic for a deeper exploration of the capability to emulate neutron irradiation with ions. Indeed, while ion irradiations are frequently used to understand phenomena observed during neutron irradiations, the question of their representativeness is often raised.

In this thesis, we aim to understand how the irradiation parameters can affect the characteristics of the a' phase in ODS steels. To do so, various ODS steels will be irradiated under different conditions (flux, dose, temperature, and type of particles, such as ions, neutrons, electrons), and subsequently analyzed at the nanoscale. The a' phase (size, chromium content) obtained for each ion irradiation condition will be compared to the one after neutron irradiation.

Applications using laser-accelerated relativistic electrons with PETAL

This PhD project focuses on the physics of plasmas generated by ultra-high-power and high-intensity lasers. The work will be carried out at the LMJ facility, using the PETAL laser which operates at intensities exceeding 10¹8 W·cm?² and enables the production of high-energy particles.
The main objective of the thesis is to investigate the generation and acceleration of relativistic electron beams in a gas jet. The potential applications of these beams will be assessed for electron–positron pair production and for electron-beam-based radiography.
The research will combine experimental and numerical approaches. The PhD candidate will take part in experimental campaigns scheduled for 2026–2027, including the implementation of diagnostics and data analysis. In parallel, Particle-In-Cell and Monte Carlo simulations will be performed to support the interpretation of the experimental results.
In a second phase, the thesis will contribute to the qualification of upgrades to the PETAL laser, focusing in particular on secondary sources of electrons, protons, and hard X-ray radiation generated by laser–matter interactions, within the framework of the PETAL-UPGRADE project.

Preconditioning of iterative schemes for the mixed finite element solution of an eigenvalue problem applied to neutronics

Neutronics is the study of the behavior of neutrons in matter and the reactions they induce, particularly the generation of power through the fission of heavy nuclei. Modeling the steady-state neutron flux in a reactor core relies on solving a generalized eigenvalue problem of the form:
Find (phi, keff) such that A phi=1/keff B phi and keff is the eigenvalue with the largest magnitude, where A is the disappearance matrix which is assumed invertible, B represents the production matrix, phi denotes the neutron flux, and keff is called the multiplication factor.

The neutronics code APOLLO3® is a joint project of CEA, Framatome, and EDF for the development of a next-generation code for reactor core physics to meet both R&D and industrial application needs [4].
The MINOS solver [2] is developed within the framework of the APOLLO3® project. This solver is based on the mixed finite element discretization of the neutron diffusion model or the simplified transport model. The strategy for solving the aforementioned generalized eigenvalue problem is iterative; it involves applying the inverse power method [6].

The convergence speed of this inverse power method algorithm depends on the spectral gap. In the context of large cores such as the EPR reactor, it is observed that the spectral gap is close to 1, which degrades the convergence of the inverse power method algorithm. It is necessary to apply acceleration techniques to reduce the number of iterations [7]. In neutron transport, the preconditioning called Diffusion Synthetic Acceleration is very popular for the so-called inner iteration [1] but has also recently been applied to the so-called outer iteration [3]. A variant of this method was introduced in [5] for solving a source problem. It is theoretically shown that this variant converges in all physical regimes.

[1] M. L. Adams, E. W. Larsen, Fast iterative methods for discrete-ordinates particle transport calculations, Progress in Nuclear Energy, Volume 40, Issue 1, 2002.

[2] A.-M. Baudron and J.-J. Lautard. MINOS: a simplified PN solver for core calculation. Nuclear Science and Engineering, volume 155(2), pp. 250–263 (2007).

[3] A. Calloo, R. Le Tellier, D. Couyras, Anderson acceleration and linear diffusion for accelerating the k-eigenvalue problem for the transport equation, Annals of Nuclear Energy, Volume 180, 2023.

[4] P. Mosca, L. Bourhrara, A. Calloo, A. Gammicchia, F. Goubioud, L. Mao, F. Madiot, F. Malouch, E. Masiello, F. Moreau, S. Santandrea, D. Sciannandrone, I. Zmijarevic, E. Y. Garcia-Cervantes, G. Valocchi, J. F. Vidal, F. Damian, P. Laurent, A. Willien, A. Brighenti, L. Graziano, and B. Vezzoni. APOLLO3®: Overview of the New Code Capabilities for Reactor Physics Analysis. Nuclear Science and Engineering, 2024.

[5] O. Palii, M. Schlottbom, On a convergent DSA preconditioned source iteration for a DGFEM method for radiative transfer, Computers & Mathematics with Applications, Volume 79, Issue 12, 2020.

[6] Y. Saad. Numerical methods for large eigenvalue problems: revised edition. Society for Industrial and Applied Mathematics, 2011.

[7] J. Willert, H. Park, and D. A. Knoll. A comparison of acceleration methods for solving the neutron transport k-eigenvalue problem. Journal of Computational Physics, 2014, vol. 274, p. 681-694.

Study of radiative decay of the nucleus using a technic like Oslo-method

Radiative neutron capture is a nuclear reaction forming a compound nucleus which decays by emitting gamma-rays at excitation energy around the neutron binding energy. This well-known reaction which we known how to accurately measure its cross section at low incident neutron energies for most stable and few unstable nuclei close the stability valley, remains difficult to measure for exotic nuclei like fission fragments. Nuclear reaction models based essentially on stable nuclei, also struggle to provide reliable predictions of cross sections for these exotic nuclei. However, in the recent years, progress made related to the models and the measurements for the radiative capture show that significant improvements in including microscopic ingredients studies. These micoscopic ingredients: gamma strength function and nuclear level density, remain accesible to the experiment. These ingredients which respectively manage the way of how the gamma cascade occurs and the nuclear structure at high excitation energy can also be measured and calculated to be compared and suggest ways to improve the predictability of models. This kind of improvements have a direct impact for instance on the cross sections for these exotic nuclei which are produced in the stellar nucleosynthesis. The subject of thie thesis is to measure these quantities for a nucleus involved in the nucleosythesis using a new setup called SFyNCS.

Elliptic Flow of Charmed Hadrons in Heavy-Ion Collisions at LHCb?

The FLOALESCENCE project explores one of the most fundamental questions in Quantum Chromodynamics (QCD): how quarks and gluons transition from a deconfined Quark–Gluon Plasma (QGP) into ordinary hadrons.?This transition, called hadronization, occurred microseconds after the Big Bang and can be recreated today in ultra-relativistic lead–lead collisions at CERN’s Large Hadron Collider (LHC).
The PhD will focus on charm quarks—excellent probes of the QGP because they are produced early in the collision and interact throughout its evolution. Using the LHCb detector, uniquely sensitive in the forward rapidity region, the project aims to measure the elliptic flow (v2) of charmed baryons (?c+) and mesons (D0) in Pb–Pb collisions.?The goal is to test whether these heavy quarks thermalize and hadronize through a coalescence mechanism, a key feature of QGP dynamics.

Objectives and tasks:
- Extract and analyze ?c+ and D0 signals in newly collected 2024–2025 Pb–Pb datasets at LHCb.
- Implement a novel flow analysis method (based on the reformulated Lee–Yang Zeros approach) for the first time at LHCb.
- Develop an event-by-event multiplicity metric to correlate flow with system energy density.
- Compare results to theoretical models and cross-check with measurements at central rapidity (ALICE).
- Publish results and present findings at international conferences.

The successful candidate will:
- Develop advanced data-analysis expertise with CERN’s LHCb software framework, ROOT, and machine learning–based signal extraction.
- Gain in-depth knowledge of QCD and relativistic heavy-ion physics, especially QGP properties and collective phenomena.
- Learn modern statistical methods for flow analysis and uncertainty estimation.
- Acquire collaborative and communication skills within a major international experiment (LHCb), including presentations in collaboration meetings and conferences.
- Build strong experience in scientific computing, big-data handling, and detector physics, valuable for both academic and industry careers.

Designing a hybrid CPU-GPU estimator for neutron transport: Advancing eco-efficient Monte Carlo simulations

Digital twins incorporating Monte Carlo simulation models are currently being developed for the design, operation, and decommissioning of nuclear facilities. These twins are capable of predicting physical quantities such as particle fluxes, gamma/neutron heating, and dose equivalent rates. However, the Monte Carlo method presents a major drawback: high computational time to achieve acceptable variance levels.
To enhance simulation efficiency, the eTLE estimator has been developed and integrated into the TRIPOLI-4® Monte Carlo code. Compared to the conventional TLE (Track Length Estimator), eTLE offers lower theoretical variance, particularly in highly absorbing media, by contributing to the detector response even when particles do not physically reach it. Nevertheless, its computational cost remains significant, especially when evaluating multiple detectors.
Two recent PhD works have proposed variants to overcome this limitation. The Forced Detection eTLE- (Guadagni, EPJ Plus 2021) employs preferential sampling that directs pseudo-particles toward the detector at each collision. It is particularly effective for small detectors and configurations with moderate shielding, especially for fast neutrons. The Split Exponential TLE (Hutinet & Antonsanti, EPJ Web 2024) is based on an asynchronous GPU approach, offloading straight-line particle transport to the graphics processor. Through multiple sampling, it maximizes GPU utilization and enables more efficient exploration of phase space.
The proposed thesis aims to combine these two approaches into a hybrid estimator named seTLE-DF. This new estimator could be used either directly or to generate importance maps without relying on auxiliary deterministic calculations. Its implementation will require dedicated GPU developments, particularly to optimize the geometry library and memory management in complex geometries.
This research topic aligns with green computing objectives, aiming to reduce the carbon footprint of high-performance computing. It relies on a hybrid CPU-GPU strategy, avoiding full porting of the Monte Carlo code to GPU. Solutions such as half-precision formats will be considered, and an energy impact assessment will be conducted before and after implementation. The future PhD student will be welcomed with the IRESNE Institute (CEA Cadarache)and will acquire strong expertise in neutron transport simulation, facilitating integration into major research institutions or companies within the nuclear sector.

Probing quantum information with the top quark at the LHC

This PhD project aims to explore the quantum nature of top-quark pair production at the Large Hadron Collider by studying spin correlations and entanglement-related observables in data recorded by the ATLAS experiment. The recent breakthrough observations of entanglement in top-antitop events have opened an entirely new window onto the quantum structure of fundamental interactions, transforming the LHC into a machine to test quantum information at the TeV scale. Building on this momentum, the thesis will focus on reconstructing the quantum state of top-quark pairs using ATLAS Run-3 data, with particular attention to the extraction of spin correlations and entanglement-sensitive observables in challenging high-momentum topologies. By improving reconstruction strategies and carefully assessing detector effects, the aim is to measure quantum properties with good precision and to contribute to understand what quantum information can bring us to our understanding of elementary particles.

SEARCH FOR DIFFUSE EMISSIONS AND SEARCHES IN VERY-HIGH-ENERGY GAMMA RAYS AND FUNDAMENTAL PHYSICS WITH H.E.S.S. AND CTAO

Observations in very-high-energy (VHE, E>100 GeV) gamma rays are crucial for understanding the most violent non-thermal phenomena at work in the Universe. The central region of the Milky Way is a complex region active in VHE gamma rays. Among the VHE gamma sources are the supermassive black hole Sagittarius A* at the heart of the Galaxy, supernova remnants and even star formation regions. The Galactic Center (GC) houses a cosmic ray accelerator up to energies of PeV, diffuse emissions from GeV to TeV including the “Galactic Center Excess” (GCE) whose origin is still unknown, potential variable sources at TeV, as well as possible populations of sources not yet resolved (millisecond pulsars, intermediate mass black holes). The GC should be the brightest source of annihilations of massive dark matter particles of the WIMPs type. Lighter dark matter candidates, axion-like particles (ALP), could convert into photons, and vice versa, in magnetic fields leaving an oscillation imprint in the gamma-ray spectra of active galactic nuclei (AGN).
The H.E.S.S. observatory located in Namibia is composed of five atmospheric Cherenkov effect imaging telescopes. It is designed to detect gamma rays from a few tens of GeV to several tens of TeV. The Galactic Center region is observed by H.E.S.S. for twenty years. These observations made it possible to detect the first Galactic Pevatron and place the strongest constraints to date on the annihilation cross section of dark matter particles in the TeV mass range. The future CTA observatory will be deployed on two sites, one in La Palma and the other one in Chile. The latter composed of more than 50 telescopes will provide an unprecedented scan of the region of the Galactic Center.
The proposed work will focus on the analysis and interpretation of H.E.S.S observations carried out in the Galactic Center region for the search for diffuse emissions (populations of unresolved sources, massive dark matter) as well as observations carried out towards a selection of active galactic nuclei for the search for ALPs constituting dark matter. These new analysis frameworks will be implemented for the CTA data analyses. An involvement in the commissioning of the first MSTs in Chile and in the data analysis for early science are expected.

Impact of magnetohydrodynamic on access and dynamics of X-point radiator regimes (XPR)

ITER and future fusion powerplants will need to operate without degrading too much the plasma facing components (PFC) in the divertor, the peripheral element with is dedicated to heat and particle exhaust in tokamaks. In this context, two key factors must be considered: heat fluxes must stay below engineering limits both in stationary conditions and during violent transient events. An operational regime recently developed can satisfy those two constraints: the X-point Radiator (XPR). Experiments on many tokamaks, in particular WEST which has the record plasma duration in this regime (> 40 seconds), have shown that it allowed to drastically reduce heat fluxes on PFCs by converting most of the plasma energy into photons and neutral particles, and that it also was able to mitigate – or even suppress – deleterious magnetohydrodynamic (MHD) edge instabilities known as ELMs (edge localised modes). The mechanisms governing these mitigation and suppression are still poorly understood. Additionally, the XPR itself can become unstable and trigger a disruption, i.e., a sudden loss of plasma confinement cause by global MHD instabilities.
The objectives of this PhD are: (i) understand the physics at play during the interaction XPR-ELMs, and (ii) optimise the access and stability of the XPR regime. To do so, the student will use the 3D non linear MHD code JOREK, the European reference code in the field. The goal is to define the operational limits of a stable XPR with small or no ELMs, and identify the main actuators (quantity and species of injected impurities, plasma geometry).
A participation to experimental campaigns of the WEST tokamak (operated by IRFM at CEA Cadarache) – and of the MAST-U tokamak operated by UKAEA – is also envisaged to confront numerical results and predictions to experimental measurements.