Abstract
The Hofstadter–Hubbard model captures the physics of strongly correlated electrons in an applied magnetic field, which is relevant to many recent experiments on Moiré materials. Few large-scale, numerically exact simulations exists for this model. In this work, we simulate the Hubbard–Hofstadter model using the determinant quantum Monte Carlo (DQMC) algorithm. We report the field and Hubbard interaction strength dependence of charge compressibility, fermion sign, local moment, magnetic structure factor, and specific heat. The gross structure of magnetic Bloch bands and band gaps determined by the non-interacting Hofstadter spectrum is preserved in the presence of U. Incompressible regions of the phase diagram have improved fermion sign. At half filling and intermediate and larger couplings, a strong orbital magnetic field delocalizes electrons and reduces the effect of Hubbard U on thermodynamic properties of the system.
Similar content being viewed by others
Introduction
Strong magnetic fields allow us to probe the phase diagram of strongly correlated materials and uncover novel phases. For example, a magnetic field induces charge/pair density wave order in cuprate superconductors1,2,3,4; magnetic field is a convenient tuning parameter for accessing quantum critical points5,6, and field-induced reentrant superconductivity also has been reported in uranium compounds7,8. With the recent proliferation of experimental evidence for fractional quantum Hall effect, superconductivity, and other correlated electron phases9,10,11,12,13,14 in graphene Moiré superlattices, there is renewed interest in studying the behavior of strongly correlated electronic systems in strong magnetic fields.
Properties of non-interacting electrons in a two-dimensional periodic lattice under the influence of a strong magnetic field are fairly well-understood. In this system, the competition between lattice and magnetic length scales leads to the fractal Hofstadter butterfly spectrum with recursive magnetic subband structure15,16, which generalizes the idea of Landau levels in a free electron gas. The Chern numbers associated with these magnetic subbands provide an elegant explanation of the integer quantum Hall effect17. Experimentally, the Hofstadter Hamiltonian has been realized in ultra-cold atoms loaded on optical lattices18, and direct observation of the Hofstadter spectrum has been reported in Moiré superlattices in graphene with high resolution19,20.
The most natural framework for understanding the simultaneous influence of magnetic field and Coulomb interaction on electrons in a periodic lattice is to take the Hofstadter Hamiltonian and add to it a Hubbard interaction term. In the literature, this is sometimes called the Hofstadter–Hubbard or Hubbard–Hofstadter model. This model has been investigated using Hartree–Fock mean-field theory21,22, exact diagonalization23,24, dynamical mean-field theory25,26, and in the large U limit via renormalized mean-field theory27. Aside from exact diagonalization, which is limited to small system sizes, all methods used to study the Hubbard–Hofstadter model have been approximate and don’t capture the full extent of quantum fluctuations. It is not conclusive, for example, whether interactions change or preserve the gap structure of the Hofstadter butterfly22,24,26.
Determinant quantum Monte Carlo (DQMC)28,29,30 is an unbiased and numerically exact algorithm for studying quantum systems at finite temperature. It employs a discrete Hubbard–Stratonovich transformation to reduce the quartic Hubbard interaction term to quadratic at the cost of introducing a fluctuating auxiliary field. This auxiliary field is then sampled using the Metropolis–Hastings algorithm. DQMC has been employed successfully in the (zero-field) Hubbard model to study spin and charge excitations31,32 and superconducting fluctuations33, as well as find evidence for fluctuating stripes34 and T-linear resistivity35. The DQMC method is especially powerful at half-filling in the absence of kinetic frustration36, where the fermion sign problem is absent due to particle-hole symmetry, even in the presence of a magnetic field. This allows simulations to be performed at much lower temperatures, providing access to properties more reflective of the ground state.
At half filling, the Fermi surface of the non-interacting Hofstadter model consists of a finite number of Dirac points at even-denominator rational fractions of magnetic flux per plaquette37. The ground state of the Hubbard–Hofstadter model is thus expected to remain a Dirac semi-metal up to some finite coupling strength Uc. At half a magnetic flux quantum per plaquette, the model also is known as the π-flux model38, and has been studied extensively numerically. As interactions are turned on, the π-flux model exhibits a quantum phase transition of the chiral Heisenberg Gross-Neveu universality class at Uc ≈ 5.6t into an antiferromagnetic Mott insulator (AFMI)39,40,41,42,43. Since the π-flux model corresponds to the Hubbard–Hofstadter Hamiltonian threaded with maximum possible flux, we may think of the zero-field Hubbard model on a half-filled square lattice as the “0-flux model”, which exibits a metal–AFMI transition with Uc = 029,30. Our simulations address intermediate field strengths between the 0-flux and π-flux Hubbard model, which, to the best of our knowledge, has not been studied via DQMC.
In this work, we study the Hubbard–Hofstadter model using DQMC and present the evolution of thermodynamic properties of correlated electrons in an orbital magnetic field B. We demonstrate that the gross structure of magnetic Bloch bands and band gaps determined by the non-interacting Hofstadter spectrum is preserved in the presence of U. Moreover, we determine that the many-body fermion sign is directly connected to electronic charge compressibility. Finally, focusing on the half-filled AFMI, we find that an orbital magnetic field tends to delocalize electrons and thus effectively lower the influence of Hubbard U.
Results and discussion
Interacting gap structure
In Fig. 1, we show the electron density 〈n〉 vs. chemical potential μ at different field strengths for U/t = 0–8. In Fig. 1a we observe an electron density plateau where there are energy gaps between magnetic Bloch bands, \(\langle n\rangle =2B\nu /{{{\Phi }}}_{0},\nu \in {\mathbb{Z}}\). As U increases in Fig. 1b–e, these plateaus are weakened and pushed outward in chemical potential, but inflections of the 〈n〉 vs. μ curves are still visible at the same density values, indicating that degeneracy of Landau levels is not modified by U. Since we are at relatively high temperature, only the most prominent band gaps (with Chern number C = ±1) 〈n〉 = 2B/Φ0 remain visible at larger U values. Additionally, at half filling, as U increases, a Mott gap appears and widens for all values of magnetic field. In Fig. 1d, e, for U/t ≥ 6, when the Mott gap is well-defined, it decreases monotonically as the magnetic field increases, consistent with previous exact diagonalization results24. The same trend can be seen in a “correlated Hofstadter butterfly” plot, as shown in Supplementary Fig. S1 and described in Supplementary Note 2. We will discuss later, and in more detail, the behavior of the Mott gap.
It is instructive to plot our data as Wannier diagrams44, i.e., color intensity plots of charge compressibility χ = ∂〈n〉/∂μ as a function of electron density 〈n〉 and magnetic field strength B. Charge compressibility, or thermodynamic density of states, is directly measurable in experiments9,13. In a non-interacting system, at zero temperature, charge compressibility is equivalent to the single-particle density of states. We measure charge compressibility in DQMC simulations as
where ni = ni↑ + ni↓. In Fig. 2, we show Wannier diagrams for U/t = 0–6. For all values of U, we observe local minima of χ (indicating incompresssible states) along straight lines satisfying the Diophantine equation
where r and s are integers, and n0 = 2 is the electron density of the completely filled system. This is consistent with what we expect from the Hofstadter spectrum in the non-interacting system44. The most prominent incompressible state with r = 1, s = 0 remains clearly visible up to U/t = 8. Less prominent incompressible states with r = 2 and r = 3 persist to U/t = 4 and U/t = 2, respectively. These results show that the integer quantum Hall states for r ≤ 3, s = 0 have no weak coupling instabilities with respect to Hubbard repulsion; the r = 1, s = 0 state remains stable up to large U. At half filling, the vertical compressibility minima indicative of the Mott gap becomes visible for U/t ≳ 4. Thus, we argue that the gross structure of magnetic Bloch bands and band gaps determined by the non-interacting Hofstadter spectrum is preserved in the presence of U, with the Mott gap at half-filling when U/t ≳ 4 superimposed as an additional feature. Due to the sign problem, our DQMC simulations are restricted to relatively high temperature βt ≤ 5, so we cannot resolve conclusively how much U changes the fine structure of the Hofstadter spectrum.
Fermion sign
An important quantity in QMC simulations of interacting fermions is the fermion sign. The Hubbard-Hofstadter model is sign-problem-free at half filling on a bipartite lattice. But the fermion sign problem45 fundamentally prevents us from obtaining high quality simulation data at low temperatures and away from half filling. Thus, any insight into factors affecting the severity of the sign problem is valuable. Since the fermion sign problem is NP-hard46, we do not expect a general solution to the fermion sign problem to exist. Nevertheless, as the sign problem is representation-dependent, it is possible to reduce or completely remove the sign problem for specific classes of non-generic Hamiltonians47.
In this work, we find a correlation between the fermion sign and the charge compressibility. In Fig. 3, we show the fermion sign 〈s〉 and charge compressibility χ, both plotted against 〈n〉, for one representative set of parameters. Local minima of charge compressibility in this interacting system exactly correspond to local maxima of the fermion sign. At these local maxima, the fermion sign may be an order of magnitude improved over its value at other electron densities and that of the standard zero-field Hubbard model. For an extended figure demonstrating that this correspondence is general across our parameter space and not a finite size artifact, see Supplementary Fig. S2 and Supplementary Note 3. Our results may mean that although the Hubbard model, in general, suffers from a sign problem, it is possible to obtain good results when we are precisely located on an integer quantum Hall plateau. Since similar sign-compressibility correspondence has been reported30,48,49,50, it appears that the improvement of fermion sign in insulating phases is quite general, consistent with our intuition that fermionic statistics become less important in localized states. Our results also relate to recent work51,52,53 suggesting that the fermion sign is not merely a coincidental barrier to accessing low-temperature physics, but may be reflective of intrinsic physics of model Hamiltonians.
Half-filling
Finally, we focus on half filling, where we believe interesting interplay between Hofstadter physics and Hubbard physics occurs. In the absence of a magnetic field, the ground state of the half-filled Hubbard model is an AFMI at any nonzero value of U, i.e., (Uc = 0)29,30, due to perfect nesting of the Fermi surface and a logarithmically divergent single-particle density of states. In the limit of strong interactions U ≫ t, the half-filled Hubbard model maps to the Heisenberg model with antiferromagnetic nearest neighbor spin exchange energy J = 4t2/U54.
In the presence of an orbital magnetic field, the non-interacting density of states at half filling is modified significantly. As can be seen in Fig. 2a, the density of states/charge compressibility does not change monotonically with field, but instead shows prominent minima at Φ/Φ0 = p/q, where p and q are co-prime and q is even, corresponding to a non-interacting ground state with q inequivalent Dirac cones37. The large-field limit corresponds to the π-flux model, in which a semimetal–AFMI transition occurs at Uc ≈ 5.6t. As the orbital magnetic field significantly changes the non-interacting density of states at half filling, we expect that critical Uc should exhibit B field dependence. It would be interesting to investigate if Uc changes monotonically with B, or if it exhibits non-monotonicity commensurate with the oscillatory behavior of the density of states. We defer the mapping of this U-B phase diagram to future work.
For the remainder of this section, we focus on the parameter region U/t ∈ [6, 10]. Here, the system is safely an AFMI at all field strengths. We examine the evolution of local magnetic moment \(\langle {m}_{z}^{2}\rangle\), antiferromagnetic structure factor S(π, π) and specific heat cv = ∂〈E〉/∂T with magnetic field strength, and see that these thermodynamic quantities all consistently show that a strong orbital magnetic field tends to modify the AFMI by delocalizing electrons and thereby reducing the effect of U on the low-energy properties of the Mott insulating phase. For finite-size analysis of thermodynamic observables at half filling, see Supplementary Fig. S3 and Supplementary Note 3.
In Fig. 4, we show the temperature and field dependence of the local moment. The local moment or sublattice magnetization
measures the degree of spin localization. It is 0.5 in the non-interacting system and approaches 1 in the U/t → ∞ limit. In the zero-field Hubbard model, \(\langle {m}_{z}^{2}\rangle\) has features at T ~ U associated with the formation of local moments, and at T ~ J, associated with short- or long-range ordering of local moments36,55. We see that at fixed U, increasing magnetic field strength reduces the local moment monotonically at all temperatures, with the effect largest below temperatures T ~ J. The zero-field and π-flux limit of local moment data are consistent with previous work43.
The magnetic structure factor is the Fourier transform of the real-space spin–spin correlation function
When the system has long-range antiferromagnetic order, the structure factor is strongly peaked at the ordering wave vector Q = (π, π), with peak height scaling linearly with lattice size36,56. In our simulations, we find that at all B and U values, the magnetic structure factor is sharply peaked at (π, π), consistent with the system being in the AFMI phase. Insets to Fig. 4 show that at all U, the magnetic field monotonically reduces S(π, π), indicating that the magnetic field reduces the AFMI ordering tendencies.
Figure 5 shows the evolution of the low temperature peak of cv with magnetic field and Hubbard U. We calculate cv numerically by measuring energy as a function of temperature 〈E(T)〉 and taking the finite difference Δ〈E〉/ΔT. In the zero-field half-filled Hubbard model, at large U, the specific heat has a “two peak” structure, with a broad high temperature peak at T ~ U associated with charge fluctuations and and a narrow low temperature peak at T ~ J associated with spin fluctuations55,57. When a magnetic field is turned on, as shown in Fig. 5b, c, U is large enough that the two peaks remain well-separated. The high temperature peak doesn’t move, while the low temperature peak shifts to higher temperatures. In Fig. 5a, the two peaks are initially well-separated at low fields, but at Φ/Φ0 ≳ 1/4, the low temperature peak shifts upwards and merges with the high-temperature peak, complicating our interpretation. This is likely due to U/t = 6 being low enough for the system to not simply map to the Heisenberg model, and for the system to be close to the AFMI phase transition at π-flux.
Since the low temperature peak in specific heat is associated with the spin exchange energy J, we are tempted to say that the orbital magnetic field increases J. However, this interpretation may be overly naive. We believe the more accurate statement is that the orbital magnetic field tends to delocalize electrons, and thus, effectively lower the influence of U on low energy properties of the system. Insets to Fig. 5 show that a magnetic field increases (in magnitude) kinetic energy in the insulating phase, which supports this interpretation. This runs contrary to our usual intuition that an orbital magnetic field localizes electrons by winding them up into Landau orbits. However, here our starting point is a correlated insulator, rather than free electrons (or a Fermi liquid). Our results in Figs. 4–5, along with the decreasing width of the Mott gap in Fig. 1d, e, suggest that in the AFMI phase, the orbital magnetic field tends to delocalize electrons, increase kinetic energy, and lower the effective influence of U. We observe that the influence of magnetic field is suppressed as U increases. As U/t → ∞, the influence from the B field will diminish and become negligible, since in the atomic limit, no hopping exists, and the orbital magnetic field cannot have an influence on the system.
Conclusions
In this work, we implemented DQMC to simulate the Hubbard–Hofstadter model and directly investigate field dependent thermodynamic properties of correlated electrons, specifically focusing on the charge compressibility, local moment, magnetic structure factor, and specific heat. By examining charge compressibility, we find that magnetic Bloch bands are smeared out by the Hubbard interaction, but the non-interacting band gaps away from half-filling persist at temperatures accessible to DQMC in the presence of local correlations. At half filling, we find that the orbital magnetic field and Hubbard potential act antagonistically. At intermediate to strong coupling U/t ∈ [6, 10], strong orbital magnetic fields reduce the apparent width of the Mott gap, reduce the magnitude of the local moment and magnetic structure factor, increase kinetic energy, and shift the low-T peak of specific heat to higher temperatures. Together, these phenomena indicate that an orbital magnetic field tends to delocalize electrons and reduce the effect of U. From the algorithmic perspective, we find that the fermion sign in DQMC simulations is improved significantly when the physical system is incompressible.
For future work, we are interested in mapping out the U-B phase diagram for U/t ∈ [0, 6] and in studying the evolution of Uc with magnetic field strength. At smaller values of U, we can achieve lower temperatures in DQMC simulations, but finite size effects also become more significant. Pinning down the precise location of the phase transition will require simulations on much larger lattices, as well as careful finite size scaling analysis.
As this work is partially motivated by understanding the large variety of exotic states in Moiré materials9,10,11,12,13,14, another potential direction for future work is to use methods developed here to study the Hubbard–Hofstadter model on the honeycomb or triangular lattice. It has been suggested that the single-band Hubbard model on a triangular lattice is directly applicable to transition metal dichalcogenide heterobilayers58,59. However, due to the linear dispersion relation of Dirac fermions, the Coulomb interaction is poorly screened in graphene. It is likely that some multi-orbital, extended Hubbard model may be required to capture the full effect of electron correlations, for example, in twisted bilayer graphene60,61.
Methods
We study the single-band Hubbard–Hofstadter model on a two dimensional square lattice
where t is the hopping integral between the nearest neighbor sites 〈ij〉, μ is chemical potential, and U is the on-site Coulomb interaction strength. \({c}_{i\sigma }^{{{{\dagger}}} }\) (ciσ) is the creation (annihilation) operator for an electron on site i with spin σ = ↑, ↓ and \({n}_{i\sigma }={c}_{i\sigma }^{{{{\dagger}}} }{c}_{i\sigma }\) measures the number of electrons of spin σ on site i. As this model only has nearest-neighbor hopping, it preserves particle-hole symmetry at half-filling with μ = 0. A uniform, orbital magnetic field is introduced by the Peierls substitution via the phase
where the integral is taken over the shortest straight line path, Φ0 = h/e is the magnetic flux quantum, and Ri is the position of site i. We choose the symmetric gauge \({{{{{{{\bf{A}}}}}}}}=(-y\hat{{{{{{{{\bf{x}}}}}}}}}+x\hat{{{{{{{{\bf{y}}}}}}}}})B/2\) and do not include any Zeeman coupling terms.
We simulate the Hamiltonian in Eq. (5) on a finite cluster with lattice constant a = 1, and Nx and Ny sites in the x and y directions, respectively. N = NxNy denotes the total number of sites. We implement modified periodic boundary conditions consistent with magnetic translation symmetry62. Requiring that the wave function be single-valued on the torus gives the flux quantization condition Φ/Φ0 = nf/N, where Φ = Ba2 is the flux through a plaquette and nf is an integer.
Allowing the hopping integral to carry a complex phase requires us to modify the standard DQMC algorithm to use complex numbers, which increases the run-time of our algorithm ~3-fold. The complexified DQMC algorithm retains the same O(M3L) scaling as the real DQMC algorithm, where M = Nx = Ny is the linear size of the lattice, and L is the number of imaginary time discretization steps. Unless otherwise specified, all Hubbard–Hofstadter DQMC simulations are performed on a Nx = Ny = 8 square cluster. Error bars in DQMC results, when shown, denote ±1 standard error of the mean, estimated by jackknife resampling. Detailed simulation parameters are listed in Supplementary Note 1. Non-interacting results, where shown, are obtained from diagonalizing the Hofstadter model on a 40 × 40 cluster in order to minimize finite size effects.
Data availability
Aggregated numerical data and analysis routines required to reproduce the figures can be found at https://doi.org/10.5281/zenodo.6383764. Raw simulation data that support the findings of this study are stored on the Sherlock cluster at Stanford University and are available from the corresponding author upon reasonable request.
Code availability
The most up-to-date version of our DQMC simulation code can be accessed at https://github.com/edwnh/dqmc.
References
Wu, T. et al. Magnetic-field-induced charge-stripe order in the high-temperature superconductor YBa2Cu3Oy. Nature 477, 191–194 (2011).
Gerber, S. et al. Three-dimensional charge density wave order in YBa2Cu3O6.67 at high magnetic fields. Science 350, 949–952 (2015).
Jang, H. et al. Ideal charge-density-wave order in the high-field state of superconducting YBCO. Proc. Natl Acad. Sci. USA 113, 14645–14650 (2016).
Edkins, S. D. et al. Magnetic field-induced pair density wave state in the cuprate vortex halo. Science 364, 976–980 (2019).
Grigera, S. A. et al. Magnetic field-tuned quantum criticality in the metallic ruthenate Sr3Ru2O7. Science 294, 329–332 (2001).
Custers, J. et al. The break-up of heavy electrons at a quantum critical point. Nature 424, 524–527 (2003).
Lévy, F., Sheikin, I., Grenier, B. & Huxley, A. D. Magnetic field-induced superconductivity in the ferromagnet URhGe. Science 309, 1343–1346 (2005).
Ran, S. et al. Extreme magnetic field-boosted superconductivity. Nat. Phys. 15, 1250–1254 (2019).
Hunt, B. et al. Massive Dirac fermions and Hofstadter butterfly in a van der Waals heterostructure. Science 340, 1427–1430 (2013).
Wang, L. et al. Evidence for a fractional fractal quantum Hall effect in graphene superlattices. Science 350, 1231–1234 (2015).
Cao, Y. et al. Unconventional superconductivity in magic-angle graphene superlattices. Nature 556, 43–50 (2018).
Spanton, E. M. et al. Observation of fractional Chern insulators in a Van der Waals heterostructure. Science 360, 62–66 (2018).
Yu, J. et al. Correlated Hofstadter spectrum and flavour phase diagram in magic-angle twisted bilayer graphene. Nat. Phys. 18, 825–831 (2022).
Saito, Y. et al. Hofstadter subband ferromagnetism and symmetry-broken Chern insulators in twisted bilayer graphene. Nat. Phys. 17, 478–481 (2021).
Hofstadter, D. R. Energy levels and wave functions of Bloch electrons in rational and irrational magnetic fields. Phys. Rev. B 14, 2239–2249 (1976).
MacDonald, A. H. Landau-level subband structure of electrons on a square lattice. Phys. Rev. B 28, 6713–6717 (1983).
Thouless, D. J., Kohmoto, M., Nightingale, M. P. & den Nijs, M. Quantized Hall conductance in a two-dimensional periodic potential. Phys. Rev. Lett. 49, 405–408 (1982).
Miyake, H., Siviloglou, G. A., Kennedy, C. J., Burton, W. C. & Ketterle, W. Realizing the Harper Hamiltonian with laser-assisted tunneling in optical lattices. Phys. Rev. Lett. 111, 185302 (2013).
Dean, C. R. et al. Hofstadter’s butterfly and the fractal quantum Hall effect in moiré superlattices. Nature 497, 598–602 (2013).
Ponomarenko, L. A. et al. Cloning of Dirac fermions in graphene superlattices. Nature 497, 594–597 (2013).
Gudmundsson, V. & Gerhardts, R. R. Effects of screening on the Hofstadter butterfly. Phys. Rev. B 52, 16744–16752 (1995).
Doh, H. & Salk, S.-H. S. Effects of electron correlations on the Hofstadter spectrum. Phys. Rev. B 57, 1312–1315 (1998).
Barelli, A., Bellissard, J., Jacquod, P. & Shepelyansky, D. L. Double butterfly spectrum for two interacting particles in the Harper model. Phys. Rev. Lett. 77, 4752–4755 (1996).
Czajka, K., Gorczyca, A., Maśka, M. M. & Mierzejewski, M. Hofstadter butterfly for a finite correlated system. Phys. Rev. B 74, 125116 (2006).
Acheche, S., Arsenault, L.-F. & Tremblay, A.-M. S. Orbital effect of the magnetic field in dynamical mean-field theory. Phys. Rev. B 96, 235135 (2017).
Markov, A. A., Rohringer, G. & Rubtsov, A. N. Robustness of the topological quantization of the Hall conductivity for correlated lattice electrons at finite temperatures. Phys. Rev. B 100, 115102 (2019).
Tu, W.-L., Schindler, F., Neupert, T. & Poilblanc, D. Competing orders in the Hofstadter t−J model. Phys. Rev. B 97, 035154 (2018).
Blankenbecler, R., Scalapino, D. J. & Sugar, R. L. Monte Carlo calculations of coupled boson-fermion systems. I. Phys. Rev. D. 24, 2278–2286 (1981).
Hirsch, J. E. Two-dimensional Hubbard model: Numerical simulation study. Phys. Rev. B 31, 4403–4419 (1985).
White, S. R. et al. Numerical study of the two-dimensional Hubbard model. Phys. Rev. B 40, 506–516 (1989).
Jia, C. J. et al. Persistent spin excitations in doped antiferromagnets revealed by resonant inelastic light scattering. Nat. Commun. 5, 3314 (2014).
Kung, Y. F. et al. Doping evolution of spin and charge excitations in the Hubbard model. Phys. Rev. B 92, 195108 (2015).
Khatami, E., Scalettar, R. T. & Singh, R. R. P. Finite-temperature superconducting correlations of the Hubbard model. Phys. Rev. B 91, 241107 (2015).
Huang, E. W. et al. Numerical evidence of fluctuating stripes in the normal state of high-tc cuprate superconductors. Science 358, 1161–1164 (2017).
Huang, E. W., Sheppard, R., Moritz, B. & Devereaux, T. P. Strange metallicity in the doped Hubbard model. Science 366, 987–990 (2019).
Varney, C. N. et al. Quantum Monte Carlo study of the two-dimensional fermion Hubbard model. Phys. Rev. B 80, 075116 (2009).
Wen, X. & Zee, A. Winding number, family index theorem, and electron hopping in a magnetic field. Nucl. Phys. B 316, 641–662 (1989).
Affleck, I. & Marston, J. B. Large-n limit of the Heisenberg–Hubbard model: Implications for high-Tc superconductors. Phys. Rev. B 37, 3774–3777 (1988).
Chang, C.-C. & Scalettar, R. T. Quantum disordered phase near the mott transition in the staggered-flux Hubbard model on a square lattice. Phys. Rev. Lett. 109, 026404 (2012).
Otsuka, Y., Yunoki, S. & Sorella, S. Mott transition in the 2D Hubbard model with π-flux. JPS Conf. Proc. 3, 013021 (2014).
Parisen Toldin, F., Hohenadler, M., Assaad, F. F. & Herbut, I. F. Fermionic quantum criticality in honeycomb and π-flux Hubbard models: Finite-size scaling of renormalization-group-invariant observables from quantum Monte Carlo. Phys. Rev. B 91, 165108 (2015).
Otsuka, Y., Yunoki, S. & Sorella, S. Universal quantum criticality in the Metal–Insulator transition of two-dimensional interacting Dirac electrons. Phys. Rev. X 6, 011029 (2016).
Guo, H. et al. Unconventional pairing symmetry of interacting Dirac fermions on a π-flux lattice. Phys. Rev. B 97, 155146 (2018).
Wannier, G. H. A result not dependent on rationality for Bloch electrons in a magnetic field. Phys. Status Solidi B Basic Res. 88, 757–765 (1978).
Loh, E. Y. et al. Sign problem in the numerical simulation of many-electron systems. Phys. Rev. B 41, 9301–9307 (1990).
Troyer, M. & Wiese, U.-J. Computational complexity and fundamental limitations to fermionic quantum Monte Carlo simulations. Phys. Rev. Lett. 94, 170201 (2005).
Li, Z.-X. & Yao, H. Sign-problem-free fermionic quantum Monte Carlo: Developments and applications. Annu. Rev. Condens. Matter Phys. 10, 337–356 (2019).
Mondaini, R., Bouadim, K., Paiva, T. & dos Santos, R. R. Finite-size effects in transport data from quantum Monte Carlo simulations. Phys. Rev. B 85, 125127 (2012).
Kung, Y. F. et al. Characterizing the three-orbital Hubbard model with determinant quantum Monte Carlo. Phys. Rev. B 93, 155166 (2016).
Huang, E. W., Vaezi, M.-S., Nussinov, Z. & Vaezi, A. Enhanced correlations and superconductivity in weakly interacting partially flat-band systems: A determinantal quantum Monte Carlo study. Phys. Rev. B 99, 235128 (2019).
Mondaini, R., Tarat, S. & Scalettar, R. T. Quantum critical points and the sign problem. Science 375, 418–424 (2022).
Wessel, S., Normand, B., Mila, F. & Honecker, A. Efficient Quantum Monte Carlo simulations of highly frustrated magnets: the frustrated spin-1/2 ladder. SciPost Phys. 3, 005 (2017).
Götz, A., Beyl, S., Hohenadler, M. & Assaad, F. F. Valence-bond solid to antiferromagnet transition in the two-dimensional Su–Schrieffer–Heeger model by langevin dynamics. Phys. Rev. B 105, 085151 (2022).
Fazekas, P. Lecture Notes on Electron Correlation and Magnetism (World Scientific, 1999).
Paiva, T., Scalettar, R. T., Huscroft, C. & McMahan, A. K. Signatures of spin and charge energy scales in the local moment and specific heat of the half-filled two-dimensional Hubbard model. Phys. Rev. B 63, 125116 (2001).
Huse, D. A. Ground-state staggered magnetization of two-dimensional quantum Heisenberg antiferromagnets. Phys. Rev. B 37, 2380–2382 (1988).
Wang, W. O., Ding, J. K., Moritz, B., Huang, E. W. & Devereaux, T. P. Magnon heat transport in a two-dimensional Mott insulator. Phys. Rev. B 105, L161103 (2022).
Wu, F., Lovorn, T., Tutuc, E. & MacDonald, A. H. Hubbard model physics in transition metal dichalcogenide Moiré bands. Phys. Rev. Lett. 121, 026402 (2018).
Tang, Y. et al. Simulation of Hubbard model physics in WSe2/WS2 Moiré superlattices. Nature 579, 353–358 (2020).
Andrews, B. & Soluyanov, A. Fractional quantum Hall states for Moiré superstructures in the Hofstadter regime. Phys. Rev. B 101, 235312 (2020).
Po, H. C., Zou, L., Senthil, T. & Vishwanath, A. Faithful tight-binding models and fragile topology of magic-angle bilayer graphene. Phys. Rev. B 99, 195455 (2019).
Assaad, F. F. Depleted kondo lattices: Quantum Monte Carlo and mean-field calculations. Phys. Rev. B 65, 115104 (2002).
Acknowledgements
We acknowledge helpful discussions with Richard Scalettar, Philip W. Phillips, Allan.H. MacDonald, Benjamin E. Feldman, Young S. Lee, and Jiachen Yu. This work was supported by the U.S. Department of Energy (DOE), Office of Basic Energy Sciences, Division of Materials Sciences and Engineering. E.W.H. was supported by the Gordon and Betty Moore Foundation’s EPiQS Initiative through grants GBMF 4305 and GBMF 8691. Y.S. was supported by the Gordon and Betty Moore Foundation’s EPiQS Initiative through grants GBMF 4302 and GBMF 8686. Computational work was performed on the Sherlock cluster at Stanford University and on resources of the National Energy Research Scientific Computing Center, supported by the U.S. DOE, Office of Science, under Contract no. DE-AC02-05CH11231.
Author information
Authors and Affiliations
Contributions
E.W.H. wrote the DQMC simulation code. J.K.D. performed simulations and data analysis. T.P.D. conceptualized the work. All authors (J.K.D., W.O.W., B.M., Y.S., E.W.H., and T.P.D.) participated in discussions and manuscript writing.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Communications Physics thanks Rok Žitko, Andrey Antipov, and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Ding, J.K., Wang, W.O., Moritz, B. et al. Thermodynamics of correlated electrons in a magnetic field. Commun Phys 5, 204 (2022). https://doi.org/10.1038/s42005-022-00968-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s42005-022-00968-2
This article is cited by
-
Interaction-driven spontaneous ferromagnetic insulating states with odd Chern numbers
npj Quantum Materials (2023)