Multiple phase transitions in extended hard-core lattice gas models in two dimensions

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.

Weakly first-order transition in an athermal lattice gas

We investigate the phase behavior of a two-dimensional athermal lattice gas in which every hard-core particle can have two or fewer nearest neighboring occupied sites on the square lattice. The ground state and close packing density are determined and it is found that at large chemical potential the model undergoes an ordering phase transition with preferential sublattice occupation. Although near the transition point the particle density and entropy exhibit an apparent discontinuity, we find that the order parameter and fluctuations of thermodynamic quantities do not scale with the system volume. These paradoxical results are reconciled by analyzing the size-dependent flow of the thermal exponent by phenomenological renormalization and the curve-crossing method, which lead to a weakly first-order phase transition scenario.

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.

First-order phase transition in a two dimensional BM3 model

The phase behavior of a Biroli–Mézard model on the two dimensional square lattice in which hard-core particles can have at most three nearest neighboring occupied sites is investigated by means of grand-canonical Monte Carlo simulations. Finite-size scaling analysis of relevant thermodynamic quantities obtained via the histogram reweighting technique reveals that at high-density, the model undergoes a first-order phase transition with preferential sublattice occupation to a crystal phase with enantiomorph ground state configurations, in close analogy to the hard-core lattice gas with the exclusion range extended up to the third shell of nearest neighbors.

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.

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.

SuSMoST: Surface Science Modeling and Simulation Toolkit

We present to the scientific community the Surface Science Modeling and Simulation Toolkit (SuSMoST), which includes a number of utilities and implementations of statistical physics algorithms and models. With SuSMoST it is possible to predict or explain the structure and thermodynamic properties of adsorption layers. SuSMoST automatically builds formal graph and tensor-network models based on atomic description of adsorption complexes and helps to do ab initio calculations of interactions between adsorbed species. Using methods of various nature SuSMoST generates representative samples of adsorption layers and computes its thermodynamic quantities such as mean energy, coverage, density, and heat capacity. From these data one can plot phase diagrams of adsorption systems, assess thermal stability of self-assembled structures, simulate thermal desorption spectra, and so forth.

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.

Tensor renormalization group study of hard-disk models on a triangular lattice

High accuracy and performance of the tensor renormalization group (TRG) method have been demonstrated for the model of hard disks on a triangular lattice. We considered a sequence of models with disk diameter ranging from a to 2sqrt[3]a, where a is the lattice constant. Practically, these models are good for approximate description of thermodynamics properties of molecular layers on crystal surfaces. Theoretically, it is interesting to analyze if and how this sequence converges to the continuous model of hard disks. The dependencies of the density and heat capacity on the chemical potential were calculated with TRG and transfer-matrix (TM) methods. We benchmarked accuracy and performance of the TRG method comparing it with TM method and with exact result for the model with nearest-neighbor exclusions (1NN). The TRG method demonstrates good convergence and turns out to be superior over TM with regard to considered models. Critical values of chemical potential (μ_{c}) have been computed for all models. For the model with next-nearest-neighbor exclusions (2NN) the TRG and TM produce consistent results (μ_{c}=1.75587 and μ_{c}=1.75398 correspondingly) that are also close to earlier Monte Carlo estimation by Zhang and Deng. We found that 3NN and 5NN models shows the first-order phase transition, with close values of μ_{c} (μ_{c}=4.4488 for 3NN and 4.4<μ_{c}<4.5 for 5NN). The 4NN model demonstrates continuous yet rapid phase transition with 2.65<μ_{c}<2.7.

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.

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.

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

A mixture of hard squares, dimers, and vacancies on a square lattice is known to undergo a transition from a low-density disordered phase to a high-density columnar ordered phase. Along the fully packed square-dimer line, the system undergoes a Kosterliz-Thouless-type transition to a phase with power law correlations. We estimate the phase boundary separating the ordered and disordered phases by calculating the interfacial tension between two differently ordered phases within two different approximation schemes. The analytically obtained phase boundary is in good agreement with Monte Carlo simulations.

Transfer-matrix study of a hard-square lattice gas with two kinds of particles and density anomaly

Using transfer matrix and finite-size scaling methods, we study the thermodynamic behavior of a lattice gas with two kinds of particles on the square lattice. Only excluded volume interactions are considered, so that the model is athermal. Large particles exclude the site they occupy and its four first neighbors, while small particles exclude only their site. Two thermodynamic phases are found: a disordered phase where large particles occupy both sublattices with the same probability and an ordered phase where one of the two sublattices is preferentially occupied by them. The transition between these phases is continuous at small concentrations of the small particles and discontinuous at larger concentrations, both transitions are separated by a tricritical point. Estimates of the central charge suggest that the critical line is in the Ising universality class, while the tricritical point has tricritical Ising (Blume-Emery-Griffiths) exponents. The isobaric curves of the total density as functions of the fugacity of small or large particles display a minimum in the disordered phase.

Columnar order and Ashkin-Teller criticality in mixtures of hard squares and dimers

We show that critical exponents of the transition to columnar order in a mixture of 2×1 dimers and 2×2 hard squares on the square lattice depends on the composition of the mixture in exactly the manner predicted by the theory of Ashkin-Teller criticality, including in the hard-square limit. This result settles the question regarding the nature of the transition in the hard-square lattice gas. It also provides the first example of a polydisperse system whose critical properties depend on composition. Our ideas also lead to some interesting predictions for a class of frustrated quantum magnets that exhibit columnar ordering of the bond energies at low temperature.

Polydispersed rods on the square lattice

We study the grand-canonical solution of a system of hard polydispersed rods placed on the square lattice using transfer matrix and finite-size scaling calculations. We determine the critical line separating an isotropic from a nematic phase. No second transition to a disordered phase is found at high density, contrary to what is observed in the monodispersed case. The estimates of critical exponents and the central charge on the critical line are consistent with the Ising universality class.

Asymptotic behavior of the isotropic-nematic and nematic-columnar phase boundaries for the system of hard rectangles on a square lattice

A system of hard rectangles of size m×mk on a square lattice undergoes three entropy-driven phase transitions with increasing density for large-enough aspect ratio k: first from a low-density isotropic to an intermediate-density nematic phase, second from the nematic to a columnar phase, and third from the columnar to a high-density sublattice phase. In this paper we show, from extensive Monte Carlo simulations of systems with m=1,2, and 3, that the transition density for the isotropic-nematic transition is ≈A(1)/k when k≫1, where A(1) is independent of m. We estimate A(1)=4.80±0.05. Within a Bethe approximation and virial expansion truncated at the second virial coefficient, we obtain A(1)=2. The critical density for the nematic-columnar transition when m=2 is numerically shown to tend to a value less than the full packing density as k(-1) when k→∞. We find that the critical Binder cumulant for this transition is nonuniversal and decreases as k(-1) for k≫1. However, the transition is shown to be in the Ising universality class.