Showing 1–50 of 63 results for author: Berkelbach, T C

Periodic Local Coupled-Cluster Theory for Insulators and Metals

Authors: Hong-Zhou Ye, Timothy C. Berkelbach

Abstract: We describe the implementation details of periodic local coupled-cluster theory with single and double excitations (CCSD) and perturbative triple excitations [CCSD(T)] using local natural orbitals (LNOs) and $k$-point symmetry. We discuss and compare several choices for orbital localization, fragmentation, and LNO construction. By studying diamond and lithium, we demonstrate that periodic LNO-CC t… ▽ More We describe the implementation details of periodic local coupled-cluster theory with single and double excitations (CCSD) and perturbative triple excitations [CCSD(T)] using local natural orbitals (LNOs) and $k$-point symmetry. We discuss and compare several choices for orbital localization, fragmentation, and LNO construction. By studying diamond and lithium, we demonstrate that periodic LNO-CC theory can be applied with equal success to both insulators and metals, achieving speedups of two to three orders of magnitude even for moderately sized $k$-point meshes. Our final predictions of the equilibrium cohesive energy, lattice constant, and bulk modulus for diamond and lithium are in good agreement with previous theoretical predictions and experimental results. △ Less

Submitted 15 July, 2024; originally announced July 2024.

Comments: 11 pages, 9 figures, 2 tables

Nonperturbative Simulation of Anharmonic Rattler Dynamics in Type-I Clathrates with Vibrational Dynamical Mean-Field Theory

Authors: Dipti Jasrasaria, Timothy C. Berkelbach

Abstract: We use vibrational dynamical mean-field theory (VDMFT) to study the vibrational structure of type-I clathrate solids, specifically X$_8$Ga$_{16}$Ge$_{30}$, where X=Ba,Sr. These materials are cage-like chemical structures hosting loosely bound guest atoms, resulting in strong anharmonicity, short phonon lifetimes, and ultra-low thermal conductivities. Presenting the methodological developments nece… ▽ More We use vibrational dynamical mean-field theory (VDMFT) to study the vibrational structure of type-I clathrate solids, specifically X$_8$Ga$_{16}$Ge$_{30}$, where X=Ba,Sr. These materials are cage-like chemical structures hosting loosely bound guest atoms, resulting in strong anharmonicity, short phonon lifetimes, and ultra-low thermal conductivities. Presenting the methodological developments necessary for this first application to three-dimensional, atomistic materials, we validate our approach through comparison to molecular dynamics simulations and show that VDMFT is extremely accurate at a fraction of the cost. Through the use of nonperturbative methods, we find that anharmonicity is dominated by four-phonon and higher-order scattering processes, and it causes rattler modes to shift up in frequency by 50% (10 cm$^{-1}$) and to have lifetimes of less than 1 ps; this behavior is not captured by traditional perturbation theory. Furthermore, we analyze the phonon self-energy and find that anharmonicity mixes guest rattling modes and cage acoustic modes, significantly changing the character of the harmonic phonons. △ Less

Submitted 31 July, 2024; v1 submitted 12 February, 2024; originally announced February 2024.

Comments: 15 pages, 12 figures

Theory of Acoustic Polarons in the Two-Dimensional SSH Model Applied to the Layered Superatomic Semiconductor Re6Se8Cl2

Authors: Petra Shih, Timothy C. Berkelbach

Abstract: Layered superatomic semiconductors, whose buildings blocks are atomically precise molecular clusters, exhibit interesting electronic and vibrational properties. In recent work [Science 382, 438 (2023)], transient reflection microscopy revealed quasi-ballistic exciton dynamics in Re6Se8Cl2, which was attributed to the formation of polarons due to coupling with acoustic phonons. Here, we characteriz… ▽ More Layered superatomic semiconductors, whose buildings blocks are atomically precise molecular clusters, exhibit interesting electronic and vibrational properties. In recent work [Science 382, 438 (2023)], transient reflection microscopy revealed quasi-ballistic exciton dynamics in Re6Se8Cl2, which was attributed to the formation of polarons due to coupling with acoustic phonons. Here, we characterize the electronic, excitonic, and phononic properties with periodic density functional theory. We further parameterize a polaron Hamiltonian with nonlocal [Su-Schrieffer-Heeger (SSH)] coupling to acoustic phonon to study the polaron ground state binding energy and dispersion relation with variational wavefunctions. We calculate a polaron binding energy of about 10 meV at room temperature, and the maximum group velocity of our polaron dispersion relation is 1.5 km/s, which is similar to the experimentally observed exciton transport velocity. △ Less

Submitted 25 January, 2024; originally announced January 2024.

Performance of periodic EOM-CCSD for band gaps of inorganic semiconductors and insulators

Authors: Ethan A. Vo, Xiao Wang, Timothy C. Berkelbach

Abstract: We calculate the band gaps of 12 inorganic semiconductors and insulators composed of atoms from the first three rows of the periodic table using periodic equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD). Our calculations are performed with atom-centered triple-zeta basis sets and up to 64 $k$-points in the Brillouin zone. We analyze the convergence behavior w… ▽ More We calculate the band gaps of 12 inorganic semiconductors and insulators composed of atoms from the first three rows of the periodic table using periodic equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD). Our calculations are performed with atom-centered triple-zeta basis sets and up to 64 $k$-points in the Brillouin zone. We analyze the convergence behavior with respect to number of orbitals and number of $k$-points sampled, using composite corrections and extrapolations to produce our final values. When accounting for electron-phonon corrections to experimental band gaps, we find that EOM-CCSD has a mean signed error of $-0.12$ eV and a mean absolute error of $0.42$ eV; the largest outliers are C (error of $-0.93$ eV), BP ($-1.00$ eV), and LiH ($+0.78$ eV). Surprisingly, we find that the more affordable partitioned EOM-MP2 theory performs as well as EOM-CCSD. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 6 pages, 5 figures, 1 table

Adsorption and Vibrational Spectroscopy of CO on the Surface of MgO from Periodic Local Coupled-Cluster Theory

Authors: Hong-Zhou Ye, Timothy C. Berkelbach

Abstract: The adsorption of CO on the surface of MgO has long been a model problem in surface chemistry. Here, we report periodic Gaussian-based calculations for this problem using second-order perturbation theory (MP2) and coupled-cluster theory with single and double excitations (CCSD) and perturbative triple excitations [CCSD(T)], with the latter two performed using a recently developed extension of the… ▽ More The adsorption of CO on the surface of MgO has long been a model problem in surface chemistry. Here, we report periodic Gaussian-based calculations for this problem using second-order perturbation theory (MP2) and coupled-cluster theory with single and double excitations (CCSD) and perturbative triple excitations [CCSD(T)], with the latter two performed using a recently developed extension of the local natural orbital approximation to problems with periodic boundary conditions. The low cost of periodic local correlation calculations allows us to calculate the full CCSD(T) binding curve of CO approaching the surface of MgO (and thus the adsorption energy) and the two-dimensional potential energy surface (PES) as a function of the distance from the surface and the CO stretching coordinate. From the PES, we obtain the fundamental vibrational frequency of CO on MgO, whose shift from the gas phase value is a common experimental probe of surface adsorption. We find that CCSD(T) correctly predicts a positive frequency shift upon adsorption of $+14.7~\textrm{cm}^{-1}$, in excellent agreement with the experimental shift of $+14.3~\textrm{cm}^{-1}$. We use our CCSD(T) results to assess the accuracy of MP2, CCSD, and several density functional theory (DFT) approximations, including exchange correlation functionals and dispersion corrections. We find that MP2 and CCSD yield reasonable binding energies and frequency shifts, whereas many DFT calculations overestimate the magnitude of the adsorption energy by $5$ -- $15$~kJ/mol and predict a negative frequency shift of about $-20~\textrm{cm}^{-1}$, which we attribute to self-interaction-induced delocalization errors that are mildly ameliorated with hybrid functionals. Our findings highlight the accuracy and computational efficiency of the periodic local correlation for the simulation of surface chemistry with accurate wavefunction methods. △ Less

Submitted 27 February, 2024; v1 submitted 26 September, 2023; originally announced September 2023.

Comments: 6 pages, 3 figures, 1 tables (excluding references and SI)

Ab initio surface chemistry with chemical accuracy

Authors: Hong-Zhou Ye, Timothy C. Berkelbach

Abstract: First-principles calculations are a cornerstone of modern surface science and heterogeneous catalysis. However, accurate reaction energies and barrier heights are frequently inaccessible due to the approximations demanded by the large number of atoms. Here we combine developments in local correlation and periodic correlated wavefunction theory to solve the many-electron Schrödinger equation for mo… ▽ More First-principles calculations are a cornerstone of modern surface science and heterogeneous catalysis. However, accurate reaction energies and barrier heights are frequently inaccessible due to the approximations demanded by the large number of atoms. Here we combine developments in local correlation and periodic correlated wavefunction theory to solve the many-electron Schrödinger equation for molecules on surfaces with chemical accuracy, commonly defined as 1~kcal/mol. As a demonstration, we study water on the surface of \ce{Al2O3} and \ce{TiO2}, two prototypical and industrially important metal oxides for which we obtain converged energies at the level of coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)], commonly known as the "gold-standard" in molecular quantum chemistry. We definitively resolve the energetics associated with water adsorption and dissociation, enabling us to address recent experiments and to analyze the errors of more commonly used approximate theories. △ Less

Submitted 7 February, 2024; v1 submitted 25 September, 2023; originally announced September 2023.

Comments: 6 pages (w/o SI), 3 figures

Toward linear scaling auxiliary field quantum Monte Carlo with local natural orbitals

Authors: Jo S. Kurian, Hong-Zhou Ye, Ankit Mahajan, Timothy C. Berkelbach, Sandeep Sharma

Abstract: We develop a local correlation variant of auxiliary field quantum Monte Carlo (AFQMC) that is based on local natural orbitals (LNO-AFQMC). In LNO-AFQMC, independent AFQMC calculations are performed for each localized occupied orbital using a truncated set of tailored orbitals. Because the size of this space does not grow with system size for a target accuracy, the method has linear scaling. Applyi… ▽ More We develop a local correlation variant of auxiliary field quantum Monte Carlo (AFQMC) that is based on local natural orbitals (LNO-AFQMC). In LNO-AFQMC, independent AFQMC calculations are performed for each localized occupied orbital using a truncated set of tailored orbitals. Because the size of this space does not grow with system size for a target accuracy, the method has linear scaling. Applying LNO AFQMC to molecular problems containing a few hundred to a thousand orbitals, we demonstrate convergence of total energies with significantly reduced costs. The savings are more significant for larger systems and larger basis sets. However, even for our smallest system studied, we find that LNO-AFQMC is cheaper than canonical AFQMC, in contrast with many other reduced-scaling methods. Perhaps most significantly, we show that energy differences converge much more quickly than total energies, making the method ideal for applications in chemistry and material science. Our work paves the way for linear scaling AFQMC calculations of strongly correlated systems, which would have a transformative effect on ab initio quantum chemistry. △ Less

Submitted 23 August, 2023; originally announced August 2023.

Comments: 9 pages, 2 figures

Can spin-component scaled MP2 achieve kJ/mol accuracy for cohesive energies of molecular crystals?

Authors: Yu Hsuan Liang, Hong-Zhou Ye, Timothy C. Berkelbach

Abstract: Achieving kJ/mol accuracy in the cohesive energy of molecular crystals, as necessary for crystal structure prediction and the resolution of polymorphism, is an ongoing challenge in computational materials science. Here, we evaluate the performance of second-order Møller-Plesset perturbation theory (MP2), including its spin-component scaled models, by calculating the cohesive energies of the 23 mol… ▽ More Achieving kJ/mol accuracy in the cohesive energy of molecular crystals, as necessary for crystal structure prediction and the resolution of polymorphism, is an ongoing challenge in computational materials science. Here, we evaluate the performance of second-order Møller-Plesset perturbation theory (MP2), including its spin-component scaled models, by calculating the cohesive energies of the 23 molecular crystals contained in the X23 dataset. Our calculations are performed with periodic boundary conditions and Brillouin zone sampling, and we converge results to the thermodynamic limit and the complete basis set limit to an accuracy of about 1 kJ/mol (0.25 kcal/mol), which is rarely achieved in previous MP2 calculations of molecular crystals. Comparing to experimental cohesive energies, we find that MP2 has a mean absolute error of 12.9 kJ/mol, which is comparable to that of DFT using the PBE functional and TS dispersion correction. Separate scaling of the opposite-spin and same-spin components of the correlation energy, with parameters previously determined for molecular interactions, reduces the mean absolute error to 9.5 kJ/mol, and reoptimizing the spin-component scaling parameters for the X23 set further reduces the mean absolute error to 7.5 kJ/mol. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Interplay between local moment and itinerant magnetism in the layered metallic antiferromagnet TaFe$_{1.14}$Te$_3$

Authors: Sae Young Han, Evan J. Telford, Asish K. Kundu, Sylvia J. Bintrim, Simon Turkel, Ren A. Wiscons, Amirali Zangiabadi, Eun-Sang Choi, Tai-De Li, Michael L. Steigerwald, Timothy C. Berkelbach, Abhay N. Pasupathy, Cory R. Dean, Colin Nuckolls, Xavier Roy

Abstract: Two-dimensional (2D) antiferromagnets have garnered considerable interest for the next generation of functional spintronics. However, many available bulk materials from which 2D antiferromagnets are isolated are limited by their sensitivity to air, low ordering temperatures, and insulating transport properties. TaFe$_{1+y}$Te$_3$ offers unique opportunities to address these challenges with increas… ▽ More Two-dimensional (2D) antiferromagnets have garnered considerable interest for the next generation of functional spintronics. However, many available bulk materials from which 2D antiferromagnets are isolated are limited by their sensitivity to air, low ordering temperatures, and insulating transport properties. TaFe$_{1+y}$Te$_3$ offers unique opportunities to address these challenges with increased air stability, metallic transport properties, and robust antiferromagnetic order. Here, we synthesize TaFe$_{1+y}$Te$_3$ ($y$ = 0.14), identify its structural, magnetic, and electronic properties, and elucidate the relationships between them. Axial-dependent high-field magnetization measurements on TaFe$_{1.14}$Te$_3$ reveal saturation magnetic fields ranging between 27-30 T with a saturation magnetic moment of 2.05-2.12 $μ_B$. Magnetotransport measurements confirm TaFe$_{1.14}$Te$_3$ is metallic with strong coupling between magnetic order and electronic transport. Angle-resolved photoemission spectroscopy measurements across the magnetic transition uncover a complex interplay between itinerant electrons and local magnetic moments that drives the magnetic transition. We further demonstrate the ability to isolate few-layer sheets of TaFe$_{1.14}$Te$_3$ through mechanical exfoliation, establishing TaFe$_{1.14}$Te$_3$ as a potential platform for 2D spintronics based on metallic layered antiferromagnets. △ Less

Submitted 3 July, 2023; originally announced July 2023.

Comments: 30 pages, 5 main figures, 23 supporting figures, and 3 supporting tables

Ab initio quantum many-body description of superconducting trends in the cuprates

Authors: Zhi-Hao Cui, Junjie Yang, Johannes Tölle, Hong-Zhou Ye, Huanchen Zhai, Raehyun Kim, Xing Zhang, Lin Lin, Timothy C. Berkelbach, Garnet Kin-Lic Chan

Abstract: Using a systematic ab initio quantum many-body approach that goes beyond low-energy models, we directly compute the superconducting pairing order of several doped cuprate materials and structures. We find that we can correctly capture two well-known trends: the pressure effect, where pairing order increases with intra-layer pressure, and the layer effect, where the pairing order varies with the nu… ▽ More Using a systematic ab initio quantum many-body approach that goes beyond low-energy models, we directly compute the superconducting pairing order of several doped cuprate materials and structures. We find that we can correctly capture two well-known trends: the pressure effect, where pairing order increases with intra-layer pressure, and the layer effect, where the pairing order varies with the number of copper-oxygen layers. From these calculations, we observe that the strength of superexchange and the covalency at optimal doping are the best descriptors of the maximal pairing order. Our microscopic analysis further identifies short-range copper spin fluctuations, together with multi-orbital charge fluctuations, as central to the pairing trends. Our work illustrates the possibility of a quantitative computational understanding of unconventional high-temperature superconducting materials. △ Less

Submitted 12 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 10 pages, 5 figures, with supplementary materials

Room temperature wavelike exciton transport in a van der Waals superatomic semiconductor

Authors: Jakhangirkhodja A. Tulyagankhodjaev, Petra Shih, Jessica Yu, Jake C. Russell, Daniel G. Chica, Michelle E. Reynoso, Haowen Su, Athena C. Stenor, Xavier Roy, Timothy C. Berkelbach, Milan Delor

Abstract: The transport of energy and information in semiconductors is limited by scattering between electronic carriers and lattice phonons, resulting in diffusive and lossy transport that curtails all semiconductor technologies. Using Re6Se8Cl2, a van der Waals (vdW) superatomic semiconductor, we demonstrate the formation of acoustic exciton-polarons, an electronic quasiparticle shielded from phonon scatt… ▽ More The transport of energy and information in semiconductors is limited by scattering between electronic carriers and lattice phonons, resulting in diffusive and lossy transport that curtails all semiconductor technologies. Using Re6Se8Cl2, a van der Waals (vdW) superatomic semiconductor, we demonstrate the formation of acoustic exciton-polarons, an electronic quasiparticle shielded from phonon scattering. We directly image polaron transport in Re6Se8Cl2 at room temperature and reveal quasi-ballistic, wavelike propagation sustained for nanoseconds and several microns. Shielded polaron transport leads to electronic energy propagation orders of magnitude greater than in other vdW semiconductors, exceeding even silicon over nanoseconds. We propose that, counterintuitively, quasi-flat electronic bands and strong exciton-acoustic phonon coupling are together responsible for the remarkable transport properties of Re6Se8Cl2, establishing a new path to ballistic room-temperature semiconductors. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Journal ref: Science, 2023, vol. 382, pp. 438-442

Highly accurate electronic structure of metallic solids from coupled-cluster theory with nonperturbative triple excitations

Authors: Verena A. Neufeld, Timothy C. Berkelbach

Abstract: Coupled-cluster theory with single, double, and perturbative triple excitations (CCSD(T)) -- often considered the "gold standard" of main-group quantum chemistry -- is inapplicable to three-dimensional metals due to an infrared divergence, preventing its application to many important problems in materials science. We study the full, nonperturbative inclusion of triple excitations (CCSDT) and propo… ▽ More Coupled-cluster theory with single, double, and perturbative triple excitations (CCSD(T)) -- often considered the "gold standard" of main-group quantum chemistry -- is inapplicable to three-dimensional metals due to an infrared divergence, preventing its application to many important problems in materials science. We study the full, nonperturbative inclusion of triple excitations (CCSDT) and propose a new, iterative method, which we call ring-CCSDT, that resums the essential triple excitations with the same $N^7$ run-time scaling as CCSD(T). CCSDT and ring-CCSDT are used to calculate the correlation energy of the uniform electron gas at metallic densities and the structural properties of solid lithium. Inclusion of connected triple excitations is shown to be essential to achieving high accuracy. We also investigate semiempirical CC methods based on spin-component scaling and the distinguishable cluster approximation and find that they enhance the accuracy of their parent ab initio methods. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Optical properties of defects in solids via quantum embedding with good active space orbitals

Authors: Bryan T. G. Lau, Brian Busemeyer, Timothy C. Berkelbach

Abstract: The study of isolated defects in solids is a natural target for classical or quantum embedding methods that treat the defect at a high level of theory and the rest of the solid at a lower level of theory. Here, in the context of active-space-based quantum embeddings, we study the performance of three active-space orbital selection schemes based on canonical (energy-ordered) orbitals, local orbital… ▽ More The study of isolated defects in solids is a natural target for classical or quantum embedding methods that treat the defect at a high level of theory and the rest of the solid at a lower level of theory. Here, in the context of active-space-based quantum embeddings, we study the performance of three active-space orbital selection schemes based on canonical (energy-ordered) orbitals, local orbitals defined in the spirit of density matrix embedding theory, and approximate natural transition orbitals. Using equation-of-motion coupled-cluster theory with single and double excitations (CCSD), we apply these active space selection schemes to the calculation of the vertical singlet excitation energy of a substitutional carbon dimer defect in hexagonal boron nitride, an oxygen vacancy in magnesium oxide, and a carbon vacancy in diamond. Especially when used in combination with a simple composite correction, we find that the best performing schemes can predict the excitation energy to about 0.1-0.2 eV of its converged value using only a few hundred orbitals, even when the full supercell has thousands of orbitals, which amounts to many-orders-of-magnitude computational savings when using correlated electronic structure theories. When compared to assigned experimental spectra and accounting for vibrational corrections, we find that CCSD predicts excitation energies that are accurate to about 0.1-0.3 eV, which is comparable to its performance in molecules and bulk solids. △ Less

Submitted 23 January, 2023; originally announced January 2023.

Comments: 8 pages, 4 figures

Machine learning potentials from transfer learning of periodic correlated electronic structure methods: Application to liquid water with AFQMC, CCSD, and CCSD(T)

Authors: Michael S. Chen, Joonho Lee, Hong-Zhou Ye, Timothy C. Berkelbach, David R. Reichman, Thomas E. Markland

Abstract: Obtaining the atomistic structure and dynamics of disordered condensed phase systems from first principles remains one of the forefront challenges of chemical theory. Here we exploit recent advances in periodic electronic structure to show that, by leveraging transfer learning starting from lower tier electronic structure methods, one can obtain machine learned potential energy surfaces for liquid… ▽ More Obtaining the atomistic structure and dynamics of disordered condensed phase systems from first principles remains one of the forefront challenges of chemical theory. Here we exploit recent advances in periodic electronic structure to show that, by leveraging transfer learning starting from lower tier electronic structure methods, one can obtain machine learned potential energy surfaces for liquid water from the higher tier AFQMC, CCSD, and CCSD(T) approaches using $\le$200 energies. By performing both classical and path integral molecular dynamics simulations on these machine learned potential energy surfaces we uncover the interplay of dynamical electron correlation and nuclear quantum effects across the entire liquid range of water while providing a general strategy for efficiently utilizing periodic correlated electronic structure methods to explore disordered condensed phase systems. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Excitons and their Fine Structure in Lead Halide Perovskite Nanocrystals from Atomistic GW/BSE Calculations

Authors: Giulia Biffi, Yeongsu Cho, Roman Krahne, Timothy C. Berkelbach

Abstract: Atomistically detailed computational studies of nanocrystals, such as those derived from the promising lead-halide perovskites, are challenging due to the large number of atoms and lack of symmetries to exploit. Here, focusing on methylammonium lead iodide nanocrystals, we combine a real-space tight binding model with the GW approximation to the self-energy and obtain exciton wavefunctions and abs… ▽ More Atomistically detailed computational studies of nanocrystals, such as those derived from the promising lead-halide perovskites, are challenging due to the large number of atoms and lack of symmetries to exploit. Here, focusing on methylammonium lead iodide nanocrystals, we combine a real-space tight binding model with the GW approximation to the self-energy and obtain exciton wavefunctions and absorption spectra via solutions of the associated Bethe-Salpeter equation. We find that the size dependence of carrier confinement, dielectric contrast, electron-hole exchange, and exciton binding energies has a strong impact on the lowest excitation energy, which can be tuned by almost 1 eV over the diameter range of 2-6 nm. Our calculated excitation energies are about 0.2 eV higher than experimentally measured photoluminescence, and they display the same qualitative size dependence. Focusing on the fine structure of the band-edge excitons, we find that the lowest-lying exciton is spectroscopically dark and about 20-30 meV lower in energy than the higher-lying triplet of bright states, whose degeneracy is slightly broken by crystal field effects. △ Less

Submitted 3 October, 2022; originally announced October 2022.

Comments: 8 pages, 4 figures

Accurate thermochemistry of covalent and ionic solids from spin-component-scaled MP2

Authors: Tamar Goldzak, Xiao Wang, Hong-Zhou Ye, Timothy C. Berkelbach

Abstract: We study the performance of spin-component-scaled second-order Møller-Plesset perturbation theory (SCS-MP2) for the prediction of the lattice constant, bulk modulus, and cohesive energy of 12 simple, three-dimensional, covalent and ionic semiconductors and insulators. We find that SCS-MP2 and the simpler scaled opposite-spin MP2 (SOS-MP2) yield predictions that are significantly improved over the… ▽ More We study the performance of spin-component-scaled second-order Møller-Plesset perturbation theory (SCS-MP2) for the prediction of the lattice constant, bulk modulus, and cohesive energy of 12 simple, three-dimensional, covalent and ionic semiconductors and insulators. We find that SCS-MP2 and the simpler scaled opposite-spin MP2 (SOS-MP2) yield predictions that are significantly improved over the already good performance of MP2. Specifically, when compared to experimental values with zero-point vibrational corrections, SCS-MP2 (SOS-MP2) yields mean absolute errors of 0.015 (0.017) Å for the lattice constant, 3.8 (3.7) GPa for the bulk modulus, and 0.06 (0.08) eV for the cohesive energy, which are smaller than those of leading density functionals by about a factor of two or more. We consider a reparameterization of the spin scaling parameters and find that the optimal parameters for these solids are very similar to those already in common use in molecular quantum chemistry, suggesting good transferability and reliable future applications to surface chemistry on insulators. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: 7 pages, 4 figures

Integral-direct Hartree-Fock and Møller-Plesset Perturbation Theory for Periodic Systems with Density Fitting: Application to the Benzene Crystal

Authors: Sylvia J. Bintrim, Timothy C. Berkelbach, Hong-Zhou Ye

Abstract: We present an algorithm and implementation of integral-direct, density-fitted Hartree-Fock (HF) and second-order Møller-Plesset perturbation theory (MP2) for periodic systems. The new code eliminates the formerly prohibitive storage requirements and allows us to study systems one order of magnitude larger than before at the periodic MP2 level. We demonstrate the significance of the development by… ▽ More We present an algorithm and implementation of integral-direct, density-fitted Hartree-Fock (HF) and second-order Møller-Plesset perturbation theory (MP2) for periodic systems. The new code eliminates the formerly prohibitive storage requirements and allows us to study systems one order of magnitude larger than before at the periodic MP2 level. We demonstrate the significance of the development by studying the benzene crystal in both the thermodynamic limit and the complete basis set limit, for which we predict an MP2 cohesive energy of $-72.8$ kJ/mol, which is about $10$--$15$ kJ/mol larger in magnitude than all previously reported MP2 calculations. Compared to the best theoretical estimate from literature, several modified MP2 models approach chemical accuracy in the predicted cohesive energy of the benzene crystal and hence may be promising cost-effective choices for future applications on molecular crystals. △ Less

Submitted 2 August, 2022; v1 submitted 3 June, 2022; originally announced June 2022.

Comments: 5 pages (w/o refs or SI), 3 figures, 1 table

Journal ref: J. Chem. Theory Comput. 2022, 18, 9, 5374--5381

Conductivity of an electron coupled to anharmonic phonons

Authors: Jonathan H. Fetherolf, Petra Shih, Timothy C. Berkelbach

Abstract: We study the impact of phonon anharmonicity on the electronic dynamics of soft materials using a nonperturbative quantum-classical approach. The method is applied to a one-dimensional model of doped organic semiconductors with low-frequency intermolecular lattice phonons. We find that anharmonicity that leads to phonon hardening increases the mobility and anharmonicity that leads to phonon softeni… ▽ More We study the impact of phonon anharmonicity on the electronic dynamics of soft materials using a nonperturbative quantum-classical approach. The method is applied to a one-dimensional model of doped organic semiconductors with low-frequency intermolecular lattice phonons. We find that anharmonicity that leads to phonon hardening increases the mobility and anharmonicity that leads to phonon softening decreases the mobility. We also test various approximations, including the use of adiabatic phonon disorder, an effective harmonic model with temperature-dependent frequencies, and the Boltzmann transport equation with second-order perturbation theory scattering rates. Overall, we find surprisingly good agreement between all methods but that accounting for phonon anharmonicity is important for accurate prediction of electronic transport including both quantitative mobility values and their qualitative temperature dependence. For the model studied, phonon lifetime effects have relatively little impact on carrier transport, but the effective frequency shift due to anharmonicity is essential. In cases with highly asymmetric, non-Gaussian disorder, an effective harmonic model cannot quantitatively reproduce mobilities or finite-frequency conductivity, and this is especially true for acoustic phonons. △ Less

Submitted 19 May, 2022; originally announced May 2022.

The normal state of attractive Fermi gases from coupled-cluster theory

Authors: James M. Callahan, John Sous, Timothy C. Berkelbach

Abstract: We introduce coupled-cluster (CC) theory for the numerical study of the normal state of two-component, dilute Fermi gases with attractive, short-range interactions at zero temperature. We focus on CC theory with double excitations (CCD) and discuss its close relationship with -- and improvement upon -- the t-matrix approximation, i.e., the resummation of ladder diagrams via a random-phase approxim… ▽ More We introduce coupled-cluster (CC) theory for the numerical study of the normal state of two-component, dilute Fermi gases with attractive, short-range interactions at zero temperature. We focus on CC theory with double excitations (CCD) and discuss its close relationship with -- and improvement upon -- the t-matrix approximation, i.e., the resummation of ladder diagrams via a random-phase approximation. We further discuss its relationship with Chevy's variational wavefunction ansatz for the Fermi polaron and argue that CCD is its natural extension to nonzero minority species concentrations. Studying normal state energetics for a range of interaction strengths below and above unitarity, we find that CCD yields good agreement with fixed-node diffusion Monte Carlo. We find that CCD does not converge for small polarizations and large interaction strengths, which we speculatively attribute to the nascent instability to a superfluid state. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: 7 pages, 3 figures

Journal ref: Phys. Rev. A 106, 033303 (2022)

Ground-state properties of metallic solids from ab initio coupled-cluster theory

Authors: Verena A. Neufeld, Hong-Zhou Ye, Timothy C. Berkelbach

Abstract: Metallic solids are a challenging target for wavefunction-based electronic structure theories and have not been studied in great detail by such methods. Here, we use coupled-cluster theory with single and double excitations (CCSD) to study the structure of solid lithium and aluminum using optimized Gaussian basis sets. We calculate the equilibrium lattice constant, bulk modulus, and cohesive energ… ▽ More Metallic solids are a challenging target for wavefunction-based electronic structure theories and have not been studied in great detail by such methods. Here, we use coupled-cluster theory with single and double excitations (CCSD) to study the structure of solid lithium and aluminum using optimized Gaussian basis sets. We calculate the equilibrium lattice constant, bulk modulus, and cohesive energy and compare them to experimental values, finding accuracy comparable to common density functionals. Because the quantum chemical "gold standard" CCSD(T) (CCSD with perturbative triple excitations) is inapplicable to metals in the thermodynamic limit, we test two approximate improvements to CCSD, which are found to improve the predicted cohesive energies. △ Less

Submitted 4 April, 2022; originally announced April 2022.

Dark-exciton driven energy funneling into dielectric inhomogeneities in two-dimensional semiconductors

Authors: Haowen Su, Ding Xu, Shan-Wen Cheng, Baichang Li, Song Liu, Kenji Watanabe, Takashi Taniguchi, Timothy C. Berkelbach, James Hone, Milan Delor

Abstract: The optoelectronic and transport properties of two-dimensional transition metal dichalcogenide semiconductors (2D TMDs) are highly susceptible to external perturbation, enabling precise tailoring of material function through post-synthetic modifications. Here we show that nanoscale inhomogeneities known as nanobubbles can be used for both strain and, less invasively, dielectric tuning of exciton t… ▽ More The optoelectronic and transport properties of two-dimensional transition metal dichalcogenide semiconductors (2D TMDs) are highly susceptible to external perturbation, enabling precise tailoring of material function through post-synthetic modifications. Here we show that nanoscale inhomogeneities known as nanobubbles can be used for both strain and, less invasively, dielectric tuning of exciton transport in bilayer tungsten disulfide (WSe2). We use ultrasensitive spatiotemporally resolved optical scattering microscopy to directly image exciton transport, revealing that dielectric nanobubbles are surprisingly efficient at funneling and trapping excitons at room temperature, even though the energies of the bright excitons are negligibly affected. Our observations suggest that exciton funneling in dielectric inhomogeneities is driven by momentum-indirect (dark) excitons whose energies are more sensitive to dielectric perturbations than bright excitons. These results reveal a new pathway to control exciton transport in 2D semiconductors with exceptional spatial and energetic precision using dielectric engineering of dark state energetic landscapes. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Full Configuration Interaction Excited-State Energies in Large Active Spaces from Subspace Iteration with Repeated Random Sparsification

Authors: Samuel M. Greene, Robert J. Webber, James E. T. Smith, Jonathan Weare, Timothy C. Berkelbach

Abstract: We present a stable and systematically improvable quantum Monte Carlo (QMC) approach to calculating excited-state energies, which we implement using our fast randomized iteration method for the full configuration interaction problem (FCI-FRI). Unlike previous excited-state quantum Monte Carlo methods, our approach, which is an asymmetric variant of subspace iteration, avoids the use of dot product… ▽ More We present a stable and systematically improvable quantum Monte Carlo (QMC) approach to calculating excited-state energies, which we implement using our fast randomized iteration method for the full configuration interaction problem (FCI-FRI). Unlike previous excited-state quantum Monte Carlo methods, our approach, which is an asymmetric variant of subspace iteration, avoids the use of dot products of random vectors and instead relies upon trial vectors to maintain orthogonality and estimate eigenvalues. By leveraging recent advances, we apply our method to calculate ground- and excited-state energies of strongly correlated molecular systems in large active spaces, including the carbon dimer with 8 electrons in 108 orbitals (8e,108o), an oxo-Mn(salen) transition metal complex (28e,28o), ozone (18e,87o), and butadiene (22e,82o). In the majority of these test cases, our approach yields total excited-state energies that agree with those from state-of-the-art methods -- including heat-bath CI, the density matrix renormalization group approach, and FCIQMC -- to within sub-milliHartree accuracy. In all cases, estimated excitation energies agree to within about 0.1 eV. △ Less

Submitted 12 October, 2022; v1 submitted 28 January, 2022; originally announced January 2022.

Comments: 16 pages, 5 figures, 3 tables

Full-frequency dynamical Bethe-Salpeter equation without frequency and a study of double excitations

Authors: Sylvia J. Bintrim, Timothy C. Berkelbach

Abstract: The Bethe-Salpeter equation (BSE) that results from the GW approximation to the self-energy is a frequency-dependent (nonlinear) eigenvalue problem due to the dynamically screened Coulomb interaction between electrons and holes. The computational time required for a numerically exact treatment of this frequency dependence is $O(N^6)$, where $N$ is the system size. To avoid the common static screen… ▽ More The Bethe-Salpeter equation (BSE) that results from the GW approximation to the self-energy is a frequency-dependent (nonlinear) eigenvalue problem due to the dynamically screened Coulomb interaction between electrons and holes. The computational time required for a numerically exact treatment of this frequency dependence is $O(N^6)$, where $N$ is the system size. To avoid the common static screening approximation, we show that the full-frequency dynamical BSE can be exactly reformulated as a frequency-independent eigenvalue problem in an expanded space of single and double excitations. When combined with an iterative eigensolver and the density fitting approximation to the electron repulsion integrals, this reformulation yields a dynamical BSE algorithm whose computational time is $O(N^5)$, which we verify numerically. Furthermore, the reformulation provides direct access to excited states with dominant double excitation character, which are completely absent in the spectrum of the statically screened BSE. We study the $2^1A_\mathrm{g}$ state of butadiene, hexatriene, and octatetraene and find that GW/BSE overestimates the excitation energy by about 1.5-2 eV and significantly underestimates the double excitation character. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: 5 pages, 2 figures, 1 table

A simplified GW/BSE approach for charged and neutral excitation energies of large molecules and nanomaterials

Authors: Yeongsu Cho, Sylvia J. Bintrim, Timothy C. Berkelbach

Abstract: Inspired by Grimme's simplified Tamm-Dancoff density functional theory approach [S. Grimme, J. Chem. Phys. \textbf{138}, 244104 (2013)], we describe a simplified approach to excited state calculations within the GW approximation to the self-energy and the Bethe-Salpeter equation (BSE), which we call sGW/sBSE. The primary simplification to the electron repulsion integrals yields the same structure… ▽ More Inspired by Grimme's simplified Tamm-Dancoff density functional theory approach [S. Grimme, J. Chem. Phys. \textbf{138}, 244104 (2013)], we describe a simplified approach to excited state calculations within the GW approximation to the self-energy and the Bethe-Salpeter equation (BSE), which we call sGW/sBSE. The primary simplification to the electron repulsion integrals yields the same structure as with tensor hypercontraction, such that our method has a storage requirement that grows quadratically with system size and computational timing that grows cubically with system size. The performance of sGW is tested on the ionization potential of the molecules in the GW100 test set, for which it differs from \textit{ab intio} GW calculations by only 0.2 eV. The performance of sBSE (based on sGW input) is tested on the excitation energies of molecules in the Thiel set, for which it differs from \textit{ab intio} GW/BSE calculations by about 0.5 eV. As examples of the systems that can be routinely studied with sGW/sBSE, we calculate the band gap and excitation energy of hydrogen-passivated silicon nanocrystals with up to 2650 electrons in 4678 spatial orbitals and the absorption spectra of two large organic dye molecules with hundreds of atoms. △ Less

Submitted 9 September, 2021; originally announced September 2021.

Comments: 7 pages, 4 figures, 2 tables

Anharmonic Lattice Dynamics from Vibrational Dynamical Mean-Field Theory

Authors: Petra Shih, Timothy C. Berkelbach

Abstract: We present a vibrational dynamical mean-field theory (VDMFT) of the dynamics of atoms in solids with anharmonic interactions. Like other flavors of DMFT, VDMFT maps the dynamics of a periodic anharmonic lattice of atoms onto those of a self-consistently defined impurity problem with local anharmonicity and coupling to a bath of harmonic oscillators. VDMFT is exact in the harmonic and molecular lim… ▽ More We present a vibrational dynamical mean-field theory (VDMFT) of the dynamics of atoms in solids with anharmonic interactions. Like other flavors of DMFT, VDMFT maps the dynamics of a periodic anharmonic lattice of atoms onto those of a self-consistently defined impurity problem with local anharmonicity and coupling to a bath of harmonic oscillators. VDMFT is exact in the harmonic and molecular limits, nonperturbative, systematically improvable through its cluster extensions, usable with classical or quantum impurity solvers (depending on the importance of nuclear quantum effects), and can be combined with existing low-level diagrammatic theories of anharmonicity. When tested on models of anharmonic optical and acoustic phonons, we find that classical VDMFT gives good agreement with classical molecular dynamics, including the temperature dependence of phonon frequencies and lifetimes. Using a quantum impurity solver, signatures of nuclear quantum effects are observed at low temperatures. We test the description of nonlocal anharmonicity via cellular VDMFT and the combination with self-consistent phonon (SCPH) theory, yielding the powerful SCPH+VDMFT approach. △ Less

Submitted 9 August, 2022; v1 submitted 31 August, 2021; originally announced September 2021.

Comments: 8 pages, 4 figures

Absorption Spectra of Solids from Periodic Equation-of-Motion Coupled-Cluster Theory

Authors: Xiao Wang, Timothy C. Berkelbach

Abstract: We present ab initio absorption spectra of six three-dimensional semiconductors and insulators calculated using Gaussian-based periodic equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD). The spectra are calculated efficiently by solving a system of linear equations at each frequency, giving access to an energy range of tens of eV without explicit enumeration o… ▽ More We present ab initio absorption spectra of six three-dimensional semiconductors and insulators calculated using Gaussian-based periodic equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD). The spectra are calculated efficiently by solving a system of linear equations at each frequency, giving access to an energy range of tens of eV without explicit enumeration of excited states. We assess the impact of Brillouin zone sampling, for which it is hard to achieve convergence due to the cost of EOM-CCSD. Although our most converged spectra exhibit lineshapes that are in good agreement with experiment, they are uniformly shifted to higher energies by about 1 eV. We tentatively attribute this discrepancy to a combination of vibrational effects and the remaining electron correlation, i.e., triple excitations and above. △ Less

Submitted 9 July, 2021; originally announced July 2021.

Comments: 7 pages, 4 figures

Improving MP2 band gaps with low-scaling approximations to EOM-CCSD

Authors: Malte F. Lange, Timothy C. Berkelbach

Abstract: Despite its reasonable accuracy for ground-state properties of semiconductors and insulators, second-order Moller-Plesset perturbation theory (MP2) significantly underestimates band gaps. Here, we evaluate the band gap predictions of partitioned equation-of-motion MP2 (P-EOM-MP2), which is a second-order approximation to equation-of-motion coupled-cluster theory with single and double excitations.… ▽ More Despite its reasonable accuracy for ground-state properties of semiconductors and insulators, second-order Moller-Plesset perturbation theory (MP2) significantly underestimates band gaps. Here, we evaluate the band gap predictions of partitioned equation-of-motion MP2 (P-EOM-MP2), which is a second-order approximation to equation-of-motion coupled-cluster theory with single and double excitations. On a test set of elemental and binary semiconductors and insulators, we find that P-EOM-MP2 overestimates band gaps by 0.3 eV on average, which can be compared to the underestimation by 0.6 eV on average exhibited by the G0W0 approximation with a PBE reference. We show that P-EOM-MP2, when interpreted as a Green's function-based theory, has a self-energy that includes all first- and second- order diagrams and a few third-order diagrams. We find that the GW approximation performs better for materials with small gaps and P-EOM-MP2 performs better for materials with large gaps, which we attribute to their superior treatment of screening and exchange, respectively. △ Less

Submitted 24 May, 2021; originally announced May 2021.

Comments: 4+ pages, 1 table, 3 figures

Approximating matrix eigenvalues by subspace iteration with repeated random sparsification

Authors: Samuel M. Greene, Robert J. Webber, Timothy C. Berkelbach, Jonathan Weare

Abstract: Traditional numerical methods for calculating matrix eigenvalues are prohibitively expensive for high-dimensional problems. Iterative random sparsification methods allow for the estimation of a single dominant eigenvalue at reduced cost by leveraging repeated random sampling and averaging. We present a general approach to extending such methods for the estimation of multiple eigenvalues and demons… ▽ More Traditional numerical methods for calculating matrix eigenvalues are prohibitively expensive for high-dimensional problems. Iterative random sparsification methods allow for the estimation of a single dominant eigenvalue at reduced cost by leveraging repeated random sampling and averaging. We present a general approach to extending such methods for the estimation of multiple eigenvalues and demonstrate its performance for several benchmark problems in quantum chemistry. △ Less

Submitted 2 March, 2022; v1 submitted 22 March, 2021; originally announced March 2021.

Comments: 31 pages, 8 figures

Dynamical correlation energy of metals in large basis sets from downfolding and composite approaches

Authors: James M. Callahan, Malte F. Lange, Timothy C. Berkelbach

Abstract: Coupled-cluster theory with single and double excitations (CCSD) is a promising ab initio method for the electronic structure of three-dimensional metals, for which second-order perturbation theory (MP2) diverges in the thermodynamic limit. However, due to the high cost and poor convergence of CCSD with respect to basis size, applying CCSD to periodic systems often leads to large basis set errors.… ▽ More Coupled-cluster theory with single and double excitations (CCSD) is a promising ab initio method for the electronic structure of three-dimensional metals, for which second-order perturbation theory (MP2) diverges in the thermodynamic limit. However, due to the high cost and poor convergence of CCSD with respect to basis size, applying CCSD to periodic systems often leads to large basis set errors. In a common "composite" method, MP2 is used to recover the missing dynamical correlation energy through a focal-point correction, but the inadequacy of MP2 for metals raises questions about this approach. Here we describe how high-energy excitations treated by MP2 can be "downfolded" into a low-energy active space to be treated by CCSD. Comparing how the composite and downfolding approaches perform for the uniform electron gas, we find that the latter converges more quickly with respect to the basis set size. Nonetheless, the composite approach is surprisingly accurate because it removes the problematic MP2 treatment of double excitations near the Fermi surface. Using the method to estimate the CCSD correlation energy in the combined complete basis set and thermodynamic limits, we find CCSD recovers over 90% of the exact correlation energy at $r_s=4$. We also test the composite and downfolding approaches with the random-phase approximation used in place of MP2, yielding a method that is more effective but more expensive. △ Less

Submitted 10 March, 2021; originally announced March 2021.

Comments: 7 pages, 4 figures

Simulations of Trions and Biexcitons in Layered Hybrid Organic-Inorganic Lead Halide Perovskites

Authors: Yeongsu Cho, Samuel M. Greene, Timothy C. Berkelbach

Abstract: Behaving like atomically-precise two-dimensional quantum wells with non-negligible dielectric contrast, the layered HOIPs have strong electronic interactions leading to tightly bound excitons with binding energies on the order of 500 meV. These strong interactions suggest the possibility of larger excitonic complexes like trions and biexcitons, which are hard to study numerically due to the comple… ▽ More Behaving like atomically-precise two-dimensional quantum wells with non-negligible dielectric contrast, the layered HOIPs have strong electronic interactions leading to tightly bound excitons with binding energies on the order of 500 meV. These strong interactions suggest the possibility of larger excitonic complexes like trions and biexcitons, which are hard to study numerically due to the complexity of the layered HOIPs. Here, we propose and parameterize a model Hamiltonian for excitonic complexes in layered HOIPs and we study the correlated eigenfunctions of trions and biexcitons using a combination of diffusion Monte Carlo and very large variational calculations with explicitly correlated Gaussian basis functions. Binding energies and spatial structures of these complexes are presented as a function of the layer thickness. The trion and biexciton of the thinnest layered HOIP have binding energies of 35 meV and 44 meV, respectively, whereas a single exfoliated layer is predicted to have trions and biexcitons with equal binding enegies of 48 meV. We compare our findings to available experimental data and to that of other quasi-two-dimensional materials. △ Less

Submitted 20 October, 2020; originally announced October 2020.

Comments: main text: 6 pages, 4 figures, 1 table; supplement: 7 pages, 2 figures

Journal ref: Phys. Rev. Lett. 126, 216402 (2021)

Regional Embedding Enables High-Level Quantum Chemistry for Surface Science

Authors: Bryan T. G. Lau, Gerald Knizia, Timothy C. Berkelbach

Abstract: Compared to common density functionals, ab initio wave function methods can provide greater reliability and accuracy, which could prove useful when modeling adsorbates or defects of otherwise periodic systems. However, the breaking of translational symmetry necessitates large supercells that are often prohibitive for correlated wave function methods. As an alternative, we introduce the regional em… ▽ More Compared to common density functionals, ab initio wave function methods can provide greater reliability and accuracy, which could prove useful when modeling adsorbates or defects of otherwise periodic systems. However, the breaking of translational symmetry necessitates large supercells that are often prohibitive for correlated wave function methods. As an alternative, we introduce the regional embedding approach, which enables correlated wave function treatments of only a target fragment of interest through small, fragment-localized orbital spaces constructed using a simple overlap criterion. Applications to the adsorption of water on lithium hydride, hexagonal boron nitride, and graphene substrates show that regional embedding combined with focal point corrections can provide converged CCSD(T) (coupled cluster) adsorption energies with very small fragment sizes. △ Less

Submitted 1 October, 2020; originally announced October 2020.

Full-Frequency GW without Frequency

Authors: Sylvia J. Bintrim, Timothy C. Berkelbach

Abstract: Efficient computer implementations of the GW approximation must approximate a numerically challenging frequency integral; the integral can be performed analytically, but doing so leads to an expensive implementation whose computational cost scales as $O(N^6)$ where $N$ is the size of the system. Here we introduce a new formulation of the full-frequency GW approximation by exactly recasting it as a… ▽ More Efficient computer implementations of the GW approximation must approximate a numerically challenging frequency integral; the integral can be performed analytically, but doing so leads to an expensive implementation whose computational cost scales as $O(N^6)$ where $N$ is the size of the system. Here we introduce a new formulation of the full-frequency GW approximation by exactly recasting it as an eigenvalue problem in an expanded space. This new formulation (1) avoids the use of time or frequency grids, (2) naturally precludes the common "diagonal" approximation, (3) enables common iterative eigensolvers that reduce the canonical scaling to $O(N^5)$, and (4) enables a density-fitted implementation that reduces the scaling to $O(N^4)$. We numerically verify these scaling behaviors and test a variety of approximations that are motivated by this new formulation. In this new formulation, the relation of the GW approximation to configuration interaction, coupled-cluster theory, and the algebraic diagrammatic construction is made especially apparent, providing a new direction for improvements to the GW approximation. △ Less

Submitted 29 September, 2020; originally announced September 2020.

Comments: 5+2 pages, 3 figures

Improved Fast Randomized Iteration Approach to Full Configuration Interaction

Authors: Samuel M. Greene, Robert J. Webber, Jonathan Weare, Timothy C. Berkelbach

Abstract: We present three modifications to our recently introduced fast randomized iteration method for full configuration interaction (FCI-FRI) and investigate their effects on the method's performance for Ne, H$_2$O, and N$_2$. The initiator approximation, originally developed for full configuration interaction quantum Monte Carlo, significantly reduces statistical error in FCI-FRI when few samples are u… ▽ More We present three modifications to our recently introduced fast randomized iteration method for full configuration interaction (FCI-FRI) and investigate their effects on the method's performance for Ne, H$_2$O, and N$_2$. The initiator approximation, originally developed for full configuration interaction quantum Monte Carlo, significantly reduces statistical error in FCI-FRI when few samples are used in compression operations, enabling its application to larger chemical systems. The semi-stochastic extension, which involves exactly preserving a fixed subset of elements in each compression, improves statistical efficiency in some cases but reduces it in others. We also developed a new approach to sampling excitations that yields consistent improvements in statistical efficiency and reductions in computational cost. We discuss possible strategies based on our findings for improving the performance of stochastic quantum chemistry methods more generally. △ Less

Submitted 20 July, 2020; v1 submitted 1 May, 2020; originally announced May 2020.

Comments: 13 pages, 5 figures

Excitons in solids from periodic equation-of-motion coupled-cluster theory

Authors: Xiao Wang, Timothy C. Berkelbach

Abstract: We present an ab initio study of electronically excited states of three-dimensional solids using Gaussian-based periodic equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD). The explicit use of translational symmetry, as implemented via Brillouin zone sampling and momentum conservation, is responsible for a large reduction in cost. Our largest system studied, wh… ▽ More We present an ab initio study of electronically excited states of three-dimensional solids using Gaussian-based periodic equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD). The explicit use of translational symmetry, as implemented via Brillouin zone sampling and momentum conservation, is responsible for a large reduction in cost. Our largest system studied, which samples the Brillouin zone using 64 k-points (a 4x4x4 mesh) corresponds to a canonical EOM-CCSD calculation of 768 electrons in 640 orbitals. We study eight simple semiconductors and insulators, with direct singlet excitation energies in the range of 3 to 15 eV. Our predicted excitation energies exhibit a mean absolute error of 0.27 eV when compared to experiment. We furthermore calculate the energy of excitons with nonzero momentum and compare the exciton dispersion of LiF with experimental data from inelastic X-ray scattering. By calculating excitation energies under strain, we extract hydrostatic deformation potentials in order to quantify the strength of interactions between excitons and acoustic phonons. Our results indicate that coupled-cluster theory is a promising method for the accurate study of a variety of exciton phenomena in solids. △ Less

Submitted 29 January, 2020; originally announced January 2020.

Comments: 8 pages, 5 figures

A unification of the Holstein polaron and dynamic disorder pictures of charge transport in organic semiconductors

Authors: Jonathan H. Fetherolf, Denis Golez, Timothy C. Berkelbach

Abstract: We present a unified and nonperturbative method for calculating spectral and transport properties of Hamiltonians with simultaneous Holstein (diagonal) and Peierls (off-diagonal) electron-phonon coupling. Our approach is motivated by the separation of energy scales in semiconducting organic molecular cystals, in which electrons couple to high-frequency intramolecular Holstein modes and to low-freq… ▽ More We present a unified and nonperturbative method for calculating spectral and transport properties of Hamiltonians with simultaneous Holstein (diagonal) and Peierls (off-diagonal) electron-phonon coupling. Our approach is motivated by the separation of energy scales in semiconducting organic molecular cystals, in which electrons couple to high-frequency intramolecular Holstein modes and to low-frequency intermolecular Peierls modes. We treat Peierls modes as quasi-classical dynamic disorder, while Holstein modes are included with a Lang-Firsov polaron transformation and no narrow-band approximation. Our method reduces to the popular polaron picture due to Holstein coupling and the dynamic disorder picture due to Peierls coupling. We derive an expression for efficient numerical evaluation of the frequency-resolved optical conductivity based on the Kubo formula and obtain the DC mobility from its zero-frequency component. We also use our method to calculate the electron-addition Green's function corresponding to the inverse photoemission spectrum. For realistic parameters, temperature-dependent DC mobility is largely determined by the Peierls-induced dynamic disorder with minor quantitative corrections due to polaronic band-narrowing, and an activated regime is not observed at relevant temperatures. In contrast, for frequency-resolved observables, a quantum mechanical treatment of the Holstein coupling is qualitatively important for capturing the phonon replica satellite structure. △ Less

Submitted 10 October, 2019; originally announced October 2019.

Comments: 12 pages, 7 figures

Journal ref: Phys. Rev. X 10, 021062 (2020)

Ab Initio Linear and Pump-Probe Spectroscopy of Excitons in Molecular Crystals

Authors: Alan M. Lewis, Timothy C. Berkelbach

Abstract: Linear and non-linear spectroscopies are powerful tools used to investigate the energetics and dynamics of electronic excited states of both molecules and crystals. While highly accurate \emph{ab initio} calculations of molecular spectra can be performed relatively routinely, extending these calculations to periodic systems is challenging. Here, we present calculations of the linear absorption spe… ▽ More Linear and non-linear spectroscopies are powerful tools used to investigate the energetics and dynamics of electronic excited states of both molecules and crystals. While highly accurate \emph{ab initio} calculations of molecular spectra can be performed relatively routinely, extending these calculations to periodic systems is challenging. Here, we present calculations of the linear absorption spectrum and pump-probe two-photon photoemission spectra of the naphthalene crystal using equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD). Molecular acene crystals are of interest due to the low-energy multi-exciton singlet states they exhibit, which have been studied extensively as intermediates involved in singlet fission. Our linear absorption spectrum is in good agreement with experiment, predicting a first exciton absorption peak at 4.4 eV, and our two-photon photoemission spectra capture the behavior of multi-exciton states, whose double-excitation character cannot be captured by current methods. The simulated pump-probe spectra provide support for existing interpretations of two-photon photoemission in closely-related acene crystals such as tetracene and pentacene. △ Less

Submitted 24 September, 2019; originally announced September 2019.

Comments: 6 pages, 2 figures

Thickness-dependent optical properties of layered hybrid organic-inorganic halide perovskites: A tight-binding GW-BSE study

Authors: Yeongsu Cho, Timothy C. Berkelbach

Abstract: We present a many-body calculation of the band structure and optical spectrum of the layered hybrid organic-inorganic halide perovskites in the Ruddlesden-Popper phase with the general formula A$^{'}_{2}$A$_{n-1}$M$_{n}$X$_{3n+1}$, focusing specifically on the lead iodide family. We calculate the mean-field band structure with spin-orbit coupling, quasiparticle corrections within the GW approximat… ▽ More We present a many-body calculation of the band structure and optical spectrum of the layered hybrid organic-inorganic halide perovskites in the Ruddlesden-Popper phase with the general formula A$^{'}_{2}$A$_{n-1}$M$_{n}$X$_{3n+1}$, focusing specifically on the lead iodide family. We calculate the mean-field band structure with spin-orbit coupling, quasiparticle corrections within the GW approximation, and optical spectra using the Bethe-Salpeter equation. The model is parameterized by first-principles calculations and classical electrostatic screening, enabling an accurate but cost-effective study of large unit cells and corresponding thickness-dependent properties. A transition of the electronic and optical properties from quasi-two-dimensional behavior to three-dimensional behavior is shown for increasing $n$ and the nonhydrogenic character of the excitonic Rydberg series is analyzed. The thickness-dependent 1s and 2s exciton energy levels are in good agreement with recently reported experiments and the 1s exciton binding energy is calculated to be 302 meV for $n=1$, 97 meV for $n=5$, and 37 meV for $n=\infty$ (bulk MAPbI3). △ Less

Submitted 25 August, 2019; originally announced August 2019.

Comments: 8 pages, 5 figures

Quantum plasmons and intraband excitons in doped nanoparticles: Failure of the Tamm-Dancoff approximation and importance of electron-hole attraction

Authors: Bryan T. G. Lau, Timothy C. Berkelbach

Abstract: We use excited-state quantum chemistry techniques to investigate the intraband absorption of doped semiconductor nanoparticles as a function of doping density, nanoparticle radius, and material properties. The excess electrons are modeled as interacting particles confined in a sphere. We compare the predictions of various single-excitation theories, including time-dependent Hartree-Fock, the rando… ▽ More We use excited-state quantum chemistry techniques to investigate the intraband absorption of doped semiconductor nanoparticles as a function of doping density, nanoparticle radius, and material properties. The excess electrons are modeled as interacting particles confined in a sphere. We compare the predictions of various single-excitation theories, including time-dependent Hartree-Fock, the random-phase approximation, and configuration interaction with single excitations. We find that time-dependent Hartree-Fock most accurately describes the character of the excitation, as compared to equation-of-motion coupled-cluster theory with single and double excitations. The excitation evolves from confinement-dominated, to excitonic, to plasmonic with increasing number of electrons at fixed density, and the threshold number of electrons to produce a plasmon increases with density due to quantum confinement. Exchange integrals (attractive electron-hole interactions) are essential to properly describe excitons, and de-excitations (i.e.~avoidance of the Tamm-Dancoff approximation) are essential to properly describe plasmons. We propose a schematic model whose analytic solutions closely reproduce our numerical calculations. Our results are in good agreement with experimental spectra of doped ZnO nanoparticles at a doping density of $1.4\times 10^{20}$ cm$^{-3}$. △ Less

Submitted 25 July, 2019; originally announced July 2019.

Comments: 10 pages, 5 figures

Coupled-cluster impurity solvers for dynamical mean-field theory

Authors: Tianyu Zhu, Carlos A. Jimenez-Hoyos, James McClain, Timothy C. Berkelbach, Garnet Kin-Lic Chan

Abstract: We describe the use of coupled-cluster theory as an impurity solver in dynamical mean-field theory (DMFT) and its cluster extensions. We present numerical results at the level of coupled-cluster theory with single and double excitations (CCSD) for the density of states and self-energies of cluster impurity problems in the one- and two-dimensional Hubbard models. Comparison to exact diagonalization… ▽ More We describe the use of coupled-cluster theory as an impurity solver in dynamical mean-field theory (DMFT) and its cluster extensions. We present numerical results at the level of coupled-cluster theory with single and double excitations (CCSD) for the density of states and self-energies of cluster impurity problems in the one- and two-dimensional Hubbard models. Comparison to exact diagonalization shows that CCSD produces accurate density of states and self-energies at a variety of values of $U/t$ and filling fractions. However, the low cost allows for the use of many bath sites, which we define by a discretization of the hybridization directly on the real frequency axis. We observe convergence of dynamical quantities using approximately 30 bath sites per impurity site, with our largest 4-site cluster DMFT calculation using 120 bath sites. We suggest coupled cluster impurity solvers will be attractive in ab initio formulations of dynamical mean-field theory. △ Less

Submitted 19 August, 2019; v1 submitted 28 May, 2019; originally announced May 2019.

Comments: 9 pages, 5 figures

Journal ref: Phys. Rev. B 100, 115154 (2019)

Beyond Walkers in Stochastic Quantum Chemistry: Reducing Error using Fast Randomized Iteration

Authors: Samuel M. Greene, Robert J. Webber, Jonathan Weare, Timothy C. Berkelbach

Abstract: We introduce a family of methods for the full configuration interaction problem in quantum chemistry, based on the fast randomized iteration (FRI) framework [L.-H. Lim and J. Weare, SIAM Rev. 59, 547 (2017)]. These methods, which we term "FCI-FRI," stochastically impose sparsity during iterations of the power method and can be viewed as a generalization of full configuration interaction quantum Mo… ▽ More We introduce a family of methods for the full configuration interaction problem in quantum chemistry, based on the fast randomized iteration (FRI) framework [L.-H. Lim and J. Weare, SIAM Rev. 59, 547 (2017)]. These methods, which we term "FCI-FRI," stochastically impose sparsity during iterations of the power method and can be viewed as a generalization of full configuration interaction quantum Monte Carlo (FCIQMC) without walkers. In addition to the multinomial scheme commonly used to sample excitations in FCIQMC, we present a systematic scheme where excitations are not sampled independently. Performing ground-state calculations on five small molecules at fixed cost, we find that the systematic FCI-FRI scheme is 11 to 45 times more statistically efficient than the multinomial FCI-FRI scheme, which is in turn 1.4 to 178 times more statistically efficient than the original FCIQMC algorithm. △ Less

Submitted 9 July, 2019; v1 submitted 2 May, 2019; originally announced May 2019.

Comments: 19 pages, 7 figures

Ab Initio Lifetime and Concomitant Double-Excitation Character of Plasmons at Metallic Densities

Authors: Alan M. Lewis, Timothy C. Berkelbach

Abstract: The accurate calculation of excited state properties of interacting electrons in the condensed phase is an immense challenge in computational physics. Here, we use state-of-the-art equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD) to calculate the dynamic structure factor, which can be experimentally measured by inelastic x-ray and electron scattering. Our cal… ▽ More The accurate calculation of excited state properties of interacting electrons in the condensed phase is an immense challenge in computational physics. Here, we use state-of-the-art equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD) to calculate the dynamic structure factor, which can be experimentally measured by inelastic x-ray and electron scattering. Our calculations are performed on the uniform electron gas at densities corresponding to Wigner-Seitz radii of $r_s=5$, 4, and 3 corresponding to the valence electron densities of common metals. We compare our results to those obtained using the random-phase approximation, which is known to provide a reasonable description of the collective plasmon excitation and which resums only a small subset of the polarizability diagrams included in EOM-CCSD. We find that EOM-CCSD, instead of providing a perturbative improvement on the RPA plasmon, predicts a many-state plasmon resonance, where each contributing state has a double-excitation character of 80\% or more. This finding amounts to an ab initio treatment of the plasmon linewidth, which is in good quantitative agreement with previous diagrammatic calculations, and highlights the strongly correlated nature of lifetime effects in condensed-phase electronic structure theory. △ Less

Submitted 8 July, 2019; v1 submitted 15 October, 2018; originally announced October 2018.

Comments: 6 pages, 4 figures; accepted version

Journal ref: Phys. Rev. Lett. 122, 226402 (2019)

Many-Body Simulation of Two-Dimensional Electronic Spectroscopy of Excitons and Trions in Monolayer Transition-Metal Dichalcogenides

Authors: Roel Tempelaar, Timothy C. Berkelbach

Abstract: We present a many-body formalism for the simulation of time-resolved nonlinear spectroscopy and apply it to study the coherent interaction between excitons and trions in doped transition-metal dichalcogenides. Although the formalism can be straightforwardly applied in a first-principles manner, for simplicity we use a parameterized band structure and a static model dielectric function, both of whi… ▽ More We present a many-body formalism for the simulation of time-resolved nonlinear spectroscopy and apply it to study the coherent interaction between excitons and trions in doped transition-metal dichalcogenides. Although the formalism can be straightforwardly applied in a first-principles manner, for simplicity we use a parameterized band structure and a static model dielectric function, both of which can be obtained from a calculation using the $GW$ approximation. Our simulation results shed light on the interplay between singlet and triplet trions in molybdenum- and tungsten-based compounds. Our two-dimensional electronic spectra are in excellent agreement with recent experiments and we accurately reproduce the beating of a cross-peak signal indicative of quantum coherence between excitons and trions. Although we confirm that the quantum beats in molybdenum-based monolayers unambigously reflect the exciton-trion coherence time, they are shown here to provide a lower-bound to the coherence time of tungsten analogues due to a destructive interference emerging from coexisting singlet and triplet trions. △ Less

Submitted 15 October, 2018; originally announced October 2018.

Vertex corrections to the polarizability do not improve the GW approximation for the ionization potential of molecules

Authors: Alan M. Lewis, Timothy C. Berkelbach

Abstract: The $GW$ approximation is based on the neglect of vertex corrections, which appear in the exact self-energy and the exact polarizability. Here, we investigate the importance of vertex corrections in the polarizability only. We calculate the polarizability with equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD), which rigorously includes a large class of diagram… ▽ More The $GW$ approximation is based on the neglect of vertex corrections, which appear in the exact self-energy and the exact polarizability. Here, we investigate the importance of vertex corrections in the polarizability only. We calculate the polarizability with equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD), which rigorously includes a large class of diagrammatically-defined vertex corrections beyond the random phase approximation (RPA). As is well-known, the frequency-dependent polarizability predicted by EOM-CCSD is quite different and generally more accurate than that predicted by the RPA. We evaluate the effect of these vertex corrections on a test set of 20 atoms and molecules. When using a Hartree-Fock reference, ionization potentials predicted by the $GW$ approximation with the RPA polarizability are typically overestimated with a mean absolute error of 0.3 eV. However, those predicted with a vertex-corrected polarizability are typically underestimated with an increased mean absolute error of 0.5 eV. This result suggests that vertex corrections in the self-energy cannot be neglected, at least for molecules. We also assess the behavior of eigenvalue self-consistency in vertex-corrected $GW$ calculations, finding a further worsening of the predicted ionization potentials. △ Less

Submitted 12 April, 2019; v1 submitted 30 September, 2018; originally announced October 2018.

Comments: 8 pages, 4 figures; results and discussion revised

On The Relation Between Equation-of-Motion Coupled-Cluster Theory and the GW Approximation

Authors: Malte F. Lange, Timothy C. Berkelbach

Abstract: We discuss the analytic and diagrammatic structure of ionization potential (IP) and electron affinity (EA) equation-of-motion coupled-cluster (EOM-CC) theory, in order to put it on equal footing with the prevalent $GW$ approximation. The comparison is most straightforward for the time-ordered one-particle Green's function, and we show that the Green's function calculated by EOM-CC with single and… ▽ More We discuss the analytic and diagrammatic structure of ionization potential (IP) and electron affinity (EA) equation-of-motion coupled-cluster (EOM-CC) theory, in order to put it on equal footing with the prevalent $GW$ approximation. The comparison is most straightforward for the time-ordered one-particle Green's function, and we show that the Green's function calculated by EOM-CC with single and double excitations (EOM-CCSD) includes fewer ring diagrams at higher order than does the $GW$ approximation, due to the former's unbalanced treatment of time-ordering. However, the EOM-CCSD Green's function contains a large number of vertex corrections, including ladder diagrams, mixed ring-ladder diagrams, and exchange diagrams. By including triple excitations, the EOM-CCSDT Green's function includes all diagrams contained in the $GW$ approximation, along with many high-order vertex corrections. In the same language, we discuss a number of common approximations to the EOM-CCSD equations, many of which can be classified as elimination of diagrams. Finally, we present numerical results by calculating the principal charged excitations energies of the molecules contained in the so-called $GW$100 test set [J. Chem. Theory Comput. 2015, 11, 5665-5687]. We argue that (in molecules) exchange is as important as screening, advocating for a Hartree-Fock reference and second-order exchange in the self-energy. △ Less

Submitted 30 April, 2018; originally announced May 2018.

Comments: 13 pages, 5 figures

Journal ref: J. Chem. Theory Comput., 14, 4224 (2018)

Random-phase approximation excitation energies from approximate equation-of-motion ring coupled-cluster doubles

Authors: Timothy C. Berkelbach

Abstract: The ground-state correlation energy calculated in the random-phase approximation (RPA) is known to be identical to that calculated using a subset of terms appearing in coupled-cluster theory with double excitations. In particular, this equivalence requires keeping only those terms that generate time-independent ring diagrams, in the Goldstone sense. Here I show that this equivalence extends to neu… ▽ More The ground-state correlation energy calculated in the random-phase approximation (RPA) is known to be identical to that calculated using a subset of terms appearing in coupled-cluster theory with double excitations. In particular, this equivalence requires keeping only those terms that generate time-independent ring diagrams, in the Goldstone sense. Here I show that this equivalence extends to neutral excitation energies, for which those calculated in the RPA are identical to those calculated using an approximation to equation-of-motion coupled-cluster theory with double excitations (EOM-CCD). The equivalence requires three approximations to EOM-CCD: first, the ground-state double-excitation amplitudes are obtained from the ring-CCD equations (the same as for the correlation energy); second, the EOM eigenvalue problem is truncated to the single-excitation (one particle + one hole) subspace; third, the similarity transformation of the Fock operator must be neglected, as it corresponds to a dressing of the single-particle propagator, which is not present in the conventional RPA. △ Less

Submitted 3 April, 2018; originally announced April 2018.

Comments: 4 pages, 2 figures

Journal ref: J. Chem. Phys. 149, 041103 (2018)

Linear and nonlinear spectroscopy from quantum master equations

Authors: Jonathan H. Fetherolf, Timothy C. Berkelbach

Abstract: We investigate the accuracy of the second-order time-convolutionless (TCL2) quantum master equation for the calculation of linear and nonlinear spectroscopies of multichromophore systems. We show that, even for systems with non-adiabatic coupling, the TCL2 master equation predicts linear absorption spectra that are accurate over an extremely broad range of parameters and well beyond what would be… ▽ More We investigate the accuracy of the second-order time-convolutionless (TCL2) quantum master equation for the calculation of linear and nonlinear spectroscopies of multichromophore systems. We show that, even for systems with non-adiabatic coupling, the TCL2 master equation predicts linear absorption spectra that are accurate over an extremely broad range of parameters and well beyond what would be expected based on the perturbative nature of the approach; non-equilibrium population dynamics calculated with TCL2 for identical parameters are significantly less accurate. For third-order (two-dimensional) spectroscopy, the importance of population dynamics and the violation of the so-called quantum regression theorem degrade the accuracy of TCL2 dynamics. To correct these failures, we combine the TCL2 approach with a classical ensemble sampling of slow microscopic bath degrees of freedom, leading to an efficient hybrid quantum-classical scheme that displays excellent accuracy over a wide range of parameters. In the spectroscopic setting, the success of such a hybrid scheme can be understood through its separate treatment of homogeneous and inhomogeneous broadening. Importantly, the presented approach has the computational scaling of TCL2, with the modest addition of an embarrassingly parallel prefactor associated with ensemble sampling. The presented approach can be understood as a generalized inhomogeneous cumulant expansion technique, capable of treating multilevel systems with non-adiabatic dynamics. △ Less

Submitted 12 September, 2017; originally announced September 2017.

Comments: 9 pages, 4 figures

Environmentally-Sensitive Theory of Electronic and Optical Transitions in Atomically-Thin Semiconductors

Authors: Yeongsu Cho, Timothy C. Berkelbach

Abstract: We present an electrostatic theory of band gap renormalization in atomically-thin semiconductors that captures the strong sensitivity to the surrounding dielectric environment. In particular, our theory aims to correct known band gaps, such as that of the three-dimensional bulk crystal. Combining our quasiparticle band gaps with an effective mass theory of excitons yields environmentally-sensitive… ▽ More We present an electrostatic theory of band gap renormalization in atomically-thin semiconductors that captures the strong sensitivity to the surrounding dielectric environment. In particular, our theory aims to correct known band gaps, such as that of the three-dimensional bulk crystal. Combining our quasiparticle band gaps with an effective mass theory of excitons yields environmentally-sensitive optical gaps as would be observed in absorption or photoluminescence. For an isolated monolayer of MoS$_2$, the presented theory is in good agreement with ab initio results based on the GW approximation and the Bethe-Salpeter equation. We find that changes in the electronic band gap are almost exactly offset by changes in the exciton binding energy, such that the energy of the first optical transition is nearly independent of the electrostatic environment, rationalizing experimental observations. △ Less

Submitted 4 September, 2017; originally announced September 2017.

Comments: 5 pages, 3 figures

Journal ref: Phys. Rev. B 97, 041409 (2018)

Coulomb engineering of the bandgap in 2D semiconductors

Authors: Archana Raja, Andrey Chaves, Jaeeun Yu, Ghidewon Arefe, Heather M. Hill, Albert F. Rigosi, Timothy C. Berkelbach, Philipp Nagler, Christian Schüller, Tobias Korn, Colin Nuckolls, James Hone, Louis E. Brus, Tony F. Heinz, David R. Reichman, Alexey Chernikov

Abstract: The ability to control the size of the electronic bandgap is an integral part of solid-state technology. Atomically-thin two-dimensional crystals offer a new approach for tuning the energies of the electronic states based on the interplay between the environmental sensitivity and unusual strength of the Coulomb interaction in these materials. By engineering the surrounding dielectric environment,… ▽ More The ability to control the size of the electronic bandgap is an integral part of solid-state technology. Atomically-thin two-dimensional crystals offer a new approach for tuning the energies of the electronic states based on the interplay between the environmental sensitivity and unusual strength of the Coulomb interaction in these materials. By engineering the surrounding dielectric environment, we are able to tune the electronic bandgap in monolayers of WS2 and WSe2 by hundreds of meV. We exploit this behavior to present an in-plane dielectric heterostructure with a spatially dependent bandgap, illustrating the feasibility of our approach for the creation of lateral junctions with nanoscale resolution. This successful demonstration of bandgap engineering based on the non-invasive modification of the Coulomb interaction should enable the design of a new class of atomically thin devices to advance the limits of size and functionality for solid-state technologies. △ Less

Submitted 7 February, 2017; v1 submitted 3 February, 2017; originally announced February 2017.

Gaussian-based coupled-cluster theory for the ground state and band structure of solids

Authors: James McClain, Qiming Sun, Garnet Kin-Lic Chan, Timothy C. Berkelbach

Abstract: We present the results of Gaussian-based ground-state and excited-state equation-of-motion coupled-cluster theory with single and double excitations for three-dimensional solids. We focus on diamond and silicon, which are paradigmatic covalent semiconductors. In addition to ground-state properties (the lattice constant, bulk modulus, and cohesive energy), we compute the quasiparticle band structur… ▽ More We present the results of Gaussian-based ground-state and excited-state equation-of-motion coupled-cluster theory with single and double excitations for three-dimensional solids. We focus on diamond and silicon, which are paradigmatic covalent semiconductors. In addition to ground-state properties (the lattice constant, bulk modulus, and cohesive energy), we compute the quasiparticle band structure and band gap. We sample the Brillouin zone with up to 64 k-points using norm-conserving pseudopotentials and polarized double- and triple-zeta basis sets, leading to canonical coupled-cluster calculations with as many as 256 electrons in 2,176 orbitals. △ Less

Submitted 17 January, 2017; originally announced January 2017.

Comments: 10 pages, 6 figures

On the accuracy of the Pade-resummed master equation approach to dissipative quantum dynamics

Authors: Hsing-Ta Chen, Timothy C. Berkelbach, David R. Reichman

Abstract: Well--defined criteria are proposed for assessing the accuracy of quantum master equations whose memory functions are approximated by Padé resummation of the first two moments in the electronic coupling. These criteria partition the parameter space into distinct levels of expected accuracy, ranging from quantitatively accurate regimes to regions of parameter space where the approach is not expecte… ▽ More Well--defined criteria are proposed for assessing the accuracy of quantum master equations whose memory functions are approximated by Padé resummation of the first two moments in the electronic coupling. These criteria partition the parameter space into distinct levels of expected accuracy, ranging from quantitatively accurate regimes to regions of parameter space where the approach is not expected to be applicable. Extensive comparison of Padé--resummed master equations with numerically exact results in the context of the spin--boson model demonstrate that the proposed criteria correctly demarcate the regions of parameter space where the Padé approximation is reliable. The applicability analysis we present is not confined to the specifics of the Hamiltonian under consideration and should provide guidelines for other classes of resummation techniques. △ Less

Submitted 1 April, 2016; v1 submitted 11 March, 2016; originally announced March 2016.