“Modeling and experimental investigation of the plasma uniformity in CF4/O2 capacitively coupled plasmas, operating in single frequency and dual frequency regime”. Zhang Y-R, Tinck S, De Schepper P, Wang Y-N, Bogaerts A, Journal of vacuum science and technology: A: vacuum surfaces and films 33, 021310 (2015). http://doi.org/10.1116/1.4906819
Abstract: A two-dimensional hybrid Monte Carlofluid model, incorporating a full-wave solution of Maxwell's equations, is employed to describe the behavior of high frequency (HF) and very high frequency capacitively coupled plasmas (CCPs), operating both at single frequency (SF) and dual frequency (DF) in a CF4/O2 gas mixture. First, the authors investigate the plasma composition, and the simulations reveal that besides CF4 and O2, also COF2, CF3, and CO2 are important neutral species, and CF+3 and F− are the most important positive and negative ions. Second, by comparing the results of the model with and without taking into account the electromagnetic effects for a SF CCP, it is clear that the electromagnetic effects are important, both at 27 and 60 MHz, because they affect the absolute values of the calculation results and also (to some extent) the spatial profiles, which accordingly affects the uniformity in plasma processing. In order to improve the plasma radial uniformity, which is important for the etch process, a low frequency (LF) source is added to the discharge. Therefore, in the major part of the paper, the plasma uniformity is investigated for both SF and DF CCPs, operating at a HF of 27 and 60 MHz and a LF of 2 MHz. For this purpose, the authors measure the etch rates as a function of position on the wafer in a wide range of LF powers, and the authors compare them with the calculated fluxes toward the wafer of the plasma species playing a role in the etch process, to explain the trends in the measured etch rate profiles. It is found that at a HF of 60 MHz, the uniformity of the etch rate is effectively improved by adding a LF power of 2 MHz and 300 W, while its absolute value increases by about 50%, thus a high etch rate with a uniform distribution is observed under this condition.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 1.374
Times cited: 3
DOI: 10.1116/1.4906819
|
“Numerical investigation of HBr/He transformer coupled plasmas used for silicon etching”. Gul B, Tinck S, De Schepper P, Aman-ur-Rehman, Bogaerts A, Journal of physics: D: applied physics 48, 025202 (2015). http://doi.org/10.1088/0022-3727/48/2/025202
Abstract: A two-dimensional hybrid Monte Carlofluid model is applied to study HBr/He inductively coupled plasmas used for etching of Si. Complete sets of gas-phase and surface reactions are presented and the effects of the gas mixing ratio on the plasma characteristics and on the etch rates are discussed. A comparison with experimentally measured etch rates is made to validate the modelling results. The etch rate in the HBr plasma is found to be quite low under the investigated conditions compared to typical etch rates of Si with F- or Cl-containing gases. This allows for a higher control and fine-tuning of the etch rate when creating ultra-small features. Our calculations predict a higher electron temperature at higher He fraction, because the electrons do not lose their energy so efficiently in vibrational and rotational excitations. As a consequence, electron impact ionization and dissociation become more important, yielding higher densities of ions, electrons and H atoms. This results in more pronounced sputtering of the surface. Nevertheless, the overall etch rate decreases upon increasing He fraction, suggesting that chemical etching is still the determining factor for the overall etch rate.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.588
Times cited: 7
DOI: 10.1088/0022-3727/48/2/025202
|
“Fluid simulation of the phase-shift effect in hydrogen capacitively coupled plasmas: 2 : radial uniformity of the plasma characteristics”. Zhang Y-R, Xu X, Bogaerts A, Wang Y-N, Journal of physics: D: applied physics 45, 015203 (2012). http://doi.org/10.1088/0022-3727/45/1/015203
Abstract: A two-dimensional fluid model, including the full set of Maxwell equations, has been developed and applied to investigate the effect of a phase shift between two power sources on the radial uniformity of several plasma characteristics in a hydrogen capacitively coupled plasma. This study was carried out at various frequencies in the range 13.56200 MHz. When the frequency is low, at 13.56 MHz, the plasma density is characterized by an off-axis peak when both power sources are in-phase (phgr = 0), and the best radial uniformity is obtained at phgr = π. This trend can be explained because the radial nonuniformity caused by the electrostatic edge effect can be effectively suppressed by the phase-shift effect at a phase difference equal to π. When the frequency rises to 60 MHz, the plasma density profiles shift smoothly from edge-peaked over uniform to centre-peaked as the phase difference increases, due to the pronounced standing-wave effect, and the best radial uniformity is reached at phgr = 0.3π. At a frequency of 100 MHz, a similar behaviour is observed, except that the maximum of the plasma density moves again towards the radial edge at the reverse-phase case (phgr = π), because of the dominant skin effect. When the frequency is 200 MHz, the bulk plasma density increases significantly with increasing phase-shift values, and a better uniformity is obtained at phgr = 0.4π. This is because the density in the centre increases faster than at the radial edge as the phase difference rises, due to the increasing power deposition Pz in the centre and the decreasing power density Pr at the radial edge. As the phase difference increases to π, the maximum near the radial edge becomes obvious again. This is because the skin effect has a predominant influence on the plasma density under this condition, resulting in a high density at the radial edge. Moreover, the axial ion flux increases monotonically with phase difference, and exhibits similar profiles to the plasma density. The calculation results illustrate that the radial uniformity of the various plasma characteristics is strongly dependent on the applied frequency and the phase shift between both power sources, which is important to realize, for controlling the uniformity of the plasma etch and deposition processes.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.588
Times cited: 15
DOI: 10.1088/0022-3727/45/1/015203
|
“Modeling ultrashort laser-induced emission from a negatively biased metal”. Wendelen W, Mueller BY, Autrique D, Bogaerts A, Rethfeld B, Applied physics letters 103, 221603 (2013). http://doi.org/10.1063/1.4830378
Abstract: A theoretical study of ultrashort laser-induced electron emission from a negatively biased metallic cathode has been performed. Classical as well as tunneling electron emission mechanisms are considered. It was found that electron emission is governed by an interplay of processes inside as well as above the cathode. A hybrid model is proposed, where the electron distribution within the target is retrieved from Boltzmann scattering integrals, while the charge distribution above it is studied by a Particle-In-Cell simulation. The results indicate that non-equilibrium effects determine the initial emission process, whereas the space charge above the target suppresses the effectively emitted charge.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 3.411
Times cited: 8
DOI: 10.1063/1.4830378
|
“Space charge corrected electron emission from an aluminum surface under non-equilibrium conditions”. Wendelen W, Mueller BY, Autrique D, Rethfeld B, Bogaerts A, Journal of applied physics 111, 113110 (2012). http://doi.org/10.1063/1.4729071
Abstract: A theoretical study has been conducted of ultrashort pulsed laser induced electron emission from an aluminum surface. Electron emission fluxes retrieved from the commonly employed Fowler-DuBridge theory were compared to fluxes based on a laser-induced non-equilibrium electron distribution. As a result, the two-and three-photon photoelectron emission parameters for the Fowler-DuBridge theory have been approximated. We observe that at regimes where photoemission is important, laser-induced electron emission evolves in a more smooth manner than predicted by the Fowler-DuBridge theory. The importance of the actual electron distribution decreases at higher laser fluences, whereas the contribution of thermionic emission increases. Furthermore, the influence of a space charge effect on electron emission was evaluated by a one dimensional particle-in-cell model. Depending on the fluences, the space charge reduces the electron emission by several orders of magnitude. The influence of the electron emission flux profiles on the effective electron emission was found to be negligible. However, a non-equilibrium electron velocity distribution increases the effective electron emission significantly. Our results show that it is essential to consider the non-equilibrium electron distribution as well as the space charge effect for the description of laser-induced photoemission. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4729071]
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.068
Times cited: 30
DOI: 10.1063/1.4729071
|
“Revisiting the interplay between ablation, collisional, and radiative processes during ns-laser ablation”. Autrique D, Gornushkin I, Alexiades V, Chen Z, Bogaerts A, Rethfeld B, Applied physics letters 103, 174102 (2013). http://doi.org/10.1063/1.4826505
Abstract: A study of ns-laser ablation is presented, which focuses on the transient behavior of the physical processes that act in and above a copper sample. A dimensionless multiphase collisional radiative model describes the interplay between the ablation, collisional, and radiative mechanisms. Calculations are done for a 6 ns-Nd:YAG laser pulse operating at 532 nm and fluences up to 15 J/cm2. Temporal intensity profiles as well as transmissivities are in good agreement with experimental results. It is found that volumetric ablation mechanisms and photo-processes both play an essential role in the onset of ns-laser induced breakdown.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 3.411
Times cited: 13
DOI: 10.1063/1.4826505
|
“Thermal conductivity of titanium nitride/titanium aluminum nitride multilayer coatings deposited by lateral rotating cathode arc”. Samani MK, Ding XZ, Khosravian N, Amin-Ahmadi B, Yi Y, Chen G, Neyts EC, Bogaerts A, Tay BK, Thin solid films : an international journal on the science and technology of thin and thick films 578, 133 (2015). http://doi.org/10.1016/j.tsf.2015.02.032
Abstract: A seriesof [TiN/TiAlN]nmultilayer coatingswith different bilayer numbers n=5, 10, 25, 50, and 100 were deposited on stainless steel substrate AISI 304 by a lateral rotating cathode arc technique in a flowing nitrogen atmosphere. The composition and microstructure of the coatings have been analyzed by using energy dispersive X-ray spectroscopy, X-ray diffraction (XRD), and conventional and high-resolution transmission electron microscopy (HRTEM). XRD analysis shows that the preferential orientation growth along the (111) direction is reduced in the multilayer coatings. TEM analysis reveals that the grain size of the coatings decreases with increasing bilayer number. HRTEMimaging of the multilayer coatings shows a high density misfit dislocation between the TiN and TiAlN layers. The cross-plane thermal conductivity of the coatings was measured by a pulsed photothermal reflectance technique. With increasing bilayer number, the multilayer coatings' thermal conductivity decreases gradually. This reduction of thermal conductivity can be ascribed to increased phonon scattering due to the disruption of columnar structure, reduced preferential orientation, decreased grain size of the coatings and present misfit dislocations at the interfaces.
Keywords: A1 Journal article; Electron microscopy for materials research (EMAT); Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 1.879
Times cited: 41
DOI: 10.1016/j.tsf.2015.02.032
|
“Rotating cylindrical magnetron sputtering: simulation of the reactive process”. Depla D, Li XY, Mahieu S, van Aeken K, Leroy WP, Haemers J, de Gryse R, Bogaerts A, Journal of applied physics 107, 113307 (2010). http://doi.org/10.1063/1.3415550
Abstract: A rotating cylindrical magnetron consists of a cylindrical tube, functioning as the cathode, which rotates around a stationary magnet assembly. In stationary mode, the cylindrical magnetron behaves similar to a planar magnetron with respect to the influence of reactive gas addition to the plasma. However, the transition from metallic mode to poisoned mode and vice versa depends on the rotation speed. An existing model has been modified to simulate the influence of target rotation on the well known hysteresis behavior during reactive magnetron sputtering. The model shows that the existing poisoning mechanisms, i.e., chemisorption, direct reactive ion implantation and knock on implantation, are insufficient to describe the poisoning behavior of the rotating target. A better description of the process is only possible by including the deposition of sputtered material on the target.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.068
Times cited: 15
DOI: 10.1063/1.3415550
|
“Development of a ReaxFF reactive force field for intrinsic point defects in titanium dioxide”. Huygh S, Bogaerts A, van Duin ACT, Neyts EC, Computational materials science 95, 579 (2014). http://doi.org/10.1016/j.commatsci.2014.07.056
Abstract: A reactive ReaxFF force field is developed for studying the influence of intrinsic point defects on the chemistry with TiO2 condensed phases. The force field parameters are optimized to ab initio data for the equations of state, relative phase stabilities for titanium and titanium dioxide, potential energy differences for (TiO2)n-clusters (n = 116). Also data for intrinsic point defects in anatase were added. These data contain formation energies for interstitial titanium and oxygen vacancies, diffusion barriers of the oxygen vacancies and molecular oxygen adsorption on a reduced anatase (101) surface. Employing the resulting force field, we study the influence of concentration of oxygen vacancies and expansion or compression of an anatase surface on the diffusion of the oxygen vacancies. Also the barrier for oxygen diffusion in the subsurface region is evaluated using this force field. This diffusion barrier of 27.7 kcal/mol indicates that the lateral redistribution of oxygen vacancies on the surface and in the subsurface will be dominated by their diffusion in the subsurface, since both this barrier as well as the barriers for diffusion from the surface to the subsurface and vice versa (17.07 kcal/mol and 21.91 kcal/mol, respectively, as calculated with DFT), are significantly lower than for diffusion on the surface (61.12 kcal/mol as calculated with DFT).
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.292
Times cited: 15
DOI: 10.1016/j.commatsci.2014.07.056
|
Bahnamiri OS, Verheyen C, Snyders R, Bogaerts A, Britun N (2021) Nitrogen fixation in pulsed microwave discharge studied by infrared absorption combined with modelling. 065007
Abstract: A pulsed microwave surfaguide discharge operating at 2.45 GHz was used for the conversion of molecular nitrogen into valuable compounds in several gas mixtures: N2 :O2 , N2 :O2 :CO2 and N2 :CO2 . The ro-vibrational absorption bands of the molecular species were monitored by a Fourier transform infrared apparatus in the post-discharge region in order to evaluate the relative number density of species, specifically NO production. The effects of specific energy input, pulse frequency, gas flow fraction, gas admixture and gas flow rate were studied for better understanding and optimization of the NO production yield and the corresponding energy cost (EC). By both the experiment and modelling, a highest NO yield is obtained at N2 :O2 (1:1) gas ratio in N2 :O2 mixture. The NO yield reveals a small growth followed by saturation when pulse repetition frequency increases. The energy efficiency start decreasing after the energy input reaches about 5 eV/molec, whereas the NO yield rises steadily at the same time. The lowest EC of about 8 MJ mol−1 corresponding to the yield and the energy efficiency of about 7% and 1% are found, respectively, in an optimum discharge condition in our case.
Keywords: A1 Journal Article;nitrogen fixation; pulsed microwave discharge; FTIR spectroscopy; discharge modelling; vibrational excitation; NO yield; energy cost; Plasma, laser ablation and surface modeling Antwerp (PLASMANT) ;
Impact Factor: 3.302
DOI: 10.1088/1361-6595/abff0e
|
“Investigating the plasma chemistry for the synthesis of carbon nanotubes/nanofibres in an inductively coupled plasma-enhanced CVD system : the effect of processing parameters”. Mao M, Bogaerts A, Journal of physics: D: applied physics 43, 315203 (2010). http://doi.org/10.1088/0022-3727/43/31/315203
Abstract: A parameter study is carried out for an inductively coupled plasma used for the synthesis of carbon nanotubes or carbon nanofibres (CNTs/CNFs), by means of the Hybrid Plasma Equipment Model. The influence of processing parameters including gas ratio for four different gas mixtures typically used for CNT/CNF growth (i.e. CH4/H2, CH4/NH3, C2H2/H2 and C2H2/NH3), inductively coupled plasma (ICP) power (501000 W), operating pressure (10 mTorr1 Torr), bias power (01000 W) and temperature of the substrate (01000 °C) on the plasma chemistry is investigated and the optimized conditions for CNT/CNF growth are analysed. Summarized, our calculations suggest that a lower fraction of hydrocarbon gases (CH4 or C2H2, i.e. below 20%) and hence a higher fraction of etchant gases (H2 or NH3) in the gas mixture result in more 'clean' conditions for controlled CNT/CNF growth. The same applies to a higher ICP power, a moderate ICP gas pressure above 100 mTorr (at least for single-walled carbon nanotubes), a high bias power (for aligned CNTs) and an intermediate substrate temperature.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.588
Times cited: 17
DOI: 10.1088/0022-3727/43/31/315203
|
“Fluid modelling of a packed bed dielectric barrier discharge plasma reactor”. Van Laer K, Bogaerts A, Plasma sources science and technology 25, 015002 (2016). http://doi.org/10.1088/0963-0252/25/1/015002
Abstract: A packed bed dielectric barrier discharge plasma reactor is computationally studied with a fluid model. Two different complementary axisymmetric 2D geometries are used to mimic the intrinsic 3D problem. It is found that a packing enhances the electric field strength and electron temperature at the contact points of the dielectric material due to polarization of the beads by the applied potential. As a result, these contact points prove to be of direct importance to initiate the plasma. At low applied potential, the discharge stays at the contact points, and shows the properties of a Townsend discharge. When a high enough potential is applied, the plasma will be able to travel through the gaps in between the beads from wall to wall, forming a kind of glow discharge. Therefore, the inclusion of a so-called ‘channel of voids’ is indispensable in any type of packed bed modelling.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 3.302
Times cited: 50
DOI: 10.1088/0963-0252/25/1/015002
|
“Influence of Gap Size and Dielectric Constant of the Packing Material on the Plasma Behaviour in a Packed Bed DBD Reactor: A Fluid Modelling Study: Influence of Gap Size and Dielectric Constant…”. Van Laer K, Bogaerts A, Plasma processes and polymers 14, 1600129 (2017). http://doi.org/10.1002/ppap.201600129
Abstract: A packed bed dielectric barrier discharge (DBD) was studied by means of fluid modelling, to investigate the influence of the dielectric constant of the packing on the plasma characteristics, for two different gap sizes. The electric field strength and electron temperature are much more enhanced in a microgap reactor than
in a mm-gap reactor, leading to more current peaks per half-cycle, but also to non-quasineutral plasma. Increasing the dielectric constant enhances the electric field further, but only up to a certain value of dielectric constant, being 9 for a microgap and 100 for a mm-gap reactor. The enhanced electric field results in a higher electron temperature, but also lower electron density. This last one strongly affects the reaction rate.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.846
Times cited: 23
DOI: 10.1002/ppap.201600129
|
“Streamer propagation in a packed bed plasma reactor for plasma catalysis applications”. Wang W, Kim H-H, Van Laer K, Bogaerts A, Chemical engineering journal 334, 2467 (2018). http://doi.org/10.1016/j.cej.2017.11.139
Abstract: A packed bed dielectric barrier discharge (DBD) is widely used for plasma catalysis applications, but the exact plasma characteristics in between the packing beads are far from understood. Therefore, we study here these plasma characteristics by means of fluid modelling and experimental observations using ICCD imaging, for packing materials with different dielectric constants. Our study reveals that a packed bed DBD reactor in dry air at atmospheric pressure may show three types of discharges, i.e. positive restrikes, filamentary microdischarges, which can also be localized between two packing beads, and surface discharges (so-called surface ionization
waves). Restrikes between the dielectric surfaces result in the formation of filamentary microdischarges, while surface charging creates electric field components parallel to the dielectric surfaces, leading to the formation of surface discharges. A transition in discharge mode occurs from surface discharges to local filamentary discharges between the packing beads when the dielectric constant of the packing rises from 5 to 1000. This may have implications for the efficiency of plasma catalytic gas treatment, because the catalyst activation may be limited by constraining the discharge to the contact points of the beads. The production of reactive species occurs most in the positive restrikes, the surface discharges and the local microdischarges in between the beads, and is less significant in the longer filamentary microdischarges. The faster streamer propagation and discharge development with higher dielectric constant of the packing beads leads to a faster production of reactive species. This study is of great interest for plasma catalysis, where packing beads with different dielectric constants are often used as supports for the catalytic materials. It allows us to better understand how different packing materials can influence the performance of packed bed plasma reactors for environmental applications.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 6.216
Times cited: 36
DOI: 10.1016/j.cej.2017.11.139
|
“Kinetic simulation of direct-current driven microdischarges in argon at atmospheric pressure”. Zhang Y, Jiang W, Bogaerts A, Journal of physics: D: applied physics 47, 435201 (2014). http://doi.org/10.1088/0022-3727/47/43/435201
Abstract: A one-dimensional, implicit particle-in-cell Monte Carlo collision model is used to simulate the plasma kinetic properties at a steady state in a parallel-plate direct current argon glow microdischarge under various operating conditions, such as driving voltage (301000 V) and gap size (101000 µm) at atmospheric pressure. First, a comparison between rf and dc modes is shown for the same pressure, driving voltage and gap spacing. Furthermore, the effect of gap size scaling (in the range of 101000 µm) on the breakdown voltage, peak electron density and peak electron current density at the breakdown voltage is examined. The breakdown voltage is lower than 150 V in all gaps considered. The microdischarge is found to have a neutral bulk plasma region and a cathode sheath region with size varying with the applied voltage and the discharge gap. In our calculations, the electron and ion densities are of the order of 10181023 m−3, which is in the glow discharge limit, as the ionization degree is lower than 1% . The electron energy distribution function shows a two-energy group distribution at a gap of 10 µm and a three-energy group distribution at larger gaps such as 200 µm and 1000 µm, emphasizing the importance of the gap spacing in dc microdischarges.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.588
Times cited: 10
DOI: 10.1088/0022-3727/47/43/435201
|
“Fluid modeling of the conversion of methane into higher hydrocarbons in an atmospheric pressure dielectric barrier discharge”. De Bie C, Verheyde B, Martens T, van Dijk J, Paulussen S, Bogaerts A, Plasma processes and polymers 8, 1033 (2011). http://doi.org/10.1002/ppap.201100027
Abstract: A one-dimensional fluid model for a dielectric barrier discharge in methane, used as a chemical reactor for gas conversion, is developed. The model describes the gas phase chemistry governing the conversion process of methane to higher hydrocarbons. The spatially averaged densities of the various plasma species as a function of time are discussed. Besides, the conversion of methane and the yields of the reaction products as a function of the residence time in the reactor are shown and compared with experimental data. Higher hydrocarbons (C2Hy and C3Hy) and hydrogen gas are typically found to be important reaction products. Furthermore, the main underlying reaction pathways are determined.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.846
Times cited: 70
DOI: 10.1002/ppap.201100027
|
“The Dominant Pathways for the Conversion of Methane into Oxygenates and Syngas in an Atmospheric Pressure Dielectric Barrier Discharge”. De Bie C, van Dijk J, Bogaerts A, The journal of physical chemistry: C : nanomaterials and interfaces 119, 22331 (2015). http://doi.org/10.1021/acs.jpcc.5b06515
Abstract: A one-dimensional fluid model for a dielectric barrier discharge in CH4/O2 and CH4/CO2 gas mixtures is developed. The model describes the gas-phase chemistry for partial oxidation and for dry reforming of methane. The spatially averaged densities of the various plasma species are presented as a function of time and initial gas mixing ratio. Besides, the conversion of the inlet gases and the selectivities of the reaction products are calculated. Syngas, higher hydrocarbons, and higher oxygenates are typically found to be important reaction products. Furthermore, the main underlying reaction pathways for the formation of syngas, methanol, formaldehyde, and other higher oxygenates are determined.
Keywords: A1 Journal article; Engineering sciences. Technology; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 4.536
Times cited: 46
DOI: 10.1021/acs.jpcc.5b06515
|
“Fluid modelling of an atmospheric pressure dielectric barrier discharge in cylindrical geometry”. Petrović, D, Martens T, van Dijk J, Brok WJM, Bogaerts A, Journal of physics: D: applied physics 42, 205206 (2009). http://doi.org/10.1088/0022-3727/42/20/205206
Abstract: A numerical parameter study has been performed for a cylindrical atmospheric pressure dielectric barrier discharge (DBD) in helium with nitrogen impurities using a two-dimensional time-dependent fluid model. The calculated electric currents and gap voltages as a function of time for a given applied potential are presented, as well as the number densities of the various plasma species. This study shows that for the geometry under consideration the applied voltage parameters have a large impact on the electric current profiles and that the discharge current is always determined by the electron and ion conduction currents while the displacement current is nearly negligible. A relative broadening of the current profiles (compared with the duration of the half cycle of the applied voltage) with an increase in the applied frequency is obtained. Nearly sinusoidal current wave forms, usually typical for radio frequency DBDs, are observed while still operating at the frequencies of tens of kilohertz. For the setup under investigation, the Townsend mode of the DBD is observed in the entire range of applied voltage amplitudes and frequencies. It is shown that the average power density dissipated in the discharge increases with rising applied voltage and frequency. An increase in applied voltage frequency leads to an increase in the electron density and a decrease in electron energy, while increasing the voltage amplitude has the opposite effect.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.588
Times cited: 29
DOI: 10.1088/0022-3727/42/20/205206
|
“The effect of O2 in a humid O2/N2/NOx gas mixture on NOx and N2O remediation by an atmospheric pressure dielectric barrier discharge”. Teodoru S, Kusano Y, Bogaerts A, Plasma processes and polymers 9, 652 (2012). http://doi.org/10.1002/ppap.201100187
Abstract: A numerical model for NxOy remediation in humid air plasma produced with a dielectric barrier discharge at atmospheric pressure is presented. Special emphasis is given to NO2 and N2O reduction with the decrease of O2 content in the feedstock gas. A detailed reaction mechanism including electronic and ionic processes, as well as the contribution of radicals and excited atomic/molecular species is proposed. The temporal evolution of the densities of NO, NO2 and N2O species, and some other by-products, is analyzed, and the major pathways for the NxOy remediation are discussed for one pulse. Subsequently, simulations are presented for a multi-pulses case, where three O2 contents are tested for optimization of the remediation process. It is found that when the gas mixture O2/N2/H2O/NOx has no initial O2 content, the best NOx and N2O remediation is achieved.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.846
Times cited: 24
DOI: 10.1002/ppap.201100187
|
“Sputter deposition of MgxAlyOz thin films in a dual-magnetron device : a multi-species Monte Carlo model”. Yusupov M, Saraiva M, Depla D, Bogaerts A, New journal of physics 14, 073043 (2012). http://doi.org/10.1088/1367-2630/14/7/073043
Abstract: A multi-species Monte Carlo (MC) model, combined with an analytical surface model, has been developed in order to investigate the general plasma processes occurring during the sputter deposition of complex oxide films in a dual-magnetron sputter deposition system. The important plasma species, such as electrons, Ar+ ions, fast Ar atoms and sputtered metal atoms (i.e. Mg and Al atoms) are described with the so-called multi-species MC model, whereas the deposition of MgxAlyOz films is treated by an analytical surface model. Targetsubstrate distances for both magnetrons in the dual-magnetron setup are varied for the purpose of growing stoichiometric complex oxide thin films. The metal atoms are sputtered from pure metallic targets, whereas the oxygen flux is only directed toward the substrate and is high enough to obtain fully oxidized thin films but low enough to avoid target poisoning. The calculations correspond to typical experimental conditions applied to grow these complex oxide films. In this paper, some calculation results are shown, such as the densities of various plasma species, their fluxes toward the targets and substrate, the deposition rates, as well as the film stoichiometry. Moreover, some results of the combined model are compared with experimental observations. Note that this is the first complete model, which can be applied for large and complicated magnetron reactor geometries, such as dual-magnetron configurations. With this model, we are able to describe all important plasma species as well as the deposition process. It can also be used to predict film stoichiometries of complex oxide films on the substrate.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 3.786
Times cited: 2
DOI: 10.1088/1367-2630/14/7/073043
|
“Monte Carlo analysis of the electron thermalization process in the afterglow of a microsecond dc pulsed glow discharge”. Martín A, Bordel N, Pereiro R, Bogaerts A, Spectrochimica acta: part B : atomic spectroscopy 63, 1274 (2008). http://doi.org/10.1016/j.sab.2008.09.012
Abstract: A Monte Carlo model is utilized for studying the behavior of electrons in the afterglow of an analytical microsecond dc pulsed glow discharge. This model uses several quantities as input data, such as electric field and potential, ion flux at the cathode, the fast argon ion and atom impact ionization rates, slow electron density, the electrical characterization of the pulse (voltage and current profiles) and temperature profile. These quantities were obtained by earlier Monte Carlo fluid calculations for a pulsed discharge. Our goal is to study the behavior of the so-called Monte Carlo electrons (i.e., those electrons created at the cathode or by ionization collisions in the plasma which are followed by using the Monte Carlo model) from their origin to the moment when they are absorbed at the cell walls or when they have lost their energy by collisions (being transferred to the group of slow electrons) in the afterglow of the pulsed discharge. The thermalization of the electrons is a phenomenon where the electron-electron Coulomb collisions acquire a special importance. Indeed, in the afterglow the cross sections of the other electron reactions taken into account in the model are very low, because of the very low electron energy. We study the electron energy distributions at several times during and after the pulse and at several positions in the plasma cell, focusing on the thermalization and on the behavior of the electrons in the afterglow. Also, the time evolution of the rates of the various collision processes, the average electron energy, the densities of Monte Carlo and slow electrons and the ionization degree are investigated.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 3.241
Times cited: 9
DOI: 10.1016/j.sab.2008.09.012
|
“Behavior of electrons in a dual-magnetron sputter deposition system : a Monte Carlo model”. Yusupov M, Bultinck E, Depla D, Bogaerts A, New journal of physics 13, 033018 (2011). http://doi.org/10.1088/1367-2630/13/3/033018
Abstract: A Monte Carlo model has been developed for investigating the electron behavior in a dual-magnetron sputter deposition system. To describe the three-dimensional (3D) geometry, different reference frames, i.e. a local and a global coordinate system, were used. In this study, the influence of both closed and mirror magnetic field configurations on the plasma properties is investigated. In the case of a closed magnetic field configuration, the calculated electron trajectories show that if an electron is emitted in (or near) the center of the cathode, where the influence of the magnetic field is low, it is able to travel from one magnetron to the other. On the other hand, when an electron is created at the race track area, it is more or less trapped in the strong magnetic field and cannot easily escape to the second magnetron region. In the case of a mirror magnetic field configuration, irrespective of where the electron is emitted from the cathode, it cannot travel from one magnetron to the other because the magnetic field lines guide the electron to the substrate. Moreover, the electron density and electron impact ionization rate have been calculated and studied in detail for both configurations.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 3.786
Times cited: 12
DOI: 10.1088/1367-2630/13/3/033018
|
“Molecular dynamics simulation of oxide thin film growth: importance of the inter-atomic interaction potential”. Georgieva V, Todorov IT, Bogaerts A, Chemical physics letters 485, 315 (2010). http://doi.org/10.1016/j.cplett.2009.12.067
Abstract: A molecular dynamics (MD) study of MgxAlyOz thin films grown by magnetron sputtering is presented using an ionic model and comparing two potential sets with formal and partial charges. The applicability of the model and the reliability of the potential sets for the simulation of thin film growth are discussed. The formal charge potential set was found to reproduce the thin film structure in close agreement with the structure of the experimentally grown thin films. Graphical abstract A molecular dynamics study of growth of MgxAlyOz thin films is presented using an ionic model and comparing two potential sets with formal and partial charges. The simulation results with the formal charge potential set showed a transition in the film from a crystalline to an amorphous structure, when the Mg metal content decreases below 50% in very close agreement with the structure of the experimentally deposited films.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 1.815
Times cited: 16
DOI: 10.1016/j.cplett.2009.12.067
|
“Computational study of the CF4 /CHF3 / H2 /Cl2 /O2 /HBr gas phase plasma chemistry”. Tinck S, Bogaerts A, Journal of physics: D: applied physics 49, 195203 (2016). http://doi.org/10.1088/0022-3727/49/19/195203
Abstract: A modelling study is performed of high-density low-pressure inductively coupled CF4/CHF3/H2/Cl2/O2/HBr plasmas under different gas mixing ratios. A reaction set describing the complete plasma chemistry is presented and discussed. The gas fraction of each component in this mixture is varied to investigate the sensitivity of the plasma properties, like electron density, plasma potential and species densities, towards the gas mixing ratios. This research is of great interest for microelectronics applications because these gases are often combined in two (or more)-component mixtures, and mixing gases or changing the fraction of a gas can sometimes yield unwanted reaction products or unexpected changes in the overall plasma properties due to the increased chemical complexity of the system. Increasing the CF4 fraction produces more F atoms for chemical etching as expected, but also more prominently lowers the density of Cl atoms, resulting in an actual drop in the etch rate under certain conditions. Furthermore, CF4 decreases the free electron density when mixed with Cl2. However, depending on the other gas components, CF4 gas can also sometimes enhance free electron density. This is the case when HBr is added to the mixture. The addition of H2 to the gas mixture will lower the sputtering process, not only due to the lower overall positive ion density at higher H2 fractions, but also because more H+, H2 + and H3 + are present and they have very low sputter yields. In contrast, a larger Cl2 fraction results in more chemical etching but also in less physical sputtering due to a smaller abundance of positive ions. Increasing the O2 fraction in the plasma will always lower the etch rate due to more oxidation of the wafer surface and due to a lower plasma density. However, it is also observed that the density of F atoms can actually increase with rising O2 gas fraction. This is relevant to note because the exact balance between fluorination and oxidation is important for fine-tuning the overall etch rate and for control of the sidewall profile. Finally, HBr is often used as a chemical etcher, but when mixed with F- or Cl-containing gases, HBr creates the same diluting effects as Ar or He, because a
higher fraction results in less chemical etching but more (physical) sputtering.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.588
Times cited: 5
DOI: 10.1088/0022-3727/49/19/195203
|
“Multi-element model for the simulation of inductively coupled plasmas : effects of helium addition to the central gas stream”. Lindner H, Bogaerts A, Spectrochimica acta: part B : atomic spectroscopy 66, 421 (2011). http://doi.org/10.1016/j.sab.2011.04.007
Abstract: A model for an atmospheric pressure inductively coupled plasma (ICP) is developed which allows rather easy extension to a variable number of species and ionisation degrees. This encompasses an easy calculation of transport parameters for mixtures, ionisation and heat capacity. The ICP is modeled in an axisymmetric geometry, taking into account the gas streaming into a flowing ambient gas. A mixture of argon and helium is applied in the injector gas stream as it is often done in laser ablation ICP spectrometry. The results show a strong influence of the added helium on the center of the ICP, which is important for chemical analysis. The length of the central channel is significantly increased and the temperature inside is significantly higher than in the case of pure argon. This means that higher gas volume flow rates can be applied by addition of helium compared to the use of pure argon. This has the advantage that the gas velocity in the transport system towards the ICP can be increased, which allows shorter washout-times. Consequently, shorter measurement times can be achieved, e.g. for spatial mapping analyses in laser ablation ICP spectrometry. Furthermore, the higher temperature and the longer effective plasma length will increase the maximum size of droplets or particles injected into the ICP that are completely evaporated at the detection site. Thus, we expect an increase of the analytical performance of the ICP by helium addition to the injector gas.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 3.241
Times cited: 28
DOI: 10.1016/j.sab.2011.04.007
|
“Elucidating the asymmetric behavior of the discharge in a dual magnetron sputter deposition system”. Yusupov M, Bultinck E, Depla D, Bogaerts A, Applied physics letters 98, 131502 (2011). http://doi.org/10.1063/1.3574365
Abstract: A magnetron discharge is characterized by drifts of the charged particles guiding center, caused by the magnetic field, in contrast to unmagnetized discharges. Because of these drifts, a pronounced asymmetry of the discharge can be observed in a dual magnetron setup. In this work, it is found that the shape of the discharge in a dual magnetron configuration depends on the magnetic field configuration. In a closed configuration, strong drifts were observed in one preferential direction, whereas in a mirror configuration the deflection of the discharge was not so pronounced. Our calculations confirm experimental observations.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 3.411
Times cited: 4
DOI: 10.1063/1.3574365
|
“Effect of small amounts of hydrogen added to argon glow discharges: hybrid Monte-Carlo-fluid model”. Bogaerts A, Gijbels R, Physical review : E : statistical, nonlinear, and soft matter physics 65, 056402 (2002). http://doi.org/10.1103/PhysRevE.65.056402
Abstract: A hybrid Monte Carlofluid modeling network is developed for an argon-hydrogen mixture, to predict the effect of small amounts of hydrogen added to a dc argon glow discharge. The species considered in the model include the Ar gas atoms, electrons, Ar+ ions and fast Ar atoms, ArH+, H+, H+2 and H+3 ions, and H atoms and H2 molecules, as well as Ar metastable atoms, sputtered Cu atoms, and the corresponding Cu+ ions. Sixty-five reactions between these species are incorporated in the model. The effect of hydrogen on various calculation results is investigated, such as the species densities, the relative role of different production and loss processes for the various species, the cathode sputtering rate and contributions by different bombarding species, and the dissociation degree of H2 and the ionization degree of Ar and Cu. The calculation results are presented and discussed for 1% H2 addition, and comparison is also made with a pure argon discharge and with only 0.1% H2 addition.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.366
Times cited: 33
DOI: 10.1103/PhysRevE.65.056402
|
“Cryogenic etching of silicon with SF6 inductively coupled plasmas: a combined modelling and experimental study”. Tinck S, Tillocher T, Dussart R, Bogaerts A, Journal of physics: D: applied physics 48, 155204 (2015). http://doi.org/10.1088/0022-3727/48/15/155204
Abstract: A hybrid Monte Carlofluid model is applied to simulate the wafer-temperature-dependent etching of silicon with SF6 inductively coupled plasmas (ICP). The bulk plasma within the ICP reactor volume as well as the surface reactions occurring at the wafer are self-consistently described. The calculated etch rates are validated by experiments. The calculations and experiments are performed at two different wafer temperatures, i.e. 300 and 173 K, resembling conventional etching and cryoetching, respectively. In the case of cryoetching, a physisorbed SFx layer (x = 06) is formed on the wafer, which is negligible at room temperature, because of fast thermal desorption, However, even in the case of cryoetching, this layer can easily be disintegrated by low-energy ions, so it does not affect the etch rates. In the investigated pressure range of 19 Pa, the etch rate is always slightly higher at cryogenic conditions, both in the experiments and in the model, and this could be explained in the model due to a local cooling of the gas above the wafer, making the gas denser and increasing the flux of reactive neutrals, like F and F2, towards the wafer.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.588
Times cited: 9
DOI: 10.1088/0022-3727/48/15/155204
|
“Numerical study of the plasma chemistry in inductively coupled SF6 and SF6/AR plasmas used for deep silicon etching applications”. Mao M, Wang YN, Bogaerts A, Journal of physics: D: applied physics 44, 435202 (2011). http://doi.org/10.1088/0022-3727/44/43/435202
Abstract: A hybrid model, called the hybrid plasma equipment model, was used to study inductively coupled SF6 plasmas used for Si etching applications. The plasma properties such as number densities of electrons, positive and negative ions, and neutrals are calculated under typical etching conditions. The electron kinetics is analysed by means of the electron energy probability function. The plasma chemistry taking place in pure SF6 and in an Ar/SF6 mixture is also discussed, and finally the effect of the argon fraction on the plasma properties is investigated.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.588
Times cited: 20
DOI: 10.1088/0022-3727/44/43/435202
|
“Simulation of an Ar/Cl2 inductively coupled plasma: study of the effect of bias, power and pressure and comparison with experiments”. Tinck S, Boullart W, Bogaerts A, Journal of physics: D: applied physics 41, 065207 (2008). http://doi.org/10.1088/0022-3727/41/6/065207
Abstract: A hybrid model, called the hybrid plasma equipment model, was used to study Ar/Cl(2) inductively coupled plasmas used for the etching of Si. The effects of substrate bias, source power and gas pressure on the plasma characteristics and on the fluxes and energies of plasma species bombarding the substrate were observed. A comparison with experimentally measured etch rates was made to investigate how the etch process is influenced and which plasma species mainly account for the etch process. First, the general plasma characteristics are investigated at the following operating conditions: 10% Ar 90% Cl(2) gas mixture, 5mTorr total gas pressure, 100 sccm gas flow rate, 250W source power, -200V dc bias at the substrate electrode and an operating frequency of 13.56MHz applied to the coil and to the substrate electrode. Subsequently, the pressure is varied from 5 to 80mTorr, the substrate bias from -100 to -300V and the source power from 250 to 1000W. Increasing the total gas pressure results in a decrease of the etch rate and a less anisotropic flux to the substrate due to more collisions of the ions in the sheath. Increasing the substrate bias has an effect on the energy of the ions bombarding the substrate and to a lesser extent on the magnitude of the ion flux. When source power is increased, it was found that, not the energy, but the magnitude of the ion flux is increased. The etch rate was more influenced by a variation of the substrate bias than by a variation of the source power, at these operating conditions. These results suggest that the etch process is mainly affected by the energy of the ions bombarding the substrate and the magnitude of the ion flux, and to a lesser extent by the magnitude of the radical flux.
Keywords: A1 Journal article; Plasma Lab for Applications in Sustainability and Medicine – Antwerp (PLASMANT)
Impact Factor: 2.588
Times cited: 31
DOI: 10.1088/0022-3727/41/6/065207
|