Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStable Quantum Monte Carlo Simulations for Entanglement Spectra of Interacting Fermions
176
0
0.0
(
0
)
Authors
Fakher F. Assaad
Ask ChatGPT about the research
No Arabic abstract
We show that the two recently proposed methods to compute Renyi entanglement entropies in the realm of determinant quantum Monte Carlo methods for fermions are in principle equivalent, but differ in sampling strategies. The analogy allows to formulate a numerically stable calculation of the entanglement spectrum at strong coupling. We demonstrate the approach by studying static and dynamical properties of the entanglement hamiltonian across the interaction driven quantum phase transition between a topological insulator and quantum antiferromagnet in the Kane-Mele Hubbard model. The formulation is not limited to fermion systems and can readily be adapted to world-line based simulations of bosonic systems.
rate research
Read More
We present a non-iterative solver based on the Schur complement method for sparse linear systems of special form which appear in Quantum Monte-Carlo (QMC) simulations of strongly interacting fermions on the lattice. While the number of floating-point operations for this solver scales as the cube of the number of lattice sites, for practically relevant lattice sizes it is still significantly faster than iterative solvers such as the Conjugate Gradient method in the regime of strong inter-fermion interactions, for example, in the vicinity of quantum phase transitions. The speed-up is even more dramatic for the solution of multiple linear systems with different right-hand sides. We present benchmark results for QMC simulations of the tight-binding models on the hexagonal graphene lattice with on-site (Hubbard) and non-local (Coulomb) interactions, and demonstrate the potential for further speed-up using GPU.
We tutorially review the determinantal Quantum Monte Carlo method for fermionic systems, using the Hubbard model as a case study. Starting with the basic ingredients of Monte Carlo simulations for classical systems, we introduce aspects such as importance sampling, sources of errors, and finite-size scaling analyses. We then set up the preliminary steps to prepare for the simulations, showing that they are actually carried out by sampling discrete Hubbard-Stratonovich auxiliary fields. In this process the Greens function emerges as a fundamental tool, since it is used in the updating process, and, at the same time, it is directly related to the quantities probing magnetic, charge, metallic, and superconducting behaviours. We also discuss the as yet unresolved minus-sign problem, and two ways to stabilize the algorithm at low temperatures.
In the context of realistic calculations for strongly-correlated materials with $d$- or $f$-electrons the efficient computation of multi-orbital models is of paramount importance. Here we introduce a set of invariants for the SU(2)-symmetric Kanamori Hamiltonian which allows to massively speed up the calculation of the fermionic trace in hybridization-expansion continuous-time quantum Monte Carlo algorithms. As an application, we show that, exploiting this set of good quantum numbers, the study of the orbital-selective Mott-transition in systems with up to seven correlated orbitals becomes feasible.
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.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
