Columnar-disorder phase boundary in a mixture of hard squares and dimers

Phases of the hard-plate lattice gas on a three-dimensional cubic lattice

We study the phase diagram of a lattice gas of 2×2×1 hard plates on the three-dimensional cubic lattice. Each plate covers an elementary plaquette of the cubic lattice, with the constraint that a site can belong to utmost one plate. We focus on the isotropic system, with equal fugacities for the three orientations of plates. We show, using grand canonical Monte Carlo simulations, that the system undergoes two phase transitions when the density of plates is increased: the first from a disordered fluid phase to a layered phase, and the second from the layered phase to a sublattice-ordered phase. In the layered phase, the system breaks up into disjoint slabs of thickness two along one spontaneously chosen Cartesian direction, corresponding to a twofold (Z_{2}) symmetry breaking of translation symmetry along the layering direction. Plates with normals perpendicular to this layering direction are preferentially contained entirely within these slabs, while plates straddling two adjacent slabs have a lower density, thus breaking the symmetry between the three types of plates. We show that the slabs exhibit two-dimensional power-law columnar order even in the presence of a nonzero density of vacancies. In contrast, interslab correlations of the two-dimensional columnar order parameter decay exponentially with the separation between the slabs. In the sublattice-ordered phase, there is twofold symmetry breaking of lattice translation symmetry along all three Cartesian directions. We present numerical evidence that the disordered to layered transition is continuous and consistent with universality class of the three-dimensional O(3) model with cubic anisotropy, while the layered to sublattice transition is first-order in nature.

Entropy-driven phases at high coverage adsorption of straight rigid rods on two-dimensional square lattices

Polymers are frequently deposited on different surfaces, which has attracted the attention of scientists from different viewpoints. In the present approach polymers are represented by rigid rods of length k (k-mers), and the substrate takes the form of an L×L square lattice whose lattice constant matches exactly the interspacing between consecutive elements of the k-mer chain. We briefly review the classical description of the nematic transition presented by this system for k≥7 observing that the high-coverage (θ) transition deserves a more careful analysis from the entropy point of view. We present a possible viewpoint for this analysis that justifies the phase transitions. Moreover, we perform Monte Carlo (MC) simulations in the grand canonical ensemble, supplemented by thermodynamic integration, to first calculate the configurational entropy of the adsorbed phase as a function of the coverage, and then to explore the different phases (and orientational transitions) that appear on the surface with increasing the density of adsorbed k-mers. In the limit of θ→1 (full coverage) the configurational entropy is obtained for values of k ranging between 2 and 10. MC data are discussed in comparison with recent analytical results [D. Dhar and R. Rajesh, Phys. Rev. E 103, 042130 (2021)2470-004510.1103/PhysRevE.103.042130]. The comparative study allows us to establish the applicability range of the theoretical predictions. Finally, the structure of the high-coverage phase is characterized in terms of the statistics of k×l domains (domains of l parallel k-mers adsorbed on the surface). A distribution of finite values of l (l≪L) is found with a predominance of k×1 (single k-mers) and k×k domains. The distribution is the same in each lattice direction, confirming that at high density the adsorbed phase goes to a state with mixed orientations and no orientational preference. An order parameter measuring the number of k×k domains in the adsorbed layer is introduced.

Entropy of fully packed hard rigid rods on d-dimensional hypercubic lattices

We determine the asymptotic behavior of the entropy of full coverings of a L×M square lattice by rods of size k×1 and 1×k, in the limit of large k. We show that full coverage is possible only if at least one of L and M is a multiple of k, and that all allowed configurations can be reached from a standard configuration of all rods being parallel, using only basic flip moves that replace a k×k square of parallel horizontal rods by vertical rods, and vice versa. In the limit of large k, we show that the entropy per site S_{2}(k) tends to Ak^{-2}lnk, with A=1. We conjecture, based on a perturbative series expansion, that this large-k behavior of entropy per site is superuniversal and continues to hold on all d-dimensional hypercubic lattices, with d≥2.

Adsorption kinetics and thermodynamic properties of a binary mixture of hard-core particles on a square lattice

The adsorption kinetics and thermodynamic properties of a binary mixture on a square lattice are studied using the random sequential adsorption with surface diffusion (RSAD). We compare the adsorption of binary species with different equilibrium rate constants and effective rates of adsorption to a surface and find that the temporal evolution of surface coverages of both species can be obtained through the use of the blocking function of a system with irreversible adsorption of highly diffusive particles. Binary mixtures, when one of the components follows the random sequential adsorption (RSA) without surface diffusion and the other follows the RSAD model, display competitive adsorption in addition to cooperative phenomena. Specifically, (i) species replacement occurs over a long period of time, while the total coverage remains unchanged after a short time, (ii) the presence of the RSAD component shifts the jamming coverage to the higher values, and (iii) the maximum jamming coverage is obtained when the effective adsorption of the RSA type components is lower than the other adsorbing particles.

Alternative characterization of the nematic transition in deposition of rods on two-dimensional lattices

We revisit the problem of excluded volume deposition of rigid rods of length k unit cells over square lattices. Two new features are introduced: (a) two new short-distance complementary order parameters, called Π and Σ, are defined, calculated, and discussed to deal with the phases present as coverage increases; (b) the interpretation is now done beginning at the high-coverage ordered phase which allows us to interpret the low-coverage nematic phase as an ergodicity breakdown present only when k≥7. In addition the data analysis invokes both mutability (dynamical information theory method) and Shannon entropy (static distribution analysis) to further characterize the phases of the system. Moreover, mutability and Shannon entropy are compared, and we report the advantages and disadvantages they present for their use in this problem.

Positional ordering of hard adsorbate particles in tubular nanopores

The phase behavior and structural properties of a monolayer of hard particles is examined in such a confinement where the adsorbed particles are constrained to the surface of a narrow hard cylindrical pore. The diameter of the pore is chosen such that only first- and second-neighbor interactions occur between the hard particles. The transfer operator method of [Percus and Zhang, Mol. Phys. 69, 347 (1990)MOPHAM0026-897610.1080/00268979000100241] is reformulated to obtain information about the structure of the monolayer. We have found that a true phase transition is not possible in the examined range of pore diameters. The monolayer of hard spheres undergoes a structural change from fluidlike order to a zigzaglike solid one with increasing surface density. The case of hard cylinders is different in the sense that a layering takes place continuously between a low-density one-row and a high-density two-row monolayer. Our results reveal a clear discrepancy with classical density functional theories, which do not distinguish smecticlike ordering in bulk from that in narrow periodic pores.

Phase transitions in a system of hard Y-shaped particles on the triangular lattice

We study the different phases and the phase transitions in a system of Y-shaped particles, examples of which include immunoglobulin-G and trinaphthylene molecules, on a triangular lattice interacting exclusively through excluded volume interactions. Each particle consists of a central site and three of its six nearest neighbors chosen alternately, such that there are two types of particles which are mirror images of each other. We study the equilibrium properties of the system using grand canonical Monte Carlo simulations that implement an algorithm with cluster moves that is able to equilibrate the system at densities close to full packing. We show that, with increasing density, the system undergoes two entropy-driven phase transitions with two broken-symmetry phases. At low densities, the system is in a disordered phase. As intermediate phases, there is a solidlike sublattice phase in which one type of particle is preferred over the other and the particles preferentially occupy one of four sublattices, thus breaking both particle symmetry as well as translational invariance. At even higher densities, the phase is a columnar phase, where the particle symmetry is restored, and the particles preferentially occupy even or odd rows along one of the three directions. This phase has translational order in only one direction, and breaks rotational invariance. From finite-size scaling, we demonstrate that both the transitions are first order in nature. We also show that the simpler system with only one type of particle undergoes a single discontinuous phase transition from a disordered phase to a solidlike sublattice phase with an increasing density of particles.