We investigate variational learning of quantum many-body ground states directly in measurement space using autoregressive neural networks. In particular, we represent quantum states via probability distributions of outcomes over a symmetric informationally complete positive operator-valued measure (SIC-POVM). The probability distribution is encoded in the parameters of an autoregressive neural-network-based on gated recurrent units (GRUs). Ground states are obtained by gradient descent that updates the probability distribution to minimize the energy with respect to a given Hamiltonian, while enforcing positivity constraints that ensure that the distribution of measurement outcomes correspond to a physical quantum state. We analyze the role of constraint enforcement (hierarchy of positivity conditions), variety of neural network architectures (multiple layers, dilation, and modifications of input data) in determining the success of this approach. We benchmark our approach on one-dimensional transverse-field Ising model and the Heisenberg model, along with gapping fields, for system sizes up to L=128, illustrating its efficacy across a wide variety of models.
We investigate variational learning of quantum many-body ground states directly in measurement space using autoregressive neural networks. In particular, we represent quantum states via probability distributions of outcomes over a symmetric informationally complete positive operator-valued measure (SIC-POVM). The probability distribution is encoded in the parameters of an autoregressive neural-network-based on gated recurrent units (GRUs). Ground states are obtained by gradient descent that updates the probability distribution to minimize the energy with respect to a given Hamiltonian, while enforcing positivity constraints that ensure that the distribution of measurement outcomes correspond to a physical quantum state. We analyze the role of constraint enforcement (hierarchy of positivity conditions), variety of neural network architectures (multiple layers, dilation, and modifications of input data) in determining the success of this approach. We benchmark our approach on one-dimensional transverse-field Ising model and the Heisenberg model, along with gapping fields, for system sizes up to L=128, illustrating its efficacy across a wide variety of models.
Solving for the lowest energy eigenstate of the many-body Schrodinger equation is a cornerstone problem that hinders understanding of a variety of quantum phenomena. The difficulty arises from the exponential nature of the Hilbert space which casts the governing equations as an eigenvalue problem of exponentially large, structured matrices. Variational methods approach this problem by searching for the best approximation within a lower-dimensional variational manifold. In this work we use graph neural networks to define a structured variational manifold and optimize its parameters to find high quality approximations of the lowest energy solutions on a diverse set of Heisenberg Hamiltonians. Using graph networks we learn distributed representations that by construction respect underlying physical symmetries of the problem and generalize to problems of larger size. Our approach achieves state-of-the-art results on a set of quantum many-body benchmark problems and works well on problems whose solutions are not positive-definite. The discussed techniques hold promise of being a useful tool for studying quantum many-body systems and providing insights into optimization and implicit modeling of exponentially-sized objects.
Recent theoretical progress has established conditions under which machine learning models can efficiently predict ground-state properties of gapped local Hamiltonians when trained on quantum-generated data. Previous experimental demonstrations in this paradigm, however, have largely been limited to small systems or highly structured states, due to the difficulty of preparing many-body ground states on quantum processors. In this work, we demonstrate learning from experimental quantum data generated from approximate ground states of the two-dimensional Heisenberg XXZ model with system sizes up to 115 qubits. We construct a dataset of single-site expectation values, two-point correlations, and 12-body loop correlations across the antiferromagnetic phase. We then train neural networks on this data and show that they can accurately predict spatially resolved observables for previously unseen Hamiltonian parameters, both within the training distribution and in an out-of-distribution regime approaching the phase boundary. Our results demonstrate the practical realization of learning from quantum data for an interacting two-dimensional many-body system at scale, motivating a path toward regimes where quantum processors could provide training data beyond the reach of classical approximation methods.
Quantum many-body computation is a central challenge across multiple domains of scientific research, and Monte Carlo method has emerged as one of the most effective numerical techniques for investigating large-scale many-body systems. In recent years, the integration of Monte Carlo methods with machine learning has led to the development of neural quantum state approach. Quantum neural networks (QNNs) offer superior expressiveness compared to classical neural networks. In this paper, we propose a novel Quantum Neural Network Quantum State (QNNQS) method, based on the variational Monte Carlo framework. This approach leverages a quantum neural network to represent the many-body quantum state, and optimizes it iteratively through gradient descent to determine the ground-state energy. To assess the performance of this method, we conduct simulations of the one-dimensional transverse-field Ising model, investigating the influence of quantum circuit architecture on circuit expressibility, entangling capability, and computational accuracy. The results demonstrate that the QNNQS method achieves comparable simulation accuracy while utilizing significantly fewer parameters and requiring less training resources than classical network quantum state methods.
Abstract In this work we apply deep neural networks to find the non-equilibrium steady state solution to correlated open quantum many-body systems. Motivated by the ongoing search to find more powerful representations of (mixed) quantum states, we design a simple prototypical convolutional neural network and show that parametrizing the density matrix directly with more powerful models can yield better variational ansatz functions and improve upon results reached by neural density operator based on the restricted Boltzmann machine. Hereby we give up the explicit restriction to positive semi-definite density matrices. However, this is fulfilled again to good approximation by optimizing the parameters. The great advantage of this approach is that it opens up the possibility of exploring more complex network architectures that can be tailored to specific physical properties. We show how translation invariance can be enforced effortlessly and reach better results with fewer parameters. We present results for the dissipative one-dimensional transverse-field Ising model and a two-dimensional dissipative Heisenberg model compared to exact values.
Thermodynamics is a theory of principles that permits a basic description of the macroscopic properties of a rich variety of complex systems from traditional ones, such as crystalline solids, gases, liquids, and thermal machines, to more intricate systems such as living organisms and black holes to name a few. Physical quantities of interest, or equilibrium state variables, are linked together in equations of state to give information on the studied system, including phase transitions, as energy in the forms of work and heat, and/or matter are exchanged with its environment, thus generating entropy. A more accurate description requires different frameworks, namely, statistical mechanics and quantum physics to explore in depth the microscopic properties of physical systems and relate them to their macroscopic properties. These frameworks also allow to go beyond equilibrium situations. Given the notably increasing complexity of mathematical models to study realistic systems, and their coupling to their environment that constrains their dynamics, both analytical approaches and numerical methods that build on these models show limitations in scope or applicability. On the other hand, machine learning, i.e., data-driven, methods prove to be increasingly efficient for the study of complex quantum systems. Deep neural networks, in particular, have been successfully applied to many-body quantum dynamics simulations and to quantum matter phase characterization. In the present work, we show how to use a variational autoencoder (VAE)—a state-of-the-art tool in the field of deep learning for the simulation of probability distributions of complex systems. More precisely, we transform a quantum mechanical problem of many-body state reconstruction into a statistical problem, suitable for VAE, by using informationally complete positive operator-valued measure. We show, with the paradigmatic quantum Ising model in a transverse magnetic field, that the ground-state physics, such as, e.g., magnetization and other mean values of observables, of a whole class of quantum many-body systems can be reconstructed by using VAE learning of tomographic data for different parameters of the Hamiltonian, and even if the system undergoes a quantum phase transition. We also discuss challenges related to our approach as entropy calculations pose particular difficulties.
Neural quantum states (NQS) provide a flexible and scalable framework for approximating quantum many-body wavefunctions. Among NQS parameterizations, autoregressive models are especially attractive because they enable exact, independent sampling from the Born distribution, avoiding the autocorrelation and mixing issues of Markov chain methods. Yet their optimization remains comparatively underexplored: Adam is a scalable method but ignores function space geometry, while stochastic reconfiguration is principled but costly and numerically fragile in large models. To address this gap, we show that variational energy minimization can be viewed as an advantage policy-gradient problem over the Born distribution, motivating trust-region optimization for NQS training. We introduce Proximal Wavefunction Optimization (PWO), a principled trust-region algorithm that clips probability-ratio changes in the amplitude channel and phase increments in the phase channel. PWO avoids explicit matrix inversion, reuses samples across multiple updates, and combines the scalability of first-order optimization with theoretical guarantees. Across Ising and frustrated J1-J2 one- and two-dimensional spin systems, PWO improves stability and wall-clock convergence over Adam, minSR, and SPRING. Finally, we fine-tune a 1.5B-parameter RWKV-7 model, demonstrating NQS optimization at a scale over three orders of magnitude beyond prior work.