In this Thesis, I study composite quantum impurities with variational techniques, both inspired by machine learning as well as fully analytic. I supplement this with exploration of other applications of machine learning, in particular artificial neural networks, in many-body physics. In Chapters 3 and 4, I study quasiparticle systems with variational approach. I derive a Hamiltonian describing the angulon quasiparticle in the presence of a magnetic field. I apply analytic variational treatment to this Hamiltonian. Then, I introduce a variational approach for non-additive systems, based on artificial neural networks. I exemplify this approach on the example of the polaron quasiparticle (Fröhlich Hamiltonian). In Chapter 5, I continue using artificial neural networks, albeit in a different setting. I apply artificial neural networks to detect phases from snapshots of two types physical systems. Namely, I study Monte Carlo snapshots of multilayer classical spin models as well as molecular dynamics maps of colloidal systems. The main type of networks that I use here are convolutional neural networks, known for their applicability to image data.
In this Thesis, I study composite quantum impurities with variational techniques, both inspired by machine learning as well as fully analytic. I supplement this with exploration of other applications of machine learning, in particular artificial neural networks, in many-body physics. In Chapters 3 and 4, I study quasiparticle systems with variational approach. I derive a Hamiltonian describing the angulon quasiparticle in the presence of a magnetic field. I apply analytic variational treatment to this Hamiltonian. Then, I introduce a variational approach for non-additive systems, based on artificial neural networks. I exemplify this approach on the example of the polaron quasiparticle (Fröhlich Hamiltonian). In Chapter 5, I continue using artificial neural networks, albeit in a different setting. I apply artificial neural networks to detect phases from snapshots of two types physical systems. Namely, I study Monte Carlo snapshots of multilayer classical spin models as well as molecular dynamics maps of colloidal systems. The main type of networks that I use here are convolutional neural networks, known for their applicability to image data.
Most interacting many-body systems in physics are not analytically solvable. Instead, numerical methods are needed for the study of these complex and high-dimensional problems. At present, there are many interesting problems in strongly correlated systems that remain unsolved with current methods. At the heart of this problem is finding an efficient representation that incorporates symmetries, correlations and general features. \n \nIn the context of computer science, machine learning techniques have had astonishing success at reducing the dimensionality of data. The leading method is through the use of artificial neural networks. These networks have been enormously successful at sifting through vast amounts of data to find patterns and regularities. In a sense, neural networks are themselves a statistical system whose properties are adjusted to mimic the features of the data. By finding an effective low-dimensional representation of the data, machine learning has greatly subdued the curse of dimensionality found in many real-world problems. \n \nIn this Thesis, we apply several machine learning techniques to the study of interacting many-body systems in classical and quantum statistical physics. We explore supervised classification of phases of matter with an emphasis on physical interpretation of the net- work. In doing so, we design a custom network architecture that possesses rotational symmetry as an inductive bias. We further investigate connections between the renormalization group and deep learning through applying a super-resolving neural network to the classical Ising model. Towards experimental efforts, we also repurpose generative machine learning to quantum state tomography for the calibration and testing of quantum devices. We conclude with a latent variable model inspired by near-term quantum algorithms. This maps to a variational Monte Carlo ansatz that produces samples efficiently for interacting quantum systems.
Over the past years, machine learning has emerged as a powerful computational tool to tackle complex problems over a broad range of scientific disciplines. In particular, artificial neural networks have been successfully deployed to mitigate the exponential complexity often encountered in quantum many-body physics, the study of properties of quantum systems built out of a large number of interacting particles. In this Article, we overview some applications of machine learning in condensed matter physics and quantum information, with particular emphasis on hands-on tutorials serving as a quick-start for a newcomer to the field. We present supervised machine learning with convolutional neural networks to learn a phase transition, unsupervised learning with restricted Boltzmann machines to perform quantum tomography, and variational Monte Carlo with recurrent neural-networks for approximating the ground state of a many-body Hamiltonian. We briefly review the key ingredients of each algorithm and their corresponding neural-network implementation, and show numerical experiments for a system of interacting Rydberg atoms in two dimensions.
Over the past few years, machine learning has emerged as a powerful computational tool to tackle complex problems in a broad range of scientific disciplines.In particular, artificial neural networks have been successfully used to mitigate the exponential complexity often encountered in quantum many-body physics, the study of properties of quantum systems built from a large number of interacting particles.In this article, we review some applications of neural networks in condensed matter physics and quantum information, with particular emphasis on hands-on tutorials serving as a quick start for a newcomer to the field.The prerequisites of this tutorial are basic probability theory and calculus, linear algebra, basic notions of neural networks, statistical physics, and quantum mechanics.The reader is introduced to supervised machine learning with convolutional neural networks to learn a phase transition, unsupervised learning with restricted Boltzmann machines to perform quantum tomography, and the variational Monte Carlo method with recurrent neural networks for approximating the ground state of a many-body Hamiltonian.For each algorithm, we briefly review the key ingredients and their corresponding neuralnetwork implementation, and show numerical experiments for a system of interacting Rydberg atoms in two dimensions.
Neural network approaches to machine learning problems rely on statistical frameworks that are designed to generalize to many situations, whereas traditional physics-based approaches to quantum problems stem from specific physical laws that govern the behavior of certain classes of systems. While quantum physics and machine learning have been largely disparate fields for most of their history, their interplay has gained considerable attention in recent years, primarily driven by the premise that the most powerful computational models are those whose foundations stem from both physical and statistical principles. Two overarching strategies in this area include building physics-based machine learning models from the ground up, or alternatively, looking for hidden ways that statistical models include aspects of physical systems. Many-body quantum systems turn out to be a useful testing ground for both of these approaches, due to the complex behavior that arises from simple fundamental rules. Here, we present several ways that neural networks can represent many-body spin states, which enables the exploration into the underlying nature of neural networks, as well as gaining insight into how effective different classes of models are for simulating physical systems. The main focus areas of this thesis are (1) connecting convolutional neural networks to maximum-entropy distributions when applied to spin systems, (2) using generative models to gain physical insights into experimental tomographic data, and (3) applying inverse reinforcement learning to probe the nature of agent-based quantum circuit design.
This thesis focuses on the development, application, and high-performance implementation of numerical methods for simulating quantum many-body systems. In the first chapter, I provide an overview of several most successful conventional approaches for simulating quantum many-body systems, along with their current limitations. Then, the basic properties of the quantum spin model, which will be the subject of study in the following chapters, and its ground state are briefly summarized. Additionally, I discuss the background and the motivations of the research projects carried out during my PhD research. Chapter 2 introduces two numerical algorithms that are extensively used in this thesis, which are exact diagonalization (ED) and variational Monte Carlo. I outline the process of constructing a Hamiltonian matrix and employing iterative methods for exact diagonalization. Subsequently, I introduce the mechanism of the four important components in the ``standard" variational Monte Carlo that are widely used in the most recent research. In Chapter 3, we propose a symmetry-projected variational ansatz in the form of linear combinations of simple restricted Boltzmann machines. This construction allows one to explore states outside of the original variational manifold and increase the representation power with moderate computational effort. Besides allowing one to restore spatial symmetries, an expansion in terms of Krylov states using a Lanczos recursion offers a solution that can further improve the quantum state accuracy. We illustrate these ideas with an application to the Heisenberg J1−J2 model on the square lattice, a paradigmatic problem under debate in condensed matter physics,and achieve state-of-the-art accuracy in the representation of the ground state. In Chapter 4, we present a variational approach based on a Chebyshev expansion of the spectral function, and each Chebyshev polynomial is represented by a neural network quantum state. The Chebyshev moments are obtained by recursively applying the Hamiltonian and projecting on the space of variational states using a modified natural gradient descent method. We compare this approach with a modified approximation of the spectral function which uses a Krylov subspace constructed from the "Chebyshev wave functions." We present results for the one-dimensional and two-dimensional Heisenberg models on the square lattice and compare them with those obtained by other methods in the literature. In Chapter 5, we describe a method to carry out thermodynamic simulations of quantum many-body models using typicality and variational representations of quantum states. Minimally entangled typical thermal states (METTS) are a construction that allows one to solve for the imaginary time evolution of quantum many-body systems. By using wave functions that are weakly entangled, one can take advantage of efficient representations in the form of matrixproduct states. We generalize these ideas to arbitrary variational wave functions and we focus, as an illustration, on the particular case of restricted Boltzmann machines. We demonstrate the application of the method to quantum spin systems in one and two spatial dimensions. In Chapter 6, we study the properties of the impurity band in heavily doped nonmagnetic semiconductors using the Jacobi-Davidson algorithm and the supervised deep learning method. The disorder averaged inverse participation ratio and Thouless number calculation show us the rich structure inside the impurity band. A convolutional neural network model, which is trained to distinguish the extended/localized phase of the Anderson model with high accuracy, shows us the results in good agreement with the conventional approach. Together, we find that there are three mobility edges in the impurity band for a specific on-site impurity potential, which means the presence of the extended states while filling the impurity band. In the final chapter, Chapter 7, we discuss the implementation of architecture-aware optimizations for spin dynamics simulations on x86 CPUs and Nvidia GPUs. The time cost of evaluating the spin-spin correlation functions in real space increases quadratically as the number of spins N, leading to significant computational effort, making the simulation of large spin systems very challenging. We propose to use a highly optimized general matrix multiply (GEMM) subroutine to calculate the dynamical spin-spin correlation function that can achieve near-optimal hardware utilization. Furthermore, we fuse the element-wise operations in the calculation of S(q,t) into the in-house GEMM kernel to reduce memory latency. On Nvidia GPU, we perform the on-the-fly calculation in the epilogue of the GEMM subroutine to avoid saving intermediate results on global memory, which significantly reduces memory usage and achieves better performance than cuBLAS.--Author's abstract
This thesis investigates emergent phenomena in strongly correlated quantum matter, with a central focus on developing first-principles deep learning methods for quantum many-body systems. A common theme throughout this work is the search for microscopic mechanisms that generate unexpected collective behavior, together with computational frameworks capable of describing such behavior in regimes where conventional approaches become inadequate. The first part of the thesis studies a doped Mott insulator in a magnetic field and uncovers a new superconducting state, which we call \emph{magnonic superconductivity}. In this state, magnons mediate pairing, and the superconducting order is carried by composite objects binding electron pairs to spin excitations. This mechanism gives rise to a partially spin-polarized superconducting phase beyond the conventional paradigm of electron pairing and further suggests the possibility of higher-charge bound states, including charge-6e superconductivity. The main thrust of the thesis is the use of deep neural networks as variational, first-principles representations of many-body wave functions. In contrast to the usual data-driven role of machine learning, the neural-network ans\"atze developed here do not learn from labeled datasets; instead, they are optimized directly from the microscopic Hamiltonian through variational Monte Carlo. Within this framework, we first develop self-attention-based neural-network wave functions for correlated electrons in solids and show that they accurately describe moir\'e systems with minimal human bias, capture interaction-driven phase transitions, and exhibit favorable scaling of representational complexity with system size. We then introduce a general physics-informed initialization strategy that reconstructs complex many-body wave functions from probability density and probability current. This approach enables highly accurate pretraining of neural wave functions, achieves near-exact overlaps for topological and paired quantum states, and allows efficient simulations of fractional quantum Hall systems with realistic interactions at unprecedented sizes. Finally, we introduce QERNEL, an efficient foundational neural variational ansatz designed to learn a family of many-body ground states within a single conditional model. By combining efficient context conditioning with modern permutation-equivariant building blocks, QERNEL achieves substantially improved expressivity at low computational cost. Applied to interacting electrons in semiconductor moir\'e heterobilayers, it captures distinct many-body phases and resolves the sharp phase transition between them in systems as large as 150 electrons. These results show that context-aware neural wave functions can move beyond Hamiltonian-by-Hamiltonian solution toward genuinely foundational models of quantum matter. Taken together, the results of this thesis show that first-principles deep learning can serve not merely as a numerical tool, but as a general framework for quantum many-body physics: one that bridges microscopic Hamiltonians, expressive wave-function ansätze, and scalable computation, while opening new paths to the discovery and characterization of strongly correlated quantum matter.