No Arabic abstract
Using Monte Carlo simulations and finite-size scaling analysis, the critical behavior of self-assembled rigid rods on triangular and honeycomb lattices at intermediate density has been studied. The system is composed of monomers with two attractive (sticky) poles that, by decreasing temperature or increasing density, polymerize reversibly into chains with three allowed directions and, at the same time, undergo a continuous isotropic-nematic (IN) transition. The determination of the critical exponents, along with the behavior of Binder cumulants, indicate that the IN transition belongs to the q=1 Potts universality class.
Monte Carlo simulations and finite-size scaling analysis have been carried out to study the critical behavior in a two-dimensional system of particles with two bonding sites that, by decreasing temperature or increasing density, polymerize reversibly into chains with discrete orientational degrees of freedom and, at the same time, undergo a continuous isotropic-nematic (IN) transition. A complete phase diagram was obtained as a function of temperature and density. The numerical results were compared with mean field (MF) and real space renormalization group (RSRG) analytical predictions about the IN transformation. While the RSRG approach supports the continuous nature of the transition, the MF solution predicts a first-order transition line and a tricritical point, at variance with the simulation results.
The critical behavior of self-assembled rigid rods on a square lattice was recently reinvestigated by Almarza et al. [Phys. Rev. E 82, 061117 (2010)]. Based on the Binder cumulants and the value of the critical exponent of the correlation length, the authors found that the isotropic-nematic phase transition occurring in the system is in the two-dimensional Ising universality class. This conclusion contrasts with that of a previous study [Lopez et al., Phys. Rev. E 80, 040105 (R) (2009)] which indicates that the transition at intermediate density belongs to the q = 1 Potts universality class. Almarza et al. attributed the discrepancy to the use of the density as the control parameter by Lopez et al. The present work shows that this suggestion is not sufficient, and that the discrepancy arises solely from the use of different statistical ensembles. Finally, the necessity of making corrections to the scaling functions in the canonical ensemble is discussed.
The statistical thermodynamics of straight rigid rods of length $k$ on triangular lattices was developed on a generalization in the spirit of the lattice-gas model and the classical Guggenheim-DiMarzio approximation. In this scheme, the Helmholtz free energy and its derivatives were written in terms of the order parameter $delta$, which characterizes the nematic phase occurring in the system at intermediate densities. Then, using the principle of minimum free energy with $delta$ as a parameter, the main adsorption properties were calculated. Comparisons with Monte Carlo simulations and experimental data were performed in order to evaluate the reaches and limitations of the theoretical model.
On Archimedean lattices, the Ising model exhibits spontaneous ordering. Three examples of these lattices of the majority-vote model with noise are considered and studied through extensive Monte Carlo simulations. The order/disorder phase transition is observed in this system. The calculated values of the critical noise parameter are q_c=0.089(5), q_c=0.078(3), and q_c=0.114(2) for honeycomb, Kagome and triangular lattices, respectively. The critical exponents beta/nu, gamma/nu and 1/nu for this model are 0.15(5), 1.64(5), and 0.87(5); 0.14(3), 1.64(3), and 0.86(6); 0.12(4), 1.59(5), and 1.08(6) for honeycomb, Kagome and triangular lattices, respectively. These results differs from the usual Ising model results and the majority-vote model on so-far studied regular lattices or complex networks. The effective dimensionalities of the system D_{eff}= 1.96(5) (honeycomb), D_{eff} =1.92(4) (Kagome), and D_{eff}= 1.83(5) (triangular) for these networks are just compatible to the embedding dimension two.
The behavior of two-dimensional Ising spin glasses at the multicritical point on triangular and honeycomb lattices is investigated, with the help of finite-size scaling and conformal-invariance concepts. We use transfer-matrix methods on long strips to calculate domain-wall energies, uniform susceptibilities, and spin-spin correlation functions. Accurate estimates are provided for the location of the multicritical point on both lattices, which lend strong support to a conjecture recently advanced by Takeda, Sasamoto, and Nishimori. Correlation functions are shown to obey rather strict conformal-invariance requirements, once suitable adaptations are made to account for geometric aspects of the transfer-matrix description of triangular and honeycomb lattices. The universality class of critical behavior upon crossing the ferro-paramagnetic phase boundary is probed, with the following estimates for the associated critical indices: $ u=1.49(2)$, $gamma=2.71(4)$, $eta_1= 0.183(3)$, distinctly different from the percolation values.