-
Understanding Discrepancies of Wavefunction Theories for Large Molecules
Authors:
Tobias Schäfer,
Andreas Irmler,
Alejandro Gallo,
Andreas Grüneis
Abstract:
Quantum mechanical many-electron calculations can predict properties of atoms, molecules and even complex materials. The employed computational methods play a quintessential role in many scientifically and technologically relevant research fields. However, a question of paramount importance is whether approximations aimed at reducing the computational complexity for solving the many-electron Schrö…
▽ More
Quantum mechanical many-electron calculations can predict properties of atoms, molecules and even complex materials. The employed computational methods play a quintessential role in many scientifically and technologically relevant research fields. However, a question of paramount importance is whether approximations aimed at reducing the computational complexity for solving the many-electron Schrödinger equation, are accurate enough. Here, we investigate recently reported discrepancies of noncovalent interaction energies for large molecules predicted by two of the most widely-trusted many-electron theories: diffusion quantum Monte Carlo and coupled-cluster theory. We are able to unequivocally pin down the source of the puzzling discrepancies and present modifications to widely-used coupled-cluster methods needed for more accurate noncovalent interaction energies of large molecules on the hundred-atom scale. This enhances the reliability of predictions from quantum mechanical many-electron theories across a wide range of critical applications, including drug design, catalysis, and the innovation of new functional materials, such as those for renewable energy technologies.
△ Less
Submitted 4 July, 2024; v1 submitted 1 July, 2024;
originally announced July 2024.
-
Sampling the reciprocal Coulomb potential in finite and anisotropic cells
Authors:
Tobias Schäfer,
William Z. Van Benschoten,
James J. Shepherd,
Andreas Grüneis
Abstract:
We present a robust strategy to numerically sample the Coulomb potential in reciprocal space for periodic Born-von Karman cells of general shape. Our approach tackles two common issues of plane-wave based implementations of Coulomb integrals under periodic boundary conditions, the treatment of the singularity at the Brillouin-zone center, as well as quadrature errors, which can cause severe conver…
▽ More
We present a robust strategy to numerically sample the Coulomb potential in reciprocal space for periodic Born-von Karman cells of general shape. Our approach tackles two common issues of plane-wave based implementations of Coulomb integrals under periodic boundary conditions, the treatment of the singularity at the Brillouin-zone center, as well as quadrature errors, which can cause severe convergence problems in anisotropic cells, necessary for the calculation of low-dimensional systems. We apply our strategy to the Hartree-Fock (HF) and coupled cluster (CC) theory and discuss the consequences of different sampling strategies on the different theories. We show that sampling the Coulomb potential via the widely used probe-charge Ewald method is unsuitable for CC calculations in anisotropic cells. To demonstrate the applicability of our developed approach, we study two representative, low-dimensional use cases: the infinite carbon chain, for which we report the first periodic CCSD(T) potential energy surface, as well as a surface slab of lithium hydride, for which we demonstrate the impact of different sampling strategies for calculating surface energies. We find that our Coulomb sampling strategy serves as a vital solution, addressing the critical need for improved accuracy in plane-wave based CC calculations for low-dimensional systems.
△ Less
Submitted 23 October, 2023; v1 submitted 20 October, 2023;
originally announced October 2023.
-
Phaseless auxiliary field quantum Monte Carlo with projector-augmented wave method for solids
Authors:
Amir Taheridehkordi,
Martin Schlipf,
Zoran Sukurma,
Moritz Humer,
Andreas Grüneis,
Georg Kresse
Abstract:
We implement the phaseless auxiliary field quantum Monte Carlo method using the plane-wave based projector augmented wave method and explore the accuracy and the feasibility of applying our implementation to solids. We use a singular value decomposition to compress the two-body Hamiltonian and thus reduce the computational cost. Consistent correlation energies from the primitive-cell sampling and…
▽ More
We implement the phaseless auxiliary field quantum Monte Carlo method using the plane-wave based projector augmented wave method and explore the accuracy and the feasibility of applying our implementation to solids. We use a singular value decomposition to compress the two-body Hamiltonian and thus reduce the computational cost. Consistent correlation energies from the primitive-cell sampling and the corresponding supercell calculations numerically verify our implementation. We calculate the equation of state for diamond and the correlation energies for a range of prototypical solid materials. A down-sampling technique along with natural orbitals accelerates the convergence with respect to the number of orbitals and crystal momentum points. We illustrate the competitiveness of our implementation in accuracy and computational cost for dense crystal momentum point meshes comparing to a well-established quantum-chemistry approach, the coupled-cluster ansatz including singles, doubles and perturbative triple particle-hole excitation operators.
△ Less
Submitted 27 April, 2023;
originally announced April 2023.
-
Coupled cluster theory for the ground and excited states of two dimensional quantum dots
Authors:
Faruk Salihbegović,
Alejandro Gallo,
Andreas Grüneis
Abstract:
We present a study of the two dimensional circular quantum dot model Hamiltonian using a range of quantum chemical ab initio methods. Ground and excited state energies are computed on different levels of perturbation theories including the coupled cluster method. We outline a scheme to compute the required Coulomb integrals in real space and utilize a semi-analytic solution to the integral over th…
▽ More
We present a study of the two dimensional circular quantum dot model Hamiltonian using a range of quantum chemical ab initio methods. Ground and excited state energies are computed on different levels of perturbation theories including the coupled cluster method. We outline a scheme to compute the required Coulomb integrals in real space and utilize a semi-analytic solution to the integral over the Coulomb kernel in the vicinity of the singularity. Furthermore, we show that the remaining basis set incompleteness error for two dimensional quantum dots scales with the inverse number of virtual orbitals, allowing us to extrapolate to the complete basis set limit energy. By varying the harmonic potential parameter we tune the correlation strength and investigate the predicted ground and excited state energies.
△ Less
Submitted 11 November, 2021;
originally announced November 2021.
-
Surface science using coupled cluster theory via local Wannier functions and in-RPA-embedding: the case of water on graphitic carbon nitride
Authors:
Tobias Schäfer,
Alejandro Gallo,
Andreas Irmler,
Felix Hummel,
Andreas Grüneis
Abstract:
A first-principles study of the adsorption of a single water molecule on a layer of graphitic carbon nitride employing an embedding approach is presented. The embedding approach involves an algorithm to obtain localized Wannier orbitals of various types expanded in a plane-wave basis and intrinsic atomic orbital projectors. The localized occupied orbitals are employed in combination with unoccupie…
▽ More
A first-principles study of the adsorption of a single water molecule on a layer of graphitic carbon nitride employing an embedding approach is presented. The embedding approach involves an algorithm to obtain localized Wannier orbitals of various types expanded in a plane-wave basis and intrinsic atomic orbital projectors. The localized occupied orbitals are employed in combination with unoccupied natural orbitals to perform many-electron perturbation theory calculations of local fragments. The fragments are comprised of a set of localized orbitals close to the adsorbed water molecule. Although the surface model contains more than 100 atoms in the simulation cell, the employed fragments are small enough to allow for calculations using high-level theories up to the coupled cluster ansatz with single, double and perturbative triple particle-hole excitation operators (CCSD(T)). To correct for the missing long-range correlation energy contributions to the adsorption energy, we embed CCSD(T) theory into the direct random phase approximation, yielding rapidly convergent adsorption energies with respect to the fragment size. Convergence of computed binding energies with respect to the virtual orbital basis set is achieved employing a number of recently developed techniques. Moreover, we discuss fragment size convergence for a range of approximate many-electron perturbation theories. The obtained benchmark results are compared to a number of density functional calculations.
△ Less
Submitted 14 October, 2021; v1 submitted 12 October, 2021;
originally announced October 2021.
-
Driven electronic bridge processes via defect states in $^{229}$Th-doped crystals
Authors:
Brenden S. Nickerson,
Martin Pimon,
Pavlo V. Bilous,
Johannes Gugler,
Georgy A. Kazakov,
Tomas Sikorsky,
Kjeld Beeks,
Andreas Gruneis,
Thorsten Schumm,
Adriana Palffy
Abstract:
The electronic defect states resulting from doping $^{229}$Th in CaF$_2$ offer a unique opportunity to excite the nuclear isomeric state $^{229m}$Th at approximately 8 eV via electronic bridge mechanisms. We consider bridge schemes involving stimulated emission and absorption using an optical laser. The role of different multipole contributions, both for the emitted or absorbed photon and nuclear…
▽ More
The electronic defect states resulting from doping $^{229}$Th in CaF$_2$ offer a unique opportunity to excite the nuclear isomeric state $^{229m}$Th at approximately 8 eV via electronic bridge mechanisms. We consider bridge schemes involving stimulated emission and absorption using an optical laser. The role of different multipole contributions, both for the emitted or absorbed photon and nuclear transition, to the total bridge rates are investigated theoretically. We show that the electric dipole component is dominant for the electronic bridge photon. In contradistinction, the electric quadrupole channel of the $^{229}$Th isomeric transition plays the dominant role for the bridge processes presented. The driven bridge rates are discussed in the context of background signals in the crystal environment and of implementation methods. We show that inverse electronic bridge processes quenching the isomeric state population can improve the performance of a solid-state nuclear clock based on $^{229m}$Th.
△ Less
Submitted 19 March, 2021;
originally announced March 2021.
-
Focal-point approach with pair-specific cusp correction for coupled-cluster theory
Authors:
Andreas Irmler,
Alejandro Gallo,
Andreas Grüneis
Abstract:
We present a basis set correction scheme for the coupled-cluster singles and doubles (CCSD) method. The scheme is based on employing frozen natural orbitals (FNOs) and diagrammatically decomposed contributions to the electronic correlation energy that dominate the basis set incompleteness error (BSIE). As recently discussed in [https://doi.org/10.1103/PhysRevLett.123.156401], the BSIE of the CCSD…
▽ More
We present a basis set correction scheme for the coupled-cluster singles and doubles (CCSD) method. The scheme is based on employing frozen natural orbitals (FNOs) and diagrammatically decomposed contributions to the electronic correlation energy that dominate the basis set incompleteness error (BSIE). As recently discussed in [https://doi.org/10.1103/PhysRevLett.123.156401], the BSIE of the CCSD correlation energy is dominated by the second-order Møller-Plesset (MP2) perturbation energy and the particle-particle ladder term. Here, we derive a simple approximation to the BSIE of the particle-particle ladder term that effectively corresponds to a rescaled pair-specific MP2 BSIE, where the scaling factor depends on the spatially averaged correlation hole depth of the coupled-cluster and first-order pair wavefunctions. The evaluation of the derived expressions is simple to implement in any existing code. We demonstrate the effectiveness of the method for the uniform electron gas. Furthermore, we apply the method to coupled-cluster theory calculations of atoms and molecules using FNOs. Employing the proposed correction and an increasing number of FNOs per occupied orbital, we demonstrate for a test set that rapidly convergent closed and open-shell reaction energies, atomization energies, electron affinities, and ionization potentials can be obtained. Moreover, we show that a similarly excellent trade-off between required virtual orbital basis set size and remaining BSIEs can be achieved for the perturbative triples contribution to the CCSD(T) energy employing FNOs and the (T*) approximation.
△ Less
Submitted 11 March, 2021;
originally announced March 2021.
-
Effective Hamiltonians for the study of real metals using quantum chemical theories
Authors:
Tina N. Mihm,
Tobias Schäfer,
Sai Kumar Ramadugu,
Andreas Grüneis,
James J. Shepherd
Abstract:
Computationally efficient and accurate quantum mechanical approximations to solve the many-electron Schrödinger equation are at the heart of computational materials science. In that respect the coupled cluster hierarchy of methods plays a central role in molecular quantum chemistry because of its systematic improvability and computational efficiency. In this hierarchy, coupled cluster singles and…
▽ More
Computationally efficient and accurate quantum mechanical approximations to solve the many-electron Schrödinger equation are at the heart of computational materials science. In that respect the coupled cluster hierarchy of methods plays a central role in molecular quantum chemistry because of its systematic improvability and computational efficiency. In this hierarchy, coupled cluster singles and doubles (CCSD) is one of the most important steps in moving towards chemical accuracy and, in recent years, its scope has successfully been expanded to the study of insulating surfaces and solids. Here, we show that CCSD theory can also be applied to real metals. In so doing, we overcome the limitation of needing extremely large supercells to capture long range electronic correlation effects. An effective Hamiltonian can be found using the transition structure factor--a map of electronic excitations from the Hartree--Fock wavefunction--which has fewer finite size effects than conventional periodic boundary conditions. This not only paves the way of applying coupled cluster methods to real metals but also reduces the computational cost by two orders of magnitude compared to previous methods. Our applications to phases of lithium and silicon show a resounding success in reaching the thermodynamic limit, taking the first step towards a truly universal quantum chemical treatment of solids.
△ Less
Submitted 10 March, 2021;
originally announced March 2021.
-
Local embedding of Coupled Cluster theory into the Random Phase Approximation using plane-waves
Authors:
Tobias Schäfer,
Florian Libisch,
Georg Kresse,
Andreas Grüneis
Abstract:
We present an embedding approach to treat local electron correlation effects in periodic environments. In a single, consistent framework, our plane-wave based scheme embeds a local high-level correlation calculation (here Coupled Cluster Theory, CC), employing localized orbitals, into a low-level correlation calculation (here the direct Random Phase Approximation, RPA). This choice allows for an a…
▽ More
We present an embedding approach to treat local electron correlation effects in periodic environments. In a single, consistent framework, our plane-wave based scheme embeds a local high-level correlation calculation (here Coupled Cluster Theory, CC), employing localized orbitals, into a low-level correlation calculation (here the direct Random Phase Approximation, RPA). This choice allows for an accurate and effcient treatment of long-range dispersion effects. Accelerated convergence with respect to the local fragment size can be observed if the low-level and high-level long-range dispersion are quantitatively similar, as is the case for CC in RPA. To demonstrate the capabilities of the introduced embedding approach, we calculate adsorption energies of molecules on a surface and in a chabazite crystal cage, as well as the formation energy of a lattice impurity in a solid at the level of highly accurate many-electron perturbation theories. The absorption energy of a methane molecule in a zeolite chabazite, for instance, is converged with an error well below 20 meV at the CC level. As our largest periodic benchmark system, we apply our scheme to the adsorption of a water molecule on titania in a supercell containing more than 1000 electrons.
△ Less
Submitted 11 December, 2020;
originally announced December 2020.
-
Structural and electronic properties of solid molecular hydrogen from many-electron theories
Authors:
Ke Liao,
Tong Shen,
Xin-Zheng Li,
Ali Alavi,
Andreas Grüneis
Abstract:
We study the structural and electronic properties of phase III of solid hydrogen using accurate many-electron theories and compare to state-of-the-art experimental findings. The atomic structures of phase III modelled by C2/c-24 crystals are fully optimized on the level of second-order perturbation theory, demonstrating that previously employed structures optimized on the level of approximate dens…
▽ More
We study the structural and electronic properties of phase III of solid hydrogen using accurate many-electron theories and compare to state-of-the-art experimental findings. The atomic structures of phase III modelled by C2/c-24 crystals are fully optimized on the level of second-order perturbation theory, demonstrating that previously employed structures optimized on the level of approximate density functionals exhibit errors in the H$_2$ bond lengths that cause significant discrepancies in the computed quasi particle band gaps and vibrational frequencies compared to experiment. Using the newly optimized atomic structures, we study the band gap closure and change in vibrational frequencies as a function of pressure. Our findings are in good agreement with recent experimental observations and may prove useful in resolving long-standing discrepancies between experimental estimates of metallization pressures possibly caused by disagreeing pressure calibrations.
△ Less
Submitted 18 November, 2020;
originally announced November 2020.
-
A periodic equation-of-motion coupled-cluster implementation applied to $F$-centers in alkaline earth oxides
Authors:
Alejandro Gallo,
Felix Hummel,
Andreas Irmler,
Andreas Grüneis
Abstract:
We present an implementation of equation of motion coupled-cluster singles and doubles (EOM-CCSD) theory using periodic boundary conditions and a plane wave basis set. Our implementation of EOM-CCSD theory is applied to study $F$-centers in alkaline earth oxides employing a periodic supercell approach. The convergence of calculated electronic excitation energies for neutral color centers in MgO, C…
▽ More
We present an implementation of equation of motion coupled-cluster singles and doubles (EOM-CCSD) theory using periodic boundary conditions and a plane wave basis set. Our implementation of EOM-CCSD theory is applied to study $F$-centers in alkaline earth oxides employing a periodic supercell approach. The convergence of calculated electronic excitation energies for neutral color centers in MgO, CaO and SrO crystals with respect to orbital basis set and system size is explored. We discuss extrapolation techniques that approximate excitation energies in the complete basis set limit and reduce finite size errors. Our findings demonstrate that EOM-CCSD theory can predict optical absorption energies of $F$-centers in good agreement with experiment. Furthermore, we discuss calculated emission energies corresponding to the decay from triplet to singlet states, responsible for the photoluminescence properties. Our findings are compared to experimental and theoretical results available in literature.
△ Less
Submitted 28 October, 2020;
originally announced October 2020.
-
Applying the Coupled-Cluster Ansatz to Solids and Surfaces in the Thermodynamic Limit
Authors:
Thomas Gruber,
Ke Liao,
Theodoros Tsatsoulis,
Felix Hummel,
Andreas Grüneis
Abstract:
Modern electronic structure theories can predict and simulate a wealth of phenomena in surface science and solid-state physics. In order to allow for a direct comparison with experiment, such ab initio predictions have to be made in the thermodynamic limit, substantially increasing the computational cost of many-electron wave-function theories. Here, we present a method that achieves thermodynamic…
▽ More
Modern electronic structure theories can predict and simulate a wealth of phenomena in surface science and solid-state physics. In order to allow for a direct comparison with experiment, such ab initio predictions have to be made in the thermodynamic limit, substantially increasing the computational cost of many-electron wave-function theories. Here, we present a method that achieves thermodynamic limit results for solids and surfaces using the "gold standard" coupled cluster ansatz of quantum chemistry with unprecedented efficiency. We study the energy difference between carbon diamond and graphite crystals, adsorption energies of water on h-BN, as well as the cohesive energy of the Ne solid, demonstrating the increased efficiency and accuracy of coupled cluster theory for solids and surfaces.
△ Less
Submitted 14 April, 2020;
originally announced April 2020.
-
Coupled cluster theory in materials science
Authors:
Igor Ying Zhang,
Andreas Grüneis
Abstract:
In this tutorial-style review we discuss basic concepts of coupled cluster theory and recent developments that increase its computational efficiency for calculations of molecules, solids and materials in general. We will touch upon the connection between coupled cluster theory and the random-phase approximation that is widely used in the field of solid-state physics. We will discuss various approa…
▽ More
In this tutorial-style review we discuss basic concepts of coupled cluster theory and recent developments that increase its computational efficiency for calculations of molecules, solids and materials in general. We will touch upon the connection between coupled cluster theory and the random-phase approximation that is widely used in the field of solid-state physics. We will discuss various approaches to improve the computational performance without compromising on accuracy. These approaches include large-scale parallel design as well as techniques that reduce the pre-factor of the computational complexity. A central part of this article discusses the convergence of calculated properties to the thermodynamic limit, which is of significant importance for reliable predictions of materials properties and constitutes an additional challenge compared to calculations of large molecules. We mention technical aspects of computer code implementations of periodic coupled cluster theories in different numerical frameworks of the one-electron orbital basis; the projector-augmented-wave formalism using a plane wave basis set and the numeric atom-centered-orbital (NAO) with resolution-of-identity. We will discuss results and the possible scope of these implementations and how they can help advance the current state of the art in electronic structure theory calculations of materials.
△ Less
Submitted 14 April, 2020;
originally announced April 2020.
-
Particle-particle ladder based basis-set corrections applied to atoms and molecules using coupled-cluster theory
Authors:
Andreas Irmler,
Andreas Grüneis
Abstract:
We investigate the basis-set convergence of electronic correlation energies calculated using coupled cluster theory and a recently proposed finite basis-set correction technique. The correction is applied to atomic and molecular systems and is based on a diagrammatically decomposed coupled cluster singles and doubles correlation energy. Only the second-order energy and the particle-particle ladder…
▽ More
We investigate the basis-set convergence of electronic correlation energies calculated using coupled cluster theory and a recently proposed finite basis-set correction technique. The correction is applied to atomic and molecular systems and is based on a diagrammatically decomposed coupled cluster singles and doubles correlation energy. Only the second-order energy and the particle-particle ladder term are corrected for their basis-set incompleteness error. We present absolute correlation energies and results for a large benchmark set. Our findings indicate that basis set reductions by two cardinal numbers are possible for atomization energies, ionization potentials and electron affinities without compromising accuracy when compared to conventional CCSD calculations. In the case of reaction energies we find that reductions by one cardinal number are possible compared to conventional CCSD calculations. The employed technique can readily be applied to other many-electron theories without the need for three- or four-electron integrals.
△ Less
Submitted 2 October, 2019; v1 submitted 13 March, 2019;
originally announced March 2019.
-
On the duality of ring and ladder diagrams and its importance for many-electron perturbation theories
Authors:
Andreas Irmler,
Felix Hummel,
Andreas Grüneis
Abstract:
We present a diagrammatic decomposition of the transition pair correlation function for the uniform electron gas. We demonstrate explicitly that ring and ladder diagrams are dual counterparts that capture significant long- and short-ranged interelectronic correlation effects, respectively. Our findings help to guide the further development of approximate many-electron theories and reveal that the…
▽ More
We present a diagrammatic decomposition of the transition pair correlation function for the uniform electron gas. We demonstrate explicitly that ring and ladder diagrams are dual counterparts that capture significant long- and short-ranged interelectronic correlation effects, respectively. Our findings help to guide the further development of approximate many-electron theories and reveal that the contribution of the ladder diagrams to the electronic correlation energy can be approximated in an effective manner using second-order perturbation theory. We employ the latter approximation to reduce the computational cost of coupled cluster theory calculations for insulators and semiconductors by two orders of magnitude without compromising accuracy.
△ Less
Submitted 14 March, 2019; v1 submitted 13 March, 2019;
originally announced March 2019.
-
On the physisorption of water on graphene: Sub-chemical accuracy from many-body electronic structure methods
Authors:
Jan Gerit Brandenburg,
Andrea Zen,
Martin Fitzner,
Benjamin Ramberger,
Georg Kresse,
Theodoros Tsatsoulis,
Andreas Grüneis,
Angelos Michaelides,
Dario Alfè
Abstract:
Molecular adsorption on surfaces plays a central role in catalysis, corrosion, desalination, and many other processes of relevance to industry and the natural world. Few adsorption systems are more ubiquitous or of more widespread importance than those involving water and carbon, and for a molecular level understanding of such interfaces water monomer adsorption on graphene is a fundamental and re…
▽ More
Molecular adsorption on surfaces plays a central role in catalysis, corrosion, desalination, and many other processes of relevance to industry and the natural world. Few adsorption systems are more ubiquitous or of more widespread importance than those involving water and carbon, and for a molecular level understanding of such interfaces water monomer adsorption on graphene is a fundamental and representative system. This system is particularly interesting as it calls for an accurate treatment of electron correlation effects, as well as posing a practical challenge to experiments. Here, we employ many-body electronic structure methodologies that can be rigorously converged and thus provide faithful references for the molecule-surface interaction. In particular, we use diffusion Monte-Carlo (DMC), coupled cluster (CCSD(T)), as well as the random phase approximation (RPA) to calculate the strength of the interaction between water and an extended graphene surface. We establish excellent, sub-chemical, agreement between the complementary high-level methodologies, and an adsorption energy estimate in the most stable configuration of approximately -100\,meV is obtained. We also find that the adsorption energy is rather insensitive to the orientation of the water molecule on the surface, despite different binding motifs involving qualitatively different interfacial charge reorganisation. In producing the first demonstrably accurate adsorption energies for water on graphene this work also resolves discrepancies amongst previously reported values for this widely studied system. It also paves the way for more accurate and reliable studies of liquid water at carbon interfaces with cheaper computational methods, such as density functional theory and classical potentials.
△ Less
Submitted 19 November, 2018;
originally announced November 2018.
-
Reaction energetics of Hydrogen on the Si(100) surface: A periodic many-electron theory study
Authors:
Theodoros Tsatsoulis,
Sung Sakong,
Axel Groß,
Andreas Grüneis
Abstract:
We report on a many-electron wavefunction theory study for the reaction energetics of hydrogen dissociation on the Si(100) surface. We demonstrate that quantum chemical wavefunction based methods using periodic boundary conditions can predict chemically accurate results for the activation barrier and the chemisorption energy in agreement with experimental findings. These highly accurate results fo…
▽ More
We report on a many-electron wavefunction theory study for the reaction energetics of hydrogen dissociation on the Si(100) surface. We demonstrate that quantum chemical wavefunction based methods using periodic boundary conditions can predict chemically accurate results for the activation barrier and the chemisorption energy in agreement with experimental findings. These highly accurate results for the reaction energetics enable a deeper understanding of the underlying physical mechanism and make it possible to benchmark widely used density functional theory methods.
△ Less
Submitted 18 September, 2018; v1 submitted 21 August, 2018;
originally announced August 2018.
-
Plasmonic enhancement of aligned semiconducting graphene nanoribbons
Authors:
M. Pfeiffer,
B. V. Senkovskiy,
D. Haberer,
F. R. Fischer,
F. Yang,
K. Meerholz,
Y. Ando,
A. Grüneis,
K. Lindfors
Abstract:
We couple photoluminescent semiconducting 7-atom wide armchair edge graphene nanoribbons to plasmonic nanoantenna arrays and demonstrate an enhancement of the photoluminescence and Raman scattering intensity of the nanoribbons by more than one order of magnitude. The increase in signal allows us to study Raman spectra with high signal-to-noise ratio. Using plasmonic enhancement we are able to dete…
▽ More
We couple photoluminescent semiconducting 7-atom wide armchair edge graphene nanoribbons to plasmonic nanoantenna arrays and demonstrate an enhancement of the photoluminescence and Raman scattering intensity of the nanoribbons by more than one order of magnitude. The increase in signal allows us to study Raman spectra with high signal-to-noise ratio. Using plasmonic enhancement we are able to detect the off-resonant Raman signals from the modified radial breathing-like mode (RBLM) due to physisorbed molecules, the 3rd order RBLM, and C-H vibrations. We find excellent agreement between data and simulations describing the spectral dependence of the enhancement and modifications of the polarization anisotropy. The strong field gradients in the optical near-field further allow us to probe the subwavelength coherence properties of the phonon modes in the nanoribbons. We theoretically model this considering a finite coherence length along the GNR direction. Our results allow estimating the coherence length in graphene nanoribbons.
△ Less
Submitted 21 December, 2017;
originally announced December 2017.
-
A comparison between quantum chemistry and quantum Monte Carlo techniques for the adsorption of water on the (001) LiH surface
Authors:
Theodoros Tsatsoulis,
Felix Hummel,
Denis Usvyat,
Martin Schütz,
George H. Booth,
Simon S. Binnie,
Michael J. Gillan,
Dario Alfè,
Angelos Michaelides,
Andreas Grüneis
Abstract:
We present a comprehensive benchmark study of the adsorption energy of a single water molecule on the (001) LiH surface using periodic coupled cluster and quantum Monte Carlo theories. We benchmark and compare different implementations of quantum chemical wave function based theories in order to verify the reliability of the predicted adsorption energies and the employed approximations. Furthermor…
▽ More
We present a comprehensive benchmark study of the adsorption energy of a single water molecule on the (001) LiH surface using periodic coupled cluster and quantum Monte Carlo theories. We benchmark and compare different implementations of quantum chemical wave function based theories in order to verify the reliability of the predicted adsorption energies and the employed approximations. Furthermore we compare the predicted adsorption energies to those obtained employing widely-used van der Waals density-functionals. Our findings show that quantum chemical approaches are becoming a robust and reliable tool for condensed phase electronic structure calculations, providing an additional tool that can also help in potentially improving currently available van der Waals density-functionals.
△ Less
Submitted 9 May, 2017; v1 submitted 21 February, 2017;
originally announced February 2017.
-
Low rank factorization of the Coulomb integrals for periodic coupled cluster theory
Authors:
Felix Hummel,
Theodoros Tsatsoulis,
Andreas Grüneis
Abstract:
We study the decomposition of the Coulomb integrals of periodic systems into a tensor contraction of six matrices of which only two are distinct. We find that the Coulomb integrals can be well approximated in this form already with small matrices compared to the number of real space grid points. The cost of computing the matrices scales as O(N^4) using a regularized form of the alternating least s…
▽ More
We study the decomposition of the Coulomb integrals of periodic systems into a tensor contraction of six matrices of which only two are distinct. We find that the Coulomb integrals can be well approximated in this form already with small matrices compared to the number of real space grid points. The cost of computing the matrices scales as O(N^4) using a regularized form of the alternating least squares algorithm. The studied factorization of the Coulomb integrals can be exploited to reduce the scaling of the computational cost of expensive tensor contractions appearing in the amplitude equations of coupled cluster methods with respect to system size. We apply the developed methodologies to calculate the adsorption energy of a single water molecule on a hexagonal boron nitride monolayer in a plane wave basis set and periodic boundary conditions.
△ Less
Submitted 10 November, 2016;
originally announced November 2016.
-
From plane waves to local Gaussians for the simulation of correlated periodic systems
Authors:
George H. Booth,
Theodoros Tsatsoulis,
Garnet Kin-Lic Chan,
Andreas Grüneis
Abstract:
We present a simple, robust and black-box approach to the implementation and use of local, periodic, atom-centered Gaussian basis functions within a plane wave code, in a computationally efficient manner. The procedure outlined is based on the representation of the Gaussians within a finite bandwidth by their underlying plane wave coefficients. The core region is handled within the projected augme…
▽ More
We present a simple, robust and black-box approach to the implementation and use of local, periodic, atom-centered Gaussian basis functions within a plane wave code, in a computationally efficient manner. The procedure outlined is based on the representation of the Gaussians within a finite bandwidth by their underlying plane wave coefficients. The core region is handled within the projected augment wave framework, by pseudizing the Gaussian functions within a cut-off radius around each nucleus, smoothing the functions so that they are faithfully represented by a plane wave basis with only moderate kinetic energy cutoff. To mitigate the effects of the basis set superposition error and incompleteness at the mean-field level introduced by the Gaussian basis, we also propose a hybrid approach, whereby the complete occupied space is first converged within a large plane wave basis, and the Gaussian basis used to construct a complementary virtual space for the application of correlated methods. We demonstrate that these pseudized Gaussians yield compact and systematically improvable spaces with an accuracy comparable to their non-pseudized Gaussian counterparts. A key advantage of the described method is its ability to efficiently capture and describe electronic correlation effects of weakly bound and low-dimensional systems, where plane waves are not sufficiently compact or able to be truncated without unphysical artefacts. We investigate the accuracy of the pseudized Gaussians for the water dimer interaction, neon solid and water adsorption on a LiH surface, at the level of second-order M\{o}ller--Plesset perturbation theory.
△ Less
Submitted 18 August, 2016; v1 submitted 21 March, 2016;
originally announced March 2016.
-
Many-body quantum chemistry for the electron gas: convergent perturbative theories
Authors:
James J. Shepherd,
Andreas Grüneis
Abstract:
We investigate the accuracy of a number of wavefunction based methods at the heart of quantum chemistry for metallic systems. Using Hartree-Fock as a reference, perturbative (Møller-Plesset, MP) and coupled cluster (CC) theories are used to study the uniform electron gas model. Our findings suggest that non-perturbative coupled cluster theories are acceptable for modelling electronic interactions…
▽ More
We investigate the accuracy of a number of wavefunction based methods at the heart of quantum chemistry for metallic systems. Using Hartree-Fock as a reference, perturbative (Møller-Plesset, MP) and coupled cluster (CC) theories are used to study the uniform electron gas model. Our findings suggest that non-perturbative coupled cluster theories are acceptable for modelling electronic interactions in metals whilst perturbative coupled cluster theories are not. Using screened interactions, we propose a simple modification to the widely-used coupled-cluster singles and doubles plus perturbative triples method (CCSD(T)) that lifts the divergent behaviour and is shown to give very accurate correlation energies for the homogeneous electron gas.
△ Less
Submitted 22 October, 2013;
originally announced October 2013.
-
Explicitly correlated plane waves: Accelerating convergence in periodic wavefunction expansions
Authors:
Andreas Grüneis,
James J. Shepherd,
Ali Alavi,
David P. Tew,
George H. Booth
Abstract:
We present an investigation into the use of an explicitly correlated plane wave basis for periodic wavefunction expansions at the level of second-order Møller-Plesset perturbation theory (MP2). The convergence of the electronic correlation energy with respect to the one-electron basis set is investigated and compared to conventional MP2 theory in a finite homogeneous electron gas model. In additio…
▽ More
We present an investigation into the use of an explicitly correlated plane wave basis for periodic wavefunction expansions at the level of second-order Møller-Plesset perturbation theory (MP2). The convergence of the electronic correlation energy with respect to the one-electron basis set is investigated and compared to conventional MP2 theory in a finite homogeneous electron gas model. In addition to the widely used Slater-type geminal correlation factor, we also derive and investigate a novel correlation factor that we term Yukawa-Coulomb. The Yukawa-Coulomb correlation factor is motivated by analytic results for two electrons in a box and allows for a further improved convergence of the correlation energies with respect to the employed basis set. We find the combination of the infinitely delocalized plane waves and local short-ranged geminals provides a complementary, and rapidly convergent basis for the description of periodic wavefunctions. We hope that this approach will expand the scope of discrete wavefunction expansions in periodic systems.
△ Less
Submitted 24 July, 2013;
originally announced July 2013.
-
Correlation Energy Divergences in Metallic Systems
Authors:
James J. Shepherd,
Andreas Grüneis
Abstract:
We numerically examine divergences of the total energy in metallic systems of approximate many-body theories using Hartree--Fock as a reference, including perturbative (M\oller-Plesset, MP), coupled cluster (CC) and configuration interaction (CI) approaches. Controlling for finite size effects and basis set incompleteness error by comparison with energies from the random phase approximation (RPA),…
▽ More
We numerically examine divergences of the total energy in metallic systems of approximate many-body theories using Hartree--Fock as a reference, including perturbative (M\oller-Plesset, MP), coupled cluster (CC) and configuration interaction (CI) approaches. Controlling for finite size effects and basis set incompleteness error by comparison with energies from the random phase approximation (RPA), we suggest convincingly that non-perturbative coupled cluster theories are acceptable for modelling electronic interactions in metals whilst perturbative coupled cluster theories are not. Data are provided from the RPA with which it is possible to test other approximate finite-basis methods for divergences with only modest computational cost.
△ Less
Submitted 30 August, 2012;
originally announced August 2012.
-
Convergence of many-body wavefunction expansions using a plane wave basis: from the homogeneous electron gas to the solid state
Authors:
James J. Shepherd,
Andreas Grüneis,
George H. Booth,
Georg Kresse,
Ali Alavi
Abstract:
Using the finite simulation-cell homogeneous electron gas (HEG) as a model, we investigate the convergence of the correlation energy to the complete basis set (CBS) limit in methods utilising plane-wave wavefunction expansions. Simple analytic and numerical results from second-order Møller-Plesset theory (MP2) suggest a 1/M decay of the basis-set incompleteness error where M is the number of plane…
▽ More
Using the finite simulation-cell homogeneous electron gas (HEG) as a model, we investigate the convergence of the correlation energy to the complete basis set (CBS) limit in methods utilising plane-wave wavefunction expansions. Simple analytic and numerical results from second-order Møller-Plesset theory (MP2) suggest a 1/M decay of the basis-set incompleteness error where M is the number of plane waves used in the calculation, allowing for straightforward extrapolation to the CBS limit. As we shall show, the choice of basis set truncation when constructing many-electron wavefunctions is far from obvious, and here we propose several alternatives based on the momentum transfer vector, which greatly improve the rate of convergence. This is demonstrated for a variety of wavefunction methods, from MP2 to coupled-cluster doubles theory (CCD) and the random-phase approximation plus second-order screened exchange (RPA+SOSEX). Finite basis-set energies are presented for these methods and compared with exact benchmarks. A transformation can map the orbitals of a general solid state system onto the HEG plane wave basis and thereby allow application of these methods to more realistic physical problems.
△ Less
Submitted 9 May, 2012; v1 submitted 22 February, 2012;
originally announced February 2012.
-
A Full Configuration Interaction Perspective on the Homogeneous Electron Gas
Authors:
James J. Shepherd,
George Booth,
Andreas Grüneis,
Ali Alavi
Abstract:
Highly accurate results for the homogeneous electron gas (HEG) have only been achieved to date within a diffusion Monte Carlo (DMC) framework. Here, we introduce a newly developed stochastic technique, Full Configuration Interaction Quantum Monte Carlo (FCIQMC), which samples the exact wavefunction expanded in plane wave Slater determinants. Despite the introduction of a basis set incompleteness e…
▽ More
Highly accurate results for the homogeneous electron gas (HEG) have only been achieved to date within a diffusion Monte Carlo (DMC) framework. Here, we introduce a newly developed stochastic technique, Full Configuration Interaction Quantum Monte Carlo (FCIQMC), which samples the exact wavefunction expanded in plane wave Slater determinants. Despite the introduction of a basis set incompleteness error, we obtain a finite-basis energy which is significantly, and variationally lower than any previously published work for the 54-electron HEG at $r_s$ = 0.5 a.u., in a Hilbert space of $10^{108}$ Slater determinants. At this value of $r_s$, as well as of 1.0 a.u., we remove the remaining basis set incompleteness error by extrapolation, yielding results comparable or better than state-of-the-art DMC backflow energies. In doing so, we demonstrate that it is possible to yield highly accurate results with the FCIQMC method in sizable periodic systems.
△ Less
Submitted 12 September, 2011;
originally announced September 2011.