Modeling of Heating and Radiative Properties of a Solid Target Irradiated by an Ultra-Intense Laser

Ultra-high-intensity (UHI) lasers can deliver petawatt (10^15 W) beams over ultra-short durations (from femtoseconds to picoseconds). The interaction of such a laser with a solid target allows the study of matter under extreme conditions, i.e., at solid densities brought to temperatures close to 10^7 K. These dense and hot plasmas have radiative properties (X-ray emission/absorption) that provide information about their state (temperature, etc.) and that can be measured experimentally. However, the physical mechanisms governing the heating and radiative properties are not yet fully understood.
The thesis aims to simulate the laser-target interaction using the CALDER code, a Particle-In-Cell (PIC) type code developed at the CEA. The PhD student will study the physics of isochoric heating with this code and seek to understand the predominant physical effects in these plasmas. He will model recent experiments aimed at studying these plasmas and conducted on academic laser facilities. Prospective simulations will help prepare future experiments, such as those planned on the PETAL laser at the CEA. Atomic physics codes may be used to refine the simulated X-ray spectra. The work will also include adapting CALDER to better model this physics. The results will contribute to a better understanding of the interaction between ultra-high-intensity laser-plasma and the radiative properties of irradiated solid targets.

Prompt fission neutron spectrum precision measurement in the spontaneous fission of 252Cf

The 252Cf(sf) prompt fission neutron spectrum (PFNS) is a reference neutron data that is widely used as a well known neutron flux for cross section measurements and neutron detector characterization. The current evaluation of the spectrum dates from the work of Mannhart in 1988. With the improvement of detection systems, the uncertainty on the spectrum has an increasingly significant impact on the uncertainty of new nuclear data measurements using it as a reference. Improving the precision on the 252Cf(sf) PFNS would therefore have a wide impact on the nuclear data community and improve the uncertainties on all data that were measured with respect to this reference spectrum. The thesis aims to measure again the 252Cf(sf) PFNS with a focus on the region below 1 MeV and the region above 8 MeV, where the uncertainties are greatest, using recoil proton detectors. The chosen candidate will have to actively participate to the design and construction of the setup, leading the technological choices through simulations, will participate to the experiment and do the data analysis. The work will then be presented in international conferences and peer-reviewed articles.

Neutron elastic and inelastic scattering measurement on 9Be with VENDETA

9Be plays a central role in fusion technology and in material testing reactor (MTR) as a neutron source moderator. However, existing nuclear data for neutron scattering on 9Be exhibit significant uncertainties, particularly in the 1.5-15 MeV energy range. In order to provide high quality data for improving nuclear reaction models and evaluated nuclear data libraries, an experiment to measure the elastic and inelastic scattering of neutrons on a 9Be target was proposed and accepted at the Neutrons For Science (NFS) facility. Neutrons will be measured using the recently developed VErsatile Neutron DETector Array (VENDETA), formed of high-resolution time-of-flight detectors which combine excellent neutron/gamma discrimination and efficiency down to 100 keV. The experiment will employ the quasi-monoenergetic neutron beams available at NFS from the p+7Li reaction, enabling a systematic investigation of scattering observables as a function of incident neutron energy.
The chosen candidate will lead the analysis to extract angular differential cross sections for both elastic and inelastic channels as a function of the incoming neutron energy. The data will have a direct impact on applications in nuclear energy and shielding design, while their interpretation in terms of partial decay width to the elastic and inelastic channels will improve our understanding of 9Be and 10Be nuclear structure.

Behavior of matter under isothermal dynamic compression: displacement of chemical reactivity; synthesis of new metastable materials; phase transition mechanisms.

The Diamond Anvil Cell equipped with piezoelectric actuators, or d-CED, is an innovative device that can generate dynamic compressions and decompressions over a wide range of pressure variation rates. The d-CED thus enables finely controlled dynamic stresses to be applied, with (de)compression rates that can vary over several orders of magnitude along isothermal paths. This paves the way for the creation of reference databases for the validation of microscopic mechanisms. Furthermore, the compression or decompression rates can be equated to ultra-fast heating or cooling rates of the sample, offering the possibility of exploring, in a highly controlled manner, certain phenomena still debated in the literature, such as the maximum stability of a solid beyond its melting point.
The objective of this thesis is to exploit the new possibilities offered by d-CED to demonstrate new phenomena or gain a detailed understanding of certain effects discussed in the literature, by performing ultra-fast temperature variations. A first application will consist of studying the nucleation kinetics of rare gases (Ar, Ne, Kr) as a function of the compression rate, and comparing them with recent measurements made at the XFEL in cryogenic jets. A second objective will be to study chemical changes, with an initial study focusing on the modification of the reactivity of nitromethane, a reference explosive. Another area of study will concern the synthesis of new molecular compounds from mixtures of dense molecular fluids (N2, H2, O2).

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.

Measurement of the speed of sound in H2 and He, key components of gas giant interiors

The goal of this thesis is to study hydrogen-helium mixtures in the fluid phase under high pressure and high temperature using Raman and Brillouin spectroscopy. The experiments will be conducted in a diamond anvil cell with laser heating, allowing exploration of a wide range of pressure and temperature conditions representative of the interiors of gas giant planets (1-300 GPa, 300-4000 K). Raman spectroscopy will be used to probe possible chemical changes occurring under extreme conditions, while Brillouin spectroscopy will provide access to the adiabatic sound velocity and the equations of state of these fluid mixtures. These data will be particularly useful for improving the modeling of Jupiter and Saturn’s interiors.

Hydrodynamic simulations of porous materials for ductile damage

The mechanical behavior of metallic materials under highly dynamical loading (schock) and especially their damage behavior is a topic of interest for the CEA-DAM. For tantalum, damage is ductile : by nucleation, growth and coalescence of voids within the material. Usual ductile damage models have been developed using the simplifying assumption that voids are isolated in the materials. However, recent studies by direct simulations explicitly describing a void population in the material (and experimental observations after failure) have shown the importance of void interaction for predicting ductile damage. Yet, the microscopical mechanisms of this interaction remain little known.
The objective of the PhD is to study the growth and coalescence phases of ductile damage through direct numerical simulations of a porous material undergoing dynamic loading. Hydrodynamic simulations, in which voids are explicitly meshed within a continuous matrix, will be used to study relevant scales of length and time. Monitoring the void population throughout the simulation will provide valuable information on the influence of void interaction during ductile damage. Firstly, the bulk behavior will be compared to the one predicted by usual models of isolated voids, showing the macroscopic effect of void interaction. Secondly, the evolution of the size distribution in the void population will be monitored. The last objective will be to understand microscopic void-to-void interaction. In order to take advantage of the wealth of simulation results, approaches based on artificial intelligence (neural networks on the graph associated with the pore population) will be used to learn the link between a void's neighborhood and its growth.
The doctoral student will have the opportunity to develop their skills in shock physics and mechanics, numerical simulations (with access to CEA-DAM supercomputers), and data science.

Multiscale modelling of twinning in tin

Twinning is a displacive deformation mechanism characterized by a continuous deformation of the material. Although widely studied for other industrial materials such as titanium alloys, this inelastic mechanism remains poorly understood and incompletely modeled for complex crystallographic structures. However, due to the reduced number of symmetries in these structures, dislocation slip is insufficient to accommodate deformation in certain loading directions, requiring the activation of twinning. This is the case for tin, which has a tetragonal structure. In particular, twinning contributes significantly to the mechanical response of tin at high strain rates and low temperatures. At intermediate temperatures and strain rates, a competition between dislocation plasticity and twinning plasticity can occur, making it crucial to describe the coupling between these two phenomena. Proposing a better description of this coupling will shed new light on the experimental data available at CEA DAM. The objective of the thesis is to develop a multiscale approach, from molecular dynamics to continuum mechanics, validated by experiments, to converge on a model that describes the behavior of tin over a wide range of temperatures and strain rates.

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.

microstructure informed kinetic model : application to solid explosives

When an explosive composition is subjected to an intense stress, such as a shock, the wave generated interacts with the microstructure and in particular with the defects it contains. Due to the nature of the defects, the energy can be localised, as when porosity is compacted, which can lead to the appearance of hot spots. Beyond a certain critical size, these hot spots grow as a result of the chemical decomposition of the explosive, and in some cases this can lead to the creation of a detonation wave. The role of these hot spots is therefore decisive in the initiation of solid explosives. The majority of macroscopic models used to study the shock-detonation transition (SDT) are phenomenological models calibrated on experiments (e.g. multi-strand gauge experiments) and therefore do not take into account the microstructural peculiarities specific to each explosive. It then becomes necessary to recalibrate a model for each composition, which limits any predictive capacity.
Microtomographic studies of real microstructures of explosive compositions have revealed that these deviate significantly from an average description based on a spherical pore. Through image segmentation, these microtomographs can provide essential ingredients for mesoscopic-scale simulation codes: these microstructures can be used directly as input for calculations or as a basis for generating virtual but realistic microstructures, thereby extending the accessible database given the experimental difficulties in generating this type of image in large numbers.
The computing power available today means that we can now envisage explicit simulations of realistic microstructures of explosive compositions. These simulations, in two or even three dimensions, will form the basis for the construction of a macroscopic kinetics model for modelling the shock-detonation transition. The results expected from this work are cross-disciplinary and can be transposed to all composite energetic materials. The effect of thermal or mechanical damage on the behaviour of an explosive or a solid propellant (vulnerability issues) could also benefit from this project. This more detailed knowledge of the role of microstructure (grain shape, porosity, etc.) could also improve filler manufacturing processes (e.g. ‘Very Insensitive’-RDX).

Top