Spectro-temporal analysis of Gamma-Ray Burst afterglows detected with SVOM

Gamma-Ray Bursts (GRB) are the most powerful explosions in the Universe. They last a few tens of seconds and emit the same amount of energy as the Sun during its entire lifetime. They gamma-ray emission is followed by a long lasting (hours to days) emission from the X-rays to the radio band. This "afterglow" emission is rich on information about the GRB nearby environnent and host galaxy. SVOM (Space based astronomical Variable Object Monitor) is a Sino-French mission, dedicated to GRB studies, and has been successfully launched in June 2024. It carries a multi-wavelength payload covering gamma-rays/X-rays/optical and includes two dedicated ground based robotic telescopes in Mexico and China.
The PHD project is focussed on the exploitation of the SVOM data for GRBs. The successful candidate will join the MXT science Teal at DAp. MXT is a new type of X-ray telescope, for which the DAp is responsible and its Instrument Centre is also hosted at DAp.
The PHD student will participate actively to the spectral and temporal analysis of MXT data. These data will be compared
to the other data acquired by the SVOM collaboration, especially in the optical an infrared domains.
This dataset will be used as a support to the physical interpretation of GRBs. More specifically, the aspects related to the modeling of the energy injection in the first phases of the afterglow will be used to determine the nature of the compact object at the origin of the relativistic flux, generating the electromagnetic emission observed.

Stochastic Neutron Noise Estimation Using a Rare-Event Simulation Approach. Application to the Monitoring of Nuclear System Reactivity

This PhD project aims to develop an innovative method to characterize the reactivity of fissile systems by analyzing their stochastic fluctuations, known as zero-power neutron noise. In a subcritical fissile medium, neutrons originating from spontaneous fission can initiate short and random chain reactions, generating a fluctuating signal. This noise carries essential information on the distance of the system to criticality, a key parameter both for the safety of nuclear installations (prevention of criticality accidents) and for the detection of undeclared fissile materials (nuclear security and non-proliferation).

Existing theoretical approaches to infer system reactivity from neutron noise are limited to idealized situations and become unsuitable in realistic configurations, particularly when the system is strongly subcritical or when significant uncertainties exist regarding its geometry or composition (as in the case of the Fukushima Daiichi corium or spent fuel storage). Monte Carlo simulations then appear as a natural alternative, but current simulations rely on variance reduction techniques that fail to correctly preserve stochastic fluctuations.

This thesis proposes to address this scientific challenge by adapting a relatively recent variance reduction method known as Adaptive Multilevel Splitting (AMS), originally developed to efficiently sample rare events while preserving their statistical properties. The goal is to extend this method to neutron transport in multiplying media and to make it a tool capable of faithfully simulating the temporal correlations characteristic of neutron noise. Following the theoretical developments, the algorithm will be implemented in Geant4, compared to analytical benchmark solutions, and experimentally validated through in situ measurements (using neutron sources or research reactors). In the long term, this work may lead to direct applications in nuclear monitoring, safety diagnostics, and detector physics, while also opening perspectives in fundamental physics and medical physics.

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.

From Cosmic Web to Galaxies: Tracing Gas Accretion at High Redshift through Observations and Simulations

This thesis aims to develop an integrated understanding of high-redshift galaxies within their large-scale structures. We will investigate how feedback and nuclear activity from these galaxies affect their environments by coupling observational data with cosmological simulations.
Our primary objectives are to:
1. Advance the diagnostic capabilities for studying diffuse gas.
2. Test and validate current paradigms of gas accretion.
Our observational work will utilize new data from Keck and the Very Large Telescope on Lyman-alpha halos around massive groups and clusters at z>2, which are already largely in hand. We will also incorporate a growing body of data from the James Webb Space Telescope (JWST) on the same targets to reveal the properties of galaxies and their active galactic nuclei (AGNs).
On the theoretical side, we will use publicly available results from the TNG100, HORIZON5, and CALIBRE simulations to understand galaxy evolution, learning from both the successes and failures in the comparison with observations. Ultimately, this will allow us to inform new, high-fidelity simulations of the circum-galactic medium, designed specifically to constrain gas accretion processes.
This research directly supports our long-term goal of preparing for the exploitation of BlueMUSE, a new instrument being built for the VLT, in which we participate. It will also address one of the key open questions in astrophysics, as highlighted by the Astro2020 Decadal Survey.

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.

Characterization and calibration of cryogenic detectors at the 100 eV scale for the detection of coherent neutrino scattering (CEvNS)

DESCRIPTIONS:

The NUCLEUS experiment [1] aims to detect reactor neutrinos via coherent elastic neutrino–nucleus scattering (CEvNS). Predicted in 1974 and first observed in 2017, this process provides a unique opportunity to test the Standard Model at low energies. Because the scattering is coherent over the entire nucleus, the cross section is enhanced by several orders of magnitude, making CEvNS also promising for reactor monitoring using neutrinos.

The NUCLEUS experimental setup is currently being installed near the EDF nuclear reactors in Chooz (Ardennes, France), which constitute an intense neutrino source. The only physical signal of a CEvNS event is the tiny recoil of the target nucleus, with an energy below 1 keV. To detect this, NUCLEUS uses CaWO4 crystals of about 1 g, placed in a cryostat cooled to 15 mK. The nuclear recoil produces vibrations in the crystal lattice, equivalent to a temperature rise of about 100 µK, measured with a Transition Edge Sensor (TES) deposited on the crystal. These detectors achieve excellent energy resolutions of only a few eV and detection thresholds on the order of ~10 eV [2]. The NUCLEUS setup was successfully tested and validated in 2024 at TU Munich [3], and data taking at Chooz is scheduled to start in summer 2026, simultaneously with the beginning of the PhD. An initial contribution will involve data acquisition and analysis at the reactor site. More specifically, the PhD student will be responsible for the characterization of the deployed cryogenic CaWO4 detectors — stability, energy resolution, calibration, and intrinsic background of the crystal.

Calibration at the sub-keV scale is a crucial challenge for CEvNS (and dark matter) experiments. Until recently, it was extremely difficult to generate nuclear recoils of known energy to characterize detector responses. The CRAB method [4, 5] addresses this issue by using thermal neutron capture (25 meV) on nuclei that constitute the cryogenic detector. The resulting compound nucleus has a well-known excitation energy — the neutron separation energy — between 5 and 8 MeV, depending on the isotope. When it de-excites by emitting a single gamma photon, the nucleus recoils with a precisely determined energy given by two-body kinematics. A calibration peak in the desired energy range of a few hundred eV then appears in the detector’s energy spectrum. A first measurement in 2022, using a NUCLEUS CaWO4 detector and a commercial ²5²Cf neutron source, validated this method [6].

The second part of the PhD will take place within the “high-precision” phase of the project, which consists in performing measurements with a pure thermal neutron beam from the TRIGA-Mark-II reactor in Vienna (TU Wien, Austria). The calibration setup was successfully installed and characterized in 2025 [7]. It consists of a cryostat housing the cryogenic detectors to be characterized, surrounded by large BaF2 crystals for coincidence detection of the nuclear recoil and the gamma ray that induced it. The whole setup is placed directly on the neutron beam axis, which provides a flux of about 450 n/cm²/s. This coincidence technique will significantly reduce background and extend the CRAB method to a wider energy range and to materials used in most cryogenic detectors. These measurements are expected to provide a unique characterization of the response of cryogenic detectors in the energy region of interest for light dark matter searches and coherent neutrino scattering. In parallel with the measurement of nuclear recoils, the installation of a low-energy X-ray source in the cryostat will generate electronic recoils, enabling a direct comparison between the detector responses to sub-keV energy deposits produced by nuclear and electronic recoils.

The arrival of the PhD student will coincide with the completion of the measurement program on CaWO4 and Al2O3 detectors of NUCLEUS and with the start of the measurement programs on Ge (TESSERACT project) and Si (BULLKID project) detectors.
The high-precision measurements will also open a new sensitivity window to subtle effects coupling nuclear physics(nuclear de-excitation times) and solid-state physics (nuclear recoil times in matter, and the creation of crystal defects induced by nuclear recoils) [8].

The PhD student will be deeply involved in all aspects of the experiment: simulation, data analysis, and interpretation of the obtained results.

WORK PLAN:

The PhD student will actively participate in data taking and in the analysis of the first results from the NUCLEUS cryogenic CaWO4 detectors at Chooz. This work will be carried out in collaboration with the Nuclear Physics Department (DPhN), the Particle Physics Department (DPhP) of CEA-Saclay, and the TU Munich team. It will begin with familiarization with the CAIT analysis framework used for cryogenic detectors. The student will focus in particular on detector calibration, studying the detector response to electronic recoils induced by optical photon pulses injected through fibers and by X-ray fluorescence generated by cosmic rays. Once this calibration is established, two types of backgrounds will be investigated: Nuclear recoils in the keV range induced by cosmogenic fast neutrons, and a low-energy background, known as the Low Energy Excess (LEE), intrinsic to the detector.
The comparison between the experimental and simulated fast neutron background spectra will be analyzed in light of the differences between nuclear and electronic recoil responses measured in the CRAB project. The long data-taking periods at the Chooz site will also be used to study the time evolution of the LEE background. This work will be conducted in collaboration with solid-state physics experts from the Institute for Applied Sciences and Simulation (CEA/ISAS) to better understand the origin of the LEE, which remains a major open question in the cryogenic detector community.
The analysis skills acquired on NUCLEUS will then be applied to the high-precision CRAB measurement campaigns planned for 2027 at the TRIGA reactor (TU Wien) with Ge and Si detectors. The student will be deeply involved in the setup, data acquisition, and analysis of results. The planned measurements on germanium, using both phonon and ionization channels, have the potential to resolve the current ambiguity in the ionization yield of low-energy nuclear recoils, a key factor for the sensitivity of future experiments.
The high calibration precision will also be exploited to study fine effects in nuclear and solid-state physics, such as timing effects and crystal defect formation induced by nuclear recoils in the detector. This study will be conducted in synergy with teams from CEA/IRESNE and CEA/ISAS, who provide detailed simulations of nuclear de-excitation gamma cascades and molecular dynamics simulations of nuclear recoil propagation in matter.

Through this work, the student will receive comprehensive training as an experimental physicist, including strong components in simulation and data analysis, as well as hands-on experience with cryogenic techniques during the commissioning of the NUCLEUS and CRAB detectors. The proposed contributions are expected to lead to several publications during the PhD, with high visibility in the CEvNS and dark matter communities. Within the CEA, the student will also benefit from the exceptionally cross-disciplinary nature of this project, which already
fosters regular interaction among the communities of nuclear physics, particle physics and condensed matter physics.

COLLABORATIONS:

NUCLEUS: Germany (TU-Munich, MPP), Austria (HEPHY, TU-Wien), Italy (INFN), France (CEA-Saclay).
CRAB: Germany (TU-Munich, MPP), Austria (HEPHY, TU-Wien), Italy (INFN), France (CEA-Saclay, CNRS-IJCLab, CNRS-IP2I, CNRS-LPSC).

BIBLIOGRAPHY:

[1] NUCLEUS Collaboration, Exploring CE?NS with NUCLEUS at the Chooz nuclear power plant, The European Physical Journal C 79 (2019) 1018.
15, 48, 160, 174
[2] R. Strauss et al., Gram-scale cryogenic calorimeters for rare-event searches, Phys. Rev. D 96 (2017) 022009. 16, 18, 78, 174
[3] H. Abele et al., Particle background characterization and prediction for the NUCLEUS reactor CE?NS experiment, https://arxiv.org/abs/2509.03559
[4] L. Thulliez, D. Lhuillier et al. Calibration of nuclear recoils at the 100 eV scale using neutron capture, JINST 16 (2021) 07, P07032
(https://arxiv.org/abs/2011.13803)
[5]https://irfu.cea.fr/dphp/Phocea/Vie_des_labos/Ast/ast.php?id_ast=4970
[6] H. Abele et al., Observation of a nuclear recoil peak at the 100 eV scale induced by neutron capture, Phys. Rev. Lett. 130, 211802 (2023) (https://arxiv.org/abs/2211.03631)
[7] H.Abele et al., The CRAB facility at the TUWien TRIGA reactor: status and related physics program, (https://arxiv.org/abs/2505.15227)
[8] G. Soum-Sidikov et al., Study of collision and ?-cascade times following neutron-capture processes in cryogenic detectors Phys. Rev. D
108, 072009 (2023) (https://arxiv.org/abs/2305.10139)

Explainable observers and interpretable AI for superconducting accelerators and radioactive isotope identification

GANIL’s SPIRAL1 and SPIRAL2 facilities produce complex data that remain hard to interpret. SPIRAL2 faces instabilities in its superconducting cavities, while SPIRAL1 requires reliable isotope identification under noisy conditions.
This PhD will develop observer-based interpretable AI, combining physics models and machine learning to detect, explain, and predict anomalies. By embedding causal reasoning and explainability tools such as SHAP and LIME, it aims to improve the reliability and transparency of accelerator operations.

Precision measurements of neutrino oscillations and search for CP violation with the T2K and Hyper-Kamiokande experiments

The study of neutrino oscillations has entered a precision era, driven by long-baseline experiments like T2K, which compare neutrino signals at near and far detectors to probe key parameters, including possible Charge-Parity Violation (CPV). Detecting CPV in neutrinos could help explain the Universe’s matter–antimatter asymmetry. T2K’s 2020 results gave first hints of CPV but remain limited by statistics. To improve sensitivity, T2K has undergone major upgrades: replacing the most upstream part of its near detector with a new target, increased accelerator power (up to 800 kW by 2025, aiming for 1.3 MW by 2030). The next-generation Hyper-Kamiokande (Hyper-K) experiment, starting in 2028, will reuse the T2K beam and near detector but with new far detector 8.4 times larger than Super-Kamiokande greatly boosting the statistics. The IRFU group has key role in the near detector upgrade and is now focusing on analysis, crucial for controlling systematic uncertainties crucial for the Hyper-K high statistics time. The proposed PhD work centers on analyzing the new near detector data: designing new sample selections taking into account for the low-momentum protons and neutrons from neutrinos, and refining neutrino–nucleus interaction models to improve energy reconstruction. The second goal is to propagate these improvements to Hyper-K, guiding future oscillation analyses. The student will also contribute to Hyper-K construction and calibration (electronics testing at CERN, installation in Japan).

Axion searches in the SuperDAWA experiment with superconducting magnets and microwave radiometry

Axions are hypothetical particles that could both explain a fundamental problem in strong interactions (the conservation of CP symmetry in QCD) and account for a significant fraction of dark matter. Their direct detection is therefore a key challenge in both particle physics and cosmology.

The SuperDAWA experiment, currently under construction at CEA Saclay, uses superconducting magnets and a microwave radiometer placed inside a cryogenic cryostat. This setup aims to convert potential axions into measurable radio waves, with frequencies directly linked to the axion mass.

The proposed PhD will combine numerical modeling with hands-on experimental work. The student will develop a detailed model of the experiment, including magnetic fields, radio signal propagation, and detector electronics, validated step by step with real measurements. Once the experiment is running, the PhD candidate will participate in data-taking campaigns and their analysis.

This project provides a unique opportunity to contribute to a state-of-the-art experiment in experimental physics, with direct implications for the global search for dark matter.

Testing the Standard Model in the Higgs-top sector in a new inclusive way with multiple leptons using the ATLAS detector at the LHC

The LHC collides protons at 13.6 TeV, producing a massive dataset to study rare processes and search for new physics. The production of a Higgs boson in association with a single top quark (tH) in the multi-lepton final state (2 same-sign leptons or 3 charged leptons) is particularly promising, but challenging to analyze due to undetected neutrinos and fake leptons. The tH process is especially interesting because its small Standard Model cross section originates from a subtle destructive interference between diagrams including the Higgs coupling to the W boson and the Higgs coupling to the top quark. This makes tH uniquely sensitive: even small deviations from the Standard Model can strongly enhance its production rate. The measurement of the tH cross section is delicate because the ttH and ttW processes have similar topologies and much larger cross sections, requiring a simultaneous extraction to obtain a reliable result and properly account for correlations between signals. ATLAS observed a moderate excess of tH using the Run 2 dataset (2.8 s), making the analysis of Run 3 data including these correlations crucial. The thesis will first exploit AI algorithms based on Transformer architectures to reconstruct event kinematics and extract observables sensitive to the CP nature of the Higgs-top coupling. In a second phase, a global approach will be adopted to analyze simultaneously the ttW, ttZ, ttH, tH, and 4-top processes, searching for anomalous couplings, including those violating CP symmetry, within the framework of the Standard Model Effective Field Theory (SMEFT). This study will provide the first complete measurement of tH in the multi-lepton channel with Run 3 data and will pave the way for a global analysis of rare processes and anomalous couplings at the LHC in this channel.