PLANETOLOGY IN A LABORATORY
A. Benuzzi-Mounaix LULI Laboratory - Ecole Polytechnique OUTLINE
ü Introduction: what is planetology in laboratory?
ü How to do planetology in a laboratory? Generation of “planetary” matter, Equation of state measurements, diagnostics .....
ü Some interesting results
ü The research program PLANETLAB (2012-2016)
ü Facilities for planetology in a laboratory, teams involved, PhD subjects Introduc on PLANETARY SCIENCE
² Study of Forma on and evolu on of planets
² Study of planet interiors
FAST GROWING SCIENCE DUE TO EXOPLANETS DISCOVERY
From 1988 to today
838 planets discovered Introduc on Key ques ons about planets inside and outside the Solar System
•What is the nature of the iron core at the center of Earth and other terrestrial planets?
•What is the interior structure of Jupiter and the other giant planets?
•What kinds of planets exist outside our solar system? Can we characterize their structures, composi ons, dynamics and evolu on?
Why the Saturn luminosity is not compa ble with its age?
To be able to answer
We have to know the equa on of state (i.e. the rela on between ρ,P,E) of the ma er exis ng in the planet interiors
But …. Introduc on What type of ma er?
Our giant planets
The important classes of materials are:
-Hydrogen and Helium
-Simple molecular
compounds H2O, NH3, CH4
-Silicates and oxides:
(Mg,Fe)2SiO4, (Mg,Fe)SiO3, SiO2, (Mg,Fe)O
-Iron and iron alloys
Introduc on At which condi ons? Introduc on This ma er is called Warm Dense Ma er
Al phase diagram (ρ,T) 104 classical plasma ! = 1 0. 1 ρsolide < ρ < 100 ρsolide 103 dense ! = 10 0. 1 eV < T < 100 eV plasma 102 WDM ! = 100 Condi ons found not only 101 high density in planetology but also in matter Iner al Confinement Fusion 100 10-4 10-2 1 102 104 Density ( g/cm3) Introduc on This ma er is called Warm Dense Ma er
Al phase diagram (ρ,T) 104 classical plasma ! = 1
103 dense ! = 10 plasma 102 WDM is the state at the intersec on WDM ! = 100 between plasma physics and condensed 101 high ma er physics. density matter WDM is correlated and degenerate 100 Correlated Γ = E / E > 1 -4 c th 10 10-2 1 102 104 Coulomb energy E Density ( g/cm3) c Thermal energy Eth Degenerate ΤF > T Perfect gas does not apply Perturba on theory is completely invalid – no small parameters Introduc on Important progresses on theore cal side: Quantum Molecular Dynamics calcula ons
ü Electrons are treated by quantum mechanics through Schrödinger equa on within the so called "Kohn-Sham" theory which is exact (ab ini o)
ü The ions mo on is deduced from Newton equa ons using quantum forces created by electrons.
ü Thermodynamical equilibrium
Electronic density of D ü Equa on of state, transport proper es, opaci es…. 2 @ 1g/cc and 29.000 K No use of empirical parameters No need of an a priori experimental knowledge
But to confirm their predic ve power Need of experimental data How do planetology in a laboratory How create this planetary ma er in a laboratory?
Sta c way
Diamond cell Isothermal Compression P ≈ 0 - a few Mbar Dynamic way
Chemical explosions, gas guns
(Nuclear explosions)
High power lasers Shock Compression P ≈ 0 - hundreds of Mbar How to do planetology in a laboratory What is a shock ?
An abrupt change in pressure, density, energy which moves through a medium
Pressure, Density, Energy Rankine Hugoniot relations
Shock velocity D Conservation of the Fluid velocity U mass
z momentum energy How to do planetology in a laboratory Laser driven shock
Cold target Hot Laser plasma
Shock
Ablation pressure 14 2 IL en 10 W/cm 2 P " 12(I /#) 3 λ en µm L P en Mbar
Today Pressures of several tens of Mbar
! How to do planetology in a laboratory EOS measurements using a shock
Rankine-Hugoniot
ρ0D = ρ (D-U) mass
ρ0DU = P-P0 momentum
2 ρ0D(E-E0+U /2)=PU energy
3 equations, 5 unknown parameters (D,U,P,E,ρ)
To have a point on f(P, E, ρ) = 0
ü 2 parameters - absolute measurement 2 velocities D and U → time and distance How to do planetology in a laboratory How do we measure shock parameters?
Usual diagnos cs
Velocity Interferometer System for Any Reflector (VISAR)
target! streak Velocities (D or U) by Doppler effect & reflectivity λ
λ0 τ D or U
Self-emission D & Temperature T target
streak Kirchoff law: I(T,λ) = α(λ) IP(T, λ) Planckian α (λ) = 1-R (λ) radiation T How to do planetology in a laboratory Recent developments to have microscopic informa on X-ray sca ering diagnos c
X-ray sca ering cross 3 contribu ons: sec on
Free e- Weakly bound e- Thightly bound e-
Inelas c features broadened by Elas c feature thermal mo ons (Rayleigh peak)
Intensity Inelastic feature Inelas c feature Elastic feature Electronic density and temperature
Elas c feature
Correla on effects, phase transi ons ECompton E0 Energy How to do planetology in a laboratory Recent developments to have microscopic informa on X-ray absorp on near edge spectroscopy
!
3s - 3p band ! F 2p 2s Absorp on X K-edge
• Probe the valence 1s electrons r • Probe the local atomic order
Powerful diagnos c to study structure changes, phase transi ons and to test approxima ons used in theories
SOME INTERESTING RESULTS Some results Some results on water
The mo va on
0.7 Metallic
0.6
Fluid 0.5 Neptune isentrope Ionic 0.4
r Main Hugoniot
0.3
c
é
l
Temperature T (eV) T Temperature Superionic Molecular Molecular
o 0.2
M
0.1 Precompressed Hugoniot Solid 0 0 0.5 1 1.5 2 2.5 3
Pressure P (Mbar) Phase diagram obtained from ab initio computations (C. Cavazzoni et al., Science 283, 1999)
Why is the magne c field of these planets so high and asymmetric?
Is it because the mantles are cons tuted of « ice » layers, fluid and conduc ng?
0.7 0.7 Metallic Metallic 0.6 0.6
Fluid
Fluid 0.5 0.5 Neptune isentrope Ionic Neptune isentrope Ionic 0.4
0.4 r
r
0.3
0.3 c
c
é
l
Temperature T (eV) T Temperature Superionic
Molecular é
l
Temperature T (eV) T Temperature o Superionic
0.2 Molecular
o 0.2
M
M
0.1 0.1 Solid
0 Solid 00 0.5 1 1.5 2 2.5 3 0 0.5 1 1.5 2 2.5 3 Pressure P (Mbar) Pressure P (Mbar) Phase diagram obtained from ab ini o computa ons (C. Cavazzoni et al., Science 283, 1999) Phase diagram obtained from ab ini o computa ons (C. Cavazzoni et al., Science 283, 1999) Some results Some results on water
0.7 Metallic
0.6
Fluid 0.5 Op cal proper es have Neptune isentrope Ionic been measured In the phase 0.4
r Main Hugoniot diagram
0.3
c
é
l
Temperature T (eV) T Temperature Superionic Molecular Molecular
o 0.2 M Hugoniot has been 0.1 Precompressed Hugoniot Solid measured 0 0 0.5 1 1.5 2 2.5 3
Pressure P (Mbar) Phase diagram obtained from ab initio computations (C. Cavazzoni et al., Science 283, 1999) And Hugoniot has been
P. Celliers et al PoP 2004
0.7 0.7 Metallic Metallic 0.6 0.6
Fluid
Fluid 0.5 0.5 Neptune isentrope Ionic Neptune isentrope Ionic 0.4
0.4 r
r
0.3
0.3 c
c
é
l
Temperature T (eV) T Temperature Superionic
Molecular é
l
Temperature T (eV) T Temperature o Superionic
0.2 Molecular
o 0.2
M
M
0.1 0.1 Solid
0 Solid 00 0.5 1 1.5 2 2.5 3 0 0.5 1 1.5 2 2.5 3 Pressure P (Mbar) Pressure P (Mbar) Phase diagram obtained from ab ini o computa ons (C. Cavazzoni et al., Science 283, 1999) Phase diagram obtained from ab ini o computa ons (C. Cavazzoni et al., Science 283, 1999) Some results Some results on iron
Mantle ! !
Outer core ! Liquid Fe!
Inner core! solid Fe! Simultaneous Temperature/ Pressure measurements across T at Earth’s Inner-Outer core boundary melting the melting region → Tmelting P=3.3 Mbar!
Critical constraint for modeling the chemical composition and energy balance of the Earthʼs core! Huser et al., PoP, 2005 Some results Study of electronic structure changes in a large WDM domain using XANES
Test WDM theories in a well-known case (Al)
Dense plasma Atomic state ! !
...
ρ 3s - 3p 3p band ! 3s F Pre-edge onset 2p 2p (1.6 g/cc) 2s 2s
density
1s 1s r r
The onset of a metal/non metal transi on has been put in evidence and confirmed by calcula ons.
A. Benuzzi-Mounaix et al PRL (2011); A. Levy et al. PRL (2012) PLANETLAB research program (ANR 2012-2016)
Coordinated by LUTH, Observatoire de Paris Funds 500.000 euros
ü Establish benchmarking values for the equa ons of Main goals: states, mel ng curves and the transport proper es in the Fe-Si-Mg-O-S complex system exis ng inside core of giant planets and super earths ü Make super earth models
THEORETICAL APPROACH EXPERIMENTAL APPROACH (ab ini o calcula ons) (using sta c and dynamic methods)
(LUTH Meudon, CEA Bruyères le (LULI Ecole Polytechnique, Chatel) IMPMC Paris 6) Planetlab research program (ANR 2012-2016)
Program based on calcula ons for exoplanets descrip on a er a careful valida on against both sta c and dynamical experiments that can be performed in a limited region of the phase diagram. Concerning silica Dynamic experiment to study dissocia on around the Hugoniot (scheduled on LULI 2000) using XANES/X-ray sca ering Ab ini o calcula ons up to 20Mbar and 50000 K: mel ng, EOS, …
Dynamic experiment using quasi-isentropic compression to detect mel ng (experiment scheduled on LIL laser at CEA)
Super-earth
Sta c experiments coupled to synchrotron radia on Laser facili es in France and around the world for planetology
High energy and long pulses
JANUS laser – LLNL - USA
PALS – Prague - Czech Republic VULCAN – RAL - UK Energy LULI2000 - Ecole Polytechnique - FRANCE GEKKO laser – Osaka University - JAPAN ORION - AWE - UK LIL - CEA - FRANCE OMEGA laser – University of Rochester - USA LMJ/PETAL laser – CEA – France NIF laser – LLNL - USA Sta c measurements and facili es
Contact: [email protected]
Located in Grenoble, This is the largest synchrotron in Europe
Located close to Orsay, this french facility is smallest + in situ X-ray than ESRF but sufficient for measurements many type of high pressure (diffrac on, experiments absorp on, …) Experimental Teams in Paris area
LULI laboratory : A. Benuzzi-Mounaix, M. Koenig, A. Ravasio, E. Brambrink, T. Vinci, A. Denoeud ( experiments with laser)
CEA (Bruyères le Chatel) : P. Loubeyre, S. Brygoo, G. Huser (experiments with DAC + laser)
IMPMC (Paris 6): G. Morard, G. Fiquet, F. Guyot (experiments with DAC)
Theore cal teams in Paris area
LUTH (Meudon) : S. Mazevet, R. Musella
CEA (Bruyères le Chatel) : V. Recoules, J. Bouchet Internship M2 (followed by PHD)
Experimental study for planetology using high power laser Supervisors : A. Ravasio Contact: [email protected]
LULI 2000: 1kJ at 1054 nm 0.5ns-3ns 100J at 1054 1-5 ps
Internship M2 + PHD subject (already funded)
Proper es of iron alloys for planetary modeling Supervisors : J. Bouchet (CEA) V Recoules (CEA) S. Mazevet (LUTH) Contact: [email protected]
Na onality requested: French Previous ab ini o results on iron
Cores are solid for planets having their
mass > 5 MT
If no liquid phase, no magne c field -> solar winds can reach the surface -> probably no possible life on these planets
Morard et al. HEDP, 2011 Internship M2 + PHD subject (already funded) Previous ab ini o results on iron
Proper es of iron alloys for planetary modeling Supervisors : J. Bouchet (CEA) V Recoules (CEA) S. Mazevet (LUTH) Contact: [email protected]
The discovery of extra solar planets (exoplanets) calls for a be er understanding of the proper es of iron and its alloys Fe-Si-S-O in temperature and pressure domains never studied un l now. In this thesis we will use the solid state simula on methods as DFT to study the alloys effects on iron proper es as the mel ng temperature and more generally on the equa on of state. Then these new data will be integrated in planetary models. This project will take place in an ANR obtained by the Observatoire de Morard et al. HEDP, 2011 Paris, the CEA, the Ecole Polytechnique et the university of Paris Jussieu and will be supervised by the Observatoire Cores are solid for planets having their
de Paris et the CEA. This subject has a strong numerical mass > 5 MT part and will call for a strong background in solid state physics, numerical simula ons, geophysics and planetary If no liquid phase, no magne c field -> modeling. solar winds can reach the surface -> probably no possible life on these Na onality requested: French planets THANK YOU FOR ATTENTION