starting in October 2016 (here are the 2018 projects)

Application deadline: 22 January 2016

Oxford Plasma Theory Group welcomes applications for DPhil studies and research in plasma physics in the areas of magnetic confinement fusion (MCF) and plasma astrophysics.

Potential supervisors:
Prof Michael Barnes (projects 1, 4), Prof Gianluca Gregori (project 5; also here), Prof Felix Parra Diaz (project 2), Prof Subir Sarkar (project 5), Prof Alexander Schekochihin (projects 3, 4, 6, 7; also here)


Projects in Fusion Plasmas


These three projects are all offered jointly with researchers at the U.K.A.E.A. Culham Centre for Fusion Energy (CCFE). They also form part of a Plasma Turbulence Collaboration between Oxford Plasma Theory Group, CCFE and York Plasma Institute (YPI) at the University of York. York also offers PhD projects within this Collaboration --- if you wish to be considered for those, please apply to York as well as to Oxford (the competitions for admission at the two Universities are separate and students will follow different training programmes, but all will be involved in the Collaboration).

At the application stage, you are not required to (although you may if you wish) indicate which project you prefer --- we will consider all applicants purely on intellectual merit. If you are made an offer of admission, we will strive to give you the opportunity to work on the project of your choice.
Note that the project descriptions given below are not set in stone and we are willing to discuss modifications and adjustments to them that might better reflect your interests and inclinations.

1. Turbulence suppression in magnetic confinement fusion plasmas
Supervisor: Prof Michael Barnes
CCFE co-supervisor: Dr Michele Romanelli

[see also this project at York]

In magnetic confinement fusion, a strong magnetic field is used to confine a hydrogenic gas as it is heated to over a hundred million degrees. At these temperatures the gas becomes ionised and the constituent ions have enough thermal energy to tunnel through the Coulomb barrier and convert mass into energy via fusion. This ionised gas, or plasma, generates and responds to electromagnetic fields and demonstrates collective behavior including unstable modes of oscillation that result in turbulent mixing. The unstable modes are driven by gradients in the equilibrium plasma density, temperature, and flow; in turn, the equilibrium is itself constrained by the resultant turbulent mixing. This nonlinear feedback loop ultimately controls the fusion energy output and exhibits a variety of rich and puzzling dynamical features that span a wide range of space-time scales. The aim of this project is to address a practical and fundamental open question regarding these multi-scale turbulent dynamics: What leads to the appearance of turbulence-free regions in some current experiments, and can we use these regions to optimise fusion energy production? The project will require a combination of analytical theory and cutting-edge numerical techniques and should benefit from comparisons with data from MAST and JET, two large fusion experiments housed nearby at the Culham Laboratory, as well as a strong collaboration with scientists working on the new JT-60SA experiment in Japan.

2. Interaction between turbulent plasma fluctuations and high-frequency plasma waves

Supervisor: Prof Felix Parra Diaz
CCFE co-supervisor: Dr Jon Hillesheim
York collaborator: Dr David Dickinson
[see also this project at York]

The performance of magnetic confinement devices is limited by electromagnetic fluctuations driven by unavoidable plasma instabilities. These fluctuations, commonly known as plasma turbulence, spoil the confinement provided by the strong external magnetic fields, and reduce the fusion energy yield. In this project, the student will study how these intrinsic turbulent fluctuations (with characteristic frequencies of the order of 10kHz), interact with high-frequency waves (of the order of 100 GHz) that are launched into the plasma. The nonlinear scattering of the high-frequency waves by the plasma fluctuations is a fundamental problem in plasma physics that requires sophisticated semi-analytical treatments to bridge the large time-scale separation. The study of the interaction between low- and high-frequency fluctuations is crucial for plasma turbulence diagnostics that use the high-frequency waves to probe the fusion plasma, which is at temperatures hotter than the centre of the sun and inaccessible to physical probes. The student will apply the developed theory to the Doppler Back-Scattering (DBS) systems in JET and MAST, the two experiments at the Culham Centre for Fusion Energy. In addition to having access to the large database accumulated in the past three years, the student will have the opportunity to work with the newly developed DBS system for MAST-Upgrade. This research is highly relevant to the future of fusion energy because DBS diagnostics are planned for ITER, the international collaboration to demonstrate the feasibility of fusion power currently under construction in France, and similar diagnostics will be necessary in future fusion reactors to monitor the plasma in real time.

3. Zonal flow-drift wave interactions in ion-scale turbulence in differentially rotating tokamak plasmas: BES measurements and gyrokinetic models
Supervisor: Prof Alexander Schekochihin
CCFE co-supervisor: Dr Anthony Field
York collaborator: Dr Istvan Czeigler
[see also this project at York]

While electromagnetic turbulence in fusion plasmas is widely accepted to be the main culprit in undermining confinement (by vigorously transporting energy from the core of the plasma to its edge), until relatively recently most of our understanding of this turbulence or even just of our knowledge of what it "looks like" had been based on indirect measures of its effects (primarily heat transport) and on gyrokinetic numerical simulations. The Beam-Emission-Spectroscopy (BES) diagnostics have allowed us to probe turbulent fluctuations directly. In particular, BES data from MAST has been successfully used to determine statistical properties of ion-Larmor-scale turbulence and thereby test theoretical predictions and gyrokinetic simulations. In particular, it appears that if the plasma is differentially rotated, i.e., if locally the turbulence is embedded in a sheared plasma flow, the nature of the turbulence changes and the heat trasport can be suppressed. We are far from complete understanding of this process or indeed the basic structure of the turbulence. Simulations show that the experimentally relevant near-marginal turbulent state is crucially regulated by the interaction between the equilibrium velocity shear, zonal flows (coherent shear flows generated by the turbulence itself) and wave-like electromagnetic fluctuations that form the bulk of the turbulence (called "drift waves"). This project will focus on this interaction by determining the scaling of the turbulent and equilibrium parameters, in particular with flow shear and collisionality (which is a control parameter for zonal flows), and by comparison with results of gyrokinetic numerical modeling. The project involves both theoretical and experimental work (which may mean developing analysis techniques and/or interacting with actual hardware, depending on your inclinations). In particular, analysis techniques are to be developed and tested for the detection of fluctuating 2D flows, with the aim of detecting localised, self-generated zonal flows. Both the MAST tokamak and its BES system are being upgraded (due to come online in 2016/17) and so the student working on this project will have a unique opportunity to extract new physics results, including planning dedicated experiments. 

Background Reading:
1. M. Barnes, F. I. Parra and A. A. Schekochihin, “Critically balanced ion temperature gradient turbulence in fusion plasmas,” Phys. Rev. Lett. 107, 115003 (2011)
2. A. A. Schekochihin, E. G. Highcock, and S. C. Cowley,
Subcritical fluctuations and suppression of turbulence in differentially rotating gyrokinetic plasmas,”  Plasma Phys. Control. Fusion 54, 055011 (2012)
3. Y.-c. Ghim, A. A. Schekochihin, A. R. Field, I. G. Abel, M. Barnes, G. Colyer, S. C. Cowley, F. I. Parra, D. Dunai, S. Zoletnik, and the MAST Team,
Experimental signatures of critically balanced turbulence in MAST,Phys. Rev. Lett. 110, 145002 (2013)
4. Y.-c. Ghim, A. R. Field, A. A. Schekochihin, E. G. Highcock, C. Michael, and the MAST Team,
Local dependence of ion temperature gradient on magnetic configuration, rotational shear and turbulent heat flux in MAST,” Nucl. Fusion 54, 042003 (2014)

Projects in Plasma Astrophysics

Candidates interested in any of these three projects or generally in plasma astrophysics, astrophysical turbulence and/or dynamo theory are welcome to get in touch with prospective supervisors for further information. A more bespoke project can be designed to align with the inclinations and interests of the student (for example how much emphasis is placed on analytical vs. numerical methods or kinetic theory vs. fluid dynamics, etc., is negotiable).  

4. Free energy flows in turbulent astrophysical plasmas.
Supervisors: Prof Michael Barnes and Prof Alexander Schekochihin

In magnetised astrophysical plasmas, there is a turbulent cascade of electromagnetic fluctuations carrying free energy from large to small scales. The energy is typically extracted from large-scale sources (e.g., in the solar wind, the violent activity in the Sun’s corona; in accretion discs, the Keplerian shear flow; in galaxy clusters, outbursts from active galactic nuclei) and deposited into heat – the internal energy of ions and electrons. In order for this dissipation of energy to happen, the energy must reach small scales – in weakly collisional plasmas, these are small scales in the 6D kinetic phase space, i.e., what emerges are large spatial gradients of electric and magnetic fields and large gradients in velocity of the particle distribution functions. This gives rise to two very intriguing questions: (1) how does the energy flow through the 6D phase space and what therefore is the structure of the fluctuations in this space: their spectra, phase-space correlation functions etc. (these fluctuations are best observed in the solar wind, but these days we can also measure density and magnetic fluctuations in galaxy clusters, via X-ray and radio observations; turbulent velocities will become observable as well when ASTRO-H launches in 2015); (2) when turbulent fluctuations are dissipated into particle heat, how is their energy partitioned between various species of particles that populate the plasma: electrons, bulk ions, fast heavy particles (e.g., alpha particles and heavy ions). The latter question is particularly important for extragalactic plasmas because all we can observe is radiation from the particles and knowing where the internal energy of each species came from is key to constructing and verifying theories both of turbulence and of macroscale dynamics and thermodynamics. This project has an analytical and a numerical dimension (which of these will dominate depends on the student’s inclinations). Analytically, we will work out a theory of phase space cascade at spatial scales between the ion and electron Larmor scales (we have done some preliminary work, so we know how to start off on this calculation, but obviously at some point we’ll be wading into unchartered waters). Numerically, we will simulate this cascade using “gyrokinetic” equations – an approach in which we average over the Larmor motion and calculate the distribution function of “Larmor rings of charge” rather than particles (this reduces the dimension of phase space to 5D, making theory more tractable and numerics more affordable). Our key objective will be to work out the energy partition between ions and electrons – this question has fascinated theoreticians for at least 15 years, but we now have the tools to sort it out. 

Background Reading:
1. A. A. Schekochihin et al., “Astrophysical gyrokinetics: kinetic and fluid turbulent cascades in magnetized weakly collisional plasmas,” Astrophys. J. Suppl. 182, 310 (2009)
2. G. G. Howes et al., “Gyrokinetic simulations of solar wind turbulence from ion to electron scales,” Phys. Rev. Lett. 107, 035004 (2011)
3. T. Tatsuno et al., “Nonlinear phase mixing and phase-space cascade of entropy in gyrokinetic plasma turbulence,” Phys. Rev. Lett. 103, 015003 (2009)
4. G. G. Howes et al., “Kinetic simulations of magnetized turbulence in astrophysical plasmas,” Phys. Rev. Lett. 100, 065004 (2008)

5. Second-order Fermi acceleration by MHD turbulence in young supernova remnants
Supervisors: Prof Subir Sarkar and Prof Gianluca Gregori
(for this project, you may also apply for a DPhil in Atomic and Laser Physics)

The origin of high-energy cosmic rays still remains a mystery, a century after their discovery. Although several different processes and astrophysical sites may be implicated, there is a consensus that diffusive shock acceleration in Galactic supernova remnants accelerates the bulk of hadronic cosmic rays, up to hundreds of TeV in energy. However, observations of radio synchrotron emission by high-energy electrons in young supernova remnants are better understood if the electrons are accelerated by turbulence behind a decelerating shock wave via the 2nd-order Fermi process. In this project, we propose to develop numerical simulations of such turbulence in order to determine its spectrum and thus calculate the time evolution of the relativistic-electron spectrum. The simulations will be validated against laboratory experiments in magnetohydrodynamic turbulence, using the largest laser systems in the world, and the predictions of the radio emission will be compared in detail with observations of Cassiopeia A as well as of other young supernova remnants such as SN1987A and G1.9+0.3.

Background Reading:
V. Petrosian Stochastic Acceleration by Turbulence,” Space Sci. Rev. 173, 535 (2012)

6. Magnetogenesis in the intergalactic plasma
Supervisor: Prof Alexander Schekochihin

In weakly collisional plasmas such as those found in galaxy clusters (and in many other places, from the heliosphere to accretion discs around black holes), pressure is a tensor rather than a scalar – and there are perpendicular and parallel pressures, with respect to the local direction of the magnetic field. As it turns out, every time such a plasma is moved around (as a fluid), the magnetic field that threads it will change (because the field moves with the plasma) and every time that happens, pressure anisotropy is created, i.e., perpendicular and parallel pressures become different. This situation turns out to be violently unstable, a sea of small-scale fluctuations emerges and the plasma thereby contrives to limit the pressure anisotropy at a certain level, which itself depends on the magnitude of the magnetic field. This means the magnetic field can only evolve in a fairly constrained way, very differently to how it behaves in conventional conducting fluids (e.g., liquid metals or very collisional plasmas) described by the equations of magnetohydrodynamics (MHD). The reason this is interesting is that (1) no one knows how exactly the field evolves (the “industry-standard” models of most astrophysical plasmas are by and large MHD, so not appropriate to the case with pressure anisotropies) – this means there is new physics to be discovered; (2) the physical regime in which this is all relevant is precisely the regime in which most plasmas in the Universe are – and so we can't understand the magnetic field structure we observe in the intergalactic medium, predict how this medium moves (more broadly, its large-scale dynamics and thermodynamics) or know where the field came from in the first place without a good model for the magnetic-field dynamics. In this project, we will aim to develop such a model.

Background Reading:
1. F. Mogavero and A. A. Schekochihin, “Models of magnetic-field evolution and effective viscosity in weakly collisional extragalactic plasmas,” Mon. Not. R. Astron. Soc. 440, 3226 (2014)

2. A. A. Schekochihin et al., “Magnetofluid dynamics of magnetized cosmic plasma: firehose and gyrothermal instabilities,” Mon. Not. R. Astron. Soc. 405, 291 (2010)
3. M. S. Rosin et al., “A nonlinear theory of the parallel firehose and gyrothermal instabilities in a weakly collisional plasma,” Mon. Not. R. Astron. Soc. 413, 7 (2011)
4. M. W. Kunz, A. A. Schekochihin, and J. M. Stone, “Firehose and mirror instabilities in a collisionless shearing plasma,” Phys. Rev. Lett. 112, 205003 (2014)
5. F. Rincon, A. A. Schekochihin, and S. C. Cowley, “Nonlinear mirror instability,” Mon. Not. R. Astron. Soc. 447, L45 (2015)

7. Turbulent cascade through the ion cyclotron resonance in space plasmas.
Supervisor: Prof Alexander Schekochihin

This project is related to the general theme of Project 1, but focuses on a single problem that requires analytical handling in the fully kinetic 6D framework and a programme of numerical simulations of the full kinetic equation for a plasma. The computational challenges involved in 6D simulations are enormous. The problem proposed for this project is therefore carefully chosen to be both of fundamental importance and solvable in a definitive way with existing resources. In environments such as the solar wind (and, as far as we can fathom, more generally in cosmic plasmas from heliospheric to extragalactic), the turbulent cascade of electromagnetic fluctuations reaches ion Larmor spatial scales while frequencies remain smaller than the ion cyclotron frequency (these are the radius and frequency of the particle gyration in a magnetic field). Under these conditions, the “gyrokinetic” approximation (5D rather than 6D) holds. However, in many parameter regimes, the ion cyclotron frequency is reached at sub-Larmor scales. Linear theory tells us that in the narrow region of frequencies around the ion cyclotron, electromagnetic fluctuations fall into resonance with Larmor motion and thus can transfer their energy into the particle motion, leading to ion heating. The nonlinear question is how much (if any) of the turbulent free energy will be sucked into this resonance and how much will, via nonlinear coupling, simply leapfrog it and continue cascading until converted into electron heat at the electron Larmor scale? This question cannot be resolved via any low-frequency model (like gyrokinetics) and will therefore be answered via 6D fully kinetic simulations restricted to a narrow window of scales around the resonance. Theory is falsifiable in this area: one can (and does) measure ion distribution functions in the solar wind (although the physics results apply more broadly). If the student is so inclined, we may get involved directly with analysis of such measurements. This project is likely to involve a collaboration with Prof Matthew Kunz at Princeton University.

Background Reading:
1. G. G. Howes et al., “A model of turbulence in magnetized plasmas: implications for the dissipation range in the solar wind,” J. Geophys. Res. 113, A05103 (2008)
2. A. A. Schekochihin et al., “Astrophysical gyrokinetics: kinetic and fluid turbulent cascades in magnetized weakly collisional plasmas,” Astrophys. J. Suppl. 182, 310 (2009)