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 atta
inable 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.
We present a technique that enables the evaluation of perturbative expansions based on one-loop-renormalized vertices up to large expansion orders. Specifically, we show how to compute large-order corrections to the random phase approximation in eith
er the particle-hole or particle-particle channels. The algorithms efficiency is achieved by the summation over contributions of all symmetrized Feynman diagram topologies using determinants, and by integrating out analytically the two-body long-range interactions in order to yield an effective zero-range interaction. Notably, the exponential scaling of the algorithm as a function of perturbation order leads to a polynomial scaling of the approximation error with computational time for a convergent series. To assess the performance of our approach, we apply it to the non-perturbative regime of the square-lattice fermionic Hubbard model away from half-filling and report, as compared to the bare interaction expansion algorithm, significant improvements of the Monte Carlo variance as well as the convergence properties of the resulting perturbative series.
We present a general formalism that allows for the computation of large-order renormalized expansions in the spacetime representation, effectively doubling the numerically attainable perturbation order of renormalized Feynman diagrams. We show that t
his formulation compares advantageously to the currently standard techniques due to its high efficiency, simplicity, and broad range of applicability. Our formalism permits to easily complement perturbation theory with non-perturbative information, which we illustrate by implementing expansions renormalized by the addition of a gap or the inclusion of Dynamical Mean-Field Theory. As a result, we present numerically-exact results for the square-lattice Fermi-Hubbard model in the low temperature non-Fermi-liquid regime and show the momentum-dependent suppression of fermionic excitations in the antinodal region.
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, w
e 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.
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 calculatio
n 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 $Sigma(omega)$ 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.
We express the recently introduced real-time diagrammatic Quantum Monte Carlo, Phys. Rev. B 91, 245154 (2015), in the Larkin-Ovchinnikov basis in Keldysh space. Based on a perturbation expansion in the local interaction $U$, the special form of the i
nteraction vertex allows to write diagrammatic rules in which vacuum Feynman diagrams directly vanish. This reproduces the main property of the previous algorithm, without the cost of the exponential sum over Keldysh indices. In an importance sampling procedure, this implies that only interaction times in the vicinity of the measurement time contribute. Such an algorithm can then directly address the long-time limit needed in the study of steady states in out-of-equilibrium systems. We then implement and discuss different variants of Monte Carlo algorithms in the Larkin-Ovchinnikov basis. A sign problem reappears, showing that the cancellation of vacuum diagrams has no direct impact on it.
We present a theoretical investigation of the effects of correlations on the electronic structure of the Mott insulator Sr$_2$IrO$_4$ upon electron doping. A rapid collapse of the Mott gap upon doping is found, and the electronic structure displays a
strong momentum-space differentiation at low doping level: The Fermi surface consists of pockets centered around $(pi/2,pi/2)$, while a pseudogap opens near $(pi,0)$. Its physical origin is shown to be related to short-range spin correlations. The pseudogap closes upon increasing doping, but a differentiated regime characterized by a modulation of the spectral intensity along the Fermi surface persists to higher doping levels. These results, obtained within the cellular dynamical mean-field theory framework, are discussed in comparison to recent photoemission experiments and an overall good agreement is found.
We introduce and compare three different Monte Carlo determinantal algorithms that allow one to compute dynamical quantities, such as the self-energy, of fermionic systems in their thermodynamic limit. We show that the most efficient approach express
es the sum of a factorial number of one-particle-irreducible diagrams as a recursive sum of determinants with exponential complexity. By comparing results for the two-dimensional Hubbard model with those obtained from state-of-the-art diagrammatic Monte Carlo, we show that we can reach higher perturbation orders and greater accuracy for the same computational effort.
We investigate the interplay of spin-orbit coupling (SOC) and electronic correlations in Sr2RuO4 using dynamical mean-field theory. We find that SOC does not affect the correlation-induced renormalizations, which validates the Hunds metal picture of
ruthenates even in the presence of the sizable SOC relevant to these materials. Nonetheless, SOC found to change significantly the electronic structure at k-points where a degeneracy applies in its absence. We explain why these two observations are consistent with one another and calculate effects of SOC on the correlated electronic structure. The magnitude of these effects is found to depend on the energy of the quasiparticle state under consideration, leading us to introduce the notion of an energy-dependent quasiparticle spin-orbit coupling. This notion is generally applicable to all materials in which both the spin-orbit coupling and electronic correlations are sizable.
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.
