Abstract: We report on the dynamical properties of interacting colloids confined to one dimension and subjected to external periodic energy landscapes. We particularly focus on the influence of hydrodynamic interactions on the mean-square displacement. Using Brownian dynamics simulations, we study colloidal systems with two types of repulsive interparticle interactions, namely, Yukawa and superparamagnetic potentials. We find that in the homogeneous case, hydrodynamic interactions lead to an enhancement of the particle mobility and the mean-square displacement at long times scales as t(alpha), with alpha = 1/2 + epsilon and epsilon being a small correction. This correction, however, becomes much more important in the presence of an external field, which breaks the homogeneity of the particle distribution along the line and, therefore, promotes a richer dynamical scenario due to the hydrodynamical coupling among particles. We provide here the complete dynamical scenario in terms of the external potential parameters: amplitude and commensurability.

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 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: 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: Diffusive properties of interacting magnetic dipoles confined in a parabolic narrow channel and in the presence of a periodic modulated (corrugated) potential along the unconfined direction are studied using Brownian dynamics simulations. We compare our simulation results with the analytical result for the effective diffusion coefficient of a single particle by Festa and d'Agliano [Physica A 90, 229 (1978)] and show the importance of interparticle interaction on the diffusion process. We present results for the diffusion of magnetic dipoles as a function of linear density, strength of the periodic modulation and commensurability factor.

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: Janus particles present an important class of building blocks for directional assembly. These are compartmentalized colloids with two different hemispheres. In this work we consider a three-dimensional model of Janus spheres that contain one hydrophobic and one charged hemisphere. Using molecular dynamics simulations, we study the morphology of these particles when confined in a channel-like environment. The interplay between the attractive and repulsive forces on each particle gives rise to a rich phase space where the relative orientation of each particle plays a dominant role in the formation of large-scale clusters. The interest in this system is primarily due to the fact that it could give a better understanding of the mechanisms of the formation of polar membranes. A variety of ordered membranelike morphologies is found consisting of single and multiple connected chain configurations. The helicity of these chains can be chosen by simply changing the salt concentration of the solution. Special attention is given to the formation of Bernal spirals. These helices are composed of regular tetrahedra and are known to exhibit nontrivial translational and rotational symmetry.

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: 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: Due to the potential application of DNA for biophysics and optoelectronics, the electronic energy states and transitions of this genetic material have attracted a great deal of attention recently. However, the fluorescence and corresponding physical process of DNA under optical excitation with photon energies below ultraviolet are still not fully clear. In this work, we experimentally investigate the photoluminescence (PL) properties of single-stranded DNA (ssDNA) samples under near-ultraviolet (NUV) and visible excitations (270 similar to 440 nm). Based on the dependence of the PL peak wavelength (lem) upon the excitation wavelength (lex), the PL behaviors of ssDNA can be approximately classified into two categories. In the relatively short excitation wavelength regime, lem is nearly constant due to exciton-like transitions associated with delocalized excitonic states and excimer states. In the relatively long excitation wavelength range, a linear relation of lem = Alex + B with A 0 or A < 0 can be observed, which comes from electronic transitions related to coupled vibrational-electronic levels. Moreover, the transition channels in different excitation wavelength regimes and the effects of strand length and base type can be analyzed on the basis of these results. These important findings not only can give a general description of the electronic energy states and transitional behaviors of ssDNA samples under NUV and visible excitations, but also can be the basis for the application of DNA in nanoelectronics and optoelectronics.

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.

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.