Abstract: Atomistic simulation methods for the quantification of free energies are in wide use. These methods operate by sampling the probability density of a system along a small set of suitable collective variables (CVs), which is, in turn, expressed in the form of a free energy surface (FES). This definition of the FES can capture the relative stability of metastable states but not that of the transition state because the barrier height is not invariant to the choice of CVs. Free energy barriers therefore cannot be consistently computed from the FES. Here, we present a simple approach to calculate the gauge correction necessary to eliminate this inconsistency. Using our procedure, the standard FES as well as its gauge-corrected counterpart can be obtained by reweighing the same simulated trajectory at little additional cost. We apply the method to a number of systems—a particle solvated in a Lennard-Jones fluid, a Diels–Alder reaction, and crystallization of liquid sodium—to demonstrate its ability to produce consistent free energy barriers that correctly capture the kinetics of chemical or physical transformations, and discuss the additional demands it puts on the chosen CVs. Because the FES can be converged at relatively short (sub-ns) time scales, a free energy-based description of reaction kinetics is a particularly attractive option to study chemical processes at more expensive quantum mechanical levels of theory.

Abstract: The evaluation of nucleation rates from molecular dynamics trajectories is hampered by the slow nucleation time scale and impact of finite size effects. Here, we show that accurate nucleation rates can be obtained in a very general fashion relying only on the free energy barrier, transition state theory, and a simple dynamical correction for diffusive recrossing. In this setup, the time scale problem is overcome by using enhanced sampling methods, in casu metadynamics, whereas the impact of finite size effects can be naturally circumvented by reconstructing the free energy surface from an appropriate ensemble. Approximations from classical nucleation theory are avoided. We demonstrate the accuracy of the approach by calculating macroscopic rates of droplet nucleation from argon vapor, spanning 16 orders of magnitude and in excellent agreement with literature results, all from simulations of very small (512 atom) systems.

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: Computer simulations based on the particle-in-cell/Monte Carlo collision method are performed to study the plasma characteristics and especially the transition in electron heating mechanisms in a hybrid direct current (dc)/dual-frequency (DF) capacitively coupled CF 4 discharge. When applying a superposed dc voltage, the plasma density first increases, then decreases, and finally increases again, which is in good agreement with experiments. This trend can be explained by the transition between the four main heating modes, i.e., DF coupling, dc and DF coupling, dc source dominant heating, and secondary electron dominant heating.

Abstract: We report on the thin film resistivity of several platinum-group metals (Ru, Pd, Ir, and Pt). Platinum-group thin films show comparable or lower resistivities than Cu for film thicknesses below about 5 nm due to a weaker thickness dependence of the resistivity. Based on experimentally determined mean linear distances between grain boundaries as well as ab initio calculations of the electron mean free path, the data for Ru, Ir, and Cu were modeled within the semiclassical Mayadas-Shatzkes model [Phys. Rev. B 1, 1382 (1970)] to assess the combined contributions of surface and grain boundary scattering to the resistivity. For Ru, the modeling results indicated that surface scattering was strongly dependent on the surrounding material with nearly specular scattering at interfaces with SiO2 or air but with diffuse scattering at interfaces with TaN. The dependence of the thin film resistivity on the mean free path is also discussed within the Mayadas-Shatzkes model in consideration of the experimental findings. Published by AIP Publishing.

Abstract: The generation of high-density plasmas on the surface of porous catalysts is very important for plasma catalysis, as it determines the active surface of the catalyst that is available for the reaction. In this work, we investigate the mechanism of surface and volume plasma streamer formation and propagation near micro-sized pores in dielectric barrier discharges operating in air at atmospheric pressure. A two-dimensional particle-in-cell/ Monte Carlo collision model is used to model the individual kinetic behavior of plasma species. Our calculations indicate that the surface discharge is enhanced on the surface of the catalyst pores compared with the microdischarge inside the catalyst pores. The reason is that the surface ionization wave induces surface charging along the catalyst pore sidewalls, leading to a strong electric field along the pore sidewalls, which in turn further enhances the surface discharge. Therefore, highly concentrated reactive species occur on the surfaces of the catalyst pores, indicating high-density plasmas on the surface of porous catalysts. Indeed, the maximum electron impact excitation and ionization rates occur on the pore surface, indicating the more pronounced production of excited state and electron-ion pairs on the pore surface than inside the pore, which may profoundly affect the plasma catalytic process. Published under license by AIP Publishing.

Abstract: We examined the vacancy trapping proficiency of Sn and P atoms in germanium using positron annihilation spectroscopy measurements, sensitive to the open-volume defects. Epitaxial Ge1 xSnx films were grown by chemical vapor deposition with different P concentrations in the 3: 0 1019-1: 5 1020 cm 3 range. We corroborate our findings with first principles simulations. Codoping of Ge with a Sn concentration of up to 9% is not an efficient method to suppress the free vacancy concentration and the formation of larger phosphorus-vacancy complexes. Experimental results confirm an increase in the number of P atoms around the monovacancy with P-doping, leading to dopant deactivation in epitaxial germanium-tin layers with similar Sn content. Vice versa, no impact on the improvement of maximum achieved P activation in Ge with increasing Sn-doping has been observed. Theoretical calculations also confirm that Pn-V (vacancy) complexes are energetically more stable than the corresponding SnmPn-V and Snm-V defect structures with the same number of alien atoms (Sn or P) around the monovacancy. he strong attraction of vacancies to the phosphorus atoms remains the dominant dopant deactivation mechanism in Ge as well as in Ge1 xSnx. Published under license by AIP Publishing.

Abstract: A systematic investigation of voltage effect on the plasma parameters, especially the species densities and chemical reaction mechanisms, in the capacitive N-2 discharges is performed by employing a two-dimensional self-consistent fluid model. The validity of the numerical model is first demonstrated by the qualitative agreement of the calculated and experimental results. Then, the densities, production mechanisms, and loss mechanisms of species from simulation are examined at various voltages. It is found that all the species densities increase monotonically with the voltage, whereas their spatial profiles at lower voltages are quite different from those at higher voltages. The electrons and Nthorn 2 ions are mainly generated by the electron impact ionization of N-2 gas, while the Nthorn ions, whose density is one or two orders of magnitude lower, are mostly formed by the ionization of N atoms. The electron impact dissociation of N-2 gas dominates the generation of N atoms, which are mostly destroyed for the Nthorn ion production. As for the excited N-2 levels, the level conversion processes play a very important role in their production and depletion mechanisms, except for the electron impact excitation of the ground state N-2 molecules. Published under license by AIP Publishing.

Abstract: The unique physical properties of carbyne, a novel carbon nanostructure, have attracted considerable interest in modern nanotechnology. While carbyne synthesis has been accomplished successfully using diverse techniques, the underlying mechanisms governing the carbon monoxide-dependent catalytic synthesis of endohedral carbyne remain poorly understood. In this simulation-based study, we investigate the synthesis of endohedral carbyne from carbon and carbon monoxide radicals in the presence of a nickel catalyst inside double-walled carbon nanotubes with a (5,5)@(10,10) structure. The outcome of our investigation demonstrates that the incorporation of the carbon atom within the Ni-n@(5,5)@(10,10) model system initiates the formation of an elongated carbon chain. In contrast, upon the introduction of carbon monoxide radicals, the growth of the carbyne chain is inhibited as a result of the oxidation of endohedral nickel clusters by oxygen atoms after the initial steps of nucleation. Our findings align with prior theoretical, simulation, and experimental investigations, reinforcing their consistency and providing valuable insights into the synthesis of carbyne-based nanodevices that hold promising potential for future advancements in nanotechnology.

Abstract: The effect of gas heating in laboratory sputter magnetrons is investigated by means of numerical modeling. The model is two-dimensional in the coordinate space and three-dimensional in the velocity space based on the particle-in-cellMonte Carlo collisions technique. It is expanded in a way that allows the inclusion of the neutral plasma particles (fast gas atoms and sputtered atoms), which makes it possible to calculate the gas temperature and its influence on the discharge behavior in a completely self-consistent way. The results of the model are compared to experimental measurements and to other existing simulation results. The results show that gas heating is pressure dependent (rising with the increase in the gas pressure) and should be taken into consideration at pressures above 10 mTorr.

Abstract: The influence of the composition on the crystallographic properties of deposited Mg(M)O (with M=Al, Cr, Ti, Y, and Zr) films is studied. For a flexible control of the composition, dual reactive magnetron sputtering was used as deposition technique. Two different approaches to predict the composition are discussed. The first is an experimental way based on the simple relationship between the deposition rate and the target-substrate distance. The second is a route using a Monte Carlo based particle trajectory code. Both methods require a minimal experimental input and enable the user to quickly predict the composition of complex thin films. Good control and flexibility allow us to study the compositional effects on the growth of Mg(M)O films. Pure MgO thin films were grown with a (111) preferential out-of-plane orientation. When adding M to MgO, two trends were noticed. The first trend is a change in the MgO lattice parameters compared to pure MgO. The second tendency is a decrease in the crystallinity of the MgO phase. The experimentally determined crystallographic properties are shown to be in correspondence with the predicted properties from molecular dynamics simulations.