Phase diagram and critical properties of Yukawa bilayers

Core-softened colloid under extreme geometrical confinement

Geometrical constraints offer a promising strategy for assembling colloidal crystal structures that are not typically observed in bulk or under 2D conditions. Core-softened colloids, in particular, have emerged as versatile chemical building blocks with applications across various scientific and technological areas. In this study, we investigate the behavior of a core-softened model confined between two parallel walls. Employing molecular dynamics simulations, we analyze the system's response under extreme confinement, where only one or two layers of colloids are permitted. The system comprises particles modeled by a ramp-like potential confined within slit nanoslits created by two flat, purely repulsive walls with a lateral side L separated by a distance Lz. Through a systematic analysis of the phase behavior as Lz increases, or as the system undergoes decompression, for different values of L, we identified a mono-to-bilayer transition associated with changes in the colloidal structure. In the monolayer regime, we observed solid phases at lower densities than those observed in the 2D case. Importantly, we demonstrated that confinement at specific Lz values, allowing particle arrangement into two layers, can lead to the emergence of the square phase, which was not observed under monolayer or 2D conditions. By correlating thermodynamic, translational, and orientational ordering, as well as the dynamics of this confined colloidal system, our findings offer valuable insights into the utilization of geometrical constraints to induce and manipulate structural changes.

Strongly coupled Yukawa trilayer liquid: Structure and dynamics

The equilibrium structure and the dispersion relations of collective excitations in trilayer Yukawa systems in the strongly coupled liquid regime are examined. The equilibrium correlations reveal a variety of structures in the liquid phase, reminiscent of the corresponding structures in the solid phase. At small layer separation substitutional disorder becomes the governing feature. Theoretical dispersion relations are obtained by applying the quasilocalized charge approximation (QLCA) formalism, while numerical data are generated by microcanonical molecular dynamics simulations. The dispersions and polarizations of the collective excitations obtained through both of these methods are compared and discussed in detail. We find that the QLCA method is, in general, very satisfactory, but that there are phenomena not covered by the QLCA. In particular, by analyzing the dynamical longitudinal and transverse current fluctuation spectra we discover the existence of a structure not related to the collective mode spectra. This also provides insight into the long-standing problem of the gap frequency discrepancy, observed in strongly coupled layered systems in earlier studies.

Buckling of two-dimensional plasma crystals with nonreciprocal interactions

Laboratory realizations of two-dimensional (2D) plasma crystals typically involve monodisperse microparticles confined into horizontal monolayers in radio-frequency (rf) plasma sheaths. This gives rise to the so-called plasma wakes beneath the microparticles. The presence of wakes renders the interactions in such systems nonreciprocal, a fact that can lead to a quite different behavior from the one expected for their reciprocal counterparts. Here we examine the buckling of a hexagonal 2D plasma crystal, occurring as the confinement strength is decreased, taking explicitly into account the nonreciprocity of the system via a well-established point-wake model. We observe that for a finite wake charge, the monolayer hexagonal crystal undergoes a transition first to a bilayer hexagonal structure, unrealizable in harmonically confined reciprocal Yukawa systems, and subsequently to a bilayer square structure. Our theoretical results are confirmed by molecular dynamics simulations for experimentally relevant parameters, indicating the potential of their observation in state-of-the-art experiments with 2D complex plasmas.

Strong-coupling theory of counterions with hard cores between symmetrically charged walls

By a combination of Monte Carlo simulations and analytical calculations, we investigate the effective interactions between highly charged planar interfaces, neutralized by mobile counterions (salt-free system). While most previous analysis have focused on pointlike counterions, we treat them as charged hard spheres. We thus work out the fate of like-charge attraction when steric effects are at work. The analytical approach partitions counterions in two subpopulations, one for each plate, and integrates out one subpopulation to derive an effective Hamiltonian for the remaining one. The effective Hamiltonian features plaquette four-particle interactions, and it is worked out by computing a Gibbs-Bogoliubov inequality for the free energy. At the root of the treatment is the fact that under strong electrostatic coupling, the system of charges forms an ordered arrangement, that can be affected by steric interactions. Fluctuations around the reference positions are accounted for. To dominant order at high coupling, it is found that steric effects do not significantly affect the interplate effective pressure, apart at small distances where hard-sphere overlap are unavoidable, and thus rule out configurations.