Ordered ground state configurations of the asymmetric Wigner bilayer system-Revisited with unsupervised learning.

We have reanalyzed the rich plethora of ground state configurations of the asymmetric Wigner bilayer system that we had recently published in a related diagram of states [Antlanger et al., Phys. Rev. Lett. 117, 118002 (2016)], comprising roughly 60 000 state points in the phase space spanned by the distance between the plates and the charge asymmetry parameter of the system. In contrast to this preceding contribution where the classification of the emerging structures was carried out "by hand," we have used for the present contribution machine learning concepts, notably based on a principal component analysis and a k-means clustering approach: using a 30-dimensional feature vector for each emerging structure (containing relevant information, such as the composition of the configuration as well as the most relevant order parameters), we were able to reanalyze these ground state configurations in a considerably more systematic and comprehensive manner than we could possibly do in the previously published classification scheme. Indeed, we were now able to identify new structures in previously unclassified regions of the parameter space and could considerably refine the previous classification scheme, thereby identifying a rich wealth of new emerging ground state configurations. Thorough consistency checks confirm the validity of the newly defined diagram of states.

The asymmetric Wigner bilayer.

We present a comprehensive discussion of the so-called asymmetric Wigner bilayer system, where mobile point charges, all of the same sign, are immersed into the space left between two parallel, homogeneously charged plates (with possibly different charge densities). At vanishing temperatures, the particles are expelled from the slab interior; they necessarily stick to one of the two plates and form there ordered sublattices. Using complementary tools (analytic and numerical), we study systematically the self-assembly of the point charges into ordered ground state configurations as the inter-layer separation and the asymmetry in the charge densities are varied. The overwhelming plethora of emerging Wigner bilayer ground states can be understood in terms of the competition of two strategies of the system: net charge neutrality on each of the plates on the one hand and particles' self-organization into commensurate sublattices on the other hand. The emerging structures range from simple, highly commensurate (and thus very stable) lattices (such as staggered structures, built up by simple motives) to structures with a complicated internal structure. The combined application of our two approaches (whose results agree within remarkable accuracy) allows us to study on a quantitative level phenomena such as over- and underpopulation of the plates by the mobile particles, the nature of phase transitions between the emerging phases (which pertain to two different universality classes), and the physical laws that govern the long-range behaviour of the forces acting between the plates. Extensive, complementary Monte Carlo simulations in the canonical ensemble, which have been carried out at small, but finite temperatures along selected, well-defined pathways in parameter space confirm the analytical and numerical predictions within high accuracy. The simple setup of the Wigner bilayer system offers an attractive possibility to study and to control complex scenarios and strategies of colloidal self-assembly, via the variation of two system parameters.

Rich Polymorphic Behavior of Wigner Bilayers.

Self-assembly into target structures is an efficient material design strategy. Combining analytical calculations and computational techniques of evolutionary and Monte Carlo types, we report about a remarkable structural variability of Wigner bilayer ground states, when charges are confined between parallel charged plates. Changing the interlayer separation, or the plate charge asymmetry, a cascade of ordered patterns emerges. At variance with the symmetric case phenomenology, the competition between commensurability features and charge neutralization leads to long range attraction, appearance of macroscopic charges, exotic phases, and nonconventional phase transitions with distinct critical indices, offering the possibility of a subtle, but precise and convenient control over patterns.

Equilibrium structures of anisometric, quadrupolar particles confined to a monolayer.

We investigate the structural properties of a two-dimensional system of ellipsoidal particles carrying a linear quadrupole moment in their center. These particles represent a simple model for a variety of uncharged, non-polar conjugated organic molecules. Using optimization tools based on ideas of evolutionary algorithms, we first examine the ground state structures as we vary the aspect ratio of the particles and the pressure. Interestingly, we find, besides the intuitively expected T-like configurations, a variety of complex structures, characterized with up to three different particle orientations. In an effort to explore the impact of thermal fluctuations, we perform constant-pressure molecular dynamics simulations within a range of rather low temperatures. We observe that ground state structures formed by particles with a large aspect ratio are in particular suited to withstand fluctuations up to rather high temperatures. Our comprehensive investigations allow for a deeper understanding of molecular or colloidal monolayer arrangements under the influence of a typical electrostatic interaction on a coarse-grained level.

Crystal phases of soft spheres systems in a slab geometry.

We have identified the ground state configurations of soft particles (interacting via inverse power potentials) confined between two hard, impenetrable walls. To this end we have used a highly reliable optimization scheme at vanishing temperature while varying the wall separation over a representative range. Apart from the expected layered triangular and square structures (which are compatible with the three-dimensional bulk fcc lattice), we have identified a cascade of highly complex intermediate structures. Taking benefit of the general scaling properties of inverse power potentials, we could identify - for a given softness value - one single master curve which relates the energy to the wall separation, irrespective of the density of the system. Via extensive Monte Carlo simulations, we have performed closer investigations of these intermediate structures at finite temperature: we could provide evidence to which extent these particle arrangements remain stable over a relatively large temperature range.

Communications: The Metropolis Monte Carlo finite element algorithm for electrostatic interactions.

The Metropolis Monte Carlo algorithm with the finite element method applied to compute electrostatic interaction energy between charge densities is described in this work. By using the finite element method to integrate numerically Poisson's equation, it is shown that the computing time to obtain the acceptance probability of an elementary trial move does not, in principle, depend on the number of charged particles present in the system.

Structure and thermodynamics of a ferrofluid bilayer.

We present extensive Monte Carlo simulations for the thermodynamic and structural properties of a planar bilayer of dipolar hard spheres for a wide range of densities, dipole moments, and layer separations. Expressions for the stress and pressure tensors of the bilayer system are derived. For all thermodynamic states considered, the interlayer energy is shown to be attractive and much smaller than the intralayer contribution to the energy. It vanishes at layer separations of the order of two hard sphere diameters. The normal pressure is negative and decays as a function of layer separation h as -1/h;{5} . Intralayer and interlayer pair distribution functions and angular correlation functions are presented. Despite the weak interlayer energy, strong positional and orientational correlations exist between particles in the two layers.