1
|
Abraham V, Harsha G, Zgid D. Relativistic Fully Self-Consistent GW for Molecules: Total Energies and Ionization Potentials. J Chem Theory Comput 2024; 20:4579-4590. [PMID: 38778459 DOI: 10.1021/acs.jctc.4c00075] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2024]
Abstract
The fully self-consistent GW (scGW) method with an iterative solution of the Dyson equation provides a consistent approach for describing the ground and excited states without any dependence on the mean-field reference. In this work, we present a relativistic version of scGW for molecules containing heavy elements using the exact two-component (X2C) Coulomb approximation. We benchmark SOC-81 data set containing closed shell heavy elements for the first ionization potential using the fully self-consistent GW as well as one-shot GW. The self-consistent GW provides superior results compared to G0W0 with PBE reference and comparable results to G0W0 with PBE0 while also removing the starting point dependence. The photoelectron spectra obtained at the X2C level demonstrate very good agreement with the experimental spectra. We also observe that scGW provides very good estimation of ionization potential for the inner d-shell orbitals. Additionally, using the well-conserved total energy, we investigate the equilibrium bond length and harmonic frequencies of a few halogen dimers using scGW. Overall, our findings demonstrate the applicability of the fully self-consistent GW method for accurate ionization potential, photoelectron spectra, and total energies in finite systems with heavy elements with a reasonable computational scaling.
Collapse
Affiliation(s)
- Vibin Abraham
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Gaurav Harsha
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Dominika Zgid
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
- Department of Physics and Astronomy, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
2
|
Schäfer T, Van Benschoten WZ, Shepherd JJ, Grüneis A. Sampling the reciprocal Coulomb potential in finite anisotropic cells. J Chem Phys 2024; 160:051101. [PMID: 38310470 DOI: 10.1063/5.0182729] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2023] [Accepted: 01/07/2024] [Indexed: 02/05/2024] Open
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 and discretization 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 and coupled cluster (CC) theories and discuss the consequences of different sampling strategies on 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, and 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.
Collapse
Affiliation(s)
- Tobias Schäfer
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria
| | | | - James J Shepherd
- Department of Chemistry, University of Iowa, Iowa City, Iowa 52242, USA
| | - Andreas Grüneis
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria
| |
Collapse
|
3
|
Monzel L, Holzer C, Klopper W. Natural virtual orbitals for the GW method in the random-phase approximation and beyond. J Chem Phys 2023; 158:144102. [PMID: 37061489 DOI: 10.1063/5.0144469] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2023] Open
Abstract
The increasingly popular GW method is becoming a convenient tool to determine vertical ionization energies in molecular systems. However, depending on the formalism used and the range of orbitals investigated, it may be hampered by a steep computational scaling. To alleviate this issue, correlated natural virtual orbitals (NVOs) based on second-order Møller-Plesset (MP2) and direct MP2 correlation energies are implemented, and the resulting correlated NVOs are tested on GW quasiparticle energies. Test cases include the popular GW variants G0W0 and evGW0 as well as more elaborate vertex corrections. We find that for increasingly larger molecular systems and basis sets, NVOs considerably improve efficiency. Furthermore, we test the performance of the truncated (frozen) NVO ansatz on the GW100 test set. For the latter, it is demonstrated that, using a carefully chosen truncation threshold, NVOs lead to a negligible loss in accuracy while providing speedups of one order of magnitude. Furthermore, we compare the resulting quasiparticle energies to very accurate vertical ionization energies obtained from coupled-cluster theory with singles, doubles, and noniterative triples [CCSD(T)], confirming that the loss in accuracy introduced by truncating the NVOs is negligible compared to the methodical errors in the GW approximation. It is also demonstrated that the choice of basis set impacts results far more than using a suitably truncated NVO space. Therefore, at the same computational expense, more accurate results can be obtained using NVOs. Finally, we provide improved reference CCSD(T) values for the GW100 test set, which have been obtained using the def2-QZVPP basis set.
Collapse
Affiliation(s)
- Laurenz Monzel
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany
| | - Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), Wolfgang-Gaede-Straße 1, 76131 Karlsruhe, Germany
| | - Wim Klopper
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany
| |
Collapse
|
4
|
Tiny (ZnO) clusters supported on graphene for solar energy trapping: A density functional theory study. J Taiwan Inst Chem Eng 2023. [DOI: 10.1016/j.jtice.2023.104769] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/27/2023]
|
5
|
Kállay M, Horváth RA, Gyevi-Nagy L, Nagy PR. Basis Set Limit CCSD(T) Energies for Extended Molecules via a Reduced-Cost Explicitly Correlated Approach. J Chem Theory Comput 2022; 19:174-189. [PMID: 36576419 PMCID: PMC9835832 DOI: 10.1021/acs.jctc.2c01031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Several approximations are introduced and tested to reduce the computational expenses of the explicitly correlated coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method for both closed and open-shell species. First, the well-established frozen natural orbital (FNO) technique is adapted to explicitly correlated CC approaches. Second, our natural auxiliary function (NAF) scheme is employed to reduce the size of the auxiliary basis required for the density fitting approximation regularly used in explicitly correlated calculations. Third, a new approach, termed the natural auxiliary basis (NAB) approximation, is proposed to decrease the size of the auxiliary basis needed for the expansion of the explicitly correlated geminals. The performance of the above approximations and that of the combined FNO-NAF-NAB approach are tested for atomization and reaction energies. Our results show that overall speedups of 7-, 5-, and 3-times can be achieved with double-, triple-, and quadruple-ζ basis sets, respectively, without any loss in accuracy. The new method can provide, e.g., reaction energies and barrier heights well within chemical accuracy for molecules with more than 40 atoms within a few days using a few dozen processor cores, and calculations with 50+ atoms are still feasible. These routinely affordable computations considerably extend the reach of explicitly correlated CCSD(T).
Collapse
Affiliation(s)
- Mihály Kállay
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,
| | - Réka A. Horváth
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - László Gyevi-Nagy
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Péter R. Nagy
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| |
Collapse
|
6
|
Schäfer T, Gallo A, Irmler A, Hummel F, Grüneis A. Surface science using coupled cluster theory via local Wannier functions and in-RPA-embedding: The case of water on graphitic carbon nitride. J Chem Phys 2021; 155:244103. [PMID: 34972356 DOI: 10.1063/5.0074936] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A first-principles study of the adsorption of a single water molecule on a layer of graphitic carbon nitride is reported employing an embedding approach for many-electron correlation methods. To this end, a plane-wave based implementation to obtain intrinsic atomic orbitals and Wannier functions for arbitrary localization potentials is presented. In our embedding scheme, the localized occupied orbitals allow for a separate treatment of short-range and long-range correlation contributions to the adsorption energy by a fragmentation of the simulation cell. In combination with unoccupied natural orbitals, the coupled cluster ansatz with single, double, and perturbative triple particle-hole excitation operators is used to capture the correlation in local fragments centered around the adsorption process. For the long-range correlation, a seamless embedding into the random phase approximation yields rapidly convergent adsorption energies with respect to the local 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.
Collapse
Affiliation(s)
- Tobias Schäfer
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria
| | - Alejandro Gallo
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria
| | - Andreas Irmler
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria
| | - Felix Hummel
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria
| | - Andreas Grüneis
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria
| |
Collapse
|
7
|
Nagy PR, Gyevi-Nagy L, Kállay M. Basis set truncation corrections for improved frozen natural orbital CCSD(T) energies. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1963495] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Péter R. Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
| | - László Gyevi-Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
| |
Collapse
|
8
|
Lee J, Feng X, Cunha LA, Gonthier JF, Epifanovsky E, Head-Gordon M. Approaching the basis set limit in Gaussian-orbital-based periodic calculations with transferability: Performance of pure density functionals for simple semiconductors. J Chem Phys 2021; 155:164102. [PMID: 34717349 PMCID: PMC8556001 DOI: 10.1063/5.0069177] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2021] [Accepted: 10/05/2021] [Indexed: 12/14/2022] Open
Abstract
Simulating solids with quantum chemistry methods and Gaussian-type orbitals (GTOs) has been gaining popularity. Nonetheless, there are few systematic studies that assess the basis set incompleteness error (BSIE) in these GTO-based simulations over a variety of solids. In this work, we report a GTO-based implementation for solids and apply it to address the basis set convergence issue. We employ a simple strategy to generate large uncontracted (unc) GTO basis sets that we call the unc-def2-GTH sets. These basis sets exhibit systematic improvement toward the basis set limit as well as good transferability based on application to a total of 43 simple semiconductors. Most notably, we found the BSIE of unc-def2-QZVP-GTH to be smaller than 0.7 mEh per atom in total energies and 20 meV in bandgaps for all systems considered here. Using unc-def2-QZVP-GTH, we report bandgap benchmarks of a combinatorially designed meta-generalized gradient approximation (mGGA) functional, B97M-rV, and show that B97M-rV performs similarly (a root-mean-square-deviation of 1.18 eV) to other modern mGGA functionals, M06-L (1.26 eV), MN15-L (1.29 eV), and Strongly Constrained and Appropriately Normed (SCAN) (1.20 eV). This represents a clear improvement over older pure functionals such as local density approximation (1.71 eV) and Perdew-Burke-Ernzerhof (PBE) (1.49 eV), although all these mGGAs are still far from being quantitatively accurate. We also provide several cautionary notes on the use of our uncontracted bases and on future research on GTO basis set development for solids.
Collapse
Affiliation(s)
- Joonho Lee
- Department of Chemistry, Columbia University, New York, New York 10027,, USA
| | | | - Leonardo A. Cunha
- Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, USA
| | - Jérôme F. Gonthier
- Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, USA
| | | | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, USA
| |
Collapse
|
9
|
Choudhary K. Quantum computation for predicting electron and phonon properties of solids. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 33:385501. [PMID: 34225258 DOI: 10.1088/1361-648x/ac1154] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/26/2021] [Accepted: 07/05/2021] [Indexed: 06/13/2023]
Abstract
Quantum chemistry is one of the most promising near-term applications of quantum computers. Quantum algorithms such as variational quantum eigen solver (VQE) and variational quantum deflation (VQD) algorithms have been mainly applied for molecular systems and there is a need to implement such methods for periodic solids. Using Wannier tight-binding Hamiltonian (WTBH) approaches, we demonstrate the application of VQE and VQD to accurately predict both electronic and phonon bandstructure properties of several elemental as well as multi-component solid-state materials. We apply VQE-VQD calculations for 307 spin-orbit coupling based electronic WTBHs and 933 finite-difference based phonon WTBHs. Also, we discuss a workflow for using VQD with lattice Green's function that can be used for solving dynamical mean-field theory problems. The WTBH model solvers can be used for testing other quantum algorithms and models also.
Collapse
Affiliation(s)
- Kamal Choudhary
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, MD 20899, United States of America
- Theiss Research, La Jolla, CA 92037, United States of America
| |
Collapse
|
10
|
Riemelmoser S, Kaltak M, Kresse G. Optimized effective potentials from the random-phase approximation: Accuracy of the quasiparticle approximation. J Chem Phys 2021; 154:154103. [PMID: 33887939 DOI: 10.1063/5.0045400] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
The optimized effective potential (OEP) method presents an unambiguous way to construct the Kohn-Sham potential corresponding to a given diagrammatic approximation for the exchange-correlation functional. The OEP from the random-phase approximation (RPA) has played an important role ever since the conception of the OEP formalism. However, the solution of the OEP equation is computationally fairly expensive and has to be done in a self-consistent way. So far, large scale solid state applications have, therefore, been performed only using the quasiparticle approximation (QPA), neglecting certain dynamical screening effects. We obtain the exact RPA-OEP for 15 semiconductors and insulators by direct solution of the linearized Sham-Schlüter equation. We investigate the accuracy of the QPA on Kohn-Sham bandgaps and dielectric constants, and comment on the issue of self-consistency.
Collapse
Affiliation(s)
- Stefan Riemelmoser
- Faculty of Physics and Center for Computational Materials Science, University of Vienna, Kolingasse 14-16, A-1090 Vienna, Austria
| | - Merzuk Kaltak
- VASP Software GmbH, Sensengasse 8/17, A-1090 Vienna, Austria
| | - Georg Kresse
- Faculty of Physics and Center for Computational Materials Science, University of Vienna, Kolingasse 14-16, A-1090 Vienna, Austria
| |
Collapse
|
11
|
Eidsvåg H, Bentouba S, Vajeeston P, Yohi S, Velauthapillai D. TiO 2 as a Photocatalyst for Water Splitting-An Experimental and Theoretical Review. Molecules 2021; 26:molecules26061687. [PMID: 33802911 PMCID: PMC8002707 DOI: 10.3390/molecules26061687] [Citation(s) in RCA: 44] [Impact Index Per Article: 14.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2021] [Revised: 03/09/2021] [Accepted: 03/10/2021] [Indexed: 11/16/2022] Open
Abstract
Hydrogen produced from water using photocatalysts driven by sunlight is a sustainable way to overcome the intermittency issues of solar power and provide a green alternative to fossil fuels. TiO2 has been used as a photocatalyst since the 1970s due to its low cost, earth abundance, and stability. There has been a wide range of research activities in order to enhance the use of TiO2 as a photocatalyst using dopants, modifying the surface, or depositing noble metals. However, the issues such as wide bandgap, high electron-hole recombination time, and a large overpotential for the hydrogen evolution reaction (HER) persist as a challenge. Here, we review state-of-the-art experimental and theoretical research on TiO2 based photocatalysts and identify challenges that have to be focused on to drive the field further. We conclude with a discussion of four challenges for TiO2 photocatalysts-non-standardized presentation of results, bandgap in the ultraviolet (UV) region, lack of collaboration between experimental and theoretical work, and lack of large/small scale production facilities. We also highlight the importance of combining computational modeling with experimental work to make further advances in this exciting field.
Collapse
Affiliation(s)
- Håkon Eidsvåg
- Department of Computing, Mathematics and Physics, Western Norway University of Applied Sciences, Inndalsveien 28, Box 5063, N-5009 Bergen, Norway;
- Correspondence: (H.E.); (D.V.); Tel.: +47-980-61-444 (H.E.); +47-55-58-77-11 (D.V.)
| | - Said Bentouba
- Department of Computing, Mathematics and Physics, Western Norway University of Applied Sciences, Inndalsveien 28, Box 5063, N-5009 Bergen, Norway;
| | - Ponniah Vajeeston
- Center for Materials Science and Nanotechnology, Department of Chemistry, University of Oslo, Box 1033 Blindern, N-0315 Oslo, Norway;
| | - Shivatharsiny Yohi
- Department of Chemistry, Faculty of Science, University of Jaffna, Sir. Pon, Ramanathan Rd, Jaffna 40000, Sri Lanka;
| | - Dhayalan Velauthapillai
- Department of Computing, Mathematics and Physics, Western Norway University of Applied Sciences, Inndalsveien 28, Box 5063, N-5009 Bergen, Norway;
- Correspondence: (H.E.); (D.V.); Tel.: +47-980-61-444 (H.E.); +47-55-58-77-11 (D.V.)
| |
Collapse
|
12
|
Ramberger B, Kresse G. New insights into the 1D carbon chain through the RPA. Phys Chem Chem Phys 2021; 23:5254-5260. [PMID: 33629671 DOI: 10.1039/d0cp06607a] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We investigated the electronic and structural properties of the infinite linear carbon chain (carbyne) using density functional theory (DFT) and the random phase approximation (RPA) to the correlation energy. The studies are performed in vacuo and for carbyne inside a carbon nano tube (CNT). In the vacuum, semi-local DFT and RPA predict bond length alternations of about 0.04 Å and 0.13 Å, respectively. The frequency of the highest optical mode at the Γ point is 1219 cm-1 and about 2000 cm-1 for DFT and the RPA. Agreement of the RPA to previous high level quantum chemistry and diffusion Monte-Carlo results is excellent. For the RPA we calculate the phonon-dispersion in the full Brillouine zone and find marked quantitative differences to DFT calculations not only at the Γ point but also throughout the entire Brillouine zone. To model carbyne inside a carbon nanotube, we considered a (10,0) CNT. Here the DFT calculations are even qualitatively sensitive to the k-points sampling. At the limes of a very dense k-points sampling, semi-local DFT predicts no bond length alternation (BLA), whereas in the RPA a sizeable BLA of 0.09 Å prevails. The reduced BLA leads to a significant red shift of the vibrational frequencies of about 350 cm-1, so that they are in good agreement with experimental estimates. Overall, the good agreement between the RPA and previously reported results from correlated wavefunction methods and experimental Raman data suggests that the RPA provides reliable results at moderate computational costs. It hence presents a useful addition to the repertoire of correlated wavefunction methods and its accuracy clearly prevails for low dimensional systems, where semi-local density functionals struggle to yield even qualitatively correct results.
Collapse
Affiliation(s)
- Benjamin Ramberger
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Kolingasse 14-16, 1090 Vienna, Austria.
| | - Georg Kresse
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Kolingasse 14-16, 1090 Vienna, Austria.
| |
Collapse
|
13
|
Gyevi-Nagy L, Kállay M, Nagy PR. Accurate Reduced-Cost CCSD(T) Energies: Parallel Implementation, Benchmarks, and Large-Scale Applications. J Chem Theory Comput 2021; 17:860-878. [PMID: 33400527 PMCID: PMC7884001 DOI: 10.1021/acs.jctc.0c01077] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2020] [Indexed: 11/28/2022]
Abstract
The accurate and systematically improvable frozen natural orbital (FNO) and natural auxiliary function (NAF) cost-reducing approaches are combined with our recent coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] implementations. Both of the closed- and open-shell FNO-CCSD(T) codes benefit from OpenMP parallelism, completely or partially integral-direct density-fitting algorithms, checkpointing, and hand-optimized, memory- and operation count effective implementations exploiting all permutational symmetries. The closed-shell CCSD(T) code requires negligible disk I/O and network bandwidth, is MPI/OpenMP parallel, and exhibits outstanding peak performance utilization of 50-70% up to hundreds of cores. Conservative FNO and NAF truncation thresholds benchmarked for challenging reaction, atomization, and ionization energies of both closed- and open-shell species are shown to maintain 1 kJ/mol accuracy against canonical CCSD(T) for systems of 31-43 atoms even with large basis sets. The cost reduction of up to an order of magnitude achieved extends the reach of FNO-CCSD(T) to systems of 50-75 atoms (up to 2124 atomic orbitals) with triple- and quadruple-ζ basis sets, which is unprecedented without local approximations. Consequently, a considerably larger portion of the chemical compound space can now be covered by the practically "gold standard" quality FNO-CCSD(T) method using affordable resources and about a week of wall time. Large-scale applications are presented for organocatalytic and transition-metal reactions as well as noncovalent interactions. Possible applications for benchmarking local CCSD(T) methods, as well as for the accuracy assessment or parametrization of less complete models, for example, density functional approximations or machine learning potentials, are also outlined.
Collapse
Affiliation(s)
- László Gyevi-Nagy
- Department of Physical Chemistry and
Materials Science, Budapest University of
Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and
Materials Science, Budapest University of
Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Péter R. Nagy
- Department of Physical Chemistry and
Materials Science, Budapest University of
Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
14
|
Schäfer T, Libisch F, Kresse G, Grüneis A. Local embedding of coupled cluster theory into the random phase approximation using plane waves. J Chem Phys 2021; 154:011101. [PMID: 33412868 DOI: 10.1063/5.0036363] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
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 (CC) theory], employing localized orbitals, into a low-level correlation calculation [here, the direct Random Phase Approximation (RPA)]. This choice allows for an accurate and efficient 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 dispersions 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 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.
Collapse
Affiliation(s)
- Tobias Schäfer
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria
| | - Florian Libisch
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria
| | - Georg Kresse
- University of Vienna, Faculty of Physics and Center for Computational Materials Science, Kolingasse 14-16, A-1090 Vienna, Austria
| | - Andreas Grüneis
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria
| |
Collapse
|
15
|
Morales MA, Malone FD. Accelerating the convergence of auxiliary-field quantum Monte Carlo in solids with optimized Gaussian basis sets. J Chem Phys 2020; 153:194111. [DOI: 10.1063/5.0025390] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Miguel A. Morales
- Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Fionn D. Malone
- Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| |
Collapse
|