Monte Carlo Renormalization Group Research Articles

Corrections to scaling in the 3D Ising model: A comparison between MC and MCRG results

Corrections to scaling in the 3D Ising model are studied based on Monte Carlo (MC) simulation results for very large lattices with linear lattice sizes up to [Formula: see text]. Our estimated values of the correction-to-scaling exponent [Formula: see text] tend to decrease below the usually accepted value about 0.83 when the smallest lattice sizes, i.e. [Formula: see text] with [Formula: see text], are discarded from the fits. This behavior apparently confirms some of the known estimates of the Monte Carlo renormalization group (MCRG) method, i.e. [Formula: see text] and [Formula: see text]. We discuss the possibilities that [Formula: see text] is either really smaller than usually expected or these values of [Formula: see text] describe some transient behavior which, eventually, turns into the correct asymptotic behavior at [Formula: see text]. We propose refining MCRG simulations and analysis to resolve this issue. Our actual MC estimations of the critical exponents [Formula: see text] and [Formula: see text] provide stable values [Formula: see text] and [Formula: see text], which well agree with those of the conformal bootstrap method, i.e. [Formula: see text] and [Formula: see text].

First-order transition in the stacked-[formula omitted]-[formula omitted] Ising model on a cubic lattice

We investigate critical properties of the stacked-J1-J2 Ising model on a cubic lattice. Using Monte Carlo simulations and renormalization group, we find a single phase transition of the first order for J2/J1>1/2. The renormgroup approach predicts that a transition can be of the second order from the universality class of the O(2) model, but the Monte Carlo results show another set of critical exponents: exponents continuously vary form the values typical for a first-order transition in the finite-size scaling theory at 1/2<J2/J1<1 to the Ising values in the limit J2/J1→∞. We also exclude the pseudo-first-order behavior observed in the J1-J2 Ising model on a square lattice for 0.67≲J2/J1≲0.9.

Magnetization of the spin- 12 Heisenberg antiferromagnet on the triangular lattice

After decades of debate, now there is a rough consensus that at zero temperature the spin-$1/2$ Heisenberg antiferromagnet on the triangular lattice is three-sublattice ${120}^{\ensuremath{\circ}}$ magnetically ordered, in contrast to a quantum spin liquid as originally proposed. However, there remains considerable discrepancy in the magnetization reported among various methods. To resolve this issue, in this work we revisit this model by the tensor-network state algorithm. The ground-state energy per bond ${E}_{b}$ and magnetization per spin ${M}_{0}$ in the thermodynamic limit are obtained with high precision. The former is estimated to be ${E}_{b}=\ensuremath{-}0.18334(10)$. This value agrees well with that from the series expansion. The three-sublattice magnetic order is firmly confirmed and the magnetization is determined as ${M}_{0}=0.161(5)$. It is about $32%$ of its classical value and slightly below the lower bound from the series expansion. In comparison with the best estimated value by Monte Carlo and density-matrix renormalization group, our result is about $20%$ smaller. This magnetic order is consistent with further analysis of the three-body correlation. Our work thus provides benchmark results for this prototypical model.

Phase diagram of a distorted kagome antiferromagnet and application to Y-kapellasite

We investigate the magnetism of a previously unexplored distorted spin-1/2 kagome model consisting of three symmetry-inequivalent nearest-neighbor antiferromagnetic Heisenberg couplings J⬡, J, and J^{prime}, and uncover a rich ground state phase diagram even at the classical level. Using analytical arguments and numerical techniques we identify a collinear overrightarrow{Q}=0 magnetic phase, two unusual non-collinear coplanar overrightarrow{Q}=(1/3,1/3) phases and a classical spin liquid phase with a degenerate manifold of non-coplanar ground states, resembling the jammed spin liquid phase found in the context of a bond-disordered kagome antiferromagnet. We further show with density functional theory calculations that the recently synthesized Y-kapellasite Y3Cu9(OH)19Cl8 is a realization of this model and predict its ground state to lie in the region of overrightarrow{Q}=(1/3,1/3) order, which remains stable even after the inclusion of quantum fluctuation effects within variational Monte Carlo and pseudofermion functional renormalization group. The presented model opens a new direction in the study of kagome antiferromagnets.

Clock model interpolation and symmetry breaking in O(2) models

Motivated by recent attempts to quantum simulate lattice models with continuous Abelian symmetries using discrete approximations, we define an extended-O(2) model by adding a $\gamma \cos(q\varphi)$ term to the ordinary O(2) model with angular values restricted to a $2\pi$ interval. In the $\gamma \rightarrow \infty$ limit, the model becomes an extended $q$-state clock model that reduces to the ordinary $q$-state clock model when $q$ is an integer and otherwise is a continuation of the clock model for noninteger $q$. By shifting the $2\pi$ integration interval, the number of angles selected can change discontinuously and two cases need to be considered. What we call case $1$ has one more angle than what we call case $2$. We investigate this class of clock models in two space-time dimensions using Monte Carlo and tensor renormalization group methods. Both the specific heat and the magnetic susceptibility show a double-peak structure for fractional $q$. In case $1$, the small-$\beta$ peak is associated with a crossover, and the large-$\beta$ peak is associated with an Ising critical point, while both peaks are crossovers in case $2$. When $q$ is close to an integer by an amount $\Delta q$ and the system is close to the small-$\beta$ Berezinskii-Kosterlitz-Thouless transition, the system has a magnetic susceptibility that scales as $\sim 1 / (\Delta q)^{1 - 1/\delta'}$ with $\delta'$ estimates consistent with the magnetic critical exponent $\delta = 15$. The crossover peak and the Ising critical point move to Berezinskii-Kosterlitz-Thouless transition points with the same power-law scaling. A phase diagram for this model in the $(\beta, q)$ plane is sketched. These results are possibly relevant for configurable Rydberg-atom arrays where the interpolations among phases with discrete symmetries can be achieved by varying continuously the distances among atoms and the detuning frequency.

Monte Carlo renormalization-group calculation for the d=3 Ising model using a modified transformation.

We present a simple approach to high-accuracy calculations of critical properties for the three-dimensional Ising model, without prior knowledge of the critical temperature. The iterative method uses a modified block-spin transformation with a tunable parameter to improve convergence in the Monte Carlo renormalization group trajectory. We found experimentally that the iterative method enables the calculation of the critical temperature simultaneously with a critical exponent.

Constraining the parameter space of a quantum spin liquid candidate in applied field with iterative optimization

The quantum spin liquid (QSL) state is an exotic state of matter featuring a high degree of entanglement and lack of long-range magnetic order in the zero-temperature limit. The triangular antiferromagnet YbMgGaO4 is a candidate QSL host, and precise determination of the Hamiltonian parameters is critical to understanding the nature of the possible ground states. However, the presence of chemical disorder has made directly measuring these parameters challenging. Here we report neutron scattering and magnetic susceptibility measurements covering a broad range of applied magnetic field at low temperature. Our data shows a field-induced crossover in YbMgGaO4, which we reproduce with complementary classical Monte Carlo and Density Matrix Renormalization Group simulations. Neutron scattering data above and below the crossover reveal a shift in scattering intensity from M to K points and, collectively, our measurements provide essential characteristics of the phase crossover that we employ to strictly constrain proposed magnetic Hamiltonian parameters despite the chemical disorder. Constrained exchange parameters further suggest the material's proximity to the QSL state in the clean limit. More broadly, our approach demonstrates a means of pursuing QSL candidates where Hamiltonian parameters might otherwise be obscured by disorder.

Pairing and Pair Superfluid Density in One-Dimensional Two-Species Fermionic and Bosonic Hubbard Models.

We use unbiased computational methods to elucidate the onset and properties of pair superfluidity in two-species fermionic and bosonic systems with onsite interspecies attraction loaded in a uniform, i.e., with no confining potential, one-dimensional optical lattice. We compare results from quantum MonteCarlo (QMC) and density matrix renormalization group (DMRG), emphasizing the one-to-one correspondence between the Drude weight tensor, calculated with DMRG, and the various winding numbers extracted from the QMC. Our results show that, for any nonvanishing attractive interaction, pairs form and are the sole contributors to superfluidity; there are no individual contributions due to the separate species. For weak attraction, the pair size diverges exponentially, i.e., Bardeen-Cooper-Schrieffer (BCS) pairing, requiring huge systems to bring out the pair-only nature of the superfluid. This crucial property is largely overlooked in many studies, thereby misinterpreting the origin and nature of the superfluid. We compare and contrast this with the repulsive case and show that the behavior is very different, contradicting previous claims about drag superfluidity and the symmetry of properties for attractive and repulsive interactions. Finally, our results show that the situation is similar for soft-core bosons: superfluidity is due only to pairs, even for the smallest attractive interaction strength compatible with the largest system sizes that we could attain.

Accurate simulation for finite projected entangled pair states in two dimensions

Based on the scheme of variational Monte Carlo sampling, we develop an accurate and efficient two-dimensional tensor-network algorithm to simulate quantum lattice models. We find that Monte Carlo sampling shows huge advantages in dealing with finite projected entangled pair states, which allows significantly enlarged system size and improves the accuracy of tensor network simulation. We demonstrate our method on the square-lattice antiferromagnetic Heisenberg model up to $32 \times 32$ sites, as well as a highly frustrated $J_1-J_2$ model up to $24\times 24$ sites. The results, including ground state energy and spin correlations, are in excellent agreement with those of the available quantum Monte Carlo or density matrix renormalization group methods. Therefore, our method substantially advances the calculation of 2D tensor networks for finite systems, and potentially opens a new door towards resolving many challenging strongly correlated quantum many-body problems.

MonteCarlo Renormalization Group for Classical Lattice Models with Quenched Disorder.

We extend to quenched-disordered systems the variational scheme for real-space renormalization group calculations that we recently introduced for homogeneous spin Hamiltonians. When disorder is present our approach gives access to the flow of the renormalized Hamiltonian distribution, from which one can compute the critical exponents if the correlations of the renormalized couplings retain finite range. Key to the variational approach is the bias potential found by minimizing a convex functional in statistical mechanics. This potential reduces dramatically the MonteCarlo relaxation time in large disordered systems. We demonstrate the method with applications to the two-dimensional dilute Ising model, the random transverse field quantum Ising chain, and the random field Ising in two- and three-dimensional lattices.

Continuous-time Monte Carlo renormalization group

We implement Monte Carlo Renormalization Group (MCRG) in the continuous-time Monte Carlo simulation of a quantum system. We demonstrate numerically the emergent isotropy between space and time at large distances for the systems that exhibit Lorentz invariance at quantum criticality. This allows us to estimate accurately the sound velocity for these quantum systems. $Q$-state Potts models in one and two space dimensions are used to illustrate the method.

Absence of Superconductivity in the Pure Two-Dimensional Hubbard Model

We study the superconducting pairing correlations in the ground state of the doped Hubbard model -- in its original form without hopping beyond nearest neighbor or other perturbing parameters -- in two dimensions at intermediate to strong coupling and near optimal doping. The nature of such correlations has been a central question ever since the discovery of cuprate high-temperature superconductors. Despite unprecedented effort and tremendous progress in understanding the properties of this fundamental model, a definitive answer to whether the ground state is superconducting in the parameter regime most relevant to cuprates has proved exceedingly difficult to establish. In this work, we employ two complementary, state-of-the-art many-body computational methods, constrained path (CP) auxiliary-field quantum Monte Carlo (AFQMC) and density matrix renormalization group (DMRG) methods, deploying the most recent algorithmic advances in each. Systematic and detailed comparisons between the two methods are performed. The DMRG is extremely reliable on small width cylinders, where we use it to validate the AFQMC. The AFQMC is then used to study wide systems as well as fully periodic systems, to establish that we have reached the thermodynamic limit. The ground state is found to be non-superconducting in the moderate to strong coupling regime in the vicinity of optimal hole doping.

Determination of the critical manifold tangent space and curvature with Monte Carlo renormalization group.

We show that the critical manifold of a statistical mechanical system in the vicinity of a critical point is locally accessible through correlation functions at that point. A practical numerical method is presented to determine the tangent space and the curvature to the critical manifold with variational Monte Carlo renormalization group. Because of the use of a variational bias potential of the coarse-grained variables, critical slowing down is greatly alleviated in the Monte Carlo simulation. In addition, this method is free of truncation error. We study the isotropic Ising model on square and cubic lattices, the anisotropic Ising model, and the tricritical Ising model on square lattices to illustrate the method.

Tangent space to the manifold of critical classical Hamiltonians representable by tensor networks.

We introduce a scheme to perform Monte Carlo renormalization group with the coupling constants of the system Hamiltonian encoded in a tensor network. With this scheme we compute the tangent space to the manifold of the critical Hamiltonians representable by a tensor network at the nearest-neighbor critical coupling for three models: the two- and three-dimensional Ising models and the two-dimensional three-state Potts model.

Study on the two-dimensional kinetic Ising model with the dynamic Monte Carlo renormalization group method

Using a modified dynamic Monte Carlo renormalization group method, the two- dimensional kinetic Ising model is studied, and the dynamic critical exponent is obtained. The critical temperature of phase transition can be obtained by the renormalization method for the correlation function. In the method we used, the correlation function is replaced with the absolute value of the magnetization, and it is found that the evolution of the absolute value of the magnetization over time satisfies the power-law form. It is found that the value of the dynamic critical exponent tends to be a stable value in the form of a power-law function as the scale of the system increases. The dynamic critical exponent obtained is z≃2.15. When the modified dynamic Monte Carlo renormalization group method is applied to the two-dimensional Glauber model, the obtained dynamic critical exponent is z≃2.25.

6a-B-2 MCRG法による格子QCDのスケーリング解析

6a-B-2 MCRG法による格子QCDのスケーリング解析

Variational Approach to MonteCarlo Renormalization Group.

We present a MonteCarlo method for computing the renormalized coupling constants and the critical exponents within renormalization theory. The scheme, which derives from a variational principle, overcomes critical slowing down, by means of a bias potential that renders the coarse grained variables uncorrelated. The two-dimensional Ising model is used to illustrate the method.

Surprising convergence of the Monte Carlo renormalization group for the three-dimensional Ising model.

We present a surprisingly simple approach to high-accuracy calculations of the critical properties of the three-dimensional Ising model. The method uses a modified block-spin transformation with a tunable parameter to improve convergence in the Monte Carlo renormalization group. The block-spin parameter must be tuned differently for different exponents to produce optimal convergence.

Nonlocal correlations in the orbital selective Mott phase of a one-dimensional multiorbital Hubbard model

We study non-local correlations in a three-orbital Hubbard model defined on an extended one-dimensional chain using determinant quantum Monte Carlo and density matrix renormalization group methods. We focus on a parameter with robust Hund's coupling, which produces an orbital selective Mott phase (OSMP) at intermediate values of the Hubbard U, as well as an orbitally ordered ferromagnetic insulating state at stronger coupling. An examination of the orbital and spin-correlation functions indicates that the orbital ordering occurs before the onset of magnetic correlations in this parameter regime as a function of temperature. In the OSMP, we find that the self-energy for the itinerant electrons is momentum dependent, indicating a degree of non-local correlations while the localized electrons have largely momentum independent self-energies. These non-local correlations also produce relative shifts of the hole-like and electron-like bands within our model. The overall momentum dependence of these quantities is strongly suppressed in the orbitally-ordered insulating phase.

Critical properties of dissipative quantum spin systems in finite dimensions

We study the critical properties of finite-dimensional dissipative quantum spin systems with uniform ferromagnetic interactions. Starting from the transverse field Ising model coupled to a bath of harmonic oscillators with Ohmic spectral density, we generalize its classical representation to classical spin systems with O(n) symmetry and then take the large-n limit to reduce the system to a spherical model. The exact solution to the resulting spherical model with long-range interactions along the imaginary time axis shows a phase transition with static critical exponents coinciding with those of the conventional short-range spherical model in dimensions, where d is the spatial dimensionality of the original quantum system. This implies that the dynamical exponent is z = 2. These conclusions are consistent with the results of Monte Carlo simulations and renormalization group calculations for dissipative transverse field Ising and O(n) models in one and two dimensions. The present approach therefore serves as a useful tool for analytically investigating the properties of quantum phase transitions of the dissipative transverse field Ising and other related models. Our method may also offer a platform to study more complex phase transitions in dissipative finite-dimensional quantum spin systems, which have recently received renewed interest in the context of quantum annealing in a noisy environment.