Abstract: We investigated the mixing and segregation of a system consisting of two different species of particles, having different charges, interacting through a pure Coulomb potential, and confined in a three-dimensional parabolic trap. The structure of the cluster and its normal mode spectrum are analyzed as a function of the relative charge and the relative number of different types of particles. We found that (a) the system can be in a mixed or segregated state depending on the relative charge ratio parameter and (b) the segregation process is mediated by a first or second order structural phase transition which strongly influences the magic cluster properties of the system.

Abstract: We expand on a recent study of a lattice model of interacting particles [Phys. Rev. Lett. 111, 110601 (2013)]. The adsorption isotherm and equilibrium fluctuations in particle number are discussed as a function of the interaction. Their behavior is similar to that of interacting particles in porous materials. Different expressions for the particle jump rates are derived from transition-state theory. Which expression should be used depends on the strength of the interparticle interactions. Analytical expressions for the self-and transport diffusion are derived when correlations, caused by memory effects in the environment, are neglected. The diffusive behavior is studied numerically with kinetic Monte Carlo (kMC) simulations, which reproduces the diffusion including correlations. The effect of correlations is studied by comparing the analytical expressions with the kMC simulations. It is found that the Maxwell-Stefan diffusion can exceed the self-diffusion. To our knowledge, this is the first time this is observed. The diffusive behavior in one-dimensional and higher-dimensional systems is qualitatively the same, with the effect of correlations decreasing for increasing dimension. The length dependence of both the self-and transport diffusion is studied for one-dimensional systems. For long lengths the self-diffusion shows a 1/L dependence. Finally, we discuss when agreement with experiments and simulations can be expected. The assumption that particles in different cavities do not interact is expected to hold quantitatively at low and medium particle concentrations if the particles are not strongly interacting.

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.

Abstract: The diffusion of charged particles interacting through a repulsive Yukawa potential, exp(-r/lambda)/r, confined by a parabolic potential in the y direction and subjected to a periodic substrate potential in the x direction is investigated. Langevin dynamic simulations are used to investigate the effect of the particle density, the amplitude of the periodic substrate, and the range of the interparticle interaction potential on the diffusive behavior of the particles. We found that in general the diffusion is suppressed with increasing the amplitude of the periodic potential, but for specific values of the strength of the substrate potential a remarkable increase of the diffusion is found with increasing the periodic potential amplitude. In addition, we found a strong dependence of the diffusion on the specific arrangement of the particles, e. g., single-chain versus multichain configuration. For certain particle configurations, a reentrant behavior of the diffusion is found as a function of the substrate strength due to structural transitions in the ordering of the particles.

Abstract: The initial stage of nanoparticle formation and growth in radiofrequency acetylene (C2H2) plasmas is investigated by means of a self-consistent one-dimensional fluid model. A detailed chemical kinetic scheme, containing electron impact, ion-neutral, and neutral-neutral reactions, has been developed in order to predict the underlying dust growth mechanisms and the most important dust precursors. The model considers 41 different species (neutrals, radicals, ions, and electrons) describing hydrocarbons (CnHm) containing up to 12 carbon atoms. Possible routes for particle growth are discussed. Both positive and negative ion reaction pathways are considered, as consecutive anion- and cation-molecule reactions seem to lead to a fast build up of the carbon skeleton.

Abstract: We analyze the combined effect of three ingredients of an aggregation model-surface tension, particle flow and particle source-representing typical characteristics of many aggregation growth processes in nature. Through extensive numerical experiments and for different underlying lattice structures we demonstrate that the location of incoming particles and their preferential direction of flow can significantly affect the resulting general shape of the aggregate, while the surface tension controls the surface roughness. Combining all three ingredients increases the aggregate shape plasticity, yielding a wider spectrum of shapes as compared to earlier works that analyzed these ingredients separately. Our results indicate that the considered combination of effects is fundamental for modeling the polymorphic growth of a wide variety of structures in confined geometries and/or in the presence of external fields, such as rocks, crystals, corals, and biominerals.

Abstract: Molecular dynamics simulations are used to study the effects of different cations on the permeation of charged polymers through flat capillaries with heights below 2 nm. Interestingly, we found that, despite being monovalent, Li+ , Na+ , and K+ cations have different effects on polymer permeation, which consequently affects their transmission speed throughout those capillaries. We attribute this phenomenon to the interplay of the cations' hydration free energies and the hydrodynamic drag in front of the polymer when it enters the capillary. Different alkali cations exhibit different surface versus bulk preferences in small clusters of water under the influence of an external electric field. This paper presents a tool to control the speed of charged polymers in confined spaces using cations.

Abstract: We study the mixing of two different kinds of particles, having different charge and/or mass, interacting through a pure Coulomb potential, and confined in a parabolic trap. The structure of the cluster and its normal mode spectrum are analyzed as a function of the ratio of the charges (mass ratio) of the two types of particles. We show that particles are not always arranged in a shell structure. Mixing of the particles goes hand in hand with a large number of metastable states. The normal modes of the system are obtained, and we find that some of the special modes can be tuned by varying the ratio between the charges (masses) of the two species. The degree of mixing of the two type of particles is summarized in a phase diagram, and an order parameter that describes quantitatively the mixing between particles is defined.

Abstract: The thermal properties of a two-dimensional classical cluster of negatively charged particles bound by a punctual positive charge are presented. The melting phenomenon is analyzed and the features which characterize such a solid-liquid transition are highlighted. We found that the presence of metastable states strongly modifies the melting scenario, and that the melting temperature of the system is determined by the height of the saddle point energy separating the ground state and the metastable state. Due to the particular type of confinement potential considered in this paper, we also found that, at sufficiently large temperature, the cluster can become thermally ionized.

Abstract: The ground state of colloidal magnetic particles in a modulated channel are investigated as a function of the tilt angle of an applied magnetic field. The particles are confined by a parabolic potential in the transversal direction while in the axial direction a periodic substrate potential is present. By using Monte Carlo simulations, we construct a phase diagram for the different crystal structures as a function of the magnetic field orientation, strength of the modulated potential, and the commensurability factor of the system. Interestingly, we found first-and second-order phase transitions between different crystal structures, which can be manipulated by the orientation of the external magnetic field. A reentrant behavior is found between two-and four-chain configurations, with continuous second-order transitions. Novel configurations are found consisting of frozen solitons of defects. By changing the orientation and/or strength of the magnetic field and/or the strength and periodicity of the substrate potential, the system transits through different phases.

Abstract: The electron distribution function (EDF) in an electron cyclotron resonance (ECR) discharge is far from Maxwellian. The self-consistent simulation of ECR discharges requires a calculation of the EDF on every magnetic line for various ion density profiles. The straightforward self-consistent simulation of ECR discharges using the Monte Carlo technique for the EDF calculation is very computer time expensive, since the electron and ion time scales are very different. An electron Boltzmann kinetic equation averaged over the fast electron bouncing and pitch-angle scattering was derived in order to develop an effective and operative tool for the fast modeling (FM) of low-pressure ECR discharges. An analytical solution for the EDF in a loss cone was derived. To check the validity of the FM, one-dimensional (in coordinate) and two-dimensional (in velocity) Monte Carlo simulation codes were developed. The validity of the fast modeling method is proved by comparison with the Monte Carlo simulations. The complete system of equations for FM is presented and ready for use in a comprehensive study of ECR discharges. The variations of plasma density and of wall and sheath potentials are analyzed by solving a self-consistent set of equations for the EDF.