Brief History of Particle in Cell Simulations
3D Potential distribution across the plasma
The following is a (very) brief history of Particle-in-cell simulation techniques developed at the Space Plasma and Plasma Processing Lab (SP3) at the Australian National University, where I did my PhD. SP3 has had a long history of involvement with computer modelling of plasma phenomena. In particular the development and use of Particle-in-cell (PIC) techniques for simulation of low pressure, low temperature, radiofrequency plasmas used in materials processing applications.
- Brief description of Particle-in-Cell techniques
- History of Simulation Usage in SP3
- My research at SP3
- Some web publications
The lab started using PIC in the mid 80s, shortly after the technique was adapted for bounded plasma systems by Prof Birdsall's group at University of California Berkley (for a history of PIC development see C.K. Birdsall (1991) IEEE Trans. Plasma Sci., 19, 65 ). The technique models charged species in the plasma, moving in the electric fields created by applied voltages and inter-particle coloumb forces. PIC codes also include collisions between charged species and background neutral gas atoms, and use boundary conditions which can include external LCR circuits. PIC is particularly adept at modelling non-equilibrium plasmas, in which the particle distributions are not solely a function of the local fields. They have been used to reproduce extremely accurately, experimental measurements of non-Maxwellian electron and ion distributions in both the bulk and at the reactor surfaces. PIC codes are widely used for modelling a variety of plasma types and conditions including rf systems, dc glow and magnetron discharges and plasma immersion ion implantation.
However, PIC is computationally intensive and are inherently
subject to noise on the density and field distributions. This currently limits simulations
to low pressure systems in which the plasma densities are less than 10^17 /m^3. Most PIC
simulations are one dimensional for the same reason (although the use of faster numerical
techniques and supercomputers has seen the development of higher dimensional codes). In
this sense PIC techniques complement fluid methods, which deal well with coherent behaviour
of bulk populations, but are not so useful for determining non-linear behaviour of small
subgroups (for example the heating of a small fraction of the electron density in the sheaths).
Hybrid models, which incorporate a combination of modelling techniques,
increasingly common for more complete simulation of plasma processing systems, and may
include plasma chemistry, plasma-wall interactions and the effect of magnetic fields.
As computers become faster some of the limitations of PIC simulations are being lifted.
PIC codes can be bounded - including walls which determine boundary conditions for the potential solver - or unbounded. In bounded codes modelling plasma processing reactors the boundary conditions can also model external voltage/current sources and LCR "matching" circuits.
Schematic of a plasma reactor with an external matching circuit
Electric fields are solved self-consistently due to the superposition of external (applied) fields and internal fields from charge distributions. This is done by using a non-physical grid across the plasma (the "cell" part of the name) and determining the charge density at each grid position by assigning particles to the grid according to their position, and a weighting scheme. Once the charge density at the grid positions is known the potentials can be calculated using Poisson's equation. Then the electric fields at the grid can be determined and finally fields at the particle positions are determined using an inverse weighting scheme. Particles can then be moved via Newton's equations, using a leap-frog finite differencing method (positions and fields are calculated at integer time-steps, velocities at half time-steps).
Schematic of leap-frog differencing technique
Schematic showing algorithm for 1 time-step of the simulation
Some useful references:
C.K. Birdsall and A.B. Langdon "Plasma Physics via Computer Simulation", McGraw-Hill Book Company, Singapore (1985).
C.K. Birdsall, IEEE Trans. Plasma Sci., 19, 65-85 (1991).
T. Tajima "Computational Plasma Physics: with Applications to Fusion and Astrophysics" Addison-Wesley Publishing Co. (1989).
History of Simulation Usage in the Lab1985
Initial development of an unbounded 1D electrostatic PIC simulation to study beam-plasma interactions
- R.W. Boswell and I.J. Morey, Appl. Phys. Lett., 52, 21 (1988)
- I.J. Morey and R.W. Boswell, Phys. Fluids B, 1, 1502 (1988)
Development of a bounded 1D simulation of a planar rf discharge modelling a hydrogen-like gas, used to study basic phenomena such as current/voltage relationships, density and potential distributions and sheath heating of the electrons.
- D. Vender and R.W. Boswell, IEEE Trans. Plasma Sci.,18, 725 (1990)
- D. Vender and R.W. Boswell, IEEE Trans. Plasma Sci., 19, 141 (1991)
- D. Vender and R.W. Boswell, J. Vac. Sci. Technol. A, 10, 1331 (1992)
- R.W. Boswell, A.J. Lichtenberg and D. Vender, IEEE Trans. Plasma Sci., 20, 62 (1992)
1D electrostatic PIC simulation of asymmetric rf discharges modelling hydrogen and argon gases with realistic cross-sections.
Use of a 2D PIC/fluid hybrid simulation with flexible boundary conditions, to model inductive/wave-sustained sources (including helicon reactors). Also used to look at diffusion of evaporant species through a plasma. (see the page on the HARE reactor at the ANU HARE )
- R.K. Porteous and D.B. Graves, IEEE Tran. Plasma Sci.,19, 204 (1991).
- R.K. Porteous, H.M. Wu and D.B. Graves, Plasma Sources Sci. Technol., 3, 25 (1994)
- B. Higgins, A. Durandet and R. Boswell, J. Vac. Sci. Technol. B 13(2), 192 (1995)
- Multiple author paper, Rev. Sci. Instrum. 66(3), 2908 (1995)
Modelling plasma breakdown and decay, in pulsed rf discharges, using both symmetric and asymmetric reactor geometries.
- D. Vender, H.B. Smith and R.W. Boswell, J. Appl. Phys., 80, 4292 (1996)
- H.B. Smith, C. Charles and R.W. Boswell, J. Appl. Phys., 82, 561 (1997).
My research: 1997-1999My main area of interest while I was a post-doctoral research scientist at ANU, was in modelling breakdown and decay in rf systems, in order to better understand pulsed discharges. The simulations were compared with results from a small experimental experimental reactor in order to make a detailed study of breakdown and approach to equilibrium in a low pressure, capacitively coupled, rf argon discharge Studies have been made of the evolution of the bias voltage in asymmetric geometry, which indicate that previously unknown behaviour of the currents and plasma potential are involved (see " Charging of the blocking capacitor in an asymmetric discharge" for details ).
The simulation was also used to investigate time-dependent effects such as the coupling of power into the electrons and ions and the growth of the charged particle densities before, and immediately after, sheath formation. Results from the simulation have been extremely important in building general analytic models of the plasma behaviour (see " Simulation of plasma breakdown in a low pressure rf plasma" ).
The effect on multiple pulses on parameters such as the plasma density, currents, potentials and power coupled into the plasma have also been investigated (see " Pulsing a low pressure rf plasma" ). Pulsed plasmas are of increasing interest in materials processing applications, for improved film properties in deposition and increased uniformity for etching. They have also been used to reduce powder formation in chemical vapour deposition reactors. In order to understand and extend their application in commercial ventures it is essential to understand the physical processes taking place during pulsing.
A model of an external "L-series" matching network was developed to include in the plasma model.
The matching network is used in the experimental system to efficiently
couple power from the generator into the plasma. In the process of doing so the network has a
large influence on the current/voltage waveforms applied to the reactor electrode and hence
on the behaviour of the plasma. In particular during breakdown the plasma impedances are
changing rapidly, which consequently effects the power supplied to the plasma.In order to
better understand how the matching network and the plasma interact, and to better model the
experimental system, it was decided to explicitly include the matching network in a PIC
simulation. Ultimately it is
hoped that by better understanding the influence of the matching network on the plasma formation
and growth, improvements to the circuit design can be made, in particular for pulsed systems
in which the plasma experiences predominately non-equilibrium conditions.