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

Exclusion statistics for structured particles on topologically correlated states. II. Multicomponent lattice gases

Statistical thermodynamics of particles having a spectrum of topological correlated states and observing statistical exclusion is developed to describe mixtures of species of arbitrary size and shape. A generalized statistical distribution is obtained through a configuration space ansatz recently introduced for single species accounting for the multiple exclusion statistical phenomena, where spatially correlated particle states can be simultaneously excluded by more than one particle. Statistical exclusion on a correlated states spectrum is characterized by exclusion statistical parameters β_{cij} which are self-consistently determined within the multiple exclusion from thermodynamic boundary conditions. Self-exclusion and cross-exclusion frequency functions e_{ij}(n) and average cumulative exclusion functions G_{ij}(n) are introduced to characterize the state exclusion spectrum as density varies. Haldane's statistics and Wu's distribution for statistically independent excluding species are recovered in the limit of uncorrelated states for single species as well as for mixtures of self- and cross-excluding species with constant mutual statistical exclusion. The multiple exclusion statistics formalism is applied to the k-mer problem on a square lattice rationalized as a mixture of two differently oriented self-excluding and cross-excluding pseudospecies. An isotropic-nematic and a high-density nematic-isotropic (disordered) phase transitions is predicted only for k≥7. The isotropic-nematic transition is continuum as expected, but the high-density transition results in a first-order one. The formalism provides phase coexistence lines and the chemical potential dependence of the low- and high-density branches in the nematic regime. The theoretical approach to lattice gases presented in this work offers a unique general framework applicable to mixtures of entropy-complex lattice gases. From this framework, k-mer phase transitions can be reproduced, and significant configuration features can be derived from the state exclusion spectrum functions.

Exclusion statistics for structured particles on topologically correlated states. I. Single species lattice gases

A statistical thermodynamics description of particles having a set of spatially correlated states with statistical exclusion is developed. A general approximation for the density of states is presented from a state-counting ansatz recently introduced accounting for the multiple state exclusion statistical phenomena as a consequence of state spatial correlations. The multiple exclusion statistics is characterized by an exclusion correlation constant g_{c} which is consistently determined within the formalism from proper thermodynamic limits. The analytical form of g_{c} is given in terms of the Lambert function from the particle-lattice geometry. A generalized statistical distribution is obtained reducing to Haldane's statistics and Wu's distribution in the limiting case of particles on a set of spatially uncorrelated states. The problem of hard rods (k-mers) on a square lattice is studied with this formalism. From the entropy density dependence of the isotropic (I) and fully oriented nematic (N) phases, the approximation predicts two transitions, I→N and high-coverage N→I (disordered), only for k≥7 with the entropy at saturation matching to the known value from a Monte Carlo (MC) simulation. Critical coverage of both transitions is given for k=7 to k=20 in the first and second orders of approximations, in qualitative and quantitative agreement with results from MC simulations. State exclusion frequency e(n) and exclusion average G(n) functions are introduced and given in terms of the chemical potential to obtain a thermodynamic characterization of the state exclusion evolution on density. Results of chemical potential and state exclusion are shown for ideal lattice gases of k-mers, squares, and rectangles on a square lattice. Analytical results are compared with fast-relaxation grand canonical MC simulations.

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.

Freezing phase transition in hard-core lattice gases on the triangular lattice with exclusion up to seventh next-nearest neighbor

Hard-core lattice-gas models are minimal models to study entropy-driven phase transitions. In the k-nearest-neighbor lattice gas, a particle excludes all sites up to the kth next-nearest neighbors from being occupied by another particle. As k increases from one, it extrapolates from nearest-neighbor exclusion to the hard-sphere gas. In this paper we study the model on the triangular lattice for k≤7 using a flat histogram algorithm that includes cluster moves. Earlier studies focused on k≤3. We show that for 4≤k≤7, the system undergoes a single phase transition from a low-density fluid phase to a high-density sublattice-ordered phase. Using partition function zeros and nonconvexity properties of the entropy, we show that the transitions are discontinuous. The critical chemical potential, coexistence densities, and critical pressure are determined accurately.

Triangles on a triangular lattice: Insights into self-assembly in two dimensions driven by shape complementarity

A series of models for reversible filling of a triangular lattice with equilateral triangles has been developed and investigated. There are eight distinct models that vary in the set of prohibitions. In zeroth approximation, these models allow one to estimate the influence of the particles' shape and complementarity of their pair configurations on the self-assembly of dense monolayers formed by reversible filling. The most symmetrical models were found to be equivalent to hard-disk models on the hexagonal lattice. When any contact of hard triangles by vertices is prohibited, the dense monolayers are disordered, and their entropy tends to the constant. If only one pair configuration is prohibited, the close-packed layer appears through the continuous phase transition. In other cases, the weak first-order transition resulting in the self-assembly of close-packed layers is observed.

Condensation and Crystal Nucleation in a Lattice Gas with a Realistic Phase Diagram

We reconsider model II of Orban et al. (J. Chem. Phys. 1968, 49, 1778−1783), a two-dimensional lattice-gas system featuring a crystalline phase and two distinct fluid phases (liquid and vapor). In this system, a particle prevents other particles from occupying sites up to third neighbors on the square lattice, while attracting (with decreasing strength) particles sitting at fourth- or fifth-neighbor sites. To make the model more realistic, we assume a finite repulsion at third-neighbor distance, with the result that a second crystalline phase appears at higher pressures. However, the similarity with real-world substances is only partial: Upon closer inspection, the alleged liquid−vapor transition turns out to be a continuous (albeit sharp) crossover, even near the putative triple point. Closer to the standard picture is instead the freezing transition, as we show by computing the free-energy barrier relative to crystal nucleation from the “liquid”.

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.

Hard core lattice gas with third next-nearest neighbor exclusion on triangular lattice: One or two phase transitions?

We obtain the phase diagram of the hard core lattice gas with third nearest neighbor exclusion on the triangular lattice using Monte Carlo simulations that are based on a rejection-free flat histogram algorithm. In a recent paper [Darjani et al., J. Chem. Phys. 151, 104702 (2019)], it was claimed that the lattice gas with third nearest neighbor exclusion undergoes two phase transitions with increasing density with the phase at intermediate densities exhibiting hexatic order with continuously varying exponents. Although a hexatic phase is expected when the exclusion range is large, it has not been seen earlier in hard core lattice gases with short range exclusion. In this paper, by numerically determining the entropies for all densities, we show that there is only a single phase transition in the system between a low-density fluid phase and a high density ordered sublattice phase and that a hexatic phase is absent. The transition is shown to be first order in nature, and the critical parameters are determined accurately.

Rejection-free cluster Wang-Landau algorithm for hard-core lattice gases

We introduce a rejection-free, flat histogram, cluster algorithm to determine the density of states of hard-core lattice gases. We show that the algorithm is able to efficiently sample low entropy states that are usually difficult to access, even when the excluded volume per particle is large. The algorithm is based on simultaneously evaporating all the particles in a strip and reoccupying these sites with a new appropriately chosen configuration. We implement the algorithm for the particular case of the hard-core lattice gas in which the first k next-nearest neighbors of a particle are excluded from being occupied. It is shown that the algorithm is able to reproduce the known results for k=1,2,3 both on the square and cubic lattices. We also show that, in comparison, the corresponding flat histogram algorithms with either local moves or unbiased cluster moves are less accurate and do not converge as the system size increases.

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.

Husimi-lattice solutions and the coherent-anomaly-method analysis for hard-square lattice gases

Although lattice gases composed of particles preventing up to their kth nearest neighbors from being occupied (the kNN models) have been widely investigated in the literature, the location and the universality class of the fluid-columnar transition in the 2NN model on the square lattice are still a topic of debate. Here, we present grand-canonical solutions of this model on Husimi lattices built with diagonal square lattices, with 2L(L+1) sites, for L⩽7. The systematic sequence of mean-field solutions confirms the existence of a continuous transition in this system, and extrapolations of the critical chemical potential μ_{2,c}(L) and particle density ρ_{2,c}(L) to L→∞ yield estimates of these quantities in close agreement with previous results for the 2NN model on the square lattice. To confirm the reliability of this approach, we employ it also for the 1NN model, where very accurate estimates for the critical parameters μ_{1,c} and ρ_{1,c}-for the fluid-solid transition in this model on the square lattice-are found from extrapolations of data for L⩽6. The nonclassical critical exponents for these transitions are investigated through the coherent anomaly method (CAM), which in the 1NN case yields β and ν differing by at most 6% from the expected Ising exponents. For the 2NN model, the CAM analysis is somewhat inconclusive, because the exponents sensibly depend on the value of μ_{2,c} used to calculate them. Notwithstanding, our results suggest that β and ν are considerably larger than the Ashkin-Teller exponents reported in numerical studies of the 2NN system.

Glassy dynamics and equilibrium state on the honeycomb lattice: Role of surface diffusion and desorption on surface crowding

The phase behavior and adsorption kinetics of hard-core particles on a honeycomb lattice are studied by means of random sequential adsorption with surface diffusion. We concentrate on reversible adsorption by introducing a desorption process into our previous model and varying the equilibrium rate constant as a control parameter. We find that an exact prediction of the temporal evolution of fractional surface coverage and the surface pressure dynamics of reversible adsorption can be achieved by use of the blocking function of a system with irreversible adsorption of highly mobile particles. For systems out of equilibrium we observe several features of glassy dynamics, such as slow relaxation dynamics, the memory effect, and aging. In particular, the analysis of our system in the limit of small desorption probability shows simple aging behavior with a power-law decay. A detailed discussion of Gibbs adsorption isotherm for nonequilibrium adsorption is given, which exhibits a hysteresis between this system and its equilibrium counterpart.

Global phase diagram of Coulomb-interacting anisotropic Weyl semimetal with disorder

Taking into account the interplay between the disorder and Coulomb interaction, the phase diagram of three-dimensional anisotropic Weyl semimetal is studied by renormalization group (RG) theory. Weak disorder is irrelevant in anisotropic Weyl semimetal, while the disorder becomes relevant and drives a quantum phase transition (QPT) from semimetal to compressible diffusive metal (CDM) phases if the disorder strength is larger than a critical value. The long-range Coulomb interaction is irrelevant in clean anisotropic Weyl semimetal. However, interestingly, we find that the long-range Coulomb interaction exerts a dramatic influence on the critical disorder strength for phase transition to CDM. Specifically, the critical disorder strength can receive a prominent change even though an arbitrarily weak Coulomb interaction is included. This novel behavior is closely related to the anisotropic screening effect of Coulomb interaction, and essentially results from the specifical energy dispersion of the fermion excitations in anisotropic Weyl semimetal. The theoretical results are helpful for understanding the physical properties of the candidates of anisotropic Weyl semimetal, such as pressured BiTeI, and some other related materials.

Phase transitions in hard-core lattice gases on the honeycomb lattice

We study lattice gas systems on the honeycomb lattice where particles exclude neighboring sites up to order k (k=1,...,5) from being occupied by another particle. Monte Carlo simulations were used to obtain phase diagrams and characterize phase transitions as the system orders at high packing fractions. For systems with first-neighbors exclusion (1NN), we confirm previous results suggesting a continuous transition in the two-dimensional Ising universality class. Exclusion up to second neighbors (2NN) lead the system to a two-step melting process where, first, a high-density columnar phase undergoes a first-order phase transition with nonstandard scaling to a solidlike phase with short-range ordered domains and, then, to fluidlike configurations with no sign of a second phase transition. 3NN exclusion, surprisingly, shows no phase transition to an ordered phase as density is increased, staying disordered even to packing fractions up to 0.98. The 4NN model undergoes a continuous phase transition with critical exponents close to the three-state Potts model. The 5NN system undergoes two first-order phase transitions, both with nonstandard scaling. We, also, propose a conjecture concerning the possibility of more than one phase transition for systems with exclusion regions further than 5NN based on geometrical aspects of symmetries.

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.

Thermodynamic behavior of binary mixtures of hard spheres: Semianalytical solutions on a Husimi lattice built with cubes

We study binary mixtures of hard particles, which exclude up to their kth nearest neighbors (kNN) on the simple cubic lattice and have activities z_{k}. In the first model analyzed, point particles (0NN) are mixed with 1NN ones. The grand-canonical solution of this model on a Husimi lattice built with cubes unveils a phase diagram with a fluid and a solid phase separated by a continuous and a discontinuous transition line which meet at a tricritical point. A density anomaly, characterized by minima in isobaric curves of the total density of particles against z_{0} (or z_{1}), is also observed in this system. Overall, this scenario is identical to the one previously found for this model when defined on the square lattice. The second model investigated consists of the mixture of 1NN particles with 2NN ones. In this case, a very rich phase behavior is found in its Husimi lattice solution, with two solid phases-one associated with the ordering of 1NN particles (S1) and the other with the ordering of 2NN ones (S2)-beyond the fluid (F) phase. While the transitions between F-S2 and S1-S2 phases are always discontinuous, the F-S1 transition is continuous (discontinuous) for small (large) z_{2}. The critical and coexistence F-S1 lines meet at a tricritical point. Moreover, the coexistence F-S1,F-S2, and S1-S2 lines meet at a triple point. Density anomalies are absent in this case.

Phase diagram of a system of hard cubes on the cubic lattice

We study the phase diagram of a system of 2×2×2 hard cubes on a three-dimensional cubic lattice. Using Monte Carlo simulations, we show that the system exhibits four different phases as the density of cubes is increased: disordered, layered, sublattice ordered, and columnar ordered. In the layered phase, the system spontaneously breaks up into parallel slabs of size 2×L×L where only a very small fraction cubes do not lie wholly within a slab. Within each slab, the cubes are disordered; translation symmetry is thus broken along exactly one principal axis. In the solidlike sublattice-ordered phase, the hard cubes preferentially occupy one of eight sublattices of the cubic lattice, breaking translational symmetry along all three principal directions. In the columnar phase, the system spontaneously breaks up into weakly interacting parallel columns of size 2×2×L, where only a very small fraction cubes do not lie wholly within a column. Within each column, the system is disordered, and thus translational symmetry is broken only along two principal directions. Using finite-size scaling, we show that the disordered-layered phase transition is continuous, while the layered-sublattice and sublattice-columnar transitions are discontinuous. We construct a Landau theory written in terms of the layering and columnar order parameters which is able to describe the different phases that are observed in the simulations and the order of the transitions. Additionally, our results near the disordered-layered transition are consistent with the O(3) universality class perturbed by cubic anisotropy as predicted by the Landau theory.