Quantitative single shot and spatially resolved plasma wakefield - - PowerPoint PPT Presentation
Quantitative single shot and spatially resolved plasma wakefield - - PowerPoint PPT Presentation
Quantitative single shot and spatially resolved plasma wakefield diagnostics Muhammad Firmansyah Kasim University of Oxford, UK PhD Supervisors: Professor Peter Norreys & Professor Philip Burrows University of Oxford, UK JAI-Fest,
JAI-Fest, 10 November 2015 2 University of Oxford, RAL, UCL
Acknowledgement
University of Oxford:
- Peter Norreys
- Philip Burrows
- James Holloway
- Matthew Levy
- Naren Ratan
- Luke Ceurvorst
- James Sadler
AWAKE Collaboration, including:
- Allen Caldwell
- Matthew Wing
- Patric Muggli
- Edda Gschwendtner
Rutherford Appleton Laboratory:
- Robert Bingham
- Raoul Trines
- Rajeev Pattathil
- Dan Symes
- Pete Brummitt
- Chris Gregory
- Steve Hawkes
OSIRIS Consortium, including:
- Luis O. Silva
- Jorge Vieira
Indonesian Endowment Fund for Education
JAI-Fest, 10 November 2015 3 University of Oxford, RAL, UCL
Introduction
JAI-Fest, 10 November 2015 4 University of Oxford, RAL, UCL
Introduction
JAI-Fest, 10 November 2015 5 University of Oxford, RAL, UCL
Introduction
JAI-Fest, 10 November 2015 6 University of Oxford, RAL, UCL
Introduction
π π§, π¨ =
π
π π¬, π¨ ππ¬
JAI-Fest, 10 November 2015 7 University of Oxford, RAL, UCL
Introduction
π π§, π¨ =
π
π π¬, π¨ ππ¬
JAI-Fest, 10 November 2015 8 University of Oxford, RAL, UCL
Introduction
JAI-Fest, 10 November 2015 9 University of Oxford, RAL, UCL
Introduction
Abel transform
JAI-Fest, 10 November 2015 10 University of Oxford, RAL, UCL
Introduction
Abel transform
Abel transform: π π§, π¨ = 2
π§ β
π(π , π¨) π ππ π 2 β π§2
JAI-Fest, 10 November 2015 11 University of Oxford, RAL, UCL
Introduction
Abel transform
Abel transform: π π§, π¨ = 2
π§ β
π(π , π¨) π ππ π 2 β π§2 Abel inversion: π π , π¨ = β 1 π
π β ππ(π§, π¨)
ππ§ ππ§ π§2 β π 2
JAI-Fest, 10 November 2015 12 University of Oxford, RAL, UCL
Moving case
JAI-Fest, 10 November 2015 13 University of Oxford, RAL, UCL
Moving case
JAI-Fest, 10 November 2015 14 University of Oxford, RAL, UCL
Moving case
ππ©π πΎ = π/π
JAI-Fest, 10 November 2015 15 University of Oxford, RAL, UCL
Moving case
ππ©π πΎ = π/π
JAI-Fest, 10 November 2015 16 University of Oxford, RAL, UCL
Moving case
ππ©π πΎ = π/π Abel transform still works!
JAI-Fest, 10 November 2015 17 University of Oxford, RAL, UCL
Moving case
ππ©π πΎ = π/π Abel transform still works!
How if π > π?
JAI-Fest, 10 November 2015 18 University of Oxford, RAL, UCL
Moving case
The probe goes through different longitudinal position, π¨, during the interaction
JAI-Fest, 10 November 2015 19 University of Oxford, RAL, UCL
Moving case
The probe goes through different longitudinal position, π¨, during the interaction Normal Abel transform does not work.
JAI-Fest, 10 November 2015 20 University of Oxford, RAL, UCL
- Normal Abel transformation
π π, π = π
π β π π, π π
ππ β ππ ππ
Forward transform
JAI-Fest, 10 November 2015 21 University of Oxford, RAL, UCL
- Normal Abel transformation
π π, π = π
π β π π, π π
ππ β ππ ππ
- Modified Abel transformation
π π, π = π
π β
ππ©π ππ ππ β ππ π π, π π ππ β ππ ππ
Forward transform
JAI-Fest, 10 November 2015 22 University of Oxford, RAL, UCL
- Normal Abel transformation
π π, π = π
π β π π, π π
ππ β ππ ππ
- Modified Abel transformation
π π, π = π
π β
ππ©π ππ ππ β ππ π π, π π ππ β ππ ππ where π = ππ©π πΎ β π/π /ππ£π¨ πΎ,
Forward transform
JAI-Fest, 10 November 2015 23 University of Oxford, RAL, UCL
- Normal Abel transformation
π π, π = π
π β π π, π π
ππ β ππ ππ
- Modified Abel transformation
π π, π = π
π β
ππ©π ππ ππ β ππ π π, π π ππ β ππ ππ where π = ππ©π πΎ β π/π /ππ£π¨ πΎ, π(π, π) and π(π, π) are Fourier Transform of π π, π and π(π, π) in longitudinal direction, respectively
Forward transform
JAI-Fest, 10 November 2015 24 University of Oxford, RAL, UCL
- Normal Abel inversion
π π, π = β π π
π β ππ π, π
ππ ππ ππ β ππ
Inverse transform
JAI-Fest, 10 November 2015 25 University of Oxford, RAL, UCL
- Normal Abel inversion
π π, π = β π π
π β ππ π, π
ππ ππ ππ β ππ
- Modified Abel inversion
π π, π = β π π
π β
ππ©ππ’ ππ ππ β ππ π π π, π ππ ππ ππ β ππ
Inverse transform
JAI-Fest, 10 November 2015 26 University of Oxford, RAL, UCL
- Objective: diagnose electron density profile, π(π, πΌ), in the
wakefield
- Can be done by sending the laser probe with oblique angle
- f incidence relative to the wakefield
Application on plasma accelerators
JAI-Fest, 10 November 2015 27 University of Oxford, RAL, UCL
- Objective: diagnose electron density profile, π(π, πΌ), in the
wakefield
- Can be done by sending the laser probe with oblique angle
- f incidence relative to the wakefield
- What can we detect?
Application on plasma accelerators
JAI-Fest, 10 November 2015 28 University of Oxford, RAL, UCL
- By photon ray theory in plasma wake (Wilks,
1989) π¬π ππ β β ππ
π
πππ
π
π ππ ππ ππΌ ππ
Theory of photon acceleration
JAI-Fest, 10 November 2015 29 University of Oxford, RAL, UCL
- By photon ray theory in plasma wake (Wilks,
1989) π¬π ππ β β ππ
π
πππ
π
π ππ ππ ππΌ ππ
- π¬π and ππ: change in frequency and the central
frequency of the EM wave
Theory of photon acceleration
JAI-Fest, 10 November 2015 30 University of Oxford, RAL, UCL
- By photon ray theory in plasma wake (Wilks,
1989) π¬π ππ β β ππ
π
πππ
π
π ππ ππ ππΌ ππ
- π¬π and ππ: change in frequency and the central
frequency of the EM wave
- ππ and π: the unperturbed and perturbed electron
density
Theory of photon acceleration
JAI-Fest, 10 November 2015 31 University of Oxford, RAL, UCL
- By photon ray theory in plasma wake (Wilks,
1989) π¬π ππ β β ππ
π
πππ
π
π ππ ππ ππΌ ππ
- π¬π and ππ: change in frequency and the central
frequency of the EM wave
- ππ and π: the unperturbed and perturbed electron
density
- πΌ: the distance in the wakefieldβs frame
Theory of photon acceleration
JAI-Fest, 10 November 2015 32 University of Oxford, RAL, UCL
- From photon ray theory, we get
π π, π = π
π β
ππ©π ππ ππ β ππ π π, π π ππ β ππ ππ
Theory of photon acceleration
JAI-Fest, 10 November 2015 33 University of Oxford, RAL, UCL
- From photon ray theory, we get
π π, π = π
π β
ππ©π ππ ππ β ππ π π, π π ππ β ππ ππ
- Where
π π, πΌ = π¬π π, πΌ ππ π π, πΌ = β π ππ ππ π, πΌ ππΌ ππ
π
πππ
π
π π ππ£π¨πΎ
Theory of photon acceleration
JAI-Fest, 10 November 2015 34 University of Oxford, RAL, UCL
- From photon ray theory, we get
π π, π = π
π β
ππ©π ππ ππ β ππ π π, π π ππ β ππ ππ
- Where
π π, πΌ = π¬π π, πΌ ππ π π, πΌ = β π ππ ππ π, πΌ ππΌ ππ
π
πππ
π
π π ππ£π¨πΎ
- π¬π π, πΌ can be detected using SPIDER, S3I, FDH, or other
interferometry method.
Theory of photon acceleration
JAI-Fest, 10 November 2015 35 University of Oxford, RAL, UCL
- From photon ray theory, we get
π π, π = π
π β
ππ©π ππ ππ β ππ π π, π π ππ β ππ ππ
- Where
π π, πΌ = π¬π π, πΌ ππ π π, πΌ = β π ππ ππ π, πΌ ππΌ ππ
π
πππ
π
π π ππ£π¨πΎ
- π π, πΌ can be obtained by the modified Abel inversion.
Theory of photon acceleration
JAI-Fest, 10 November 2015 36 University of Oxford, RAL, UCL
- Plasma and probe parameters
β π = π. π Γ ππππππ§βπ β ππππππ = πππ π¨π§ (plane wave) β ππ = ππ. π π π
- Driver parameters (electron beam)
β ππππ = π. π Γ ππππππ§βπ β π = π. π ππ§ β π = ππ πππ β πΎ = πππ©
- Repetitive boundaries
3D simulation
JAI-Fest, 10 November 2015 37 University of Oxford, RAL, UCL
3D simulation
JAI-Fest, 10 November 2015 38 University of Oxford, RAL, UCL
- Transform the density profile to the frequency-
shift profile
Results from density profile
Probeβs frequency profile
JAI-Fest, 10 November 2015 39 University of Oxford, RAL, UCL
- Transform the density profile to the frequency-
shift profile
π π, π = β π π
π β
ππ©ππ’ ππ ππ β ππ π π π, π ππ ππ ππ β ππ
Results from density profile
Probeβs frequency profile
JAI-Fest, 10 November 2015 40 University of Oxford, RAL, UCL
- Transform the density profile to the frequency-
shift profile
π π, π = β π π
π β
ππ©ππ’ ππ ππ β ππ π π π, π ππ ππ ππ β ππ
Results from density profile
Inverse Modified Abel transform
Probeβs frequency profile Retrieved density profile
JAI-Fest, 10 November 2015 41 University of Oxford, RAL, UCL
- Qualitative comparison with the retrieved density profile
from the probe pulse
Results
Obtained from the electron density profile Retrieved from the probe pulseβs frequency profile
JAI-Fest, 10 November 2015 42 University of Oxford, RAL, UCL
- Quantitative comparison between the density data and the
measured data.
Results
Less than 15% relative error
JAI-Fest, 10 November 2015 43 University of Oxford, RAL, UCL
- Quantitative comparison between the density data and the
measured data.
Results
Less than 15% relative error
JAI-Fest, 10 November 2015 44 University of Oxford, RAL, UCL
- Quantitative comparison between the density data and the
measured data.
Results
Less than 15% relative error
JAI-Fest, 10 November 2015 45 University of Oxford, RAL, UCL
- Quantitative comparison between the density data and the
measured data.
Results
Less than 15% relative error diffraction
JAI-Fest, 10 November 2015 46 University of Oxford, RAL, UCL
- Frequency-shift profile after interacting with the
wakefield
Limitations
Just after the interaction with the wakefield
JAI-Fest, 10 November 2015 47 University of Oxford, RAL, UCL
- Frequency-shift profile after interacting with the
wakefield
Limitations
Propagates 12 ΞΌm
JAI-Fest, 10 November 2015 48 University of Oxford, RAL, UCL
- Frequency-shift profile after interacting with the
wakefield
Limitations
Propagates 30 ΞΌm
JAI-Fest, 10 November 2015 49 University of Oxford, RAL, UCL
- Frequency-shift profile after interacting with the
wakefield
Limitations
Propagates 30 ΞΌm Initial profile
JAI-Fest, 10 November 2015 50 University of Oxford, RAL, UCL
- Frequency-shift profile after interacting with the
wakefield
Limitations
Propagates 30 ΞΌm Initial profile
Lower value
JAI-Fest, 10 November 2015 51 University of Oxford, RAL, UCL
- As πΎ decreases, interaction length increases.
Limitations
JAI-Fest, 10 November 2015 52 University of Oxford, RAL, UCL
- As πΎ decreases, interaction length increases.
- If the interaction length β₯ diffraction length, it
reads lower value.
Limitations
JAI-Fest, 10 November 2015 53 University of Oxford, RAL, UCL
- As πΎ decreases, interaction length increases.
- If the interaction length β₯ diffraction length, it
reads lower value.
- Diffraction limit:
ππ£π¨ πΎ β₯ ππ πππ
- ππ is the probeβs wavelength
- ππ is the wakefieldβs radius
Limitations
JAI-Fest, 10 November 2015 54 University of Oxford, RAL, UCL
- We have developed a mathematical model for diagnostics of cylindrically
symmetric & relativistic moving objects.
- Using photon acceleration method, we could diagnose the plasma
wakefield density quantitatively.
- More information:
- We are planning to try the diagnostic on laser plasma wakefield
accelerator at Rutherford Appleton Lab and on AWAKE at CERN.
Conclusions
JAI-Fest, 10 November 2015 55 University of Oxford, RAL, UCL
THANK YOU!
Image credit to: http://sf.co.ua/13/07/wallpaper-2937024.jpg
JAI-Fest, 10 November 2015 56 University of Oxford, RAL, UCL
- Frequency and phase shift relation:
π¬π = π¬π ππ
- Phase modulation
π¬π β ππ ππ π¬π ππ
- Detecting π¬π can be done using Spectral Interferometry
method
- Then π¬π can be obtained by π¬π = ππ¬π/ππ
Detecting π¬π and π¬π
JAI-Fest, 10 November 2015 57 University of Oxford, RAL, UCL
- Uses 2 chirped pulses: (1) reference and (2) probe pulses.
- Separated by time separation π.
- The reference isnβt modulated, but the probe is.
- The reference can be obtained from unmodulated part of a
single probe (self-referenced).
- Successfully implemented in some experiments:
β Single-shot supercontinuum spectral interferometry (Kim, et al., 2002) β Frequency domain holography (Matlis, et al., 2006) β Frequency domain streak camera (Li, et al., 2010) β etc.
Spectral interferometry
JAI-Fest, 10 November 2015 58 University of Oxford, RAL, UCL
- If the reference electric field has the profile: ππ(π)
- And the probe has the profile:
ππ π = ππ π β π πππ¬π(πβπ)
- where:
β π¬π π is the probeβs phase modulation
- The spectrum of the pulses:
π± π = π | ππ π |π π + ππ©π π¬ π π + ππ
- Using FFT, we can extract π¬
π π , hence π¬π(π)
Spectral interferometry
- scillating part
JAI-Fest, 10 November 2015 59 University of Oxford, RAL, UCL
- The spectrum of the pulses:
π± π = π | ππ π |π π + ππ©π π¬ π π + ππ
Spectral interferometry
FFT HPF IFFT
Phase
JAI-Fest, 10 November 2015 60 University of Oxford, RAL, UCL
- If the modulation too large β the sideband of the
sideband will overlap the central element and/or the sideband gets out of range.
- If the modulation too small β the noise will spoil
the modulation information.
Spectral interferometry limitations
π ππ
JAI-Fest, 10 November 2015 61 University of Oxford, RAL, UCL
- If the modulation too large β adjust the delay
between two pulses, π, and/or get spectrometer with higher resolution.
- If the modulation too small β reduce the noise