Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPhase diagram of self-assembled rigid rods on two-dimensional lattices: Theory and Monte Carlo simulations
103
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
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.
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.
Numerical simulations and finite-size scaling analysis have been carried out to study the percolation behavior of straight rigid rods of length $k$ ($k$-mers) on two-dimensional square lattices. The $k$-mers, containing $k$ identical units (each one occupying a lattice site), were adsorbed at equilibrium on the lattice. The process was monitored by following the probability $R_{L,k}(theta)$ that a lattice composed of $L times L$ sites percolates at a concentration $theta$ of sites occupied by particles of size $k$. A nonmonotonic size dependence was observed for the percolation threshold, which decreases for small particles sizes, goes through a minimum, and finally asymptotically converges towards a definite value for large segments. This striking behavior has been interpreted as a consequence of the isotropic-nematic phase transition occurring in the system for large values of $k$. Finally, the universality class of the model was found to be the same as for the random percolation model.
In simple ferromagnetic quantum Ising models characterized by an effective double-well energy landscape the characteristic tunneling time of path-integral Monte Carlo (PIMC) simulations has been shown to scale as the incoherent quantum-tunneling time, i.e., as $1/Delta^2$, where $Delta$ is the tunneling gap. Since incoherent quantum tunneling is employed by quantum annealers (QAs) to solve optimization problems, this result suggests there is no quantum advantage in using QAs w.r.t. quantum Monte Carlo (QMC) simulations. A counterexample is the recently introduced shamrock model, where topological obstructions cause an exponential slowdown of the PIMC tunneling dynamics with respect to incoherent quantum tunneling, leaving the door open for potential quantum speedup, even for stoquastic models. In this work, we investigate the tunneling time of projective QMC simulations based on the diffusion Monte Carlo (DMC) algorithm without guiding functions, showing that it scales as $1/Delta$, i.e., even more favorably than the incoherent quantum-tunneling time, both in a simple ferromagnetic system and in the more challenging shamrock model. However a careful comparison between the DMC ground-state energies and the exact solution available for the transverse-field Ising chain points at an exponential scaling of the computational cost required to keep a fixed relative error as the system size increases.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
