No Arabic abstract
Lattice Monte Carlo calculations of interacting systems on non-bipartite lattices exhibit an oscillatory imaginary phase known as the phase or sign problem, even at zero chemical potential. One method to alleviate the sign problem is to analytically continue the integration region of the state variables into the complex plane via holomorphic flow equations. For asymptotically large flow times the state variables approach manifolds of constant imaginary phase known as Lefschetz thimbles. However, flowing such variables and calculating the ensuing Jacobian is a computationally demanding procedure. In this paper we demonstrate that neural networks can be trained to parameterize suitable manifolds for this class of sign problem and drastically reduce the computational cost. We apply our method to the Hubbard model on the triangle and tetrahedron, both of which are non-bipartite. At strong interaction strengths and modest temperatures the tetrahedron suffers from a severe sign problem that cannot be overcome with standard reweighting techniques, while it quickly yields to our method. We benchmark our results with exact calculations and comment on future directions of this work.
The repulsive fermionic Hubbard model is a typical model describing correlated electronic systems. Although it is a simple model with only a kinetic term and a local interaction term, their competition generates rich phases. When the interaction part is significant, usually in many strongly correlated, flat or narrow band systems, lots of novel correlated phases may emerge. One way to understand the possible correlated phases is to go beyond finite interaction and solve the infinite-$U$ Hubbard model. Solving infinite-$U$ Hubbard model is usually extremely hard, and a large-scale unbiased numerical study is still missing. In this Letter, we propose a projection approach, such that a controllable quantum Monte Carlo (QMC) simulation on infinite-$U$ Hubbard model may be done at some integer fillings where either it is sign problem free or surprisingly has an algebraic sign structure -- a power law dependence of average sign on system size. We demonstrate our scheme on the infinite-$U$ $SU(2N)$ fermionic Hubbard model on both square and honeycomb lattice at half-filling, where it is sign problem free, and suggest possible correlated ground states. The method can be generalized to study certain extended Hubbard models applying to cluster Mott insulators or 2D Morie systems, among one of them at certain non-half integer filling, the sign has an algebraic behavior such that it can be numerically solved within a polynomial time. Further, our projection scheme can also be generalized to implement the Gutzwiller projection to spin basis such that $SU(2N)$ quantum spin models and Kondo lattice models may be studied in the framework of fermionic QMC simulations.
Quantum gas microscopes for ultracold atoms can provide high-resolution real-space snapshots of complex many-body systems. We implement machine learning to analyze and classify such snapshots of ultracold atoms. Specifically, we compare the data from an experimental realization of the two-dimensional Fermi-Hubbard model to two theoretical approaches: a doped quantum spin liquid state of resonating valence bond type, and the geometric string theory, describing a state with hidden spin order. This approach considers all available information without a potential bias towards one particular theory by the choice of an observable and can therefore select the theory which is more predictive in general. Up to intermediate doping values, our algorithm tends to classify experimental snapshots as geometric-string-like, as compared to the doped spin liquid. Our results demonstrate the potential for machine learning in processing the wealth of data obtained through quantum gas microscopy for new physical insights.
We study the Hubbard model with non-Hermitian asymmetric hopping terms. The conjugate hopping terms are introduced for two spin components so that the negative sign is canceled out. This ensures that the quantum Monte Carlo simulation is free from the negative sign problem. We analyze the antiferromagnetic order and its suppression by the non-Hermiticity.
The solution of complex many-body lattice models can often be found by defining an energy functional of the relevant density of the problem. For instance, in the case of the Hubbard model the spin-resolved site occupation is enough to describe the system total energy. Similarly to standard density functional theory, however, the exact functional is unknown and suitable approximations need to be formulated. By using a deep-learning neural network trained on exact-diagonalization results we demonstrate that one can construct an exact functional for the Hubbard model. In particular, we show that the neural network returns a ground-state energy numerically indistinguishable from that obtained by exact diagonalization and, most importantly, that the functional satisfies the two Hohenberg-Kohn theorems: for a given ground-state density it yields the external potential and it is fully variational in the site occupation.
Complex behavior poses challenges in extracting models from experiment. An example is spin liquid formation in frustrated magnets like Dy$_2$Ti$_2$O$_7$. Understanding has been hindered by issues including disorder, glass formation, and interpretation of scattering data. Here, we use a novel automated capability to extract model Hamiltonians from data, and to identify different magnetic regimes. This involves training an autoencoder to learn a compressed representation of three-dimensional diffuse scattering, over a wide range of spin Hamiltonians. The autoencoder finds optimal matches according to scattering and heat capacity data and provides confidence intervals. Validation tests indicate that our optimal Hamiltonian accurately predicts temperature and field dependence of both magnetic structure and magnetization, as well as glass formation and irreversibility in Dy$_2$Ti$_2$O$_7$. The autoencoder can also categorize different magnetic behaviors and eliminate background noise and artifacts in raw data. Our methodology is readily applicable to other materials and types of scattering problems.