Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAlleviating the Sign Problem in Quantum Monte Carlo Simulations of Spin-Orbit-Coupled Multi-Orbital Hubbard Models
88
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a strategy to alleviate the sign problem in continuous-time quantum Monte Carlo (CTQMC) simulations of the dynamical-mean-field-theory (DMFT) equations for the spin-orbit-coupled multiorbital Hubbard model. We first identify the combinations of rotationally invariant Hund coupling terms present in the relativistic basis which lead to a severe sign problem. Exploiting the fact that the average sign in CTQMC depends on the choice of single-particle basis, we propose a bonding-antibonding basis $V_{j3/2mathrm{BA}}$ which shows an improved average sign compared to the widely used relativistic basis for most parameter sets investigated. We then generalize this procedure by introducing a stochastic optimization algorithm that exploits the space of single-particle bases and show that $V_{j3/2mathrm{BA}}$ is very close to optimal within the parameter space investigated. Our findings enable more efficient DMFT simulations of materials with strong spin-orbit coupling.
rate research
Read More
Building on a recent investigation of the Shastry-Sutherland model [S. Wessel et al., Phys. Rev. B 98, 174432 (2018)], we develop a general strategy to eliminate the Monte Carlo sign problem near the zero temperature limit in frustrated quantum spin models. If the Hamiltonian of interest and the sign-problem-free Hamiltonian---obtained by making all off-diagonal elements negative in a given basis---have the same ground state and this state is a member of the computational basis, then the average sign returns to one as the temperature goes to zero. We illustrate this technique by studying the triangular and kagome lattice Heisenberg antiferrromagnet in a magnetic field above saturation, as well as the Heisenberg antiferromagnet on a modified Husimi cactus in the dimer basis. We also provide detailed appendices on using linear programming techniques to automatically generate efficient directed loop updates in quantum Monte Carlo simulations.
We characterize the three-orbital Hubbard model using state-of-the-art determinant quantum Monte Carlo (DQMC) simulations with parameters relevant to the cuprate high-temperature superconductors. The simulations find that doped holes preferentially reside on oxygen orbitals and that the ({pi},{pi}) antiferromagnetic ordering vector dominates in the vicinity of the undoped system, as known from experiments. The orbitally-resolved spectral functions agree well with photoemission spectroscopy studies and enable identification of orbital content in the bands. A comparison of DQMC results with exact diagonalization and cluster perturbation theory studies elucidates how these different numerical techniques complement one another to produce a more complete understanding of the model and the cuprates. Interestingly, our DQMC simulations predict a charge-transfer gap that is significantly smaller than the direct (optical) gap measured in experiment. Most likely, it corresponds to the indirect gap that has recently been suggested to be on the order of 0.8 eV, and demonstrates the subtlety in identifying charge gaps.
Exciton-polaron formation in one-dimensional lattice models with short or long-range carrier-phonon interaction is studied by quantum Monte Carlo simulations. Depending on the relative sign of electron and hole-phonon coupling, the exciton-polaron size increases or decreases with increasing interaction strength. Quantum phonon fluctuations determine the (exciton-)polaron size and yield translation-invariant states at all finite couplings.
We discuss a projector Monte Carlo method for quantum spin models formulated in the valence bond basis, using the S=1/2 Heisenberg antiferromagnet as an example. Its singlet ground state can be projected out of an arbitrary basis state as the trial state, but a more rapid convergence can be obtained using a good variational state. As an alternative to first carrying out a time consuming variational Monte Carlo calculation, we show that a very good trial state can be generated in an iterative fashion in the course of the simulation itself. We also show how the properties of the valence bond basis enable calculations of quantities that are difficult to obtain with the standard basis of Sz eigenstates. In particular, we discuss quantities involving finite-momentum states in the triplet sector, such as the dispersion relation and the spectral weight of the lowest triplet.
We derive the equations for calculating the high-frequency asymptotics of the local two-particle vertex function for a multi-orbital impurity model. These relate the asymptotics for a general local interaction to equal-time two-particle Greens functions, which we sample using continuous-time quantum Monte Carlo simulations with a worm algorithm. As specific examples we study the single-orbital Hubbard model and the three $t_{2g}$ orbitals of SrVO$_3$ within dynamical mean field theory (DMFT). We demonstrate how the knowledge of the high-frequency asymptotics reduces the statistical uncertainties of the vertex and further eliminates finite box size effects. The proposed method benefits the calculation of non-local susceptibilities in DMFT and diagrammatic extensions of DMFT.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
