Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSlow Mixing of Glauber Dynamics for the Six-Vertex Model in the Ordered Phases
142
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The six-vertex model in statistical physics is a weighted generalization of the ice model on $mathbb{Z}^2$ (i.e., Eulerian orientations) and the zero-temperature three-state Potts model (i.e., proper three-colorings). The phase diagram of the model depicts its physical properties and suggests where local Markov chains will be efficient. In this paper, we analyze the mixing time of Glauber dynamics for the six-vertex model in the ordered phases. Specifically, we show that for all Boltzmann weights in the ferroelectric phase, there exist boundary conditions such that local Markov chains require exponential time to converge to equilibrium. This is the first rigorous result bounding the mixing time of Glauber dynamics in the ferroelectric phase. Our analysis demonstrates a fundamental connection between correlated random walks and the dynamics of intersecting lattice path models (or routings). We analyze the Glauber dynamics for the six-vertex model with free boundary conditions in the antiferroelectric phase and significantly extend the region for which local Markov chains are known to be slow mixing. This result relies on a Peierls argument and novel properties of weighted non-backtracking walks.
rate research
Read More
We prove an optimal $Omega(n^{-1})$ lower bound on the spectral gap of Glauber dynamics for anti-ferromagnetic two-spin systems with $n$ vertices in the tree uniqueness regime. This spectral gap holds for all, including unbounded, maximum degree $Delta$. Consequently, we have the following mixing time bounds for the models satisfying the uniqueness condition with a slack $deltain(0,1)$: $bullet$ $C(delta) n^2log n$ mixing time for the hardcore model with fugacity $lambdale (1-delta)lambda_c(Delta)= (1-delta)frac{(Delta - 1)^{Delta - 1}}{(Delta - 2)^Delta}$; $bullet$ $C(delta) n^2$ mixing time for the Ising model with edge activity $betainleft[frac{Delta-2+delta}{Delta-delta},frac{Delta-delta}{Delta-2+delta}right]$; where the maximum degree $Delta$ may depend on the number of vertices $n$, and $C(delta)$ depends only on $delta$. Our proof is built upon the recently developed connections between the Glauber dynamics for spin systems and the high-dimensional expander walks. In particular, we prove a stronger notion of spectral independence, called the complete spectral independence, and use a novel Markov chain called the field dynamics to connect this stronger spectral independence to the rapid mixing of Glauber dynamics for all degrees.
We study numerically the two-point correlation functions of height functions in the six-vertex model with domain wall boundary conditions. The correlation functions and the height functions are computed by the Markov chain Monte-Carlo algorithm. Particular attention is paid to the free fermionic point ($Delta=0$), for which the correlation functions are obtained analytically in the thermodynamic limit. A good agreement of the exact and numerical results for the free fermionic point allows us to extend calculations to the disordered ($|Delta|<1$) phase and to monitor the logarithm-like behavior of correlation functions there. For the antiferroelectric ($Delta<-1$) phase, the exponential decrease of correlation functions is observed.
We show that the height function of the six-vertex model, in the parameter range $mathbf a=mathbf b=1$ and $mathbf cge1$, is delocalized with logarithmic variance when $mathbf cle 2$. This complements the earlier proven localization for $mathbf c>2$. Our proof relies on Russo--Seymour--Welsh type arguments, and on the local behaviour of the free energy of the cylindrical six-vertex model, as a function of the unbalance between the number of up and down arrows.
We develop an efficient method to compute the torus partition function of the six-vertex model exactly for finite lattice size. The method is based on the algebro-geometric approach to the resolution of Bethe ansatz equations initiated in a previous work, and on further ingredients introduced in the present paper. The latter include rational $Q$-system, primary decomposition, algebraic extension and Galois theory. Using this approach, we probe new structures in the solution space of the Bethe ansatz equations which enable us to boost the efficiency of the computation. As an application, we study the zeros of the partition function in a partial thermodynamic limit of $M times N$ tori with $N gg M$. We observe that for $N to infty$ the zeros accumulate on some curves and give a numerical method to generate the curves of accumulation points.
We consider a system of classical particles, interacting via a smooth, long-range potential, in the mean-field regime, and we optimally analyze the propagation of chaos in form of sharp estimates on many-particle correlation functions. While approaches based on the BBGKY hierarchy are doomed by uncontrolled losses of derivatives, we propose a novel non-hierarchical approach that focusses on the empirical measure of the system and exploits a Glauber type calculus with respect to initial data in form of higher-order Poincare inequalities for cumulants. This main result allows to rigorously truncate the BBGKY hierarchy to an arbitrary precision on the mean-field timescale, thus justifying the Bogolyubov corrections to mean field. As corollaries, we also deduce a quantitative central limit theorem for fluctuations of the empirical measure, and we partially justify the Lenard-Balescu limit for a spatially homogeneous system away from thermal equilibrium.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
