Simulating large, strongly interacting fermionic systems remains a major challenge for existing numerical methods. In this Letter, we introduce Gutzwiller projected hidden fermion determinant states (G-HFDS) to simulate the strongly interacting limit of the Fermi-Hubbard model, namely the t-J model, across the entire doping regime. We demonstrate that the G-HFDS achieve energies competitive with matrix product states on lattices as large as 10×10 sites while using several orders of magnitude fewer parameters, suggesting the potential for efficient application to even larger system sizes. This remarkable efficiency enables us to probe low-energy physics across the full doping range, providing new insights into the competition between kinetic and magnetic interactions and the nature of emergent quasiparticles. Starting from the low-doping regime, where magnetic polarons dominate the low energy physics, we track their evolution with increasing doping and different next-nearest neighbor hopping amplitudes through analyses of spin and polaron correlation functions as well as the Fermi surface. Our findings demonstrate the potential of determinant-based neural quantum states with an inherent fermionic sign structure, opening the way for simulating large-scale fermionic systems at any particle filling.
Simulating large, strongly interacting fermionic systems remains a major challenge for existing numerical methods. In this Letter, we introduce Gutzwiller projected hidden fermion determinant states (G-HFDS) to simulate the strongly interacting limit of the Fermi-Hubbard model, namely the t-J model, across the entire doping regime. We demonstrate that the G-HFDS achieve energies competitive with matrix product states on lattices as large as 10×10 sites while using several orders of magnitude fewer parameters, suggesting the potential for efficient application to even larger system sizes. This remarkable efficiency enables us to probe low-energy physics across the full doping range, providing new insights into the competition between kinetic and magnetic interactions and the nature of emergent quasiparticles. Starting from the low-doping regime, where magnetic polarons dominate the low energy physics, we track their evolution with increasing doping and different next-nearest neighbor hopping amplitudes through analyses of spin and polaron correlation functions as well as the Fermi surface. Our findings demonstrate the potential of determinant-based neural quantum states with an inherent fermionic sign structure, opening the way for simulating large-scale fermionic systems at any particle filling.
Ultracold alkaline-earth atoms and molecules now enable experimental realizations of SU(N)-symmetric Fermi-Hubbard models, yet theoretical understanding of these systems, particularly at finite doping remains limited. Here we investigate the strong-coupling limit of the SU(3) symmetric Fermi-Hubbard model on the triangular lattice across the full doping range. Using a three-flavor extension of Gutzwiller-projected hidden fermion determinant states (G-HFDS), a neural network based variational ansatz, we analyze two- and three-point spin-spin and spin-spin-hole correlations of the SU(3) Cartan generators. We further study the structure of a pair of doped holes for large periodic systems, and compare our results to the paradigmatic SU(2) square lattice equivalent, finding strikingly similar magnetic correlations, non-s-wave pairing symmetry, and enhanced binding energies. Our results provide a foundation for future exploration of doped SU(N) Mott insulators, providing insights for both theoretical developments and quantum simulation experiments.
Whether the ground state of the square lattice Hubbard model exhibits superconductivity remains a major open question, central to understanding high temperature cuprate superconductors and ultra-cold fermions in optical lattices. Numerical studies have found evidence for stripe-ordered states and superconductivity at strong coupling but the phase diagram remains controversial. Here, we show that one can resolve the subtle energetics of metallic, superconducting, and stripe phases using a new class of neural quantum state (NQS) wavefunctions that extend hidden fermion determinant states to Pfaffians. We simulate several hundred electrons using fast Pfaffian algorithms allowing us to measure off-diagonal long range order. At strong coupling and low hole-doping, we find that a non-superconducting filled stripe phase prevails, while superconductivity coexisting with partially-filled stripes is stabilized by a negative next neighbor hopping t-prime, with |t-prime| > 0.1. At larger doping levels, we introduce momentum-space correlation functions to mitigate finite size effects that arise from weakly-bound pairs. These provide evidence for uniform d-wave superconductivity at U = 4, even when t-prime = 0. Our results highlight the potential of NQS approaches, and provide a fresh perspective on superconductivity in the square lattice Hubbard model.
Understanding the physics underlying high-temperature superconductivity in cuprates and, more recently, infinite-layer nickelates has remained a central challenge in condensed-matter physics. We establish neural quantum states (NQS), specifically Hidden Fermion Determinant States (HFDS), as a scalable variational approach to the three-band Emery model of the copper- and nickel-oxide layers in these materials. After benchmarking HFDS against matrix product states on width-two geometries, we study ground states of fully two-dimensional (2D) systems of up to 10×10 unit cells (300 sites). We characterize the momentum-space distribution of dopants and find a pronounced electron-hole dichotomy similar to cuprate experiments. We further downfold the three-band model to effective single-band descriptions by constructing interacting Wannier functions. We consider a wide range of parameters -- from the charge-transfer regime relevant to cuprates to the Hubbard-Mott regime of nickelates, as well as systematic scans of the charge-transfer gap that has been demonstrated to impact the critical superconducting temperatures. Across all regimes, the effective model significantly deviates from the usual Fermi-Hubbard model: The typical ratio U/t is enhanced, some parameters experience a significant doping dependence, and sizable terms beyond the conventional Hubbard model are present, most notably a density-assisted hopping tn. Notably, in all effective models, tn has the largest contribution to particle-hole asymmetry, rather than next-nearest-neighbor hopping contributions. The effective parameters sensitively depend on the charge-transfer energy, doping, and interaction ratios. Our results establish HFDS as an efficient tool for studying the 2D Emery model and demonstrate that single-band descriptions can require interaction terms generated by the underlying multi-band models.
The two-dimensional Hubbard model was originally proposed as the minimal model of uniform d-wave superconductivity in the cuprate high-temperature superconductors, but recent calculations with relevant parameters in a realistic doping regime have often found other ground states. Further, uncertainty remains over such results due to the intrinsic bias of different computational methods. Here we address this uncertainty through state-of-the-art neural tensor network and neural backflow state simulations, enabled by large-scale GPU computation, that together reach an unprecedented variational accuracy on lattices with up to 1024 sites. At 1/8 doping in the pure Hubbard model, as well as with next-nearest-neighbor hopping relevant to the hole-doped cuprates, we reveal a near degeneracy among a low-energy uniform d-wave state, superconducting stripes, and non-superconducting stripe states on the milli-t energy scale. Our work rebalances the narrative that the Hubbard model supports a single dominant order in relevant parameter regimes towards a picture of extreme competition between stripes and superconductivity.
Abstract Numerically simulating large, spinful, fermionic systems is of great interest in condensed matter physics. However, the exponential growth of the Hilbert space dimension with system size renders exact quantum state parameterizations impractical. Owing to their representative power, neural networks often allow to overcome this exponential scaling. Here, we investigate the ability of neural quantum states (NQS) to represent the bosonic and fermionict − Jmodel – the high interaction limit of the Hubbard model – on various 1D and 2D lattices. Using autoregressive, tensorized recurrent neural networks (RNNs), we study ground state representations upon hole doping the half-filled system. Additionally, we propose a method to calculate quasiparticle dispersions, applicable to any network architecture or lattice geometry, and allowing to infer the low-energy physics from NQS. By analyzing the strengths and weaknesses of the RNN ansatz we shed light on the challenges and promises of NQS for simulating bosonic and fermionic systems.
The recent development of Neural Quantum States (NQS) has established them as one of the most accurate methods for studying strongly correlated many-fermion systems, outperforming existing many-body approaches for large systems. However, NQS calculations are currently highly resource intensive. In this work, we introduce an improved Pareto frontier of efficiency and accuracy for NQS in the simulation of strongly correlated lattice fermions. This frontier is defined by two complementary backflow-related architectures: the Sparse Convolutional Ansatz for Lattice Electrons (SCALE), which achieves high computational efficiency, and the Accurate Convolutional ansatz for lattice Electrons (ACE), which delivers high variational accuracy, based on benchmark results on the Hubbard model and the t – J model on large lattice sizes. SCALE utilizes a tailored convolutional design to enable an efficient local update based on the low-rank update of determinants. This structure reduces the computational scaling from O ( N 4 ) to O ( N 3 ) in backflow methods, resulting in a practical speedup exceeding 40 × in our test cases while retaining accurate variational accuracy. As an application, we study the 1 / 8 -doped pure Hubbard model up to size 32 × 32 and find no significant energy difference between the horizontal and vertical unit-filled stripe states, in contrast to the half-filled stripe state when next-nearest-neighbor hoppings are included. ACE, conversely, employs a deep convolutional stack to maximize expressive power, achieving high accuracy on large systems. We perform extensive benchmark calculations on the Hubbard and t – J models with the two architectures. On the challenging repulsive Hubbard model, SCALE achieves variational energies competitive with leading methods at a fraction of the computational cost. Meanwhile, ACE achieves lower variational energies than previously reported results while requiring only one-sixth of the run-time for system size 16 × 4 . These NQS approaches provide , , and tools for exploring the rich physics of strongly correlated fermionic systems.