Dynamical Cluster Approximation Research Articles

Spin-singlet topological superconductivity in the attractive Rashba-Hubbard model

Fully gapped, spin-singlet superconductors with antisymmetric spin-orbit coupling in a Zeeman magnetic field provide a promising route to realize superconducting states with non-Abelian topological order and therefore fault-tolerant quantum computation. Here we use a quantum Monte Carlo dynamical cluster approximation to study the superconducting properties of a doped two-dimensional attractive Hubbard model with Rashba spin-orbit coupling in a Zeeman magnetic field. We generally find that the Rashba coupling has a beneficial effect towards $s$-wave superconductivity. In the presence of a finite Zeeman field, when superconductivity is suppressed by Pauli pair breaking, the Rashba coupling counteracts the spin imbalance created by the Zeeman field by mixing the spins, and thus restores superconductivity at finite temperatures. We show that this favorable effect of the spin-orbit coupling is traced to a spin-flip driven enhancement of the amplitude for the propagation of a pair of electrons in time-reversed states. Moreover, by inspecting the Fermi surface of the interacting model, we show that for sufficiently large Rashba coupling and Zeeman field, the superconducting state is expected to be topologically nontrivial.

Robust charge-density-wave correlations in the electron-doped single-band Hubbard model

There is growing evidence that the hole-doped single-band Hubbard and t − J models do not have a superconducting ground state reflective of the high-temperature cuprate superconductors but instead have striped spin- and charge-ordered ground states. Nevertheless, it is proposed that these models may still provide an effective low-energy model for electron-doped materials. Here we study the finite temperature spin and charge correlations in the electron-doped Hubbard model using quantum Monte Carlo dynamical cluster approximation calculations and contrast their behavior with those found on the hole-doped side of the phase diagram. We find evidence for a charge modulation with both checkerboard and unidirectional components decoupled from any spin-density modulations. These correlations are inconsistent with a weak-coupling description based on Fermi surface nesting, and their doping dependence agrees qualitatively with resonant inelastic x-ray scattering measurements. Our results provide evidence that the single-band Hubbard model describes the electron-doped cuprates.

Mott transition, Widom line, and pseudogap in the half-filled triangular lattice Hubbard model

The Mott transition is observed experimentally in materials that are magnetically frustrated so that long-range order does not hide the Mott transition at finite temperature. The Hubbard model on the triangular lattice at half-filling is a paradigmatic model to study the interplay of interactions and frustration on the normal-state phase diagram. We use the dynamical cluster approximation with continuous time auxiliary field quantum Monte Carlo to solve this model for 1, 4, 6, 12, and 16 site clusters with detailed analysis performed for the 6 site cluster. We show that a) for every cluster there is an inflection point in the double occupancy as a function of interaction, defining a Widom line that extends above the critical point of the first-order Mott transition; b) the presence of this line and the cluster size dependence argue for the observability of the Mott transition at finite temperature in the thermodynamic limit; c) the loss of spectral weight in the metal to Mott insulator transition as a function of temperature and for strong interactions is momentum dependent, the hallmark of a pseudogap. That pseudogap spans a large region of the phase diagram near the Mott transition.

Non-Fermi Liquid Behavior in the Three-Dimensional Hubbard Model

We present a numerical study on the non-Fermi liquid behavior of a three-dimensional strongly correlated system. The Hubbard model in a simple cubic lattice is simulated by the dynamical cluster approximation; in particular, the quasi-particle weight is calculated at finite dopings for a range of temperatures. By fitting the quasi-particle weight to the marginal Fermi liquid form at finite doping near the putative quantum critical point, we find evidence of a separatrix between Fermi liquid and non-Fermi liquid regions. Our results suggest that a marginal Fermi liquid and possibly a quantum critical point exist in the non-symmetry broken solution of the three-dimensional interacting electron systems. We also calculate the spectral function, close to the half-filling, and we obtain evidence of pseudogap.

Single- and two-particle finite size effects in interacting lattice systems

Simulations of extended quantum systems are typically performed by extrapolating results of a sequence of finite-system-size simulations to the thermodynamic limit. In the quantum Monte Carlo community, twist-averaging was pioneered as an efficient strategy to eliminate one-body finite size effects. In the dynamical mean field community, cluster generalizations of the dynamical mean field theory were formulated to study systems with nonlocal correlations. In this work, we put the twist-averaging and the dynamical cluster approximation variant of the dynamical mean field theory onto equal footing, discuss commonalities and differences, and compare results from both techniques to the standard periodic boundary technique. At the example of Hubbard-type models with local, short-range and Yukawa-like longer range interactions we show that all methods converge to the same limit, but that the convergence speed differs in practice. We show that embedding theories are an effective tool for managing both one-body and two-body finite size effects, in particular if interactions are averaged over twist angles.

Quantifying the role of antiferromagnetic fluctuations in the superconductivity of the doped Hubbard model

We study the contribution of the electron-spin fluctuation coupling to the superconducting state of the two dimensional Hubbard model within dynamical cluster approximation (DCA) using a numerical exact continuous time Monte Carlo solver. By analyzing the frequency dependence of the self energy, we show that only about half of the superconductivity can be attributed to a "pairing glue" arising from treating spin fluctuations as a pairing boson in the standard one-loop theory.

Enhancing Tc in a composite superconductor/metal bilayer system: A dynamical cluster approximation study

It has been proposed that the superconducting transition temperature ${T}_{\text{c}}$ of an unconventional superconductor with a large pairing scale but strong phase fluctuations can be enhanced by coupling it to a metal. However, the general efficacy of this approach across different parameter regimes remains an open question. Using the dynamical cluster approximation, we study this question in a system composed of an attractive Hubbard layer in the intermediate coupling regime, where the magnitude of the attractive Coulomb interaction $|U|$ is slightly larger than the bandwidth $W$, hybridized with a noninteracting metallic layer. We find that while the superconducting transition becomes more mean-field-like with increasing interlayer hopping, the superconducting transition temperature ${T}_{\text{c}}$ exhibits a nonmonotonic dependence on the strength of the hybridization ${t}_{\ensuremath{\perp}}$. This behavior arises from a reduction of the effective pairing interaction in the correlated layer that outcompetes the growth in the intrinsic pair-field susceptibility induced by the coupling to the metallic layer. We find that the largest ${T}_{\text{c}}$ inferred here for the composite system is comparable to the maximum value currently estimated for the isolated negative-$U$ Hubbard model.

Superconductivity in the bilayer Hubbard model: Two Fermi surfaces are better than one

Fully occupied or unoccupied bands in a solid are often considered inert and irrelevant to a material's low-energy properties. But the discovery of enhanced superconductivity in heavily electron-doped FeSe-derived superconductors poses questions about the possible role of incipient bands (those laying close to but not crossing the Fermi level) in pairing. To answer this question, researchers have studied pairing correlations in the bilayer Hubbard model, which has an incipient band for large interlayer hopping ${t}_{\ensuremath{\perp}}$, using many-body perturbation theory and variational methods. They have generally found that superconductivity is enhanced as one of the bands approaches the Lifshitz transition and even when it becomes incipient. Here we address this question using the nonperturbative quantum Monte Carlo (QMC) dynamical cluster approximation (DCA) to study the bilayer Hubbard model's pairing correlations. We find that the model has robust ${s}^{\ifmmode\pm\else\textpm\fi{}}$ pairing correlations in the large ${t}_{\ensuremath{\perp}}$ limit, which can become stronger as one band is made incipient. While this behavior is linked to changes in the effective interaction, we further find that it is counteracted by a suppression of the intrinsic pair-field susceptibility and does not translate to an increased ${T}_{c}$. Our results demonstrate that the highest achievable transition temperatures in the bilayer Hubbard model occur when the system has two bands crossing the Fermi level.

Dynamical cluster approximation study of electron localization in the extended Hubbard model

We perform a detailed study of the phase transitions and mechanisms of electron localization in the extended Hubbard model using the dynamical cluster approximation on a $2\times 2$ cluster. We explore the interplay of charge order and Mott physics. We find that a nearest-neighbor Coulomb interaction $V$ causes "screening" effects close to the Mott phase transition, pushing the phase boundary to larger values of $U$. We also demonstrate the different effects of $V$ on correlations in metallic and insulating regimes and document the different correlation aspects of charge order and Mott states.

Pairing correlations in the cuprates: A numerical study of the three-band Hubbard model

We study the three-band Hubbard model for the copper oxide plane of the high-temperature superconducting cuprates using determinant quantum Monte Carlo and the dynamical cluster approximation (DCA) and provide a comprehensive view of the pairing correlations in this model using these methods. Specifically, we compute the pair-field susceptibility and study its dependence on temperature, doping, interaction strength, and charge-transfer energy. Using the DCA, we also solve the Bethe-Salpeter equation for the two-particle Green's function in the particle-particle channel to determine the transition temperature to the superconducting phase on smaller clusters. Our calculations reproduce many aspects of the cuprate phase diagram and indicate that there is an "optimal" value of the charge-transfer energy for the model where $T_c$ is largest. These results have implications for our understanding of superconductivity in both the cuprates and other doped charge-transfer insulators.

Spin density wave order in interacting type-I and type-II Weyl semimetals

Weyl semimetals, featuring massless linearly dispersing chiral fermions in three dimensions, provide an excellent platform for studying the interplay of electronic interactions and topology, and exploring new correlated states of matter. Here, we examine the effect of a local repulsive interaction on an inversion-symmetry breaking Weyl semimetal model, using cluster dynamical mean field theory and variational cluster approximation methods. Our analysis reveals a continuous transition from the gapless Weyl semimetal phase to a gapped spin density wave ordered phase at a critical value of the interaction, which is determined by the band structure parameters. Further, we introduce a finite tilt in the linear dispersion and examine the corresponding behavior for a type-II Weyl semimetal model, where the critical interaction strength is found to be significantly diminished, indicating a greater susceptibility towards interactions. The behavior of different physical quantities, such as the double occupancy, the spectral function and the Berry curvature associated with the Weyl nodes are obtained in both the semimetallic and the magnetically ordered states. Finally, we provide an interaction-induced phase diagram for the Weyl semimetal model, as a function of the tilt parameter.

Effects of frustration on the nonequilibrium dynamics of photoexcited lattice systems

We theoretically investigate the effects of the lattice geometry on the nonequilibrium dynamics of photo-excited carriers in a half-filled two-dimensional Hubbard model. Using a nonequilibrium generalization of the dynamical cluster approximation, we compare the relaxation dynamics in lattices which interpolate between the triangular lattice and square lattice configuration and thus reveal the role of the geometric frustration in these strongly correlated nonequilibrium systems. In particular, we show that the cooling effect resulting from the disordering of the spin background is less effective in the triangular case because of the frustration. This manifests itself in a longer relaxation time of the photo-doped population, as measured by the time-resolved photo-emission signal, and a higher effective temperature of the photo-doped carriers in the non-thermal steady state after the intra-Hubbard-band thermalization.

Effect of short-ranged spatial correlations on the Anderson localization of phonons in mass-disordered systems

We investigate the effect of spatially correlated disorder on the Anderson transition of phonons in three dimensions using a Green’s function-based approach, namely, the typical medium dynamical cluster approximation, in mass-disordered systems. We numerically demonstrate that correlated disorder with pairwise correlations mitigates the localization of the vibrational modes. A correlation-driven localization–delocalization transition can emerge in a three-dimensional disordered system with an increase in the strength of correlations.

Overdoped end of the cuprate phase diagram

Studying the disappearance of superconductivity at the end of the overdoped region of the cuprate phase diagram offers a different approach for investigating the interaction which is responsible for pairing in these materials. In the underdoped region this question is complicated by the presence of charge and stripe ordered phases as well as the pseudogap. In the overdoped region the situation appears simpler with only a normal phase, a superconducting phase and impurity scattering. Here, for the overdoped region, we report the results of a combined dynamic cluster approximation (DCA) and a weak Born impurity scattering calculation for a $t-t'-U$ Hubbard model. We find that a decrease in the $d$-wave pairing strength of the two-particle scattering vertex is closely coupled to changes in the momentum and frequency structure of the magnetic spin fluctuations as the system is overdoped. Treating the impurity scattering within a disordered BCS $d$-wave approximation, we see how the combined effects of the decreasing $d$-wave pairing strength and weak impurity scattering lead to the end of the $T_c$ dome.

Continuous momentum dependence in the dynamical cluster approximation

The dynamical cluster approximation (DCA) is a quantum cluster extension to the single-site dynamical mean-field theory that incorporates spatially nonlocal dynamic correlations systematically and nonperturbatively. The DCA$^+$ algorithm addresses the cluster shape dependence of the DCA and improves the convergence with cluster size by introducing a lattice self-energy with continuous momentum dependence. However, we show that the DCA$^+$ algorithm is plagued by a fundamental problem when its self-consistency equations are formulated using the bare Green's function of the cluster. This problem is most severe in the strongly correlated regime at low doping, where the DCA$^+$ self-energy becomes overly metallic and local, and persists to cluster sizes where the standard DCA has long converged. In view of the failure of the DCA$^+$ algorithm, we propose to complement DCA simulations with a post-interpolation procedure for single-particle and two-particle correlation functions to preserve continuous momentum dependence and the associated benefits in the DCA. We demonstrate the effectiveness of this practical approach with results for the half-filled and hole-doped two-dimensional Hubbard model.

Short-range charge fluctuations in the two-dimensional Hubbard model

We investigate charge fluctuations in the two-dimensional Hubbard model as a function of doping, interaction strength, next-nearest-neighbor hopping, and temperature within the eight-site dynamical cluster approximation. In the regime of intermediate interaction strengths, we find that $d$ density wave fluctuations, which had previously been postulated using analytical arguments, are present and strong but cannot be interpreted as the cause of the pseudogap in the model, due to their evolution with doping and interaction strength. For all parameters away from half filling, the charge fluctuations investigated, including $d$ density wave fluctuations, are weaker than $d$-wave superconducting fluctuations.

Ferromagnetic spin correlations in the two-dimensional Hubbard model

We analyze the dynamical nearest-neighbor and next-nearest-neighbor spin correlations in the 4-site and 8-site dynamical cluster approximation to the two-dimensional Hubbard model. Focusing on the robustness of these correlations at long imaginary times, we reveal enhanced ferromagnetic correlations on the lattice diagonal, consistent with the emergence of composite spin-1 moments at a temperature scale that essentially coincides with the pseudo-gap temperature $T^*$. We discuss these results in the context of the spin-freezing theory of unconventional superconductivity.

Photoenhanced excitonic correlations in a Mott insulator with nonlocal interactions

We investigate the effect of nonlocal interactions on the photo-doped Mott insulating state of the two-dimensional Hubbard model using a nonequilibrium generalization of the dynamical cluster approximation. In particular, we compare the situation where the excitonic states are lying within the continuum of doublon-holon excitations to a set-up where the excitons appear within the Mott gap. In the first case, the creation of nearest-neighbor doublon-holon pairs by excitations across the Mott gap results in enhanced excitonic correlations, but these excitons quickly decay into uncorrelated doublons and holons. In the second case, photo-excitation results in long-lived excitonic states. While in a low-temperature equilibrium state, excitonic features are usually not evident in single-particle observables such as the photoemission spectrum, we show that the photo-excited nonequilibrium system can exhibit in-gap states associated with the excitons. The comparison with exact-diagonalization results for small clusters allows us to identify the signatures of the excitons in the photo-emission spectrum.

Dual parquet scheme for the two-dimensional Hubbard model: Modeling low-energy physics of high- Tc cuprates with high momentum resolution

We present a method to treat the two-dimensional (2D) Hubbard model for parameter regimes which are relevant for the physics of the high-${T}_{c}$ superconducting cuprates. Unlike previous attempts to attack this problem, our approach takes into account all fluctuations in different channels on equal footing and is able to treat reasonable large lattice sizes up to $32\ifmmode\times\else\texttimes\fi{}32$. This is achieved by the following three-step procedure: (i) We transform the original problem to a representation (dual fermions) in which all purely local correlation effects from the dynamical mean field theory are already considered in the bare propagator and bare interaction of the new problem. (ii) The strong $1/{(i\ensuremath{\nu})}^{2}$ decay of the bare propagator allows us to integrate out all higher Matsubara frequencies aside from the lowest using low-order diagrams. The new effective action depends only on the two lowest Matsubara frequencies which allows us to (iii) apply the two-particle self-consistent parquet formalism, which takes into account the competition between different low-energy bosonic modes in an unbiased way, on much finer momentum grids than usual. In this way, we were able to map out the phase diagram of the 2D Hubbard model as a function of temperature and doping. Consistently with the experimental evidence for hole-doped cuprates and previous dynamical cluster approximation calculations, we find an antiferromagnetic region at low doping and a superconducting dome at higher doping. Our results also support the role of the van Hove singularity as an important ingredient for the high value of ${T}_{c}$ at optimal doping. At small doping, the destruction of antiferromagnetism is accompanied by an increase of charge fluctuations supporting the scenario of a phase-separated state driven by quantum critical fluctuations.

Extended Crossover from a Fermi Liquid to a Quasiantiferromagnet in the Half-Filled 2D Hubbard Model.

The ground state of the Hubbard model with nearest-neighbor hopping on the square lattice at half filling is known to be that of an antiferromagnetic (AFM) band insulator for any on-site repulsion. At finite temperature, the absence of long-range order makes the question of how the interaction-driven insulator is realized nontrivial. We address this problem with controlled accuracy in the thermodynamic limit using self-energy diagrammatic determinant MonteCarlo and dynamical cluster approximation methods and show that development of long-range AFM correlations drives an extended crossover from Fermi liquid to insulating behavior in the parameter regime that precludes a metal-to-insulator transition. The intermediate crossover state is best described as a non-Fermi liquid with a partially gapped Fermi surface.