Strongly correlated electron systems host a variety of poorly understood correlations in their high-temperature normal state. Unlike ordered phases defined by order parameters, regions of the normal state are often defined through unconventional properties such as strange metallic transport or spectroscopic pseudogaps. Characterizing the microscopic correlations in the normal state is necessary to elucidate mechanisms that lead to these properties and their connection to ground-state orders. Here we establish the presence of intertwined charge and spin stripes in the strange metal normal state of the Hubbard model using determinant quantum Monte Carlo calculations. The charge and spin density waves constituting the stripes are fluctuating and short ranged; yet they obey a mutual commensurability relation and remain microscopically interlocked, as evidenced through measurements of three-point spin-spin-hole correlation functions. Our findings demonstrate the ability of many-body numerical simulations to unravel the microscopic correlations that define quantum states of matter.
Strongly correlated electron systems host a variety of poorly understood correlations in their high-temperature normal state. Unlike ordered phases defined by order parameters, regions of the normal state are often defined through unconventional properties such as strange metallic transport or spectroscopic pseudogaps. Characterizing the microscopic correlations in the normal state is necessary to elucidate mechanisms that lead to these properties and their connection to ground-state orders. Here we establish the presence of intertwined charge and spin stripes in the strange metal normal state of the Hubbard model using determinant quantum Monte Carlo calculations. The charge and spin density waves constituting the stripes are fluctuating and short ranged; yet they obey a mutual commensurability relation and remain microscopically interlocked, as evidenced through measurements of three-point spin-spin-hole correlation functions. Our findings demonstrate the ability of many-body numerical simulations to unravel the microscopic correlations that define quantum states of matter.
The high-temperature superconducting cuprates host unidirectional spin- and charge-density-wave orders that can intertwine with superconductivity in nontrivial ways. While the charge components of these stripes have now been observed in nearly all cuprate families, their detailed evolution with doping varies across different materials and at high and low temperatures. We address this problem using nonperturbative determinant quantum Monte Carlo calculations for the three-band Hubbard model. Using an efficient implementation, we can resolve the model's fluctuating spin and charge modulations and map their evolution as a function of the charge transfer energy and doping. We find that the incommensurability of the charge modulations is decoupled from the spin modulations and decreases with hole doping, consistent with experimental measurements at high temperatures. These findings support the proposal that the high-temperature charge correlations are distinct from the intertwined stripe order observed at low-temperature and in the single-band Hubbard model.
Upon doping, Mott insulators often exhibit symmetry breaking where charge carriers and their spins organize into patterns known as stripes. For high-transition temperature cuprate superconductors, stripes are widely suspected to exist in a fluctuating form. We used numerically exact determinant quantum Monte Carlo calculations to demonstrate dynamical stripe correlations in the three-band Hubbard model, which represents the local electronic structure of the copper-oxygen plane. Our results, which are robust to varying parameters, cluster size, and boundary conditions, support the interpretation of experimental observations such as the hourglass magnetic dispersion and the Yamada plot of incommensurability versus doping in terms of the physics of fluctuating stripes. These findings provide a different perspective on the intertwined orders emerging from the cuprates' normal state.
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.
The single-band Hubbard model is one of the most extensively studied models in condensed matter physics, giving rise to intertwined spin and charge stripes that coexist with, or lie in the vicinity of, superconductivity in the phase diagram. However, whether the low energy physics of the single-band Hubbard model is fully equivalent to the multi-band (multi-orbital) Emery model remains an unsettled question. While the intertwined stripes and nematicity have been studied in the single-band Hubbard model, a comprehensive picture in the Emery model is lacking. In this paper, we focus on the less investigated intertwined charge stripes using complementary density matrix renormalization group (DMRG) and determinant quantum Monte Carlo (DQMC) techniques. Our ground state DMRG confirms the presence of the oxygen-centered charge stripes at a reduced amplitude in the Emery model parameter regime widely used in the study of intertwined stripes. Close analysis of the oxygen orbital structure of the static charge correlation function from DQMC reveals the charge stripe pattern in real-space, showcasing stronger charge density modulation on p-orbitals pointing along ``the rivers of charge'', consistent with DMRG. For the parameter set with the largest fermion signs, we managed to reach a temperature where the system first demonstrated tendencies to form purely unidirectional spin and charge stripes, and the B1g component becomes dominant in the bond-charge nematic susceptibility. This observation correlates with the doping dependence of the kinetic energy anisotropy, suggesting a close relation between the nematicity and charge stripes in the Emery model.
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 for the first time that they survive in the doped Hubbard model in the thermodynamic limit. This discovery also provides a new opportunity to study the influence of fluctuating stripe correlations on the model's pairing correlations within a unified numerical framework.
We revisit the effective theory for fluctuating spin stripes coupled to a Fermi surface, and consider the parameter regime where a spin nematic phase intervenes between the spin density wave state and the symmetric state. It is shown that adding potential disorder to this theory, which acts as an unconventional type of random-field disorder, naturally gives rise to a phase diagram containing a quantum critical point that is described by the universal theory of strange metals with spatial disorder in both the magnitude and sign of the electron-boson coupling term [A.A. Patel, H. Guo, I. Esterlis and S. Sachdev, Science 381, 790 (2023)]. One difference compared to the original theory, however, is that at non-zero temperatures the disordered spin-stripe model automatically self-averages over the sign of the coupling. We also study the effects of thermal fluctuations in a phenomenological model for the disordered spin density wave state, and find from Monte Carlo simulations that a short anti-ferromagnetic correlation length (order 4-5 lattice constants) already leads to pronounced Fermi arcs in the electronic spectral weight.