No Arabic abstract
We present a method for computing resonant inelastic x-ray scattering (RIXS) spectra in one-dimensional systems using the density matrix renormalization group (DMRG) method. By using DMRG to address the problem, we shift the computational bottleneck from the memory requirements associated with exact diagonalization (ED) calculations to the computational time associated with the DMRG algorithm. This approach is then used to obtain RIXS spectra on cluster sizes well beyond state-of-the-art ED techniques. Using this new procedure, we compute the low-energy magnetic excitations observed in Cu $L$-edge RIXS for the challenging corner shared CuO$_4$ chains, both for large multi-orbital clusters and downfolded $t$-$J$ chains. We are able to directly compare results obtained from both models defined in clusters with identical momentum resolution. In the strong coupling limit, we find that the downfolded $t$-$J$ model captures the features of the magnetic excitations probed by RIXS after a uniform scaling of the spectra is taken into account.
The ladder compound Sr$_{14}$Cu$_{24}$O$_{41}$ is of interest both as a quasi-one-dimensional analog of the superconducting cuprates and as a superconductor in its own right when Sr is substituted by Ca. In order to model resonant inelastic x-ray scattering (RIXS) spectra for this compound, we investigate the simpler SrCu$_{2}$O$_{3}$ system in which the crystal structure contains very similar ladder planes. We approximate the LDA dispersion of SrCu$_{2}$O$_{3}$ by a Cu only two-band tight-binding model. Strong correlation effects are incorporated by assuming an anti-ferromagnetic ground state. The available angle-resolved photoemission (ARPES) and RIXS data on the ladder compound are found to be in reasonable accord with our theoretical predictions.
We analyze the resonant inelastic x-ray scattering (RIXS) spectra at the K edge of Mn in the antiferromagnetic insulating manganite LaMnO3. We make use of the Keldysh-type Green-function formalism, in which the RIXS intensity is described by a product of an incident-photon-dependent factor and a density-density correlation function in the 3d states. We calculate the former factor using the 4p density of states given by an ab initio band structure calculation and the latter using a multi-orbital tight-binding model. The ground state of the model Hamiltonian is evaluated within the Hartree-Fock approximation. Correlation effects are treated within the random phase approximation (RPA). We obtain the RIXS intensity in a wide range of energy-loss 2-15 eV. The spectral shape is strongly modified by the RPA correlation, showing good agreement with the experiments. The incident-photon-energy dependence also agrees well with the experiments. The present mechanism that the RIXS spectra arise from band-to-band transitions to screen the core-hole potential is quite different from the orbiton picture previously proposed, enabling a comprehensive understanding of the RIXS spectra.
The control and detection of crystallographic chirality is an important and challenging scientific problem. Chirality has wide ranging implications from medical physics to cosmology including an intimate but subtle connection in magnetic systems, for example Mn$_{1-x}$Fe$_{x}$Si. X-ray diffraction techniques with resonant or polarized variations of the experimental setup are currently utilized to characterize lattice chirality. We demonstrate using theoretical calculations the feasibility of indirect $K$ -edge bimagnon resonant inelastic X-ray scattering (RIXS) spectrum as a viable experimental technique to distinguish crystallographic handedness. We apply spin wave theory to the recently discovered $sqrt {5}timessqrt {5}$ vacancy ordered chalcogenide Rb$_{0.89}$Fe$_{1.58}$Se$_{2}$ for realistic X-ray experimental set up parameters (incoming energy, polarization, and Bragg angle) to show that the computed RIXS spectrum is sensitive to the underlying handedness (right or left) of the lattice. A Flack parameter definition that incorporates the right- and left- chiral lattice RIXS response is introduced. It is shown that the RIXS response of the multiband magnon system RbFeSe arises both from inter- and intra- band scattering processes. The extinction or survival of these RIXS peaks are sensitive to the underlying chiral lattice orientation. This in turn allows for the identification of the two chiral lattice orientations.
We introduce the transcorrelated Density Matrix Renormalization Group (tcDMRG) theory for the efficient approximation of the energy for strongly correlated systems. tcDMRG encodes the wave function as a product of a fixed Jastrow or Gutzwiller correlator and a matrix product state. The latter is optimized by applying the imaginary-time variant of time-dependent (TD) DMRG to the non-Hermitian transcorrelated Hamiltonian. We demonstrate the efficiency of tcDMRG at the example of the two-dimensional Fermi-Hubbard Hamiltonian, a notoriously difficult target for the DMRG algorithm, for different sizes, occupation numbers, and interaction strengths. We demonstrate fast energy convergence of tcDMRG, which indicates that tcDMRG could increase the efficiency of standard DMRG beyond quasi-monodimensional systems and provides a generally powerful approach toward the dynamic correlation problem of DMRG.
We present an infinite density-matrix renormalization group (DMRG) study of an interacting continuum model of twisted bilayer graphene (tBLG) near the magic angle. Because of the long-range Coulomb interaction and the large number of orbital degrees of freedom, tBLG is difficult to study with standard DMRG techniques -- even constructing and storing the Hamiltonian already poses a major challenge. To overcome these difficulties, we use a recently developed compression procedure to obtain a matrix product operator representation of the interacting tBLG Hamiltonian which we show is both efficient and accurate even when including the spin, valley and orbital degrees of freedom. To benchmark our approach, we focus mainly on the spinless, single-valley version of the problem where, at half-filling, we find that the ground state is a nematic semimetal. Remarkably, we find that the ground state is essentially a k-space Slater determinant, so that Hartree-Fock and DMRG give virtually identical results for this problem. Our results show that the effects of long-range interactions in magic angle graphene can be efficiently simulated with DMRG, and opens up a new route for numerically studying strong correlation physics in spinful, two-valley tBLG, and other moire materials, in future work.