A recently introduced one-dimensional three-orbital Hubbard model displays orbital-selective Mott phases with exotic spin arrangements such as spin block states [J. Rincón et al., Phys. Rev. Lett. 112, 106405 (2014)PRLTAO0031-900710.1103/PhysRevLett.112.106405]. In this publication we show that the constrained-path quantum Monte Carlo (CPQMC) technique can accurately reproduce the phase diagram of this multiorbital one-dimensional model, paving the way to future CPQMC studies in systems with more challenging geometries, such as ladders and planes. The success of this approach relies on using the Hartree-Fock technique to prepare the trial states needed in CPQMC. We also study a simplified version of the model where the pair-hopping term is neglected and the Hund coupling is restricted to its Ising component. The corresponding phase diagrams are shown to be only mildly affected by the absence of these technically difficult-to-implement terms. This is confirmed by additional density matrix renormalization group and determinant quantum Monte Carlo calculations carried out for the same simplified model, with the latter displaying only mild fermion sign problems. We conclude that these methods are able to capture quantitatively the rich physics of the several orbital-selective Mott phases (OSMP) displayed by this model, thus enabling computational studies of the OSMP regime in higher dimensions, beyond static or dynamic mean-field approximations.
Using the density matrix renormalization group technique we study the effect of spin-orbit coupling on a three-orbital Hubbard model in the
(t2g)4 sector and in one dimension. Fixing the Hund coupling to a robust value compatible with some multiorbital materials, we present the phase diagram varying the Hubbard
U and spin-orbit coupling $\ensuremath{\lambda}$, at zero temperature. Our results are shown to be qualitatively similar to those recently reported using the dynamical mean-field theory in higher dimensions, providing a robust basis to approximate many-body techniques. Among many results, we observe an interesting transition from an orbital-selective Mott phase to an excitonic insulator with increasing $\ensuremath{\lambda}$ at intermediate
U. In the strong
U coupling limit, we find a nonmagnetic insulator with an effective angular momentum $\ensuremath{\langle}{({\mathbf{J}}^{\mathrm{eff}})}^{2}\ensuremath{\rangle}\ensuremath{\ne}0$ near the excitonic phase, smoothly connected to the $\ensuremath{\langle}{({\mathbf{J}}^{\mathrm{eff}})}^{2}\ensuremath{\rangle}=0$ regime. We also provide a list of quasi-one-dimensional materials where the physics discussed in this paper could be realized.
We study nonlocal correlations in a three-orbital Hubbard model defined on an extended one-dimensional chain using determinant quantum Monte Carlo and density matrix renormalization group methods. We focus on a parameter regime with robust Hund's coupling, which produces an orbital selective Mott phase (OSMP) at intermediate values of the Hubbard
U, as well as an orbitally ordered ferromagnetic insulating state at stronger coupling. An examination of the orbital- and spin-correlation functions indicates that the orbital ordering occurs before the onset of magnetic correlations in this parameter regime as a function of temperature. In the OSMP, we find that the self-energy for the itinerant electrons is momentum dependent, indicating a degree of nonlocal correlations while the localized electrons have largely momentum independent self-energies. These nonlocal correlations also produce relative shifts of the holelike and electronlike bands within our model. The overall momentum dependence of these quantities is strongly suppressed in the orbitally ordered insulating phase.
Using the dynamical mean field theory with Lanczos method as its impurity solver, we study the orbital-selective Mott transition (OSMT) in the two-orbital J model and Jz model. In the multi-orbital systems, the Mott metal-insulator transition occurs successively when the widths of the bands are different. As the narrow orbital becomes Mott insulator while the wide orbital is still in metallic phase, we find an orbital-selective Mott phase (OSMP). There are two different Hubbard models that are frequently used to describe the OSMT, which are named J model and Jz model, respectively. The J Model is composed of the whole Hund's rule coupling terms, including the spin-flip term, the pair-hopping term and the Ising type Hund's rule coupling term. However, there is only Ising type Hund's rule coupling term in the Jz model.#br#We study the ratio of bandwidth W2/W1 on the OSMT by analyzing the results of the density of states and quasi-particle weight. Comparing the phase diagrams obtained from the J and Jz models with the Hund's rule coupling J(Jz)=U/4, we find that the OSMP region of the J model is much larger than that of the Jz model when W2/W1=0.5 or W2/W1=0.2. When the ratio of bandwidth increases to W2/W1=0.8, the OSMP disappears completely in the Jz model. However in the J model, we can still find the OSMT but the area of the OSMP shrinks significantly. Therefore, the OSMT happens more easily in the J model than in the Jz model.#br#In order to discuss the cooperative effect of the bandwidth and Hund's rule coupling on the OSMT, we compare the phase diagrams for different Hund's rule couplings J(Jz)=U/4 and J(Jz)=U/2. We find that when the bandwidth W2/W1≥q 0.7, the OSMT disappears in Jz model in the case of either Jz=U/4 or Jz=U/2. However, the OSMP always exists in the J model if the bandwidths of the two orbitals are different, suggesting that the rotation invariances of the Hund's rule couplings can protect the OSMP. Therefore, one should be more careful when using the Jz model instead of the J model to study the OSMP.
The recent detailed study of quasi-one-dimensional iron-based ladders, with the
3d iron electronic density
n=6, has unveiled surprises, such as orbital-selective phases. However, similar studies for
n=6 iron chains are still rare. Here a three-orbital electronic Hubbard model was constructed to study the magnetic and electronic properties of the quasi-one-dimensional
n=6 iron chain
Ce2O2FeSe2, with focus on the effect of doping. Specifically, introducing the Hubbard
U and Hund
JH couplings and studying the model via the density matrix renormalization group, we report the ground-state phase diagram varying the electronic density away from
n=6. For the realistic Hund coupling
JH/U=1/4, several electronic phases were obtained, including a metal, orbital-selective Mott, and Mott insulating phases. Doping away from the parent phase, the competition of many tendencies leads to a variety of magnetic states, such as ferromagnetism, as well as several antiferromagnetic and magnetic ``block'' phases. In the hole-doping region, two different interesting orbital-selective Mott phases were found: OSMP1 (with one localized orbital and two itinerant orbitals) and OSMP2 (with two localized orbitals and one itinerant orbital). Moreover, charge disproportionation phenomena were found in special doping regions. We argue that our predictions can be tested by simple modifications in the original chemical formula of
Ce2O2FeSe2.
The canonical one-band Hubbard model is studied using a computational method that mixes the Monte Carlo procedure with the mean field approximation. This technique allows us to incorporate thermal fluctuations and the development of short-range magnetic order above ordering temperatures, contrary to the crude finite-temperature Hartree-Fock approximation, which incorrectly predicts a N\'eel temperature
TN that grows linearly with the Hubbard
U/t. The effective model studied here contains quantum and classical degrees of freedom. It thus belongs to the ``spin fermion'' model family widely employed in other contexts. Using exact diagonalization, supplemented by the traveling cluster approximation, for the fermionic sector, and classical Monte Carlo for the classical fields, the Hubbard
U/t vs temperature
T/t phase diagram is studied employing large three- and two-dimensional clusters. We demonstrate that the method is capable of capturing the formation of local moments in the normal state without long-range order, the nonmonotonicity of
TN with increasing
U/t, the development of gaps and pseudogaps in the density of states, and the two-peak structure in the specific heat. Extensive comparisons with determinant quantum Monte Carlo results suggest that the present approach is qualitatively, and often quantitatively, accurate, particularly at intermediate and high temperatures. Finally, we study the Hubbard model including plaquette diagonal hopping (i.e., the $t\text{\ensuremath{-}}{t}^{\ensuremath{'}}$ Hubbard model) in two dimensions and show that our approach allows us to study low-temperature properties where determinant quantum Monte Carlo fails due to the fermion sign problem. Future applications of this method include multiorbital Hubbard models such as those needed for iron-based superconductors.