The simulation of interacting many-body quantum systems is a central chal- lenge in condensed matter physics due to the exponential growth of the Hilbert space with system size. In this project, we implement a Neural Network Quantum State (NQS) ansatz based on a Complex Restricted Boltzmann Machine (RBM) to solve for the ground state of the 1D Transverse Field Ising Model (TFIM). Using Variational Monte Carlo (VMC) optimization via Stochastic Reconfiguration (SR), we successfully train the network to approximate the ground state wave function across the paramagnetic-ferromagnetic phase transition. We benchmark our results against Exact Diagonalization (ED) for small systems (N = 12), demonstrating high accuracy in ground state energy estimation. Furthermore, we investigate quantum correlations and order parameters, specifically the transverse magnetization mx and longitudinal spin-spin correlations Czz(r), highlighting the RBM’s ability to capture long-range order and critical behavior near the quantum critical point hx/J = 1.
The simulation of interacting many-body quantum systems is a central chal- lenge in condensed matter physics due to the exponential growth of the Hilbert space with system size. In this project, we implement a Neural Network Quantum State (NQS) ansatz based on a Complex Restricted Boltzmann Machine (RBM) to solve for the ground state of the 1D Transverse Field Ising Model (TFIM). Using Variational Monte Carlo (VMC) optimization via Stochastic Reconfiguration (SR), we successfully train the network to approximate the ground state wave function across the paramagnetic-ferromagnetic phase transition. We benchmark our results against Exact Diagonalization (ED) for small systems (N = 12), demonstrating high accuracy in ground state energy estimation. Furthermore, we investigate quantum correlations and order parameters, specifically the transverse magnetization mx and longitudinal spin-spin correlations Czz(r), highlighting the RBM’s ability to capture long-range order and critical behavior near the quantum critical point hx/J = 1.
The Transverse Field Ising Model (TFIM) is an important quantum many-body model used to study the competition between spin-spin interactions and quantum fluctuations generated by a transverse field. Due to the exponential growth of the Hilbert space with increasing system size, numerical methods quickly become computationally infeasible motivating the use of approximate approaches. In this work Neural Quantum States were used as a variational method to approximate the ground-state properties of the one-dimensional TFIM. The NQS was represented using a feedforward Neural Network and optimized through Monte Carlo sampling together with stochastic reconfiguration. The implementation was done through the NetKet framework. The obtained results showed stable convergence of the variational energy and good agreement with exact diagonalization for the smaller system. The quantum phase transition was observed through the order parameter and spin-spin correlations. Finite-size effects were also observed when comparing the two system sizes. The smaller system more closely approaches the symmetric cat-like ground-state. In contrast, the larger system tends to remain near one symmetry-broken configuration. As a result, the smaller system showed stronger spin expectation fluctuations and a smoother behavior near the phase transition. The result demonstrates that the NQS approach is a good approximate method in studying the ground-state physics of the one-dimensional TFIM.
Recently, quantum-state representation using artificial neural networks has started to be recognized as a powerful tool. However, due to the black-box nature of machine learning, it is difficult to analyze what machine learns or why it is powerful. Here, by applying one of the simplest neural networks, the restricted Boltzmann machine (RBM), to the ground-state representation of the one-dimensional (1D) transverse-field Ising (TFI) model, we make an attempt to directly analyze the optimized network parameters. In the RBM optimization, a zero-temperature quantum state is mapped onto a finite-temperature classical state of the extended Ising spins that constitute the RBM. We find that the quantum phase transition from the ordered phase to the disordered phase in the 1D TFI model by increasing the transverse field is clearly reflected in the behaviors of the optimized RBM parameters and hence in the finite-temperature phase diagram of the classical RBM Ising system. The present finding of a correspondence between the neural-network parameters and quantum phases suggests that a careful investigation of the neural-network parameters may provide a new route to extracting nontrivial physical insights from the neural-network wave functions.
Accuracy of restricted Boltzmann machines for the one-dimensional J1−J2 Heisenberg model
Neural networks have been recently proposed as variational wave functions for quantum many-body systems [G. Carleo and M. Troyer, Science 355, 602 (2017)]. In this work, we focus on a specific architecture, known as restricted Boltzmann machine (RBM), and analyse its accuracy for the spin-1/2 J_1-J_2 J 1 − J 2 antiferromagnetic Heisenberg model in one spatial dimension. The ground state of this model has a non-trivial sign structure, especially for J_2/J_1>0.5 J 2 / J 1 > 0.5 , forcing us to work with complex-valued RBMs. Two variational Ansätze are discussed: one defined through a fully complex RBM, and one in which two different real-valued networks are used to approximate modulus and phase of the wave function. In both cases, translational invariance is imposed by considering linear combinations of RBMs, giving access also to the lowest-energy excitations at fixed momentum k. We perform a systematic study on small clusters to evaluate the accuracy of these wave functions in comparison to exact results, providing evidence for the supremacy of the fully complex RBM. Our calculations show that this kind of Ansätze is very flexible and describes both gapless and gapped ground states, also capturing the incommensurate spin-spin correlations and low-energy spectrum for J_2/J_1>0.5 J 2 / J 1 > 0.5 . The RBM results are also compared to the ones obtained with Gutzwiller-projected fermionic states, often employed to describe quantum spin models [F. Ferrari, A. Parola, S. Sorella and F. Becca, Phys. Rev. B 97, 235103 (2018)]. Contrary to the latter class of variational states, the fully-connected structure of RBMs hampers the transferability of the wave function from small to large clusters, implying an increase in the computational cost with the system size.
Accuracy of restricted Boltzmann machines for the one-dimensional J1−J2 Heisenberg model
Neural networks have been recently proposed as variational wave functions for quantum many-body systems [G. Carleo and M. Troyer, Science 355, 602 (2017)]. In this work, we focus on a specific architecture, known as Restricted Boltzmann Machine (RBM), and analyse its accuracy for the spin-1/2 J1−J2 antiferromagnetic Heisenberg model in one spatial dimension. The ground state of this model has a non-trivial sign structure, especially for J2/J1>0.5, forcing us to work with complex-valued RBMs. Two variational Ans\"atze are discussed: one defined through a fully complex RBM, and one in which two different real-valued networks are used to approximate modulus and phase of the wave function. In both cases, translational invariance is imposed by considering linear combinations of RBMs, giving access also to the lowest-energy excitations at fixed momentum k. We perform a systematic study on small clusters to evaluate the accuracy of these wave functions in comparison to exact results, providing evidence for the supremacy of the fully complex RBM. Our calculations show that this kind of Ans\"atze is very flexible and describes both gapless and gapped ground states, also capturing the incommensurate spin-spin correlations and low-energy spectrum for J2/J1>0.5. The RBM results are also compared to the ones obtained with Gutzwiller-projected fermionic states, often employed to describe quantum spin models [F. Ferrari, A. Parola, S. Sorella and F. Becca, Phys. Rev. B 97, 235103 (2018)]. Contrary to the latter class of variational states, the fully-connected structure of RBMs hampers the transferability of the wave function from small to large clusters, implying an increase of the computational cost with the system size.
Neural network quantum states provide a novel representation of the many-body states of interacting quantum systems and open up a promising route to solve frustrated quantum spin models that evade other numerical approaches. Yet its capacity to describe complex magnetic orders with large unit cells has not been demonstrated, and its performance in a rugged energy landscape has been questioned. Here we apply restricted Boltzmann machines (RBMs) and stochastic gradient descent to seek the ground states of a compass spin model on the honeycomb lattice, which unifies the Kitaev model, Ising model and the quantum 120° model with a single tuning parameter. We report calculation results on the variational energy, order parameters and correlation functions. The phase diagram obtained is in good agreement with the predictions of tensor network ansatz, demonstrating the capacity of RBMs in learning the ground states of frustrated quantum spin Hamiltonians. The limitations of the calculation are discussed. A few strategies are outlined to address some of the challenges in machine learning frustrated quantum magnets.
Neural network quantum states (NQS) have been widely applied to spin-1/2 systems, where they have proven to be highly effective. The application to systems with larger on-site dimension, such as spin-1 or bosonic systems, has been explored less and predominantly using spin-1/2 Restricted Boltzmann Machines (RBMs) with a one-hot/unary encoding. Here, we propose a more direct generalization of RBMs for spin-1 that retains the key properties of the standard spin-1/2 RBM, specifically trivial product states representations, labeling freedom for the visible variables and gauge equivalence to the tensor network formulation. To test this new approach, we present variational Monte Carlo (VMC) calculations for the spin-1 anti-ferromagnetic Heisenberg (AFH) model and benchmark it against the one-hot/unary encoded RBM demonstrating that it achieves the same accuracy with substantially fewer variational parameters. Furthermore, we investigate how the hidden unit complexity of NQS depend on the local single-spin basis used. Exploiting the tensor network version of our RBM we construct an analytic NQS representation of the Affleck-Kennedy-Lieb-Tasaki (AKLT) state in the xyz spin-1 basis using only M=2N hidden units, compared to M∼O(N2) required in the Sz basis. Additional VMC calculations provide strong evidence that the AKLT state in fact possesses an exact compact NQS representation in the xyz basis with only M=N hidden units. These insights help to further unravel how to most effectively adapt the NQS framework for more complex quantum systems.