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

Recurrence solution of monomer-polymer models on two-dimensional rectangular lattices

The problem of counting polymer coverings on rectangular lattices is investigated. In this model, a linear rigid polymer covers k adjacent lattice sites such that no two polymers occupy a common site. Those unoccupied lattice sites are considered as monomers. We prove that for a given number of polymers (k-mers), the number of arrangements for the polymers on two-dimensional rectangular lattices satisfies simple recurrence relations. These recurrence relations are quite general and apply for arbitrary polymer length (k) and the width of the lattices (n). The well-studied monomer-dimer problem is a special case of the monomer-polymer model when k=2. It is known the enumeration of monomer-dimer configurations in planar lattices is #P complete. The recurrence relations shown here have the potential for hints for the solution of long-standing problems in this class of computational complexity.

Entropy-driven phases at high coverage adsorption of straight rigid rods on three-dimensional cubic lattices

Combining Monte Carlo simulations and thermodynamic integration method, we study the configurational entropy per site of straight rigid rods of length k (k-mers) adsorbed on three-dimensional (3D) simple cubic lattices. The process is monitored by following the dependence of the lattice coverage θ on the chemical potential μ (adsorption isotherm). Then, we perform the integration of μ(θ) over θ to calculate the configurational entropy per site of the adsorbed phase s(k,θ) as a function of the coverage. Based on the behavior of the function s(k,θ), different phase diagrams are obtained according to the k values: k≤4, disordered phase; k=5,6, disordered and layered-disordered phases; and k≥7, disordered, nematic and layered-disordered phases. In the limit of θ→1 (full coverage), the configurational entropy per site is determined for values of k ranging between 2 and 8. For k≥6, MC data coincide (within the statistical uncertainty) with recent analytical predictions [D. Dhar and R. Rajesh, Phys. Rev. E 103, 042130 (2021)2470-004510.1103/PhysRevE.103.042130] for very large rods. This finding represents the first numerical validation of the expression obtained by Dhar and Rajesh for d-dimensional lattices with d>2. In addition, for k≥5, the values of s(k,θ→1) for simple cubic lattices are coincident with those values reported in [P. M. Pasinetti et al., Phys. Rev. E 104, 054136 (2021)2470-004510.1103/PhysRevE.104.054136] for two-dimensional (2D) square lattices. This is consistent with the picture that at high densities and k≥5, the layered-disordered phase is formed on the lattice. Under these conditions, the system breaks to 2D layers, and the adsorbed phase becomes essentially 2D. The 2D behavior of the fully covered lattice reinforces the conjecture that the large-k behavior of entropy per site is superuniversal, and holds on d-dimensional hypercubical lattices for all d≥2.

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 of rigid k-mers on a square lattice

Using the transfer matrix technique, we estimate the entropy for a gas of rods of sizes equal to k (named k-mers), which cover completely a square lattice. Our calculations were made considering three different constructions, using periodical and helical boundary conditions. One of those constructions, which we call profile method, was based on the calculations performed by Dhar and Rajesh to obtain a lower limit to the entropy of very large chains placed on the square lattice. This method, so far as we know, was never used before to define the transfer matrix, but turned out to be very useful, since it produces matrices with smaller dimensions than those obtained using the usual approach. Our results were obtained for chain sizes ranging from k=2 to k=10 and they are compared with results already available in the literature. In the case of dimers (k=2) our results are compatible with the exact result. For trimers (k=3), recently investigated by Ghosh et al., also our results were compatible, with the same happening for the simulational estimates obtained by Pasinetti et al. in the whole range of rod sizes. Our results are also consistent with the asymptotic expression for the behavior of the entropy as a function of the size k, proposed by Dhar and Rajesh for very large rods (k≫1).

Orientational phase transition in monolayers of multipolar straight rigid rods: The case of 2-thiophene molecule adsorption on the Au(111) surface

Monte Carlo simulations and finite-size scaling theory have been carried out to study the critical behavior and universality for the isotropic-nematic (IN) phase transition in a system of straight rigid pentamers adsorbed on a triangular lattice with polarized nonhomogeneous intermolecular interactions. The model was inspired by the deposition of 2-thiophene molecules over the Au(111) surface, which was previously characterized by experimental techniques and density functional theory. A nematic phase, observed experimentally by the formation of a self-assembled monolayer of parallel molecules, is separated from the isotropic state by a continuous transition occurring at a finite density. The precise determination of the critical exponents indicates that the transition belongs to the three-state Potts universality class. The finite-size scaling analysis includes the study of mutability and diversity. These two quantities are derived from information theory and they have not previously been considered as part of the conventional treatment of critical phenomena for the determination of critical exponents. The results obtained here contribute to the understanding of formation processes of self-assembled monolayers, phase transitions, and critical phenomena from novel compression algorithms for studying mutual information in sequences of data.

Phase transition from nematic to high-density disordered phase in a system of hard rods on a lattice

A system of hard rigid rods of length k on hypercubic lattices is known to undergo two phase transitions when chemical potential is increased: from a low density isotropic phase to an intermediate density nematic phase, and on further increase to a high-density phase with no orientational order. In this paper, we argue that, for large k, the second phase transition is a first-order transition with a discontinuity in density in all dimensions greater than 1. We show that the chemical potential at the transition is ≈kln[k/lnk] for large k, and that the density of uncovered sites drops from a value ≈(lnk)/k^{2} to a value of order exp(-ak), where a is some constant, across the transition. We conjecture that these results are asymptotically exact, in all dimensions d≥2. We also present evidence of coexistence of nematic and disordered phases from Monte Carlo simulations for rods of length 9 on the square lattice.

Entropy of fully packed rigid rods on generalized Husimi trees: A route to the square-lattice limit

Although hard rigid rods (k-mers) defined on the square lattice have been widely studied in the literature, their entropy per site, s(k), in the full-packing limit is only known exactly for dimers (k=2) and numerically for trimers (k=3). Here, we investigate this entropy for rods with k≤7, by defining and solving them on Husimi lattices built with diagonal and regular square-lattice clusters of effective lateral size L, where L defines the level of approximation to the square lattice. Due to an L-parity effect, by increasing L we obtain two systematic sequences of values for the entropies s_{L}(k) for each type of cluster, whose extrapolations to L→∞ provide estimates of these entropies for the square lattice. For dimers, our estimates for s(2) differ from the exact result by only 0.03%, while that for s(3) differs from best available estimates by 3%. In this paper, we also obtain a new estimate for s(4). For larger k, we find that the extrapolated results from the Husimi tree calculations do not lie between the lower and upper bounds established in the literature for s(k). In fact, we observe that, to obtain reliable estimates for these entropies, we should deal with levels L that increase with k. However, it is very challenging computationally to advance to solve the problem for large values of L and for large rods. In addition, the exact calculations on the generalized Husimi trees provide strong evidence for the fully packed phase to be disordered for k≥4, in contrast to the results for the Bethe lattice wherein it is nematic.

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.