Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDiagrammatic Monte Carlo study of the Frohlich polaron dispersion in 2D and 3D
103
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present results for the solution of the large polaron Frohlich Hamiltonian in 3-dimensions (3D) and 2-dimensions (2D) obtained via the Diagrammatic Monte Carlo (DMC) method. Our implementation is based on the approach by Mishchenko [A.S. Mishchenko et al., Phys. Rev. B 62, 6317 (2000)]. Polaron ground state energies and effective polaron masses are successfully benchmarked with data obtained using Feynmans path integral formalism. By comparing 3D and 2D data, we verify the analytically exact scaling relations for energies and effective masses from 3D$to$2D, which provides a stringent test for the quality of DMC predictions. The accuracy of our results is further proven by providing values for the exactly known coefficients in weak- and strong coupling expansions. Moreover, we compute polaron dispersion curves which are validated with analytically known lower and upper limits in the small coupling regime and verify the first order expansion results for larger couplings, thus disproving previous critiques on the apparent incompatibility of DMC with analytical results and furnishing useful reference for a wide range of coupling strengths.
rate research
Read More
We consider two large polaron systems that are described by a Fr{o}hlich type of Hamiltonian, namely the Bose-Einstein condensate (BEC) polaron in the continuum and the acoustic polaron in a solid. We present ground-state energies of these two systems calculated with the Diagrammatic Monte Carlo (DiagMC) method and with a Feynman all-coupling approach. The DiagMC method evaluates up to very high order a diagrammatic series for the polaron Greens function. The Feynman all-coupling approach is a variational method that has been used for a wide range of polaronic problems. For the acoustic and BEC polaron both methods provide remarkably similar non-renormalized ground-state energies that are obtained after introducing a finite momentum cutoff. For the renormalized ground-state energies of the BEC polaron, there are relatively large discrepancies between the DiagMC and the Feynman predictions. These differences can be attributed to the renormalization procedure for the contact interaction.
We present the first approximation free diagrammatic Monte Carlo study of a lattice polaron interacting with an acoustic phonon branch through the deformation potential. Weak and strong coupling regimes are separated by a self-trapping region where quantum resonance between various possible lattice deformations is seen in the ground state properties, spectral function, and optical conductivity. The unique feature of such polaron is the interplay between long- and short wavelength acoustic vibrations creating a composite phonon cloud and leading to persistent self-trapping due to the existence of multiple quasi-stable states. This results in a spectral response whose structure is much more complex than in any of the previously considered polaron models.
We present a simple trick that allows to consider the sum of all connected Feynman diagrams at fixed position of interaction vertices for general fermionic models. With our approach one achieves superior performance compared to Diagrammatic Monte Carlo, while rendering the algorithmic part dramatically simpler. As we consider the sum of all connected diagrams at once, we allow for cancellations between diagrams with different signs, alleviating the sign problem. Moreover, the complexity of the calculation grows exponentially with the order of the expansion, which should be constrasted with the factorial growth of the standard diagrammatic technique. We illustrate the efficiency of the technique for the two-dimensional Fermi-Hubbard model.
Over the last several years, a new generation of quantum simulations has greatly expanded our understanding of charge density wave phase transitions in Hamiltonians with coupling between local phonon modes and the on-site charge density. A quite different, and interesting, case is one in which the phonons live on the bonds, and hence modulate the electron hopping. This situation, described by the Su-Schrieffer-Heeger (SSH) Hamiltonian, has so far only been studied with quantum Monte Carlo in one dimension. Here we present results for the 2D SSH model, and show that a bond ordered wave (BOW) insulator is present in the ground state at half-filling, and argue that a critical value of the electron-phonon coupling is required for its onset, in contradistinction with the 1D case where BOW exists for any nonzero coupling. We determine the precise nature of the bond ordering pattern, which has hitherto been controversial, and the critical transition temperature, which is associated with a spontaneous breaking of ${cal Z}_4$ symmetry.
Two-dimensional (2D) post-transition metal chalcogenides (PTMC) have attracted attention due to their suitable band gaps and lower exciton binding energies, making them more appropriate for electronic, optical and water-splitting devices than graphene and monolayer transition metal dichalcogenides (TMDs). Of the predicted 2D PTMCs, GaSe has been reliably synthesized and experimentally characterized. Despite this fact, quantities such as lattice parameters and band character vary significantly depending on which density functional theory (DFT) functional is used. Although many-body perturbation theory (GW approximation) has been used to correct the electronic structure and obtain the excited state properties of 2D GaSe, and solving the Bethe-Salpeter equation (BSE) has been used to find the optical gap, we find that the results depend strongly on the starting wavefunction. In attempt to correct these discrepancies, we employed the many-body Diffusion Monte Carlo (DMC) method to calculate the ground and excited state properties of GaSe because DMC has a weaker dependence on the trial wavefunction. We benchmark these results with available experimental data, DFT [local-density approximation, Perdew-Burke-Ernzerhof (PBE), strongly constrained and appropriately normed (SCAN) meta-GGA, and hybrid (HSE06) functionals] and GW-BSE (using PBE and SCAN wavefunctions) results. Our findings confirm monolayer GaSe is an indirect gap semiconductor (Gamma-M) with a quasiparticle electronic gap in close agreement with experiment and low exciton binding energy. We also benchmark the optimal lattice parameter, cohesive energy and ground state charge density with DMC and various DFT methods. We aim to present a terminal theoretical benchmark for pristine monolayer GaSe, which will aide in the further study of 2D PTMCs using DMC methods.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
