Dynamics in a triangular-lattice Hubbard model is studied at half filling on the basis of the cluster dynamical mean field theory. Numerical calculations use a solver of the continuous-time quantum Monte Carlo method based on the strong-coupling expansion. We examine the change in the nearest-neighbor dynamical correlations of doublon and holon in the time domain with varying the Coulomb repulsion near the Mott transition. We demonstrate that the nearest-neighbor doublon-holon pair shows a strong attractive correlation, particularly in the insulating phase, while the nearest-neighbor doublon-doublon pair shows a repulsive correlation. Useful information is provided by the trajectories in the complex plane of the dynamical correlation functions. The short-time dynamics of the nearest-neighbor doublon and holon in the metallic phase indicates larger fluctuations than in the insulating phase. Their trajectories in the complex plane show that the short-time dynamics of the doublon-holon pair has an opposite behavior to that of the doublon-doublon pair. We find that their time scale of dynamics can be characterized by the period in which the phase rotates π around the long-time limit in both phases. In the long-time region, fluctuations persist up to a very long time in the metallic phase, while they quickly vanish in the insulating phase.
Dynamics in a triangular-lattice Hubbard model is studied at half filling on the basis of the cluster dynamical mean field theory. Numerical calculations use a solver of the continuous-time quantum Monte Carlo method based on the strong-coupling expansion. We examine the change in the nearest-neighbor dynamical correlations of doublon and holon in the time domain with varying the Coulomb repulsion near the Mott transition. We demonstrate that the nearest-neighbor doublon-holon pair shows a strong attractive correlation, particularly in the insulating phase, while the nearest-neighbor doublon-doublon pair shows a repulsive correlation. Useful information is provided by the trajectories in the complex plane of the dynamical correlation functions. The short-time dynamics of the nearest-neighbor doublon and holon in the metallic phase indicates larger fluctuations than in the insulating phase. Their trajectories in the complex plane show that the short-time dynamics of the doublon-holon pair has an opposite behavior to that of the doublon-doublon pair. We find that their time scale of dynamics can be characterized by the period in which the phase rotates π around the long-time limit in both phases. In the long-time region, fluctuations persist up to a very long time in the metallic phase, while they quickly vanish in the insulating phase.
We study the dynamics of doublons in the half-filled Hubbard model on the triangular lattice by using cellular dynamical mean-field theory. Investigating the nearest-neighbor dynamical correlations, we demonstrate that a nearest-neighbor doublon-holon pair shows a strong attraction, in particular in the insulating phase. The nearest-neighbor doublon and holon dynamics obtained indicate that the characteristic time scales inversely with the Coulomb interaction in the short-time region in both the metallic and insulating phases. We find that the metallic phase has larger fluctuations in their dynamics. In the long-time region, the metallic phase has persistent fluctuations in various nearest-neighbor configurations, while the fluctuations are vanishingly small in the insulating phase. We also calculate the on-site dynamical correlation of doublons and find that the lifetime of the doublon is longer in the metallic phase than in the insulating phase.
We investigate the Hubbard model on the triangular lattice at half filling using the dynamical cluster approximation (DCA) and dual fermion (DF) methods in combination with continuous-time quantum Monte Carlo (CT QMC) and semiclassical approximation (SCA) methods. We study the one-particle properties and nearest-neighbor spin correlations using the DCA method. We calculate the spectral functions using the CT QMC and SCA methods. The spectral function in the SCA and obtained by analytic continuation of the Pad\'e approximation in CT QMC are in good agreement. We determine the metal-insulator transition and the hysteresis associated with a first-order transition in the double-occupancy and nearest-neighbor spin-correlation functions as functions of temperature. As a further check, we employ the DF method and discuss the advantages and limitation of the dynamical mean-field theory, DCA, and recently developed DF methods by comparing the Green's functions. We find an enhancement of antiferromagnetic correlations and provide evidence for magnetically ordered phases by calculating the spin susceptibility.
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.
We employ dynamical mean-field theory (DMFT) with a quantum Monte Carlo atomic solver to investigate the finite-temperature Mott transition in the Hubbard model with nearest-neighbor hopping on a triangular lattice at half-filling. We estimate the value of the critical interaction to be ${U}_{c}=12.0\ifmmode\pm\else\textpm\fi{}0.5$ in units of the hopping amplitude t through the evolution of the magnetic moment, spectral function, internal energy, and specific heat as the interaction U and temperature T are varied. This work also presents a comparison between DMFT and finite-size determinant quantum Monte Carlo calculations and a discussion of the advantages and limitations of both methods.
The Mott metal-insulator transition in the Hubbard model is studied by constructing a dynamical slave-boson mean-field theory in the limit of large lattice coordination number z that incorporates the binding between doubly occupied (doublon) and empty (holon) sites. On the Mott insulating side where all doublons and holons bond in real space into excitonic pairs leading to the charge gap, the theory simplifies considerably to leading order in 1/z and becomes exact on the infinite-z Bethe lattice. An asymptotic solution is obtained for a continuous Mott transition associated with the closing of the charge gap at a critical value of the Hubbard Uc and the corresponding doublon density ndc, hopping ${\ensuremath{\chi}}_{d}^{c}$ and doublon-holon pairing ${\mathrm{\ensuremath{\Delta}}}_{d}^{c}$ amplitudes. We find ${U}_{c}={U}_{\mathrm{BR}}[1\ensuremath{-}2{n}_{d}^{c}\ensuremath{-}\sqrt{z}({\ensuremath{\chi}}_{d}^{c}+{\mathrm{\ensuremath{\Delta}}}_{d}^{c})]\ensuremath{\simeq}0.8{U}_{\mathrm{BR}}$, where UBR is the critical value for the Brinkman-Rice transition in the Gutzwiller approximation captured in the static mean-field solution of the slave-boson formulation of Kotliar and Ruckenstein. Thus, the Mott transition can be viewed as the quantum correction to the Brinkman-Rice transition due to doublon-holon binding. Quantitative comparisons are made to the results of the dynamical mean-field theory, showing good agreement. In the absence of magnetic order, the Mott insulator is a U(1) quantum spin liquid with nonzero intersite spinon hopping that survives the large-z limit and lifts the 2N-fold degeneracy of the local moments. We show that the spinons are coupled to the doublons/holons by a dissipative compact U(1) gauge field in the deconfined phase, realizing the spin-charge separated gapless spin liquid Mott insulator.