Simulating the Two-Dimensional <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mrow><mml:mi>t</mml:mi><mml:mtext>−</mml:mtext><mml:mi>J</mml:mi></mml:mrow></mml:math> Model at Finite Doping with Neural Quantum States

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.