1
|
Akritidis M, Fytas NG, Weigel M. Geometric clusters in the overlap of the Ising model. Phys Rev E 2023; 108:044145. [PMID: 37978672 DOI: 10.1103/physreve.108.044145] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2023] [Accepted: 09/22/2023] [Indexed: 11/19/2023]
Abstract
We study the percolation properties of geometrical clusters defined in the overlap space of two statistically independent replicas of a square-lattice Ising model that are simulated at the same temperature. In particular, we consider two distinct types of clusters in the overlap, which we dub soft- and hard-constraint clusters, and which are subsets of the regions of constant spin overlap. By means of Monte Carlo simulations and a finite-size scaling analysis we estimate the transition temperature as well as the set of critical exponents characterizing the percolation transitions undergone by these two cluster types. The results suggest that both soft- and hard-constraint clusters percolate at the critical temperature of the Ising model and their critical behavior is governed by the correlation-length exponent ν=1 found by Onsager. At the same time, they exhibit nonstandard and distinct sets of exponents for the average cluster size and percolation strength.
Collapse
Affiliation(s)
- Michail Akritidis
- Centre for Fluid and Complex Systems, Coventry University, Coventry, CV1 5FB, United Kingdom
| | - Nikolaos G Fytas
- Institut für Physik, Technische Universität Chemnitz, 09107 Chemnitz, Germany
- Department of Mathematical Sciences, University of Essex, Colchester CO4 3SQ, United Kingdom
| | - Martin Weigel
- Institut für Physik, Technische Universität Chemnitz, 09107 Chemnitz, Germany
| |
Collapse
|
2
|
Macêdo ARS, Vasilopoulos A, Akritidis M, Plascak JA, Fytas NG, Weigel M. Two-dimensional dilute Baxter-Wu model: Transition order and universality. Phys Rev E 2023; 108:024140. [PMID: 37723817 DOI: 10.1103/physreve.108.024140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2023] [Accepted: 08/07/2023] [Indexed: 09/20/2023]
Abstract
We investigate the critical behavior of the two-dimensional spin-1 Baxter-Wu model in the presence of a crystal-field coupling Δ with the goal of determining the universality class of transitions along the second-order part of the transition line as one approaches the putative location of the multicritical point. We employ extensive Monte Carlo simulations using two different methodologies: (i) a study of the zeros of the energy probability distribution, closely related to the Fisher zeros of the partition function, and (ii) the well-established multicanonical approach employed to study the probability distribution of the crystal-field energy. A detailed finite-size scaling analysis in the regime of second-order phase transitions in the (Δ,T) phase diagram supports previous claims that the transition belongs to the universality class of the four-state Potts model. For positive values of Δ, we observe the presence of strong finite-size effects, indicative of crossover effects due to the proximity of the first-order part of the transition line. Finally, we demonstrate how a combination of cluster and heat-bath updates allows one to equilibrate larger systems, and we demonstrate the potential of this approach for resolving the ambiguities observed in the regime of Δ≳0.
Collapse
Affiliation(s)
- A R S Macêdo
- Departamento de Física, Instituto de Ciências Exatas, Universidade Federal de Minas Gerais, C.P. 702, Belo Horizonte 65919-050, MG, Brazil
- Instituto Federal do Maranhão - Campus Imperatriz, Imperatriz 65919-050, MA, Brazil
| | - A Vasilopoulos
- Centre for Fluid and Complex Systems, Coventry University, Coventry CV1 5FB, United Kingdom
| | - M Akritidis
- Centre for Fluid and Complex Systems, Coventry University, Coventry CV1 5FB, United Kingdom
| | - J A Plascak
- Departamento de Física, Instituto de Ciências Exatas, Universidade Federal de Minas Gerais, C.P. 702, Belo Horizonte 65919-050, MG, Brazil
- Universidade Federal da Paraíba, Centro de Ciências Exatas e da Natureza - Campus I, Departamento de Física - CCEN Cidade Universitária, João Pessoa 58051-970, PB, Brazil
- Center for Simulational Physics, University of Georgia, Athens, Georgia 30602, USA
| | - N G Fytas
- Centre for Fluid and Complex Systems, Coventry University, Coventry CV1 5FB, United Kingdom
- Department of Mathematical Sciences, University of Essex, Colchester CO4 3SQ, United Kingdom
| | - M Weigel
- Institut für Physik, Technische Universität Chemnitz, 09107 Chemnitz, Germany
| |
Collapse
|
3
|
Vasilopoulos A, Fytas NG, Vatansever E, Malakis A, Weigel M. Universality in the two-dimensional dilute Baxter-Wu model. Phys Rev E 2022; 105:054143. [PMID: 35706204 DOI: 10.1103/physreve.105.054143] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2022] [Accepted: 05/02/2022] [Indexed: 06/15/2023]
Abstract
We study the question of universality in the two-dimensional spin-1 Baxter-Wu model in the presence of a crystal field Δ. We employ extensive numerical simulations of two types, providing us with complementary results: Wang-Landau sampling at fixed values of Δ and a parallelized variant of the multicanonical approach performed at constant temperature T. A detailed finite-size scaling analysis in the regime of second-order phase transitions in the (Δ,T) phase diagram indicates that the transition belongs to the universality class of the four-state Potts model. Previous controversies with respect to the nature of the transition are discussed and attributed to the presence of strong finite-size effects, especially as one approaches the pentacritical point of the model.
Collapse
Affiliation(s)
| | - Nikolaos G Fytas
- Centre for Fluid and Complex Systems, Coventry University, Coventry CV1 5FB, United Kingdom
- Institut für Physik, Technische Universität Chemnitz, D-09107 Chemnitz, Germany
| | - Erol Vatansever
- Centre for Fluid and Complex Systems, Coventry University, Coventry CV1 5FB, United Kingdom
- Department of Physics, Dokuz Eylül University, TR-35160 Izmir, Turkey
| | - Anastasios Malakis
- Centre for Fluid and Complex Systems, Coventry University, Coventry CV1 5FB, United Kingdom
- Department of Physics, University of Athens, Panepistimiopolis, GR-15784 Zografou, Greece
| | - Martin Weigel
- Institut für Physik, Technische Universität Chemnitz, D-09107 Chemnitz, Germany
| |
Collapse
|
4
|
Parisen Toldin F. Finite-size scaling at fixed renormalization-group invariant. Phys Rev E 2022; 105:034137. [PMID: 35428125 DOI: 10.1103/physreve.105.034137] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2022] [Accepted: 03/07/2022] [Indexed: 06/14/2023]
Abstract
Finite-size scaling at fixed renormalization-group invariant is a powerful and flexible technique to analyze Monte Carlo data at a critical point. It consists in fixing a given renormalization-group invariant quantity to a given value, thereby trading its statistical fluctuations with those of a parameter driving the transition. One remarkable feature is the observed significant improvement of statistical accuracy of various quantities, as compared to a standard analysis. We review the method, discussing in detail its implementation, the error analysis, and a previously introduced covariance-based optimization. Comprehensive benchmarks on the Ising model in two and three dimensions show large gains in the statistical accuracy, which are due to cross-correlations between observables. As an application, we compute an accurate estimate of the inverse critical temperature of the improved O(2) ϕ^{4} model on a three-dimensional cubic lattice.
Collapse
Affiliation(s)
- Francesco Parisen Toldin
- Institut für Theoretische Physik und Astrophysik, Universität Würzburg, Am Hubland, D-97074 Würzburg, Germany
| |
Collapse
|
5
|
Torstensson E, Goyal G, Johnning A, Westerlund F, Ambjörnsson T. Combining dense and sparse labeling in optical DNA mapping. PLoS One 2021; 16:e0260489. [PMID: 34843574 PMCID: PMC8629184 DOI: 10.1371/journal.pone.0260489] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2021] [Accepted: 11/10/2021] [Indexed: 11/19/2022] Open
Abstract
Optical DNA mapping (ODM) is based on fluorescent labeling, stretching and imaging of single DNA molecules to obtain sequence-specific fluorescence profiles, DNA barcodes. These barcodes can be mapped to theoretical counterparts obtained from DNA reference sequences, which in turn allow for DNA identification in complex samples and for detecting structural changes in individual DNA molecules. There are several types of DNA labeling schemes for ODM and for each labeling type one or several types of match scoring methods are used. By combining the information from multiple labeling schemes one can potentially improve mapping confidence; however, combining match scores from different labeling assays has not been implemented yet. In this study, we introduce two theoretical methods for dealing with analysis of DNA molecules with multiple label types. In our first method, we convert the alignment scores, given as output from the different assays, into p-values using carefully crafted null models. We then combine the p-values for different label types using standard methods to obtain a combined match score and an associated combined p-value. In the second method, we use a block bootstrap approach to check for the uniqueness of a match to a database for all barcodes matching with a combined p-value below a predefined threshold. For obtaining experimental dual-labeled DNA barcodes, we introduce a novel assay where we cut plasmid DNA molecules from bacteria with restriction enzymes and the cut sites serve as sequence-specific markers, which together with barcodes obtained using the established competitive binding labeling method, form a dual-labeled barcode. All experimental data in this study originates from this assay, but we point out that our theoretical framework can be used to combine data from all kinds of available optical DNA mapping assays. We test our multiple labeling frameworks on barcodes from two different plasmids and synthetically generated barcodes (combined competitive-binding- and nick-labeling). It is demonstrated that by simultaneously using the information from all label types, we can substantially increase the significance when we match experimental barcodes to a database consisting of theoretical barcodes for all sequenced plasmids.
Collapse
Affiliation(s)
- Erik Torstensson
- Department of Astronomy and Theoretical Physics, Lund University, Lund, Sweden
| | - Gaurav Goyal
- Department of Biology and Biological Engineering, Chalmers University of Technology, Gothenburg, Sweden
| | - Anna Johnning
- Department of Mathematical Sciences, Chalmers University of Technology and the University of Gothenburg, Gothenburg, Sweden
- Systems and Data Analysis, Fraunhofer-Chalmers Centre, Gothenburg, Sweden
- Centre for Antibiotic Resistance Research, CARe, University of Gothenburg, Gothenburg, Sweden
| | - Fredrik Westerlund
- Department of Biology and Biological Engineering, Chalmers University of Technology, Gothenburg, Sweden
| | - Tobias Ambjörnsson
- Department of Astronomy and Theoretical Physics, Lund University, Lund, Sweden
- * E-mail:
| |
Collapse
|
6
|
Weigel M, Barash L, Shchur L, Janke W. Understanding population annealing Monte Carlo simulations. Phys Rev E 2021; 103:053301. [PMID: 34134227 DOI: 10.1103/physreve.103.053301] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2021] [Accepted: 04/06/2021] [Indexed: 11/07/2022]
Abstract
Population annealing is a recent addition to the arsenal of the practitioner in computer simulations in statistical physics and it proves to deal well with systems with complex free-energy landscapes. Above all else, it promises to deliver unrivaled parallel scaling qualities, being suitable for parallel machines of the biggest caliber. Here we study population annealing using as the main example the two-dimensional Ising model, which allows for particularly clean comparisons due to the available exact results and the wealth of published simulational studies employing other approaches. We analyze in depth the accuracy and precision of the method, highlighting its relation to older techniques such as simulated annealing and thermodynamic integration. We introduce intrinsic approaches for the analysis of statistical and systematic errors and provide a detailed picture of the dependence of such errors on the simulation parameters. The results are benchmarked against canonical and parallel tempering simulations.
Collapse
Affiliation(s)
- Martin Weigel
- Centre for Fluid and Complex Systems, Coventry University, Coventry CV1 5FB, United Kingdom.,Institut für Physik, Technische Universität Chemnitz, 09107 Chemnitz, Germany
| | - Lev Barash
- Landau Institute for Theoretical Physics, 142432 Chernogolovka, Russia
| | - Lev Shchur
- Landau Institute for Theoretical Physics, 142432 Chernogolovka, Russia.,National Research University Higher School of Economics, 101000 Moscow, Russia
| | - Wolfhard Janke
- Institut für Theoretische Physik, Universität Leipzig, IPF 231101, 04081 Leipzig, Germany
| |
Collapse
|
7
|
Salas J. Phase diagram for the bisected-hexagonal-lattice five-state Potts antiferromagnet. Phys Rev E 2020; 102:032124. [PMID: 33076030 DOI: 10.1103/physreve.102.032124] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2020] [Accepted: 08/28/2020] [Indexed: 11/07/2022]
Abstract
In this paper we study the phase diagram of the five-state Potts antiferromagnet on the bisected-hexagonal lattice. This question is important since Delfino and Tartaglia recently showed that a second-order transition in a five-state Potts antiferromagnet is allowed, and the bisected-hexagonal lattice had emerged as a candidate for such a transition on numerical grounds. By using high-precision Monte Carlo simulations and two complementary analysis methods, we conclude that there is a finite-temperature first-order transition point. This one separates a paramagnetic high-temperature phase, and a low-temperature phase where five phases coexist. This phase transition is very weak in the sense that its latent heat (per edge) is two orders of magnitude smaller than that of other well-known weak first-order phase transitions.
Collapse
Affiliation(s)
- Jesús Salas
- Departamento de Matemáticas, Escuela Politécnica Superior, Universidad Carlos III de Madrid, 28911 Leganés, Spain and Grupo de Teorías de Campos y Física Estadística, Associate Unit Instituto Gregorio Millán (UC3M)-Instituto de Estructura de la Materia (CSIC), Madrid, Spain
| |
Collapse
|
8
|
Xu J, Ferrenberg AM, Landau DP. High-resolution Monte Carlo study of the order-parameter distribution of the three-dimensional Ising model. Phys Rev E 2020; 101:023315. [PMID: 32168706 DOI: 10.1103/physreve.101.023315] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2019] [Accepted: 01/29/2020] [Indexed: 11/07/2022]
Abstract
We apply extensive Monte Carlo simulations to study the probability distribution P(m) of the order parameter m for the simple cubic Ising model with periodic boundary condition at the transition point. Sampling is performed with the Wolff cluster flipping algorithm, and histogram reweighting together with finite-size scaling analyses are then used to extract a precise functional form for the probability distribution of the magnetization, P(m), in the thermodynamic limit. This form should serve as a benchmark for other models in the three-dimensional Ising universality class.
Collapse
Affiliation(s)
- Jiahao Xu
- Center for Simulational Physics, University of Georgia, Athens, Georgia 30602, USA
| | - Alan M Ferrenberg
- Department of Computer Science and Software Engineering, Miami University of Ohio, Oxford, Ohio 45056, USA
| | - David P Landau
- Center for Simulational Physics, University of Georgia, Athens, Georgia 30602, USA
| |
Collapse
|
9
|
Xu J, M. Ferrenberg A, P. Landau D. 92 Years of the Ising Model: A High Resolution Monte Carlo Study. ACTA ACUST UNITED AC 2018. [DOI: 10.1088/1742-6596/1012/1/012002] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
10
|
Ferrenberg AM, Xu J, Landau DP. Pushing the limits of Monte Carlo simulations for the three-dimensional Ising model. Phys Rev E 2018; 97:043301. [PMID: 29758673 DOI: 10.1103/physreve.97.043301] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2017] [Indexed: 06/08/2023]
Abstract
While the three-dimensional Ising model has defied analytic solution, various numerical methods like Monte Carlo, Monte Carlo renormalization group, and series expansion have provided precise information about the phase transition. Using Monte Carlo simulation that employs the Wolff cluster flipping algorithm with both 32-bit and 53-bit random number generators and data analysis with histogram reweighting and quadruple precision arithmetic, we have investigated the critical behavior of the simple cubic Ising Model, with lattice sizes ranging from 16^{3} to 1024^{3}. By analyzing data with cross correlations between various thermodynamic quantities obtained from the same data pool, e.g., logarithmic derivatives of magnetization and derivatives of magnetization cumulants, we have obtained the critical inverse temperature K_{c}=0.221654626(5) and the critical exponent of the correlation length ν=0.629912(86) with precision that exceeds all previous Monte Carlo estimates.
Collapse
Affiliation(s)
- Alan M Ferrenberg
- Information Technology Services and Department of Chemical, Paper & Biomedical Engineering, Miami University, Oxford, Ohio 45056, USA
| | - Jiahao Xu
- Center for Simulational Physics, University of Georgia, Athens, Georgia 30602, USA
| | - David P Landau
- Center for Simulational Physics, University of Georgia, Athens, Georgia 30602, USA
| |
Collapse
|
11
|
Elçi EM, Weigel M. Efficient simulation of the random-cluster model. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:033303. [PMID: 24125381 DOI: 10.1103/physreve.88.033303] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/27/2013] [Indexed: 06/02/2023]
Abstract
The simulation of spin models close to critical points of continuous phase transitions is heavily impeded by the occurrence of critical slowing down. A number of cluster algorithms, usually based on the Fortuin-Kasteleyn representation of the Potts model, and suitable generalizations for continuous-spin models have been used to increase simulation efficiency. The first algorithm making use of this representation, suggested by Sweeny in 1983, has not found widespread adoption due to problems in its implementation. However, it has been recently shown that it is indeed more efficient in reducing critical slowing down than the more well-known algorithm due to Swendsen and Wang. Here, we present an efficient implementation of Sweeny's approach for the random-cluster model using recent algorithmic advances in dynamic connectivity algorithms.
Collapse
Affiliation(s)
- Eren Metin Elçi
- Applied Mathematics Research Centre, Coventry University, Coventry, CV1 5FB, United Kingdom and Institut für Physik, Johannes Gutenberg-Universität Mainz, Staudinger Weg 7, D-55099 Mainz, Germany
| | | |
Collapse
|
12
|
Romero-Vargas Castrillón S, Matysiak S, Stillinger FH, Rossky PJ, Debenedetti PG. Thermal stability of hydrophobic helical oligomers: a lattice simulation study in explicit water. J Phys Chem B 2012; 116:9963-70. [PMID: 22877080 DOI: 10.1021/jp305134w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We investigate the thermal stability of helical hydrophobic oligomers using a three-dimensional, water-explicit lattice model and the Wang-Landau Monte Carlo method. The degree of oligomer helicity is controlled by the parameter ε(mm) < 0, which mimics monomer-monomer hydrogen bond interactions leading to the formation of helical turns in atomistic proteins. We vary |ε(mm)| between 0 and 4.5 kcal/mol and therefore investigate systems ranging from flexible homopolymers (i.e., those with no secondary structure) to helical oligomers that are stable over a broad range of temperatures. We find that systems with |ε(mm)| ≤ 2.0 kcal/mol exhibit a broad thermal unfolding transition at high temperature, leading to an ensemble of random coils. In contrast, the structure of conformations involved in a second, low-temperature, transition is strongly dependent on |ε(mm)|. Weakly helical oligomers are observed when |ε(mm)| ≤ 1.0 kcal/mol and exhibit a low-temperature, cold-unfolding-like transition to an ensemble of strongly water-penetrated globular conformations. For higher |ε(mm)| (1.7 kcal/mol ≤ |ε(mm)| ≤ 2.0 kcal/mol), cold unfolding is suppressed, and the low-temperature conformational transition becomes a "crystallization", in which a "molten" helix is transformed into a defect-free helix. The molten helix preserves ≥50% of the helical contacts observed in the "crystal" at a lower temperature. When |ε(mm)| = 4.5 kcal/mol, we find that conformational transitions are largely suppressed within the range of temperatures investigated.
Collapse
|
13
|
Romero-Vargas Castrillón S, Matysiak S, Stillinger FH, Rossky PJ, Debenedetti PG. Phase Behavior of a Lattice Hydrophobic Oligomer in Explicit Water. J Phys Chem B 2012; 116:9540-8. [DOI: 10.1021/jp3039237] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
| | - Silvina Matysiak
- Fischell Department of Bioengineering, University of Maryland, College Park, Maryland 20742,
United States
| | - Frank H. Stillinger
- Department
of Chemistry, Princeton University, Princeton,
New Jersey 08544,
United States
| | - Peter J. Rossky
- Institute for Computational Engineering & Sciences and Department of Chemistry and Biochemistry, University of Texas at Austin, Austin, Texas 78712, United States
| | - Pablo G. Debenedetti
- Department of Chemical and Biological
Engineering, Princeton University, Princeton,
New Jersey 08544, United States
| |
Collapse
|
14
|
Parisen Toldin F. Improvement of Monte Carlo estimates with covariance-optimized finite-size scaling at fixed phenomenological coupling. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 84:025703. [PMID: 21929057 DOI: 10.1103/physreve.84.025703] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/13/2011] [Revised: 06/20/2011] [Indexed: 05/31/2023]
Abstract
In the finite-size scaling analysis of Monte Carlo data, instead of computing the observables at fixed Hamiltonian parameters, one may choose to keep a renormalization-group invariant quantity, also called phenomenological coupling, fixed at a given value. Within this scheme of finite-size scaling, we exploit the statistical covariance between the observables in a Monte Carlo simulation in order to reduce the statistical errors of the quantities involved in the computation of the critical exponents. This method is general and does not require additional computational time. This approach is demonstrated in the Ising model in two and three dimensions, where large gain factors in CPU time are obtained.
Collapse
Affiliation(s)
- Francesco Parisen Toldin
- Max-Planck-Institut für Physik Komplexer Systeme, Nöthnitzer Strasse 38, D-01187 Dresden, Germany.
| |
Collapse
|