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 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: A classical two-dimensional (2D) model for an artificial atom is used to make a numerical exact study of elastic and nonelastic scattering. Interesting differences in the scattering angle distribution between this model and the well-known Rutherford scattering are found in the small energy and/or small impact parameter scattering regime. For scattering off a classical 2D hydrogen atom different phenomena such as ionization, exchange of particles, and inelastic scattering can occur. A scattering regime diagram is constructed as function of the impact parameter (b) and the initial velocity (v) of the incoming particle. In a small regime of the (b,v) space the system exhibits chaos, which is studied in more detail. Analytic expressions for the scattering angle are given in the high impact parameter asymptotic limit.

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: A two-dimensional (2D) binary colloidal system consisting of interacting dipoles is investigated using an analytical approach. Within the harmonic approximation we obtain the phonon spectrum of the system as a function of the composition, dipole-moment ratio, and mass ratio between the small and big particles. Through a systematic analysis of the phonon spectra we are able to determine the stability region of the different lattice structures of the colloidal alloys. The gaps in the phonon frequency spectrum, the optical frequencies in the long-wavelength limit, and the sound velocity are discussed as well. Using the modified Lindemann criterion and within the harmonic approximation we estimate the melting temperature of the sublattice generated by the big 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 ordered configurations of a monolayer of interacting magnetic dipoles confined in a circular parabolic potential are investigated as a function of the dipole moment of the particles. Despite the circular confinement, we find very asymmetric ordered structures like chains and Y-shaped configurations when a magnetic field is applied parallel to the plane of the particles. The normal-mode spectrum of the particles and its dependence on the magnetic field and the strength of the dipole moment of the particles are studied. The vibrational and rotational modes of the spectrum, which are associated with the stability of the system, are investigated in detail. The number of particles is varied and we found different ordering of the particles for different values of the dipole moment and the magnetic field. A ring structure with a large number of particles is observed for high values of the dipole moment of the particles.

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: Molecular dynamics simulations are used to investigate the structural properties of a two-dimensional ensemble of magnetic rods, which are modeled as aligned single dipolar beads. The obtained self-assembled configurations can be characterized as (1) clusters, (2) percolated, and (3) ordered structures, and their structural properties are investigated in detail. By increasing the aspect ratio of the magnetic rods, we show that the percolation transition is suppressed due to the reduced mobility of the rods in two dimensions. Such a behavior is opposite to the one observed in three dimensions. A magnetic bulk phase is found with local ferromagnetic order and an unusual nonmonotonic behavior of the nematic order is observed.

Abstract: Similar to the spontaneous formation of colonies of bacteria, flocks of birds, or schools of fish, “living crystals” can be formed by artificial self-propelled particles such as Janus colloids. Unlike usual solids, these “crystals” are far from thermodynamic equilibrium. They fluctuate in time forming a crystalline structure, breaking apart and re-forming again. We propose a method to stabilize living crystals by applying a weak confinement potential that does not suppress the ability of the particles to perform self-propelled motion, but it stabilizes the structure and shape of the dynamical clusters. This gives rise to such configurations of living crystals as “living shells” formed by Janus colloids. Moreover, the shape of the stable living clusters can be controlled by tuning the potential strength. Our proposal can be verified experimentally with either artificial microswimmers such as Janus colloids, or with living active matter.

Abstract: We report on a molecular dynamics study of small classical two-dimensional clusters with ringlike configurations. We focus on the particles motion at low temperatures before the radial and angular melting sets in. It is shown that in magic number configurations a local radial melting of subshells occur, which is related to the intershell rotation.