We discuss a nonequilibrium dynamical mean-field framework for simulating inhomogeneous Hubbard models with local disorders. Our approach treats electron interactions and disorders on equal footing, by considering only local dynamical fluctuations. The theory reduces to nonequilibrium dynamical mean-field theory in the presence of only electron-electron interactions and to the coherent potential approximation in noninteracting systems with disorders. Both time-dependent and steady-state problems are treated by implementing the theory on the three-branch Kadanoff-Baym contour and two-branch Schwinger-Keldysh contour, respectively. Benchmarks on a eight-site cube show that the method yields rather accurate spectral functions in both the weakly and strongly interacting regimes. In a cubic lattice, we demonstrate energy conservation after an interaction quench and thermalization after just a few hopping times in both clean and disordered systems. As an application, we study transport through a serial double quantum-dot sandwiched between two leads, focusing on the current and dot occupations after a voltage quench.
We discuss a nonequilibrium dynamical mean-field framework for simulating inhomogeneous Hubbard models with local disorders. Our approach treats electron interactions and disorders on equal footing, by considering only local dynamical fluctuations. The theory reduces to nonequilibrium dynamical mean-field theory in the presence of only electron-electron interactions and to the coherent potential approximation in noninteracting systems with disorders. Both time-dependent and steady-state problems are treated by implementing the theory on the three-branch Kadanoff-Baym contour and two-branch Schwinger-Keldysh contour, respectively. Benchmarks on a eight-site cube show that the method yields rather accurate spectral functions in both the weakly and strongly interacting regimes. In a cubic lattice, we demonstrate energy conservation after an interaction quench and thermalization after just a few hopping times in both clean and disordered systems. As an application, we study transport through a serial double quantum-dot sandwiched between two leads, focusing on the current and dot occupations after a voltage quench.
We present a solution for the nonequilibrium dynamics of an interacting disordered system. The approach adapts the combination of the equilibrium dynamical mean-field theory and the equilibrium coherent potential approximation methods to the nonequilibrium many-body formalism, using the Kadanoff-Baym-Keldysh complex time contour, for the dynamics of interacting disordered systems away from equilibrium. We use our time domain solution to obtain the equilibrium density of states of the disordered interacting system described by the Anderson-Hubbard model, bypassing the necessity for the cumbersome analytical continuation process. We further apply the nonequilibrium solution to the interaction quench problem for an isolated disordered system. Here, the interaction is abruptly changed from zero (noninteracting system) to another constant (finite) value at which it is subsequently kept. We observe via the time dependence of the potential, kinetic, and total energies the effect of disorder on the relaxation of the system as a function of final interaction strength. The real-time approach has the potential to shed light on the fundamental role of disorder in the nonequilibrium dynamics of interacting quantum systems.
Dynamical mean-field theory (DMFT) is a widely used, nonperturbative computational framework for strongly correlated electron systems. This thesis extends DMFT in two directions: treating disorder and higher-order correlators away from equilibrium, and using near-term quantum computers as impurity solvers. In the first part, we develop two methodological advances for the nonequilibrium DMFT treatment of disordered, interacting systems. The first overcomes the central computational bottleneck of the recently developed nonequilibrium DMFT plus coherent potential approximation (DMFT+CPA) by replacing its explicit numerical disorder averaging with an exact analytical integration, exploiting a normalization property of the Schwinger-Keldysh contour; applied to a disorder-and-interaction quench of the Anderson-Hubbard model, this reveals that weak random disorder promotes thermalization. The second demonstrates that the out-of-time-order correlator (OTOC) can be successfully computed within the existing DMFT+CPA framework, by generalizing the Schwinger-Keldysh contour to a double-folded structure; as a proof-of-principle demonstration, weak disorder is found to accelerate the OTOC's decay in the half-filled Anderson-Hubbard model in the intermediate-coupling regime. In the second part, we explore the use of near-term quantum computers as DMFT impurity solvers. A stable real-time DMFT iteration scheme, designed for compatibility with quantum hardware, is first developed and validated classically against the metal-to-insulator transition. A symmetry-adapted variational quantum eigensolver (VQE) is then benchmarked as the Anderson impurity solver for the converged DMFT Hamiltonian, and is found to reproduce ground-state energies accurately and to extract the single-particle Green's function reliably at intermediate to strong coupling, while accurate ground-state energies are shown not to guarantee accurate dynamics at weak coupling, a key challenge for extending such approaches to correlated metals. Together, we demonstrate that DMFT can be methodologically extended to study disordered and interacting systems away from equilibrium, and to incorporate near-term quantum hardware as a genuinely new computational resource to enhance its capabilities.
Correlated, disordered lattice fermion systems out of equilibrium have been understudied, at least in part for lack of tools to study them. The development of a new method, the NE-DMFT+CPA, is explicated, for investigating the properties of these systems. This method combines DMFT and the CPA using the Kadanoff-Baym-Keldysh contour-time formalism to study the dynamics of interacting disordered systems out of equilibrium. The method is detailed, then applied to the case of the Anderson-Hubbard model: first to replicate well-known equilibrium results, then to study the case of the model under the effect of an interaction quench. Here the effect of disorder on the relaxation of the system is studied by the time-dependence of the energies and momentum distribution function, which display the screening of the Coulomb interaction by disorder. It is then used to examine the thermalization of the same system by using the fluctuation-dissipation theorem to extract the distribution function and thereby the effective temperature after relaxation. This allows the characterization of the effect of the interplay of the disorder strength and the Coulomb interaction on the thermalization behavior. It is found that the heating of the system due to the quench is mitigated by stronger disorder. Finally, current progress of the work to apply the method to binary alloy disorder and an interaction ramp is summarized. Binary disorder appears to have the effect of damping the initial transient of the potential energy more strongly than the comparable uniformly-disordered system. The relaxation of the energies under the effect of an interaction ramp suggests that the crossing of the kinetic energies may be affected by the form of the time-dependent interaction.
This work is motivated by the fact that the investigation of non-equilibrium phenomena in strongly correlated electron systems has developed into one of the most active and exciting branches of condensed matter physics as it provides rich new insights that could not be obtained from the study of equilibrium situations. However, a theoretical description of those phenomena is missing. Therefore, in this thesis, we develop a numerical method that can be used to study two minimal models -- the Hubbard model and the Anderson impurity model with general parameter range and time dependence. We begin by introducing the theoretical framework and the general features of the Hubbard model. We then describe the dynamical mean field theory (DMFT), which was first invented by Georges in 1992. It provides a feasible way to approach strongly correlated electron systems and reduces the complexity of the calculations via a mapping of lattice models onto quantum impurity models subject to a self-consistency condition. We employ the non-equilibrium extension of DMFT and map the Hubbard model to the single impurity Anderson model (SIAM). Since the fundamental component of the DMFT method is a solver of the single impurity Anderson model, we continue with a description of the formalism to study the real-time dynamics of the impurity model staring at its thermal equilibrium state. We utilize the non-equilibrium strong-coupling perturbation theory and derive semi-analytical approximation methods such as the non-crossing approximation (NCA) and the one-crossing approximation (OCA). We then use the Quantum Monte-Carlo method (QMC) as a numerically exact method and present proper measurements of local observables, current and Green's functions. We perform simulations of the current after a quantum quench from equilibrium by rapidly applying a bias voltage in a wide range of initial temperatures. The current exhibits short equilibrium times and saturates upon the decrease of temperature at all times, indicating Kondo behavior both in the transient regime and in the steady state. However, this bare QMC solver suffers from a dynamical sign problem for long time propagations. To overcome the limitations of this bare treatment, we introduce the ``Inchworm algorithm'', based on iteratively reusing the information obtained in previous steps to extend the propagation to longer times and stabilize the calculations. We show that this algorithm greatly reduces the required order for each simulation and re-scales the exponential challenge to quadratic in time. We introduce a method to compute Green's functions, spectral functions, and currents for inchworm Monte Carlo and show how systematic error assessments in real time can be obtained. We illustrate the capabilities of the algorithm with a study of the behavior of quantum impurities after an instantaneous voltage quench from a thermal equilibrium state. We conclude with the applications of the unbiased inchworm impurity solver to DMFT calculations. We employ the methods for a study of the one-band paramagnetic Hubbard model on the Bethe lattice in equilibrium, where the DMFT approximation becomes exact. We begin with a brief introduction of the Mott metal insulator phase diagram. We present the results of both real time Green's functions and spectral functions from our nonequilibrium calculations. We observe the metal-insulator crossover as the on-site interaction is increased and the formation of a quasi-particle peak as the temperature is lowered. We also illustrate the convergence of our algorithms in different aspects.
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.
We study the paramagnetic Anderson-Hubbard model using an extension of dynamical mean field theory (DMFT), known as statistical DMFT, that allows us to treat disorder and strong electronic correlations on equal footing. An approximate nonlocal Green's function is found for individual disorder realizations and then configuration averaged. We apply this method to two-dimensional lattices with up to 1000 sites in the strong disorder limit, where an atomic-limit approximation is made for the self-energy. We investigate the scaling of the inverse participation ratio at quarter- and half-filling, and find a nonmonotonic dependence of the localization length on the interaction strength. For strong disorder, we do not find evidence for an insulator-metal transition, and the disorder potential becomes unscreened near the Mott transition. Furthermore, strong correlations suppress the Altshuler-Aronov density of states anomaly near half-filling.