The three-dimensional (3D) Ising model is a foundational model in statistical physics and critical phenomena, yet its analytical intractability has long impeded the precise determination of universal critical exponents. While high-precision estimates have been obtained through classical numerical methods and conformal bootstrap techniques, a direct quantum simulation of the 3D Ising criticality remains challenging, requiring nontrivial connectivity, sufficient system size, and high spectral resolution. In this work, assisted by the state-operator correspondence of conformal field theory, we perform a digital quantum simulation of the 3D Ising critical exponents using a multiply-connected 9-qubit superconducting quantum processor with a Platonic lattice geometry. Employing an extended variational quantum eigensolver equipped with a phase-based loss function, we variationally prepare the low-energy eigenstates of the transverse-field Ising model on a cubic Platonic lattice encoded in an 8-qubit register. The four lowest eigenenergies are extracted via Fourier-transform analysis and high-precision numerical fitting, agreeing with the exact diagonalization values up to +/- 0.001. The resulting scaling dimension Delta_epsilon = 1.5850 and critical exponent nu = 0.7067 match well with theory.
Comments
19 pages, 11 figures, 5 tables; supplementary information included
The three-dimensional (3D) Ising model is a foundational model in statistical physics and critical phenomena, yet its analytical intractability has long impeded the precise determination of universal critical exponents. While high-precision estimates have been obtained through classical numerical methods and conformal bootstrap techniques, a direct quantum simulation of the 3D Ising criticality remains challenging, requiring nontrivial connectivity, sufficient system size, and high spectral resolution. In this work, assisted by the state-operator correspondence of conformal field theory, we perform a digital quantum simulation of the 3D Ising critical exponents using a multiply-connected 9-qubit superconducting quantum processor with a Platonic lattice geometry. Employing an extended variational quantum eigensolver equipped with a phase-based loss function, we variationally prepare the low-energy eigenstates of the transverse-field Ising model on a cubic Platonic lattice encoded in an 8-qubit register. The four lowest eigenenergies are extracted via Fourier-transform analysis and high-precision numerical fitting, agreeing with the exact diagonalization values up to +/- 0.001. The resulting scaling dimension Delta_epsilon = 1.5850 and critical exponent nu = 0.7067 match well with theory.
The critical dynamics of `model A'' of Hohenberg and Halperin has been studied by the Monte Carlo method. Simulations have been carried out in the three-dimensional (3d) simple cubic Ising model for lattices of sizes L=16 to L=512. Using Wolff's cluster algorithm, the critical temperature is precisely found as ${\ensuremath{\beta}}_{c}=0.221\phantom{\rule{0.16em}{0ex}}654\phantom{\rule{0.16em}{0ex}}68(5)$. By Fourier transform of the lattice configurations, the critical scattering intensities $I(\stackrel{P\vec}{q}\phantom{\rule{0.16em}{0ex}})$ can be obtained. After circular averaging, the static simulations with L=256 and L=512 provide an estimate of the critical exponent $\ensuremath{\gamma}/\ensuremath{\nu}=2\ensuremath{-}\ensuremath{\eta}=1.9640(5)$. The $|\stackrel{P\vec}{q}\phantom{\rule{0.16em}{0ex}}|$-dependent distribution of $I(\stackrel{P\vec}{q}\phantom{\rule{0.16em}{0ex}})$ showed an exponential distribution, corresponding to a Gaussian distribution of the scattering amplitudes for a large q domain. The time-dependent intensities were then used for the study of the critical dynamics of 3d lattices at the critical point. To simulate results of an x-ray photon correlation spectroscopy experiment, the time-dependent correlation function of the intensities was studied for each $|\stackrel{P\vec}{q}\phantom{\rule{0.16em}{0ex}}|$-value. In the q region where $I(\stackrel{P\vec}{q}\phantom{\rule{0.16em}{0ex}})$ had an exponential distribution, the time correlations can be fit to a stretched exponential, where the exponent $\ensuremath{\mu}=\ensuremath{\gamma}/\ensuremath{\nu}z\ensuremath{\simeq}0.975$ provides an estimate of the dynamic exponent z. This corresponds to z=2.0145, in agreement with the observed variations of the characteristic fluctuation time of the intensity: $\ensuremath{\tau}(q)\ensuremath{\propto}{q}^{\ensuremath{-}z}$, which gives z=2.015. These results agree with the $\ensuremath{\epsilon}$ expansion of field-theoretical methods (2.017). In this paper, the need to take account of the anomalous time behavior ($\ensuremath{\mu}<1$) in the dynamics is exemplified. This dynamics reflects a nonlinear time behavior of model A, and its large time extension is discussed in detail.
We extend the Prometheus framework for unsupervised phase transition discovery from two-dimensional classical systems to three-dimensional classical systems and quantum many-body systems. Building upon preliminary observations from a 2D Ising model student abstract [Yee et al., 2026], we address two fundamental questions: (1) Does the framework scale to higher dimensions where exact solutions are unavailable? (2) Can it generalize to quantum phase transitions driven by quantum fluctuations rather than thermal fluctuations? For the 3D Ising model on lattices up to L=32, we achieve critical temperature detection within 0.01\% of literature values ($\Tc/J = 4.511 \pm 0.005$) and extract critical exponents with ≥70% accuracy, with statistical analysis correctly identifying the 3D Ising universality class (p=0.72). For quantum systems, we develop quantum-aware VAE (Q-VAE) architectures operating on complex-valued wavefunctions with fidelity-based loss functions, achieving 2\% accuracy in quantum critical point detection for the transverse field Ising model. For the disordered TFIM, we perform a consistency check of activated dynamical scaling $\ln ξ\sim |h - \hc|^{-ψ}$, extracting tunneling exponent ψ=0.48±0.08 consistent with theoretical predictions (ψ=0.5, Δχ2=12.3, p<0.001). This demonstrates that unsupervised learning can identify qualitatively different \emph{types} of critical behavior, serving as a consistency check on known IRFP physics.
The quantum phase transitions provide a paradigm for studying collective quantum phenomena that are a result of competing non-commuting interactions. This paper will study the ground state properties and quantum critical dynamics of the one-dimensional transverse field Ising model through a combined perspective that includes exact diagonalisation, variational quantum eigensolver (VQE) simulations, and simulations on realistic physical quantum devices. We focus on a lattice of four spins, where we calculate the ground-state energies, magnetic order parameters and correlation functions at uniformly applied conditions, which is repeated by all systems. Precise diagonalisation provides both a benchmark, which is symmetry-conserving, and a depth-two, physics inspired variational approximation, which provides simulations accessible to hardware. The circuits that have been optimised identically are then placed on the IQM Garnet quantum processor, using a resource-efficient batched protocol. We find that the ground-state energies of shallow variational circuits are reliably captured by the circuit over the entire parameter space; the magnetic arrangement parameters and observables sensitive to correlation signal significantly more noise. The error analysis of quantitative analysis reveals a strong broadening of critical crossover on hardware, which is consistent with the noise attenuation of long-range correlations. These findings highlight the current capabilities as well as the fundamental limitations of noisy intermediate-scale quantum systems in modelling quantum critical phenomena as a benchmark to future enhancements in obtaining quantum hardware and quantum algorithms development.
Quantum phase transitions driven by electronic correlations are central to understanding the physics of graphene and related two-dimensional materials. A paradigmatic example is the semimetal-to-Mott-insulator transition on the honeycomb lattice, governed by the Gross-Neveu-Heisenberg universality class, yet consensus on its critical exponents has remained elusive for over a decade due to severe finite-size effects and the absence of rigorous conformal bootstrap benchmarks. Here we try to resolve this long-standing controversy by performing projector determinant quantum Monte Carlo simulations on lattices of unprecedented size, reaching 10,368 sites. By developing a novel projected submatrix update algorithm, we achieve a significant algorithmic speedup that enables us to access the thermodynamic limit with high precision. We observe that the fermion anomalous dimension and the correlation length exponent converge rapidly, while the boson anomalous dimension exhibits a systematic size dependence that we resolve via linear extrapolation. To validate our analysis, we perform parallel large-scale simulations of the spinless t-V model on the honeycomb lattice, which belongs to the Gross-Neveu-Ising class. Our results for the t-V model show agreement with conformal bootstrap predictions, thereby corroborating the robustness of our methodology. Our work provides state-of-the-art critical exponents for the honeycomb Hubbard model and establishes a systematic finite-size scaling workflow applicable to a broad class of strongly correlated quantum systems, paving the way for resolving other challenging fermionic quantum critical phenomena.
Herdeiro and Doyon [Phys. Rev. E 94, 043322 (2016)2470-004510.1103/PhysRevE.94.043322] introduced a numerical recipe, dubbed uv sampler, offering precise estimations of the conformal field theory (CFT) data of the planar two-dimensional (2D) critical Ising model. It made use of scale invariance emerging at the critical point in order to sample finite sublattice marginals of the infinite plane Gibbs measure of the model by producing holographic boundary distributions. The main ingredient of the Markov chain Monte Carlo sampler is the invariance under dilation. This paper presents a generalization to higher dimensions with the critical 3D Ising model. This leads to numerical estimations of a subset of the CFT data-scaling weights and structure constants-through fitting of measured correlation functions. The results are shown to agree with the recent most precise estimations from numerical bootstrap methods [Kos, Poland, Simmons-Duffin, and Vichi, J. High Energy Phys. 08 (2016) 03610.1007/JHEP08(2016)036].
We perform a high-statistics simulation of the three-dimensional randomly dilute Ising model on cubic lattices L3 with L< or =256. We choose a particular value of the density, x=0.8, for which the leading scaling corrections are suppressed. We determine the critical exponents, obtaining nu=0.683(3), eta=0.035(2), beta=0.3535(17), and alpha=-0.049(9), in agreement with previous numerical simulations. We also estimate numerically the fixed-point values of the four-point zero-momentum couplings that are used in field-theoretical fixed-dimension studies. Although these results somewhat differ from those obtained using perturbative field theory, the field-theoretical estimates of the critical exponents do not change significantly if the Monte Carlo result for the fixed point is used. Finally, we determine the six-point zero-momentum couplings, relevant for the small-magnetization expansion of the equation of state, and the invariant amplitude ratio R(+)(xi) that expresses the universality of the free-energy density per correlation volume. We find R(+)(xi)=0.2885(15).