The collective spin and charge excitations of doped cuprates and their relationship to superconductivity are not yet fully understood, particularly in the case of the charge excitations. Here, we study the doping-dependent dynamical spin and charge s
tructure factors of single and multi-orbital models for the one-dimensional corner shared spin-chain cuprates using several numerically exact methods. We find that the singleband Hubbard model can describe the spin and charge excitations of the $pd$-model in the low-energy region, including the particle-hole asymmetry in the spin response. However, our results also reveal that the weight of the interorbital spin excitations between Cu and O orbitals is comparable to the weight of the spin excitations between two Cu orbitals. This finding elucidates the microscopic nature of the spin excitations in the 1D cuprates and sheds light on the spin properties of other oxides. Importantly, we find a particle-hole asymmetry in the orbital-resolved charge excitations, which cannot be described by the singleband Hubbard model and is relevant to resonant inelastic x-ray scattering experiments. Our results imply that the explicit inclusion of the oxygen degrees of freedom may be required to understand experimental observations.
Due to the strong coupling between magnetism and ferroelectricity, $(mathrm{ND}_4)_2mathrm{FeCl}_5cdotmathrm{D}_2mathrm{O}$ exhibits several intriguing magnetic and electric phases. In this letter, we include high-order onsite spin anisotropic intera
ctions in a spin model that successfully captures the ferroelectric phase transitions of $(mathrm{ND}_4)_2mathrm{FeCl}_5cdotmathrm{D}_2mathrm{O}$ under a magnetic field and produces the large weights of high-order harmonic components in the cycloid structure that are observed from neutron diffraction experiments. Moreover, we predict a new ferroelectric phase sandwiched between the FE II and FE III phases in a magnetic field. Our results emphasize the importance of the high-order spin anisotropic interactions and provide a guideline to understand multiferroic materials with rich phase diagrams.
Motivated by a recent inelastic neutron scattering experiment on $mathrm{YbMgGaO}_4$ cite{William2019}, we reinvestigate the homogeneous spin model on the triangular lattice. Using the cluster mean-field theory, we study the phase diagram and the mag
netic-field-induced phase transition. We find that the phase boundary between the stripe state and the $120^{circ}$ antiferromagnetic state is broadened by the magnetic field, leading to a field-induced phase transition. This phase transition is suppressed by the next-nearest neighbor exchange interaction $J_2/J_1$ and vanishes as $J_2/J_1>0.13$. We find a parameter space at $J_2/J_1=0.1$, in which the field-induce transition can be achieved and the deviation of theoretical spin excitation energies from experimental data is only $5.4%$. Our results imply that an effective homogeneous spin model still works in $mathrm{YbMgGaO}_4$.
We study magnetic and charge susceptibilities in the half-filled two-dimensional triangular Hubbard model within the dual fermion approximation in the metallic, Mott insulating, and crossover regions of parameter space. In the textcolor{black}{insula
ting state}, we find strong spin fluctuations at the K point at low energy corresponding to the textcolor{black}{120$^{circ}$} antiferromagnetic order. These spin fluctuations persist into the metallic phase and move to higher energy. We also present data for simulated neutron spectroscopy and textcolor{black}{spin-lattice} relaxation times, and perform direct comparisons to inelastic neutron spectroscopy experiments on the triangular material Ba$_8$CoNb$_6$O$_{24}$ and to the relaxation times on $kappa$-(ET)$_2$Cu$_2$(CN)$_3$. Finally, we present charge susceptibilities in different areas of parameter space, which should correspond to momentum-resolved electron-loss spectroscopy measurements on triangular compounds.
Monte Carlo (MC) simulations are essential computational approaches with widespread use throughout all areas of science. We present a method for accelerating lattice MC simulations using fully connected and convolutional artificial neural networks th
at are trained to perform local and global moves in configuration space, respectively. Both networks take local spacetime MC configurations as input features and can, therefore, be trained using samples generated by conventional MC runs on smaller lattices before being utilized for simulations on larger systems. This new approach is benchmarked for the case of determinant quantum Monte Carlo (DQMC) studies of the two-dimensional Holstein model. We find that both artificial neural networks are capable of learning an unspecified effective model that accurately reproduces the MC configuration weights of the original Hamiltonian and achieve an order of magnitude speedup over the conventional DQMC algorithm. Our approach is broadly applicable to many classical and quantum lattice MC algorithms.
We study a three-orbital Su-Schrieffer-Heeger model defined on a two-dimensional Lieb lattice and in the negative charge transfer regime using determinant quantum Monte Carlo. At half-filling (1 hole/unit cell), we observe a bipolaron insulating phas
e, where the ligand oxygen atoms collapse and expand about alternating cation atoms to produce a bond-disproportionated state. This phase is robust against moderate hole doping but is eventually suppressed at large hole concentrations, leading to a metallic polaron-liquid-like state with fluctuating patches of local distortions. Our results suggest that the polarons are highly disordered in the metallic state and freeze into a periodic array across the metal-to-insulator transition. We also find an $s$-wave superconducting state at finite doping that primarily appears on the oxygen sublattices. Our approach provides an efficient, non-perturbative way to treat bond phonons in higher dimensions and our results have implications for many materials where coupling to bond phonons is the dominant interaction.
We study the interplay between the electron-phonon (e-ph) and on-site electron-electron (e-e) interactions in a three-orbital Hubbard-Holstein model on an extended one-dimensional lattice using determinant quantum Monte Carlo. For weak e-e and e-ph i
nteractions, we observe a competition between an orbital-selective Mott phase (OSMP) and a (multicomponent) charge-density-wave (CDW) insulating phase, with an intermediate metallic phase located between them. For large e-e and e-ph couplings, the OSMP and CDW phases persist, while the metallic phase develops short-range orbital correlations and becomes insulating when both the e-e and e-ph interactions are large but comparable. Many of our conclusions are in line with those drawn from a prior dynamical mean field theory study of the two-orbital Hubbard-Holstein model [Phys. Rev. B 95, 12112(R) (2017)] in infinite dimension, suggesting that the competition between the e-ph and e-e interactions in multiorbital Hubbard-Holstein models leads to rich physics, regardless of the dimension of the system.
We study the interplay between the electron-electron (e-e) and the electron-phonon (e-ph) interactions in the two-orbital Hubbard-Holstein model at half filling using the dynamical mean field theory. We find that the e-ph interaction, even at weak co
uplings, strongly modifies the phase diagram of this model and introduces an orbital-selective Peierls insulating phase (OSPI) that is analogous to the widely studied orbital-selective Mott phase (OSMP). At small e-e and e-ph coupling, we find a competition between the OSMP and the OSPI, while at large couplings, a competition occurs between Mott and charge-density-wave (CDW) insulating phases. We further demonstrate that the Hunds coupling influences the OSPI transition by lowering the energy associated with the CDW. Our results explicitly show that one must be cautious when neglecting the e-ph interaction in multiorbital systems, where multiple electronic interactions create states that are readily influenced by perturbing interactions.
We use determinant quantum Monte Carlo to study the single particle properties of quasiparticles and phonons in a variant of the two-dimensional Holstein model that includes an additional non-linear electron-phonon (e-ph) interaction. We find that a
small positive non-linear interaction reduces the effective coupling between the electrons and the lattice, suppresses charge-density wave (CDW) correlations, and hardens the effective phonon frequency. Conversely, a small negative non-linear interaction can enhance the e-ph coupling resulting in heavier quasiparticles, an increased tendency towards a CDW phase at all fillings, and a softened phonon frequency. An effective linear model with a renormalized interaction strength and phonon frequency can qualitatively capture this physics; however, the quantitative effects of the non-linearity on both the electronic and phononic degrees of freedom cannot be captured by such a model. These results are significant for typical non-linear coupling strengths found in real materials, indicating that non-linearity can have a significant influence on the physics of many e-ph coupled systems.
Determinant quantum Monte Carlo (DQMC) simulations are used to study non-linear electron-phonon interactions in a two-dimensional Holstein-like model on a square lattice. We examine the impact of non-linear electron-lattice interactions on supercondu
ctivity and on Peierls charge-density-wave (CDW) correlations at finite temperatures and carrier concentrations. We find that the CDW correlations are dramatically suppressed with the inclusion of even a small non-linear interaction. Conversely, the effect of the non-linearity on superconductivity is found to be less dramatic at high temperatures; however, we find evidence that the non-linearity is ultimately detrimental to superconductivity. These effects are attributed to the combined hardening of the phonon frequency and a renormalization of the effective linear electron-phonon coupling towards weaker values. These results demonstrate the importance of non-linear interactions at finite carrier concentrations when one is addressing CDW and superconducting order and have implications for experiments that drive the lattice far from equilibrium.
