BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Youjin Deng

Youjin Deng contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.stat-mechcond-mat.dis-nncond-mat.softphysics.soc-phArtificial Intelligencecond-mat.quant-gascond-mat.str-elMachine Learningmath-phmath.MPphysics.data-anPopulations and Evolution

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
12topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

22 published item(s)

preprint2026arXiv

Inverse Knowledge Search over Verifiable Reasoning: Synthesizing a Scientific Encyclopedia from a Long Chains-of-Thought Knowledge Base

Most scientific materials compress reasoning, presenting conclusions while omitting the derivational chains that justify them. This compression hinders verification by lacking explicit, step-wise justifications and inhibits cross-domain links by collapsing the very pathways that establish the logical and causal connections between concepts. We introduce a scalable framework that decompresses scientific reasoning, constructing a verifiable Long Chain-of-Thought (LCoT) knowledge base and projecting it into an emergent encyclopedia, SciencePedia. Our pipeline operationalizes an endpoint-driven, reductionist strategy: a Socratic agent, guided by a curriculum of around 200 courses, generates approximately 3 million first-principles questions. To ensure high fidelity, multiple independent solver models generate LCoTs, which are then rigorously filtered by prompt sanitization and cross-model answer consensus, retaining only those with verifiable endpoints. This verified corpus powers the Brainstorm Search Engine, which performs inverse knowledge search -- retrieving diverse, first-principles derivations that culminate in a target concept. This engine, in turn, feeds the Plato synthesizer, which narrates these verified chains into coherent articles. The initial SciencePedia comprises approximately 200,000 fine-grained entries spanning mathematics, physics, chemistry, biology, engineering, and computation. In evaluations across six disciplines, Plato-synthesized articles (conditioned on retrieved LCoTs) exhibit substantially higher knowledge-point density and significantly lower factual error rates than an equally-prompted baseline without retrieval (as judged by an external LLM). Built on this verifiable LCoT knowledge base, this reasoning-centric approach enables trustworthy, cross-domain scientific synthesis at scale and establishes the foundation for an ever-expanding encyclopedia.

0 citations0 reviewsNo code linkedOpen access
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
preprint2022arXiv

Backbone and shortest-path exponents of the two-dimensional $Q$-state Potts model

We present a Monte Carlo study of the backbone and the shortest-path exponents of the two-dimensional $Q$-state Potts model in the Fortuin-Kasteleyn bond representation. We first use cluster algorithms to simulate the critical Potts model on the square lattice and obtain the backbone exponents $d_{\rm B} = 1.732 \, 0(3)$ and $1.794(2)$ for $Q=2,3$ respectively. However, for large $Q$, the study suffers from serious critical slowing down and slowly converging finite-size corrections. To overcome these difficulties, we consider the O$(n)$ loop model on the honeycomb lattice in the densely packed phase, which is regarded to correspond to the critical Potts model with $Q=n^2$. With a highly efficient cluster algorithm, we determine from domains enclosed by the loops $d_{\rm B} =1.643\,39(5), 1.732\,27(8), 1.793\,8(3), 1.838\,4(5), 1.875\,3(6)$ for $Q = 1, 2, 3, 2 \! + \! \sqrt{3}, 4$, respectively, and $d_{\rm min} = 1.094\,5(2), 1.067\,5(3), 1.047\,5(3), 1.032\,2(4)$ for $Q=2,3, 2+\sqrt{3}, 4$ respectively. Our estimates significantly improve over the existing results for both $d_{\rm B}$ and $d_{\rm min}$. Finally, by studying finite-size corrections in backbone-related quantities, we conjecture an exact formula as a function of $n$ for the leading correction exponent.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Geometric Upper Critical Dimensions of the Ising Model

The upper critical dimension of the Ising model is known to be $d_c=4$, above which critical behavior is regarded as trivial. We hereby argue from extensive simulations that, in the random-cluster representation, the Ising model simultaneously exhibits two upper critical dimensions at $(d_c= 4, d_p=6)$, and critical clusters for $d \geq d_p$, except the largest one, are governed by exponents from percolation universality. We predict a rich variety of geometric properties and then provide strong evidence in dimensions from 4 to 7 and on complete graphs. Our findings significantly advance the understanding of the Ising model, which is a fundamental system in many branches of physics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Machine learning for percolation utilizing auxiliary Ising variables

Machine learning for phase transition has received intensive research interest in recent years. However, its application in percolation still remains challenging. We propose an auxiliary Ising mapping method for machine learning study of the standard percolation as well as a variety of statistical mechanical systems in correlated percolation representations. We demonstrate that unsupervised machine learning is able to accurately locate the percolation threshold, independent of the spatial dimension of system or the type of phase transition, which can be first order or continuous. Moreover, we show that, by neural network machine learning, auxiliary Ising configurations for different universalities can be classified with high confidence level. Our results indicate that the auxiliary Ising mapping method, despite of it simplicity, can advance the application of machine learning in statistical and condensed-matter physics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Monte Carlo study of duality and the Berezinskii-Kosterlitz-Thouless phase transitions of the two-dimensional $q$-state clock model in flow representations

The two-dimensional $q$-state clock model for $q \geq 5$ undergoes two Berezinskii-Kosterlitz-Thouless (BKT) phase transitions as temperature decreases. Here we report an extensive worm-type simulation of the square-lattice clock model for $q=$5--9 in a pair of flow representations, from the high- and low-temperature expansions, respectively. By finite-size scaling analysis of susceptibility-like quantities, we determine the critical points with a precision improving over the existing results. Due to the dual flow representations, each point in the critical region is observed to simultaneously exhibit a pair of anomalous dimensions, which are $η_1=1/4$ and $η_2 = 4/q^2$ at the two BKT transitions. Further, the approximate self-dual points $β_{\rm sd}(L)$, defined by the stringent condition that the susceptibility like quantities in both flow representations are identical, are found to be nearly independent of system size $L$ and behave as $β_{\rm sd} \simeq q/2π$ asymptotically at the large-$q$ limit. The exponent $η$ at $β_{\rm sd}$ is consistent with $1/q$ within statistical error as long as $q \geq 5$. Based on this, we further conjecture that $η(β_{\rm sd}) = 1/q$ holds exactly and is universal for systems in the $q$-state clock universality class. Our work provides a vivid demonstration of rich phenomena associated with the duality and self-duality of the clock model in two dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multiple and Asymmetric Scalings in Explosive Percolation

Explosive percolation in the Achlioptas process has recently attracted much research attention. From extensive simulations in an event-based ensemble, we find that, in dimensions from $2$ to $6$ and on random graphs, the Achlioptas processes all have two scaling windows and multiple fractal structures. The mixing of these multiple scalings successfully explains the previously observed anomalous phenomena in the conventional ensemble, and, moreover, correct critical exponents are now determined with a high precision by the event-based method. The multiple scalings and the ensemble inequivalence may bring new insights for other statistical systems.

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
preprint2022arXiv

Universal Critical Behavior of Percolation in Orientationally Ordered Janus Particles and Other Anisotropic Systems

We combine percolation theory and Monte Carlo simulation to study in two dimensions the connectivity of an equilibrium lattice model of interacting Janus disks which self-assemble into an orientationally ordered stripe phase at low temperature. As the patch size is increased or the temperature is lowered, clusters of patch-connected disks grow, and a percolating cluster emerges at a threshold. In the stripe phase, the critical clusters extend longer in the direction parallel to the stripes than in the perpendicular direction, and percolation is thus anisotropic. It is found that the critical behavior of percolation in the Janus system is consistent with that of standard isotropic percolation, when an appropriate spatial rescaling is made. The rescaling procedure can be applied to understand other anisotropic systems, such as the percolation of aligned rigid rods and of the $q$-state Potts model with anisotropic interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Unwrapped two-point functions on high-dimensional tori

We study unwrapped two-point functions for the Ising model, the self-avoiding walk and a random-length loop-erased random walk on high-dimensional lattices with periodic boundary conditions. While the standard two-point functions of these models have been observed to display an anomalous plateau behaviour, the unwrapped two-point functions are shown to display standard mean-field behaviour. Moreover, we argue that the asymptotic behaviour of these unwrapped two-point functions on the torus can be understood in terms of the standard two-point function of a random-length random walk model on Zd. A precise description is derived for the asymptotic behaviour of the latter. Finally, we consider a natural notion of the Ising walk length, and show numerically that the Ising and SAW walk lengths on high-dimensional tori show the same universal behaviour known for the SAW walk length on the complete graph.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Nested Closed Paths in Two-Dimensional Percolation

For two-dimensional percolation on a domain with the topology of a disc, we introduce a nested-path operator (NP) and thus a continuous family of one-point functions $W_k \equiv \langle \mathcal{R} \cdot k^\ell \rangle $, where $\ell$ is the number of independent nested closed paths surrounding the center, $k$ is a path fugacity, and $\mathcal{R}$ projects on configurations having a cluster connecting the center to the boundary. At criticality, we observe a power-law scaling $W_k \sim L^{X_{\rm NP}}$, with $L$ the linear system size, and we determine the exponent $X_{\rm NP}$ as a function of $k$. On the basis of our numerical results, we conjecture an analytical formula, $X_{\rm NP} (k) = \frac{3}{4}ϕ^2 -\frac{5}{48}ϕ^2/ (ϕ^2-\frac{2}{3})$ where $k = 2 \cos(πϕ)$, which reproduces the exact results for $k=0,1$ and agrees with the high-precision estimate of $X_{\rm NP}$ for other $k$ values. In addition, we observe that $W_2(L)=1$ for site percolation on the triangular lattice with any size $L$, and we prove this identity for all self-matching lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Percolation effects in the Fortuin-Kasteleyn Ising model on the complete graph

The Fortuin-Kasteleyn (FK) random cluster model, which can be exactly mapped from the $q$-state Potts spin model, is a correlated bond percolation model. By extensive Monte Carlo simulations, we study the FK bond representation of the critical Ising model ($q=2$) on a finite complete graph, i.e. the mean-field Ising model. We provide strong numerical evidence that the configuration space for $q=2$ contains an asymptotically vanishing sector in which quantities exhibit the same finite-size scaling as in the critical uncorrelated bond percolation ($q=1$) on the complete graph. Moreover, we observe that in the full configuration space, the power-law behaviour of the cluster-size distribution for the FK Ising clusters except the largest one is governed by a Fisher exponent taking the value for $q=1$ instead of $q=2$. This demonstrates the percolation effects in the FK Ising model on the complete graph.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Two-scale scenario of rigidity percolation of sticky particles

In the presence of attraction, the jamming transition of packings of frictionless particles corresponds to the rigidity percolation. When the range of attraction is long, the distribution of the size of rigid clusters, $P(s)$, is continuous and shows a power-law decay. For systems with short-range attractions, however, $P(s)$ appears discontinuous. There is a power-law decay for small cluster sizes, followed by a low probability gap and a peak near the system size. We find that this appearing ``discontinuity'' does not mean that the transition is discontinuous. In fact, it signifies the coexistence of two distinct length scales, associated with the largest cluster and smaller ones, respectively. The comparison between the largest and second largest clusters indicates that their growth rates with system size are rather different. However, both cluster sizes tend to diverge in the large system size limit, suggesting that the jamming transition of systems with short-range attractions is still continuous. In the framework of the two-scale scenario, we also derive a generalized hyperscaling relation. With robust evidence, our work challenges the former single-scale view of the rigidity percolation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

$N$-cluster correlations in four- and five-dimensional percolation

We study $N$-cluster correlation functions in four- and five-dimensional (4D, 5D) bond percolation by extensive Monte Carlo simulation. We reformulate the transfer Monte Carlo algorithm for percolation [Phys. Rev. E {\bf 72}, 016126 (2005)] using the disjoint-set data structure, and simulate a cylindrical geometry $L^{d-1}\times \infty$, with the linear size up to $L=512$ for 4D and $128$ for 5D. We determine with a high precision all possible $N$-cluster exponents, for $N \! =\!2$ and $3$, and the universal amplitude for a logarithmic correlation function. From the symmetric correlator with $N \! = \!2$, we obtain the correlation-length critical exponent as $1/ν\! =\! 1.4610(12)$ for 4D and $1/ν\! =\! 1.737 (2)$ for 5D, significantly improving over the existing results. Estimates for the other exponents and the universal logarithmic amplitude have not been reported before to our knowledge. Our work demonstrates the validity of logarithmic conformal field theory and adds to the growing knowledge for high-dimensional percolation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Complete graph and Gaussian fixed point asymptotics in the five-dimensional Fortuin-Kasteleyn Ising model with periodic boundaries

We present an extensive Markov-chain Monte Carlo study of the finite-size scaling behavior of the Fortuin-Kasteleyn Ising model on five-dimensional hypercubic lattices with periodic boundary conditions. We observe that physical quantities, which include the contribution of the largest cluster, exhibit complete graph asymptotics. However, for quantities, where the contribution of the largest cluster is removed, we observe that the scaling behavior is mainly controlled by the Gaussian fixed point. Our results therefore suggest that \textit{both} scaling predictions, i.e. the complete graph \textit{and} the Gaussian fixed point asymptotics, are needed to provide a complete description for the five-dimensional finite-size scaling behavior on the torus.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Critical polynomials in the nonplanar and continuum percolation models

Exact or precise thresholds have been intensively studied since the introduction of the percolation model. Recently the critical polynomial $P_{\rm B}(p,L)$ was introduced for planar-lattice percolation models, where $p$ is the occupation probability and $L$ is the linear system size. The solution of $P_{\rm B} = 0$ can reproduce all known exact thresholds and leads to unprecedented estimates for thresholds of unsolved planar-lattice models. In two dimensions, assuming the universality of $P_{\rm B}$, we use it to study a nonplanar lattice model, i.e., the equivalent-neighbor lattice bond percolation, and the continuum percolation of identical penetrable disks, by Monte Carlo simulations and finite-size scaling analysis. It is found that, in comparison with other quantities, $P_{\rm B}$ suffers much less from finite-size corrections. As a result, we obtain a series of high-precision thresholds $p_c(z)$ as a function of coordination number $z$ for equivalent-neighbor percolation with $z$ up to O$(10^5)$, and clearly confirm the asymptotic behavior $zp_c-1 \sim 1/\sqrt{z}$ for $z \rightarrow \infty$. For the continuum percolation model, we surprisingly observe that the finite-size correction in $P_{\rm B}$ is unobservable within uncertainty O$(10^{-5})$ as long as $L \geq 3$. The estimated threshold number density of disks is $ρ_c = 1.436 325 05(10)$, slightly below the most recent result $ρ_c = 1.436 325 45(8)$ of Mertens and Moore obtained by other means. Our work suggests that the critical polynomial method can be a powerful tool for studying nonplanar and continuum systems in statistical mechanics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Extended crossover from Fermi liquid to quasi-antiferromagnet in the half-filled 2D Hubbard model

The ground state of the Hubbard model with nearest-neighbor hopping on the square lattice at half filling is known to be that of an antiferromagnetic (AFM) band insulator for any on-site repulsion. At finite temperature, the absence of long-range order makes the question of how the interaction-driven insulator is realized nontrivial. We address this problem with controlled accuracy in the thermodynamic limit using self-energy diagrammatic determinant Monte Carlo and dynamical cluster approximation methods and show that development of long-range AFM correlations drives an extended crossover from Fermi liquid to insulating behavior in the parameter regime that precludes a metal-to-insulator transition. The intermediate crossover state is best described as a non-Fermi liquid with a partially gapped Fermi surface.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

History-dependent percolation on multiplex networks

The structure of interconnected systems and its impact on the system dynamics is a much-studied cross-disciplinary topic. Although various critical phenomena have been found in different models, the study on the connections between different percolation transitions is still lacking. Here we propose a unified framework to study the origins of the discontinuous transitions of the percolation process on interacting networks. The model evolves in generations with the result of the present percolation depending on the previous state and thus is history-dependent. Both theoretical analysis and Monte Carlo simulations reveal that the nature of the transition remains the same at finite generations but exhibits an abrupt change for the infinite generation. We use brain functional correlation and morphological similarity data to show that our model also provides a general method to explore the network structure and can contribute to many practical applications, such as detecting the abnormal structures of human brain networks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Oscillatory-like Expansion of a Fermionic Superfluid

We study the expansion behaviours of a Fermionic superfluid in a cigar-shaped optical dipole trap for the whole BEC-BCS crossover and various temperatures. At low temperature ($0.06(1) T_F$), the atom cloud undergoes an anisotropic hydrodynamic expansion over 30~ms, which behaves like oscillation in the horizontal plane. By analyzing the expansion dynamics according to the superfluid hydrodynamic equation, the effective polytropic index $\barγ$ of Equation-of-State of Fermionic superfluid is extracted. The $\barγ$ values show a non-monotonic behavior over the BEC-BCS crossover, and have a good agreement with the theoretical results in the unitarity and BEC side. The normalized quasi-frequencies of the oscillatory expansion are measured, which drop significantly from the BEC side to the BCS side and reach a minimum value of 1.73 around $1/k_Fa=-0.25$. Our work improves the understanding of the dynamic properties of strongly interacting Fermi gas.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Scaling features in the spreading of COVID-19

Since the outbreak of COVID-19, many data analyses have been done. Some of them are based on the classical epidemiological approach that assumes an exponential growth, but a few studies report that a power-law scaling may provide a better fit to the currently available data. Hereby, we examine the data in China (01/20/2020--02/24/2020), and indeed find that the growth closely follows a power-law kinetics over a significantly wide time period. The exponents are $2.48(20)$, $2.21(6)$ and $4.26(12)$ for the number of confirmed infections, deaths and cured cases, respectively, indicating an underlying small-world network structure in the pandemic. While no obvious deviations from the power-law growth can be seen yet for the number of deaths and cured cases, negative deviations have clearly appeared in the number of infections, particularly that for the region outside Hubei. This suggests the beginning of the slowing-down of the virus spreading due to the huge containment effort. Meanwhile, we find that despite the dramatic difference in magnitudes, the growth kinetics of the infection number exhibits much similarity for Hubei province and the region outside Hubei. On this basis, in log-log plot, we rescale the infection number for the region outside Hubei such that it overlaps as much as possible with the total infection number in China, from which an approximate extrapolation yields the maximum of the pandemic around March 3, 2020, with the number of infections about $83,000$. Further, by analyzing the kinetics of the mortality in log-log scale, we obtains a rough estimate that near March 3, the death rate of COVID-19 would be about $4.7\%\thicksim 5.0\%$ for Hubei province and $0.7\%\thicksim1.0\%$ for the region outside Hubei. We emphasize that our predictions may be quantitatively unreliable, since the data analysis is purely empirical and various assumptions are used.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Worm-algorithm-type Simulation of Quantum Transverse-Field Ising Model

We apply a worm algorithm to simulate the quantum transverse-field Ising model in a path-integral representation of which the expansion basis is taken as the spin component along the external-field direction. In such a representation, a configuration can be regarded as a set of non-intersecting loops constructed by &#34;kinks&#34; for pairwise interactions and spin-down (or -up) imaginary-time segments. The wrapping probability for spin-down loops, a dimensionless quantity characterizing the loop topology on a torus, is observed to exhibit small finite-size corrections and yields a high-precision critical point in two dimensions (2D) as $h_c \! =\! 3.044\, 330(6)$, significantly improving over the existing results and nearly excluding the best one $h_c \! =\! 3.044\, 38 (2)$. At criticality, the fractal dimensions of the loops are estimated as $d_{\ell \downarrow} (1{\rm D}) \! = \! 1.37(1) \! \approx \! 11/8 $ and $d_{\ell \downarrow} (2{\rm D}) \! = \! 1.75 (3)$, consistent with those for the classical 2D and 3D O(1) loop model, respectively. An interesting feature is that in one dimension (1D), both the spin-down and -up loops display the critical behavior in the whole disordered phase ($ 0 \! \leq \! h \! < \! h_c$), having a fractal dimension $d_{\ell} \! = \! 1.750 (7)$ that is consistent with the hull dimension $d_{\rm H} \! = \! 7/4$ for critical 2D percolation clusters. The current worm algorithm can be applied to simulate other quantum systems like hard-core boson models with pairing interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

The elastic and directed percolation backbone

We argue that the elastic backbone (EB) (union of shortest paths) on a cylindrical system, recently studied by Sampaio Filho et al. [Phys. Rev. Lett. 120, 175701 (2018)], is in fact the backbone of two-dimensional directed percolation (DP). We simulate the EB on the same system as considered by these authors, and also study the DP backbone directly using an algorithm that allows backbones to be generated in a completely periodic manner. We find that both the EB in the bulk and the DP backbone have a fractal dimension of $d_{b} = d_{B,\rm DP} = 1.681\,02(15)$ at the identical critical point $p_{c,\rm{DP}} \approx 0.705\,485\,22$. We also measure the fractal dimension at the edge of the EB system and for the full DP clusters, and find $d_e = d_{\rm DP} = 1.840 \, 54 (4)$. We argue that those two fractal dimensions follow from the DP exponents as $d_{B,\rm DP} = 2-2β/ν_\parallel = 1.681 \, 07 2 (12) $ and $d_{\rm DP} = 2-β/ν_\parallel = 1.840\, 536 (6)$. Our fractal dimensions differ from the value 1.750(3) found by Sampaio Filho et al., whose value may represent a crossover effect.

0 citations0 reviewsNo code linkedOpen access