Quantum Material’s Hidden Switch Finally Understood after Decades of Debate

Researchers are actively investigating the quantum critical transition between a semimetal and a Mott insulator within the honeycomb Hubbard model, a system governed by the Gross-Neveu-Heisenberg universality class. Fo-Hong Wang, Fanjie Sun, and Chenghao He, alongside Xiao Yan Xu from the Tsung-Dao Lee Institute and Shanghai Jiao Tong University, present new findings that address long-standing discrepancies in determining the precise critical exponents of this transition. Their work is significant because it employs projector determinant Monte Carlo simulations on lattices of unprecedented size, reaching 10,368 sites, and a novel algorithmic approach to achieve high precision and access the thermodynamic limit. By establishing a robust finite-size scaling workflow and validating their methodology with the spinless Fermi-Hubbard model, the team delivers state-of-the-art critical exponents and offers a pathway towards resolving other challenging fermionic critical phenomena.

Critical exponents define the semimetal-to-Mott-insulator transition in the honeycomb Hubbard model and characterize its universality class

Scientists have achieved a breakthrough in understanding quantum materials by resolving a long-standing controversy surrounding the precise critical exponents governing the semimetal-to-Mott-insulator transition in the honeycomb Hubbard model. This transition, belonging to the Gross-Neveu-Heisenberg (GNH) universality class, is fundamental to understanding how massless Dirac fermions, relativistic quasiparticles found in materials like graphene, behave under strong interactions.
Despite its importance, determining the precise values of the critical exponents has remained elusive due to challenges posed by finite-size effects in numerical simulations and a lack of robust benchmarks. Researchers have now employed projector determinant quantum Monte Carlo (QMC) simulations on lattices of unprecedented size, reaching 10,368 sites, to address this issue.

This work introduces a novel projected submatrix update algorithm that significantly accelerates computations, enabling access to the thermodynamic limit with high precision. By performing these large-scale simulations, the fermion anomalous dimension and the correlation length exponent were observed to converge rapidly, while the boson anomalous dimension exhibited a systematic size dependence.

This dependence was resolved through careful linear extrapolation, providing definitive values for all three critical exponents. To validate their methodology, the team also conducted parallel large-scale simulations of the spinless t-V model on the honeycomb lattice, which belongs to the Gross-Neveu-Ising class.

Notably, the results for the t-V model, including the first QMC determination of the fermion anomalous dimension, align with predictions from conformal bootstrap theory, confirming the robustness of the approach. This research delivers state-of-the-art critical exponents for the honeycomb Hubbard model and establishes a systematic finite-size scaling workflow applicable to a broad range of strongly correlated quantum systems. The findings pave the way for resolving other challenging fermionic quantum critical phenomena and offer insights into the behaviour of materials with strong electron correlations.

Computational methods for correlated fermions on the honeycomb lattice are rapidly advancing

Projector determinant quantum Monte Carlo simulations formed the core of this research, performed on lattices reaching an unprecedented 10,368 sites. A novel projected submatrix update algorithm was developed to significantly accelerate computations, enabling access to the thermodynamic limit with enhanced precision.

This algorithmic speedup was crucial for mitigating severe finite-size effects previously hindering accurate determination of critical exponents in strongly correlated systems. The study investigated the interplay between Dirac fermions and electronic correlations on the honeycomb lattice, focusing on the semimetal to Mott insulator transition governed by the Gross-Neveu-Heisenberg universality class.

To achieve high statistical accuracy, simulations were conducted using large system sizes and long thermalisation times, carefully monitoring convergence of key observables. The fermion anomalous dimension and correlation length exponent were observed to converge rapidly with increasing lattice size, indicating a well-defined critical regime.

A systematic size dependence was identified in the boson anomalous dimension, which was then resolved through linear extrapolation to the thermodynamic limit. To validate the methodology and ensure robustness, parallel large-scale simulations of the spinless t-V model on the honeycomb lattice, belonging to the Gross-Neveu-Ising universality class, were also performed.

Results for this model, including the first quantum Monte Carlo determination of the fermion anomalous dimension, demonstrated agreement with conformal bootstrap predictions, corroborating the reliability of the computational approach. This work establishes a systematic finite-size scaling workflow applicable to a broad range of strongly correlated fermionic systems, facilitating the resolution of other challenging quantum critical phenomena.

Fermion and boson anomalous dimensions in honeycomb lattice models via large-scale quantum Monte Carlo simulations reveal universal behavior

Researchers achieved a fermion anomalous dimension of 0.235 and a correlation length exponent of 0.805 through projector determinant Monte Carlo simulations on lattices reaching 10,368 sites. These calculations were performed on the honeycomb lattice to investigate the semimetal to Mott insulator transition governed by the Gross-Neveu-Heisenberg universality class.

A novel projected submatrix update algorithm was developed, significantly accelerating computations and enabling access to the thermodynamic limit with increased precision. The boson anomalous dimension exhibited systematic size dependence, which was resolved via linear extrapolation to obtain a final value.

Simulations were also conducted on the spinless t-V model, belonging to the Gross-Neveu-Ising class, to validate the methodology. Results for this model, including the first quantum Monte Carlo determination of the fermion anomalous dimension, demonstrated agreement with conformal bootstrap predictions, confirming the robustness of the approach.

The study established a systematic finite-size scaling workflow applicable to a broad range of strongly correlated systems. Measurements of the fermion anomalous dimension and correlation length exponent converged rapidly, indicating reliable results even with moderate system sizes. The achieved precision in determining critical exponents for the honeycomb Hubbard model provides a benchmark for future investigations of challenging fermionic critical phenomena. This work represents a substantial advancement in understanding quantum phase transitions in strongly correlated materials.

Critical exponent values established for semimetal to Mott insulator transition are consistent with theoretical predictions

Researchers have achieved high-precision determination of critical exponents for a fundamental phase transition occurring in a honeycomb lattice system. This transition involves the shift from a semimetal to a Mott insulator, a phenomenon governed by the Gross-Neveu-Heisenberg universality class, and has long been subject to debate regarding the precise values of its critical exponents.

By employing projector determinant Monte Carlo simulations on lattices containing up to 10,368 sites, a significantly larger scale than previously possible, they have obtained robust estimates for these exponents. The study utilized a novel algorithm to accelerate computations, enabling access to the thermodynamic limit with improved accuracy.

Analysis revealed distinct convergence patterns for the critical exponents, allowing for their definitive extraction through systematic finite-size scaling. Validation was performed through parallel large-scale simulations of a related spinless model, confirming the methodology and yielding the first quantum Monte Carlo determination of the fermion anomalous dimension for that system.

These results align with predictions from conformal bootstrap theory, bolstering confidence in the approach. The authors acknowledge that finite-size effects remain a challenge in numerical studies of this type, and their analysis involved careful extrapolation to mitigate these limitations. Future research may focus on applying this systematic finite-size scaling workflow to other strongly correlated fermionic systems, potentially resolving long-standing questions in that field. This work establishes a benchmark for critical exponents in the honeycomb Hubbard model and provides a valuable tool for investigating challenging quantum critical phenomena.