We show how to define the appropriate network measure that characterizes the nonmagnetic Mott transition starting from the Hubbard model on a square lattice. To simulate the Mott transition, we employ the variational Monte Carlo method based on the Gutzwiller wave function with nearest-neighbor doublon-holon binding correlations. In order to search the Mott transition point, we use the doublon-holon correlation function as the weighted network representation for analyzing the half-filled electronic system. We then show that a network quantity, small-worldness, can be applied to defining the order parameter describing a hidden symmetry broken in the network representation, which allows us to capture the Mott transition point. As a result, we suggest that the picture of the Mott transition can be obtained using a mass of data entirely contained in topologies of the proper weighted networks.
We investigate the Mott transition in a two-dimensional two-orbital Hubbard model on a square lattice at half filling by using the variational Monte Carlo method. As a variational trial wave function, we consider Gutzwiller's on-site correlations and nearest-neighbor doublon-holon correlations. We calculate the ground-state energy and the other physical quantities, all of which clearly signal a first-order Mott transition. We thus obtain the phase diagram at zero temperature, where we find the coexistence region of the metal and the Mott insulator. We discuss how Hund's coupling affects the nature of the Mott transition in this system.
The mechanism of Mott transition in the Hubbard model on the square lattice is studied without explicit introduction of magnetic and superconducting correlations, using a variational Monte Carlo method. In the trial wave functions, we consider various types of binding factors between a doubly-occupied site (doublon, D) and an empty site (holon, H), like a long-range type as well as a conventional nearest-neighbor type, and add independent long-range D-D (H-H) factors. It is found that a wide choice of D-H binding factor leads to Mott transitions at critical values near the band width. We renew the D-H binding picture of Mott transitions by introducing two characteristic length scales, the D-H binding length l_{DH} and the minimum D-D distance l_{DD}, which we appropriately estimate. A Mott transition takes place at l_{DH}=l_{DD}. In the metallic regime (l_{DH}>l_{DD}), the domains of D-H pairs overlap with one another, thereby doublons and holons can move independently by exchanging the partners one after another. In contrast, the D-D factors give only a minor contribution to the Mott transition.
To understand the mechanism of Mott transitions in case of no magnetic influence, superfluid-insulator (Mott) transitions in the S=0 Bose Hubbard model at unit filling are studied on the square and triangular lattices, using a variational Monte Carlo method. In trial many-body wave functions, we introduce various types of attractive correlation factors between a doubly-occupied site (doublon, D) and an empty site (holon, H), which play a central role for Mott transitions, in addition to the onsite repulsive (Gutzwiller) factor. By optimizing distance-dependent parameters, we study various properties of this type of wave functions. With a hint from the Mott transition arising in a completely D-H bound state, we propose an improved picture of Mott transitions, by introducing two characteristic length scales, the D-H binding length
ξdh and the minimum D-D exclusion length
ξdd. Generally, a Mott transition occurs when
ξdh becomes comparable to
ξdd. In the conductive (superfluid) state, domains of D-H pairs overlap with each other (
ξdh>ξdd); thereby D and H can propagate independently as density carriers by successively exchanging the partners. In contrast, intersite repulsive Jastrow (D-D and H-H) factors have little importance for the Mott transition.
Mechanisms of Mott transitions and d x 2 Ày 2 -wave superconductivity (SC) are studied in the half-filledband Hubbard model on square lattices with a diagonal hopping term (t 0 ), using an optimization (or correlated) variational Monte Carlo method.In the trial wave functions, a doublon-holon binding effect is introduced in addition to the onsite Gutzwiller projection.We mainly treat a d-wave singlet state and a projected Fermi sea.In both wave functions, first-order Mott transitions without direct relevance to magnetic orders occur at U ¼ U c , which is approximately the bandwidth, for arbitrary t 0 =t.These transitions originate in the binding or unbinding of a doublon to a holon.d-wave SC appears in a narrow range immediately below U c .The robust d-wave superconducting correlation is necessarily accompanied by enhanced antiferromagnetic correlation; the strength of SC decreases, as t 0 =t increases.
We use a convolutional restricted Boltzmann machine (CRBM) neural network to construct a variational wave function (WF) for the Hubbard model on a square lattice and study it using the variational Monte Carlo (VMC) method. In the wave function, the CRBM acts as a correlation factor to a mean-field BCS state. The number of variational parameters in the WF does not grow automatically with the lattice size and it is computationally much more efficient compared to other neural network based WFs. We find that in the intermediate to strong coupling regime of the model at half-filling, the wave function outperforms even the highly accurate long range backflow-Jastrow correlated wave function. Using the WF, we study the ground state of the half-filled model as a function of onsite Coulomb repulsion
U. We consider two cases for the next-nearest-neighbor hopping parameter, e.g.,
t′=0 as well as a frustrated model case with
t′=0. By examining several quantities, e.g., double occupancy, charge gap, momentum distribution, and spin-spin correlations, we find that the weekly correlated phase in both cases is paramagnetic metallic (PM). As
U is increased, the system undergoes a first-order Mott transition to an insulating state at a critical
Uc, the value of which depends upon
t′. The Mott state in both cases is spin gapped with long range antiferromagnetic (AF) order. Remarkably, the AF order emerges spontaneously from the wave function which does not have any explicitly broken symmetry in it. Apart from some quantitative differences in the results for the two values of
t′, we find some interesting qualitative differences in the way the Mott transition takes place in the two cases.