QDB: a new database of plasma chemistries and reactions

TL;DR

QDB is a comprehensive database platform that collates, validates, and identifies gaps in plasma chemistry data, facilitating research and industrial applications by providing complete and validated reaction datasets.

Contribution

The paper introduces QDB, a novel platform with a new data model for collating and validating plasma chemistry datasets, including initial chemistry sets for specific gas mixtures.

Findings

Validated chemistry sets for SF₆/CF₄/O₂ and SF₆/CF₄/N₂/H₂ are presented.

QDB identifies key unreported reactions and data gaps in plasma chemistry.

The platform enables data exchange and gap filling using theoretical methods.

Abstract

One of the most challenging and recurring problems when modelling plasmas is the lack of data on key atomic and molecular reactions that drive plasma processes. Even when there are data for some reactions, complete and validated datasets of chemistries are rarely available. This hinders research on plasma processes and curbs development of industrial applications. The QDB project aims to address this problem by providing a platform for provision, exchange, and validation of chemistry datasets. A new data model developed for QDB is presented. QDB collates published data on both electron scattering and heavy-particle reactions. These data are formed into reaction sets, which are then validated against experimental data where possible. This process produces both complete chemistry sets and identifies key reactions that are currently unreported in the literature. Gaps in the datasets can be filled using established theoretical methods. Initial validated chemistry sets for SF$_6$/CF$_4$/O$_2$ and SF$_6$/CF$_4$/N$_2$/H$_2$ are presented as examples.

Tables

Table 1. Table 1: Classification of electron collision processes considered in QDB.

Abbreviation	Types of Reaction	Description
EDX	deexcitation	$e+A^{\ast }\to e+A$
EEL	elastic scattering	$e+A\to e+A$
EIN	ionization	$e+A\to e+A^{+}+e$
EMT	momentum transfer
ERR	radiative recombination	$e+A^{+}\to A+h\nu$
EDR	dissociative recombination	$e+A{B}^{+}\to A+B$
EDS	dissociation	$e+AB\to e+A+B$
EDA	dissociative attachment	$e+AB\to A+B^{-}$
EDE	dissociative excitation	$e+AB\to A^{\ast }+B+e$
EDI	dissociative ionization	$e+AB\to A^{+}+B+2e$
EEX	electron-impact electronic excitation	$e+A\to e+A^{\ast }$
ECX	change of excitation	$e+A^{\ast }\to e+A^{\ast \ast }$
ERC	recombination (general)	$e+A^{+z}\to A^{+\left(z-1\right)}$
EDT	electron attachment	$e+A+B\to A+B^{-}$
EVX	electron-impact vibrational excitation	$e+A\to e+A\left[v=\ast \right]$
ETS	electron total scattering	$e+A\to e+\Sigma A$
ETI	electron total ionisation	$e+A\to e+e+\Sigma A^{+}$
ETA	electron total attachment	$e+A\to \Sigma A^{-}$

Table 2. Table 2: Classification of heavy-particle processes considered in QDB.

Abbreviation	Types of Reaction	Description
HGN	associative electron detachment	$A^{-}+B\to AB+e$
HCX	charge transfer	$A^{+}+B\to A+B^{+}$
HIR	heavy-particle interchange	$A+BC\to AB+C$
HAS	association	$A+B\to AB$
HIN	heavy-particle collisional ionization	$A+B\to A+B^{+}+e$
HIA	heavy-particle association and ionization	$A+B\to A{B}^{+}+e$
HPI	Penning ionization	$A+B^{\ast }\to A^{+}+B+e$
HNE	neutralization	$e+B^{-}\to B+2e$
HMM	ions recombination	$A^{-}+B^{+}\to A+B$
HDS	heavy-particle collisional dissociation	$AB+C\to A+B+C$
HDX	heavy-particle collisional deexcitation	$A+B^{\ast }\to A+B$
HDN	heavy-particle dissociative neutralization	$A{B}^{-}+C^{+}\to A+B+C$
HDC	heavy-particle dissociation and charge transfer	$AB+C^{+}\to A^{+}+B+C$
HDI	heavy-particle dissociation and ionization	$AB+C^{\ast }\to A^{+}+B+C$
HEX	heavy-particle excitation	$A+B\to A+B^{\ast }$
HED	heavy-particle electron detachment	$A^{-}+B\to A+B+e$

Table 3. Table 3: Classification of photon processes considered in QDB.

Abbreviation	Types of Reaction	description
PDS	photodissociation	$AB+h\nu \to A+B$
PEX	photoexcitation	$A+h\nu \to A^{\ast }$
PRD	radiative decay	$A^{\ast }\to A+h\nu$

Table 4. Table 4: Data compilations used as input to QDB.

Lead author	System
Itikawa	N2 [39], H2O [40], CH4 [16], CO [41], CO2 [42], H2 [43], O2 [44, 45], O [46].
Christophorou	CF4 [47, 48], CHF3 [49, 47], C2F6 [47], C3F8 [47], CCl2F2 [50], Cl2 [47], SF6 [51, 52], C2F6 [47], BCl3 [47], CF3I [53], C4F8 [54]
Janev	H2 [55, 56], CXHY [17, 57, 58], SiH4 [59],
Phelps	N2 [60, 61], SF6 [62], O2 [63]
Bartschat	Ar [64], F [65], B [66]
Quantemol	BF3 [67], C2H2, C3, C3H4 [68], C3N, CF4, CH, CH4 [69], HCN/HNC [70],CONH3, COS, CS [71], CaF, F2O,
	H2S, HBr, HCHO, HCP, Kr, NH3,O3, PH3, SO2, SiF2, SiH4, SiO [72]
Tennyson	C2 [73], C2OH6 [74], CF [75, 76], CF2 [77, 76], CO [78, 79, 80], CO2 [81], H2 [82, 83, 84], H2O [85], N2 [86, 87],
	N2O [88], NO2 [89], O2 [90, 91, 92, 93]

Table 5. Table 5: Data sources for electron collision processes included in QDB classified by reaction type.

Reaction	Data source
Deexcitation (EDX)	[94, 95, 96, 31, 97, 32, 98, 27, 99, 100, 22, 101, 102]
Elastic Scattering (EEL)	[103, 48, 46, 43, 45, 42, 96, 104, 52, 100, 105, 106, 21, 107, 26, 27, 28, 108, 31, 109, 110, 111]
	[32, 112, 98, 113, 114, 39, 115, 116, 117, 118, 65, 119, 120, 121, 49, 122, 40, 123, 124, 125, 126, 16]
Ionization (EIN)	[103, 48, 127, 45, 94, 42, 96, 97, 104, 128, 129, 106, 107, 130, 26, 131, 132, 95, 133, 134, 135, 31, 51, 32, 136]
	[124, 113, 137, 98, 138, 139, 140, 141, 119, 142, 121, 120, 33, 122, 40, 123, 143, 144, 65, 145, 39, 146, 102]
Momentum Transfer (EMT)	[116, 47, 48, 61, 45, 43, 95, 51, 40, 65, 16, 42]
Dissociative Recombination (EDR)	[19, 94, 147, 148, 149, 150, 105, 151, 152, 21, 153, 154, 108, 112, 155, 156, 117, 157]
	[34, 158, 125]
Dissociation (EDS)	[103, 48, 96, 159, 128, 129, 106, 21, 153, 130, 51, 110, 160, 113, 137, 155, 161, 162, 61]
	[119, 49, 163, 122, 40, 123, 124, 62]
Dissociative Attachment (EDA)	[103, 48, 164, 45, 165, 42, 43, 106, 21, 107, 166, 51, 167, 168, 110, 119, 49, 122, 158, 40, 169, 123, 124, 16]
Dissociative Excitation (EDE)	[170, 171, 21]
Dissociative Ionization (EDI)	[103, 172, 48, 45, 43, 53, 42, 173, 174, 96, 128, 21, 131, 51, 135, 175, 176, 54, 113, 137, 138]
	[139, 49, 163, 16, 40, 102, 141, 39, 146, 177]

Table 6. Table 6: Data sources for chemical reactions included in QDB classified by reaction type.

Reaction	Data source
Associative Electron Detachment (HGN)	[178, 105, 18, 122, 158, 179, 180]
Charge Transfer (HCX)	[19, 181, 94, 182, 147, 183, 184, 185, 186, 187, 179, 188, 180, 189, 190, 191, 192, 18]
	[22, 101, 193, 23, 29, 194, 166, 110, 195, 196, 197, 198, 33, 199, 200, 157, 122, 158]
Electron-Impact Electronic Excitation (EEX)	[103, 127, 45, 43, 201, 202, 203, 204, 147, 60, 107, 26, 27, 95, 31]
	[110, 168, 167, 98, 118, 61, 119, 121, 40, 62, 205, 124, 65, 206, 207, 39, 16, 102]
Heavy-Particle Interchange (HIR)	[147, 208, 184, 209, 210, 151, 152, 211, 212, 213, 18, 214, 215, 216]
	[217, 218, 219, 220, 110, 198, 221, 222, 223, 224, 225, 226, 19, 227, 150, 189, 228, 21]
	[229, 230, 24, 231, 232, 233, 25, 122, 234, 179, 190, 192, 235, 236, 237, 238, 29, 239]
	[240, 241, 242, 243, 244, 245, 246, 247, 248, 196, 33, 249, 157, 250, 251, 158, 252, 253]
Association (HAS)	[19, 181, 94, 208, 184, 129, 150, 99, 254, 180, 255, 256, 257, 258, 21, 192, 235, 18, 236]
	[259, 260, 24, 261, 232, 262, 242, 30, 198, 33, 249, 263, 264, 122, 234, 225, 252, 158, 35]
Heavy-Particle Collisional Ionization (HIN)	[265, 94, 266, 99, 267]
Heavy-Particle Association And Ionization (HIA)	[268, 269, 122]
Penning Ionization (HPI)	[24, 268, 270, 122, 179, 22, 101, 180]
Neutralization (HNE)	[119, 221, 18, 128, 26, 271, 180]
Ions Recombination (HMM)	[19, 94, 208, 186, 272, 150, 273, 21, 22, 18, 101, 24, 274, 110, 178, 275, 276, 122, 158]
Heavy-Particle Collisional Dissociation (HDS)	[19, 150, 181, 147, 208, 184, 277, 185, 278, 99, 187, 179, 279, 280, 281, 282]
	[283, 271, 235, 22, 101, 23, 262, 284, 33, 249, 122, 234, 252]
Heavy-Particle Collisional Deexcitation (HDX)	[29, 19, 181, 249, 285, 33, 30, 122, 208, 185, 150, 286]
Heavy-Particle Dissociative Neutralization (HDN)	[24, 178, 273, 151, 287, 288, 275, 21, 122, 150, 18]
Heavy-Particle Dissociation & Charge Transfer (HDC)	[19, 184, 187, 179, 289, 189, 20, 290, 21, 190, 22, 18, 101, 193]
	[25, 110, 98, 198, 33, 122, 158, 234]
Heavy-Particle Dissociation And Ionization (HDI)	[291]

Table 7. Table 7: List and validation status of QDB Chemistries, November 2016

ID	# reactions	Mixture	Validated?
C3	146	N2/H2
C4	63	Ar/H2
C5	155	O2/H2
C6	194	SF6/O2	Yes
C7	208	CF4/O2	Yes
C8	61	SF6	Yes
C9	81	CF4	Yes
C10	409	CF4/O2/H2/N2
C11	193	C4F8	Yes
C13	49	SiH4
C14	73	SiH4/NH3
C15	59	Ar/O2
C16	412	Ar/O2/C4F8	Yes
C17	207	SiH4/Ar/O2
C18	26	Ar/Cu
C19	70	Cl2/O2/Ar
C20	55	Ar/BCl3/Cl2
C21	238	Ar/NH3
C22	187	CH4/H2
C23	81	C2H2/H2
C24	286	CH4/NH3
C25	180	C2H2/NH3
C26	128	He/O2
C27	562	CF4/CHF3/H2/Cl2/O2/HBr
C28	590	CH4/N2
C29	334	SF6 / CF4/ O2	Yes
C30	20	Ar/Cu/He
C31	104	Ar/NF3
C32	192	SF6/CF4/N2/H2	Yes

Table 8. Table 8: Rating scheme for QDB chemistry sets

Ranking	Description of comparison conditions
1	Self-consistent but behaviour differs from available measurements.
2	Not yet compared; no suitable measurements found.
3	Comparison with measurements for the same process conditions versus one variable:
	power, gas flow or pressure
4	Comparison with some measurements for different process conditions
	(more than one comparison) e.g. validation for different pressure regimes.
5	Chemistry tested using more than a program:
	E.g. Quantemol-P and Quantemol-VT or another plasma simulation model.

Table 9. Table 9: Conditions assumed in the model of SF 6 -CF 4 -O 2 etching.

Parameter	Value
Gas Pressure	500 and 700 mTorr
Gas Flow Rate
SF6	800 sccm
CF4	150 sccm
O2	10 sccm
Power	1000 to 2000 W
Substrate	Si

Table 10. Table 10: The conditions for processing in the model for SF 6 /CF 4 /N 2 /H 2 .

Parameter	Value
Gas Pressure	500 mTorr
Gas Flow Rate
SF6	800 sccm
CF4	150 sccm
N2/H2	10 sccm
Power	1500 W
Substrate	Si

Table 11. Table 11: Recognised quantum state designations in QDB

Generic excited state	*, **
Arbitary key-value pairs	n=2
Atomic electronic configuration	1s2.2s1, [Ar].3d4
Atomic term symbol	3P0, 2Po_1/2
Molecular term symbol	X(3$\Sigma$-g) C(1E"g)
Diatomic vibrational state	0, 1
Polyatomic vibrational state	v1+v2
Racah notation symbol	3d[3/2]_2
Energy, frequency or wavelength	$\lambda$=532 nm, E=12.4 eV

Equations

$$A{\left(\frac{T_e}{1\mathrm{eV}}\right)}^n\exp \left(-\frac{E}{T_e}\right)$$

$$A(\frac{T_g}{300}{)}^n\exp (\frac{-E}{T_g}),$$

Full text

††: Plasma Sources Sci. Technol.

1 Introduction

Realistic plasma models of many processes rest on the availability of reliable atomic and molecular data, so that the models are able to replicate the processes that drive the plasma at the submicroscopic level. Particularly for low-temperature plasmas, which are substantially molecular in composition, the set of possible processes, which we refer to as reactions below, can be very large. For low-temperature plasmas, accurate and comprehensive reaction datasets enable complex modelling of plasma-using technologies that empower our technology-based society [1]. Assembling appropriate datasets is therefore of critical importance.

For a given plasma composition, there are sets of species that are present in the plasma and a set of processes, generally called reactions, that will link the species or different states of the species. This reaction set is described as the “chemistry” for that plasma. For anything but the simplest molecular plasma, the number of possible reactions that could make up a chemistry can be very large [2]. The important reactions in a given plasma will be a subset of all these possible reactions, although it is not always possible to say in advance precisely which these reactions will be. In this context it is appropriate to characterise a useful chemistry as one which has three attributes: 1) The chemistry should be complete, that is contain all the important reactions for the given plasma. 2) It should be consistent, that is the reactions should not be unbalanced, thus resulting in the plasma composition being driven away from the true composition. 3) Finally, the plasma chemistry should be correct; this criterion cannot be demonstrated on theoretical grounds alone and requires validation against experimental measurements made in plasmas.

Assembling plasma chemistries is far from straightforward. While there may be several chemistries available for relatively simple systems such as molecular nitrogen plasmas [3, 4, 5, 6, 7], they generally do not exist for more complex problems such as the chemical mixtures typically used in etching and other technological plasmas. Indeed, given that reactions involving molecular radicals frequently remain completely uncharacterised [8], it is often a challenge to assemble a complete reaction set for these chemistries.

Here we present the Quantemol database (QDB). There are a growing number of databases aimed at supplying the needs of plasma modellers. For example, the recent LXCat project of Pitchford et al [9] aims to provide a web-based platform for data needed to model low-temperature plasmas. In practice LXCat considers electron collision processes but not heavy-particle (chemical) reactions. While both QDB and LXCat are set up to accept and provide multiple datasets for a single process if they are available, QDB aims to recommend a dataset for a particular application while LXCat leaves this choice to its users. The Phys4Entry database provides (ro)vibrationally resolved collisional data, including heterogeneous processes, for modeling re-entry plasmas [10]. For low-temperature, astronomical plasmas KIDA [11, 12] and BASECOL [13] provide data on chemical reactions and collisional excitation, respectively.

QDB aims to provide a repository for cross sections and/or rates for key reactions needed for models of low-temperature, i.e. molecular, plasmas. QDB collects data on both electron scattering and heavy-particle reactions and aims to facilitate and encourage peer-to-peer data sharing by its users. At present the data provided are largely for two-body reactions and hence are appropriate for low-pressure plasmas, but this will change in the future. Given sets of reactions, QDB then assembles these sets in chemistries for important plasma mixtures. If there are suitable experimental data available, these chemistries can be validated.

The following section gives an overview of QDB, with a technical specification of the data model given in the appendix. Section 3 catagorises the process types included in the database while Section 4 summarises the data sources used. A list of the reactions with a complete set of references is given as supplementary data to this article. Section 5 explains our chemistry construction and validation procedure; this is illustrated for two chemistries, those comprising SF6/CF4/O2 and SF6/CF4/N2/H2, respectively. These chemistries were selected due to their importance in silicon etching. Section 6 discusses future developments planned for the database, and the last section provides a summary and conclusions.

2 Overview

QDB provides reaction rates, cross sections and chemistries. The basic data item is the species, which can be state-specified, e.g. N(4S), or not, e.g. N2. At present QDB considers three generic species: the electron, the photon, and $M$, the third body in three-body reactions, plus 405 other atomic and molecular species. This total rises to 904 species when state-specified species are counted separately.

Species can undergo a series of processes called reactions. This includes, for example, elastic scattering or momentum transfer that are not generally regarded as chemical reactions. These processes are considered in more detail below. At present the database contains data on 4096 distinct reactions, comprising 2861 energy-dependent cross sections and 2259 temperature-dependent rate coefficients in Arrhenius form. Note that QDB allows multiple datasets for the same reaction, and for some reactions we have distinct data which are available as either cross sections or rates coefficients.

Many of these reactions are compiled into chemistries. Currently QDB has 28 chemistries, which are tabulated and discussed below. These chemistries can be validated, provided appropriate experimental data are available. Currently QDB contains 8 chemistries with some degree of validation; two of these chemistries are considered in detail below. Notes on the validation procedure are provided in the form of a datasheet for each validated chemistry. Chemistries are awarded a star rating that reflects how far they have been shown to satisfy the criteria of being complete, consistent, and correct. QDB is structured as a MySQL relational database; the data model used is discussed in the Appendix to this article.

Users can upload new data using the interface on the QDB website (www.quantemoldb.com) and download data using a choice of file formats, which is being expanded. Currently supported formats are comma-separated text for each reaction, provided as a zip file or in qdat format, which facilitates input for Kushner’s Hybrid Plasma Equipment Model (HPEM) [14, 15] code. The zip format contains all the cross sections, as individual comma-separated text files, and the rate coefficients that are needed for the specific chemistry. It also includes a manifest file listing all the files provided in the zip archive, a readme file, and a set of citations in bibtex format. The online view of each dataset contains data sheets and a form allowing users to provide feedback. Figure 1 shows a sample screen shot from QDB giving a comparison between various cross sections for electron-impact ionisation of methane. Note that for some reactions, cross sections for the process are obtained from two different sources; in these cases the agreement is excellent.

Development of QDB is performed with input from the Advisory Board whose members are co-authors of this paper.

3 Process types

Each reaction dataset in QDB is classified as containing cross sections or rate coefficients. The latter can be generated from the cross sections. The rate coefficients are expressed in Arrhenius form stored in the form of three parameters ($A$, $n$, and $E$), which can be used to compute the rate coefficients at the desired temperatures. For electron-impact reactions the Arrhenius formula is

[TABLE]

where $T_{e}$ is the electron temperature in eV and E is the activation energy in eV. For heavy-particle reactions the Arrhenius formula employed is

[TABLE]

where $T_{g}$ is the gas temperature in K and $E$ is the activation energy in K. In both cases, $A$ is the Arrhenius coefficient whose units depend on the order of the reactions. First-order reactions such as photodissociation and photoexcitation are expressed in s*-1*; second-order reactions, such as electron-impact reactions or two-body heavy-particle reactions, are expressed in cm3s*-1*; and three-body reactions use cm6s*-1*. Cross sections, e.g. for electron-neutral-molecule scattering, are given in units of cm2 as a function of electron energy in eV.

Each reaction is classified according to the process considered. These processes are listed in Tables 1, 2, and 3 which consider electron collision processes, heavy-particle reactions and processes involving photons, respectively. Note that some heavy-particle processes, in particular HAS and HIR, also involve a third body, generically denoted $M$ in the database. The process label does not depend on the presence of $M$ which is therefore not included in the process description.

4 Reactions in QDB

The scientific literature contains many measurements and calculations of reaction data that provide potentially useful input to plasma models. However, the task of extracting these data is far from straightforward. So far our strategy has been to focus on major data compilations and data sources. A list of those included so far is given in Table 4. In addition a variety of data was taken from models performed by Kushner and co-workers [18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35]. These sources were augmented with individual reactions taken directly from the original scientific literature. Where no suitable data could be found, internal (Quantemol) electron collision cross sections were generated using the Quantemol-N [36] implementation of the UK Molecular R-matrix code (UKRMol) [37]. As implied by Table 4, only a few of these cross sections have been published, although some have already been made available via LXCat [9] and the Virtual Atomic and Molecular Data Centre (VAMDC) [38]. The process of adding new data to QDB is a continuous one. The present results represent a snapshot of the situation as of November 2016.

A complete list of all reactions currently given in QDB with appropriate bibliographic references is provided in the supplementary data. Tables 5 and 6 summarise the sources of the data currently available in QDB by process type. At present there are relatively few radiative process in the database; the only ones involve radiative decay (PRD) in atoms [22, 101, 292].

5 Chemistry construction and validation

The chemistry sets are assembled starting from reactions already present in QDB; missing reactions are then extracted from the literature and added to QDB where possible. In cases where important reactions have not been previously studied, the missing reaction data are calculated using appropriate methods, such as Quantemol-N [36] for electron-molecule scattering reactions, or by scaling laws, or estimated to provide the necessary data. These data are added to QDB. This allows us to provide complete and self-consistent chemistry sets that form the starting point for validation. Figure 2 illustrates the network of 196 reactions assembled to characterise the chemistry of SF6 and O2. Table 7 lists the chemistries avaible in QDB as of November 2016.

The self-consistency of each chemistry set is checked using a range of models including Kushner’s zero-dimensional GlobalKin [293, 294, 30, 33] as implemented in Quantemol-P. GlobalKin couples molecular data (i.e. reaction probabilities) with plasma models to determine plasma properties, such as equilibrium concentrations of species. A variety of plasmas can be simulated using this software, such as etching and atmospheric pressure plasma reactors. HPEM, as implemented in Quantemol-VT, is also used.

Validation is achieved by using a chemistry set as the basis for the modelling of different industrial reactors. Comparison of the model output with measurement is the principal means by which validation is achieved. For higher-dimensional simulation, the behaviour of the species and the surface parameters across the wafer, such as etching or deposition rates, can also be used for comparison. Chemistry sets are given a reliability rating on the basis of these comparisons; see Table 8.

We now illustrate this process using two chemistry sets: those for SF6/CF4/O2 and SF6/CF4/N2/H2 gas mixture etching Si. As discussed below, these sets comprise subsets of mixtures, which provide important chemistries as well. These were also validated as part of the validation process. The chemistries were validated using GlobalKin and compared with experimental data provided by Infineon. Unfortunately the data available for complex reactions sets such as the ones considered is often very limited making the validation at best partial for these cases. We note that the primary goal of these examples is to illustrate the verification process in a simple form, not to reflect the highest verification rating. In fact, we consider the examples to only achieve the lowest level of validation – that is, agreement in the general trends between simulation and experiment for one discharge.

The Infineon tool consists of two parts. The first part is a coaxial microwave plasma discharge. In this part the electromagnetic wave propagates along the interface between the plasma column and the surrounding dielectric tube and the plasma column is sustained by electromagnetic energy. Free radicals are formed in this chamber with high efficiency. This chamber is connected to a larger vessel where the flux of particles propagates and where the remote wafer to be etched is located. Our GlobalKin models only attempted to model the second, larger chamber. GlobalKin performs a spatially homogeneous plasma chemistry simulation which are coupled with surface reaction modules. The model uses a Boltzmann solver to obtain electron impact reaction rate coefficients. These models assumed a plasma volume of 90 000 cm3, an area around the plasma of 10 700 cm2. The models, which used an assumed diffusion length of 8.3 cm, were initiated using the feedstock gases. They were run for 500 iterations, corresponding to a total of 1 s, which proved sufficient to reach steady state.

5.1 SF6/CF4/O2

Initially, distinct sets of chemistries for SF6/O2 and CF4/O2 were constructed and validated separately. These chemistries were then merged and missing reactions, such as SF6 + CF${}_{3}^{+}\rightarrow$ SF${}_{5}^{+}$ + CF4  [295] or SF${}_{x}^{-}$ + CF${}_{y}^{+}\rightarrow$ SFx + CFy, were identified and added. We then set up separate surface chemistries for SF6 and CF4, with a focus on silicon etching by F-radicals in the case of SF6 and CF4. Surface chemistry parameters were taken from Kokkoris et al. [296]. Since CF4 formed a smaller percentage of the mixture, see Table 9, we did not include the polymer deposition by CFx radicals. In addition to the F atom reactions, we added reactions of SFx radicals using the same reaction scheme as Kokkoris et al. [296], see also comments on this work by Nelson et al. [297].

Only just over 1% of oxygen is added to the mixture, as this increases the dissociation of SF6 and CF4. Due to its low density, the oxygen-related surface reactions were not included in the model. We also assumed that no significant concentration of CS molecules is formed during plasma processing with a mixture of SF6 and CF4, due to the large concentration of the F radicals in the mixture. This results in a much higher probability of formation of CxFy and SFx species. The conditions of the experiments are presented in Table 9.

To validate the chemistry set, we tested power variation between 1000 and 2000 W at a fixed pressure of 500 mTorr. As a second test we varied the pressure with power fixed at 2000 W. We used 40% of the experimental power in our simulation, in order to simulate the energy dissipation. This is almost certainly an over estimate of the power reaching the actual plasma in the Infineon device which explains why the simulations give faster etch rates than the measurements, see below.

Our GlobalKin simulations only provide global average values and therefore cannot provide an absolute quantitative comparison. For validation of our global model simulation we therefore compare trends. In particular, we compared the trends in the etching rate with measurements provided by Infineon Technologies, who studied the effect of both power variation and pressure variation. According to these measurements, the etch rate increases with increasing power but decreases with increasing pressure. Figure 3 illustrates the effect of varying power and pressure on the silicon etch rate for a mixture of SF6/CF4/O2. Good agreement in the trends is observed between the results of our simulation and the experimental data of Infineon. By increasing the power in the measurements and simulation we observe an increase in the Si etching rate but the rate drops as the pressure is raised.

Given the limited nature of the validation tests we have been able to perform for this chemistry we give this chemistry set a 3 star rating which is the lowest rating for a validated chemistry, see Table 8.

5.2 SF6/CF4/N2/H2

Initially, two sets of chemistries, for SF6/CF4, from the previous validation task, and for N2/H2, were constructed and validated separately. These chemistries were then merged. Since there is only a small proportion of N2/H2 in the mixture, we excluded species like NFx and CHxFy, as well as the reactions that lead to them. Rates for ion recombination (HMM) reactions such as ${\rm N}^{+}+{\rm F}^{-}\rightarrow{\rm N}+{\rm F}$ and ${\rm N}_{2}^{+}{\rm F}^{-}\rightarrow{\rm N}_{2}+{\rm F}$ were estimated from Moseley et al [272] and added to the chemistry list. The SF6/CF4 chemistry was generated using the same assumptions used for the SF6/CF4/O2 mixture. We also dealt with the surface reactions in a similar fashion to SF6/CF4/O2.

Table 10 summarises the experimental conditions assumed in the model. The experimental data from Infineon we tested against compared etch rates, under similar conditions, for the SF6/CF4/N2/H2 and SF6/CF4/O2 mixtures. According to the measurements, a higher etch rate is found for SF6/CF4/O2 than for SF6/CF4/N2/H2. As shown in Fig. 5, due to the lower dissociation in presence of N2/H2 compared with O2, both simulation and experimental results show a higher etch rate for the SF6/CF4/O2 mixtures compared to SF6/CF4/N2/H2 due to the lower production of F.

Given the limited nature of the validation tests we have been able to perform for this chemistry we give this chemistry set a 3 star rating which is the lowest rating for a validated chemistry, see Table 8.

6 Future developments

The process of adding both more reactions and chemistries to QDB is continuous and ongoing. We will also progressively improve the validation status of current chemistries and validate more chemistries, although these activities require appropriate experimental data to be available for us to validate against. We have developed an initial rating system for these chemistries, and we plan to develop this further by allowing users to also submit ratings. At the same time we plan to implement more formal uncertainty quantification (UQ) procedures [298].

The processes currently covered by QDB are listed in Tables 1, 2, and 3. These lists do cover all possible low-pressure gas phase processes. For example, inclusion of vibrationally-resolved reactions for molecules, such as in electron collisions with CO2 [81], is important for a number of plasma studies. QDB has the capability to hold such data but more work on data input and processes considered will be required to make it fully functional. At present QDB does not include processes that occur on surfaces and has only limited data for processes involving a third body. Both of these will be included in the database in the future. The inclusion of three-body reactions will extend the coverage to atmospheric pressure plasmas.

At present the data can be downloaded in two formats: a generic one and one that is appropriate for HPEM. With increasingly large datasets and sophisticated modelling programs, it is desirable for data to be transferred directly from the database to the model using an application program interface (API). We plan to develop APIs for commonly used plasma modelling programs in order to facilitate the use of QDB. We are also currenly implementing facilities for users to self-assemble chemistries in their own basket; in the longer term we plan to facilitate this with an automated chemistry generation tool.

7 Conclusions

One of the challenging problems when modelling plasmas is the lack of reliable chemistry data. For this purpose, we have developed the Quantemol Database (QDB), which aims to provide a platform for the exchange and validation of reactions that are important in plasmas and plasma chemistry datasets. The database provides data on both electron scattering and heavy-particle reactions, and it aims to facilitate and encourage peer-to-peer data sharing by its users. QDB currently includes almost 5000 reactions and 28 complete sets of chemistries; so far 8 of these sets have undergone some sort of validation. The set of reactions includes more than 2800 cross sections and more than 2100 sets of reaction rate coefficients in Arrhenius format for more than 980 species in different states.

Acknowledgements

We thank Prof. Uwe Czarnetzki, Dr. Nobuyuki Kuboi, Prof. Satoshi Hamaguchi, Prof. Mark Kushner and Prof. Leanne Pitchford for helpful discussions during this project. Partial funding was received from the UK STFC under grants ST/M007774 and ST/K004069, and from the Powerbase project. Powerbase has received funding from the Electronic Component Systems for European Leadership Joint Undertaking under grant agreement No 662133. This Joint Undertaking receives support from the European Union’s Horizon 2020 research and innovation programme and Austria, Belgium, Germany, Italy, Netherlands, Norway, Slovakia, Spain, United Kingdom.