Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userHybridization expansion Monte Carlo simulation of multi-orbital quantum impurity problems: matrix product formalism and improved Monte Carlo sampling
168
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We explore two complementary modifications of the hybridization-expansion continuous-time Monte Carlo method, aiming at large multi-orbital quantum impurity problems. One idea is to compute the imaginary-time propagation using a matrix product states representation. We show that bond dimensions considerably smaller than the dimension of the Hilbert space are sufficient to obtain accurate results, and that this approach scales polynomially, rather than exponentially with the number of orbitals. Based on scaling analyses, we conclude that a matrix product state implementation will outperform the exact-diagonalization based method for quantum impurity problems with more than 12 orbitals. The second idea is an improved Monte Carlo sampling scheme which is applicable to all variants of the hybridization expansion method. We show that this so-called sliding window sampling scheme speeds up the simulation by at least an order of magnitude for a broad range of model parameters, with the largest improvements at low temperature.
rate research
Read More
Polaron tunneling is a prominent example of a problem characterized by different energy scales, for which the standard quantum Monte Carlo methods face a slowdown problem. We propose a new quantum-tunneling Monte Carlo (QTMC) method which is free from this issue and can be used for a wide range of tunneling phenomena. We apply it to study an impurity interacting with a one-dimensional Bose-Einstein condensate and simultaneously trapped in an external double-well potential. Our scheme works for an arbitrary coupling between the particle and condensate and, at the same time, allows for an account of tunneling effects. We discover two distinct quasi-particle peaks associated, respectively, with the phonon-assisted tunneling and the self-trapping of the impurity, which are in a crossover regime for the system modeled. We observe and analyze changes in the weights and spectral positions of the peaks (or, equally, effective masses of the quasi-particles) when the coupling strength is increased. Possible experimental realizations using cold atoms are discussed.
We develop an energy density matrix that parallels the one-body reduced density matrix (1RDM) for many-body quantum systems. Just as the density matrix gives access to the number density and occupation numbers, the energy density matrix yields the energy density and orbital occupation energies. The eigenvectors of the matrix provide a natural orbital partitioning of the energy density while the eigenvalues comprise a single particle energy spectrum obeying a total energy sum rule. For mean-field systems the energy density matrix recovers the exact spectrum. When correlation becomes important, the occupation energies resemble quasiparticle energies in some respects. We explore the occupation energy spectrum for the finite 3D homogeneous electron gas in the metallic regime and an isolated oxygen atom with ground state quantum Monte Carlo techniques implemented in the QMCPACK simulation code. The occupation energy spectrum for the homogeneous electron gas can be described by an effective mass below the Fermi level. Above the Fermi level evanescent behavior in the occupation energies is observed in similar fashion to the occupation numbers of the 1RDM. A direct comparison with total energy differences shows a quantitative connection between the occupation energies and electron addition and removal energies for the electron gas. For the oxygen atom, the association between the ground state occupation energies and particle addition and removal energies becomes only qualitative. The energy density matrix provides a new avenue for describing energetics with quantum Monte Carlo methods which have traditionally been limited to total energies.
The recently developed density matrix quantum Monte Carlo (DMQMC) algorithm stochastically samples the N -body thermal density matrix and hence provides access to exact properties of many-particle quantum systems at arbitrary temperatures. We demonstrate that moving to the interaction picture provides substantial benefits when applying DMQMC to interacting fermions. In this first study, we focus on a system of much recent interest: the uniform electron gas in the warm dense regime. The basis set incompleteness error at finite temperature is investigated and extrapolated via a simple Monte Carlo sampling procedure. Finally, we provide benchmark calculations for a four-electron system, comparing our results to previous work where possible.
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.
Frustrated spin systems generically suffer from the negative sign problem inherent to Monte Carlo methods. Since the severity of this problem is formulation dependent, optimization strategies can be put forward. We introduce a phase pinning approach in the realm of the auxiliary field quantum Monte Carlo algorithm. If we can find an anti-unitary operator that commutes with the one body Hamiltonian coupled to the auxiliary field, then the phase of the action is pinned to $0$ and $pi$. For generalized Kitaev models, we can successfully apply this strategy and observe a remarkable improvement of the average sign. We use this method to study thermodynamical and dynamical properties of the Kitaev-Heisenberg model down to temperatures corresponding to half of the exchange coupling constant. Our dynamical data reveals finite temperature properties of ordered and spin-liquid phases inherent to this model.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
