Search

We present a generalization of the discrete Lehmann representation (DLR) to three-point correlation and vertex functions in imaginary time and Matsubara frequency. The representation takes the form of a linear combination of judiciously chosen exponentials in imaginary time, and products of simple poles in Matsubara frequency, which are universal for a given temperature and energy cutoff. We present a systematic algorithm to generate compact sampling grids, from which the coefficients of such an expansion can be obtained by solving a linear system. We show that the explicit form of the representation can be used to evaluate diagrammatic expressions involving infinite Matsubara sums, such as polarization functions or self-energies, with controllable, high-order accuracy. This collection of techniques establishes a framework through which methods involving three-point objects can be implemented robustly, with a substantially reduced computational cost and memory footprint.

We introduce cppdlr, a C++ library implementing the discrete Lehmann representation (DLR) of functions in imaginary time and Matsubara frequency, such as Green's functions and self-energies. The DLR is based on a low-rank approximation of the analytic continuation kernel, and yields a compact and explicit basis consisting of exponentials in imaginary time and simple poles in Matsubara frequency. cppdlr constructs the DLR basis and associated interpolation grids, and implements standard operations. It provides a flexible yet high-level interface, facilitating the incorporation of the DLR into both small-scale applications and existing large-scale software projects.

The inchworm expansion is a promising approach to solving strongly correlated quantum impurity models due to its reduction of the sign problem in real and imaginary time. However, inchworm Monte Carlo is computationally expensive, converging as $1/\sqrt{N}$ where $N$ is the number of samples. We show that the imaginary time integration is amenable to quasi Monte Carlo, with enhanced $1/N$ convergence, by mapping the Sobol low-discrepancy sequence from the hypercube to the simplex with the so-called Root transform. This extends the applicability of the inchworm method to, e.g., multi-orbital Anderson impurity models with off-diagonal hybridization, relevant for materials simulation, where continuous time hybridization expansion Monte Carlo has a severe sign problem.

We present a deterministic algorithm for the efficient evaluation of imaginary time diagrams based on the recently introduced discrete Lehmann representation (DLR) of imaginary time Green's functions. In addition to the efficient discretization of diagrammatic integrals afforded by its approximation properties, the DLR basis is separable in imaginary time, allowing us to decompose diagrams into linear combinations of nested sequences of one-dimensional products and convolutions. Focusing on the strong coupling bold-line expansion of generalized Anderson impurity models, we show that our strategy reduces the computational complexity of evaluating an $M$th-order diagram at inverse temperature $\beta$ and spectral width $\omega_{\max}$ from $\mathcal{O}((\beta \omega_{\max})^{2M-1})$ for a direct quadrature to $\mathcal{O}(M (\log (\beta \omega_{\max}))^{M+1})$, with controllable high-order accuracy. We benchmark our algorithm using third-order expansions for multi-band impurity problems with off-diagonal hybridization and spin-orbit coupling, presenting comparisons with exact diagonalization and quantum Monte Carlo approaches. In particular, we perform a self-consistent dynamical mean-field theory calculation for a three-band Hubbard model with strong spin-orbit coupling representing a minimal model of Ca$_2$RuO$_4$, demonstrating the promise of the method for modeling realistic strongly correlated multi-band materials. For both strong and weak coupling expansions of low and intermediate order, in which diagrams can be enumerated, our method provides an efficient, straightforward, and robust black-box evaluation procedure. In this sense, it fills a gap between diagrammatic approximations of the lowest order, which are simple and inexpensive but inaccurate, and those based on Monte Carlo sampling of high-order diagrams.

Dynamical mean-field theory describes the impact of strong local correlation effects in many-electron systems. While the single-particle spectral function is directly obtained within the formalism, two-particle susceptibilities can also be obtained by solving the Bethe-Salpeter equation. The solution requires handling infinite matrices in Matsubara frequency space. This is commonly treated using a finite frequency cut-off, resulting in slow linear convergence. A decomposition of the two-particle response in local and non-local contributions enables a reformulation of the Bethe-Salpeter equation inspired by the dual boson formalism. The re-formulation has a drastically improved cubic convergence with respect to the frequency cut-off, facilitating the calculation of susceptibilities in multi-orbital systems considerably. This improved convergence arises from the fact that local contributions can be measured in the impurity solver. The dual Bethe-Salpeter equation uses the fully reducible vertex which is free from vertex divergences. We benchmark the approach on several systems including the spin susceptibility of strontium ruthenate Sr$_2$RuO$_4$, a strongly correlated Hund's metal with three active orbitals.

Many-electron systems undergo a collective Larmor precession in the presence of a magnetic field. In a paramagnetic metal, the resulting spin wave provides insight into the correlation effects generated by the electron-electron interaction. Here, we use dynamical mean-field theory to investigate the collective Larmor precession in the strongly correlated regime, where dynamical correlation effects such as quasiparticle lifetimes and non-quasiparticle states are essential. We study the spin excitation spectrum, which includes a dispersive Larmor mode as well as electron-hole excitations that lead to Stoner damping. We also extract the momentum-resolved damping of slow spin waves. The accurate theoretical description of these phenomena relies on the Ward identity, which guarantees a precise cancellation of self-energy and vertex corrections at long wavelengths. Our findings pave the way towards a better understanding of spin wave damping in correlated materials.

The real-time contour formalism for Green's functions provides time-dependent information of quantum many-body systems. In practice, the long-time simulation of systems with a wide range of energy scales is challenging due to both the storage requirements of the discretized Green's function and the computational cost of solving the Dyson equation. In this manuscript, we apply a real-time discretization based on a piece-wise high-order orthogonal-polynomial expansion to address these issues. We present a superconvergent algorithm for solving the real-time equilibrium Dyson equation using the Legendre spectral method and the recursive algorithm for Legendre convolution. We show that the compact high order discretization in combination with our Dyson solver enables long-time simulations using far fewer discretization points than needed in conventional multistep methods. As a proof of concept, we compute the molecular spectral functions of H$_2$, LiH, He$_2$ and C$_6$H$_4$O$_2$ using self-consistent second-order perturbation theory and compare the results with standard quantum chemistry methods as well as the auxiliary second-order Green's function perturbation theory method.

We introduce libdlr, a library implementing the recently introduced discrete Lehmann representation (DLR) of imaginary time Green's functions. The DLR basis consists of a collection of exponentials chosen by the interpolative decomposition to ensure stable and efficient recovery of Green's functions from imaginary time or Matsbuara frequency samples. The library provides subroutines to build the DLR basis and grids, and to carry out various standard operations. The simplicity of the DLR makes it straightforward to incorporate into existing codes as a replacement for less efficient representations of imaginary time Green's functions, and libdlr is intended to facilitate this process. libdlr is written in Fortran, provides a C header interface, and contains a Python module pydlr. We also introduce a stand-alone Julia implementation, Lehmann.jl.

We consider the numerical solution of the real time equilibrium Dyson equation, which is used in calculations of the dynamical properties of quantum many-body systems. We show that this equation can be written as a system of coupled, nonlinear, convolutional Volterra integro-differential equations, for which the kernel depends self-consistently on the solution. As is typical in the numerical solution of Volterra-type equations, the computational bottleneck is the quadratic-scaling cost of history integration. However, the structure of the nonlinear Volterra integral operator precludes the use of standard fast algorithms. We propose a quasilinear-scaling FFT-based algorithm which respects the structure of the nonlinear integral operator. The resulting method can reach large propagation times, and is thus well-suited to explore quantum many-body phenomena at low energy scales. We demonstrate the solver with two standard model systems: the Bethe graph, and the Sachdev-Ye-Kitaev model.

We study the superconducting gap function of Sr$_2$RuO$_4$. By solving the linearized Eliashberg equation with a correlated pairing vertex extracted from a dynamical mean-field calculation we identify the dominant pairing channels. An analysis of the candidate gap functions in orbital and quasiparticle band basis reveals that an inter-orbital singlet pairing of even parity is in agreement with experimental observations. It reconciles in particular the occurrence of a two-component order parameter with the presence of line-nodes of quasiparticles along the c-axis in the superconducting phase. The strong angular dependence of the gap along the Fermi surface is in stark contrast to its quasi-locality when expressed in the orbital basis. We identify local inter-orbital spin correlations as the driving force for the pairing and thus reveal the continuation of Hund's physics into the superconducting phase., Comment: 20 pages, 12 figures

Quantum many-body systems in thermal equilibrium can be described by the imaginary time Green's function formalism. However, the treatment of large molecular or solid ab inito problems with a fully realistic Hamiltonian in large basis sets is hampered by the storage of the Green's function and the precision of the solution of the Dyson equation. We present a Legendre-spectral algorithm for solving the Dyson equation that addresses both of these issues. By formulating the algorithm in Legendre coefficient space, our method inherits the known faster-than-exponential convergence of the Green's function's Legendre series expansion. In this basis, the fast recursive method for Legendre polynomial convolution, enables us to develop a Dyson equation solver with quadratic scaling. We present benchmarks of the algorithm by computing the dissociation energy of the helium dimer He$_2$ within dressed second-order perturbation theory. For this system, the application of the Legendre spectral algorithm allows us to achieve an energy accuracy of $10^{-9} E_h$ with only a few hundred expansion coefficients.

The nonequilibrium dynamics of correlated many-particle systems is of interest in connection with pump-probe experiments on molecular systems and solids, as well as theoretical investigations of transport properties and relaxation processes. Nonequilibrium Green's functions are a powerful tool to study interaction effects in quantum many-particle systems out of equilibrium, and to extract physically relevant information for the interpretation of experiments. We present the open-source software package NESSi (The Non-Equilibrium Systems Simulation package) which allows to perform many-body dynamics simulations based on Green's functions on the L-shaped Kadanoff-Baym contour. NESSi contains the library libcntr which implements tools for basic operations on these nonequilibrium Green's functions, for constructing Feynman diagrams, and for the solution of integral and integro-differential equations involving contour Green's functions. The library employs a discretization of the Kadanoff-Baym contour into time $N$ points and a high-order implementation of integration routines. The total integrated error scales up to $\mathcal{O}(N^{-7})$, which is important since the numerical effort increases at least cubically with the simulation time. A distributed-memory parallelization over reciprocal space allows large-scale simulations of lattice systems. We provide a collection of example programs ranging from dynamics in simple two-level systems to problems relevant in contemporary condensed matter physics, including Hubbard clusters and Hubbard or Holstein lattice models. The libcntr library is the basis of a follow-up software package for nonequilibrium dynamical mean-field theory calculations based on strong-coupling perturbative impurity solvers., Comment: 86 pages, 15 figures

The crossover from fluctuating atomic constituents to a collective state as one lowers temperature or energy is at the heart of the dynamical mean-field theory description of the solid state. We demonstrate that the numerical renormalization group is a viable tool to monitor this crossover in a real-materials setting. The renormalization group flow from high to arbitrarily small energy scales clearly reveals the emergence of the Fermi-liquid state of Sr$_2$RuO$_4$. We find a two-stage screening process, where orbital fluctuations are screened at much higher energies than spin fluctuations, and Fermi-liquid behavior, concomitant with spin coherence, below a temperature of 25 K. By computing real-frequency correlation functions, we directly observe this spin--orbital scale separation and show that the van Hove singularity drives strong orbital differentiation. We extract quasiparticle interaction parameters from the low-energy spectrum and find an effective attraction in the spin-triplet sector.

We study the magnetic susceptibility in the normal state of Sr$_2$RuO$_4$ using dynamical mean-field theory including dynamical vertex corrections. Besides the well known incommensurate response, our calculations yield quasi-local spin fluctuations which are broad in momentum and centered around the $\Gamma$ point, in agreement with recent inelastic neutron scattering experiments [P. Steffens, et al., Phys. Rev. Lett. 122, 047004 (2019)]. We show that these quasi-local fluctuations are controlled by the Hund's coupling and account for the dominant contribution to the momentum-integrated response. While all orbitals contribute equally to the incommensurate response, the enhanced $\Gamma$ point response originates from the planar xy orbital., Comment: 6 pages, 5 figures

We study equilibrium and nonequilibrium properties of electron-phonon systems described by the Hubbard-Holstein model using the dynamical mean-field theory. In equilibrium, we benchmark the results for impurity solvers based on the one-crossing approximation and slave-rotor approximation against non-perturbative numerical renormalization group reference data. We also examine how well the low energy properties of the electron-boson coupled systems can be reproduced by an effective static electron-electron interaction. The one-crossing and slave-rotor approximations are then used to simulate insulator-to-metal transitions induced by a sudden switch-on of the electron-phonon interaction. The slave-rotor results suggest the existence of a critical electron-phonon coupling above which the system is transiently trapped in a non-thermal metallic state with coherent quasiparticles. The same quench protocol in the one-crossing approximation results in a bad metallic state., Comment: 16 pages, 19 figures

Photoinduced hidden phases are often observed in materials with intertwined orders. Understanding the formation of these non-thermal phases is challenging and requires a resolution of the cooperative interplay between different orders on the ultrashort timescale. In this work, we demonstrate that non-equilibrium photo-excitations can induce a state with spin-orbital orders entirely different from the equilibrium state in the three-quarter-filled two-band Hubbard model. We identify a general mechanism governing the transition to the hidden state, which relies on a non-thermal partial melting of the intertwined orders mediated by photoinduced charge excitations in the presence of strong spin-orbital exchange interactions. Our study theoretically confirms the crucial role played by orbital degrees of freedom in the light-induced dynamics of strongly correlated materials and it shows that the switching to hidden states can be controlled already on the fs timescale of the electron dynamics.

Using a nonequilibrium implementation of the extended dynamical mean field theory (EDMFT) we simulate the relaxation after photo excitation in a strongly correlated electron system with antiferromagnetic spin interactions. We consider the $t$-$J$ model and focus on the interplay between the charge- and spin-dynamics in different excitation and doping regimes. The appearance of string states after a weak photo excitation manifests itself in a nontrivial scaling of the relaxation time with the exchange coupling and leads to a correlated oscillatory evolution of the kinetic energy and spin-spin correlation function. A strong excitation of the system, on the other hand, suppresses the spin correlations and results in a relaxation that is controlled by hole scattering. We discuss the possibility of detecting string states in optical and cold atom experiments.

We extend the self-energy functional theory (SFT) to the case of interacting lattice bosons in the presence of symmetry breaking and quenched disorder. The self-energy functional we derive depends only on the self-energies of the disorder-averaged propagators, allowing for the construction of general non-perturbative approximations. Using a simple single-site reference system with only three variational parameters, we are able to reproduce numerically exact quantum Monte Carlo (QMC) results on local observables of the Bose-Hubbard model with box disorder with high accuracy. At strong interactions, the phase boundaries are reproduced qualitatively but shifted with respect to the ones observed with QMC due to the extremely low condensate fraction in the superfluid phase. Deep in the strongly-disordered weakly-interacting regime, the simple reference system employed is insufficient and no stationary solutions can be found within its restricted variational subspace. By systematically analyzing thermodynamical observables and the spectral function, we find that the strongly-interacting Bose glass is characterized by different regimes, depending on which local occupations are activated as a function of the disorder strength. We find that the particles delocalize into isolated superfluid lakes over a strongly localized background around maximally-occupied sites whenever these sites are particularly rare. Our results indicate that the transition from the Bose glass to the superfluid phase around unit filling at strong interactions is driven by the percolation of superfluid lakes which form around doubly occupied sites., Comment: 34 pages, 5 figures

We study the relaxation dynamics of photo-carriers in the paramagnetic Mott insulating phase of the half-filled two-band Hubbard model. Using nonequilibrium dynamical mean field theory, we excite charge carriers across the Mott gap by a short hopping modulation, and simulate the evolution of the photo-doped population within the Hubbard bands. We observe an ultrafast charge-carrier relaxation driven by emission of local spin excitations with an inverse relaxation time proportional to the Hund's coupling. The photo-doping generates additional side-bands in the spectral function, and for strong Hund's coupling, the photo-doped population also splits into several resonances. The dynamics of the local many-body states reveals two effects, thermal blocking and kinetic freezing, which manifest themselves when the Hund's coupling becomes of the order of the temperature or the bandwidth, respectively. These effects, which are absent in the single-band Hubbard model, should be relevant for the interpretation of experiments on correlated materials with multiple active orbitals. In particular, the features revealed in the non-equilibrium energy distribution of the photo-carriers are experimentally accessible, and provide information on the role of the Hund's coupling in these materials., Comment: 15 pages, 11 figures

We derive the reciprocal cluster mean-field method to study the strongly-interacting bosonic Harper-Hofstadter-Mott model. The system exhibits a rich phase diagram featuring band insulating, striped superfluid, and supersolid phases. Furthermore, for finite hopping anisotropy we observe gapless uncondensed liquid phases at integer fillings, which are analyzed by exact diagonalization. The liquid phases at fillings 1 and 3 exhibit the same band fillings as the fermionic integer quantum Hall effect, while the phase at filling 2 is CT-symmetric with zero charge response. We discuss how these phases become gapped on a quasi-one-dimensional cylinder, leading to a quantized Hall response, which we characterize by introducing a suitable measure for non-trivial many-body topological properties. Incompressible metastable states at fractional filling are also observed, indicating competing fractional quantum Hall phases. The combination of reciprocal cluster mean-field and exact diagonalization yields a promising method to analyze the properties of bosonic lattice systems with non-trivial unit cells in the thermodynamic limit., Comment: 16 pages, 11 figures

We derive the self-energy functional theory for bosonic lattice systems with broken $U(1)$ symmetry by parametrizing the bosonic Baym-Kadanoff effective action in terms of one- and two-point self-energies. The formalism goes beyond other approximate methods such as the pseudoparticle variational cluster approximation, the cluster composite boson mapping, and the Bogoliubov+U theory. It simplifies to bosonic dynamical-mean field theory when constraining to local fields, whereas when neglecting kinetic contributions of non-condensed bosons it reduces to the static mean-field approximation. To benchmark the theory we study the Bose-Hubbard model on the two- and three-dimensional cubic lattice, comparing with exact results from path integral quantum Monte Carlo. We also study the frustrated square lattice with next-nearest neighbor hopping, which is beyond the reach of Monte Carlo simulations. A reference system comprising a single bosonic state, corresponding to three variational parameters, is sufficient to quantitatively describe phase-boundaries, and thermodynamical observables, while qualitatively capturing the spectral functions, as well as the enhancement of kinetic fluctuations in the frustrated case. On the basis of these findings we propose self-energy functional theory as the omnibus framework for treating bosonic lattice models, in particular, in cases where path integral quantum Monte Carlo methods suffer from severe sign problems (e.g. in the presence of non-trivial gauge fields or frustration). Self-energy functional theory enables the construction of diagrammatically sound approximations that are quantitatively precise and controlled in the number of optimization parameters, but nevertheless remain computable by modest means., Comment: 23 pages, 13 figures

We investigate features in the single-particle spectral function beyond the Hubbard bands in the strongly correlated normal phase of the Bose-Hubbard model. There are two distinct classes of additional peaks generated by the bosonic statistics. The first type is thermally activated Hubbard "sidebands", with the same physical origin as the zero-temperature Hubbard bands, but generated by excitations from thermally activated local occupation number states. The second class are two-particle fluctuation resonances driven by the lattice dynamics. In the unity filling Mott insulator, this takes the form of a localized triplon combined with a dispersing holon. Both types of resonances also manifest themselves in the structure factor and the interaction modulation spectra obtained from nonequilibrium bosonic dynamical mean-field theory calculations. Our findings explain experimental lattice modulation and Bragg spectroscopy results, and they predict a strong temperature dependence of the first sideband, thereby opening the door to precise thermometry of strongly correlated lattice bosons., Comment: 9 pages, 3 figures, including appendices

We develop the nonequilibrium extension of bosonic dynamical mean field theory (BDMFT) and a Nambu real-time strong-coupling perturbative impurity solver. In contrast to Gutzwiller mean-field theory and strong coupling perturbative approaches, nonequilibrium BDMFT captures not only dynamical transitions, but also damping and thermalization effects at finite temperature. We apply the formalism to quenches in the Bose-Hubbard model, starting both from the normal and Bose-condensed phases. Depending on the parameter regime, one observes qualitatively different dynamical properties, such as rapid thermalization, trapping in metastable superfluid or normal states, as well as long-lived or strongly damped amplitude oscillations. We summarize our results in non-equilibrium "phase diagrams" which map out the different dynamical regimes., Comment: 18 pages, 8 figures

We argue that, because of the quantum-entanglement, the local physics of the strongly-correlated materials at zero temperature is described in very good approximation by a simple generalized Gibbs distribution, which depends on a relatively small number local quantum thermodynamical potentials. We demonstrate that our statement is exact in certain limits, and we perform numerical calculations of the iron compounds FeSe and FeTe and of the elemental cerium by employing the Gutzwiller Approximation (GA) that strongly support our theory in general., Comment: Manuscript: 5 pages, 3 figures. Supplemental material: 1 page, 1 figure

We show that repulsive local Coulomb interaction alone can drive valence-skipping charge disproportionation in the degenerate d-band, resulting in effective negative-U. This effect is shown to originate from anisotropic orbital-multipole scattering, and occurs only for $d^1$, $d^4$, $d^6$ and $d^9$ fillings (and their immediate surroundings). Explicit boundaries for valence-skipping are derived and the paramagnetic phase diagram for $d^4$ and $d^6$ is calculated. We also establish that the valence-skipping metal is very different, in terms of its local valence distribution, compared to the atomic-like Hund's metal. These findings explains why transition metal compounds with the aforementioned d-band fillings are more prone to valence-skipping charge order and anomalous superconductivity., Comment: 5 pages, 5 figures

We develop a method for calculating the self-energy of a quantum impurity coupled to a continuous bath by stochastically generating a distribution of finite Anderson models that are solved by exact diagonalization, using the noninteracting local spectral function as a probability distribution for the sampling. The method enables calculation of the full analytic self-energy and single-particle Green's function in the complex frequency plane, without analytic continuation, and can be used for both finite and zero temperature at arbitrary fillings. Results are in good agreement with imaginary frequency data from continuous-time quantum Monte Carlo calculations for the single impurity Anderson model and the two-orbital Hubbard model within dynamical mean field theory (DMFT) as well as real frequency data for self energy of the single band Hubbard model within DMFT using numerical renormalization group. The method should be applicable to a wide range of quantum impurity models and particularly useful when high-precision real frequency results are sought., Comment: Two band calculation updated, 4 pages, 4 figures

We present a self-consistent numerical approach to solve the Gutzwiller variational problem for general multi-band models with arbitrary on-site interaction. The proposed method generalizes and improves the procedure derived by Deng et al., Phys. Rev. B. 79 075114 (2009), overcoming the restriction to density-density interaction without increasing the complexity of the computational algorithm. Our approach drastically reduces the problem of the high-dimensional Gutzwiller minimization by mapping it to a minimization only in the variational density matrix, in the spirit of the Levy and Lieb formulation of DFT. For fixed density the Gutzwiller renormalization matrix is determined as a fixpoint of a proper functional, whose evaluation only requires ground-state calculations of matrices defined in the Gutzwiller variational space. Furthermore, the proposed method is able to account for the symmetries of the variational function in a controlled way, reducing the number of variational parameters. After a detailed description of the method we present calculations for multi-band Hubbard models with full (rotationally invariant) Hund's rule on-site interaction. Our analysis shows that the numerical algorithm is very efficient, stable and easy to implement. For these reasons this method is particularly suitable for first principle studies -- e.g., in combination with DFT -- of many complex real materials, where the full intra-atomic interaction is important to obtain correct results., Comment: 19 pages, 7 figures

We present a spectral weight conserving formalism for Fermionic thermal Green's functions that are discretized in imaginary time and thus periodic in imaginary ("Matsubara") frequency. The formalism requires a generalization of the Dyson equation and the Baym-Kadanoff-Luttinger-Ward functional for the free energy. A conformal transformation is used to analytically continue the periodized Matsubara Green's function to the continuous real axis in a way that conserves the discontinuity at t=0 of the corresponding real-time Green's function. For given discretization the method allows numerical Green's function calculations of very high precision and it appears to give a well controlled convergent approximation as we decrease the discretization interval. The ideas are tested on Dynamical Mean Field Theory calculations of the paramagnetic Hubbard model.

We extend the time-dependent Gutzwiller variational approach, recently introduced by Schir\`o and Fabrizio, Phys. Rev. Lett. 105 076401 (2010), to impurity problems. Furthermore, we derive a consistent theory for the steady state, and show its equivalence with the previously introduced nonequilibrium steady-state extension of the Gutzwiller approach. The method is shown to be able to capture dissipation in the leads, so that a steady state is reached after a sufficiently long relaxation time. The time-dependent method is applied to the single orbital Anderson impurity model at half-filling, modeling a quantum dot coupled to two leads. In these first exploratory calculations the Gutzwiller projector is limited to act only on the impurity. The strengths and the limitations of this approximation are assessed via comparison with state of the art continuous time quantum Monte Carlo results. Finally, we discuss how the method can be systematically improved by extending the region of action of the Gutzwiller projector., Comment: 13 pages, 6 figures

We consider the finite temperature metal-insulator transition in the half filled paramagnetic Hubbard model on the infinite dimensional Bethe lattice. A new method for calculating the Dynamical Mean Field Theory fixpoint surface in the phase diagram is presented and shown to be free from the convergence problems of standard forward recursion. The fixpoint equation is then analyzed using dynamical systems methods. On the fixpoint surface the eigenspectra of its Jacobian is used to characterize the hysteresis boundaries of the first order transition line and its second order critical end point. The critical point is shown to be a cusp catastrophe in the parameter space, opening a pitchfork bifurcation along the first order transition line, while the hysteresis boundaries are shown to be saddle-node bifurcations of two merging fixpoints. Using Landau theory the properties of the critical end point is determined and related to the critical eigenmode of the Jacobian. Our findings provide new insights into basic properties of this intensively studied transition., Comment: 11 pages, 12 figures, 1 table