Source author record

Youjin Deng

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

79works

25topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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 "kinks" 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.

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.

preprint2017arXiv

Irreversible Markov chain Monte Carlo algorithm for self-avoiding walk

We formulate an irreversible Markov chain Monte Carlo algorithm for the self-avoiding walk (SAW), which violates the detailed balance condition and satisfies the balance condition. Its performance improves significantly compared to that of the Berretti-Sokal algorithm, which is a variant of the Metropolis-Hastings method. The gained efficiency increases with the spatial dimension (D), from approximately $10$ times in 2D to approximately $40$ times in 5D. We simulate the SAW on a 5D hypercubic lattice with periodic boundary conditions, for a system with a linear size up to $L=128$, and confirm that as for the 5D Ising model, the finite-size scaling of the SAW is governed by renormalized exponents $ν^*=2/d$ and $γ/ν^*=d/2$. The critical point is determined, which is approximately $8$ times more precise than the best available estimate.

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.

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.

preprint2016arXiv

No-enclave percolation corresponds to holes in the cluster backbone

The no-enclave percolation (NEP) model introduced recently by Sheinman et al. can be mapped to a problem of holes within a standard percolation backbone, and numerical measurements of these holes gives the size-distribution exponent $τ= 1.82(1)$ of the NEP model. An argument is given that $τ=1 + d_B/2 \approx 1.822$ where $d_B$ is the backbone dimension. On the other hand, a model of simple holes within a percolation cluster implies $τ= 1 + d_f/2 = 187/96 \approx 1.948$, where $d_f$ is the fractal dimension of the cluster, and this value is consistent with Sheinman et al.'s experimental results of gel collapse which gives $τ= 1.91(6)$. Both models yield a discontinuous maximum hole size at $p_c$, signifying explosive percolation behavior. At $p_c$, the largest hole fills exactly half the system, due to symmetry. Extensive numerical simulations confirm our results.

preprint2016arXiv

Spin-Ice State of the Quantum Heisenberg Antiferromagnet on the Pyrochlore Lattice

We study the low-temperature physics of the $SU(2)$-symmetric spin-$1/2$ Heisenberg antiferromagnet on a pyrochlore lattice and find fingerprint evidence for the thermal spin-ice state in this frustrated quantum magnet. Our conclusions are based on the results of bold diagrammatic Monte Carlo simulations, with good convergence of the skeleton series down to the temperature $T/J=1/6$. The identification of the spin-ice state is done through a remarkably accurate microscopic correspondence for the static structure factor between the quantum Heisenberg, classical Heisenberg, and Ising models at all accessible temperatures, and the characteristic bowtie pattern with pinch points observed at $T/J=1/6$. The dynamic structure factor at real frequencies (obtained by the analytic continuation of numerical data) is consistent with diffusive spinon dynamics at the pinch points.

preprint2016arXiv

Trapping Centers at the Superfluid-Mott-insulator Criticality: Transition between Charge-quantized States

Under the conditions of superfluid-Mott-insulator criticality in two dimensions, the trapping centers--i.e., local potential wells and bumps--are generically characterized by an integer charge corresponding to the number of trapped particles (if positive) or holes (if negative). Varying the strength of the center leads to a transition between two competing ground states with charges differing by $\pm 1$. The hallmark of the transition scenario is a splitting of the number density distortion, $δn(r)$, into a half-integer core and a large halo carrying the complementary charge of $\pm 1/2$. The sign of the halo changes across the transition and the radius of the halo, $r_0$, diverges on the approach to the critical strength of the center, $V = V_c$, by the law $r_0 \propto |V-V_c|^{-\tildeν}$, with $\tildeν \approx 2.33(5)$.

preprint2015arXiv

Adaptive Multi-GPU Exchange Monte Carlo for the 3D Random Field Ising Model

We present an adaptive multi-GPU Exchange Monte Carlo method designed for the simulation of the 3D Random Field Model. The algorithm design is based on a two-level parallelization scheme that allows the method to scale its performance in the presence of faster and GPUs as well as multiple GPUs. The set of temperatures is adapted according to the exchange rate observed from short trial runs, leading to an increased exchange rate at zones where the exchange process is sporadic. Performance results show that parallel tempering is an ideal strategy for being implemented on the GPU, and runs between one to two orders of magnitude with respect to a single-core CPU version, with multi-GPU scaling being approximately $99\%$ efficient. The results obtained extend the possibilities of simulation to sizes of $L = 32, 64$ for a workstation with two GPUs.

preprint2015arXiv

Clique percolation in random graphs

As a generation of the classical percolation, clique percolation focuses on the connection of cliques in a graph, where the connection of two $k$-cliques means that they share at least $l<k$ vertices. In this paper, we develop a theoretical approach to study clique percolation in Erdős-Rényi graphs, which gives not only the exact solutions of the critical point, but also the corresponding order parameter. Based on this, we prove theoretically that the fraction $ψ$ of cliques in the giant clique cluster always makes a continuous phase transition as the classical percolation. However, the fraction $ϕ$ of vertices in the giant clique cluster for $l>1$ makes a step-function-like discontinuous phase transition in the thermodynamic limit and a continuous phase transition for $l=1$. More interesting, our analysis shows that at the critical point, the order parameter $ϕ_c$ for $l>1$ is neither $0$ nor $1$, but a constant depending on $k$ and $l$. All these theoretical findings are in agreement with the simulation results, which give theoretical support and clarification for previous simulation studies of clique percolation.

preprint2015arXiv

Coulomb liquid phases of bosonic cluster Mott insulators on a pyrochlore lattice

Employing large-scale quantum Monte Carlo simulations, we reveal the full phase diagram of the extended Hubbard model of hard-core bosons on the pyrochlore lattice with partial fillings. When the intersite repulsion is dominant, the system is in a cluster Mott insulator phase with an integer number of bosons localized inside the tetrahedral units of the pyrochlore lattice. We show that the full phase diagram contains three cluster Mott insulator phases with 1/4, 1/2, and 3/4 boson fillings, respectively. We further demonstrate that all three cluster Mott insulators are Coulomb liquid phases and its low-energy property is described by the emergent compact U(1) quantum electrodynamics. In addition to measuring the specific heat and entropy of the cluster Mott insulators, we investigate the correlation function of the emergent electric field and verify it is consistent with the compact U(1) quantum electrodynamics description. Our result sheds light on the magnetic properties of various pyrochlore systems, as well as the charge physics of the cluster magnets.

preprint2015arXiv

Emergent BCS regime of the two-dimensional fermionic Hubbard model: ground-state phase diagram

A significant part of the phase diagram of the two-dimensional fermionic Hubbard model for moderate interactions and filling factors ($U < 4, \, n<0.7$) is governed by effective Fermi liquid physics with weak BCS-type instabilities. We access this regime in a controlled way by a combination of the bold-line diagrammatic Monte Carlo method with an additional ladder-diagram summation trick and semi-analytic treatment of the weak instability in the Cooper channel. We obtain the corresponding ground-state phase diagram in the $(n,U)$ plane describing the competition between the $p-$ and $d-$wave superfluid states. We also claim the values of the dimensionless BCS coupling constants controlling the superfluid $T_c$ at the phase boundaries, which prove to be very small up to $U=4, n = 0.6$.

preprint2015arXiv

Ground state phase diagram of the repulsive fermionic $t-t^{\prime}$ Hubbard model on the square lattice from weak-coupling

We obtain a complete and exact in the weak-coupling limit ($U \rightarrow 0$) ground state phase diagram of the repulsive fermionic Hubbard model on the square lattice for filling factors $0 < n < 2$ and next-nearest-neighbour hopping amplitudes $0 \le t^{\prime} \le 0.5$. Phases are distinguished by the symmetry and the number of nodes of the superfluid order parameter. The phase diagram is richer than may be expected and typically features states with a high --- higher than that of the fundamental mode of the corresponding irreducible representation --- number of nodes. The effective coupling strength in the Cooper channel $λ$, which determines the critical temperature $T_c$ of the superfluid transition, is calculated in the whole parameter space and regions with high values of $λ$ are identified. It is shown that besides the expected increase of $λ$ near the Van Hove singularity line, joining the ferromagnetic and antiferromagnetic points, another region with high values of $λ$ can be found at quarter filling and $t^{\prime}=0.5$ due to the presence of a line of nesting at $t^{\prime} \ge 0.5$. The results can serve as benchmarks for controlled non-perturbative methods and guide the ongoing search for high-$T_c$ superconductivity in the Hubbard model.

preprint2015arXiv

Leaf-excluded percolation in two and three dimensions

We introduce the \emph{leaf-excluded} percolation model, which corresponds to independent bond percolation conditioned on the absence of leaves (vertices of degree one). We study the leaf-excluded model on the square and simple-cubic lattices via Monte Carlo simulation, using a worm-like algorithm. By studying wrapping probabilities, we precisely estimate the critical thresholds to be $0.355\,247\,5(8)$ (square) and $0.185\,022(3)$ (simple-cubic). Our estimates for the thermal and magnetic exponents are consistent with those for percolation, implying that the phase transition of the leaf-excluded model belongs to the standard percolation universality class.

preprint2015arXiv

Realization of Two-Dimensional Spin-orbit Coupling for Bose-Einstein Condensates

Cold atoms with laser-induced spin-orbit (SO) interactions provide intriguing new platforms to explore novel quantum physics beyond natural conditions of solids. Recent experiments demonstrated the one-dimensional (1D) SO coupling for boson and fermion gases. However, realization of 2D SO interaction, a much more important task, remains very challenging. Here we propose and experimentally realize, for the first time, 2D SO coupling and topological band with $^{87}$Rb degenerate gas through a minimal optical Raman lattice scheme, without relying on phase locking or fine tuning of optical potentials. A controllable crossover between 2D and 1D SO couplings is studied, and the SO effects and nontrivial band topology are observed by measuring the atomic cloud distribution and spin texture in the momentum space. Our realization of 2D SO coupling with advantages of small heating and topological stability opens a broad avenue in cold atoms to study exotic quantum phases, including the highly-sought-after topological superfluid phases.

preprint2015arXiv

Recursive Percolation

We introduce a simple lattice model in which percolation is constructed on top of critical percolation clusters, and show that it can be repeated recursively any number $n$ of generations. In two dimensions, we determine the percolation thresholds up to $n=5$. The corresponding critical clusters become more and more compact as $n$ increases, and define universal scaling functions of the standard two-dimensional form and critical exponents that are distinct for any $n$. This family of exponents differs from any previously known universality class, and cannot be accommodated by existing analytical methods. We confirm that recursive percolation is well defined also in three dimensions.

preprint2015arXiv

Short-range correlations in percolation at criticality

We derive the critical nearest-neighbor connectivity $g_n$ as $3/4$, $3(7-9p_c^{tri})/[4(5-4p_c^{tri})]$, and $3(2+7p_c^{tri})/[4(5-p_c^{tri})]$ for bond percolation on the square, honeycomb and triangular lattice respectively, where $p_c^{tri}=2\sin(π/18)$ is the percolation threshold for the triangular lattice; and confirm these values via Monte Carlo simulations. On the square lattice, we also numerically determine the critical next-nearest-neighbor connectivity as $g_{nn}=0.687\;500\;0(2)$, which confirms a conjecture by Mitra and Nienhuis in J. Stat. Mech. P10006 (2004), implying the exact value $g_{nn}=11/16$. We also determine the connectivity on a free surface as $g_n^{surf}=0.625\;000\;1(13)$ and conjecture that this value is exactly equal to $5/8$. In addition, we find that at criticality, the connectivities depend on the linear finite size L as $\sim L^{y_t-d}$, and the associated specific-heat-like quantities $C_n$ and $C_{nn}$ scale as $\sim L^{2y_t-d} \ln (L/L_0)$, where $d$ is the lattice dimensionality, $y_t=1/ν$ the thermal renormalization exponent, and $L_0$ a non-universal constant. We provide an explanation of this logarithmic factor in the theoretical framework reported recently by Vasseur et al. in J. Stat. Mech. L07001 (2012).

preprint2015arXiv

Solutions of the Two Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms

Numerical results for ground state and excited state properties (energies, double occupancies, and Matsubara-axis self energies) of the single-orbital Hubbard model on a two-dimensional square lattice are presented, in order to provide an assessment of our ability to compute accurate results in the thermodynamic limit. Many methods are employed, including auxiliary field quantum Monte Carlo, bare and bold-line diagrammatic Monte Carlo, method of dual fermions, density matrix embedding theory, density matrix renormalization group, dynamical cluster approximation, diffusion Monte Carlo within a fixed node approximation, unrestricted coupled cluster theory, and multi-reference projected Hartree-Fock. Comparison of results obtained by different methods allows for the identification of uncertainties and systematic errors. The importance of extrapolation to converged thermodynamic limit values is emphasized. Cases where agreement between different methods is obtained establish benchmark results that may be useful in the validation of new approaches and the improvement of existing methods.

preprint2015arXiv

The Massive Goldstone (Higgs) mode in two-dimensional ultracold atomic lattice systems

We discuss how to reveal the massive Goldstone mode, often referred to as the Higgs amplitude mode, near the Superfluid-to-Insulator quantum critical point (QCP) in a system of two-dimensional ultracold bosonic atoms in optical lattices. The spectral function of the amplitude response is obtained by analytic continuation of the kinetic energy correlation function calculated by Monte Carlo methods. Our results enable a direct comparison with the recent experiment [M. Endres, T. Fukuhara, D. Pekker, M. Cheneau, P. Schauß, C. Gross, E. Demler, S. Kuhr, and I. Bloch, Nature 487, 454-458 (2012)], and demonstrate a good agreement for temperature shifts induced by lattice modulation. Based on our numerical analysis, we formulate the necessary conditions in terms of homogeneity, detuning from the QCP and temperature in order to reveal the massive Goldstone resonance peak in spectral functions experimentally. We also propose to apply a local modulation at the trap center to overcome the inhomogeneous broadening caused by the parabolic trap confinement.

preprint2015arXiv

Universal Critical Wrapping Probabilities in the Canonical Ensemble

Universal dimensionless quantities, such as Binder ratios and wrapping probabilities, play an important role in the study of critical phenomena. We study the finite-size scaling behavior of the wrapping probability for the Potts model in the random-cluster representation, under the constraint that the total number of occupied bonds is fixed, so that the canonical ensemble applies. We derive that, in the limit $L \rightarrow \infty$, the critical values of the wrapping probability are different from those of the unconstrained model, i.e. the model in the grand-canonical ensemble, but still universal, for systems with $2y_t - d > 0$ where $y_t = 1/ν$ is the thermal renormalization exponent and $d$ is the spatial dimension. Similar modifications apply to other dimensionless quantities, such as Binder ratios. For systems with $2y_t-d \le 0$, these quantities share same critical universal values in the two ensembles. It is also derived that new finite-size corrections are induced. These findings apply more generally to systems in the canonical ensemble, e.g. the dilute Potts model with a fixed total number of vacancies. Finally, we formulate an efficient cluster-type algorithm for the canonical ensemble, and confirm these predictions by extensive simulations.

preprint2014arXiv

Bond and Site Percolation in Three Dimensions

We simulate the bond and site percolation models on a simple-cubic lattice with linear sizes up to L=512, and estimate the percolation thresholds to be $p_c ({\rm bond})=0.248\,811\,82(10)$ and $p_c ({\rm site})=0.311\,607\,7(2)$. By performing extensive simulations at these estimated critical points, we then estimate the critical exponents $1/ν=1.141\,0(15)$, $β/ν=0.477\,05(15)$, the leading correction exponent $y_i =-1.2(2)$, and the shortest-path exponent $d_{\rm min}=1.375\,6(3)$. Various universal amplitudes are also obtained, including wrapping probabilities, ratios associated with the cluster-size distribution, and the excess cluster number. We observe that the leading finite-size corrections in certain wrapping probabilities are governed by an exponent $\approx -2$, rather than $y_i \approx -1.2$.

preprint2014arXiv

Geometric structure of percolation clusters

We investigate the geometric properties of percolation clusters, by studying square-lattice bond percolation on the torus. We show that the density of bridges and nonbridges both tend to 1/4 for large system sizes. Using Monte Carlo simulations, we study the probability that a given edge is not a bridge but has both its loop arcs in the same loop, and find that it is governed by the two-arm exponent. We then classify bridges into two types: branches and junctions. A bridge is a {\em branch} iff at least one of the two clusters produced by its deletion is a tree. Starting from a percolation configuration and deleting the branches results in a {\em leaf-free} configuration, while deleting all bridges produces a bridge-free configuration. Although branches account for $\approx 43%$ of all occupied bonds, we find that the fractal dimensions of the cluster size and hull length of leaf-free configurations are consistent with those for standard percolation configurations. By contrast, we find that the fractal dimensions of the cluster size and hull length of bridge-free configurations are respectively given by the backbone and external perimeter dimensions. We estimate the backbone fractal dimension to be $1.643\,36(10)$.

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}$.

preprint2014arXiv

Softening of Roton and Phonon Modes in a Bose-Einstein Condensate with Spin-Orbit Coupling

Roton-type excitations usually emerge from strong correlations or long-range interactions, as in superfluid helium or dipolar ultracold atoms. However, in weakly short-range interacting quantum gas, the recently synthesized spin-orbit (SO) coupling can lead to various unconventional phases of superfluidity, and give rise to an excitation spectrum of roton-maxon character. Using Bragg spectroscopy we study a SO coupled Bose-Einstein condensate of $^{87}$Rb atoms, and show that the excitation spectrum in a "magnetized" phase clearly possesses a two-branch and roton-maxon structure. As Raman coupling strength $Ω$ is decreased, a roton-mode softening is observed, as a precursor of the phase transition to a stripe phase that spontaneously breaks spatially translational symmetry. The measured roton gaps agree well with theoretical calculations. Further, we determine sound velocities both in the magnetized and the non-magnetized phase, and a phonon-mode softening is observed around the phase transition in between. The validity of the $f$-sum rule is examined.

preprint2014arXiv

Two-species hard-core bosons on the triangular lattice: A quantum Monte Carlo study

Using worm-type quantum Monte Carlo simulations, we investigate bosonic mixtures on the triangular lattice of two species of bosons, which interact via nearest-neighbour intraspecies ($V$) and onsite interspecies ($U$) repulsions. For the case of symmetric hopping amplitude ($t_A/V=t_B/V$) and $U/V=1$, we determine a rich ground-state phase diagram that contains double solid, double superfluid (2SF), supersolid (SS), solid-superfluid (Solid-SF) and counterflow supersolid (CSS) states. The SS, Solid-SF and CSS states exhibit spontaneous symmetry breaking among the three sublattices of the triangular lattice and between the two species, which leads to nonzero crystalline density wave order in each species. We furthermore show that the CSS and the SS states are present for $t_A/V \neq t_B/V$, and the latter even survives up to $t_A/V \rightarrow \infty$ or $t_B/V \rightarrow \infty$ limit. The effects induced by the variation of $U/V$ and by the imbalance of particle numbers of the two species are also explored.

preprint2014arXiv

Universal Conductivity in a Two-Dimensional Superfluid-to-Insulator Quantum Critical System

We compute the universal conductivity of the (2+1)-dimensional XY universality class, which is realized for a superfluid-to-Mott insulator quantum phase transition at constant density. Based on large-scale Monte Carlo simulations of the classical (2+1)-dimensional $J$-current model and the two-dimensional Bose-Hubbard model, we can precisely determine the conductivity on the quantum critical plateau, $σ(\infty)=0.359(4)σ_Q$ with $σ_Q$ the conductivity quantum. The universal conductivity curve is the standard example with the lowest number of components where the bottoms-up AdS/CFT correspondence from string theory can be tested and made to use [R. C. Myers, S. Sachdev, and A. Singh, Phys. Rev. D83, 066017 (2011)]. For the first time, the shape of the $σ(iω_n)- σ(\infty)$ function in the Matsubara representation is accurate enough for a conclusive comparison and establishes the particle-like nature of charge transport. We find that the holographic gauge/gravity duality theory for transport properties can be made compatible with the data if temperature of the horizon of the black brane is different from the temperature of the conformal field theory. The requirements for measuring the universal conductivity in a cold gas experiment are also determined by our calculation.

preprint2013arXiv

Deconfined criticality flow in the Heisenberg model with ring-exchange interactions

Quantum transition points in the J -Q model---the test bed of the deconfined critical point theory---and the SU(2)-symmetric discrete noncompact CP^1 representation of the deconfined critical action are directly compared by the flowgram method. We find that the flows of two systems coincide in a broad region of linear system sizes (10 < L < 50 for the J -Q model), implying that the deconfined critical point theory correctly captures the mesoscopic physics of competition between the antiferromagnetic and valence-bond orders in quantum spin systems. At larger sizes, however, we observe significant deviations between the two flows which both demonstrate strong violations of scale invariance. This reliably rules out the second-order transition scenario in at least one of the two models and suggests the most likely explanation for the nature of the transition in the J-Q model.

preprint2013arXiv

High-precision Monte Carlo study of directed percolation in (d+1) dimensions

We present a Monte Carlo study of the bond and site directed (oriented) percolation models in $(d+1)$ dimensions on simple-cubic and body-centered-cubic lattices, with $2 \leq d \leq 7$. A dimensionless ratio is defined, and an analysis of its finite-size scaling produces improved estimates of percolation thresholds. We also report improved estimates for the standard critical exponents. In addition, we study the probability distributions of the number of wet sites and radius of gyration, for $1 \leq d \leq 7$.

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.

preprint2013arXiv

Percolation in the canonical ensemble

We study the bond percolation problem under the constraint that the total number of occupied bonds is fixed, so that the canonical ensemble applies. We show via an analytical approach that at criticality, the constraint can induce new finite-size corrections with exponent y_{can}=2y_t-d both in energy-like and magnetic quantities, where y_t=1/ν is the thermal renormalization exponent and d is the spatial dimension. Furthermore, we find that while most of universal parameters remain unchanged, some universal amplitudes, like the excess cluster number, can be modified and become non-universal. We confirm these predictions by extensive Monte Carlo simulations of the two-dimensional percolation problem which has y_{can}=-1/2.

preprint2013arXiv

Scattering and effective interactions of ultracold atoms with spin-orbit coupling

We derive an analytical expression for the scattering amplitude of two ultracold atoms of arbitrary spin and with general spin-orbit (SO) coupling, on the basis of our recent work (Phys. Rev. A \textbf{86}, 053608 (2012)). As an application, we demonstrate that SO coupling can induce scattering resonance in the case with finite scattering length. The same approach can be applied to calculate the two-body bound state of SO-coupled ultracold atoms. For the ultracold spin-1/2 fermi gases in three- or two- dimensional systems with SO coupling, we also obtain the renormalization relation of effective contact interaction with momentum cutoff, as well as the applicability of Huang-Yang pseudo-potential.

preprint2013arXiv

Simultaneous analysis of three-dimensional percolation models

We simulate the bond and site percolation models on several three-dimensional lattices, including the diamond, body-centered cubic, and face-centered cubic lattices. As on the simple-cubic lattice [Phys. Rev. E, \textbf{87} 052107 (2013)], it is observed that in comparison with dimensionless ratios based on cluster-size distribution, certain wrapping probabilities exhibit weaker finite-size corrections and are more sensitive to the deviation from percolation threshold $p_c$, and thus provide a powerful means for determining $p_c$. We analyze the numerical data of the wrapping probabilities simultaneously such that universal parameters are shared by the aforementioned models, and thus significantly improved estimates of $p_c$ are obtained.

preprint2013arXiv

Stability of Excited Dressed States with Spin-Orbit Coupling

We study the decay behaviors of ultracold atoms in metastable states with spin-orbit coupling (SOC), and demonstrate that there are two SOC-induced decay mechanisms. One arises from the trapping potential and the other is due to interatomic collision. We present general schemes for calculating decay rates from these two mechanisms, and illustrate how the decay rates can be controlled by experimental parameters.We experimentally measure the decay rates over a broad parameter region, and the results agree well with theoretical calculations. This work provides an insight for both quantum simulation involving metastable dressed states and studies on few-body problems with SO coupling.

preprint2013arXiv

Two-dimensional Potts antiferromagnets with a phase transition at arbitrarily large q

We exhibit infinite families of two-dimensional lattices (some of which are triangulations or quadrangulations of the plane) on which the q-state Potts antiferromagnet has a finite-temperature phase transition at arbitrarily large values of q. This unexpected result is proven rigorously by using a Peierls argument to measure the entropic advantage of sublattice long-range order. Additional numerical data are obtained using transfer matrices, Monte Carlo simulation, and a high-precision graph-theoretic method.

preprint2013arXiv

Universal Properties of the Higgs Resonance in (2+1)-Dimensional U(1) Critical Systems

We present spectral functions for the magnitude squared of the order parameter in the scaling limit of the two-dimensional superfluid to Mott insulator quantum phase transition at constant density, which has emergent particle-hole symmetry and Lorentz invariance. The universal functions for the superfluid, Mott insulator, and normal liquid phases reveal a low-frequency resonance which is relatively sharp and is followed by a damped oscillation (in the first two phases only) before saturating to the quantum critical plateau. The counterintuitive resonance feature in the insulating and normal phases calls for deeper understanding of collective modes in the strongly coupled (2+1)-dimensional relativistic field theory. Our results are derived from analytically continued correlation functions obtained from path-integral Monte Carlo simulations of the Bose-Hubbard model.

preprint2012arXiv

Collective Dipole Oscillation of a Spin-Orbit Coupled Bose-Einstein Condensate

We present an experimental study of the collective dipole oscillation of a spin-orbit coupled Bose-Einstein condensate in a harmonic trap. Dynamics of the center-of-mass dipole oscillation is studied in a broad parameter region, as a function of spin-orbit coupling parameters as well as oscillation amplitude. Anharmonic properties beyond effective-mass approximation are revealed, such as amplitude-dependent frequency and finite oscillation frequency at place with divergent effective mass. These anharmonic behaviors agree quantitatively with variational wave-function calculations. Moreover, we experimentally demonstrate a unique feature of spin-orbit coupled system predicted by a sum-rule approach, stating that spin polarization susceptibility--a static physical quantity--can be measured via dynamics of dipole oscillation. The divergence of polarization susceptibility is observed at the quantum phase transition that separates magnetic nonzero-momentum condensate from nonmagnetic zero-momentum phase. The good agreement between the experimental and theoretical results provides a bench mark for recently developed theoretical approaches.

preprint2012arXiv

Computational Study of a Multistep Height Model

An equilibrium random surface multistep height model proposed in [Abraham and Newman, EPL, 86, 16002 (2009)] is studied using a variant of the worm algorithm. In one limit, the model reduces to the two-dimensional Ising model in the height representation. When the Ising model constraint of single height steps is relaxed, the critical temperature and critical exponents are continuously varying functions of the parameter controlling height steps larger than one. Numerical estimates of the critical exponents can be mapped via a single parameter-- the Coulomb gas coupling-- to the exponents of the O(n) loop model on the honeycomb lattice with n <= 1.

preprint2012arXiv

Crossover from Isotropic to Directed Percolation

We generalize the directed percolation (DP) model by relaxing the strict directionality of DP such that propagation can occur in either direction but with anisotropic probabilities. We denote the probabilities as $p_{\downarrow}= p \cdot p_d$ and $p_{\uparrow}=p \cdot (1-p_d)$, with $p $ representing the average occupation probability and $p_d$ controlling the anisotropy. The Leath-Alexandrowicz method is used to grow a cluster from an active seed site. We call this model with two main growth directions {\em biased directed percolation} (BDP). Standard isotropic percolation (IP) and DP are the two limiting cases of the BDP model, corresponding to $p_d=1/2$ and $p_d=0,1$ respectively. In this work, besides IP and DP, we also consider the $1/2<p_d<1$ region. Extensive Monte Carlo simulations are carried out on the square and the simple-cubic lattices, and the numerical data are analyzed by finite-size scaling. We locate the percolation thresholds of the BDP model for $p_d=0.6$ and 0.8, and determine various critical exponents. These exponents are found to be consistent with those for standard DP. We also determine the renormalization exponent associated with the asymmetric perturbation due to $p_d -1/2 \neq 0$ near IP, and confirm that such an asymmetric scaling field is relevant at IP.

preprint2012arXiv

Experimental demonstration of topological error correction

Topological error correction--a novel method to actively correct errors based on cluster states with topological properties--has the highest order of tolerable error rates known to date (10^{-2}). Moreover, the scheme requires only nearest-neighbour interaction, particularly suitable for most physical systems. Here we report the first experimental demonstration of topological error correction with an 8-qubit optical cluster state. In the experiment, it is shown that a correlation can be protected against a single error on any single qubit. In addition, when all qubits are simultaneously subjected to errors with equal probability, the effective error rate is significantly reduced, clearly verifying the advantage of topological error correction. The quantum gate with the error rate below the threshold is within the current experimental technology. We believe topological error correction should be a critical ingredient for the future large-scale quantum computation.

preprint2012arXiv

Modified Bethe-Peierls boundary condition for ultracold atoms with Spin-Orbit coupling

We show that the Bethe-Peierls (BP) boundary condition should be modified for ultracold atoms with spin-orbit (SO) coupling. Moreover, we derive a general form of the modified BP boundary condition, which is applicable to a system with arbitrary kind of SO coupling. In the modified BP condition, an anisotropic term appears and the inter-atomic scattering length is normally SO-coupling dependent. For the special system in the current experiments, however, it can be proved that the scattering length is SO-coupling independent, and takes the same value as in the case without SO coupling. Our result is helpful for the study of both few-body and many-body physics in SO-coupled ultracold gases.

preprint2012arXiv

Novel phase transitions in XY Antiferromagnets on Plane Triangulations

Using Monte Carlo simulations and finite-size scaling, we investigate the XY antiferromagnet on the triangular, Union Jack and bisected-hexagonal lattices, and in each case find both Ising and Kosterlitz-Thouless transitions. As is well-known, on the triangular lattice, as the temperature decreases the system develops chiral order for temperatures $T < \Tc$, and then quasi-long-range magnetic order on its sublattices when $T < \Ts$, with $\Ts < \Tc$. The behavior $\Ts<\Tc$ is predicted by theoretical arguments due to Korshunov, based on the unbinding of kink-antikink pairs. On the Union Jack and bisected-hexagonal lattices, by contrast, we find that as $T$ decreases the magnetizations on some of the sublattices become quasi-long-range ordered at a temperature $\Ts > \Tc$, before chiral order develops. In some cases, the sublattice spins then undergo a second transition, of Ising type, separating two quasi-long-range ordered phases. On the Union Jack lattice, the magnetization on the degree-4 sublattice remains disordered until $\Tc$ and then undergoes an Ising transition to a quasi-long-range ordered phase.

preprint2012arXiv

Phase diagram of the toric code model in a parallel magnetic field

Ground-state phase diagram of the toric code model in a parallel magnetic field has three distinct phases: topological, charge-condensed, and vortex-condensed states. To study it we consider an implicit local order parameter characterizing the transition between the topological and charge-condensed phases, and sample it using continuous-time Monte Carlo simulations. The corresponding second-order transition line is obtained by finite-size scaling analysis of this order parameter. Symmetry breaking between charges and vortices along the first-order transition line is also observed, and our numerical result shows that the end point of the first-order transition line is located at $h_x^{(c)}=h_z^{(c)}=0.418(2)$.

preprint2012arXiv

Shortest-Path Fractal Dimension for Percolation in Two and Three Dimensions

We carry out a high-precision Monte Carlo study of the shortest-path fractal dimension $\dm$ for percolation in two and three dimensions, using the Leath-Alexandrowicz method which grows a cluster from an active seed site. A variety of quantities are sampled as a function of the chemical distance, including the number of activated sites, a measure of the radius, and the survival probability. By finite-size scaling, we determine $\dm = 1.130 77(2)$ and $1.375 6(6)$ in two and three dimensions, respectively. The result in 2D rules out the recently conjectured value $\dm=217/192$ [Phys. Rev. E 81, 020102(R) (2010)].

preprint2012arXiv

The Hintermann-Merlini-Baxter-Wu and the Infinite-Coupling-Limit Ashkin-Teller Models

We show how the Hintermann-Merlini-Baxter-Wu model (which is a generalization of the well-known Baxter-Wu model to a general Eulerian triangulation) can be mapped onto a particular infinite-coupling-limit of the Ashkin-Teller model. We work out some mappings among these models, also including the standard and mixed Ashkin-Teller models. Finally, we compute the phase diagram of the infinite-coupling-limit Ashkin-Teller model on the square, triangular, hexagonal, and kagome lattices.

preprint2012arXiv

The O(n) loop model on a three-dimensional lattice

We study a class of loop models, parameterized by a continuously varying loop fugacity n, on the hydrogen-peroxide lattice, which is a three-dimensional cubic lattice of coordination number 3. For integer n > 0, these loop models provide graphical representations for n-vector models on the same lattice, while for n = 0 they reduce to the self-avoiding walk problem. We use worm algorithms to perform Monte Carlo studies of the loop model for n = 0, 0.5, 1, 1.5, 2, 3, 4, 5 and 10 and obtain the critical points and a number of critical exponents, including the thermal exponent yt, magnetic exponent yh, and loop exponent yl. For integer n, the estimated values of yt and yh are found to agree with existing estimates for the three-dimensional O(n) universality class. The efficiency of the worm algorithms is reflected by the small value of the dynamic exponent z, determined from our analysis of the integrated autocorrelation times.

preprint2011arXiv

Deterministic spin-wave interferometer based on Rydberg blockade

The spin-wave (SW) NOON state is an $N$-particle Fock state with two atomic spin-wave modes maximally entangled. Attributed to the property that the phase is sensitive to collective atomic motion, the SW NOON state can be utilized as a novel atomic interferometer and has promising application in quantum enhanced measurement. In this paper we propose an efficient protocol to deterministically produce the atomic SW NOON state by employing Rydberg blockade. Possible errors in practical manipulations are analyzed. A feasible experimental scheme is suggested. Our scheme is far more efficient than the recent experimentally demonstrated one, which only creates a heralded second-order SW NOON state.

preprint2011arXiv

Finite-temperature phase transition in a class of 4-state Potts antiferromagnets

We argue that the 4-state Potts antiferromagnet has a finite-temperature phase transition on any Eulerian plane triangulation in which one sublattice consists of vertices of degree 4. We furthermore predict the universality class of this transition. We then present transfer-matrix and Monte Carlo data confirming these predictions for the cases of the union-jack and bisected hexagonal lattices.

preprint2011arXiv

Quench dynamics of topological quantum phase transition in Wen-plaquette model

We study the quench dynamics of the topological quantum phase transition in the two-dimensional transverse Wen-plaquette model, which has a phase transition from a Z2 topologically ordered to a spin-polarized state. By mapping the Wen-plaquette model onto a one-dimensional quantum Ising model, we calculate the expectation value of the plaquette operator Fi during a slowly quenching process from a topologically ordered state. A logarithmic scaling law of quench dynamics near the quantum phase transition is found, which is analogous to the well-known static critical behavior of the specific heat in the one-dimensional quantum Ising model.

preprint2011arXiv

Universal critical properties of the Eulerian bond-cubic model

We investigate the Eulerian bond-cubic model on the square lattice by means of Monte Carlo simulations, using an efficient cluster algorithm and a finite-size scaling analysis. The critical points and four critical exponents of the model are determined for several values of $n$. Two of the exponents are fractal dimensions, which are obtained numerically for the first time. Our results are consistent with the Coulomb gas predictions for the critical O($n$) branch for $n < 2$ and the results obtained by previous transfer matrix calculations. For $n=2$, we find that the thermal exponent, the magnetic exponent and the fractal dimension of the largest critical Eulerian bond component are different from those of the critical O(2) loop model. These results confirm that the cubic anisotropy is marginal at $n=2$ but irrelevant for $n<2$.

preprint2011arXiv

Worm Monte Carlo study of the honeycomb-lattice loop model

We present a Markov-chain Monte Carlo algorithm of "worm"type that correctly simulates the O(n) loop model on any (finite and connected) bipartite cubic graph, for any real n>0, and any edge weight, including the fully-packed limit of infinite edge weight. Furthermore, we prove rigorously that the algorithm is ergodic and has the correct stationary distribution. We emphasize that by using known exact mappings when n=2, this algorithm can be used to simulate a number of zero-temperature Potts antiferromagnets for which the Wang-Swendsen-Kotecky cluster algorithm is non-ergodic, including the 3-state model on the kagome-lattice and the 4-state model on the triangular-lattice. We then use this worm algorithm to perform a systematic study of the honeycomb-lattice loop model as a function of n<2, on the critical line and in the densely-packed and fully-packed phases. By comparing our numerical results with Coulomb gas theory, we identify the exact scaling exponents governing some fundamental geometric and dynamic observables. In particular, we show that for all n<2, the scaling of a certain return time in the worm dynamics is governed by the magnetic dimension of the loop model, thus providing a concrete dynamical interpretation of this exponent. The case n>2 is also considered, and we confirm the existence of a phase transition in the 3-state Potts universality class that was recently observed via numerical transfer matrix calculations.

preprint2011arXiv

Worm-type Monte Carlo simulation of the Ashkin-Teller model on the triangular lattice

We investigate the symmetric Ashkin-Teller (AT) model on the triangular lattice in the antiferromagnetic two-spin coupling region ($J<0$). In the $J \rightarrow -\infty$ limit, we map the AT model onto a fully-packed loop-dimer model on the honeycomb lattice. On the basis of this exact transformation and the low-temperature expansion, we formulate a variant of worm-type algorithms for the AT model, which significantly suppress the critical slowing-down. We analyze the Monte Carlo data by finite-size scaling, and locate a line of critical points of the Ising universality class in the region $J<0$ and $K>0$, with K the four-spin interaction. Further, we find that, in the $J \rightarrow -\infty$ limit, the critical line terminates at the decoupled point $K=0$. From the numerical results and the exact mapping, we conjecture that this `tricritical' point ($J \rightarrow -\infty, K=0$) is Berezinsky-Kosterlitz-Thouless-like and the logarithmic correction is absent. The dynamic critical exponent of the worm algorithm is estimated as $z=0.28(1)$ near $(J \rightarrow -\infty, K=0)$.

preprint2010arXiv

Interference of overlap-free entangled photons with a Mach-Zehnder-like interferometer

Using spontaneous parametric down conversion and a 50:50 beam splitter, we generate coaxial polarization-entangled photon pairs, of which the two photons are far separated from each other. The photons are then sent one by one through one port of a modified Mach-Zehnder interferometer. We observe interference fringes with a periodicity half of the single-photon wavelength, independent of the distance between the photons. This feature can find applications in quantum-enhanced measurement.

preprint2010arXiv

Light pulse in $Λ$-type cold atomic gases

We investigate the behavior of the light pulse in $Lambda$-type cold atomic gases with two counterpropagating control lights with equal strength by directly simulating the dynamic equations and exploring the dispersion relation. Our analysis shows that, depending on the length $L_0$ of the stored wave packet and the decay rate $γ$ of ground-spin coherence, the recreated light can behave differently. For long $L_0$ and/or large $γ$, a stationary light pulse is produced, while two propagating light pulses appear for short $L_0$ and/or small $γ$. In the $γ\to 0$ limit, the light always splits into two propagating pulses for sufficiently long time. This scenario agrees with a recent experiment [Y.-W.Lin, et al., Phys. Rev. Lett. 102, 213601(2009)] where two propagating light pulses are generated in laser-cooled cold atomic ensembles.

preprint2010arXiv

Quantum interface between frequency-uncorrelated down-converted entanglement and atomic-ensemble quantum memory

Photonic entanglement source and quantum memory are two basic building blocks of linear-optical quantum computation and long-distance quantum communication. In the past decades, intensive researches have been carried out, and remarkable progress, particularly based on the spontaneous parametric down-converted (SPDC) entanglement source and atomic ensembles, has been achieved. Currently, an important task towards scalable quantum information processing (QIP) is to efficiently write and read entanglement generated from a SPDC source into and out of an atomic quantum memory. Here we report the first experimental realization of a quantum interface by building a 5 MHz frequency-uncorrelated SPDC source and reversibly mapping the generated entangled photons into and out of a remote optically thick cold atomic memory using electromagnetically induced transparency. The frequency correlation between the entangled photons is almost fully eliminated with a suitable pump pulse. The storage of a triggered single photon with arbitrary polarization is shown to reach an average fidelity of 92% for 200 ns storage time. Moreover, polarization-entangled photon pairs are prepared, and one of photons is stored in the atomic memory while the other keeps flying. The CHSH Bell's inequality is measured and violation is clearly observed for storage time up to 1 microsecond. This demonstrates the entanglement is stored and survives during the storage. Our work establishes a crucial element to implement scalable all-optical QIP, and thus presents a substantial progress in quantum information science.

preprint2010arXiv

Some geometric critical exponents for percolation and the random-cluster model

We introduce several infinite families of new critical exponents for the random-cluster model and present scaling arguments relating them to the k-arm exponents. We then present Monte Carlo simulations confirming these predictions. These new exponents provide a convenient way to determine k-arm exponents from Monte Carlo simulations. An understanding of these exponents also leads to a radically improved implementation of the Sweeny Monte Carlo algorithm. In addition, our Monte Carlo data allow us to conjecture an exact expression for the shortest-path fractal dimension d_min in two dimensions: d_min = (g+2)(g+18)/(32g) where g is the Coulomb-gas coupling, related to the cluster fugacity q via q = 2 + 2 cos(gπ/2) with 2 \le g \le 4.

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.

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

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.

preprint2009arXiv

Single-cluster dynamics for the random-cluster model

We formulate a single-cluster Monte Carlo algorithm for the simulation of the random-cluster model. This algorithm is a generalization of the Wolff single-cluster method for the $q$-state Potts model to non-integer values $q>1$. Its results for static quantities are in a satisfactory agreement with those of the existing Swendsen-Wang-Chayes-Machta (SWCM) algorithm, which involves a full cluster decomposition of random-cluster configurations. We explore the critical dynamics of this algorithm for several two-dimensional Potts and random-cluster models. For integer $q$, the single-cluster algorithm can be reduced to the Wolff algorithm, for which case we find that the autocorrelation functions decay almost purely exponentially, with dynamic exponents $z_{\rm exp} =0.07 (1), 0.521 (7)$, and $1.007 (9)$ for $q=2, 3$, and 4 respectively. For non-integer $q$, the dynamical behavior of the single-cluster algorithm appears to be very dissimilar to that of the SWCM algorithm. For large critical systems, the autocorrelation function displays a range of power-law behavior as a function of time. The dynamic exponents are relatively large. We provide an explanation for this peculiar dynamic behavior.

Youjin Deng

What is connected

Connect this record

See the researcher in context

Building this map preview

79 published item(s)

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

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

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

Geometric Upper Critical Dimensions of the Ising Model

Machine learning for percolation utilizing auxiliary Ising variables

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

Multiple and Asymmetric Scalings in Explosive Percolation

Surface criticality of antiferromagnetic Potts model

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

Unwrapped two-point functions on high-dimensional tori

Nested Closed Paths in Two-Dimensional Percolation

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

Two-scale scenario of rigidity percolation of sticky particles

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

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

Critical polynomials in the nonplanar and continuum percolation models

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

History-dependent percolation on multiplex networks

Oscillatory-like Expansion of a Fermionic Superfluid

Scaling features in the spreading of COVID-19

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

The elastic and directed percolation backbone

Irreversible Markov chain Monte Carlo algorithm for self-avoiding walk

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

Equivalent-neighbor Potts models in two dimensions

No-enclave percolation corresponds to holes in the cluster backbone

Spin-Ice State of the Quantum Heisenberg Antiferromagnet on the Pyrochlore Lattice

Trapping Centers at the Superfluid-Mott-insulator Criticality: Transition between Charge-quantized States

Adaptive Multi-GPU Exchange Monte Carlo for the 3D Random Field Ising Model

Clique percolation in random graphs

Coulomb liquid phases of bosonic cluster Mott insulators on a pyrochlore lattice

Emergent BCS regime of the two-dimensional fermionic Hubbard model: ground-state phase diagram

Ground state phase diagram of the repulsive fermionic $t-t^{\prime}$ Hubbard model on the square lattice from weak-coupling

Leaf-excluded percolation in two and three dimensions

Realization of Two-Dimensional Spin-orbit Coupling for Bose-Einstein Condensates

Recursive Percolation

Short-range correlations in percolation at criticality

Solutions of the Two Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms

The Massive Goldstone (Higgs) mode in two-dimensional ultracold atomic lattice systems

Universal Critical Wrapping Probabilities in the Canonical Ensemble

Bond and Site Percolation in Three Dimensions

Geometric structure of percolation clusters

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

Softening of Roton and Phonon Modes in a Bose-Einstein Condensate with Spin-Orbit Coupling

Two-species hard-core bosons on the triangular lattice: A quantum Monte Carlo study

Universal Conductivity in a Two-Dimensional Superfluid-to-Insulator Quantum Critical System

Deconfined criticality flow in the Heisenberg model with ring-exchange interactions

High-precision Monte Carlo study of directed percolation in (d+1) dimensions

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

Percolation in the canonical ensemble

Scattering and effective interactions of ultracold atoms with spin-orbit coupling

Simultaneous analysis of three-dimensional percolation models

Stability of Excited Dressed States with Spin-Orbit Coupling

Two-dimensional Potts antiferromagnets with a phase transition at arbitrarily large q

Universal Properties of the Higgs Resonance in (2+1)-Dimensional U(1) Critical Systems

Collective Dipole Oscillation of a Spin-Orbit Coupled Bose-Einstein Condensate

Computational Study of a Multistep Height Model

Crossover from Isotropic to Directed Percolation

Experimental demonstration of topological error correction

Modified Bethe-Peierls boundary condition for ultracold atoms with Spin-Orbit coupling

Novel phase transitions in XY Antiferromagnets on Plane Triangulations

Phase diagram of the toric code model in a parallel magnetic field

Shortest-Path Fractal Dimension for Percolation in Two and Three Dimensions

The Hintermann-Merlini-Baxter-Wu and the Infinite-Coupling-Limit Ashkin-Teller Models

The O(n) loop model on a three-dimensional lattice

Deterministic spin-wave interferometer based on Rydberg blockade

Finite-temperature phase transition in a class of 4-state Potts antiferromagnets

Quench dynamics of topological quantum phase transition in Wen-plaquette model

Universal critical properties of the Eulerian bond-cubic model

Worm Monte Carlo study of the honeycomb-lattice loop model

Worm-type Monte Carlo simulation of the Ashkin-Teller model on the triangular lattice

Interference of overlap-free entangled photons with a Mach-Zehnder-like interferometer

Light pulse in $Λ$-type cold atomic gases

Quantum interface between frequency-uncorrelated down-converted entanglement and atomic-ensemble quantum memory