Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2512.04663

Fermionic neural Gibbs states

Jannes Nys, Juan Carrasquilla

quant-phcond-mat.str-elcs.LGphysics.comp-phSubmitted Dec 4, 2025• Announced Dec 5, 2025• 0 citations
arXivHTMLPDF

Abstract

We introduce fermionic neural Gibbs states (fNGS), a variational framework for modeling finite-temperature properties of strongly interacting fermions. fNGS starts from a reference mean-field thermofield-double state and uses neural-network transformations together with imaginary-time evolution to systematically build strong correlations. Applied to the doped Fermi-Hubbard model, a minimal lattice model capturing essential features of strong electronic correlations, fNGS accurately reproduces thermal energies over a broad range of temperatures, interaction strengths, even at large dopings, for system sizes beyond the reach of exact methods. These results demonstrate a scalable route to studying finite-temperature properties of strongly correlated fermionic systems beyond one dimension with neural-network representations of quantum states.

Related papers

Thermalization Dynamics in the Two-Dimensional Hubbard Model with Neural-Network Quantum States

Alessandro Sinibaldi, Luciano Loris Viteritti, Riccardo Rende +2 more

Thermalization in strongly correlated fermionic systems remains a central open problem in quantum many-body physics. In this work, we investigate the real-time dynamics and the approach to thermalization in the two-dimensional Hubbard model, a paradigmatic framework for correlated electrons, relevant to high-temperature superconductivity and ultracold quantum simulation. Focusing on the half-filled square lattice, we monitor the time evolution of the double occupancy following a quench in the on-site interaction $U$, and assess whether its long-time value is captured by a canonical thermal ensemble. We employ time-dependent variational Monte Carlo methods combined with transformer-based Neural-Network Quantum States to accurately describe the nonequilibrium dynamics of fermions, especially for the behavior at long times, thereby accessing regimes that were previously inaccessible to numerical simulations. Our results reveal two dynamical behaviors: for weak to intermediate interactions, the double occupancy rapidly approaches the thermal prediction, consistent with ergodic evolution; beyond a critical interaction $U_{C}$, the dynamics remains distinct from the thermal expectation on the timescales investigated, revealing signatures of a prethermal plateau delaying fast relaxation. These results establish numerical simulation as a powerful tool to probe nonequilibrium quantum phenomena in correlated fermionic matter.

cond-mat.str-elCross-listedarXiv:2606.05293•Jun 5, 2026•0 citations

Thermal tensor network simulations of lattice fermions with fixed filling

Qiaoyi Li, Dai-Wei Qu, Bin-Bin Chen +2 more

Numerical simulations of strongly correlated fermions at finite temperature are essential for studying high-temperature superconductivity and other quantum many-body phenomena. The recently developed tangent-space tensor renormalization group (tanTRG) provides an efficient and accurate framework by representing thermal density operators as matrix product operators. However, the particle number generally varies during the cooling process. The conventional strategy of fine-tuning chemical potentials to reach a target filling is computationally demanding. Here we propose a fixed-$N$ tanTRG algorithm that stabilizes the average particle number by adaptively tuning the chemical potential within the imaginary-time evolution. We benchmark its accuracy on exactly solvable free fermions, and further apply it to the square-lattice Hubbard model. For hole-doped cases, we study the temperature evolution of charge and spin correlations, identifying several characteristic temperature scales for stripe formation. Our results establish fixed-$N$ tanTRG as an efficient and reliable tool for finite-temperature studies of correlated fermion systems.

arXiv:2511.07303•Physical review. B./Physical review. B•Feb 2, 2026•2 citations

Purifying Deep Boltzmann Machines for Thermal Quantum States

Yusuke Nomura, Nobuyuki Yoshioka, Franco Nori

We develop two cutting-edge approaches to construct deep neural networks representing the purified finite-temperature states of quantum many-body systems. Both methods commonly aim to represent the Gibbs state by a highly expressive neural-network wave function, exemplifying the idea of purification. The first method is an entirely deterministic approach to generate deep Boltzmann machines representing the purified Gibbs state exactly. This strongly assures the remarkable flexibility of the ansatz which can fully exploit the quantum-to-classical mapping. The second method employs stochastic sampling to optimize the network parameters such that the imaginary time evolution is well approximated within the expressibility of neural networks. Numerical demonstrations for transverse-field Ising models and Heisenberg models show that our methods are powerful enough to investigate the finite-temperature properties of strongly correlated quantum many-body systems, even when the problematic effect of frustration is present.

arXiv:2103.04791•Physical Review Letters•Aug 4, 2021•35 citations

Gaussian variational method to Fermi-Hubbard model in one and two dimensions

Yue-Ran Shi, Yuan-Yao He, Ruijin Liu +1 more

Abstract The study of ground-state properties of the Fermi-Hubbard model is a long-lasting task in the research of strongly correlated systems, and is crucial for the understanding of notable quantum phenomena including superconductivity and magnetism. Owing to the exponentially growing complexity of the system, a quantitative analysis usually demands high computational cost and is restricted to small samples, especially in two or higher dimensions. Here, we introduce a variational method in the frame of fermionic Gaussian states, and obtain the ground states of one- and two-dimensional attractive Hubbard models via imaginary-time evolution. We calculate the total energy and benchmark the results in a wide range of interaction strength and filling factor with those obtained via exact two-body results, the density matrix renormalization group based on matrix product states (MPS), and projector Quantum Monte Carlo (QMC) method. For both 1D and 2D cases, the Gaussian variational method presents accurate results for total energy with a maximum systematic error ∼ 4 % in the intermediate interaction region. The accuracy of these results has negligible dependence on the system size. We further calculate the double occupancy and find excellent agreement with MPS and QMC, as well as the experimental results of cold quantum gases in optical lattices. The results suggest that the Gaussian pairing state is a good approximation to the ground states of attractive Hubbard model, in particular in the strong and weak coupling limits. Moreover, we generalize the method to the attractive Hubbard model with a finite spin-polarization, which can be mapped to the repulsive interaction case via particle-hole transformation, and obtain accurate results of ground state energy and double occupancy. Our work demonstrates the ability of the Gaussian variational method to extract ground state properties of strongly correlated many-body systems with negligible computational cost, especially of large size and in higher dimensions.

arXiv:2308.06700•New Journal of Physics•Sep 1, 2023•2 citations

Variational Simulation of Lattice Fermions with Neural Quantum States

Imelda Romero

Neural quantum states provide a flexible and scalable framework for strongly correlated systems. This thesis develops and benchmarks variational neural ansÀtze for fermionic lattice models, targeting ground states and low-lying excited states, with a particular focus on architectural performance. The constructions embed known structure such as fermionic antisymmetry, lattice translations, and total momentum while learning the remaining correlations efficiently in log amplitude. In the first part, we introduce a single dimensionless, intensive accuracy indicator, the V-score. Defined from the mean energy and energy variance of a state, it enables objective, method-agnostic comparisons across system sizes, interaction strengths, and even distinct momentum sectors, including regimes without exact references. We use it as a unifying measure of variational quality throughout the thesis. In the second part, we introduce a symmetry-aware Slaterâ Backflowâ Jastrow ansatz for spinless fermions in the t-V model. A convolutional backflow is built to be translation equivariant and consistent with particle exchange, allowing quantum-number projection to a fixed momentum without additional forward passes. This state is used to target ground states and momentum-resolved excitations, and to probe ordering tendencies via structure factors. In the third part, we develop a momentum-projected vision-transformer ansatz (DefeViT) for the one-hole t-J family. The state is factorized in a fixed-hole frame with a neural amplitude for the mobile defect, then restored to a chosen total momentum by summing lattice translations with the appropriate fermionic phase. On square lattices and their diagonal extensions (t'-J'), the approach captures ground-state properties and single-hole dispersions and tracks how spin correlators and hole-centered bond textures evolve under frustration. Within this setting, physics-informed neural states provide compact and accurate descriptions of ground states and low-lying excited states.

•Jan 1, 2025•0 citations

Neural Network Backflow with Low-Rank Multi-Determinant Updates

Tianshu Huang, Xiaowei Ou, Vidvuds Ozolins

Simulating strongly correlated fermions remains a long-standing challenge due to the exponential complexity of the Hilbert space and the intricate sign structure of many-body wavefunctions. We introduce a variational framework centered on a neural network backflow transformation that combines deep learning with variational Monte Carlo. The proposed ansatz employs a multi-determinant expansion with low-rank shifts to capture non-local correlations and complex sign structures. Applied to the two-dimensional Hubbard model at both half-filling and $1/8$ doping, the method achieves energies within $0.45\%$ of auxiliary-field quantum Monte Carlo at half-filling and captures intertwined charge- and spin-density stripe patterns at $1/8$ doping. These results demonstrate the potential of this framework as a scalable and interpretable approach to variational simulations of strongly correlated fermionic systems.

cond-mat.str-elCross-listedarXiv:2609.26544•Sep 23, 2026•0 citations
We introduce fermionic neural Gibbs states (fNGS), a variational framework for modeling finite-temperature properties of strongly interacting fermions. fNGS starts from a reference mean-field thermofield-double state and uses neural-network transformations together with imaginary-time evolution to systematically build strong correlations. Applied to the doped Fermi-Hubbard model, a minimal lattice model capturing essential features of strong electronic correlations, fNGS accurately reproduces thermal energies over a broad range of temperatures, interaction strengths, even at large dopings, for system sizes beyond the reach of exact methods. These results demonstrate a scalable route to studying finite-temperature properties of strongly correlated fermionic systems beyond one dimension with neural-network representations of quantum states.
Thermalization in strongly correlated fermionic systems remains a central open problem in quantum many-body physics. In this work, we investigate the real-time dynamics and the approach to thermalization in the two-dimensional Hubbard model, a paradigmatic framework for correlated electrons, relevant to high-temperature superconductivity and ultracold quantum simulation. Focusing on the half-filled square lattice, we monitor the time evolution of the double occupancy following a quench in the on-site interaction UUU, and assess whether its long-time value is captured by a canonical thermal ensemble. We employ time-dependent variational Monte Carlo methods combined with transformer-based Neural-Network Quantum States to accurately describe the nonequilibrium dynamics of fermions, especially for the behavior at long times, thereby accessing regimes that were previously inaccessible to numerical simulations. Our results reveal two dynamical behaviors: for weak to intermediate interactions, the double occupancy rapidly approaches the thermal prediction, consistent with ergodic evolution; beyond a critical interaction UCU_{C}UC​, the dynamics remains distinct from the thermal expectation on the timescales investigated, revealing signatures of a prethermal plateau delaying fast relaxation. These results establish numerical simulation as a powerful tool to probe nonequilibrium quantum phenomena in correlated fermionic matter.
Numerical simulations of strongly correlated fermions at finite temperature are essential for studying high-temperature superconductivity and other quantum many-body phenomena. The recently developed tangent-space tensor renormalization group (tanTRG) provides an efficient and accurate framework by representing thermal density operators as matrix product operators. However, the particle number generally varies during the cooling process. The conventional strategy of fine-tuning chemical potentials to reach a target filling is computationally demanding. Here we propose a fixed-NNN tanTRG algorithm that stabilizes the average particle number by adaptively tuning the chemical potential within the imaginary-time evolution. We benchmark its accuracy on exactly solvable free fermions, and further apply it to the square-lattice Hubbard model. For hole-doped cases, we study the temperature evolution of charge and spin correlations, identifying several characteristic temperature scales for stripe formation. Our results establish fixed-NNN tanTRG as an efficient and reliable tool for finite-temperature studies of correlated fermion systems.
We develop two cutting-edge approaches to construct deep neural networks representing the purified finite-temperature states of quantum many-body systems. Both methods commonly aim to represent the Gibbs state by a highly expressive neural-network wave function, exemplifying the idea of purification. The first method is an entirely deterministic approach to generate deep Boltzmann machines representing the purified Gibbs state exactly. This strongly assures the remarkable flexibility of the ansatz which can fully exploit the quantum-to-classical mapping. The second method employs stochastic sampling to optimize the network parameters such that the imaginary time evolution is well approximated within the expressibility of neural networks. Numerical demonstrations for transverse-field Ising models and Heisenberg models show that our methods are powerful enough to investigate the finite-temperature properties of strongly correlated quantum many-body systems, even when the problematic effect of frustration is present.
Abstract The study of ground-state properties of the Fermi-Hubbard model is a long-lasting task in the research of strongly correlated systems, and is crucial for the understanding of notable quantum phenomena including superconductivity and magnetism. Owing to the exponentially growing complexity of the system, a quantitative analysis usually demands high computational cost and is restricted to small samples, especially in two or higher dimensions. Here, we introduce a variational method in the frame of fermionic Gaussian states, and obtain the ground states of one- and two-dimensional attractive Hubbard models via imaginary-time evolution. We calculate the total energy and benchmark the results in a wide range of interaction strength and filling factor with those obtained via exact two-body results, the density matrix renormalization group based on matrix product states (MPS), and projector Quantum Monte Carlo (QMC) method. For both 1D and 2D cases, the Gaussian variational method presents accurate results for total energy with a maximum systematic error ∼ 4 % in the intermediate interaction region. The accuracy of these results has negligible dependence on the system size. We further calculate the double occupancy and find excellent agreement with MPS and QMC, as well as the experimental results of cold quantum gases in optical lattices. The results suggest that the Gaussian pairing state is a good approximation to the ground states of attractive Hubbard model, in particular in the strong and weak coupling limits. Moreover, we generalize the method to the attractive Hubbard model with a finite spin-polarization, which can be mapped to the repulsive interaction case via particle-hole transformation, and obtain accurate results of ground state energy and double occupancy. Our work demonstrates the ability of the Gaussian variational method to extract ground state properties of strongly correlated many-body systems with negligible computational cost, especially of large size and in higher dimensions.
Neural quantum states provide a flexible and scalable framework for strongly correlated systems. This thesis develops and benchmarks variational neural ansÀtze for fermionic lattice models, targeting ground states and low-lying excited states, with a particular focus on architectural performance. The constructions embed known structure such as fermionic antisymmetry, lattice translations, and total momentum while learning the remaining correlations efficiently in log amplitude. In the first part, we introduce a single dimensionless, intensive accuracy indicator, the V-score. Defined from the mean energy and energy variance of a state, it enables objective, method-agnostic comparisons across system sizes, interaction strengths, and even distinct momentum sectors, including regimes without exact references. We use it as a unifying measure of variational quality throughout the thesis. In the second part, we introduce a symmetry-aware Slaterâ Backflowâ Jastrow ansatz for spinless fermions in the t-V model. A convolutional backflow is built to be translation equivariant and consistent with particle exchange, allowing quantum-number projection to a fixed momentum without additional forward passes. This state is used to target ground states and momentum-resolved excitations, and to probe ordering tendencies via structure factors. In the third part, we develop a momentum-projected vision-transformer ansatz (DefeViT) for the one-hole t-J family. The state is factorized in a fixed-hole frame with a neural amplitude for the mobile defect, then restored to a chosen total momentum by summing lattice translations with the appropriate fermionic phase. On square lattices and their diagonal extensions (t'-J'), the approach captures ground-state properties and single-hole dispersions and tracks how spin correlators and hole-centered bond textures evolve under frustration. Within this setting, physics-informed neural states provide compact and accurate descriptions of ground states and low-lying excited states.
Simulating strongly correlated fermions remains a long-standing challenge due to the exponential complexity of the Hilbert space and the intricate sign structure of many-body wavefunctions. We introduce a variational framework centered on a neural network backflow transformation that combines deep learning with variational Monte Carlo. The proposed ansatz employs a multi-determinant expansion with low-rank shifts to capture non-local correlations and complex sign structures. Applied to the two-dimensional Hubbard model at both half-filling and 1/81/81/8 doping, the method achieves energies within 0.45%0.45\%0.45% of auxiliary-field quantum Monte Carlo at half-filling and captures intertwined charge- and spin-density stripe patterns at 1/81/81/8 doping. These results demonstrate the potential of this framework as a scalable and interpretable approach to variational simulations of strongly correlated fermionic systems.