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

Chengxiang Ding

Chengxiang Ding 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.softquant-phhep-th

Catalog footprint

What is connected

18works
6topics
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

18 published item(s)

preprint2026arXiv

Phase transitions and crtical exponents in the six-vertex model on kagome lattices

Inspired by the experimental realization of direct kagome spin ice [Yue et al., Nat. Nanotechnol. 19, 1101 (2024)], the theoretical six-vertex model on the kagome lattice is systematically simulated using the directed loop Monte Carlo method. Four distinct vortex lattice phases are identified: (i) antiferromagnetic leg states and vortex lattice order on both triangular and honeycomb faces, with a winding number $k=1$. (ii) ferromagnetic leg states and vortex lattice order on both types of faces, with $k=-2$ on the honeycomb faces and $k=1$ on the triangular faces. (iii) paramagnetic leg states and vortex lattice order on the triangular faces with $k=1$; and (iv) paramagnetic leg states and vortex lattice order on the honeycomb faces with $k=1$. As for ferromagnetic to different types of paramagnetic phase, besides the Ising universality class with $y_t=1$, varying critical exponents have also been found with different values of vertex weights. The transition between the third type and fourth type of vortex lattice phases occurs with the new exponent $y_t=1.340(3)$. The third and fourth types of the vortex lattice phase to the vortex disorder phase are found to be of the Berezinskii-Kosterlitz-Thouless type. These findings contribute to the search for and understanding of ice on complex lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Universal Quench Dynamics of an Open Quantum System

Taking the quantum Kitaev chain as an example, we have studied the universal dynamical behaviors resulting from quantum criticality under the condition of environmental temperature quench. Our findings reveal that when the quantum parameter is at its critical value, both the excess excitation density at the end of linear quench and the subsequent free relaxation behavior exhibit universal scaling behaviors. The scaling laws observed upon quenching to the zero-temperature quantum critical point and non-zero temperature points exhibit distinct scaling exponents, which are all intimately related to the dynamical critical exponents of the quantum phase transition. Additionally, for the case of linear quench to finite temperatures, we have also discovered an intrinsic universal dynamical behavior that is independent of quantum criticality. Our research offers profound insights into the relationship between quantum criticality and nonequilibrium dynamics from two perspectives: Kibble-Zurek-like scaling behavior and free relaxation dynamics. Notably, the Kibble-Zurek-like scaling behavior in this context differs from the standard Kibble-Zurek mechanism. These two aspects jointly open up a new avenue for us to understand quantum criticality through real-time dynamical behavior, even at finite temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Finite-Size Scaling Theory at a Self-Dual Quantum Critical Point

The nondivergence of the generalized Grüneisen ratio (GR) at a quantum critical point (QCP) has been proposed to be a universal thermodynamic signature of self-duality. In this work, we study how the Kramers-Wannier-type self-duality manifests itself in the finite-size scaling behavior of thermodynamic quantities in the quantum critical regime. While the self-duality cannot be realized as a unitary transformation in the total Hilbert space for the Hamiltonian with the periodic boundary condition, it can be implemented in certain symmetry sectors with proper boundary conditions. Therefore, the GR and the transverse magnetization of the one-dimensional transverse-field Ising model exhibit different finite-size scaling behaviors in different sectors. This implies that the numerical diagnosis of self-dual QCP requires identifying the proper symmetry sectors.

0 citations0 reviewsNo code linkedOpen access
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 criticality of antiferromagnetic Potts model

We study the three-state antiferromagnetic Potts model on the simple-cubic lattice, paying attention to the surface critical behaviors. When the nearest neighboring interactions of the surface is tuned, we obtain a phase diagram similar to the XY model, owing to the emergent O(2) symmetry of the bulk critical point. For the ordinary transition, we get $y_{h1}=0.780(3)$, $η_\parallel=1.44(1)$, and $η_\perp=0.736(6)$; for the special transition, we get $y_s=0.59(1)$, $y_{h1}=1.693(2)$, $η_\parallel=-0.391(4)$, and $η_\perp=-0.179(5)$; in the extraordinary-log phase, the surface correlation function $C_\parallel(r)$ decays logarithmically, with decaying exponent $q=0.60(2)$, however, the correlation $C_\perp(r)$ still decays algebraically, with critical exponent $η_\perp=-0.442(5)$. If the ferromagnetic next nearest neighboring surface interactions are added, we find two transition points, the first one is a special point between the ordinary phase and the extraordinary-log phase, the second one is a transition between the extraordinary-log phase and the $Z_6$ symmetry-breaking phase, with critical exponent $y_{\rm s}=0.41(2)$. The scaling behaviors of the second transition is very interesting, the surface spin correlation function $C_\parallel(r)$ and the surface squared staggered magnetization at this point decays logarithmically, with exponent $q=0.37(1)$; however, the surface structure factor with the smallest wave vector and the correlation function $C_\perp(r)$ satisfy power-law decaying, with critical exponents $η_\parallel=-0.69(1)$ and $η_\perp=-0.37(1)$, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Phase transitions in 3D Ising model with cluster weight by Monte Carlo method

A cluster weight Ising model is proposed by introducing an additional cluster weight in the partition function of the traditional Ising model. It is equivalent to the O($n$) loop model or $n$-component face cubic loop model on the two-dimensional lattice, but on the three-dimensional lattice, it is still not very clear whether or not these models have the same universality. In order to simulate the cluster weight Ising model and search for new universality class, we apply a cluster algorithm, by combining the color-assignation and the Swendsen-Wang methods. The dynamical exponent for the absolute magnetization is estimated to be $z=0.45(3)$ at $n=1.5$, consistent with that of the traditional Swendsen-Wang methods. The numerical estimation of the thermal exponent $y_t$ and magnetic exponent $y_m$, show that the universalities of the two models on the three-dimensional lattice are different. We obtain the global phase diagram containing paramagnetic and ferromagnetic phases. The phase transition between the two phases are second order at $1\leq n< n_c$ and first order at $n\geq n_c$, where $n_c\approx 2$. The scaling dimension $y_t$ equals to the system dimension $d$ when the first-order transition occurs. Our results are helpful in the understanding of some traditional statistical mechanics models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamical quantum phase transition from a critical quantum quench

We study the dynamical quantum phase transition of the critical quantum quench, in which the prequenched Hamiltonian, or the postquenched Hamiltonian, or both of them are set to be the critical points of equilibrium quantum phase transitions, we find half-quantized or unquantized dynamical topological order parameter and dynamical Chern number; these results and also the existence of dynamical quantum phase transition are all closely related to the singularity of the Bogoliubov angle at the gap-closing momentum. The effects of the singularity may also be canceled out if both the prequenched and postquenched Hamiltonians are critical, then the dynamical topological order parameter and dynamical Chern number restore to integer ones. Our findings show that the widely accepted definitions of dynamical topological order parameter and dynamical Chern number are problematic for the critical quenches in the perspective of topology, which call for new definitions of them.

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
preprint2016arXiv

Emergent O(n) Symmetry in a series of three-dimensional Potts Models

We study the q-state Potts model on the simple cubic lattice with ferromagnetic interactions in one lattice direction, and antiferromagnetic interactions in the two other directions. As the temperature T decreases, the system undergoes a second-order phase transition that fits in the universality class of the 3D O(n) model with n=q-1. This conclusion is based on the estimated critical exponents, and histograms of the order parameter. At even smaller T we find, for q=4 and 5, a first-order transition to a phase with a different type of long-range order. This long-range order dissolves at T=0, and the system effectively reduces to a disordered two-dimensional Potts antiferromagnet. These results are obtained by means of Monte Carlo simulations and finite-size scaling.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Percolation of the Site Random-Cluster Model by Monte Carlo Method

Herein, we propose a site random cluster model by introducing an additional cluster weight in the partition function of the traditional site percolation. To simulate the model on a square lattice, we combine the color-assignation and the Swendsen-Wang methods together to design a highly efficient cluster algorithm with a small critical slowing-down phenomenon. To verify whether or not it is consistent with the bond random cluster model, we measure several quantities such as the wrapping probability $R_e$, the percolation strength $P_\infty$, and the magnetic susceptibility per site $χ_p$ as well as two exponents such as the thermal exponent $y_t$ and the fractal dimension $y_h$ of the largest cluster. We find that for different exponents of cluster weight q=1.5, 2, 2.5, 3, 3.5 and 4, the numerical estimation of the exponents $y_t$ and $y_h$ are consistent with the theoretical values. The universalities of the site random cluster model and the bond random cluster model are completely identical. For larger values of $q$, we find obvious signatures of the first-order percolation transition by the histograms and the hysteresis loops of the percolation strength and the energy per site. Our results are helpful for the understanding of the percolation of traditional statistical models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Trimer superfluid and supersolid on two-dimensional optical lattices

By the photoassociation method, the trimer superfluid phase disappears in the one dimensional state-dependent optical lattice if the ratio of the three-body interaction $W$ to the trimer tunneling $J$is kept at $W/J=12$ [Phys Rev A. {\bf 90}, 033622(2014)]. To search for a trimer superfluid and trimer supersolid, we load the cold atom into two-dimensional lattices, whose coordinate number $z$ and kinetic energy $-zJ$ are respectively larger and lower than those of a one dimensional lattice. Herein, we study the Bose-Hubbard model, which has an additional trimer tunneling term, a three-body interaction and a next-nearest repulsion. The on-site trimer and trimer superfluid exist if the on-site two-body repulsion and three-body repulsion are smaller than some thresholds. With atom-tunneling terms, the phase transitions from trimer superfluid phase to both atom superfluid and atom supersolid phases are first order. With $W/J=12$, in a one dimensional lattice, the trimer superfluid phase does not exist at all. In contrast, the trimer superfluid phase, exists in the lower density regions $0 \textless ρ\textless2$ on the square lattice if $J$ is not very large. The trimer superfluid phase emerges in a wider range $0 \textless ρ\textless3$ in the triangular lattice, or in the cubic lattice ($z=6$). When the three-body interaction is turned on, a trimer supersolid phase emerges due to the classical degeneracy between the quasi trimer solid and the trimer solid being broken by the quantum fluctuation. The phase transitions from the trimer supersolid phase to quasi trimer solid are first order and the phase transition from the trimer supersolid phase to trimer solid is continuous. Our results, obtained by mean-field and quantum Monte Carlo methods, will be helpful in realizing the trimer superfluid and supersolid by cold atom experiments.

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
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

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
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

Potts and percolation models on bowtie lattices

We give the exact critical frontier of the Potts model on bowtie lattices. For the case of $q=1$, the critical frontier yields the thresholds of bond percolation on these lattices, which are exactly consistent with the results given by Ziff et al [J. Phys. A 39, 15083 (2006)]. For the $q=2$ Potts model on the bowtie-A lattice, the critical point is in agreement with that of the Ising model on this lattice, which has been exactly solved. Furthermore, we do extensive Monte Carlo simulations of Potts model on the bowtie-A lattice with noninteger $q$. Our numerical results, which are accurate up to 7 significant digits, are consistent with the theoretical predictions. We also simulate the site percolation on the bowtie-A lattice, and the threshold is $s_c=0.5479148(7)$. In the simulations of bond percolation and site percolation, we find that the shape-dependent properties of the percolation model on the bowtie-A lattice are somewhat different from those of an isotropic lattice, which may be caused by the anisotropy of the lattice.

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
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