Ultrafast irradiation of correlated electronic systems triggers complex dynamics involving quasi-particle excitations, doublons, charge carriers, and spin fluctuations. To describe these effects, we develop an efficient non-equilibrium approach, dubbed D-$GW$, that enables a self-consistent treatment of local correlations within dynamical mean-field theory (DMFT) and spatial charge and spin fluctuations, that are accounted for simultaneously within a diagrammatic framework. The method is formulated in the real-time domain and provides direct access to single- and two-particle momentum- and energy-dependent response functions without the need for analytical continuation, which is required in Matsubara frequency-based approaches. We apply the D-$GW$ method to investigate the dynamics of a photo-excited extended Hubbard model, the minimal system that simultaneously hosts strong charge and spin fluctuations. Focusing on the challenging parameter regime near the Mott transition, we demonstrate that correlated metals and narrow-gap Mott insulators undergo distinct thermalization processes involving complex energy transfer between single-particle and collective electronic excitations.
Ultrafast irradiation of correlated electronic systems triggers complex dynamics involving quasi-particle excitations, doublons, charge carriers, and spin fluctuations. To describe these effects, we develop an efficient non-equilibrium approach, dubbed D-GW, that enables a self-consistent treatment of local correlations within dynamical mean-field theory (DMFT) and spatial charge and spin fluctuations, that are accounted for simultaneously within a diagrammatic framework. The method is formulated in the real-time domain and provides direct access to single- and two-particle momentum- and energy-dependent response functions without the need for analytical continuation, which is required in Matsubara frequency-based approaches. We apply the D-GW method to investigate the dynamics of a photo-excited extended Hubbard model, the minimal system that simultaneously hosts strong charge and spin fluctuations. Focusing on the challenging parameter regime near the Mott transition, we demonstrate that correlated metals and narrow-gap Mott insulators undergo distinct thermalization processes involving complex energy transfer between single-particle and collective electronic excitations.
Recent advances in time-resolved spectroscopies have enabled direct access to the momentum-selective nonequilibrium dynamics of correlated quantum materials, revealing a strongly momentum-dependent response of electrons and collective excitations. Interpreting these observations requires a real-time theoretical framework that consistently captures the interplay between strong local electronic correlations and nonlocal collective fluctuations, a capability that remains beyond state-of-the-art nonequilibrium approaches. Using a recently developed real-time many-body framework, we resolve the momentum-selective ultrafast dynamics of a photoexcited correlated electron system. We predict a transient nodal-antinodal anisotropy in electronic heating, providing a microscopic explanation for the momentum-dependent response debated in time-resolved photoemission and Raman experiments, and identify the nonthermal spectral-weight transfer responsible for the transient antinodal in-gap states observed in ultrafast photoemission. We further uncover a momentum-selective magnetic response, in which antiferromagnetic fluctuations undergo a strongly nonthermal, quench-like excitation far above the electronic temperature while preserving their correlation length, before relaxing through a momentum-space magnon cascade toward lower-momentum modes. Finally, by tracking the real-time local spin susceptibility, we identify a dynamical, experimentally accessible signature of local-moment formation and its photoinduced melting. Our results establish a unified microscopic picture of ultrafast electronic and magnetic dynamics, providing a framework for interpreting momentum-resolved pump-probe experiments.
We investigate the two-dimensional Hubbard model using a real-frequency implementation of the TPSC+DMFT approach. This hybrid method combines the nonlocal correlations captured by the two-particle self-consistent (TPSC) approach with the local dynamical correlations of dynamical mean field theory (DMFT). The results demonstrate that TPSC+DMFT effectively describes pseudogap physics and nonlocal fluctuations in the moderately correlated regime, while also reproducing the Mott insulating state at larger interaction strengths. For doped Mott insulators, we find that TPSC+DMFT captures the evolution of Fermi pockets into Fermi arcs, consistent with the results from cluster DMFT and photoemission studies. These findings highlight the capability of TPSC+DMFT to bridge the gap between weak and strong coupling physics in Hubbard models, providing insights into spin and charge fluctuations, as well as their role in the pseudogap formation.
The main focus of the present thesis lies on the development and application of numerical approaches based on quantum lattice models to describe collective excitations in strongly correlated materials. To this end, an efficient and balanced dual triply irreducible local expansion (D-TRILEX) approach that captures the major effects of spatial correlations in strongly correlated fermionic systems is presented, which is based on partial bosonization of the collective electronic fluctuations. Thereby, a specific choice of the bare interaction in the corresponding channel of instability allows a simultaneous description of collective excitations in different channels without facing the Fierz ambiguity problem which usually arises in methods based on partial bosonization. Compared to more elaborate methods, D-TRILEX does not require the computation of the numerically expensive four-point vertex. Furthermore, taking into account only the leading longitudinal contribution of the bosonic modes to the self-energy reduces the complexity of the diagrammatic structure drastically making the method applicable to more complex realistic multi-orbital systems. As the main result of this thesis, the D-TRILEX approach is introduced and its limit of applicability is studied by applying it to the two-dimensional single-band Hubbard model. Performing a comprehensive diagrammatic analysis by comparing D-TRILEX with methods, which consider different sets of diagrams and different approximations for the four-point vertex, we find that contributions that are not considered by the partially bosonized approximation of the four-point vertex have only a minor effect on electronic degrees of freedom and the longitudinal bosonic modes taking into account in D-TRILEX self-energy are the leading ones in a broad range of control parameters. After investigating the applicability of the D-TRILEX approach we apply the method to study collective electronic effects in the hole-doped InSe monolayers using a realistic electronic-structure model. Thereby, we find that due to the weakly screened long-range Coulomb interaction the system shows a charge density waves phase for the broad range of considered parameters, however, inside this regime we have found a coexisting ferromagnetic phase. D-TRILEX can be considered as an approximation of the more elaborate dual boson (DB) theory, which accounts for local correlations within the dynamical mean-field theory (DMFT) and nonlocal correlations perturbatively. As a further result of this thesis, we investigate the weak coupling regime of the two-dimensional single-band Hubbard model within the ladder DB approximation and compare our results to the numerically exact diagrammatic quantum Monte Carlo (DiagMC) method. We find that the DB method qualitatively captures all the different regimes. Even quantitatively DB shows a good agreement with the exact results for single- and two-particle quantities at high and moderate temperatures. However, the ladder approximation slightly underestimates the strong antiferromagnetic fluctuations that appear in the system at low temperatures. As the final result of this thesis, we perform a parametrization of the local four-point vertex in terms of the (three-point) Hedin vertex and screened interactions in the framework of single-boson exchange (SBE) decomposition. Similar to the D-TRILEX approach the computational costs are reduced by neglecting the interaction-irreducible part of the four-point vertex. The approximated vertex is used to calculate the DMFT susceptibility for different interaction strengths at different temperatures. We find that the approximation shows a good agreement in the limits of weak and strong couplings. Thereby, our study highlights the importance of the Hedin vertex for local vertex corrections.
In strongly correlated materials, quantum fluctuations taking place at the level of atomic sites play a central role. These local processes are captured nonperturbatively in dynamical mean field theory (DMFT). Nonlocal mechanisms, however, are also suspected to be essential for the proper description of some of these materials, including the well-known high-temperature cuprate superconductors and the more exotic systems of atoms adsorbed on semiconducting surfaces.In this thesis, the effect of nonlocal, potentially long-range fluctuations beyond the DMFT framework is investigated from two angles.First, long-ranged interactions in adatom systems, computed from first principles within the constrained random-phase approximation (cRPA), are treated within the extended DMFT (EDMFT) framework and its combination with the diagrammatic GW method (GW +EDMFT). They are shown to explain some material trends observed in these compounds. In addition, thanks to the spatial and temporal information delivered by GW+EDMFT, they solve existing puzzles raised by different experimental conclusions drawn from different probes.Second, a new approach to the Hubbard model including nonlocal correlation effects is proposed. Formally based on the local expansion of the triply-irreducible functional, this approach, dubbed TRILEX, approximates the renormalized three-leg electron-boson vertex by a local but frequency-dependent vertex computed from an impurity model with retarded density-density and spin-spin interactions. This local vertex is used to construct frequency- and momentum-dependent self-energies capturing nonlocal correlation effects. By construction, this method reconciles two theoretical pictures of cuprate superconductors, namely the spin fluctuation theory at weak interaction strengths and the doped-Mott-insulator picture at strong coupling.This novel approach is applied to the single-band Hubbard model. It is shown to reduce to DMFT at strong coupling and high temperatures, and to fluctuation-exchange theories at weak coupling, while substantial deviations from both limits are observed at intermediate coupling. There, upon doping, antiferromagnetic spin fluctuations are observed to lead to a momentum-differentiated spectral intensity on the Fermi surface, which is reminiscent of the formation of Fermi arcs in photoemission measurements on cuprate superconductors.These two methodological developments are intrinsically linked to algorithmic aspects which are reviewed in the last part of this work. In particular, an open-source library dedicated to massively parallel computations in the field of quantum many-body physics, TRIQS, is presented, with a focus on error estimation through the binning and jackknife methods. The hybridization-expansion continuous-time quantum Monte-Carlo algorithm in the segment picture -- central to the TRILEX method -- and its state-of-the-art implementation are explained. Finally, as a further illustration, a new lightweight method relying on this algorithm and designed for realistic systems, the combined screened-exchange and dynamical DMFT method (SEx+DDMFT), is presented together with an application to the BaCo2As2 compound.
Charge excitations in Mott insulators (MIs) are distinct from their band-insulator counterparts and they can provide a mechanism for energy harvesting in solar cells based on strongly correlated electronic materials. In this paper, we study the real-time dynamics of holon-doublon pairs in an MI connected to metallic leads using the time-dependent density-matrix renormalization-group method. The transfer of charge across the MI-metal interface is controlled by both the electron-electron interaction strength within the MI and the voltage difference between the leads. We find an overall enhancement of the charge transfer as compared to the case of a (noninteracting) band insulator-metal interface with a similar band gap. Moreover, the propagation of holon-doublon excitations within the MI dynamically changes the spin-spin correlations, introducing time-dependent phase shifts in the spin-structure factor.
Optimizing various properties of materials have been under intensive focus for many decades. In recent years, the exposure of ultrafast techniques enhance the possibility of seizing this goal. As by employing these techniques, one can expose the system under strong short fields, and investigate the transient response of various low-lying excitations. Such studies will strengthen our understanding of interplays among various degrees of freedom in correlated materials, and thereby facilitate the chance to manipulate different properties of these systems. Of our particular attention is the electron-lattice dynamics, where modulating a particular vibrational mode of the crystal lattice, will lead to the emergence of new quasi stable'' lattice structures as well as novel transient phenomena. Aside from these experimental prospects, the above mentioned goals were also the subject of an extensive theoretical studies. In this thesis, we investigated the in and out of equilibrium dynamics of electronic systems which are interacting with lattice degrees of freedom. We employed the two-time Keldysh Green's function within the framework of the dynamical mean field theory to conduct our researches on this subject. In our first attempt, we studied the dynamical response of a one-electron system under a linear electron-lattice interaction. We explored the real-time formation of a bound quasiparticle, known as a polaron, under various parameter regimes of the model. We presented that when the carrier is suddenly coupled to the lattice degrees of freedom, the long-time response of the system in the strong coupling limit exhibits both bound, polaronic, and unbound, delocalized, states. We characterized the nature of such a mixed state using the adiabatic picture of the model. We also showed that the formed polaron is dressed by excited phononic states, in contrast to its equilibrium counterpart, where the ground-state phononic states are bounding the electron. We also addressed the response of a many-body electronic system, weakly coupled to a bosonic bath, under various time dependent excitations, which brings the system to a bad-metallic phase. We employed a model which treated the Coulomb repulsive interaction locally. We presented that the relaxation of this excited state varies by the strength of the electron-electron interaction. We demonstrated that, for both insulating and metallic regimes of the system, a slow relaxation dynamics can be observed close to the phase transition from the metallic to insulating regimes of the system. Our analysis revealed that in the metallic side, this slow evolution is rooted in the spin-related physics of the problem. However, in the insulating regimes of the model, our results showed that presence of a gap in the spectral density of the insulating system prohibits the relaxation of the excited particles. We also studied the associated timescale for the adiabatic dynamics near the Mott transition. Our results exhibited that there is a nonadiabatic window'' in the metallic side of the phase transition, in which the adiabatic timescale is extremely long. We pointed out that the presence of such a window is intertwined with the slow relaxation timescale of the spin-related dynamics close to the Mott transition. We also captured an induced metallic phase near metal-to-insulator transition out of equilibrium. We initiated this study using the slave-rotor impurity solver, which is been coupled to the self-energies obtained from the first-order Luttinger-Ward functional for the Phonon degrees of freedom. We presented the emergence of the larger spectral density at small frequencies both in metallic and insulating phases close to the phase transition, which is surrounded by pseudo gaps. We explored the physics of such spectral densities as well as their dependencies on the phononic parameters. We postulated that this response of the system is due to the evolution of the system, in a non-thermal mechanism, towards a temperature which is less than the initial temperature of the system, before coupling to the bath. We also pointed out that the associated very slow relaxation dynamics is the result of the formation of psuedo-gaps around the quasiparticle peak.