Search

We consider a class of Hubbard-Stratonovich transformations suitable for treating Hubbard interactions in the context of quantum Monte Carlo simulations. A tunable parameter $p$ allows us to continuously vary from a discrete Ising auxiliary field ($p=\infty$) to a compact auxiliary field that couples to electrons sinusoidally ($p=0$). In tests on the single-band Hubbard model, we find that the severity of the sign problem decreases systematically with increasing $p$. Selecting $p$ finite, however, enables continuous sampling methods like the Langevin or Hamiltonian Monte Carlo methods. We explore the tradeoffs between various simulation methods through numerical benchmarks.

The high-temperature superconducting cuprates are governed by intertwined spin, charge, and superconducting orders. While various state-of-the-art numerical methods have demonstrated that these phases also manifest themselves in doped Hubbard models, they differ on which is the actual ground state. Finite-cluster methods typically indicate that stripe order dominates, while embedded quantum-cluster methods, which access the thermodynamic limit by treating long-range correlations with a dynamical mean field, conclude that superconductivity does. Here, we report the observation of fluctuating spin and charge stripes in the doped single-band Hubbard model using a quantum Monte Carlo dynamical cluster approximation (DCA) method. By resolving both the fluctuating spin and charge orders using DCA, we demonstrate that they survive in the doped Hubbard model in the thermodynamic limit. This discovery also provides an opportunity to study the influence of fluctuating stripe correlations on the model's pairing correlations within a unified numerical framework. Using this approach, we also find evidence for pair-density-wave correlations whose strength is correlated with that of the stripes., Proceedings of the National Academy of Sciences of the United States of America, 119 (7), ISSN:0027-8424, ISSN:1091-6490

In ultrathin films of FeSe grown on SrTiO3 (FeSe/STO), the superconducting transition temperature Tc is increased by almost an order of magnitude, raising questions on the pairing mechanism. As in other superconductors, antiferromagnetic spin fluctuations have been proposed to mediate SC making it essential to study the evolution of the spin dynamics of FeSe from the bulk to the ultrathin limit. Here, we investigate the spin excitations in bulk and monolayer FeSe/STO using resonant inelastic x-ray scattering (RIXS) and quantum Monte Carlo (QMC) calculations. Despite the absence of long-range magnetic order, bulk FeSe displays dispersive magnetic excitations reminiscent of other Fe-pnictides. Conversely, the spin excitations in FeSe/STO are gapped, dispersionless, and significantly hardened relative to its bulk counterpart. By comparing our RIXS results with simulations of a bilayer Hubbard model, we connect the evolution of the spin excitations to the Fermiology of the two systems revealing a remarkable reconfiguration of spin excitations in FeSe/STO, essential to understand the role of spin fluctuations in the pairing mechanism., Here, Pelliciari et al. present resonant inelastic X-ray scattering on monolayer samples of unconventional superconductor FeSe, finding evidence for gapped and dispersionless spin excitations. These experiments are very difficult due to the extremely small scattering volume of the FeSe monolayer.

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., Comment: 5 pages, 5 figures

We show that, by an appropriate choice of auxiliary fields and exact integration of the phonon degrees of freedom, it is possible to define a 'sign-free' path integral for the so called Hubbard-Holstein model at half-filling. We use a statistical method, based on an accelerated and efficient Langevin dynamics, for evaluating all relevant correlation functions of the model. Preliminary calculations at $U/t=4$ and $U/t=1$, for $\omega_0/t=1$, indicate a quite extended region around $U \simeq {g^2 \over \omega_0}$ without either antiferromagnetic or charge-density-wave orders, separating two quantum critical points at zero temperature. The elimination of the sign problem in a model without explicit particle-hole symmetry may open new perspectives for strongly correlated models, even away from the purely attractive or particle-hole symmetric cases., Comment: 5 pages, 2 figures

By using variational Monte Carlo and auxiliary-field quantum Monte Carlo methods, we perform an accurate finite-size scaling of the $s$-wave superconducting order parameter and the pairing correlations for the negative-$U$ Hubbard model at zero temperature in the square lattice. We show that the twist-averaged boundary conditions (TABCs) are extremely important to control finite-size effects and to achieve smooth and accurate extrapolations to the thermodynamic limit. We also show that TABCs are much more efficient in the grand-canonical ensemble rather than in the standard canonical ensemble with fixed number of electrons. The superconducting order parameter as a function of the doping is presented for several values of $|U|/t$ and is found to be significantly smaller than the mean-field BCS estimate already for moderate couplings. This reduction is understood by a variational ansatz able to describe the low-energy behavior of the superconducting phase by means of a suitably chosen Jastrow factor including long-range density-density correlations.

By using variational wave functions and quantum Monte Carlo techniques, we investigate the interplay between electron-electron and electron-phonon interactions in the two-dimensional Hubbard-Holstein model. Here, the ground-state phase diagram is triggered by several energy scales, i.e., the electron hopping $t$, the on-site electron-electron interaction $U$, the phonon energy $\omega_0$, and the electron-phonon coupling $g$. At half filling, the ground state is an antiferromagnetic insulator for $U \gtrsim 2g^2/\omega_0$, while it is a charge-density-wave (or bi-polaronic) insulator for $U \lesssim 2g^2/\omega_0$. In addition to these phases, we find a superconducting phase that intrudes between them. For $\omega_0/t=1$, superconductivity emerges when both $U/t$ and $2g^2/t\omega_0$ are small; then, by increasing the value of the phonon energy $\omega_0$, it extends along the transition line between antiferromagnetic and charge-density-wave insulators. Away from half filling, phase separation occurs when doping the charge-density-wave insulator, while a uniform (superconducting) ground state is found when doping the superconducting phase. In the analysis of finite-size effects, it is extremely important to average over twisted boundary conditions, especially in the weak-coupling limit and in the doped case., Comment: 11 pages, 10 figures

We present a systematic investigation of the special twist method introduced by Rajagopal $\textit{et al.}$ [ Phys. Rev. B 51, 10591 (1995) ] for reducing finite-size effects in correlated calculations of periodic extended systems with Coulomb interactions and Fermi statistics. We propose a procedure for finding special twist values which, at variance with previous applications of this method, reproduce the energy of the mean-field infinite-size limit solution within an adjustable (arbitrarily small) numerical error. This choice of the special twist is shown to be the most accurate single-twist solution for curing one-body finite-size effects in correlated calculations. For these reasons we dubbed our procedure "exact special twist" (EST). EST only needs a fully converged independent-particles or mean-field calculation within the primitive cell and a simple fit to find the special twist along a specific direction in the Brillouin zone. We first assess the performances of EST in a simple correlated model such as the 3D electron gas. Afterwards, we test its efficiency within $\textit{ab initio}$ quantum Monte Carlo simulations of metallic elements of increasing complexity. We show that EST displays an overall good performance in reducing finite-size errors comparable to the widely used twist average technique but at a much lower computational cost, since it involves the evaluation of just one wavefunction. We also demonstrate that the EST method shows similar performances in the calculation of correlation functions, such as the ionic forces for structural relaxation and the pair radial distribution function in liquid hydrogen. Our conclusions point to the usefulness of EST for correlated supercell calculations, our method will be particularly relevant when the physical problem under consideration requires large periodic cells., 15 pages, 8 figures