Abstract: We apply density-functional theory to study the adsorption of water clusters on the surface of a graphene sheet and find i) graphene is highly hydrophobic and ii) adsorbed water has very little effect on the electronic structure of graphene. A single water cluster on graphene has a very small average dipole moment which is in contrast with an ice layer that exhibits a strong dipole moment.

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: Using Brownian dynamics simulations, we investigate the dynamics of colloids confined in two-dimensional narrow channels driven by a nonuniform force Fdr(y). We considered linear-gradient, parabolic, and deltalike driving-force profiles. This driving force induces melting of the colloidal solid (i.e., shear-induced melting), and the colloidal motion experiences a transition from elastic to plastic regime with increasing Fdr. For intermediate Fdr (i.e., in the transition region) the response of the system, i.e., the distribution of the velocities of the colloidal chains õi(y), in general does not coincide with the profile of the driving force Fdr(y), and depends on the magnitude of Fdr, the width of the channel, and the density of colloids. For example, we show that the onset of plasticity is first observed near the boundaries while the motion in the central region is elastic. This is explained by: (i) (in)commensurability between the chains due to the larger density of colloids near the boundaries, and (ii) the gradient in Fdr. Our study provides a deeper understanding of the dynamics of colloids in channels and could be accessed in experiments on colloids (or in dusty plasma) with, e.g., asymmetric channels or in the presence of a gradient potential field.

Abstract: Single-file diffusion (SFD) of an infinite one-dimensional chain of interacting particles has a long-time mean-square displacement ∝t1/2, independent of the type of interparticle repulsive interaction. This behavior is also observed in finite-size chains, although only for certain intervals of time t depending on the chain length L, followed by the ∝t for t→∞, as we demonstrate for a closed circular chain of diffusing interacting particles. Here, we show that spatial correlation of noise slows down SFD and can result, depending on the amount of correlated noise, in either subdiffusive behavior ∝tα, where 0<α<1/2, or even in a total suppression of diffusion (in the limit N→∞). Spatial correlation can explain the subdiffusive behavior in recent SFD experiments in circular channels.

Abstract: The influence of an ellipsoid on the angular order of two-dimensional classical clusters is investigated through Brownian dynamics simulations. We found the following: (1) The presence of an ellipsoid does not influence the start of the angular melting, but reduces the rate at which the inner rings can rotate with respect to each other. (2) Even a small eccentricity of the ellipsoid leads to a stabilization of the angular order of the system. (3) Depending on the position of the ellipsoid in the cluster, a reentrant behavior in the angular order is observed before full radial melting of the cluster sets in. (4) The ellipsoid can lead to a two-step angular melting process: First, the rotation of the inner rings with respect to each other is hindered by the ellipsoid, but on further increasing the kinetic energy of the system, the ellipsoid just starts to behave as a spherical particle with different mobility. The effect of an ellipsoid on the molten system does not depend crucially on the interparticle interaction, but a softer parabolic confinement reduces the angular stabilization.

Abstract: Molecular dynamic simulations are performed to investigate the melting process of small three-dimensional clusters (i.e., systems with one and two shells) of classical charged particles trapped in an isotropic parabolic potential. The confined particles interact through a repulsive potential. We find that the ground-state configurations for systems with N=6, 12, 13, and 38 particles interacting through a Coulomb potential are magic clusters. Such magic clusters have an octahedral or icosahedral symmetry and are found to have a large stability against intrashell diffusion leading to an intershell melting transition prior to the intrashell and radial melting process. For systems with two shells a local radial melting of subshells is found at low temperatures resulting in a structural transition leading to an increased symmetry of the ordered system. Using Lindemanns criterion the different melting temperatures are determined and the influence of the screening of the interparticle interaction was investigated. A normal mode analysis is performed and some of the normal modes are found to be determinantal for the melting process.

Abstract: We studied equally charged particles, suspended in a complex plasma, which move in a plane and interact with a screened Coulomb potential (Yukawa type) and with an additional external confining parabolic potential in one direction, which makes the system quasi-one-dimensional (Q1D). The normal modes of the system are studied in the presence of dissipation. We also investigated how a perpendicular magnetic field couples the phonon modes with each other. Two different ways of exciting the normal modes are discussed: (1) a uniform excitation of the Q1D lattice, and (2) a local forced excitation of the system in which one particle is driven by, e.g., a laser. Our results are in very good agreement with recent experimental findings on a finite single chain system [Liu , Phys. Rev. Lett. 91, 255003 (2003)]. Predictions are made for the normal modes of multichain structures in the presence of damping.

Abstract: Strong terahertz (1THz=1012Hz) radiation can be generated by the electron oscillation in fs-laser-induced wake fields. The interaction of a fs-laser pulse with a low-density plasma layer is studied in detail using numerical simulations. The spatial distribution and temporal evolution of terahertz electron current developed in a low-density plasma layer are presented, which enables us to calculate the intensity distribution of THz radiation. It is shown that laser and plasma parameters, such as laser intensity, pulse width, and background plasma density, are of key importance to the process. The optimum condition for wake-field excitation and terahertz emission is discussed upon the simulation results. Radiation peaked at 6.4 THz, with 900 fs duration and 9% bandwidth, can be generated in a plasma of density 5×1017cm−3. It turns out that the maximum radiation intensity scales as n03a04 when wake field is resonantly excited, where n0 and a0 are, respectively, the plasma density and the normalized field amplitude of the laser pulse.

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: 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.