We develop a Fermi-Bose bootstrap embedding framework for the ground state of interacting electrons coupled to a phonon mean field. The method combines bootstrap embedding for correlated electrons with a self-consistent coherent-state mean-field treatment for phonons. This method models the interacting electron-phonon problem as a system of correlated electrons traveling in a self-consistently specified potential landscape, allowing for efficient treatment of large lattice systems. Convergence of the methods for fragment size and total system size is demonstrated for the one-dimensional Hubbard-Holstein model for up to 350 sites. Finite-size scaling is performed to extrapolate to the infinite system size. Benchmarking against the density matrix renormalization group for a small 8-site system at half- and quarter-filling shows an orders-of-magnitude runtime advantage. The comparison further reveals that the method performs best in regimes dominated by localization, such as the Mott insulating phase and the strong-coupling tiny polaron regime, where the local embedding ansatz is still valid. However, due to the mean-field treatment for phonons, we find limitations of our methods in the weakly coupled delocalized region and at the Peierls transition, where quantum phonon fluctuations and long-range kinetic correlations become substantial.
We develop a Fermi-Bose bootstrap embedding framework for the ground state of interacting electrons coupled to a phonon mean field. The method combines bootstrap embedding for correlated electrons with a self-consistent coherent-state mean-field treatment for phonons. This method models the interacting electron-phonon problem as a system of correlated electrons traveling in a self-consistently specified potential landscape, allowing for efficient treatment of large lattice systems. Convergence of the methods for fragment size and total system size is demonstrated for the one-dimensional Hubbard-Holstein model for up to 350 sites. Finite-size scaling is performed to extrapolate to the infinite system size. Benchmarking against the density matrix renormalization group for a small 8-site system at half- and quarter-filling shows an orders-of-magnitude runtime advantage. The comparison further reveals that the method performs best in regimes dominated by localization, such as the Mott insulating phase and the strong-coupling tiny polaron regime, where the local embedding ansatz is still valid. However, due to the mean-field treatment for phonons, we find limitations of our methods in the weakly coupled delocalized region and at the Peierls transition, where quantum phonon fluctuations and long-range kinetic correlations become substantial.
The Hubbard-Holstein model describes fermions on a discrete lattice, with on-site repulsion between fermions and a coupling to phonons that are localized on sites. Generally, at half-filling, increasing the coupling g to the phonons drives the system towards a Peierls charge density wave state, whereas increasing the electron-electron interaction U drives the fermions into a Mott antiferromagnet. At low g and U, or when doped, the system is metallic. In one dimension, using quantum Monte Carlo simulations, we study the case where fermions have a long-range coupling to phonons, with characteristic range $\ensuremath{\xi}$, interpolating between the Holstein and Fr\"ohlich limits. Without electron-electron interaction, the fermions adopt a Peierls state when the coupling to the phonons is strong enough. This state is destabilized by a small coupling range $\ensuremath{\xi}$ and leads to a collapse of the fermions, and, consequently, phase separation. Increasing interaction U will drive any of these three phases (metallic, Peierls, phase separation) into a Mott insulator phase. The phase separation region is once again present in the $U\ensuremath{\ne}0$ case, even for small values of the coupling range.
We develop a generalization of the ghost-rotationally-invariant slave-boson (ghost-RISB) method that incorporates local phonon modes coupled to arbitrary on-site electronic degrees of freedom, enabling a nonperturbative treatment of electron-electron and electron-phonon interactions within an efficient variational framework. The method extends the ghost-orbital construction to capture dynamical self-energy effects and phonon-induced renormalizations beyond static slave-boson approaches. Benchmarking against dynamical mean-field theory (DMFT) results for the Hubbard-Holstein model, we find excellent quantitative agreement for electron quasiparticle weights and phonon properties across a wide range of coupling strengths, and it captures accurately the competition between electron-electron and electron-phonon interactions. We show that the inclusion of the ghost orbitals is crucial to accurately describe the regime of low-frequency, strongly dynamical, phonons. The extended ghost-RISB achieves this accuracy at a fraction of the computational cost of DMFT, due to its self-consistency rooted in static observables instead of dynamical ones, enabling rapid exploration of correlated electron-phonon phase diagrams. We exploit this advantage to characterize the most demanding regime of strong coupling and adiabatic phonons. Our analysis shows a suppression of the superconducting order parameter, which is interpreted as a Franck-Condon-like reduction of the overlap between the phonon wavefunctions associated with empty and doubly-occupied sites in the bipolaronic regime.
We describe the extension of the density matrix embedding theory framework to coupled interacting fermion-boson systems. This provides a frequency-independent, entanglement embedding formalism to treat bulk fermion-boson problems. We illustrate the concepts within the context of the one-dimensional Hubbard-Holstein model, where the phonon bath states are obtained from the Schmidt decomposition of a self-consistently adjusted coherent state. We benchmark our results against accurate density matrix renormalization group calculations.
Strong correlation poses a difficult problem for electronic structure theory, with computational cost scaling quickly with system size. Fragment embedding is an attractive approach to this problem. By dividing a large complicated system into smaller manageable fragments "embedded" in an approximate description of the rest of the system, we can hope to ameliorate the steep cost of correlated calculations. While appealing, these methods often converge slowly with fragment size because of small errors at the boundary between fragment and bath. We describe a new electronic embedding method, dubbed "Bootstrap Embedding," a self-consistent wavefunction-in-wavefunction embedding theory that uses overlapping fragments to improve the description of fragment edges. We apply this method to the one dimensional Hubbard model and a translationally asymmetric variant, and find that it performs very well for energies and populations. We find Bootstrap Embedding converges rapidly with embedded fragment size, overcoming the surface-area-to-volume-ratio error typical of many embedding methods. We anticipate that this method may lead to a low-scaling, high accuracy treatment of electron correlation in large molecular systems.
We present a density-matrix embedding theory (DMET) study of the one-dimensional Hubbard-Holstein model, which is paradigmatic for the interplay of electron-electron and electron-phonon interactions. Analyzing the single-particle excitation gap, we find a direct Peierls insulator to Mott insulator phase transition in the adiabatic regime of slow phonons in contrast to a rather large intervening metallic phase in the anti-adiabatic regime of fast phonons. We benchmark the DMET results for both on-site energies and excitation gaps against density-matrix renormalization group (DMRG) results and find good agreement of the resulting phase boundaries. We also compare the full quantum treatment of phonons against the standard Born-Oppenheimer (BO) approximation. The BO approximation gives qualitatively similar results to DMET in the adiabatic regime but fails entirely in the anti-adiabatic regime, where BO predicts a sharp direct transition from Mott to Peierls insulator, whereas DMET correctly shows a large intervening metallic phase. This highlights the importance of quantum fluctuations in the phononic degrees of freedom for metallicity in the one-dimensional Hubbard-Holstein model.
This thesis investigates the interplay between electron-electron and electron-phonon interactions in low-dimensional lattice models, specifically within the Hubbard model incorporating Su-Schrieffer-Heeger (SSH) coupling with phonons. SSH interactions, which modulate the electron hopping amplitude through lattice displacements, have recently attracted renewed interest due to their potential relevance in strongly correlated materials. The main objective of this work is to determine whether the inclusion of SSH couplings in the Hubbard model can stabilize superconducting or spin-gapped metallic states, as well as other exotic phases of matter. While the half-filled SSH-Hubbard model has been extensively studied in one dimension and only partially examined in two dimensions, its phase diagram in the doped regime remains largely unexplored. This thesis aims to fill this gap, focusing on the emergence of correlated phases away from half-filling in both one and two spatial dimensions. We use the variational Monte Carlo method with a trial wave function defined in a hybrid electron-phonon Hilbert space. In addition to Jastrow factors that account for correlation effects, the wave function includes an antisymmetric component that captures backflow-inspired electron-phonon correlations. This approach enables us to accurately model both electron-electron and electron-phonon interactions. Moreover, since the variational method avoids the sign problem, we can simulate the system in the doped regime, an area typically inaccessible to exact quantum Monte Carlo methods (QMC), which are limited to half-filling. In one dimension, we observe the stabilization of a Luther-Emery liquid, a metallic phase with gapless charge excitations and a finite spin gap, upon doping a Peierls insulating phase. This finding resolves a long-standing debate regarding its realization in SSH chains. Benchmarking with density-matrix renormalization group results confirms the reliability of our approach. In two dimensions, we investigate the SSH-Hubbard model on the square lattice, obtaining results that are consistent with previous exact QMC studies at half-filling. Upon doping, we observe that SSH phonons can induce s-wave superconductivity in the weakly correlated regime. However, we suggest that the enhancement of d-wave superconductivity by SSH coupling in the strongly correlated regime may significantly depend on the specific implementation of the SSH interaction. In one particular case, we find that d-wave superconductivity is slightly suppressed by SSH coupling. These results are further confirmed by a Hartree-Fock analysis performed in the anti-adiabatic limit. Overall, this thesis highlights the crucial role of lattice effects in correlated electron systems, shedding light on the distinct contributions of SSH-type couplings in one and two dimensions. Our findings provide new insights into the microscopic mechanisms behind exotic metallic phases and set the stage for further studies of phonon-mediated phenomena in correlated systems.