Solution of an associating lattice-gas model with density anomaly on a Husimi lattice

Semianalytical solutions of Ising-like and Potts-like magnetic polymers on the Bethe lattice

We study magnetic polymers, defined as self-avoiding walks where each monomer i carries a "spin" s_{i} and interacts with its first neighbor monomers, let us say j, via a coupling constant J(s_{i},s_{j}). Ising-like [s_{i}=±1, with J(s_{i},s_{j})=ɛs_{i}s_{j}] and Potts-like [s_{i}=1,...,q, with J(s_{i},s_{j})=ɛ(s_{i})δ(s_{i},s_{j})] models are investigated. Some particular cases of these systems have recently been studied in the continuum and on regular lattices and are related to interesting applications. Here, we solve these models on Bethe lattices of branching number σ, focusing on the ferromagnetic case in zero external magnetic field. In most cases, the phase diagrams present a nonpolymerized (NP) and two polymerized phases: a paramagnetic (PP) and a ferromagnetic (FP) one. However, quite different thermodynamic properties are found depending on q in the Potts-like polymers and on whether one uses the Ising or Potts coupling in the two-state systems. For instance, for the standard Potts model [where ɛ(1)=⋯=ɛ(q)] with q=2, beyond the θ-point (where the critical and discontinuous NP-PP transition lines meet), a second tricritical point connecting a critical and a discontinuous transition line between the PP-FP phases is found in the system. A triple point where the NP-PP, NP-FP, and PP-FP first-order transition lines meet is also present in the phase diagram. For q≥3 the PP-FP transition is always discontinuous, and the scenario with the triple and θ points appears for q≤6. Interestingly, for q≥7, as well as for the Ising-like model the θ point becomes metastable and the critical NP-PP transition line ends at a critical end point (CEP), where it meets the NP-FP and PP-FP coexistence lines. Importantly, these results indicate that when q≤6 the spin ordering transition is preceded by the polymer collapse transition, whereas for q≥7 and in the Ising case these transitions happen together at the CEPs. Some interesting nonstandard Potts models are also studied, such as the lattice version of the model for epigenetic marks in the chromatin introduced by Michieletto et al. [Phys. Rev. X 6, 041047 (2016)2160-330810.1103/PhysRevX.6.041047]. In addition, the solution of the dilute Ising and dilute Potts models on the Bethe lattice are also presented here, since they are important to understand the PP-FP transitions.

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.

Phase diagram and critical properties of a two-dimensional associating lattice gas

We revisit the associating lattice gas (ALG) introduced by Henriques et al. [Phys. Rev. E 71, 031504 (2005)PLEEE81539-375510.1103/PhysRevE.71.031504] in its symmetric version. In this model, defined on the triangular lattice, interaction between molecules occupying nearest-neighbor sites depends on their relative orientation, mimicking the formation of hydrogen bonds in network-forming fluids. Although all previous studies of this model agree that it has a disordered fluid (DF), a low-density liquid (LDL), and a high-density liquid (HDL) phase, quite different forms have been reported for its phase diagram. Here, we present a thorough investigation of its phase behavior using both transfer matrix calculations and Monte Carlo (MC) simulations, along with finite-size scaling extrapolations. Results in striking agreement are found using these methods. The critical point associated with the DF-HDL transition at full occupancy, identified by Furlan et al. [Phys. Rev. E 100, 022109 (2019)2470-004510.1103/PhysRevE.100.022109] is shown to be one terminus of a critical line separating these phases. In opposition to previous simulation studies, we find that the transition between the DF and LDL phases is always discontinuous, similar to the LDL-HDL transition. The associated coexistence lines meet at the point where the DF-HDL critical line ends, making it a critical-end-point. Overall, the form of the phase diagram observed in our simulations is very similar to that found in the exact solution of the model on a Husimi lattice. Our results confirm that, despite the existence of some waterlike anomalies in this model, it is unable to reproduce key features of the phase behavior of liquid water.

Fluid-fluid demixing and density anomaly in a ternary mixture of hard spheres

We report the grand-canonical solution of a ternary mixture of discrete hard spheres defined on a Husimi lattice built with cubes, which provides a mean-field approximation for this system on the cubic lattice. The mixture is composed of pointlike particles (0NN) and particles which exclude up to their first (1NN) and second neighbors (2NN), with activities z_{0}, z_{1}, and z_{2}, respectively. Our solution reveals a very rich thermodynamic behavior, with two solid phases associated with the ordering of 1NN (S1) or 2NN particles (S2), and two fluid phases, one being regular (RF) and the other characterized by a dominance of 0NN particles (F0 phase). However, in most part of the phase diagram these fluid (F) phases are indistinguishable. Discontinuous transitions are observed between all the four phases, yielding several coexistence surfaces in the system, among which a fluid-fluid and a solid-solid demixing surface. The former one is limited by a line of critical points and a line of triple points (where the phases RF-F0-S2 coexist), both meeting at a special point, after which the fluid-fluid coexistence becomes metastable. Another line of triple points is found, connecting the F-S1, F-S2, and S1-S2 coexistence surfaces. A critical F-S1 surface is also observed meeting the F-S1 coexistence one at a line of tricritical points. Furthermore, a thermodynamic anomaly characterized by minima in isobaric curves of the total density of particles is found, yielding three surfaces of minimal density in the activity space, depending on which activity is kept fixed during its calculation.

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.

Order-disorder transition in a two-dimensional associating lattice gas

We study an associating lattice gas (ALG) using Monte Carlo simulation on the triangular lattice and semianalytical solutions on Husimi lattices. In this model, the molecules have an orientational degree of freedom and the interactions depend on the relative orientations of nearest-neighbor molecules, mimicking the formation of hydrogen bonds. We focus on the transition between the high-density liquid (HDL) phase and the isotropic phase in the limit of full occupancy, corresponding to chemical potential μ→∞, which has not yet been studied systematically. Simulations yield a continuous phase transition at τ_{c}=k_{B}T_{c}/γ=0.4763(1) (where -γ is the bond energy) between the low-temperature HDL phase, with a nonvanishing mean orientation of the molecules, and the high-temperature isotropic phase. Results for critical exponents and the Binder cumulant indicate that the transition belongs to the three-state Potts model universality class, even though the ALG Hamiltonian does not have the full permutation symmetry of the Potts model. In contrast with simulation, the Husimi lattice analyses furnish a discontinuous phase transition, characterized by a discontinuity of the nematic order parameter. The transition temperatures (τ_{c}=0.51403 and 0.51207 for trees built with triangles and hexagons, respectively) are slightly higher than that found via simulation. Since the Husimi lattice studies show that the ALG phase diagram features a discontinuous isotropic-HDL line for finite μ, three possible scenarios arise for the triangular lattice. The first is that in the limit μ→∞ the first-order line ends in a critical point; the second is a change in the nature of the transition at some finite chemical potential; the third is that the entire line is one of continuous phase transitions. Results from other ALG models and the fact that mean-field approximations show a discontinuous phase transition for the three-state Potts model (known to possess a continuous transition) lends some weight to the third alternative.

Grand-canonical solution of semiflexible self-avoiding trails on the Bethe lattice

We consider a model of semiflexible interacting self-avoiding trails (sISATs) on a lattice, where the walks are constrained to visit each lattice edge at most once. Such models have been studied as an alternative to the self-attracting self-avoiding walks (SASAWs) to investigate the collapse transition of polymers, with the attractive interactions being on site as opposed to nearest-neighbor interactions in SASAWs. The grand-canonical version of the sISAT model is solved on a four-coordinated Bethe lattice, and four phases appear: non-polymerized (NP), regular polymerized (P), dense polymerized (DP), and anisotropic nematic (AN), the last one present in the phase diagram only for sufficiently stiff chains. The last two phases are dense, in the sense that all lattice sites are visited once in the AN phase and twice in the DP phase. In general, critical NP-P and DP-P transition surfaces meet with a NP-DP coexistence surface at a line of bicritical points. The region in which the AN phase is stable is limited by a discontinuous critical transition to the P phase, and we study this somewhat unusual transition in some detail. In the limit of rods, where the chains are totally rigid, the P phase is absent and the three coexistence lines (NP-AN, AN-DP, and NP-DP) meet at a triple point, which is the endpoint of the bicritical line.

Nature of the collapse transition in interacting self-avoiding trails

We study the interacting self-avoiding trail (ISAT) model on a Bethe lattice of general coordination q and on a Husimi lattice built with squares and coordination q=4. The exact grand-canonical solutions of the model are obtained, considering that up to K monomers can be placed on a site and associating a weight ω_{i} with an i-fold visited site. Very rich phase diagrams are found with nonpolymerized, regular polymerized, and dense polymerized phases separated by lines (or surfaces) of continuous and discontinuous transitions. For a Bethe lattice with q=4 and K=2, the collapse transition is identified with a bicritical point and the collapsed phase is associated with the dense polymerized (solidlike) phase instead of the regular polymerized (liquidlike) phase. A similar result is found for the Husimi lattice, which may explain the difference between the collapse transition for ISATs and for interacting self-avoiding walks on the square lattice. For q=6 and K=3 (studied on the Bethe lattice only), a more complex phase diagram is found, with two critical planes and two coexistence surfaces, separated by two tricritical and two critical end-point lines meeting at a multicritical point. The mapping of the phase diagrams in the canonical ensemble is discussed and compared with simulational results for regular lattices.

Phase transitions of the p-spin model on pure Husimi lattices

We consider the p-spin model with spin 1/2 on all pure Husimi lattices. Using an effective representation of the recursion relations, the phase transitions of the model on all pure Husimi lattices are investigated. First, the nonexistence of the second order phase transitions in the model on all pure Husimi lattices is proven exactly. Then the existence and properties of the first order phase transitions in a zero external magnetic field are studied in detail. An implicit polynomial equation for determining temperatures below which the paramagnetic and ferromagnetic phases coexist in the model on all pure Husimi lattices is found. In addition, an implicit equation for exactly determining the transition temperatures of the first order phase transitions in a zero external magnetic field on all pure Husimi lattices is derived and discussed.

Entropy of polydisperse chains: solution on the Husimi lattice

We consider the entropy of polydisperse chains placed on a lattice. In particular, we study a model for equilibrium polymerization, where the polydispersity is determined by two activities, for internal and endpoint monomers of a chain. We solve the problem exactly on a Husimi lattice built with squares and with arbitrary coordination number, obtaining an expression for the entropy as a function of the density of monomers and mean molecular weight of the chains. We compare this entropy with the one for the monodisperse case, and find that the excess of entropy due to polydispersity is identical to the one obtained for the one-dimensional case. Finally, we obtain a distribution of molecular weights with a rather complex behavior, but which becomes exponential for very large mean molecular weight of the chains, as required by scaling properties, which should apply in this limit.

Solution on the Bethe lattice of a hard core athermal gas with two kinds of particles

Athermal lattice gases of particles with first neighbor exclusion have been studied for a long time as simple models exhibiting a fluid-solid transition. At low concentration the particles occupy randomly both sublattices, but as the concentration is increased one of the sublattices is occupied preferentially. Here, we study a mixed lattice gas with excluded volume interactions only in the grand-canonical formalism with two kinds of particles: small ones, which occupy a single lattice site and large ones, which, when placed on a site, do not allow other particles to occupy its first neighbors also. We solve the model on a Bethe lattice of arbitrary coordination number q. In the parameter space defined by the activities of both particles, at low values of the activity of small particles (z(1)) we find a continuous transition from the fluid to the solid phase as the activity of large particles (z(2)) is increased. At higher values of z(1) the transition becomes discontinuous, both regimes are separated by a tricritical point. The critical line has a negative slope at z(1) = 0 and displays a minimum before reaching the tricritical point, so that a re-entrant behavior is observed for constant values of z(2) in the region of low density of small particles. The isobaric curves of the total density of particles as a function of the density or the activity of small particles show a minimum in the fluid phase.