Search

Self-consistent dynamical approximations for strongly correlated fermion systems are particularly successful in capturing the dynamical competition of local correlations. In these, the effect of spatially extended degrees of freedom is usually only taken into account in a mean field fashion or as a secondary effect. As a result, critical exponents associated with phase transitions have a mean field character. Here we demonstrate that diagrammatic multiscale methods anchored around local approximations are indeed capable of capturing the non-mean-field nature of the critical point of the lattice model encoded in a nonvanishing anomalous dimension and of correctly describing the transition to mean-field-like behavior as the number of spatial dimensions increases. [ABSTRACT FROM AUTHOR]

The nonequilibrium spectral properties of the Anderson impurity model with a chemical potential bias are investigated within a numerically exact real-time quantum Monte Carlo formalism. The two-time correlation function is computed in a form suitable for nonequilibrium dynamical mean field calculations. Additionally, the evolution of the model's spectral properties are simulated in an alternative representation, defined by a hypothetical but experimentally realizable weakly coupled auxiliary lead. The voltage splitting of the Kondo peak is confirmed and the dynamics of its formation after a coupling or gate quench are studied. This representation is shown to contain additional information about the dot's population dynamics. Further, we show that the voltage-dependent differential conductance gives a reasonable qualitative estimate of the equilibrium spectral function, but significant qualitative differences are found including incorrect trends and spurious temperature dependent effects. [ABSTRACT FROM AUTHOR]

We present results for the equation of state of the two-dimensional Hubbard model on an isotropic square lattice as obtained from a controlled and numerically exact large-cluster dynamical mean field simulation. Our results are obtained for large but finite systems and are extrapolated to infinite system size using a known finite-size scaling relation, and are supplemented by reliable error bars accounting for all sources of errors. We establish the importance of examining the decay of spatial spin correlations to determine whether a sufficiently large cluster has been used and with this in mind we present the energy, entropy, double occupancy, and nearest-neighbor spin correlations extrapolated to the thermodynamic limit. We discuss the implications of these calculations on pseudogap physics of the 2D Hubbard model away from half filling, where we find a strong behavioral shift in energy below a temperature T* which becomes more pronounced for larger clusters. Finally, we provide reference calculations and tables for the equation of state for values of doping away from half filling which are of interest to cold-atom experiments. [ABSTRACT FROM AUTHOR]

Recently developed numerical methods have enabled the explicit construction of the superconducting state of the Hubbard model of strongly correlated electrons in parameter regimes where die model also exhibits a pseudogap and a Mott insulating phase. dx²-y² symmetry superconductivity is found to occur in proximity to the Mott insulator, but separated from it by a pseudogapped nonsuperconducting phase. The superconducting transition temperature and order parameter amplitude are found to be maximal at the onset of the normal-state pseudogap. The emergence of superconductivity from the normal state pseudogap leads to a decrease in the excitation gap. All of these features are consistent witii die observed behavior of the copper-oxide superconductors. [ABSTRACT FROM AUTHOR]

We investigate the dynamical and steady-state spin response of the nonequilibrium Anderson model to magnetic fields, bias voltage, and temperature using a numerically exact method combining a bold-line quantum Monte Carlo technique with a memory function formalism. We obtain converged results in a range of previously inaccessible regimes, in particular calculating the spin dynamics for a range of temperatures down to the crossover to the Kondo domain. We provide predictions for nonequilibrium phenomena, including nonmonotonic temperature dependence of observables at high bias voltage and oscillatory quench dynamics at high magnetic fields. [ABSTRACT FROM AUTHOR]

We show that for multiorbital quantum impurity models the noncrossing approximation and one-crossing approximation versions of the self-consistent hybridization expansions violate the sum rules relating the coefficients of the high-frequency expansion of the self-energy and the product of the self-energy and Green's function to thermodynamic expectation values. Comparison of noncrossing/one-crossing results to numerically exact quantum Monte Carlo calculations shows that the consistency with sum rules provides a useful estimate of the reliability of the approximations. The sum rule violations are more pronounced, and therefore the quality of the noncrossing/one-crossing approximation is poorer, in situations with multiple orbitals and away from particle-hole symmetry but becomes less severe as the correlation strength increases. The one-crossing approximation is markedly superior to the noncrossing approximation. [ABSTRACT FROM AUTHOR]

The development of polynomial cost solvers for correlated quantum impurity models, with controllable errors, is a central challenge in quantum many-body physics, where these models find applications ranging from nanoscience to the dynamical mean-field theory (DMFT). Here, we describe how configuration interaction (CI) approximations to exact diagonalization (ED) may be used as solvers in DMFT. CI approximations retain the main advantages of ED, such as the ability to treat general interactions and off-diagonal hybridizations and to obtain real spectral information, but are of polynomial cost. Furthermore, their errors can be controlled by monitoring the convergence of physical quantities as a function of the CI hierarchy. Using benchmark DMFT applications, such as single-site DMFT of the one-dimensional Hubbard model and 2x2 cluster DMFT of the two-dimensional Hubbard model, we show that CI approximations allow us to obtain near-exact ED results for a tiny fraction of the cost. This is true over the entire range of interaction strengths including "difficult" regimes, such as in the pseudogap phase of the two-dimensional Hubbard model. We use the ability of CI approximations to treat large numbers of orbitals to demonstrate convergence of the bath representation in the 2 x 2 cluster DMFT using a 24-bath orbital representation. CI approximations thus form a promising route to extend ED to problems that are currently difficult to study using other solvers such as continuous-time quantum Monte Carlo, including impurity models with large numbers of orbitals and general interactions. [ABSTRACT FROM AUTHOR]

A method is presented for the unbiased numerical computation of two-particle response functions of correlated electron materials via a solution of the dynamical mean-field equations in the presence of a perturbing field. The power of the method is demonstrated via a computation of the Raman B1g and B2g scattering intensities of the two-dimensional Hubbard model in parameter regimes believed to be relevant to high-temperature superconductivity. The theory reproduces the "two-magnon" peak characteristic of the Raman intensity of insulating parent compounds of high-Tc, copper oxide superconductors, and shows how it evolves to a quasiparticle response, as carriers are added. The method can be applied in any situation where a solution of equilibrium dynamical mean-field equations is feasible. [ABSTRACT FROM AUTHOR]

A Monte Carlo sampling of diagrammatic corrections to the noncrossing approximation is shown to provide numerically exact estimates of the long-time dynamics and steady-state properties of nonequilibrium quantum impurity models. This "bold" expansion converges uniformly in time and significantly ameliorates the sign problem that has heretofore limited the power of real-time Monte Carlo approaches to strongly interacting real-time quantum problems. The approach enables the study of previously intractable problems ranging from generic long-time nonequilibrium transport characteristics in systems with large on-site repulsion to the direct description of spectral functions on the real frequency axis in dynamical mean field theory. [ABSTRACT FROM AUTHOR]

We consider the effects of local impurities in highly frustrated diamond-lattice antiferromagnets, which exhibit large but nonextensive ground-state degeneracies. Such models are appropriate to many A-site magnetic spinels. We argue very generally that sufficiently dilute impurities induce an ordered magnetic ground state and provide a mechanism of degeneracy breaking. The states that are selected can be determined by a "swiss cheese model" analysis, which we demonstrate numerically for a particular impurity model in this case. Moreover, we present criteria for estimating the stability of the resulting ordered phase to a competing frozen (spin glass) one. The results may explain the contrasting finding of frozen and ordered ground states in CoAl2O4 and MnSc2S4, respectively. [ABSTRACT FROM AUTHOR]

We present momentum-resolved single-particle spectra for the three-dimensional Hubbard model for the paramagnetic and antiferromagnetically ordered phase obtained within the dynamical cluster approximation. The effective cluster problem is solved by continuous-time quantum Monte Carlo simulations. The absence of a time discretization error and the ability to perform Monte Carlo measurements directly in Matsubara frequencies enable us to analytically continue the self-energies by maximum entropy, which is essential to obtaining momentum-resolved spectral functions for the Néel state. We investigate the dependence on temperature and interaction strength and the effect of magnetic frustration introduced by a next-nearest-neighbor hopping. One particular question we address here is the influence of the frustrating interaction on the metal-insulator transition of the three-dimensional Hubbard model. [ABSTRACT FROM AUTHOR]

We present a submatrix update algorithm for the continuous-time auxiliary field method that allows the simulation of large lattice and impurity problems. The algorithm takes optimal advantage of modern CPU architectures by consistently using matrix instead of vector operations, resulting in a speedup of a factor of ≈8 and thereby allowing access to larger systems and lower temperature. We illustrate the power of our algorithm at the example of a cluster dynamical mean field simulation of the Néel transition in the three-dimensional Hubbard model, where we show momentum dependent self-energies for clusters with up to 100 sites. [ABSTRACT FROM AUTHOR]

We present a numerically exact inchworm Monte Carlo method for equilibrium multiorbital quantum impurity problems with general interactions and hybridizations. We show that the method, originally developed to overcome the dynamical sign problem in certain real-time propagation problems, can also overcome the sign problem as a function of temperature for equilibrium quantum impurity models. This is shown in several cases where the current method of choice, the continuous-time hybridization expansion, fails due to the sign problem. Our method therefore enables simulations of impurity problems as they appear in embedding theories without further approximations, such as the truncation of the hybridization or interaction structure or a discretization of the impurity bath with a set of discrete energy levels, and eliminates a crucial bottleneck in the simulation of ab initio embedding problems. [ABSTRACT FROM AUTHOR]

The steadily increasing size of scientific Monte Carlo simulations and the desire for robust, correct, and reproducible results necessitates rigorous testing procedures for scientific simulations in order to detect numerical problems and programming bugs. However, the testing paradigms developed for deterministic algorithms have proven to be ill suited for stochastic algorithms. In this paper we demonstrate explicitly how the technique of statistical hypothesis testing, which is in wide use in other fields of science, can be used to devise automatic and reliable tests for Monte Carlo methods, and we show that these tests are able to detect some of the common problems encountered in stochastic scientific simulations. We argue that hypothesis testing should become part of the standard testing toolkit for scientific simulations. [ABSTRACT FROM AUTHOR]

The development of numerical methods capable of simulating realistic materials with strongly correlated electrons, with controllable errors, is a central challenge in quantum many-body physics. Here we describe a framework for a general multiscale method based on embedding a self-energy of a strongly correlated subsystem into a self-energy generated by a method able to treat large weakly correlated systems approximately. As an example, we present the embedding of an exact diagonalization self-energy into a self-energy generated from self-consistent second-order perturbation theory. Using a quantum impurity model, generated from a cluster dynamical mean field approximation to the two-dimensional Hubbard model, as a benchmark, we illustrate that our method allows us to obtain accurate results at a fraction of the cost of typical Monte Carlo calculations. We test the method in multiple regimes of interaction strengths and dopings of the model. The general embedding framework we present avoids difficulties such as double counting corrections, frequency-dependent interactions, or vertex functions. As it is solely formulated at the level of the single-particle Green's function, it provides a promising route for the simulation of realistic materials that are currently difficult to study with other methods. [ABSTRACT FROM AUTHOR]

We compute the two-particle quantities relevant for superconducting correlations in the two-dimensional Hubbard model within the dynamical cluster approximation. In the normal state we identify the parameter regime in density, interaction, and second-nearest-neighbor hopping strength that maximizes the dx2-y2superconducting transition temperature. We find in all cases that the optimal transition temperature occurs at intermediate coupling strength, and is suppressed at strong and weak interaction strengths. Similarly, superconducting fluctuations are strongest at intermediate doping and suppressed towards large doping and half filling. We find a change in sign of the vertex contributions to dxy superconductivity from repulsive near half filling to attractive at large doping, p -wave superconductivity is not found at the parameters we study, and s-wave contributions are always repulsive. For negative second-nearest-neighbor hopping the optimal transition temperature shifts towards the electron-doped side in opposition to the van Hove singularity, which moves towards hole doping. We surmise that an increase of the local interaction of the electron-doped compounds would increase Tc. [ABSTRACT FROM AUTHOR]

Heterostructures of transition-metal oxides have emerged as a new route to engineer electronic systems with desired functionalities. Motivated by these developments, we study a two-orbital Hubbard model in a thin-film geometry confined along the cubic [001] direction using the dynamical mean-field theory. We contrast the results of two approximate impurity solvers (exact diagonalization and one-crossing approximation) to the results of the numerically exact continuous-time quantum Monte Carlo solver. Consistent with earlier studies, we find that the one-crossing approximation performs well in the insulating regime, while the advantage of the exact-diagonalization-based solver is more pronounced in the metallic regime. We then investigate various aspects of strongly correlated eg-orbital systems in thin-film geometries. In particular, we show how the interfacial orbital polarization dies off quickly a few layers from the interface and how the film thickness affects the location of the interaction-driven Mott transition. In addition, we explore the changes in the electronic structure with varying carrier concentration and identify large variations of the orbital polarization in the strongly correlated regime. [ABSTRACT FROM AUTHOR]

We study the continuous bandwidth-controlled Mott transition in the two-dimensional single-band Hubbard model with a focus on the critical scaling behavior of charge and spin degrees of freedom. Using plaquette cluster dynamical mean-field theory, we find charge and spin criticality consistent with experimental results for organic conductors. In particular, the charge degree of freedom calculated via the local density of states at the Fermi level shows a smoother transition than expected for the Ising universality class and in single-site dynamical mean-field theory, revealing the importance of short-ranged nonlocal correlations in two spatial dimensions. The spin criticality obtained from the local spin susceptibility agrees quantitatively with nuclear magnetic resonance measurements of the spin-lattice relaxation rate. [ABSTRACT FROM AUTHOR]

The two-dimensional Hubbard model exhibits superconductivity with d-wave symmetry even at half-filling in the presence of a next-nearest neighbor hopping. Using plaquette cluster dynamical mean-field theory with a continuous-time quantum Monte Carlo impurity solver, we reveal the non-Fermi liquid character of the metallic phase in proximity to the superconducting state. Specifically, the low-frequency scattering rate for momenta near (π,0) varies nonmonotonically at low temperatures, and the dc conductivity is T linear at elevated temperatures with an upturn upon cooling. Evidence is provided that pairing fluctuations dominate the normal-conducting state even considerably above the superconducting transition temperature. [ABSTRACT FROM AUTHOR]

Simulations of finite temperature quantum systems provide imaginary frequency Green's functions that correspond one to one to experimentally measurable real-frequency spectral functions. However, due to the bad conditioning of the continuation transform from imaginary to real frequencies, established methods tend to either wash out spectral features at high frequencies or produce spectral functions with unphysical negative parts. Here, we show that explicitly respecting the analytic "Nevanlinna" structure of the Green's function leads to intrinsically positive and normalized spectral functions, and we present a continued fraction expansion that yields all possible functions consistent with the analytic structure. Application to synthetic trial data shows that sharp, smooth, and multipeak data is resolved accurately. Application to the band structure of silicon demonstrates that high energy features are resolved precisely. Continuations in a realistic correlated setup reveal additional features that were previously unresolved. By substantially increasing the resolution of real frequency calculations our work overcomes one of the main limitations of finite-temperature quantum simulations. [ABSTRACT FROM AUTHOR]

We study the half-filled extended Hubbard model on a two-dimensional square lattice using cluster dynamical mean-field theory on clusters of size 8-20. We show that the model exhibits metallic, Mott-insulating, and charge-ordered phases, and determine the location of the charge-ordering phase-transition line and the properties of the phases as a function of temperature, local interaction, and nearest-neighbor interaction. We find strong nonlocal correlations outside the charge-ordered phase and a pronounced screening effect in the vicinity of the phase transition, where nonlocal interactions push the system towards metallic behavior. In contrast, correlations in the charge-ordered phase are mostly local to the unit cell. Finally, we demonstrate how strong nonlocal electron-electron interactions can increase electron mobility by turning a charge-ordered insulator into a metal. We analyze finite-size effects and the convergence of our data to the thermodynamic limit. Control of all sources of errors allows us to assess the regime of applicability of simpler approximation schemes for systems with nonlocal interactions. [ABSTRACT FROM AUTHOR]

We introduce a numerical algorithm to stochastically sample the dual fermion perturbation series around the dynamical mean field theory, generating all topologies of two-particle interaction vertices. We show results in the weak and strong coupling regime of the half-filled Hubbard model in two dimensions, illustrating that the method converges quickly where dynamical mean field theory is a good approximation, and show that corrections are large in the strong correlation regime at intermediate interaction. The fast convergence of dual corrections to dynamical mean field results illustrates the power of the approach and opens a practical avenue towards the systematic inclusion of nonlocal correlations in correlated materials simulations. An analysis of the frequency scale shows that only low-frequency propagators contribute substantially to the diagrams, putting the inclusion of higher order vertices within reach. [ABSTRACT FROM AUTHOR]

We examine the dynamics of a correlated quantum dot in the mixed valence regime. We perform numerically exact calculations 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 equilibration times and saturates upon the decrease of temperature at all times, indicating Kondo behavior both in the transient regime and in the steady state. The time-dependent current saturation temperature connects the equilibrium Kondo temperature to a substantially increased value at voltages outside of the linear response. These signatures are directly observable by experiments in the time domain. [ABSTRACT FROM AUTHOR]

We examine the nature of the transition to the antiferromagnetically ordered state in the half-filled three-dimensional Hubbard model using the dual-fermion multiscale approach. Consistent with analytics, in the weak-coupling regime we find that spin-flip excitations across the Fermi surface are important and that the strong-coupling regime is described by Heisenberg physics. In the intermediate-interaction, strong-correlation regime we find aspects of both local and nonlocal correlations. We analyze the critical exponents of the transition in the strong-coupling regime and find them to be consistent with Heisenberg physics down to an interaction of U/t = 10. [ABSTRACT FROM AUTHOR]

We study the anisotropic 3D Hubbard model with increased nearest-neighbor tunneling amplitudes along one direction using the dynamical cluster approximation and compare the results to a quantum simulation experiment of ultracold fermions in an optical lattice. We find that the short-range spin correlations are significantly enhanced in the direction with stronger tunneling amplitudes. Our results agree with the experimental observations and show that the experimental temperature is lower than the strong tunneling amplitude. We characterize the system by examining the spin correlations beyond neighboring sites and determine the distribution of density, entropy, and spin correlation in the trapped system. We furthermore investigate the dependence of the critical entropy at the Néel transition on anisotropy. [ABSTRACT FROM AUTHOR]

We present two methods for computing two-time correlation functions or Green's functions from real-time bold-line continuous-time quantum Monte Carlo. One method is a formally exact generalized auxiliary lead formalism by which spectral properties may be obtained from single-time observables. The other involves the evaluation of diagrams contributing to two-time observables directly on the Keldysh contour. Additionally, we provide a detailed description of the bold-line Monte Carlo method. Our methods are general and numerically exact, and able to reliably resolve high-energy features such as band edges. We compare the spectral functions obtained from real-time methods to analytically continued spectral functions obtained from imaginary-time Monte Carlo, thus probing the limits of analytic continuation. [ABSTRACT FROM AUTHOR]

We present real-time inchworm quantum Monte Carlo results for single-site dynamical mean field theory on an infinite coordination number Bethe lattice. Our numerically exact results are obtained on the L-shaped Keldysh contour and, being evaluated in real time, avoid the analytic continuation issues typically encountered in Monte Carlo calculations. Our results show that inchworm Monte Carlo methods have now reached a state where they can be used as dynamical mean field impurity solvers and the dynamical sign problem can be overcome. As nonequilibrium problems can be simulated at the same cost, we envisage the main use of these methods as dynamical mean field solvers for time-dependent problems far from equilibrium. [ABSTRACT FROM AUTHOR]

We present a theoretical framework and implementation details for self-energy embedding theory (SEET) with the GW approximation for the treatment of weakly correlated degrees of freedom and configuration interactions solver for handing the strongly correlated degrees. On a series of molecular examples, for which the exact results are known within a given basis, we demonstrate that SEET(CI/GW) is a systematically improvable and well controlled method capable of giving accurate results and well behaved causal self-energies and Green's functions. We compare the theoretical framework of SEET(CI/GW) to that of GW+DMFT and comment on differences between these to approaches that aim to treat both strongly and weakly correlated simultaneously. [ABSTRACT FROM AUTHOR]

We generalize the recently developed inchworm quantum Monte Carlo method to the full Keldysh contour with forward, backward, and equilibrium branches to describe the dynamics of strongly correlated impurity problems with time-dependent parameters. 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 then 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. [ABSTRACT FROM AUTHOR]

Correlated electron systems may give rise to multiple effective interactions whose combined impact on quasiparticle properties can be difficult to disentangle. We introduce an unambiguous decomposition of the electronic self-energy which allows us to quantify the contributions of various effective interactions simultaneously. We use this tool to revisit the hole-doped Hubbard model within the dynamical cluster approximation, where commonly spin fluctuations are considered to be the origin of the pseudogap. While our fluctuation decomposition confirms that spin fluctuations indeed suppress antinodal electronic spectral weight, we show that they alone cannot capture the pseudogap self-energy quantitatively. Nonlocal multiboson Feynman diagrams yield substantial contributions and are needed for a quantitative description of the pseudogap.

A description of long-lived photodoped states in Mott insulators is challenging, as it needs to address exponentially separated timescales. We demonstrate how properties of such states can be computed using numerically exact steady state techniques, in particular, the quantum Monte Carlo algorithm, by using a time-local ansatz for the distribution function with separate Fermi functions for the electron and hole quasiparticles. The simulations show that the Mott gap remains robust to large photodoping, and the photodoped state has hole and electron quasiparticles with strongly renormalized properties.

Response functions of quantum systems, such as electron Green's functions, magnetic, or charge susceptibilities, describe the response of a system to an external perturbation. They are the central objects of interest in field theories and quantum computing and measured directly in experiment. Response functions are intrinsically causal. In equilibrium and steady-state systems, they correspond to a positive spectral function in the frequency domain. Since response functions define an inner product on a Hilbert space and thereby induce a positive definite function, the properties of this function can be used to reduce noise in measured data and, in equilibrium and steady state, to construct positive definite extensions for data known on finite time intervals, which are then guaranteed to correspond to positive spectra.

We study bulk particle transport in a Fermi-Hubbard model on an infinite-dimensional Bethe lattice, driven by a constant electric field. Previous numerical studies showed that one dimensional analogs of this system exhibit a breakdown of diffusion due to Stark many-body localization at least up to time that scales exponentially with the system size. Here, we consider systems initially in a spin density wave state using a combination of numerically exact and approximate techniques. We show that for sufficiently weak electric fields, the wave's momentum component decays exponentially with time in a way consistent with normal diffusion. By studying different wavelengths, we extract the dynamical exponent and the generalized diffusion coefficient at each field strength. Interestingly, we find a nonmonotonic dependence of the dynamical exponent on the electric field. As the field increases toward a critical value proportional to the Hubbard interaction strength, transport slows down, becoming subdiffusive. At large interaction strengths, however, transport speeds up again with increasing field, exhibiting superdiffusive characteristics when the electric field is comparable to the interaction strength. Eventually, at the large field limit, localization occurs and the current through the system is suppressed.

Simulations of finite temperature quantum systems provide imaginary frequency Green's functions that correspond one to one to experimentally measurable real-frequency spectral functions. However, due to the bad conditioning of the continuation transform from imaginary to real frequencies, established methods tend to either wash out spectral features at high frequencies or produce spectral functions with unphysical negative parts. Here, we show that explicitly respecting the analytic "Nevanlinna" structure of the Green's function leads to intrinsically positive and normalized spectral functions, and we present a continued fraction expansion that yields all possible functions consistent with the analytic structure. Application to synthetic trial data shows that sharp, smooth, and multipeak data is resolved accurately. Application to the band structure of silicon demonstrates that high energy features are resolved precisely. Continuations in a realistic correlated setup reveal additional features that were previously unresolved. By substantially increasing the resolution of real frequency calculations our work overcomes one of the main limitations of finite-temperature quantum simulations.

Current nonequilibrium Monte Carlo methods suffer from a dynamical sign problem that makes simulating real-time dynamics for long times exponentially hard. We propose a new "inchworm algorithm," based on iteratively reusing information obtained in previous steps to extend the propagation to longer times. The algorithm largely overcomes the dynamical sign problem, changing the scaling from exponential to quadratic. We use the method to solve the Anderson impurity model in the Kondo and mixed valence regimes, obtaining results both for quenches and for spin dynamics in the presence of an oscillatory magnetic field.

We compute the two-particle quantities relevant for superconducting correlations in the two-dimensional Hubbard model within the dynamical cluster approximation. In the normal state we identify the parameter regime in density, interaction, and second-nearest-neighbor hopping strength that maximizes the d_{x^{2}-y^{2}} superconducting transition temperature. We find in all cases that the optimal transition temperature occurs at intermediate coupling strength, and is suppressed at strong and weak interaction strengths. Similarly, superconducting fluctuations are strongest at intermediate doping and suppressed towards large doping and half filling. We find a change in sign of the vertex contributions to d_{xy} superconductivity from repulsive near half filling to attractive at large doping. p-wave superconductivity is not found at the parameters we study, and s-wave contributions are always repulsive. For negative second-nearest-neighbor hopping the optimal transition temperature shifts towards the electron-doped side in opposition to the van Hove singularity, which moves towards hole doping. We surmise that an increase of the local interaction of the electron-doped compounds would increase T_{c}.

We study the thermodynamic properties of the 3D Hubbard model for temperatures down to the Néel temperature by using cluster dynamical mean-field theory. In particular, we calculate the energy, entropy, density, double occupancy, and nearest-neighbor spin correlations as a function of chemical potential, temperature, and repulsion strength. To make contact with cold-gas experiments, we also compute properties of the system subject to an external trap in the local density approximation. We find that an entropy per particle S/N ≈ 0.65(6) at U/t = 8 is sufficient to achieve a Néel state in the center of the trap, substantially higher than the entropy required in a homogeneous system. Precursors to antiferromagnetism can clearly be observed in nearest-neighbor spin correlators.

We present the effective action and self-consistency equations for the bosonic dynamical mean field approximation to the bosonic Hubbard model and show that it provides remarkably accurate phase diagrams and correlation functions. To solve the bosonic dynamical mean field equations, we use a continuous-time Monte Carlo method for bosonic impurity models based on a diagrammatic expansion in the hybridization and condensate coupling. This method is readily generalized to bosonic mixtures, spinful bosons, and Bose-Fermi mixtures.

A single-site dynamical mean-field study of a three band model with the rotationally invariant interactions appropriate to the t_(2g) levels of a transition metal oxide reveals a quantum phase transition between a paramagnetic metallic phase and an incoherent metallic phase with frozen moments. The Mott transitions occurring at electron densities n=2, 3 per site take place inside the frozen moment phase. The critical line separating the two phases is characterized by a self-energy with the frequency dependence Sigma(omega) approximately sqrt[omega] and a broad quantum critical regime. The findings are discussed in the context of the power law observed in the optical conductivity of SrRuO3.