Abstract:The physical processes at the interface of a low-temperature plasma and a solid are extremely complex. They involve a huge number of elementary processes in the plasma, in the solid as well as charge, momentum and energy transfer across the interface. In the majority of plasma simulations these surface processes are either neglected or treated via phenomenological parameters such as sticking coefficients, sputter rates or secondary electron emission coefficients. However, those parameters are known only in some cases, so such an approach is very inaccurate and does not have predictive capability. Therefore, improvements are highly needed. In this paper we briefly summarize relevant theoretical methods from solid state and surface physics that are able to contribute to an improved simulation of plasma-surface interaction in the near future. Even though the (quantum–mechanical) equations of motion for the participating charged and neutral particles are known, in principle, full ab initio quantum simulations are feasible only for extremely short times and/or small system sizes. A substantial simplification is achieved when electronic quantum effects are not treated explicitly. Then one arrives at much simpler semi-classical molecular dynamics (MD) simulations for the heavy particles that have become the main workhorse in surface science simulations. Using microscopically (i.e., density functional theory) founded potentials and force fields as an input, these MD simulations approach the quality of ab initio simulations, in many cases. However, despite their simplified nature, these simulations require a time step that is of the order or below one femtosecond making it prohibitive to reach experimentally relevant scales of seconds or minutes and system sizes of micrometers. To bridge this gap in length and time scales without compromising the first principles character and predictive power of the simulations, many physical and computational strategies have been put forward in surface science. This paper presents a brief overview on different methods and their underlying physical ideas, and we compare their strengths and weaknesses. Finally, we discuss their potential relevance for future plasma-surface simulations. The first class are 'acceleration' techniques that include metadynamics, hyperdynamics, temperature accelerated dynamics, collective variable driven hyperdynamics and others. Recently we have presented a novel approach: selective process acceleration (Abraham et al 2016 J. Appl. Phys. 119 185301) which we discuss in some more detail. The second promising route to longer accurate simulations is dynamical freeze out of dominant modes which we have introduced recently for the simulation of neutral atom sticking on a metal surface (Filinov et al 2018 Plasma Sources Sci. Technol. 27 064003). In this article we give a more general view on this method that allows to accurately combine first principles MD simulations with semi-analytical models and discuss possible applications that are of potential relevance for plasma physics.

Kinetic Monte Carlo simulations of plasma-surface reactions on heterogeneous surfaces

A Practical Guide to Surface Kinetic Monte Carlo Simulations

Monte Carlo and Kinetic Monte Carlo Models for Deposition Processes: A Review of Recent Works

Modelling and Simulation of Heterogeneous Reactions with Statistical Particle Methods

Molecular Dynamics for Low Temperature Plasma-Surface Interaction Studies

Ab Initio-Based Stochastic Simulations of Kinetic Processes on Surfaces

Cluster Diffusion and Coalescence on Metal Surfaces: applications of a Self-learning Kinetic Monte-Carlo method

Beyond Mean-Field Microkinetics: Toward Accurate and Efficient Theoretical Modeling in Heterogeneous Catalysis

Towards an integrated modeling of the plasma-solid interface

Acceleration and sensitivity analysis of lattice kinetic Monte Carlo simulations using parallel processing and rate constant rescaling.

Monte Carlo hydrodynamic simulation of neutral radical transport in low pressure remote plasma activated chemical vapor deposition

Monte Carlo Simulation Study of Electron Interaction with Solids and Surfaces

Diffusion Limited Kinetics In Reactive Systems

Kinetic modelling of heterogeneous catalytic systems

Kinetic Monte Carlo Simulation of Thin Film Growth

Numerical Simulation of Electrochemical Adsorption Process by Kinetic Monte Carlo Methods

Highly detailed computational study of a surface reaction model with diffusion: four algorithms analyzed via time-dependent and steady-state Monte Carlo simulations

Extending first principle plasma-surface simulations to experimentally relevant scales

Reactions at surfaces studied by ab initio dynamics calculations

General concepts, assumptions, drawbacks, and misuses in kinetic Monte Carlo and microkinetic modeling simulations applied to computational heterogeneous catalysis

Plasma Catalysis Modeling: How Ideal Is Atomic Hydrogen for Eley-Rideal?