We present a modification of the continuous-time hybridization-expansion quantum Monte Carlo method, which allows one to take into account spin rotational symmetry. At present, studies of real compounds by means of dynamical mean-field theory employ an approximate form of the Coulomb interaction which breaks the rotational symmetry, leading, in particular, to considerable overestimation of the Curie-Weiss temperature. The obtained temperature dependence of magnetic susceptibility and Curie-Weiss temperature for the three-band model are in good agreement with the exact solution. Owing to its computational efficiency, the presented method is suitable for investigation of magnetic properties of real compounds with strong electronic correlations.
The microscopic mechanism of itinerant ferromagnetism is a long-standing problem due to the lack of non-perturbative methods to handle strong magnetic fluctuations of itinerant electrons. We have non-pertubatively studied thermodynamic properties and magnetic phase transitions of a two-dimensional multi-orbital Hubbard model exhibiting ferromagnetic ground states. Quantum Monte-Carlo simulations are employed, which are proved in a wide density region free of the sign problem usually suffered by simulations for fermions. Both Hund's coupling and electron itinerancy are essential for establishing the ferromagnetic coherence. No local magnetic moments exist in the system as a priori, nevertheless, the spin channel remains incoherent showing the Curie-Weiss type spin magnetic susceptibility down to very low temperatures at which the charge channel is already coherent exhibiting a weakly temperature-dependent compressibility. For the SU(2) invariant systems, the spin susceptibility further grows exponentially as approaching zero temperature in two dimensions. In the paramagnetic phase close to the Curie temperature, the momentum space Fermi distributions exhibit strong resemblance to those in the fully polarized state. The long-range ferromagnetic ordering appears when the symmetry is reduced to the Ising class, and the Curie temperature is accurately determined. These simulations provide helpful guidance to searching for novel ferromagnetic materials in both strongly correlated
d-orbital transition metal oxide layers and the
p-orbital ultra-cold atom optical lattice systems.
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.
Magnetic and electronic properties of the Hubbard model on the Bethe and fcc lattices in infinite dimensions have been investigated numerically on the basis of the dynamical coherent potential approximation (CPA) theory combined with the harmonic approximation (HA) in order to clarify the quantitative aspects of the theory. It is shown that the dynamical CPA+HA reproduces well the sublattice magnetization, the magnetizations, susceptibilities, and the Néel temperatures (
TN) as well as the Curie temperatures calculated by the Quantum Monte-Carlo (QMC) method. The critical Coulomb interactions (
Uc) for the metal-insulator (MI) transition are also shown to agree with the QMC results above
TN. Below
TN,
Uc deviate from the QMC values by about 30% at low temperature regime. These results indicate that the dynamical CPA+HA is applicable to the quantitative description of the magnetic properties in high dimensional systems, but one needs to take into account higher-order dynamical corrections in order to describe the MI transition quantitatively at low temperatures.
Magnetic and electronic properties of the Hubbard model on the Bethe and fcc lattices in infinite dimensions have been investigated numerically on the basis of the dynamical coherent potential approximation (CPA) theory combined with the harmonic approximation (HA) in order to clarify the quantitative aspects of the theory. It is shown that the dynamical CPA+HA reproduces well the sublattice magnetizations, the magnetizations, susceptibilities, and the Néel temperatures ( T N ) as well as the Curie temperatures calculated by the Quantum Monte-Carlo (QMC) method. The critical Coulomb interactions ( U c ) for the metal–insulator (MI) transition are also shown to agree with the QMC results above T N . Below T N , U c deviate from the QMC values by about 30% at low temperature regime. These results indicate that the dynamical CPA+HA is applicable to the quantitative description of the magnetic properties in high dimensional systems, but one needs to take into account higher-order dynamical corrections in order to describe the MI transition quantitatively at low temperatures.
We develop a projective quantum Monte Carlo algorithm of the Hirsch-Fye type for obtaining ground state properties of the Anderson impurity model. This method is employed to solve the self-consistency equations of dynamical mean field theory. It is shown that the approach converges rapidly to the ground state so that reliable zero-temperature results are obtained. As a first application, we study the Mott-Hubbard metal-insulator transition of the frustrated one-band Hubbard model, reconfirming the numerical renormalization group results.