BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Wenan Guo

Wenan Guo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.str-elcond-mat.quant-gascond-mat.softcond-matcond-mat.dis-nncond-mat.othercond-mat.supr-conhep-lat

Catalog footprint

What is connected

43works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

43 published item(s)

preprint2022arXiv

Exotic surface behaviors induced by geometrical settings of the two-dimensional dimerized quantum XXZ model

We study the surface behavior of the two-dimensional columnar dimerized quantum antiferromagnetic XXZ model with easy-plane anisotropy, with particular emphasis on the surface critical behaviors of the (2+1)-dimensional quantum critical points of the model that belong to the classical three-dimensional O(2) universality class, for both $S=1/2$ and $S=1$ spins using quantum Monte Carlo simulations. We find completely different surface behaviors on two different surfaces of geometrical settings: the dangling-ladder surface, which is exposed by cutting a row of weak bonds, and the dangling-chain surface, which is formed by cutting a row of strong bonds along the direction perpendicular to the strong bonds of a periodic system. Similar to the Heisenberg limit, we find an ordinary transition on the dangling-ladder surface for both $S=1$ and $S=1/2$ spin systems. However, the dangling-chain surface shows much richer surface behaviors than in the Heisenberg limit. For the $S=1/2$ easy-plane model, at the bulk critical point, we provide evidence supporting an extraordinary surface transition with a long-range order established by effective long-range interactions due to bulk critical fluctuations. The possibility that the state is an extraordinary-log state seems unlikely. For the $S=1$ system, we find surface behaviors similar to that of the three-dimensional classical XY model with sufficiently enhanced surface coupling, suggesting an extraordinary-log state at the bulk critical point.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Surface critical properties of the three-dimensional clock model

Using Monte Carlo simulations and finite-size scaling analysis, we show that the $q$-state clock model with $q=6$ on the simple cubic lattice with open surfaces has a rich phase diagram; in particular, it has an extraordinary-log phase, besides the ordinary and extraordinary transitions at the bulk critical point. We prove numerically that the presence of the intermediate extraordinary-log phase is due to the emergence of an O(2) symmetry in the surface state before the surface enters the $Z_{q}$ symmetry-breaking region as the surface coupling is increased at the bulk critical point, while O(2) symmetry emerges for the bulk. The critical behaviors of the extraordinary-log transition, as well as the ordinary and the special transition separating the ordinary and the extraordinary-log transition are obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Thermodynamic and Dynamical Signatures of a Quantum Spin-Hall Insulator to Superconductor Transition

Thermodynamic and dynamical properties of a model of Dirac fermions with a deconfined quantum critical point (DQCP) separating an interaction-generated quantum spin-Hall insulator from an s-wave superconductor [Nature Comm.~{\bf 10}, 2658 (2019)] are studied by quantum Monte Carlo simulations. Inside the deconfined quantum critical region bound by the single-particle gap, spinons and spinless charge-2e skyrmions emerge. Since the model conserves total spin and charge, and has a single length scale, these excitations lead to a characteristic linear temperature dependence of the uniform spin and charge susceptibilities. At the DQCP, the order parameter dynamic structure factors show remarkable similarities that support emergent Lorentz symmetry. Above a critical temperature, superconductivity is destroyed by the proliferation of spin-1/2 vortices.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Extreme Suppression of Antiferromagnetic Order and Critical Scaling in a Two-Dimensional Random Quantum Magnet

Sr$_2$CuTeO$_6$ is a square-lattice Néel antiferromagnet with superexchange between first-neighbor $S=1/2$ Cu spins mediated by plaquette centered Te ions. Substituting Te by W, the affected impurity plaquettes have predominantly second-neighbor interactions, thus causing local magnetic frustration. Here we report a study of Sr$_2$CuTe$_{1-x}$W$_x$O$_6$ using neutron diffraction and $μ$SR techniques, showing that the Néel order vanishes already at $x = 0.025 \pm 0.005$. We explain this extreme order suppression using a two-dimensional Heisenberg spin model, demonstrating that a W-type impurity induces a deformation of the order parameter that decays with distance as $1/r^2$ at temperature $T=0$. The associated logarithmic singularity leads to loss of order for any $x>0$. Order for small $x>0$ and $T>0$ is induced by weak interplane couplings. In the nonmagnetic phase of Sr$_2$CuTe$_{1-x}$W$_x$O$_6$, the $μ$SR relaxation rate exhibits quantum critical scaling with a large dynamic exponent, $z \approx 3$, consistent with a random-singlet state.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Monte Carlo Renormalization Flows in the Space of Relevant and Irrelevant Operators: Application to Three-Dimensional Clock Models

We present a way to visualize and quantify renormalization group flows in a space of observables computed using Monte Carlo simulations. We apply the method to classical three-dimensional clock models, i.e., the planar (XY) spin model perturbed by a $Z_q$ symmetric anisotropy field. The method performs significantly better than standard techniques for determining the scaling dimension $y_q$ of the $Z_q$ field at the critical point if it is irrelevant ($q\ge4$). Furthermore, we analyze all stages of the complex renormalization flow, including the cross-over from the U(1) Nambu-Goldstone fixed point to the ultimate $Z_q$ symmetry-breaking fixed point due to the relevance of the $Z_q$ field inside the ordered phase. We expect our method to be particularly useful in the context of quantum-critical points with inherent dangerously irrelevant operators that cannot be tuned away microscopically but whose renormalization flows can be analyzed exactly as we do here for the clock models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantum-critical scaling properties of the two-dimensional random-singlet state

We use QMC simulations to study effects of disorder on the $S=1/2$ Heisenberg model with exchange constant $J$ on the square lattice supplemented by multispin interactions $Q$. It was found recently [L. Lu et al., Phys. Rev. X 8, 041040 (2018)] that the ground state of this $J$-$Q$ model with random couplings undergoes a quantum phase transition from the Néel state into a randomness-induced spin-liquid-like state that is a close analogue to the well known random-singlet (RS) state of the random Heisenberg chain. The 2D RS state arises from spinons localized at topological defects. The interacting spinons form a critical state with mean spin-spin correlations decaying with distance $r$ as $r^{-2}$, as in the 1D RS state. The dynamic exponent $z \ge 2$, varying continuously with the model parameters. We here further investigate the properties of the RS state in the $J$-$Q$ model with random $Q$ couplings. We study the temperature dependence of the specific heat and various susceptibilities for large enough systems to reach the thermodynamic limit and also analyze the size dependence of the critical magnetic order parameter and its susceptibility in the ground state. For all these quantities, we find consistency with conventional quantum-critical scaling when the condition implied by the $r^{-2}$ form of the spin correlations is imposed. All quantities can be explained by the same value of the dynamic exponent $z$ at fixed model parameters. We argue that the RS state identified in the $J$-$Q$ model corresponds to a generic renormalization group fixed point that can be reached in many quantum magnets with random couplings, and may already have been observed experimentally.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Surface critical behaviors of coupled Haldane chains

The special surface transition at (2+1)-dimensional quantum critical point is precluded in corresponding classical critical point. The mechanism of such behavior which is only found in dimerized Heisenberg models so far is still under debate. To illuminate the role of symmetry protected topological (SPT) phase in inducing such nonordinary behaviors, we study a system on a two-dimensional square lattice consisted by interacting spin-1 Haldane chains, which has a genuine SPT phase--the Haldane phase--at weak interchain interactions and a quantum critical point belonging to the classical 3D O(3) universality class to the Néel phase. Different from models studied previously, there is no dimerization in the current model. Cutting the system along the chain direction or perpendicular to the chain direction exposes two different surfaces. Using unbiased quantum Monte Carlo simulations, we find that the two different types of surface show completely different surface critical behaviors at the bulk critical point, resulted from different surface states in the SPT phase. For the system with surfaces along the chain direction, the surface critical behavior is of ordinary type of the bulk 3D O(3) critical point, while for the surfaces perpendicular to the chain direction, the surface critical behavior is nonordinary, consistent with special transitions found in dimerized Heisenberg models. Our numerical results demonstrate that the gapless surface state in the gapped SPT phase together with the gapless mode of critical point is a pure quantum scenario that leads to the nonordinary transition.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The three-state Potts model on the centered triangular lattice

We study phase transitions of the Potts model on the centered-triangular lattice with two types of couplings, namely $K$ between neighboring triangular sites, and $J$ between the centered and the triangular sites. Results are obtained by means of a finite-size analysis based on numerical transfer matrix calculations and Monte Carlo simulations. Our investigation covers the whole $(K, J)$ phase diagram, but we find that most of the interesting physics applies to the antiferromagnetic case $K<0$, where the model is geometrically frustrated. In particular, we find that there are, for all finite $J$, two transitions when K is varied. Their critical properties are explored. In the limits $J\to \pm \infty$ we find algebraic phases with infinite-order transitions to the ferromagnetic phase.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Equivalent-neighbor Potts models in two dimensions

We investigate the two-dimensional $q=3$ and 4 Potts models with a variable interaction range by means of Monte Carlo simulations. We locate the phase transitions for several interaction ranges as expressed by the number $z$ of equivalent neighbors. For not too large $z$, the transitions fit well in the universality classes of the short-range Potts models. However, at longer ranges the transitions become discontinuous. For $q=3$ we locate a tricritical point separating the continuous and discontinuous transitions near $z=80$, and a critical fixed point between $z=8$ and 12. For $q=4$ the transition becomes discontinuous for $z > 16$. The scaling behavior of the $q=4$ model with $z=16$ approximates that of the $q=4$ merged critical-tricritical fixed point predicted by the renormalization scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quantum criticality with two length scales

The theory of deconfined quantum critical points describes phase transitions at temperature T = 0 outside the standard paradigm, predicting continuous transformations between certain ordered states where conventional theory requires discontinuities. Numerous computer simulations have offered no proof of such transitions, however, instead finding deviations from expected scaling relations that were neither predicted by the DQC theory nor conform to standard scenarios. Here we show that this enigma can be resolved by introducing a critical scaling form with two divergent length scales. Simulations of a quantum magnet with antiferromagnetic and dimerized ground states confirm the form, proving a continuous transition with deconfined excitations and also explaining anomalous scaling at T > 0. Our findings revise prevailing paradigms for quantum criticality, with potentially far-reaching implications for many strongly-correlated materials.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Special transitions in an O($n$) loop model with an Ising-like constraint

We investigate the O($n$) nonintersecting loop model on the square lattice under the constraint that the loops consist of ninety-degree bends only. The model is governed by the loop weight $n$, a weight $x$ for each vertex of the lattice visited once by a loop, and a weight $z$ for each vertex visited twice by a loop. We explore the $(x,z)$ phase diagram for some values of $n$. For $0<n<1$, the diagram has the same topology as the generic O($n$) phase diagram with $n<2$, with a first-order line when $z$ starts to dominate, and an O($n$)-like transition when $x$ starts to dominate. Both lines meet in an exactly solved higher critical point. For $n>1$, the O($n$)-like transition line appears to be absent. Thus, for $z=0$, the $(n,x)$ phase diagram displays a line of phase transitions for $n\le 1$. The line ends at $n=1$ in an infinite-order transition. We determine the conformal anomaly and the critical exponents along this line. These results agree accurately with a recent proposal for the universal classification of this type of model, at least in most of the range $-1 \leq n \leq 1$. We also determine the exponent describing crossover to the generic O($n$) universality class, by introducing topological defects associated with the introduction of `straight' vertices violating the ninety-degree-bend rule. These results are obtained by means of transfer-matrix calculations and finite-size scaling.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Anomalous quantum glass of bosons in a random potential in two dimensions

We present a quantum Monte Carlo study of the "quantum glass" phase of the 2D Bose-Hubbard model with random potentials at filling $ρ=1$. In the narrow region between the Mott and superfluid phases the compressibility has the form $κ\sim {\rm exp}(-b/T^α)+c$ with $α<1$ and $c$ vanishing or very small. Thus, at $T=0$ the system is either incompressible (a Mott glass) or nearly incompressible (a Mott-glass-like anomalous Bose glass). At stronger disorder, where a glass reappears from the superfluid, we find a conventional highly compressible Bose glass. On a path connecting these states, away from the superfluid at larger Hubbard repulsion, a change of the disorder strength by only $10\%$ changes the low-temperature compressibility by more than four orders of magnitude, lending support to two types of glass states separated by a phase transition or a sharp cross-over.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Emergent topological excitations in a two-dimensional quantum spin system

We study the mechanism of decay of a topological (winding-number) excitation due to finite-size effects in a two-dimensional valence-bond solid state, realized in an $S=1/2$ spin model ($J$-$Q$ model) and studied using projector Monte Carlo simulations in the valence bond basis. A topological excitation with winding number $|W|>0$ contains domain walls, which are unstable due to the emergence of long valence bonds in the wave function, unlike in effective descriptions with the quantum dimer model. We find that the life time of the winding number in imaginary time diverges as a power of the system length $L$. The energy can be computed within this time (i.e., it converges toward a "quasi-eigenvalue" before the winding number decays) and agrees for large $L$ with the domain-wall energy computed in an open lattice with boundary modifications enforcing a domain wall. Constructing a simplified two-state model and using the imaginary-time behavior from the simulations as input, we find that the real-time decay rate out of the initial winding sector is exponentially small in $L$. Thus, the winding number rapidly becomes a well-defined conserved quantum number for large systems, supporting the conclusions reached by computing the energy quasi-eigenvalues. Including Heisenberg exchange interactions which brings the system to a quantum-critical point separating the valence-bond solid from an antiferromagnetic ground state (the putative "deconfined" quantum-critical point), we can also converge the domain wall energy here and find that it decays as a power-law of the system size. Thus, the winding number is an emergent quantum number also at the critical point, with all winding number sectors becoming degenerate in the thermodynamic limit. This supports the description of the critical point in terms of a U(1) gauge-field theory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Topological properties of a Valence-Bond-Solid

We present a projector quantum Monte Carlo study of the topological properties of the valence-bond-solid ground state in the $J$-$Q_3$ spin model on the square lattice. The winding number is a topological number counting the number of domain walls in the system and is a good quantum number in the thermodynamic limit. We study the finite-size behaviour and obtain the domain wall energy density for a topological nontrivial valence-bond-solid state.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Completely packed O($n$) loop models and their relation with exactly solved coloring models

We explore the physical properties of the completely packed O($n$) loop model on the square lattice, and its generalization to an Eulerian graph model, which follows by including cubic vertices which connect the four incoming loop segments. This model includes crossing bonds as well. Our study of the properties of this model involve transfer-matrix calculations and finite-size scaling. The numerical results are compared to existing exact solutions, including solutions of special cases of a so-called coloring model, which are shown to be equivalent with our generalized loop model. The latter exact solutions correspond with seven one-dimensional branches in the parameter space of our generalized loop model. One of these branches, describing the case of nonintersecting loops, is already known to correspond with the ordering transition of the Potts model. We find that another exactly solved branch, which describes a model with nonintersecting loops and cubic vertices, corresponds with a first-order Ising-like phase transition for $n>2$. For $1<n<2$, this branch can be interpreted in terms of a low-temperature O($n$) phase with corner-cubic anisotropy. For $n>2$ this branch is the locus of a first-order phase boundary between a phase with a hard-square lattice-gas like ordering, and a phase dominated by cubic vertices. The first-order character of this transition is in agreement with a mean-field argument.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Exact finite-size corrections and corner free energies for the c=-2 universality class

We consider (a) the partition functions of the anisotropic dimer model on the rectangular (2M-1) x (2N-1) lattice with free and cylindrical boundary conditions with a single monomer residing on the boundary and (b) the partition function of the anisotropic spanning tree on an M x N rectangular lattice with free boundary conditions. We express (a) and (b) in terms of a principal partition function with twisted boundary conditions. Based on these expressions, we derive the exact asymptotic expansions of the free energy for both cases (a) and (b). We confirm the conformal field theory prediction for the corner free energy of these models, and find the central charge is c = - 2. We also show that the dimer model on the cylinder with an odd number of sites on the perimeter exhibits the same finite-size corrections as on the plane.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Reentrance of Berezinskii-Kosterlitz-Thouless-like transitions in three-state Potts antiferromagnetic thin film

Using Monte Carlo simulations and finite-size scaling, we study three-state Potts antiferromagnet on layered square lattice with two and four layers $L_z=2$ and $4$. As temperature decreases, the system develops quasi-long-range order via a Berezinskii-Kosterlitz-Thouless transition at finite temperature $T_{c1}$. For $L_z=4$, as temperature is further lowered, a long-range order breaking the $Z_6$ symmetry develops at a second transition at $T_{c2} < T_{c1}$. The transition at $T_{c2}$ is also Berezinskii-Kosterlitz-Thouless-like, but has magnetic critical exponent $η=1/9$ instead of the conventional value $η= 1/4$. The emergent $U(1)$ symmetry is clearly demonstrated in the quasi-long-range ordered region $T_{c2} \leq T \leq T_{c1}$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rényi Information flow in the Ising model with single-spin dynamics

The $n$-index Rényi mutual information and transfer entropy for the two-dimensional kinetic Ising model with arbitrary single-spin dynamics in the thermodynamic limit are derived as functions of thermodynamic quantities. By means of Monte Carlo simulations with the Wolff algorithm, we calculate the information flows in the Ising model with the Metropolis dynamics and the Glauber dynamics. We find that, not only the global Rényi transfer entropy, but also the pairwise Rényi transfer entropy peaks in the disorder phase. Therefore, the Rényi information flows may be used as better tools than the Shannon counterparts in the study of phase transitions in complex dynamical systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Accurate expansions of internal energy and specific heat of critical two-dimensional Ising model with free boundaries

The bond-propagation (BP) algorithm for the specific heat of the two dimensional Ising model is developed and that for the internal energy is completed. Using these algorithms, we study the critical internal energy and specific heat of the model on the square lattice and triangular lattice with free boundaries. Comparing with previous works [X.-T. Wu {\it et al} Phys. Rev. E {\bf 86}, 041149 (2012) and Phys. Rev. E {\bf 87}, 022124 (2013)], we reach much higher accuracy ($10^{-26}$) of the internal energy and specific heat, compared to the accuracy $10^{-11}$ of the internal energy and $10^{-9}$ of the specific heat reached in the previous works. This leads to much more accurate estimations of the edge and corner terms. The exact values of some edge and corner terms are therefore conjectured. The accurate forms of finite-size scaling for the internal energy and specific heat are determined for the rectangle-shaped square lattice with various aspect ratios and various shaped triangular lattice. For the rectangle-shaped square and triangular lattices and the triangle-shaped triangular lattice, there is no logarithmic correction terms of order higher than 1/S, with S the area of the system. For the triangular lattice in rhombus, trapezoid and hexagonal shapes, there exist logarithmic correction terms of order higher than 1/S for the internal energy, and logarithmic correction terms of all orders for the specific heat.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Critical properties of the Hintermann-Merlini model

Many critical properties of the Hintermann-Merlini model are known exactly through the mapping to the eight-vertex model. Wu [J. Phys. C {\bf 8}, 2262 (1975)] calculated the spontaneous magnetizations of the model on two sublattices by relating them to the conjectured spontaneous magnetization and polarization of the eight-vertex model, respectively. The latter conjecture remains unproved. In this paper, we numerically study the critical properties of the model by means of a finite-size scaling analysis based on transfer matrix calculations and Monte Carlo simulations. All analytic predictions for the model are confirmed by our numerical results. The central charge $c=1$ is found for the critical manifold investigated. In addition, some unpredicted geometry properties of the model are studied. Fractal dimensions of the largest Ising clusters on two sublattices are determined. The fractal dimension of the largest Ising cluster on the sublattice A takes fixed value $D_{\rm a}=1.888(2)$, while that for sublattice B varies continuously with the parameters of the model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Extended Bose-Hubbard model with pair hopping on the triangular lattice

We study systematically an extended Bose-Hubbard model on the triangular lattice by means of a meanfield method based on the Gutzwiller ansatz. Pair hopping terms are explicitly included and a three-body constraint is applied. The ground-state phase diagram and a variety of quantum phase transitions are investigated in great detail. In particular, we show the existence and stability of the pair supersolid phase.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Ising-like phase transition of an n-component Eulerian face-cubic model

By means of Monte Carlo simulations and a finite-size scaling analysis, we find a critical line of an n-component Eulerian face-cubic model on the square lattice and the simple cubic lattice in the region v>1, where v is the bond weight. The phase transition belongs to the Ising universality class independent of n. The critical properties of the phase transition can also be captured by the percolation of the complement of the Eulerian graph.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Ising-like transitions in the O($n$) loop model on the square lattice

We explore the phase diagram of the O($n$) loop model on the square lattice in the $(x,n)$ plane, where $x$ is the weight of a lattice edge covered by a loop. These results are based on transfer-matrix calculations and finite-size scaling. We express the correlation length associated with the staggered loop density in the transfer-matrix eigenvalues. The finite-size data for this correlation length, combined with the scaling formula, reveal the location of critical lines in the diagram. For $n>>2$ we find Ising-like phase transitions associated with the onset of a checkerboard-like ordering of the elementary loops, i.e., the smallest possible loops, with the size of an elementary face, which cover precisely one half of the faces of the square lattice at the maximum loop density. In this respect, the ordered state resembles that of the hard-square lattice gas with nearest-neighbor exclusion, and the finiteness of $n$ represents a softening of its particle-particle potentials. We also determine critical points in the range $-2\leq n\leq 2$. It is found that the topology of the phase diagram depends on the set of allowed vertices of the loop model. Depending on the choice of this set, the $n>2$ transition may continue into the dense phase of the $n \leq 2$ loop model, or continue as a line of $n \leq 2$ O($n$) multicritical points.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Phase transitions of the q-state Potts model on multiply-laced Sierpinski gaskets

We present an exact solution of the q-state Potts model on a class of generalized Sierpinski fractal lattices. The model is shown to possess an ordered phase at low temperatures and a continuous transition to the high temperature disordered phase at any q>=1. Multicriticality is observed in the presence of a symmetry-breaking field. Exact renormalization group analysis yields the phase diagram of the model and a complete set of critical exponents at various transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Skyrmion crystals in the pseudo-spin-1/2 Bose-Einstein condensates

Exact two-dimensional solutions are constructed for the pseudo-spin-1/2 Bose-Einstein condensates which are described by the coupled nonlinear Gross-Pitaevskii equations where the intraspecies and interspecies coupling constants are assumed to be equal. The equations are decoupled by means of re-combinations of the nonlinear terms of the hyperfine states according to the spatial dimensions. These stationary solutions form various spin textures which are identified as skyrmion crystals. In a special case, the crystal of skyrmion-antiskyrmion pairs is formed in the soliton limit.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Critical manifold of the Potts model: Exact results and homogeneity approximation

The $q$-state Potts model has stood at the frontier of research in statistical mechanics for many years. In the absence of a closed-form solution, much of the past efforts have focused on locating its critical manifold, trajectory in the parameter $\{q, e^J\}$ space where $J$ is the reduced interaction, along which the free energy is singular. However, except in isolated cases, antiferromagnetic (AF) models with $J<0$ have been largely neglected. In this paper we consider the Potts model with AF interactions focusing on deducing its critical manifold in exact and/or closed-form expressions. We first re-examine the known critical frontiers in light of AF interactions. For the square lattice we confirm the Potts self-dual point to be the sole critical point for $J>0$. We also locate its critical frontier for $J<0$ and find it to coincide with a solvability condition observed by Baxter in 1982. For the honeycomb lattice we show that the known critical point holds for {all} $J$, and determine its critical $q_c = \frac 1 2 (3+\sqrt 5) = 2.61803$ beyond which there is no transition. For the triangular lattice we confirm the known critical point to hold only for $J>0$. More generally we consider the centered-triangle (CT) and Union-Jack (UJ) lattices consisting of mixed $J$ and $K$ interactions, and deduce critical manifolds under homogeneity hypotheses. For K=0 the CT lattice is the diced lattice, and we determine its critical manifold for all $J$ and find $q_c = 3.32472$. For K=0 the UJ lattice is the square lattice and from this we deduce both the $J>0$ and $J<0$ critical manifolds and find $q_c=3$ for the square lattice. Our theoretical predictions are compared with recent tensor-based numerical results and Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Exact critical points of the O($n$) loop model on the martini and the 3-12 lattices

We derive the exact critical line of the O($n$) loop model on the martini lattice as a function of the loop weight $n$.A finite-size scaling analysis based on transfer matrix calculations is also performed.The numerical results coincide with the theoretical predictions with an accuracy up to 9 decimal places. In the limit $n\to 0$, this gives the exact connective constant $μ=1.7505645579...$ of self-avoiding walks on the martini lattice. Using similar numerical methods, we also study the O($n$) loop model on the 3-12 lattice. We obtain similarly precise agreement with the exact critical points given by Batchelor [J. Stat. Phys. 92, 1203 (1998)].

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Finite size behaviors of critical Ising model on a rectangle with free boundaries

Using the bond-propagation algorithm, we study the Ising model on a rectangle of size $M \times N$ with free boundaries. For five aspect ratios $ρ=M/N=1,2,4,8,16$, the critical free energy, internal energy and specific heat are calculated. The largest size reached is $M \times N=64\times 10^6$. The accuracy of the free energy reaches $10^{-26}$. Basing on these accurate data, we determine exact expansions of the critical free energy, internal energy and specific heat. With these expansions, we extract the bulk, surface and corner parts of free energy, internal energy and specific heat. The fitted bulk free energy density is given by $f_{\infty}=0.92969539834161021499(1)$, comparing with Onsager's exact result $f_{\infty}=0.92969539834161021506...$. We prove the conformal field theory(CFT) prediction of the corner free energy, in which the central charge of the Ising model is found to be $c=0.5\pm 1\times 10^{-10}$ comparing with the CFT result $c=0.5$. We find that not only the corner free energy but also the corner internal energy and specific heat are geometry independent, i.e., independent of aspect ratio. The implication of this finding on the finite scaling is discussed. In the second order correction of the free energy, we prove the geometry dependence predicted by CFT and find out a geometry independent constant beyond CFT. High order corrections are also obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Pair supersolid of the extended Bose-Hubbard model with atom-pair hopping on the triangular Lattice

We systematically study an extended Bose-Hubbard model with atom hopping and atom-pair hopping in the presence of a three-body constraint on the triangular lattice. By means of large-scale Quantum Monte Carlo simulations, the ground-state phase diagram are studied. We find a continuous transition between the atomic superfluid phase and the pair superfluid when the ratio of the atomic hopping and the atom-pair hopping is adapted. We then focus on the interplay among the atom-pair hopping, the on-site repulsion and the nearest-neighbor repulsion. With on-site repulsion present, we observe first order transitions between the Mott Insulators and pair superfluid driven by the pair hopping. With the nearest-neighbor repulsion turning on, three typical solid phases with 2/3, 1 and 4/3-filling emerge at small atom-pair hopping region. A stable pair supersolid phase is found at small on-site repulsion. This is due to the three-body constraint and the pair hopping, which essentially make the model a quasi hardcore boson system. Thus the pair supersolid state emerges basing on the order-by-disorder mechanism, by which hardcore bosons avoid classical frustration on the triangular lattice. The transition between the pair supersolid and the pair superfluid is first order, except for the particle-hole symmetric point. We compare the results with those obtained by means of mean-field analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Phase transitions in the frustrated Ising model on the square lattice

We consider the thermal phase transition from a paramagnetic to stripe-antiferromagnetic phase in the frustrated two-dimensional square-lattice Ising model with competing interactions J1<0 (nearest neighbor, ferromagnetic) and J2 >0 (second neighbor, antiferromagnetic). The striped phase breaks a Z4 symmetry and is stabilized at low temperatures for g=J2/|J1|>1/2. Despite the simplicity of the model, it has proved difficult to precisely determine the order and the universality class of the phase transitions. This was done convincingly only recently by Jin et al. [PRL 108, 045702 (2012)]. Here, we further elucidate the nature of these transitions and their anomalies by employing a combination of cluster mean-field theory, Monte Carlo simulations, and transfer-matrix calculations. The J1-J2 model has a line of very weak first-order phase transitions in the whole region 1/2<g<g*, where g* = 0.67(1). Thereafter, the transitions from g above g* are continuous and can be fully mapped, using universality arguments, to the critical line of the well known Ashkin-Teller model from its 4-state Potts point to the decoupled Ising limit. We also comment on the pseudo-first-order behavior at the Potts point and its neighborhood in the Ashkin-Teller model on finite lattices, which in turn leads to the appearance of similar effects in the vicinity of the multicritical point g* in the J1-J2 model. The continuous transitions near g* can therefore be mistaken to be first-order transitions, and this realization was the key to understanding the paramagnetic-striped transition for the full range of g>1/2. Most of our results are based on Monte Carlo calculations, while the cluster mean-field and transfer-matrix results provide useful methodological bench-marks for weakly first-order behaviors and Ashkin-Teller criticality.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Shape dependent finite-size effect of critical two-dimensional Ising model on a triangular lattice

Using the bond-propagation algorithm, we study the finite-size behavior of the critical two-dimensional Ising model on a finite triangular lattice with free boundaries in five shapes: triangle, rhombus, trapezoid, hexagon and rectangle. The critical free energy, internal energy and specific heat are calculated. The accuracy of the free energy reaches $10^{-26}$. Based on accurate data on several finite systems with linear size up to N=2000, we extract the bulk, surface and corner parts of the free energy, internal energy and specific heat accurately. We confirm the conformal field theory prediction of the corner free energy to be universal and find logarithmic corrections in higher order terms in the critical free energy for the rhombus, trapezoid, and hexagon shaped systems, which are absent for the triangle and rectangle shaped systems. The logarithmic edge corrections due to edges parallel or perpendicular to the bond directions in the internal energy are found to be identical, while the logarithmic edge corrections due to corresponding edges in the free energy and the specific heat are different. The corner internal energy and corner specific heat for angles $π/3$, $π/2$ and $2π/3$ are obtained, as well as higher order corrections. Comparing with the corner internal energy and corner specific heat previously found on a rectangle of the square lattice (Phys. Rev. E. 86 041149 (2012)), we conclude that the corner internal energy and corner specific heat for the rectangle shape are not universal.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Critical frontier for the Potts and percolation models on triangular-type and kagome-type lattices II: Numerical analysis

In a recent paper (arXiv:0911.2514), one of us (FYW) considered the Potts model and bond and site percolation on two general classes of two-dimensional lattices, the triangular-type and kagome-type lattices, and obtained closed-form expressions for the critical frontier with applications to various lattice models. For the triangular-type lattices Wu's result is exact, and for the kagome-type lattices Wu's expression is under a homogeneity assumption. The purpose of the present paper is two-fold: First, an essential step in Wu's analysis is the derivation of lattice-dependent constants $A, B, C$ for various lattice models, a process which can be tedious. We present here a derivation of these constants for subnet networks using a computer algorithm. Secondly, by means of a finite-size scaling analysis based on numerical transfer matrix calculations, we deduce critical properties and critical thresholds of various models and assess the accuracy of the homogeneity assumption. Specifically, we analyze the $q$-state Potts model and the bond percolation on the 3-12 and kagome-type subnet lattices $(n\times n):(n\times n)$, $n\leq 4$, for which the exact solution is not known. To calibrate the accuracy of the finite-size procedure, we apply the same numerical analysis to models for which the exact critical frontiers are known. The comparison of numerical and exact results shows that our numerical determination of critical thresholds is accurate to 7 or 8 significant digits. This in turn infers that the homogeneity assumption determines critical frontiers with an accuracy of 5 decimal places or higher. Finally, we also obtained the exact percolation thresholds for site percolation on kagome-type subnet lattices $(1\times 1):(n\times n)$ for $1\leq n \leq 6$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Crossover phenomena involving the dense O($n$) phase

We explore the properties of the low-temperature phase of the O($n$) loop model in two dimensions by means of transfer-matrix calculations and finite-size scaling. We determine the stability of this phase with respect to several kinds of perturbations, including cubic anisotropy, attraction between loop segments, double bonds and crossing bonds. In line with Coulomb gas predictions, cubic anisotropy and crossing bonds are found to be relevant and introduce crossover to different types of behavior. Whereas perturbations in the form of loop-loop attractions and double bonds are irrelevant, sufficiently strong perturbations of these types induce a phase transition of the Ising type, at least in the cases investigated. This Ising transition leaves the underlying universal low-temperature O($n$) behavior unaffected.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Effects of edge disorder in nano-scale antiferromagnetic clusters

We study the distribution of local magnetic susceptibilities in the two-dimensional antiferromagnetic S=1/2 Heisenberg model on various random clusters, in order to determine whether effects of edge disorder could be detected in NMR experiments (through the line shape, as given by the distribution of local Knight shifts). Although the effects depend strongly on the nature of the edge and the cluster size, our results indicate that line widths broader than the average shift should be expected even in clusters as large as $\approx 1000$ lattice spacing in diameter. Experimental investigations of the NMR line width should give insights into the magnetic structure of the edges.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Hardcore bosons on the dual of the bowtie lattice

We study the zero temperature phase diagram of hardcore bosons on the dual of the bowtie lattice. Two types of striped diagonal long-range order (striped order I and striped order II) are discussed. A state with type-II striped order and superfluidity is found, even without nearest-neighbor repulsion. The emergence of such a state is due to the inhomogeneity and the anisotropy of the lattice structure. However, neither the translational symmetry nor the symmetry between sublattices of the original lattice is broken. In this paper, we restrict a 'solid state' of lattice bosons as a diagonal long-range ordered state breaking either the translational symmetry of the original lattice or the symmetry of different sublattices. We thus name such a phase a striped superfluid phase (SSF). In the presence of a nearest-neighbor repulsion, we find two striped charge density wave phases(SCDW I and II) with boson density $ρ=1/2$ (with striped order I) and $ρ=2/3$ (with striped order II) respectively, when the hopping amplitude is small compared with the repulsion. The SCDW I state is a solid, in which the translational symmetry of the original lattice is broken. We observe a rather special first-order phase transition showing an interesting multi-loop hysteresis phenomenon between the two SCDW phases when the hopping term is small enough. This can be accounted for by the special degeneracy of the ground states near the classical limit. The SSF re-appears outside the two SCDW phases. The transition between the SCDW I and SSF phases is first order, while the transition between SCDW II and SSF phases is continuous. We find that the superfluid stiffness is anisotropic in the SSF states with and without repulsion. In the SSF with repulsion, the superfluid stiffness is subject to different types of anisotropy in the region near half filling and above 2/3-filling.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Conducting-angle-based percolation in the XY model

We define a percolation problem on the basis of spin configurations of the two dimensional XY model. Neighboring spins belong to the same percolation cluster if their orientations differ less than a certain threshold called the conducting angle. The percolation properties of this model are studied by means of Monte Carlo simulations and a finite-size scaling analysis. Our simulations show the existence of percolation transitions when the conducting angle is varied, and we determine the transition point for several values of the XY coupling. It appears that the critical behavior of this percolation model can be well described by the standard percolation theory. The critical exponents of the percolation transitions, as determined by finite-size scaling, agree with the universality class of the two-dimensional percolation model on a uniform substrate. This holds over the whole temperature range, even in the low-temperature phase where the XY substrate is critical in the sense that it displays algebraic decay of correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Crossing bonds in the random-cluster model

We derive the scaling dimension associated with crossing bonds in the random-cluster representation of the two-dimensional Potts model, by means of a mapping on the Coulomb gas. The scaling field associated with crossing bonds appears to be irrelevant, on the critical as well as on the tricritical branch. The latter result stands in a remarkable contrast with the existing result for the tricritical O(n) model that crossing bonds are relevant. In order to obtain independent confirmation of the Coulomb gas result for the crossing-bond exponent, we perform a finite-size-scaling analysis based on numerical transfer-matrix calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Percolation and critical O($n$) loop configurations

We study a percolation problem based on critical loop configurations of the O($n$) loop model on the honeycomb lattice. We define dual clusters as groups of sites on the dual triangular lattice that are not separated by a loop, and investigate the the bond-percolation properties of these dual clusters. The universal properties at the percolation threshold are argued to match those of Kasteleyn-Fortuin random clusters in the critical Potts model. This relation is checked numerically by means of cluster simulations of several O($n$) models in the range $1\leq n \leq 2$. The simulation results include the percolation threshold for several values of $n$, as well as the universal exponents associated with bond dilution and the size distribution of the diluted clusters at the percolation threshold. Our numerical results for the exponents are in agreement with existing Coulomb gas results for the random-cluster model, which confirms the relation between both models. We discuss the renormalization flow of the bond-dilution parameter $p$ as a function of $n$, and provide an expression that accurately describes a line of unstable fixed points as a function of $n$, corresponding with the percolation threshold. Furthermore, the renormalization scenario indicates the existence, in $p$ versus $n$ diagram, another line of fixed points at $p=1$, which is stable with respect to $p$.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Phase transitions in self-dual generalizations of the Baxter-Wu model

We study two types of generalized Baxter-Wu models, by means of transfer-matrix and Monte Carlo techniques. The first generalization allows for different couplings in the up- and down triangles, and the second generalization is to a $q$-state spin model with three-spin interactions. Both generalizations lead to self-dual models, so that the probable locations of the phase transitions follow. Our numerical analysis confirms that phase transitions occur at the self-dual points. For both generalizations of the Baxter-Wu model, the phase transitions appear to be discontinuous.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Critical properties of a dilute O($n$) model on the kagome lattice

A critical dilute O($n$) model on the kagome lattice is investigated analytically and numerically. We employ a number of exact equivalences which, in a few steps, link the critical O($n$) spin model on the kagome lattice to the exactly solvable critical $q$-state Potts model on the honeycomb lattice with $q=(n+1)^2$. The intermediate steps involve the random-cluster model on the honeycomb lattice, and a fully packed loop model with loop weight $n'=\sqrt{q}$ and a dilute loop model with loop weight $n$, both on the kagome lattice. This mapping enables the determination of a branch of critical points of the dilute O($n$) model, as well as some of its critical properties. For $n=0$, this model reproduces the known universal properties of the $θ$ point describing the collapse of a polymer. For $n\neq 0$ it displays a line of multicritical points, with the same universal properties as a branch of critical behavior that was found earlier in a dilute O($n$) model on the square lattice. These findings are supported by a finite-size-scaling analysis in combination with transfer-matrix calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Tricritical O(n) models in two dimensions

We show that the exactly solved low-temperature branch of the two-dimensional O($n$) model is equivalent with an O($n$) model with vacancies and a different value of $n$. We present analytic results for several universal parameters of the latter model, which is identified as a tricritical point. These results apply to the range $n \leq 3/2$, and include the exact tricritical point, the conformal anomaly and a number of scaling dimensions, among which the thermal and magnetic exponent, the exponent associated with crossover to ordinary critical behavior, and to tricritical behavior with cubic symmetry. We describe the translation of the tricritical model in a Coulomb gas. The results are verified numerically by means of transfer-matrix calculations. We use a generalized ADE model as an intermediary, and present the expression of the one-point distribution function in that language. The analytic calculations are done both for the square and the hexagonal lattice.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

First and second order transitions in dilute O(n) models

We explore the phase diagram of an O(n) model on the honeycomb lattice with vacancies, using finite-size scaling and transfer-matrix methods. We make use of the loop representation of the O(n) model, so that $n$ is not restricted to positive integers. For low activities of the vacancies, we observe critical points of the known universality class. At high activities the transition becomes first order. For n=0 the model includes an exactly known theta point, used to describe a collapsing polymer in two dimensions. When we vary $n$ from 0 to 1, we observe a tricritical point which interpolates between the universality classes of the theta point and the Ising tricritical point.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Phase Diagram of a Loop on the Square Lattice

The phase diagram of the O(n) model, in particular the special case $n=0$, is studied by means of transfer-matrix calculations on the loop representation of the O(n) model. The model is defined on the square lattice; the loops are allowed to collide at the lattice vertices, but not to intersect. The loop model contains three variable parameters that determine the loop density or temperature, the energy of a bend in a loop, and the interaction energy of colliding loop segments. A finite-size analysis of the transfer-matrix results yields the phase diagram in a special plane of the parameter space. These results confirm the existence of a multicritical point and an Ising-like critical line in the low-temperature O(n) phase.

0 citations0 reviewsNo code linkedOpen access