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

arXiv:1104.3854 · Physical Review B

Dual fermion dynamical cluster approach for strongly correlated systems

Shuxiang Yang, Herbert F. Fotso, Hartmut Hafermann, Ka-Ming Tam, Juana Moreno, Thomas Pruschke, Mark Jarrell

Published Oct 11, 2011• 30 citations
arXivHTMLPDF

Abstract

We have designed a multiscale approach for strongly correlated systems by combining the dynamical cluster approximation (DCA) and the recently introduced dual fermion formalism. This approach employs an exact mapping from a real lattice to a DCA cluster of linear size ${L}_{c}$ embedded in a dual fermion lattice. Short-length-scale physics is addressed by the DCA cluster calculation, while longer-length-scale physics is addressed diagrammatically using dual fermions. The bare and dressed dual fermionic Green functions scale as $\mathcal{O}(1/{L}_{c})$, so perturbation theory on the dual lattice converges very quickly, e.g., the dual Fermion self-energy calculated with simple second-order perturbation theory is of order $\mathcal{O}(1/{L}_{c}^{3})$ with third-order and three-body corrections down by an additional factor of $\mathcal{O}(1/{L}_{c})$.

DOI
10.1103/physrevb.84.155106

Related papers

Momentum-space cluster dual-fermion method

Sergei Iskakov, Hanna Terletska, Emanuel Gull

Recent years have seen the development of two types of nonlocal extensions to the single-site dynamical mean field theory. On one hand, cluster approximations, such as the dynamical cluster approximation, recover short-range momentum-dependent correlations nonperturbatively. On the other hand, diagrammatic extensions, such as the dual-fermion theory, recover long-ranged corrections perturbatively. The correct treatment of both strong short-ranged and weak long-ranged correlations within the same framework is therefore expected to lead to a quick convergence of results, and offers the potential of obtaining smooth self-energies in nonperturbative regimes of phase space. In this paper, we present an exact cluster dual-fermion method based on an expansion around the dynamical cluster approximation. Unlike previous formulations, our method does not employ a coarse-graining approximation to the interaction, which we show to be the leading source of error at high temperature, and converges to the exact result independently of the size of the underlying cluster. We illustrate the power of the method with results for the second-order cluster dual-fermion approximation to the single-particle self-energies and double occupancies.

arXiv:1712.03295•Physical review. B./Physical review. B•Mar 12, 2018•14 citations

Next generation multi-scale quantum simulations for strongly correlated materials

Shuxiang Yang

This thesis represents our effort to develop the next generation multi-scale quantum simulation methods suitable for strongly-correlated materials, where complicated phase-diagrams prevail, suggesting complicated underlying physics. We first give a detailed description of the parquet formalism. With its help, different approximate methods can be unified and a hierarchy of approximate methods with different accuracies and computational complexity can thus be designed. Next, we present a numerical solution of the parquet approximation. Results on the Hubbard model are compared to those obtained from Determinant Quantum Monte Carlo (DQMC), FLuctuation EXchange (FLEX), and self-consistent second-order approximation methods. The comparison shows a satisfactory agreement with DQMC and a significant improvement over the FLEX or the self-consistent second-order approximation. The parquet formalism can also be used to analyze the superconducting mechanism of the high-temperature superconductors. The dynamical cluster approximation (DCA) method is used to understand the proximity of the superconducting dome to the quantum critical point in the 2-D Hubbard model. At optimal doping, where Vd is revealed to be featureless, we find a power-law behavior of chi_0d(w=0), replacing the BCS logarithm behavior, and strongly enhanced T_c. After that we propose another multi-scale approach by combining the DCA and the recently introduced dual-fermion formalism. Within this approach, short and long length scale physics is addressed by the DCA cluster calculation, while intermediate length scale physics is addressed diagrammatically using dual fermions. The bare and dressed dual fermionic Green functions scale as O(1/Lc), so perturbation theory on the dual lattice converges very quickly. Lastly, we study the responses to dynamical modulation of the optical lattice potential by analyzing properties of the repulsive fermionic Hubbard model in an optical lattice. We provide numerical evidence showing the modulations by on-site local interaction cannot be ignored, and instead can even strongly contribute to the dynamical behaviors of the system in highly-doped cases.

Civil War Book Review•Aug 22, 2011•2 citations

Beyond quantum cluster theories: multiscale approaches for strongly correlated systems

Herbert F. Fotso, Ka-Ming Tam, Juana Moreno

Abstract The degrees of freedom that confer to strongly correlated systems their many intriguing properties also render them fairly intractable through typical perturbative treatments. For this reason, the mechanisms responsible for their technologically promising properties remain mostly elusive. Computational approaches have played a major role in efforts to fill this void. In particular, dynamical mean field theory and its cluster extension, the dynamical cluster approximation have allowed significant progress. However, despite all the insightful results of these embedding schemes, computational constraints, such as the minus sign problem in quantum Monte Carlo (QMC), and the exponential growth of the Hilbert space in exact diagonalization (ED) methods, still limit the length scale within which correlations can be treated exactly in the formalism. A recent advance aiming to overcome these difficulties is the development of multiscale many body approaches whereby this challenge is addressed by introducing an intermediate length scale between the short length scale where correlations are treated exactly using a cluster solver such QMC or ED, and the long length scale where correlations are treated in a mean field manner. At this intermediate length scale correlations can be treated perturbatively. This is the essence of multiscale many-body methods. We will review various implementations of these multiscale many-body approaches, the results they have produced, and the outstanding challenges that should be addressed for further advances.

arXiv:2011.05522•Quantum Science and Technology•Apr 14, 2022•6 citations

Dynamical cluster approximation with continuous lattice self-energy

Peter Staar, Thomas Maier, T. C. Schulthess

The dynamical cluster approximation (DCA) is a systematic extension beyond the single-site approximation in dynamical mean field theory, to include spatially nonlocal correlations in quantum many-body simulations of strongly correlated systems. We extend the DCA with a continuous lattice self-energy in order to achieve better convergence with cluster size. This method, which we call DCA${}^{+}$, cures the cluster-shape dependence problems of the DCA, without suffering from causality violations of previous attempts to interpolate the cluster self-energy. A practical approach based on standard inference techniques is given to deduce the continuous lattice self-energy from an interpolated cluster self-energy. We study the pseudogap region of a hole-doped two-dimensional Hubbard model and find that, in the DCA${}^{+}\phantom{\rule{0.222222em}{0ex}}$algorithm, the self-energy and pseudogap temperature ${T}^{*}$ converge monotonously with cluster size. Introduction of a continuous lattice self-energy eliminates artificial long-range correlations and thus significantly reduces the sign problem of the quantum Monte Carlo cluster solver in the DCA${}^{+}\phantom{\rule{0.222222em}{0ex}}$algorithm compared to the normal DCA. Simulations with much larger cluster sizes thus become feasible, which, along with the improved convergence in cluster size, raises hope that precise extrapolations to the exact infinite cluster size limit can be reached for other physical quantities as well.

arXiv:1304.3624•Physical Review B•Sep 3, 2013•36 citations

Interlaced coarse-graining for the dynamic cluster approximation

Peter Staar, Mi Jiang, Urs R. Hähner +2 more

The dynamical cluster approximation (DCA) and its ${\mathrm{DCA}}^{+}$ extension use coarse-graining of the momentum space to reduce the complexity of quantum many-body problems, thereby mapping the bulk lattice to a cluster embedded in a dynamical mean-field host. Here, we introduce a new form of an interlaced coarse-graining and compare it with the traditional coarse-graining. While it gives a more localized self-energy for a given cluster size, we show that it leads to more controlled results with weaker cluster shape and smoother cluster size dependence, which converge to the results obtained from the standard coarse-graining with increasing cluster size. Most importantly, the new coarse-graining reduces the severity of the fermionic sign problem of the underlying quantum Monte Carlo cluster solver and thus allows for calculations on larger clusters. This enables the treatment of correlations longer ranged than those accessible with the standard coarse-graining and thus can allow for the evaluation of the exact infinite cluster size result via finite size scaling. As a demonstration, we study the hole-doped two-dimensional Hubbard model and show that the interlaced coarse-graining in combination with the extended ${\mathrm{DCA}}^{+}$ algorithm permits the determination of the superconducting ${T}_{c}$ on cluster sizes for which the results can be fit with a Kosterlitz-Thouless scaling law.

arXiv:1601.03838•Physical review. B./Physical review. B•Apr 27, 2016•12 citations

Dual Fermion Approach to Disordered Correlated Systems

Patrick Haase

Disorder is ubiquitous in real materials and influences the physical properties like the conductivity to varying degrees. \nIf electron-electron interactions are strong, theoretical and numerical treatment of these systems becomes challenging. In this thesis a numerical approach is developed to address these systems, \ntreating both interactions and disorder on equal footing. \nThe approach is based on the dual fermion approach for interacting systems developed by Rubtsov et al. \nTerletska et al. applied the ideas of the dual fermion approach to disordered non-interacting systems. In this approach, the replica trick is used to integrate out the disorder in favor of an effective electron-electron interaction. \nWe extended the approach from Terletska et al. to treat disordered interacting systems. \n \nDual Fermions allow to take into account non-local fluctuations by means of a perturbative expansion around an impurity problem. The impurity reference system is determined self-consistently, analogously to the dynamical mean-field \ntheory. \nThe perturbative expansion is expected to yield good results for small and large values of interaction strength and disorder. A priori, it is not clear what to expect for intermediate values, \nbut experience shows that oftentimes good results are obtained for this region. \n \nAn advantage of the dual fermion approach is that there is no sign-problem for a single orbital model if quantum Monte Carlo is used to solve the interacting reference system. \nAdditionally, perturbation theory is usually numerically much cheaper than fully solving an interacting lattice or cluster problem. \nThus, the dual fermion approach allows to address regions of parameter space that are not accessible to lattice quantum Monte Carlo calculations or cluster extension of dynamical mean-field theory. \nCluster extensions of the dynamical mean-field theory are for example the dynamical cluster approximation or the cellular dynamical mean-field theory. \n \nThe new approach is benchmarked with respect to the Anderson-Falicov-Kimball model. The Anderson-Hubbard model is the main application. \nFor both models, the algorithm is tested in one dimension for the single particle Green function and compared to the dynamical cluster approximation to estimate the quality of the results. \nThe Anderson-Falicov-Kimball model is studied further in two dimensions and the phase diagram on the disorder-interaction plane is calculated from both the single particle Green function and the dc conductivity. \n \nFor the Anderson-Hubbard model, the antiferromagnetic phase transition, the Mott transition and the Anderson transition are investigated in three dimensions. The dynamical mean-field theory results are \nqualitatively confirmed for the most part. Quantitative corrections are in agreement with dynamical cluster approximation results, thus confirming the validity of the results. \nWe found a qualitative difference between the dual fermion approach and dynamical mean-field theory for the temperature dependence of the hysteresis of the double occupancy. The dual fermion result \nshows a decreasing lower critical interaction strength for decreasing temperature, which is in contrast with the dynamical mean-field result that shows an increase of the lower critical interaction strength. \nThis behavior happens for small temperatures that are not accessible using for example the dynamical cluster approximation, thus highlighting the usefulness of the dual fermion approach. \nWe found this behavior for all values of the disorder that we considered.

•Nov 5, 2015•0 citations
We have designed a multiscale approach for strongly correlated systems by combining the dynamical cluster approximation (DCA) and the recently introduced dual fermion formalism. This approach employs an exact mapping from a real lattice to a DCA cluster of linear size Lc{L}_{c}Lc​ embedded in a dual fermion lattice. Short-length-scale physics is addressed by the DCA cluster calculation, while longer-length-scale physics is addressed diagrammatically using dual fermions. The bare and dressed dual fermionic Green functions scale as O(1/Lc)\mathcal{O}(1/{L}_{c})O(1/Lc​), so perturbation theory on the dual lattice converges very quickly, e.g., the dual Fermion self-energy calculated with simple second-order perturbation theory is of order O(1/Lc3)\mathcal{O}(1/{L}_{c}^{3})O(1/Lc3​) with third-order and three-body corrections down by an additional factor of O(1/Lc)\mathcal{O}(1/{L}_{c})O(1/Lc​).
Recent years have seen the development of two types of nonlocal extensions to the single-site dynamical mean field theory. On one hand, cluster approximations, such as the dynamical cluster approximation, recover short-range momentum-dependent correlations nonperturbatively. On the other hand, diagrammatic extensions, such as the dual-fermion theory, recover long-ranged corrections perturbatively. The correct treatment of both strong short-ranged and weak long-ranged correlations within the same framework is therefore expected to lead to a quick convergence of results, and offers the potential of obtaining smooth self-energies in nonperturbative regimes of phase space. In this paper, we present an exact cluster dual-fermion method based on an expansion around the dynamical cluster approximation. Unlike previous formulations, our method does not employ a coarse-graining approximation to the interaction, which we show to be the leading source of error at high temperature, and converges to the exact result independently of the size of the underlying cluster. We illustrate the power of the method with results for the second-order cluster dual-fermion approximation to the single-particle self-energies and double occupancies.
This thesis represents our effort to develop the next generation multi-scale quantum simulation methods suitable for strongly-correlated materials, where complicated phase-diagrams prevail, suggesting complicated underlying physics. We first give a detailed description of the parquet formalism. With its help, different approximate methods can be unified and a hierarchy of approximate methods with different accuracies and computational complexity can thus be designed. Next, we present a numerical solution of the parquet approximation. Results on the Hubbard model are compared to those obtained from Determinant Quantum Monte Carlo (DQMC), FLuctuation EXchange (FLEX), and self-consistent second-order approximation methods. The comparison shows a satisfactory agreement with DQMC and a significant improvement over the FLEX or the self-consistent second-order approximation. The parquet formalism can also be used to analyze the superconducting mechanism of the high-temperature superconductors. The dynamical cluster approximation (DCA) method is used to understand the proximity of the superconducting dome to the quantum critical point in the 2-D Hubbard model. At optimal doping, where Vd is revealed to be featureless, we find a power-law behavior of chi_0d(w=0), replacing the BCS logarithm behavior, and strongly enhanced T_c. After that we propose another multi-scale approach by combining the DCA and the recently introduced dual-fermion formalism. Within this approach, short and long length scale physics is addressed by the DCA cluster calculation, while intermediate length scale physics is addressed diagrammatically using dual fermions. The bare and dressed dual fermionic Green functions scale as O(1/Lc), so perturbation theory on the dual lattice converges very quickly. Lastly, we study the responses to dynamical modulation of the optical lattice potential by analyzing properties of the repulsive fermionic Hubbard model in an optical lattice. We provide numerical evidence showing the modulations by on-site local interaction cannot be ignored, and instead can even strongly contribute to the dynamical behaviors of the system in highly-doped cases.
Abstract The degrees of freedom that confer to strongly correlated systems their many intriguing properties also render them fairly intractable through typical perturbative treatments. For this reason, the mechanisms responsible for their technologically promising properties remain mostly elusive. Computational approaches have played a major role in efforts to fill this void. In particular, dynamical mean field theory and its cluster extension, the dynamical cluster approximation have allowed significant progress. However, despite all the insightful results of these embedding schemes, computational constraints, such as the minus sign problem in quantum Monte Carlo (QMC), and the exponential growth of the Hilbert space in exact diagonalization (ED) methods, still limit the length scale within which correlations can be treated exactly in the formalism. A recent advance aiming to overcome these difficulties is the development of multiscale many body approaches whereby this challenge is addressed by introducing an intermediate length scale between the short length scale where correlations are treated exactly using a cluster solver such QMC or ED, and the long length scale where correlations are treated in a mean field manner. At this intermediate length scale correlations can be treated perturbatively. This is the essence of multiscale many-body methods. We will review various implementations of these multiscale many-body approaches, the results they have produced, and the outstanding challenges that should be addressed for further advances.
The dynamical cluster approximation (DCA) is a systematic extension beyond the single-site approximation in dynamical mean field theory, to include spatially nonlocal correlations in quantum many-body simulations of strongly correlated systems. We extend the DCA with a continuous lattice self-energy in order to achieve better convergence with cluster size. This method, which we call DCA+{}^{+}+, cures the cluster-shape dependence problems of the DCA, without suffering from causality violations of previous attempts to interpolate the cluster self-energy. A practical approach based on standard inference techniques is given to deduce the continuous lattice self-energy from an interpolated cluster self-energy. We study the pseudogap region of a hole-doped two-dimensional Hubbard model and find that, in the DCA+{}^{+}\phantom{\rule{0.222222em}{0ex}}+algorithm, the self-energy and pseudogap temperature T∗{T}^{*}T∗ converge monotonously with cluster size. Introduction of a continuous lattice self-energy eliminates artificial long-range correlations and thus significantly reduces the sign problem of the quantum Monte Carlo cluster solver in the DCA+{}^{+}\phantom{\rule{0.222222em}{0ex}}+algorithm compared to the normal DCA. Simulations with much larger cluster sizes thus become feasible, which, along with the improved convergence in cluster size, raises hope that precise extrapolations to the exact infinite cluster size limit can be reached for other physical quantities as well.
The dynamical cluster approximation (DCA) and its DCA+{\mathrm{DCA}}^{+}DCA+ extension use coarse-graining of the momentum space to reduce the complexity of quantum many-body problems, thereby mapping the bulk lattice to a cluster embedded in a dynamical mean-field host. Here, we introduce a new form of an interlaced coarse-graining and compare it with the traditional coarse-graining. While it gives a more localized self-energy for a given cluster size, we show that it leads to more controlled results with weaker cluster shape and smoother cluster size dependence, which converge to the results obtained from the standard coarse-graining with increasing cluster size. Most importantly, the new coarse-graining reduces the severity of the fermionic sign problem of the underlying quantum Monte Carlo cluster solver and thus allows for calculations on larger clusters. This enables the treatment of correlations longer ranged than those accessible with the standard coarse-graining and thus can allow for the evaluation of the exact infinite cluster size result via finite size scaling. As a demonstration, we study the hole-doped two-dimensional Hubbard model and show that the interlaced coarse-graining in combination with the extended DCA+{\mathrm{DCA}}^{+}DCA+ algorithm permits the determination of the superconducting Tc{T}_{c}Tc​ on cluster sizes for which the results can be fit with a Kosterlitz-Thouless scaling law.
Disorder is ubiquitous in real materials and influences the physical properties like the conductivity to varying degrees. \nIf electron-electron interactions are strong, theoretical and numerical treatment of these systems becomes challenging. In this thesis a numerical approach is developed to address these systems, \ntreating both interactions and disorder on equal footing. \nThe approach is based on the dual fermion approach for interacting systems developed by Rubtsov et al. \nTerletska et al. applied the ideas of the dual fermion approach to disordered non-interacting systems. In this approach, the replica trick is used to integrate out the disorder in favor of an effective electron-electron interaction. \nWe extended the approach from Terletska et al. to treat disordered interacting systems. \n \nDual Fermions allow to take into account non-local fluctuations by means of a perturbative expansion around an impurity problem. The impurity reference system is determined self-consistently, analogously to the dynamical mean-field \ntheory. \nThe perturbative expansion is expected to yield good results for small and large values of interaction strength and disorder. A priori, it is not clear what to expect for intermediate values, \nbut experience shows that oftentimes good results are obtained for this region. \n \nAn advantage of the dual fermion approach is that there is no sign-problem for a single orbital model if quantum Monte Carlo is used to solve the interacting reference system. \nAdditionally, perturbation theory is usually numerically much cheaper than fully solving an interacting lattice or cluster problem. \nThus, the dual fermion approach allows to address regions of parameter space that are not accessible to lattice quantum Monte Carlo calculations or cluster extension of dynamical mean-field theory. \nCluster extensions of the dynamical mean-field theory are for example the dynamical cluster approximation or the cellular dynamical mean-field theory. \n \nThe new approach is benchmarked with respect to the Anderson-Falicov-Kimball model. The Anderson-Hubbard model is the main application. \nFor both models, the algorithm is tested in one dimension for the single particle Green function and compared to the dynamical cluster approximation to estimate the quality of the results. \nThe Anderson-Falicov-Kimball model is studied further in two dimensions and the phase diagram on the disorder-interaction plane is calculated from both the single particle Green function and the dc conductivity. \n \nFor the Anderson-Hubbard model, the antiferromagnetic phase transition, the Mott transition and the Anderson transition are investigated in three dimensions. The dynamical mean-field theory results are \nqualitatively confirmed for the most part. Quantitative corrections are in agreement with dynamical cluster approximation results, thus confirming the validity of the results. \nWe found a qualitative difference between the dual fermion approach and dynamical mean-field theory for the temperature dependence of the hysteresis of the double occupancy. The dual fermion result \nshows a decreasing lower critical interaction strength for decreasing temperature, which is in contrast with the dynamical mean-field result that shows an increase of the lower critical interaction strength. \nThis behavior happens for small temperatures that are not accessible using for example the dynamical cluster approximation, thus highlighting the usefulness of the dual fermion approach. \nWe found this behavior for all values of the disorder that we considered.