Source author record

Michel Ferrero

Michel Ferrero appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2021arXiv

Fast Principal Minor Algorithms for Quantum Many Body Systems

The computation of determinants plays a central role in diagrammatic Monte Carlo algorithms for strongly correlated systems. The evaluation of large numbers of determinants can often be the limiting computational factor determining the number of attainable diagrammatic expansion orders. In this work we build upon the algorithm presented in [\emph{Linear Algebra and its applications} 419.1 (2006), 107-124] which computes all principal minors of a matrix in $O(2^n)$ operations. We present multiple generalizations of the algorithm to the efficient evaluation of certain subsets of all principal minors with immediate applications to Connected Determinant Diagrammatic Monte Carlo within the normal and symmetry-broken phases as well as Continuous-time Quantum Monte Carlo. Additionally, we improve the asymptotic scaling of diagrammatic Monte Carlo formulated in real-time to $O(2^n)$ and report speedups of up to a factor $25$ at computationally realistic expansion orders. We further show that all permanent-principal-minors, corresponding to sums of bosonic Feynman diagrams, can be computed in $O(3^n)$, making diagrammatic Monte Carlo for bosonic and mixed systems a viable path worth exploring.

preprint2021arXiv

The Weak, the Strong and the Long Correlation Regimes of the Two-Dimensional Hubbard Model at Finite Temperature

We investigate the momentum-resolved spin and charge susceptibilities, as well as the chemical potential and double occupancy in the two-dimensional Hubbard model as functions of doping, temperature and interaction strength. Through these quantities, we identify a weak-coupling regime, a strong-coupling regime with short-range correlations and an intermediate-coupling regime with long magnetic correlation lengths. In the spin channel, we observe an additional crossover from commensurate to incommensurate correlations. In contrast, we find charge correlations to be only short ranged for all studied temperatures, which suggests that the spin and charge responses are decoupled. These findings were obtained by a novel connected determinant diagrammatic Monte Carlo algorithm for the computation of double expansions, which we introduce in this paper. This permits us to obtain numerically exact results at unprecedentedly low temperatures $T\geq 0.067$ for interactions up to $U\leq 8$, while working on arbitrarily large lattices. Our method also allows us to gain physical insights from investigating the analytic structure of perturbative series. We connect to previous work by studying smaller lattice geometries and report substantial finite-size effects.

preprint2020arXiv

Real-frequency Diagrammatic Monte Carlo at Finite Temperature

Diagrammatic expansions are a central tool for treating correlated electron systems. At thermal equilibrium, they are most naturally defined within the Matsubara formalism. However, extracting any dynamic response function from a Matsubara calculation ultimately requires the ill-defined analytical continuation from the imaginary- to the real-frequency domain. It was recently proposed [Phys. Rev. B 99, 035120 (2019)] that the internal Matsubara summations of any interaction-expansion diagram can be performed analytically by using symbolic algebra algorithms. The result of the summations is then an analytical function of the complex frequency rather than Matsubara frequency. Here we apply this principle and develop a diagrammatic Monte Carlo technique which yields results directly on the real-frequency axis. We present results for the self-energy $Σ(ω)$ of the doped 32x32 cyclic square-lattice Hubbard model in a non-trivial parameter regime, where signatures of the pseudogap appear close to the antinode. We discuss the behavior of the perturbation series on the real-frequency axis and in particular show that one must be very careful when using the maximum entropy method on truncated perturbation series. Our approach holds great promise for future application in cases when analytical continuation is difficult and moderate-order perturbation theory may be sufficient to converge the result.

preprint2019arXiv

Finite temperature density matrix embedding theory

We describe a formulation of the density matrix embedding theory at finite temperature. We present a generalization of the ground-state bath orbital construction that embeds a mean-field finite-temperature density matrix up to a given order in the Hamiltonian, or the Hamiltonian up to a given order in the density matrix. We assess the performance of the finite-temperature density matrix embedding on the 1D Hubbard model both at half-filling and away from it, and the 2D Hubbard model at half-filling, comparing to exact data where available, as well as results from finite-temperature density matrix renormalization group, dynamical mean-field theory, and dynamical cluster approximations. The accuracy of finite-temperature density matrix embedding appears comparable to that of the ground-state theory, with at most a modest increase in bath size, and competitive with that of cluster dynamical mean-field theory.

preprint2019arXiv

Not all doped Mott insulators have a pseudogap: key role of van Hove singularities

The Mott insulating phase of the parent compounds is frequently taken as starting point for the underdoped high-$T_c$ cuprate superconductors. In particular, the pseudogap state is often considered as deriving from the Mott insulator. In this work, we systematically investigate different weakly-doped Mott insulators on the square and triangular lattice to clarify the relationship between the pseudogap and Mottness. We show that doping a two-dimensional Mott insulator does not necessarily lead to a pseudogap phase. Despite its inherent strong-coupling nature, we find that the existence or absence of a pseudogap depends sensitively on non-interacting band parameters and identify the crucial role played by the van Hove singularities of the system. Motivated by a SU(2) gauge theory for the pseudogap state, we propose and verify numerically a simple equation that governs the evolution of characteristic features in the electronic scattering rate.

preprint2016arXiv

Enhancing superconductivity in A3C60 fullerides

Motivated by the recent experimental report of a possible light-induced superconductivity in A3C60 at high temperature [Mitrano et al., Nature 530, 451 (2016)], we investigate theoretical mechanisms for enhanced superconductivity in A3C60 fullerenes. We find that an `interaction imbalance' corresponding to a smaller value of the Coulomb matrix element for two of the molecular orbitals in comparison to the third one, efficiently enhances superconductivity. Furthermore, we perform first-principle calculations of the changes in the electronic structure and in the screened Coulomb matrix elements of A3C60, brought in by the deformation associated with the pumped T1u intra-molecular mode. We find that an interaction imbalance is indeed induced, with a favorable sign and magnitude for superconductivity enhancement. The physical mechanism responsible for this enhancement consists in a stabilisation of the intra-molecular states containing a singlet pair, while preserving the orbital fluctuations allowing for a coherent inter-orbital delocalization of the pair. Other perturbations have also been considered and found to be detrimental to superconductivity. The light-induced deformation and ensuing interaction imbalance is shown to bring superconductivity further into the strong-coupling regime.

preprint2016arXiv

Large effects of subtle electronic correlations on the energetics of vacancies in alpha-Fe

We study the effect of electronic Coulomb correlations on the vacancy formation energy in paramagnetic alpha-Fe within ab initio dynamical mean-field theory. The calculated value for the formation energy is substantially lower than in standard density-functional calculations and in excellent agreement with experiment. The reduction is caused by an enhancement of electronic correlations at the nearest neighbors of the vacancy. This effect is explained by subtle changes in the corresponding spectral function of the d-electrons. The local lattice relaxations around the vacancy are substantially increased by many-body effects.

preprint2016arXiv

TRIQS/DFTTools: A TRIQS application for ab initio calculations of correlated materials

We present the TRIQS/DFTTools package, an application based on the TRIQS library that connects this toolbox to realistic materials calculations based on density functional theory (DFT). In particular, TRIQS/DFTTools together with TRIQS allows an efficient implementation of DFT plus dynamical mean-field theory (DMFT) calculations. It supplies tools and methods to construct Wannier functions and to perform the DMFT self-consistency cycle in this basis set. Post-processing tools, such as band-structure plotting or the calculation of transport properties are also implemented. The package comes with a fully charge self-consistent interface to the Wien2k band structure code, as well as a generic interface that allows to use TRIQS/DFTTools together with a large variety of DFT codes. It is distributed under the GNU General Public License (GPLv3).

preprint2015arXiv

Non-existence of the Luttinger-Ward functional and misleading convergence of skeleton diagrammatic series for Hubbard-like models

The Luttinger-Ward functional $Φ[\mathbf{G}]$, which expresses the thermodynamic grand potential in terms of the interacting single-particle Green's function $\mathbf{G}$, is found to be ill-defined for fermionic models with the Hubbard on-site interaction. In particular, we show that the self-energy $\mathbfΣ[\mathbf{G}] \propto δΦ[\mathbf{G}]/δ\mathbf{G}$ is not a single-valued functional of $\mathbf{G}$: in addition to the physical solution for $\mathbfΣ[\mathbf{G}]$, there exists at least one qualitatively distinct unphysical branch. This result is demonstrated for several models: the Hubbard atom, the Anderson impurity model, and the full two-dimensional Hubbard model. Despite this pathology, the skeleton Feynman diagrammatic series for $\mathbfΣ$ in terms of $\mathbf{G}$ is found to converge at least for moderately low temperatures. However, at strong interactions, its convergence is to the unphysical branch. This reveals a new scenario of breaking down of diagrammatic expansions. In contrast, the bare series in terms of the non-interacting Green's function $\mathbf{G}_0$ converges to the correct physical branch of $\mathbfΣ$ in all cases currently accessible by diagrammatic Monte Carlo. Besides their conceptual importance, these observations have important implications for techniques based on the explicit summation of diagrammatic series.

preprint2015arXiv

Solutions of the Two Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms

Numerical results for ground state and excited state properties (energies, double occupancies, and Matsubara-axis self energies) of the single-orbital Hubbard model on a two-dimensional square lattice are presented, in order to provide an assessment of our ability to compute accurate results in the thermodynamic limit. Many methods are employed, including auxiliary field quantum Monte Carlo, bare and bold-line diagrammatic Monte Carlo, method of dual fermions, density matrix embedding theory, density matrix renormalization group, dynamical cluster approximation, diffusion Monte Carlo within a fixed node approximation, unrestricted coupled cluster theory, and multi-reference projected Hartree-Fock. Comparison of results obtained by different methods allows for the identification of uncertainties and systematic errors. The importance of extrapolation to converged thermodynamic limit values is emphasized. Cases where agreement between different methods is obtained establish benchmark results that may be useful in the validation of new approaches and the improvement of existing methods.

preprint2015arXiv

TRIQS: A Toolbox for Research on Interacting Quantum Systems

We present the TRIQS library, a Toolbox for Research on Interacting Quantum Systems. It is an open-source, computational physics library providing a framework for the quick development of applications in the field of many-body quantum physics, and in particular, strongly-correlated electronic systems. It supplies components to develop codes in a modern, concise and efficient way: e.g. Green's function containers, a generic Monte Carlo class, and simple interfaces to HDF5. TRIQS is a C++/Python library that can be used from either language. It is distributed under the GNU General Public License (GPLv3). State-of-the-art applications based on the library, such as modern quantum many-body solvers and interfaces between density-functional-theory codes and dynamical mean-field theory (DMFT) codes are distributed along with it.

preprint2015arXiv

TRIQS/CTHYB: A Continuous-Time Quantum Monte Carlo Hybridization Expansion Solver for Quantum Impurity Problems

We present TRIQS/CTHYB, a state-of-the art open-source implementation of the continuous-time hybridisation expansion quantum impurity solver of the TRIQS package. This code is mainly designed to be used with the TRIQS library in order to solve the self-consistent quantum impurity problem in a multi-orbital dynamical mean field theory approach to strongly-correlated electrons, in particular in the context of realistic calculations. It is implemented in C++ for efficiency and is provided with a high-level Python interface. The code is ships with a new partitioning algorithm that divides the local Hilbert space without any user knowledge of the symmetries and quantum numbers of the Hamiltonian. Furthermore, we implement higher-order configuration moves and show that such moves are necessary to ensure ergodicity of the Monte Carlo in common Hamiltonians even without symmetry-breaking.

preprint2014arXiv

Dynamical correlations and screened exchange on the experimental bench: spectral properties of the cobalt pnictide BaCo2As2

Understanding the Fermi surface and low-energy excitations of iron or cobalt pnictides is crucial for assessing electronic instabilities such as magnetic or superconducting states. Here, we propose and implement a new approach to compute the low-energy properties of correlated electron materials, taking into account both screened exchange beyond the local density approximation and local dynamical correlations. The scheme allows us to resolve the puzzle of BaCo2As2, for which standard electronic structure techniques predict a ferromagnetic instability not observed in nature.

preprint2014arXiv

Orbital Polarization in Strained LaNiO$_{3}$: Structural Distortions and Correlation Effects

Transition-metal heterostructures offer the fascinating possibility of controlling orbital degrees of freedom via strain. Here, we investigate theoretically the degree of orbital polarization that can be induced by epitaxial strain in LaNiO$_3$ films. Using combined electronic structure and dynamical mean-field theory methods we take into account both structural distortions and electron correlations and discuss their relative influence. We confirm that Hund's rule coupling tends to decrease the polarization and point out that this applies to both the $d^8\underline{L}$ and $d^7$ local configurations of the Ni ions. Our calculations are in good agreement with recent experiments, which revealed sizable orbital polarization under tensile strain. We discuss why full orbital polarization is hard to achieve in this specific system and emphasize the general limitations that must be overcome to achieve this goal.

preprint2013arXiv

How bad metals turn good: spectroscopic signatures of resilient quasiparticles

We investigate transport in strongly correlated metals. Within dynamical mean-field theory, we calculate the resistivity, thermopower, optical conductivity and thermodynamic properties of a hole-doped Mott insulator. Two well-separated temperature scales are identified: T_FL below which Landau Fermi liquid behavior applies, and T_MIR above which the resistivity exceeds the Mott-Ioffe-Regel value and `bad-metal' behavior is found. We show that quasiparticle excitations remain well-defined above T_FL and dominate transport throughout the intermediate regime T_FL < T_MIR. The lifetime of these `resilient quasiparticles' is longer for electron-like excitations, and this pronounced particle-hole asymmetry has important consequences for the thermopower. The crossover into the bad-metal regime corresponds to the disappearance of these excitations, and has clear signatures in optical spectroscopy.

preprint2013arXiv

Theoretical prediction and spectroscopic fingerprints of an orbital transition in CeCu2Si2

We show that the heavy-fermion compound CeCu2Si2 undergoes a transition between two regimes dominated by different crystal-field states. At low pressure P and low temperature T the Ce 4f electron resides in the atomic crystal-field ground state, while at high P or T the electron occupancy and spectral weight is transferred to an excited crystal-field level that hybridizes more strongly with itinerant states. These findings result from first-principles dynamical-mean-field-theory calculations. We predict experimental signatures of this orbital transition in X-ray spectroscopy. The corresponding fluctuations may be responsible for the second high-pressure superconducting dome observed in this and similar materials.

preprint2011arXiv

Hallmark of strong electronic correlations in LaNiO$_3$: photoemission kink and broadening of fully occupied bands

Recent angular-resolved photoemission experiments on LaNiO$_3$ reported a renormalization of the Fermi velocity of $e_g$ quasiparticles, a kink in their dispersion at $ -0.2$ eV and a large broadening and weakened dispersion of the occupied $t_{2g}$ states. We show here that all these features result from electronic correlations and are quantitatively reproduced by calculations combining density-functional theory and dynamical mean-field theory. The importance and general relevance of correlation effects in filled bands coupled by inter-orbital interactions to a partially-filled band are pointed out.

preprint2011arXiv

Orthogonal Polynomial Representation of Imaginary-Time Green's Functions

We study the expansion of single-particle and two-particle imaginary-time Matsubara Green's functions of quantum impurity models in the basis of Legendre orthogonal polynomials. We discuss various applications within the dynamical mean-field theory (DMFT) framework. The method provides a more compact representation of the Green's functions than standard Matsubara frequencies and therefore significantly reduces the memory-storage size of these quantities. Moreover, it can be used as an efficient noise filter for various physical quantities within the continuous-time quantum Monte Carlo impurity solvers recently developed for DMFT and its extensions. In particular, we show how to use it for the computation of energies in the context of realistic DMFT calculations in combination with the local density approximation to the density functional theory (LDA+DMFT) and for the calculation of lattice susceptibilities from the local irreducible vertex function.

preprint2010arXiv

Momentum space anisotropy and pseudogaps: a comparative cluster dynamical mean field analysis of the doping-driven metal-insulator transition in the two dimensional Hubbard model

Cluster dynamical mean field calculations based on 2, 4, 8 and 16 site clusters are used to analyze the doping-driven metal-insulator transition in the two dimensional Hubbard model. Comparison of results obtained on different clusters enables a determination of those aspects of the physics that are common to all clusters and permits identification of artifacts associated with particular cluster geometries. A modest particle-hole asymmetry in the underlying band structure is shown to lead to qualitatively different behavior on the hole doped side than on the electron doped side. For particle-hole asymmetry of the sign and magnitude appropriate to high-$T_c$ cuprates, the approach to the insulator from the hole-doping side is found to proceed in two stages from a high-doping region where the properties are those of a Fermi liquid with moderately renormalized parameters and very weak momentum dependence. As doping is reduced the system first enters an intermediate doping regime where the Fermi liquid renormalizations are larger and the electron self energy varies significantly around the Fermi surface and then passes to a small doping regime characterized by a gap in some regions of the Fermi surface but gapless behavior in other regions. On the electron doped side the partially gapped regime does not occur, and the momentum dependence of the electron self energy is less pronounced. Implications for the high-$T_c$ cuprates and for the use of cluster dynamical mean field methods in wider classes of problems are discussed.