Quantum Monte Carlo Calculation of the Binding Energy of Bilayer Graphene

We report diffusion quantum Monte Carlo calculations of the interlayer binding energy of bilayer graphene. We find the binding energies of the AA- and AB-stacked structures at the equilibrium separation to be 11.5(9) and 17.7(9) meV/atom, respectively. The out-of-plane zone-center optical phonon frequency predicted by our binding-energy curve is consistent with available experimental results. As well as assisting the modeling of interactions between graphene layers, our results will facilitate the development of van der Waals exchange-correlation functionals for density functional theory calculations.Comment: 5 pages and 3 figures, submitted to Phys. Rev. Lett.; supplemental material is available on arXiv via the ancillary files attached to this submissio

Binding energies of trions and biexcitons in two-dimensional semiconductors from diffusion quantum Monte Carlo calculations

Excitonic effects play a particularly important role in the optoelectronic behavior of two-dimensional (2D) semiconductors. To facilitate the interpretation of experimental photoabsorption and photoluminescence spectra we provide statistically exact diffusion quantum Monte Carlo binding-energy data for Mott-Wannier models of excitons, trions, and biexcitons in 2D semiconductors. We also provide contact pair densities to allow a description of contact (exchange) interactions between charge carriers using first-order perturbation theory. Our data indicate that the binding energy of a trion is generally larger than that of a biexciton in 2D semiconductors. We provide interpolation formulas giving the binding energy and contact density of 2D semiconductors as functions of the electron and hole effective masses and the in-plane polarizability.</p

Diffusion quantum Monte Carlo study of excitonic complexes in two-dimensional transition-metal dichalcogenides

Excitonic effects play a particularly important role in the optoelectronic behavior of two-dimensional semiconductors. To facilitate the interpretation of experimental photoabsorption and photoluminescence spectra we provide (i) statistically exact diffusion quantum Monte Carlo binding-energy data for a Mott-Wannier model of (donor/acceptor-bound) excitons, trions, and biexcitons in two-dimensional semiconductors in which charges interact via the Keldysh potential, (ii) contact pair-distribution functions to allow a perturbative description of contact interactions between charge carriers, and (iii) an analysis and classification of the different types of bright trions and biexcitons that can be seen in single-layer molybdenum and tungsten dichalcogenides. We investigate the stability of biexcitons in which two charge carriers are indistinguishable, finding that they are only bound when the indistinguishable particles are several times heavier than the distinguishable ones. Donor/acceptor-bound biexcitons have similar binding energies to the experimentally measured biexciton binding energies. We predict the relative positions of all stable free and bound excitonic complexes of distinguishable charge carriers in the photoluminescence spectra of WSe2 and MoSe2

Charge-carrier complexes in monolayer semiconductors

The photoluminescence (PL) spectra of monolayer (1L) semiconductors feature peaks ascribed to different charge-carrier complexes. We perform diffusion quantum Monte Carlo simulations of the binding energies of these complexes and examine their response to electric and magnetic fields. We focus on quintons (charged biexcitons), since they are the largest free charge-carrier complexes in transition-metal dichalcogenides (TMDs). We examine the accuracy of the Rytova-Keldysh interaction potential between charges by comparing the binding energies of charge-carrier complexes in 1L-TMDs with results obtained using $\textit{ab initio}$ interaction potentials. Magnetic fields$<8$T change the binding energies (BEs) by$\sim0.2$ meV,T$^{-1}$, in agreement with experiments, with the BE variations of different complexes being very similar. Our results will help identify charge complexes in the PL spectra of 1L-semiconductor

Charge-tuneable biexciton complexes in monolayer WSe2

Monolayer transition metal dichalcogenides have strong Coulomb-mediated many-body interactions. Theoretical studies have predicted the existence of numerous multi-particle excitonic states. Two-particle excitons and three-particle trions have been identified by their optical signatures. However, more complex states such as biexcitons have been elusive due to limited spectral quality of the optical emission. Here, we report direct evidence of two biexciton complexes in monolayer tungsten diselenide: the four-particle neutral biexciton and the five-particle negatively charged biexciton. We distinguish these states by power-dependent photoluminescence and demonstrate full electrical switching between them. We determine the band states of the elementary particles comprising the biexcitons through magneto-optical spectroscopy. We also resolve a splitting of 2.5 meV for the neutral biexciton, which we attribute to the fine structure, providing reference for subsequent studies. Our results unveil the nature of multi-exciton complexes in transitionmetal dichalcogenides and offer direct routes towards deterministic control in many-body quantum phenomena

Noncovalent Interactions by QMC: Speedup by One-Particle Basis-Set Size Reduction

While it is empirically accepted that the fixed-node diffusion Monte-Carlo (FN-DMC) depends only weakly on the size of the one-particle basis sets used to expand its guiding functions, limits of this observation are not settled yet. Our recent work indicates that under the FN error cancellation conditions, augmented triple zeta basis sets are sufficient to achieve a benchmark level of 0.1 kcal/mol in a number of small noncovalent complexes. Here we report on a possibility of truncation of the one-particle basis sets used in FN-DMC guiding functions that has no visible effect on the accuracy of the production FN-DMC energy differences. The proposed scheme leads to no significant increase in the local energy variance, indicating that the total CPU cost of large-scale benchmark noncovalent interaction energy FN-DMC calculations may be reduced.Comment: ACS book chapter, accepte

Density of states of magnetic substitutional impurity-doped graphene in the paramagnetic and ferromagnetic phases

We reveal the effects of magnetic substitutional impurities on the density of states (DOS) of a graphene monolayer through the s-f model as well as coherent potential approximation. We show that the magnetic exchange interaction between the itinerant electrons in graphene and magnetic moments leads to ferromagnetic order and spin-splitting band below Curie temperature, which substantially affect the DOS. Furthermore, ferromagnetic graphene exhibits metallic behavior owing to the strong s-f exchange leading to the appearance of a sharp quasiparticle peak near the Fermi level. These phenomena, along with the gaps appearing in the DOS open prospects for new applications in spintronics and optics. © 2013 Elsevier BY. All rights reserved

Quantum Monte Carlo Calculation of the Binding Energy of Bilayer Graphene

We report diffusion quantum Monte Carlo calculations of the interlayer binding energy of bilayer graphene. We find the binding energies of the AA-and AB-stacked structures at the equilibrium separation to be 11.5(9) and 17.7(9)meV/atom, respectively. The out-of-plane zone-center optical phonon frequency predicted by our binding-energy curve is consistent with available experimental results. As well as assisting the modeling of interactions between graphene layers, our results will facilitate the development of van der Waals exchange-correlation functionals for density functional theory calculations

Quasiparticle and excitonic gaps of one-dimensional carbon chains

We report diffusion quantum Monte Carlo (DMC) calculations of the quasiparticle and excitonic gaps of hydrogen-terminated oligoynes and extended polyyne. The electronic gaps are found to be very sensitive to the atomic structure in these systems. We have therefore optimised the geometry of polyyne by directly minimising the DMC energy with respect to the lattice constant and the Peierls-induced carbon-carbon bond-length alternation. We find the bond-length alternation of polyyne to be 0.136(2) Å and the excitonic and quasiparticle gaps to be 3.30(7) and 3.4(1) eV, respectively. The DMC zone-centre longitudinal optical phonon frequency of polyyne is 2084(5) cm-1, which is consistent with Raman spectroscopic measurements for large oligoynes