Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Journal of Experimental and Theoretical Physics Letters

Preserving spin rotational symmetry in quantum Monte Carlo for description of magnetic properties of strongly correlated compounds

A. G. Dyachenkо, Alexander S. Belozerov, Vladimir I. Anisimov

Published Mar 1, 2015• 0 citations
Publisher

Abstract

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.

DOI
10.1134/s0021364015050082

Related papers

Magnetism of iron and nickel from rotationally invariant Hirsch-Fye quantum Monte Carlo calculations

Alexander S. Belozerov, I. Leonov, Vladimir I. Anisimov

We present a rotationally invariant Hirsch-Fye quantum Monte Carlo algorithm in which the spin rotational invariance of Hund's exchange is approximated by averaging over all possible directions of the spin quantization axis. We employ this technique to perform benchmark calculations for the two- and three-band Hubbard models on the infinite-dimensional Bethe lattice. Our results agree quantitatively well with those obtained using the continuous-time quantum Monte Carlo method with rotationally invariant Coulomb interaction. The proposed approach is employed to compute the electronic and magnetic properties of paramagnetic $\ensuremath{\alpha}$ iron and nickel. The obtained Curie temperatures agree well with experiment. Our results indicate that the magnetic transition temperature is significantly overestimated by using the density-density type of Coulomb interaction.

arXiv:1210.0911•Physical Review B•Mar 25, 2013•49 citations

Quantum Monte Carlo study for multiorbital systems with preserved spin and orbital rotational symmetries

志朗 酒井, Ryotaro Arita, Karsten Held +1 more

We propose to combine the Trotter decomposition and a series expansion of the partition function for Hund's exchange coupling in a quantum Monte Carlo (QMC) algorithm for multiorbital systems that preserves spin and orbital rotational symmetries. This enables us to treat the Hund's (spin-flip and pair-hopping) terms, which is difficult in the conventional QMC method. To demonstrate this, we first apply the algorithm to study ferromagnetism in the two-orbital Hubbard model within the dynamical mean-field theory (DMFT). The result reveals that the preservation of the SU(2) symmetry in Hund's exchange is important, where the Curie temperature is grossly overestimated when the symmetry is degraded, as is often done, to Ising $({\mathrm{Z}}_{2})$. We then calculate the ${t}_{2g}$ spectral functions of ${\mathrm{Sr}}_{2}{\mathrm{RuO}}_{4}$ by a three-band DMFT calculation with tight-binding parameters taken from the local density approximation with proper rotational symmetry.

arXiv:cond-mat/0605526•Physical Review B•Oct 4, 2006•41 citations

Benchmarkings for a semiclassical impurity solver for dynamical-mean-field theory: Self-energies and magnetic transitions of the single-orbital Hubbard model

Satoshi Okamoto, Andreas Fuhrmann, Armin Comanac +1 more

An investigation is presented of the utility of semiclassical approximations for solving the quantum-impurity problems arising in the dynamical-mean-field approach to correlated-electron models. The method is based on performing an exact numerical integral over the zero-Matsubara-frequency component of the spin part of a continuous Hubbard-Stratonovich field, along with a spin-field-dependent steepest descent treatment of the charge part. We test this method by applying it to one- or two-site approximations to the single-band Hubbard model with different band structures, and comparing the results to quantum Monte Carlo and simplified exact diagonalization calculations. The resulting electron self-energies, densities of states, and magnetic transition temperatures show reasonable agreement with the quantum Monte Carlo simulation over wide parameter ranges, suggesting that the semiclassical method is useful for obtaining a reasonable picture of the physics in situations where other techniques are too expensive.

arXiv:cond-mat/0502067•Physical Review B•Jun 24, 2005•17 citations

Realistic theory of electronic correlations in nanoscopic systems

Malte Schüler, Stefan Barthel, Tim Oliver Wehling +3 more

Nanostructures with open shell transition metal or molecular constituents host often strong electronic correlations and are highly sensitive to atomistic material details. This tutorial review discusses method developments and applications of theoretical approaches for the realistic description of the electronic and magnetic properties of nanostructures with correlated electrons. First, the implementation of a flexible interface between density functional theory and a variant of dynamical mean field theory (DMFT) highly suitable for the simulation of complex correlated structures is explained and illustrated. On the DMFT side, this interface is largely based on recent developments of quantum Monte Carlo and exact diagonalization techniques allowing for efficient descriptions of general four fermion Coulomb interactions, reduced symmetries and spin-orbit coupling, which are explained here. With the examples of the Cr (001) surfaces, magnetic adatoms, and molecular systems it is shown how the interplay of Hubbard U and Hund’s J determines charge and spin fluctuations and how these interactions drive different sorts of correlation effects in nanosystems. Non-local interactions and correlations present a particular challenge for the theory of low dimensional systems. We present our method developments addressing these two challenges, i.e., advancements of the dynamical vertex approximation and a combination of the constrained random phase approximation with continuum medium theories. We demonstrate how non-local interaction and correlation phenomena are controlled not only by dimensionality but also by coupling to the environment which is typically important for determining the physics of nanosystems.

arXiv:1707.08333•The European Physical Journal Special Topics•Jul 1, 2017•25 citations

Continuous-time impurity solvers

James E. Gubernatis, Naoki Kawashima, Philipp Werner

This chapter, like the previous one, describes quantum Monte Carlo methods for the simulation of a system of interacting electrons at nonzero temperature. Here, we target a particular class of problems, characterized by a small number of correlated sites or orbitals, that we call impurity problems. While we could apply the methods of the previous chapter to such impurity models, the algorithms discussed here are more advantageous not only because they work in continuous imaginary time and hence lack the Trotter approximation error inherent to the previous methods, but also because they manipulate determinants of matrices of smaller size, which can be handled more efficiently. At the same time, the continuous-time approach is applicable to broader classes of impurity problems. Driving the development of these continuous-time impurity solvers was the desire to perform more efficient simulations of correlated lattice models within the framework of dynamical mean-field theory. This formalism maps the lattice problem onto an impurity problem, whose parameters are determined by a self-consistency condition. In some cases, the self-consistent solution gives a good approximation of the properties of the original lattice problem. When we combine this dynamical mean-field approximation with the local-density approximation for electronic structure calculations, we obtain a powerful scheme for electronic structure calculations of strongly correlated materials. In this chapter we discuss the most important classes of impurity models, sketch the basics of the dynamical mean-field approximation, and detail the different variants of continuous-time impurity solvers. Quantum impurity models A quantum impurity model describes an atom or molecule embedded in some host with which it exchanges electrons, spin, and energy. This exchange allows the impurity to make transitions between different quantum states, and in the presence of interactions within or between the impurity orbitals, these transitions lead to nontrivial dynamical properties. Quantum impurity models play a prominent role, for example, in the theoretical description of the magnetic and electric properties of dilute metal alloys and in theoretical studies of quantum dots and molecular conductors. These models also appear as an auxiliary problem whose solution yields the dynamical mean-field description of correlated lattice models.

Cambridge University Press eBooks•Jan 31, 2016•1 citation

Continuous-time quantum Monte Carlo solver for dynamical mean field theory in the compact Legendre representation

Evan Sheridan, Cédric Weber, Evgeny A. Plekhanov +1 more

Dynamical mean-field theory (DMFT) is one of the most widely-used methods to treat accurately electron correlation effects in ab-initio real material calculations. Many modern large-scale implementations of DMFT in electronic structure codes involve solving a quantum impurity model with a continuous-time quantum Monte Carlo (CT-QMC) solver [Rubtsov et al., Phys. Rev. B 72, 035122 (2005); Werner et al., Phys. Rev. Lett. 97, 076405 (2006); Werner and Millis, Phys. Rev. B 74, 155107 (2006); Gull et al., Rev. Mod. Phys. 83, 349 (2011)]. The main advantage of CT-QMC is that, unlike standard quantum Monte Carlo approaches, it is able to generate the local Green's functions $G(\ensuremath{\tau})$ of the correlated system on an arbitrarily fine imaginary time $\ensuremath{\tau}$ grid, and is free of any systematic errors. In this work, we extend a hybrid QMC solver proposed by Khatami et al. [Phys. Rev. E 81, 056703 (2010)] and Rost et al. [Phys. Rev. E 87, 053305 (2013)] to a multiorbital context. This has the advantage of enabling impurity solver QMC calculations to scale linearly with inverse temperature $\ensuremath{\beta}$, and permit its application to $d$- and $f$-band materials. In addition, we present a Green's-function processing scheme which generates accurate quasicontinuous imaginary time solutions of the impurity problem which overcome errors inherent to standard QMC approaches. This solver and processing scheme are incorporated into a full DFT+DMFT calculation using the castep DFT code [Clark et al., Z. Kristallogr. 220, 567 (2005)]. Benchmark calculations for ${\mathrm{SrVO}}_{3}$ properties are presented. The computational efficiency of this method is also demonstrated.

arXiv:1811.12869•Physical review. B./Physical review. B•May 30, 2019•14 citations
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.
We present a rotationally invariant Hirsch-Fye quantum Monte Carlo algorithm in which the spin rotational invariance of Hund's exchange is approximated by averaging over all possible directions of the spin quantization axis. We employ this technique to perform benchmark calculations for the two- and three-band Hubbard models on the infinite-dimensional Bethe lattice. Our results agree quantitatively well with those obtained using the continuous-time quantum Monte Carlo method with rotationally invariant Coulomb interaction. The proposed approach is employed to compute the electronic and magnetic properties of paramagnetic $\ensuremath{\alpha}$ iron and nickel. The obtained Curie temperatures agree well with experiment. Our results indicate that the magnetic transition temperature is significantly overestimated by using the density-density type of Coulomb interaction.
We propose to combine the Trotter decomposition and a series expansion of the partition function for Hund's exchange coupling in a quantum Monte Carlo (QMC) algorithm for multiorbital systems that preserves spin and orbital rotational symmetries. This enables us to treat the Hund's (spin-flip and pair-hopping) terms, which is difficult in the conventional QMC method. To demonstrate this, we first apply the algorithm to study ferromagnetism in the two-orbital Hubbard model within the dynamical mean-field theory (DMFT). The result reveals that the preservation of the SU(2) symmetry in Hund's exchange is important, where the Curie temperature is grossly overestimated when the symmetry is degraded, as is often done, to Ising (Z2)({\mathrm{Z}}_{2})(Z2​). We then calculate the t2g{t}_{2g}t2g​ spectral functions of Sr2RuO4{\mathrm{Sr}}_{2}{\mathrm{RuO}}_{4}Sr2​RuO4​ by a three-band DMFT calculation with tight-binding parameters taken from the local density approximation with proper rotational symmetry.
An investigation is presented of the utility of semiclassical approximations for solving the quantum-impurity problems arising in the dynamical-mean-field approach to correlated-electron models. The method is based on performing an exact numerical integral over the zero-Matsubara-frequency component of the spin part of a continuous Hubbard-Stratonovich field, along with a spin-field-dependent steepest descent treatment of the charge part. We test this method by applying it to one- or two-site approximations to the single-band Hubbard model with different band structures, and comparing the results to quantum Monte Carlo and simplified exact diagonalization calculations. The resulting electron self-energies, densities of states, and magnetic transition temperatures show reasonable agreement with the quantum Monte Carlo simulation over wide parameter ranges, suggesting that the semiclassical method is useful for obtaining a reasonable picture of the physics in situations where other techniques are too expensive.
Nanostructures with open shell transition metal or molecular constituents host often strong electronic correlations and are highly sensitive to atomistic material details. This tutorial review discusses method developments and applications of theoretical approaches for the realistic description of the electronic and magnetic properties of nanostructures with correlated electrons. First, the implementation of a flexible interface between density functional theory and a variant of dynamical mean field theory (DMFT) highly suitable for the simulation of complex correlated structures is explained and illustrated. On the DMFT side, this interface is largely based on recent developments of quantum Monte Carlo and exact diagonalization techniques allowing for efficient descriptions of general four fermion Coulomb interactions, reduced symmetries and spin-orbit coupling, which are explained here. With the examples of the Cr (001) surfaces, magnetic adatoms, and molecular systems it is shown how the interplay of Hubbard U and Hund’s J determines charge and spin fluctuations and how these interactions drive different sorts of correlation effects in nanosystems. Non-local interactions and correlations present a particular challenge for the theory of low dimensional systems. We present our method developments addressing these two challenges, i.e., advancements of the dynamical vertex approximation and a combination of the constrained random phase approximation with continuum medium theories. We demonstrate how non-local interaction and correlation phenomena are controlled not only by dimensionality but also by coupling to the environment which is typically important for determining the physics of nanosystems.
This chapter, like the previous one, describes quantum Monte Carlo methods for the simulation of a system of interacting electrons at nonzero temperature. Here, we target a particular class of problems, characterized by a small number of correlated sites or orbitals, that we call impurity problems. While we could apply the methods of the previous chapter to such impurity models, the algorithms discussed here are more advantageous not only because they work in continuous imaginary time and hence lack the Trotter approximation error inherent to the previous methods, but also because they manipulate determinants of matrices of smaller size, which can be handled more efficiently. At the same time, the continuous-time approach is applicable to broader classes of impurity problems. Driving the development of these continuous-time impurity solvers was the desire to perform more efficient simulations of correlated lattice models within the framework of dynamical mean-field theory. This formalism maps the lattice problem onto an impurity problem, whose parameters are determined by a self-consistency condition. In some cases, the self-consistent solution gives a good approximation of the properties of the original lattice problem. When we combine this dynamical mean-field approximation with the local-density approximation for electronic structure calculations, we obtain a powerful scheme for electronic structure calculations of strongly correlated materials. In this chapter we discuss the most important classes of impurity models, sketch the basics of the dynamical mean-field approximation, and detail the different variants of continuous-time impurity solvers. Quantum impurity models A quantum impurity model describes an atom or molecule embedded in some host with which it exchanges electrons, spin, and energy. This exchange allows the impurity to make transitions between different quantum states, and in the presence of interactions within or between the impurity orbitals, these transitions lead to nontrivial dynamical properties. Quantum impurity models play a prominent role, for example, in the theoretical description of the magnetic and electric properties of dilute metal alloys and in theoretical studies of quantum dots and molecular conductors. These models also appear as an auxiliary problem whose solution yields the dynamical mean-field description of correlated lattice models.
Dynamical mean-field theory (DMFT) is one of the most widely-used methods to treat accurately electron correlation effects in ab-initio real material calculations. Many modern large-scale implementations of DMFT in electronic structure codes involve solving a quantum impurity model with a continuous-time quantum Monte Carlo (CT-QMC) solver [Rubtsov et al., Phys. Rev. B 72, 035122 (2005); Werner et al., Phys. Rev. Lett. 97, 076405 (2006); Werner and Millis, Phys. Rev. B 74, 155107 (2006); Gull et al., Rev. Mod. Phys. 83, 349 (2011)]. The main advantage of CT-QMC is that, unlike standard quantum Monte Carlo approaches, it is able to generate the local Green's functions $G(\ensuremath{\tau})$ of the correlated system on an arbitrarily fine imaginary time $\ensuremath{\tau}$ grid, and is free of any systematic errors. In this work, we extend a hybrid QMC solver proposed by Khatami et al. [Phys. Rev. E 81, 056703 (2010)] and Rost et al. [Phys. Rev. E 87, 053305 (2013)] to a multiorbital context. This has the advantage of enabling impurity solver QMC calculations to scale linearly with inverse temperature $\ensuremath{\beta}$, and permit its application to ddd- and fff-band materials. In addition, we present a Green's-function processing scheme which generates accurate quasicontinuous imaginary time solutions of the impurity problem which overcome errors inherent to standard QMC approaches. This solver and processing scheme are incorporated into a full DFT+DMFT calculation using the castep DFT code [Clark et al., Z. Kristallogr. 220, 567 (2005)]. Benchmark calculations for SrVO3{\mathrm{SrVO}}_{3}SrVO3​ properties are presented. The computational efficiency of this method is also demonstrated.