Abstract: Molecular-dynamic simulations were performed in order to investigate the melting processes of isotropically confined binary systems. We considered two species of particles, which differ by their amount of electric charge. A Lindemann type of criterion was used to determine the angular melting temperature. We demonstrate that the magic-to-normal cluster transition can evolve in two distinct ways, that is, through a structural phase transition of the first order or via a smooth transition where an increase of the shells' width leads to a continuous decreasing mechanical stability of the system. Moreover, for large systems, we demonstrate that the internal cluster exerts a minor effect on the mechanical stability of the external shell. Furthermore, we show that highly symmetric configurations, such as those found for multiple ring structures, have large mechanical stability, i.e., high angular melting temperature.

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: 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: Janus particles present an important class of building blocks for directional assembly. These are compartmentalized colloids with two different hemispheres. Here, we consider a two-dimensional model of Janus disks consisting of a hydrophobic semicircle and an electro-negatively charged one. Placed in a solution, the hydrophobic sides will attract each other while the charged sides will give rise to a repulsive force. Using molecular dynamics simulations, we study the morphology of these particles when confined in a channel-like environment using a one dimensional harmonic confinement potential. The interest to this system is first of all due to the fact that it could serve as a simple model for membrane formation. Indeed, the recently synthesized new class of artificial amphiphiles, known as Janus dendrimers, were shown to self-assemble in bilayer structures mimicking biological membranes. In turn, Janus particles that combine the amphiphilicity and colloidal rigidity serve as a good model for Janus dendrimers. A variety of ordered membrane-like morphologies are found consisting of single and multiple chain configurations with different orientations of the particles with respect to each other that we summarize in a phase diagram.

Abstract: The diffusion of a system of ferromagnetic dipoles confined in a quasi-one-dimensional parabolic trap is studied using Brownian dynamics simulations. We show that the dynamics of the system is tunable by an in-plane external homogeneous magnetic field. For a strong applied magnetic field, we find that the mobility of the system, the exponent of diffusion, and the crossover time among different diffusion regimes can be tuned by the orientation of the magnetic field. For weak magnetic fields, the exponent of diffusion in the subdiffusive regime is independent of the orientation of the external field. DOI: 10.1103/PhysRevE.87.012307

Abstract: We study the structure and phonon spectrum of a two-dimensional bilayer system of classical charged dipoles oriented perpendicular to the plane of the layers for equal density in each layer. This system can be tuned through six different crystalline phases by changing the interlayer separation or the charge and/or dipole moment of the particle. The presence of the charge on the dipole particles is responsible for the nucleation of five staggered phases and a disordered phase which are not found in the magnetic dipole bilayer system. These extra phases are a consequence of the competition between the repulsive Coulomb and the attractive dipole interlayer interaction. We present the phase diagram and determine the order of the phase transitions. The phonon spectrum of the system was calculated within the harmonic approximation, and a nonmonotonic behavior of the phonon spectrum is found as a function of the effective strength of the interparticle interaction. The stability of the different phases is determined.

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 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 ordering of N equally charged particles (-e) moving in two dimensions and confined by a Coulomb potential, resulting from a displaced positive charge Ze is discussed. This is a classical model system for atoms. We obtain the configurations of charged particles which, depending on the value of N and Z, may result in ring structures, hexagonal-type configurations, and for N/Z approximate to 1 in an inner structure of particles which is separated by an outer ring of particles. For N/Z << 1, the Hamiltonian of the parabolic confinement case is recovered. For N/Z approximate to 1, the configurations are very different from those found in the case of a parabolic confinement potential. A hydrodynamic analysis is presented in order to highlight the correlations effects.

Abstract: The purpose of this Brief Report is to point out the mistake in a formula for anisotropic electron scattering, previously published in Phys. Rev. A 41, 1112 (1990), which is widely used in Monte Carlo models of gas discharges. Anisotropic electron scattering is investigated based on the screened Coulomb potential between electrons and neutral atoms. The approach is also applied for electron scattering by nonpolar neutral molecules. Differential cross sections for electron scattering by Ar, N2, and CH4 are constructed on the basis of momentum and integrated cross sections. The formula derived in this paper is useful for Monte Carlo simulations of gas discharges.

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.