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: Nonequilibrium self-organized patterns formed by particles interacting through competing range interaction are driven over a substrate by an external force. We show that, with increasing driving force, the preexisted static patterns evolve into dynamic patterns either via disordered phase or depinned patterns or via the formation of nonequilibrium stripes. Strikingly, the stripes are formed either in the direction of the driving force or in the transverse direction, depending on the pinning strength. The revealed dynamical patterns are summarized in a dynamical phase diagram.

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 autoregressive neural networks are emerging as a powerful computational tool to solve relevant problems in classical and quantum mechanics. One of their appealing functionalities is that, after they have learned a probability distribution from a dataset, they allow exact and efficient sampling of typical system configurations. Here we employ a neural autoregressive distribution estimator (NADE) to boost Markov chain Monte Carlo (MCMC) simulations of a paradigmatic classical model of spin-glass theory, namely, the two-dimensional Edwards-Anderson Hamiltonian. We show that a NADE can be trained to accurately mimic the Boltzmann distribution using unsupervised learning from system configurations generated using standard MCMC algorithms. The trained NADE is then employed as smart proposal distribution for the Metropolis-Hastings algorithm. This allows us to perform efficient MCMC simulations, which provide unbiased results even if the expectation value corresponding to the probability distribution learned by the NADE is not exact. Notably, we implement a sequential tempering procedure, whereby a NADE trained at a higher temperature is iteratively employed as proposal distribution in a MCMC simulation run at a slightly lower temperature. This allows one to efficiently simulate the spin-glass model even in the low-temperature regime, avoiding the divergent correlation times that plague MCMC simulations driven by local-update algorithms. Furthermore, we show that the NADE-driven simulations quickly sample ground-state configurations, paving the way to their future utilization to tackle binary optimization problems.

Abstract: The dielectric properties of confined water is of fundamental interest and is still controversial. For water confined in channels with height smaller than h = 8 angstrom, we found a commensurability effect and an extraordinary decrease in the out-of-plane dielectric constant down to the limit of the dielectric constant of optical water. Spatial resolved polarization density data obtained from molecular dynamics simulations are found to be antisymmetric across the channel and are used as input in a mean-field model for the dielectric constant as a function of the height of the channel for h > 15 angstrom. Our results are in excellent agreement with a recent experiment [L. Fumagalli et al., Science 360, 1339 (2018)].

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 global relaxation process after pulsed laser induced photodetachment in a rf electronegative SIH4 discharge is studied by a self-consistent kinetic one-dimensional particle-in-cell-Monte Carlo model. Our results reveal a comprehensive physical picture of the relaxation process, including the main plasma variables, after a perturbation up to the full recovery of the steady state. A strong influence of the photodetachment on the discharge is found, which results from an increase of the electron density, leading to a weaker bulk field, and hence to a drop in the high energy tail of the electron energy distribution function (EEDF), a reduction of the reaction rates of electron impact attachment and ionization, and a subsequent decrease of the positive and negative ion densities. All the plasma quantities related to electrons recover synchronously. The recovery time of the ion densities is about 1-2 orders of magnitude longer than that of the electrons due to different recovery mechanisms. The modeled behavior of all the charged particles agrees very well with experimental results from the literature. In addition, our work clarifies some unclear processes assumed in the literature, such as the relaxation of the EEDF, the evolution of the electric field, and the recovery of negative ions.

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

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.