Abstract: In this contribution, we present an overview of the various techniques for combining atomistic molecular dynamics with Monte Carlo simulations, mainly in the context of condensed matter systems, as well as a brief summary of the main accelerated dynamics techniques. Special attention is given to the force bias Monte Carlo technique and its combination with molecular dynamics, in view of promising recent developments, including a definable timescale. Various examples of the application of combined molecular dynamics / Monte Carlo simulations are given, in order to demonstrate the enhanced simulation efficiency with respect to either pure molecular dynamics or Monte Carlo.

Abstract: A hybrid model is used to investigate the fragmentation of C4F8 inductive discharges. Indeed, the resulting reactive species are crucial for the optimization of the Si-based etching process, since they determine the mechanisms of fluorination, polymerization, and sputtering. In this paper, we present the dissociation degree, the density ratio of F vs. CxFy (i.e., fluorocarbon (fc) neutrals), the neutral vs. positive ion density ratio, details on the neutral and ion components, and fractions of various fc neutrals (or ions) in the total fc neutral (or ion) density in a C4F8 inductively coupled plasma source, as well as the effect of pressure and power on these results. To analyze the fragmentation behavior, the electron density and temperature and electron energy probability function (EEPF) are investigated. Moreover, the main electron-impact generation sources for all considered neutrals and ions are determined from the complicated C4F8 reaction set used in the model. The C4F8 plasma fragmentation is explained, taking into account many factors, such as the EEPF characteristics, the dominance of primary and secondary processes, and the thresholds of dissociation and ionization. The simulation results are compared with experiments from literature, and reasonable agreement is obtained. Some discrepancies are observed, which can probably be attributed to the simplified polymer surface kinetics assumed in the model.

Abstract: We apply an implicit particle-in-cell Monte-Carlo (PIC-MC) method to study a radio-frequency argon microdischarge at steady state in the glow discharge limit, in which the microdischarge is sustained by secondary electron emission from the electrodes. The plasma density, electron energy distribution function (EEDF), and electron temperature are calculated in a wide range of operating conditions, including driving voltage, microdischarge gap, and pressure. Also, the effect of gap size scaling (in the range of 50-1000 μm) on the plasma sustaining voltage and peak electron density at atmospheric pressure is examined, which has not been explored before. In our simulations, three different EEDFs, i.e., a so-called three temperature hybrid mode, a two temperature α mode, and a two temperature γ mode distribution, are identified at different gaps and voltages. The maximum sustaining voltage to avoid a transition from the glow mode to an arc is predicted, as well as the minimum sustaining voltage for a steady glow discharge. Our calculations elucidate that secondary electrons play an essential role in sustaining the discharge, and as a result the relationship between breakdown voltage and gap spacing is far away from the Paschen law at atmospheric pressure.

Abstract: In this paper, the negative ion behavior in a C4F8 inductively coupled plasma (ICP) is investigated using a hybrid model. The model predicts a non-monotonic variation of the total negative ion density with power at low pressure (1030 mTorr), and this trend agrees well with experiments that were carried out in many fluorocarbon (fc) ICP sources, like C2F6, CHF3, and C4F8. This behavior is explained by the availability of feedstock C4F8 gas as a source of the negative ions, as well as by the presence of low energy electrons due to vibrational excitation at low power. The maximum of the negative ion density shifts to low power values upon decreasing pressure, because of the more pronounced depletion of C4F8 molecules, and at high pressure (∼50 mTorr), the anion density continuously increases with power, which is similar to fc CCP sources. Furthermore, the negative ion composition is identified in this paper. Our work demonstrates that for a clear understanding of the negative ion behavior in radio frequency C4F8 plasma sources, one needs to take into account many factors, like the attachment characteristics, the anion composition, the spatial profiles, and the reactor configuration. Finally, a detailed comparison of our simulation results with experiments is conducted.

Abstract: Laser-induced plasmas are transient systems rapidly aging in few nanoseconds of evolution. Time-of-flight spectrometry allowed studying initial plasma characteristics based on frozen translational degrees of freedom, hence overcoming intrinsic limitations of optical spectroscopy. Experimental ion velocity distributions were reconstructed as developed during the longitudinal plasma expansion. The obtained onset plasma temperatures are in the range of similar to 18-45 eV depending on the ablated metals. Also the ion angular spreads were found to be a function of ablated metal, e.g., the narrowest for Fe, the broadest for Al, due to different collisional coupling in the plasma population. (c) 2007 American Institute of Physics.

Abstract: We study the properties of a capacitive 13.56 MHz discharge with a mixture of Ar/C₂H₂ taking into account the plasmochemistry and growth of heavy hydrocarbons. A hybrid model was developed to combine the kinetic description for electron motion and the fluid approach for negative and positive ion transports and plasmochemical processes. A significant change in plasma parameters related to injection of 5.8% portion of acetylene in argon was observed and analyzed. We found that the electronegativity of the mixture is about 30%. The densities of negatively and positively charged heavy hydrocarbons are sufficiently large to be precursors for the formation of nanoparticles in the discharge volume.

Abstract: Molecular dynamics simulations using improved TersoffBrenner potential parameters were performed to investigate Cl+ etching of a {2×1} reconstructed Si(100) surface. Steady-state Si etching accompanying the Cl coverage of the surface is observed. Furthermore, a steady-state chlorinated reaction layer is formed. The thickness of this reaction layer is found to increase with increasing energy. The stoichiometry of SiClx species in the reaction layer is found to be SiCl:SiCl2:SiCl3 = 1.0:0.14:0.008 at 50 eV. These results are in excellent agreement with available experimental data. While elemental Si products are created by physical sputtering, most SiClx (0<x<4) etch products are produced by chemical-enhanced physical sputtering.

Abstract: Particle-in-cell/Monte Carlo collision simulations, coupled with an external circuit, are used to investigate the behavior of pulsed dual-frequency (DF) capacitively coupled plasmas (CCPs). It is found that the phase shift θ between the high (or low) frequency source and the pulse modulation has a great influence on the ion density and the ionization rate. By pulsing the high frequency source, the time-averaged ion density shows a maximum when θ = 90∘. The time-averaged ion energy distribution functions (IEDFs) at the driven electrode, however, keep almost unchanged, illustrating the potential of pulsed DF-CCP for independent control of ion density (and flux) and ion energy. A detailed investigation of the temporal evolution of the plasma characteristics indicates that several high frequency harmonics can be excited at the initial stage of a pulse period by tuning the phase shift θ, and this gives rise to strong sheath oscillations, and therefore high ionization rates. For comparison, the pulsing of the low frequency source is also studied. In this case, the ion density changes slightly as a function of time, and the time-averaged ion density shows the same trend as in the HF modulation for different phase shifts θ. Moreover, the time-averaged IEDFs at the driven electrode can be modulated, showing the potential to reduce the maximum ion bombardment energy.

Abstract: The present study focuses on the role of mass removal mechanisms in ns-laser ablation. A copper sample is placed in argon, initially set at standard pressure and temperature. Calculations are performed for a 6 ns laser pulse with a wavelength of 532 nm and laser fluences up to 10 J/cm2. The transient behavior in and above the copper target is described by a hydrodynamic model. Transmission profiles and ablation depths are compared with experimental results and similar trends are found. Our calculations reveal an interesting self-inhibiting mechanism: volumetric mass removal in the supercritical region triggers plasma shielding and therefore stops proceeding. This self-limiting process indicates that volumetric mass removal does not necessarily result in large ablation depths.