1
|
Li J, Zhu T. Restoring translational symmetry in periodic all-orbital dynamical mean-field theory simulations. Faraday Discuss 2024. [PMID: 39076013 DOI: 10.1039/d4fd00068d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/31/2024]
Abstract
Dynamical mean-field theory (DMFT) and its cluster extensions provide an efficient Green's function formalism to simulate spectral properties of periodic systems at the quantum many-body level. However, traditional cluster DMFT breaks translational invariance in solid-state materials, and the best strategy to capture non-local correlation effects within cluster DMFT remains elusive. In this work, we investigate the use of overlapping atom-centered impurity fragments in recently-developed ab initio all-orbital DMFT, where all local orbitals within the impurity are treated with high-level quantum chemistry impurity solvers. We demonstrate how the translational symmetry of the lattice self-energy can be restored by designing symmetry-adapted embedding problems, which results in an improved description of spectral functions in two-dimensional boron nitride monolayers and graphene at the levels of many-body perturbation theory (GW) and coupled-cluster theory. Furthermore, we study the convergence of self-energy and density of states as the embedding size is systematically expanded in one-shot and self-consistent DMFT calculations.
Collapse
Affiliation(s)
- Jiachen Li
- Department of Chemistry, Yale University, New Haven, CT, 06520, USA.
| | - Tianyu Zhu
- Department of Chemistry, Yale University, New Haven, CT, 06520, USA.
| |
Collapse
|
2
|
Venturella C, Hillenbrand C, Li J, Zhu T. Machine Learning Many-Body Green's Functions for Molecular Excitation Spectra. J Chem Theory Comput 2024; 20:143-154. [PMID: 38150268 DOI: 10.1021/acs.jctc.3c01146] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2023]
Abstract
We present a machine learning (ML) framework for predicting Green's functions of molecular systems, from which photoemission spectra and quasiparticle energies at quantum many-body level can be obtained. Kernel ridge regression is adopted to predict self-energy matrix elements on compact imaginary frequency grids from static and dynamical mean-field electronic features, which gives direct access to real-frequency many-body Green's functions through analytic continuation and Dyson's equation. Feature and self-energy matrices are represented in a symmetry-adapted intrinsic atomic orbital plus projected atomic orbital basis to enforce rotational invariance. We demonstrate good transferability and high data efficiency of the proposed ML method across molecular sizes and chemical species by showing accurate predictions of density of states (DOS) and quasiparticle energies at the level of many-body perturbation theory (GW) or full configuration interaction. For the ML model trained on 48 out of 1995 molecules randomly sampled from the QM7 and QM9 data sets, we report the mean absolute errors of ML-predicted highest occupied and lowest unoccupied molecular orbital energies to be 0.13 and 0.10 eV, respectively, compared to GW@PBE0. We further showcase the capability of this method by applying the same ML model to predict DOS for significantly larger organic molecules with up to 44 heavy atoms.
Collapse
Affiliation(s)
- Christian Venturella
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | | | - Jiachen Li
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Tianyu Zhu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| |
Collapse
|
3
|
Banerjee S, Sokolov AY. Algebraic Diagrammatic Construction Theory for Simulating Charged Excited States and Photoelectron Spectra. J Chem Theory Comput 2023. [PMID: 37191264 DOI: 10.1021/acs.jctc.3c00251] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/17/2023]
Abstract
Charged excitations are electronic transitions that involve a change in the total charge of a molecule or material. Understanding the properties and reactivity of charged species requires insights from theoretical calculations that can accurately describe orbital relaxation and electron correlation effects in open-shell electronic states. In this Review, we describe the current state of algebraic diagrammatic construction (ADC) theory for simulating charged excitations and its recent developments. We start with a short overview of ADC formalism for the one-particle Green's function, including its single- and multireference formulations and extension to periodic systems. Next, we focus on the capabilities of ADC methods and discuss recent findings about their accuracy for calculating a wide range of excited-state properties. We conclude our Review by outlining possible directions for future developments of this theoretical approach.
Collapse
Affiliation(s)
- Samragni Banerjee
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Alexander Yu Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
4
|
Pathak H, Panyala A, Peng B, Bauman NP, Mutlu E, Rehr JJ, Vila FD, Kowalski K. Real-Time Equation-of-Motion Coupled-Cluster Cumulant Green's Function Method: Heterogeneous Parallel Implementation Based on the Tensor Algebra for Many-Body Methods Infrastructure. J Chem Theory Comput 2023; 19:2248-2257. [PMID: 37096369 DOI: 10.1021/acs.jctc.3c00045] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2023]
Abstract
We report the implementation of the real-time equation-of-motion coupled-cluster (RT-EOM-CC) cumulant Green's function method [ J. Chem. Phys. 2020, 152, 174113] within the Tensor Algebra for Many-body Methods (TAMM) infrastructure. TAMM is a massively parallel heterogeneous tensor library designed for utilizing forthcoming exascale computing resources. The two-body electron repulsion matrix elements are Cholesky-decomposed, and we imposed spin-explicit forms of the various operators when evaluating the tensor contractions. Unlike our previous real algebra Tensor Contraction Engine (TCE) implementation, the TAMM implementation supports fully complex algebra. The RT-EOM-CC singles (S) and doubles (D) time-dependent amplitudes are propagated using a first-order Adams-Moulton method. This new implementation shows excellent scalability tested up to 500 GPUs using the Zn-porphyrin molecule with 655 basis functions, with parallel efficiencies above 90% up to 400 GPUs. The TAMM RT-EOM-CCSD was used to study core photoemission spectra in the formaldehyde and ethyl trifluoroacetate (ESCA) molecules. Simulations of the latter involve as many as 71 occupied and 649 virtual orbitals. The relative quasiparticle ionization energies and overall spectral functions agree well with available experimental results.
Collapse
Affiliation(s)
- Himadri Pathak
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Ajay Panyala
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Bo Peng
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Nicholas P Bauman
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Erdal Mutlu
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - John J Rehr
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - Fernando D Vila
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - Karol Kowalski
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| |
Collapse
|
5
|
Backhouse OJ, Booth GH. Constructing “Full-Frequency” Spectra via Moment Constraints for Coupled Cluster Green’s Functions. J Chem Theory Comput 2022; 18:6622-6636. [DOI: 10.1021/acs.jctc.2c00670] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
| | - George H. Booth
- Department of Physics, King’s College London, Strand, London WC2R 2LS, U.K
| |
Collapse
|
6
|
Laughon K, Yu JM, Zhu T. Periodic Coupled-Cluster Green's Function for Photoemission Spectra of Realistic Solids. J Phys Chem Lett 2022; 13:9122-9128. [PMID: 36162126 DOI: 10.1021/acs.jpclett.2c02534] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
We present an efficient implementation of the coupled-cluster Green's function (CCGF) method for simulating photoemission spectra of periodic systems. We formulate the periodic CCGF approach with Brillouin zone sampling in the Gaussian basis at the coupled-cluster singles and doubles (CCSD) level. To enable CCGF calculations of realistic solids, we propose an active-space self-energy correction scheme by combining CCGF with the cheaper many-body perturbation theory (GW) and implement the model order reduction (MOR) frequency interpolation technique. We find that the active-space self-energy correction and MOR techniques significantly reduce the computational cost of CCGF while maintaining the high accuracy. We apply the developed CCGF approaches to compute spectral properties and band structure of silicon (Si) and zinc oxide (ZnO) crystals using triple-ζ Gaussian basis sets and medium-size k-point sampling and find good agreement with experimental measurements.
Collapse
Affiliation(s)
- Katelyn Laughon
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Jason M Yu
- Department of Chemistry, University of California, Irvine, Irvine, California 92697, United States
| | - Tianyu Zhu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| |
Collapse
|
7
|
Vila FD, Rehr JJ, Pathak H, Peng B, Panyala A, Mutlu E, Bauman NP, Kowalski K. Real-time equation-of-motion CC cumulant and CC Green's function simulations of photoemission spectra of water and water dimer. J Chem Phys 2022; 157:044101. [DOI: 10.1063/5.0099192] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses results obtained with the real-time equation-of-motion CC cumulant approach (RT-EOM-CC), and CC Green's function (CCGF) approaches in applications to the water and water dimer molecules. We compare the ionization potentials obtained with these methods for the valence region with the results obtained with the CCSD(T) formulation as a difference of energies for N and N-1 electron systems. All methods show good agreement with each other. They also agree well with experiment, with errors usually below 0.1 eV for the ionization potentials.We also analyze unique features of the spectral functions, associated with the position of satellite peaks, obtained with the RT-EOM-CC and CCGF methods employing single and double excitations, as a function of the monomer OH bond length and the proton transfer coordinate in the dimer. Finally, we analyze the impact of the basis set effects on the quality of calculated ionization potentials and find that the basis set effects are less pronounced for the augmented-type sets.
Collapse
Affiliation(s)
| | - John J. Rehr
- Department of Physics, University of Washington College of Arts and Sciences, United States of America
| | - Himadri Pathak
- Pacific Northwest National Laboratory, Pacific Northwest National Laboratory, United States of America
| | - Bo Peng
- Pacific Northwest National Laboratory, United States of America
| | - Ajay Panyala
- Pacific Northwest National Laboratory, United States of America
| | - Erdal Mutlu
- Pacific Northwest National Laboratory, United States of America
| | | | - Karol Kowalski
- Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, United States of America
| |
Collapse
|
8
|
Vila FD, Kas JJ, Rehr JJ, Kowalski K, Peng B. Equation-of-Motion Coupled-Cluster Cumulant Green's Function for Excited States and X-Ray Spectra. Front Chem 2021; 9:734945. [PMID: 34631660 PMCID: PMC8493088 DOI: 10.3389/fchem.2021.734945] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2021] [Accepted: 09/06/2021] [Indexed: 11/13/2022] Open
Abstract
Green’s function methods provide a robust, general framework within many-body theory for treating electron correlation in both excited states and x-ray spectra. Conventional methods using the Dyson equation or the cumulant expansion are typically based on the GW self-energy approximation. In order to extend this approximation in molecular systems, a non-perturbative real-time coupled-cluster cumulant Green’s function approach has been introduced, where the cumulant is obtained as the solution to a system of coupled first order, non-linear differential equations. This approach naturally includes non-linear corrections to conventional cumulant Green’s function techniques where the cumulant is linear in the GW self-energy. The method yields the spectral function for the core Green’s function, which is directly related to the x-ray photoemission spectra (XPS) of molecular systems. The approach also yields very good results for binding energies and satellite excitations. The x-ray absorption spectrum (XAS) is then calculated using a convolution of the core spectral function and an effective, one-body XAS. Here this approach is extended to include the full coupled-cluster-singles (CCS) core Green’s function by including the complete form of the non-linear contributions to the cumulant as well as all single, double, and triple cluster excitations in the CC amplitude equations. This approach naturally builds in orthogonality and shake-up effects analogous to those in the Mahan-Noizeres-de Dominicis edge singularity corrections that enhance the XAS near the edge. The method is illustrated for the XPS and XAS of NH3.
Collapse
Affiliation(s)
- F D Vila
- Department of Physics, University of Washington, Seattle, WA, United States
| | - J J Kas
- Department of Physics, University of Washington, Seattle, WA, United States
| | - J J Rehr
- Department of Physics, University of Washington, Seattle, WA, United States
| | - K Kowalski
- Physical and Computational Science Directorate, Pacific Northwest National Laboratory, Richland, WA, United States
| | - B Peng
- Physical and Computational Science Directorate, Pacific Northwest National Laboratory, Richland, WA, United States
| |
Collapse
|
9
|
Banerjee S, Sokolov AY. Efficient implementation of the single-reference algebraic diagrammatic construction theory for charged excitations: Applications to the TEMPO radical and DNA base pairs. J Chem Phys 2021; 154:074105. [PMID: 33607870 DOI: 10.1063/5.0040317] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
We present an efficient implementation of the second- and third-order single-reference algebraic diagrammatic construction (ADC) theory for electron attachment and ionization energies and spectra [EA/IP-ADC(n), n = 2, 3]. Our new EA/IP-ADC program features spin adaptation for closed-shell systems, density fitting for efficient handling of the two-electron integral tensors, and vectorized and parallel implementation of tensor contractions. We demonstrate capabilities of our efficient implementation by applying the EA/IP-ADC(n) (n = 2, 3) methods to compute the photoelectron spectrum of the (2,2,6,6-tetramethylpiperidin-1-yl)oxyl (TEMPO) radical, as well as the vertical and adiabatic electron affinities of TEMPO and two DNA base pairs (guanine-cytosine and adenine-thymine). The spectra and electron affinities computed using large diffuse basis sets with up to 1028 molecular orbitals are found to be in good agreement with the best available results from the experiment and theoretical simulations.
Collapse
Affiliation(s)
- Samragni Banerjee
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Alexander Yu Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
10
|
Bauman NP, Peng B, Kowalski K. Coupled Cluster Green's function formulations based on the effective Hamiltonians. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1725669] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Nicholas P. Bauman
- William R. Wiley Environmental Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory, Richland, WA, USA
| | - Bo Peng
- Advanced Computing, Mathematics, and Data Division, Battelle, Pacific Northwest National Laboratory, Richland, WA, USA
| | - Karol Kowalski
- William R. Wiley Environmental Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory, Richland, WA, USA
| |
Collapse
|
11
|
Banerjee S, Sokolov AY. Third-order algebraic diagrammatic construction theory for electron attachment and ionization energies: Conventional and Green’s function implementation. J Chem Phys 2019; 151:224112. [DOI: 10.1063/1.5131771] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Samragni Banerjee
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Alexander Yu. Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
12
|
Zhu T, Cui ZH, Chan GKL. Efficient Formulation of Ab Initio Quantum Embedding in Periodic Systems: Dynamical Mean-Field Theory. J Chem Theory Comput 2019; 16:141-153. [DOI: 10.1021/acs.jctc.9b00934] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Tianyu Zhu
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Zhi-Hao Cui
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
13
|
Teke NK, Pavošević F, Peng C, Valeev EF. Explicitly correlated renormalized second-order Green’s function for accurate ionization potentials of closed-shell molecules. J Chem Phys 2019; 150:214103. [DOI: 10.1063/1.5090983] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Nakul K. Teke
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Fabijan Pavošević
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Chong Peng
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Edward F. Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| |
Collapse
|
14
|
Peng B, Van Beeumen R, Williams-Young DB, Kowalski K, Yang C. Approximate Green's Function Coupled Cluster Method Employing Effective Dimension Reduction. J Chem Theory Comput 2019; 15:3185-3196. [PMID: 30951302 DOI: 10.1021/acs.jctc.9b00172] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The Green's function coupled cluster (GFCC) method, originally proposed in the early 1990s, is a powerful many-body tool for computing and analyzing the electronic structure of molecular and periodic systems, especially when electrons of the system are strongly correlated. However, in order for the GFCC to become a method that may be routinely used in the electronic structure calculations, robust numerical techniques and approximations must be employed to reduce its extremely high computational overhead. In our recent studies, it has been demonstrated that the GFCC equations can be solved directly in the frequency domain using iterative linear solvers, which can be easily distributed in a massively parallel environment. In the present work, we demonstrate a successful application of model-order-reduction (MOR) techniques in the GFCC framework. Briefly speaking, for a frequency regime of interest that requires high-resolution descriptions of spectral function, instead of solving the GFCC linear equation of full dimension for every single frequency point of interest, an efficiently solvable linear system model of a reduced dimension may be built upon projecting the original GFCC linear system onto a subspace. From this reduced order model is obtained a reasonable approximation to the full dimensional GFCC linear equations in both interpolative and extrapolative spectral regions. Here, we show that the subspace can be properly constructed in an iterative manner from the auxiliary vectors of the GFCC linear equations at some selected frequencies within the spectral region of interest. During the iterations, the quality of the subspace, as well as the linear system model, can be systematically improved. The method is tested in this work in terms of the efficiency and accuracy of computing spectral functions for some typical molecular systems such as carbon monoxide, 1,3-butadiene, benzene, and adenine. To reach the same level of accuracy as that of the original GFCC method, the application of MOR in the GFCC method is able to significantly lower the original computational cost for the aforementioned molecules in designated frequency regimes. As a byproduct, the reduced order model obtained by this method is found to provide a high-quality initial guess, which improves the convergence rate for the existing iterative linear solver.
Collapse
Affiliation(s)
- Bo Peng
- William R. Wiley Environmental Molecular Sciences Laboratory, Battelle , Pacific Northwest National Laboratory , K8-91, P.O. Box 999, Richland , Washington 99352 , United States
| | - Roel Van Beeumen
- Computational Research Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| | - David B Williams-Young
- Computational Research Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| | - Karol Kowalski
- William R. Wiley Environmental Molecular Sciences Laboratory, Battelle , Pacific Northwest National Laboratory , K8-91, P.O. Box 999, Richland , Washington 99352 , United States
| | - Chao Yang
- Computational Research Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| |
Collapse
|
15
|
Kosugi T, Matsushita YI. Wannier interpolation of one-particle Green's functions from coupled-cluster singles and doubles (CCSD). J Chem Phys 2019; 150:114104. [PMID: 30902011 DOI: 10.1063/1.5079474] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We propose two schemes for interpolation of the one-particle Green's function (GF) calculated within a coupled-cluster singles and doubles (CCSD) method for a periodic system. These schemes use Wannier orbitals for circumventing huge cost for a large number of sampled k points. One of the schemes is the direct interpolation, which obtains the GF straightforwardly by using Fourier transformation. The other is the self-energy-mediated interpolation, which obtains the GF via the Dyson equation. We apply the schemes to a LiH chain and trans-polyacetylene and examine their validity in detail. It is demonstrated that the direct-interpolated GFs suffer from numerical artifacts stemming from slow convergence of CCSD GFs in real space, while the self-energy-mediated interpolation provides more physically appropriate GFs due to the localized nature of CCSD self-energies. Our schemes are also applicable to other correlated methods capable of providing GFs.
Collapse
Affiliation(s)
- Taichi Kosugi
- Department of Physics, University of Tokyo, Tokyo 113-0033, Japan
| | - Yu-Ichiro Matsushita
- Laboratory for Materials and Structures, Institute of Innovative Research, Tokyo Institute of Technology, Yokohama 226-8503, Japan
| |
Collapse
|
16
|
Peng B, Kowalski K. Green's function coupled cluster formulations utilizing extended inner excitations. J Chem Phys 2018; 149:214102. [PMID: 30525725 DOI: 10.1063/1.5046529] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper, we analyze new approximations of the Green's function coupled cluster (GFCC) method where locations of poles are improved by extending the excitation level of inner auxiliary operators. These new GFCC approximations can be categorized as the GFCC-i(n, m) method, where the excitation level of the inner auxiliary operators (m) used to describe the ionization potential and electron affinity effects in the N - 1 and N + 1 particle spaces is higher than the excitation level (n) used to correlate the ground-state coupled cluster wave function for the N-electron system. Furthermore, we reveal the so-called "n + 1" rule in this category [or the GFCC-i(n, n + 1) method], which states that in order to maintain size-extensivity of the Green's function matrix elements, the excitation level of inner auxiliary operators X p (ω) and Y q (ω) cannot exceed n + 1. We also discuss the role of the moments of coupled cluster equations that in a natural way assures these properties. Our implementation in the present study is focused on the first approximation in this GFCC category, i.e., the GFCC-i(2,3) method. As our first practice, we use the GFCC-i(2,3) method to compute the spectral functions for the N2 and CO molecules in the inner and outer valence regimes. In comparison with the Green's function coupled cluster singles, doubles results, the computed spectral functions from the GFCC-i(2,3) method exhibit better agreement with the experimental results and other theoretical results, particularly in terms of providing higher resolution of satellite peaks and more accurate relative positions of these satellite peaks with respect to the main peak positions.
Collapse
Affiliation(s)
- Bo Peng
- William R. Wiley Environmental Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory, K8-91, P.O. Box 999, Richland, Washington 99352, USA
| | - Karol Kowalski
- William R. Wiley Environmental Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory, K8-91, P.O. Box 999, Richland, Washington 99352, USA
| |
Collapse
|
17
|
Sokolov AY. Multi-reference algebraic diagrammatic construction theory for excited states: General formulation and first-order implementation. J Chem Phys 2018; 149:204113. [DOI: 10.1063/1.5055380] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Alexander Yu. Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
18
|
Nishi H, Kosugi T, Furukawa Y, Matsushita YI. Quasiparticle energy spectra of isolated atoms from coupled-cluster singles and doubles (CCSD): Comparison with exact CI calculations. J Chem Phys 2018; 149:034106. [DOI: 10.1063/1.5029536] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Hirofumi Nishi
- Department of Applied Physics, The University of Tokyo, Tokyo 113-8656, Japan
| | - Taichi Kosugi
- Department of Applied Physics, The University of Tokyo, Tokyo 113-8656, Japan
| | - Yoritaka Furukawa
- Department of Applied Physics, The University of Tokyo, Tokyo 113-8656, Japan
| | | |
Collapse
|
19
|
Lange MF, Berkelbach TC. On the Relation between Equation-of-Motion Coupled-Cluster Theory and the GW Approximation. J Chem Theory Comput 2018; 14:4224-4236. [DOI: 10.1021/acs.jctc.8b00455] [Citation(s) in RCA: 45] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Malte F. Lange
- Department of Chemistry and James Franck Institute, University of Chicago, Chicago, Illinois 60637, United States
| | - Timothy C. Berkelbach
- Department of Chemistry and James Franck Institute, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
20
|
Peng B, Kowalski K. Green's Function Coupled-Cluster Approach: Simulating Photoelectron Spectra for Realistic Molecular Systems. J Chem Theory Comput 2018; 14:4335-4352. [PMID: 29957945 DOI: 10.1021/acs.jctc.8b00313] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this paper, we present an efficient implementation for the analytical energy-dependent Green's function coupled-cluster with singles and doubles (GFCCSD) approach with our first practice being computing spectral functions of realistic molecular systems. Because of its algebraic structure, the presented method is highly scalable and is capable of computing spectral function for a given molecular system in any energy region. Several typical examples have been given to demonstrate its capability of computing spectral functions not only in the valence band but also in the core-level energy region. Satellite peaks have been observed in the inner valence band and core-level energy region where a many-body effect becomes significant and the single particle picture of ionization often breaks down. The accuracy test has been carried out by extensively comparing the computed spectral functions by our GFCCSD method with experimental photoelectron spectra as well as the theoretical ionization potentials obtained from other methods. It turns out the GFCCSD method is able to provide a qualitative or semiquantitative level of description of ionization processes in both the core and valence regimes. To significantly improve the GFCCSD results for the main ionic states, a larger basis set can usually be employed, whereas the improvement of the GFCCSD results for the satellite states needs higher-order many-body terms to be included in the GFCC implementation.
Collapse
Affiliation(s)
- Bo Peng
- William R. Wiley Environmental Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory , K8-91, P.O. Box 999, Richland , Washington 99352 , United States
| | - Karol Kowalski
- William R. Wiley Environmental Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory , K8-91, P.O. Box 999, Richland , Washington 99352 , United States
| |
Collapse
|
21
|
Kosugi T, Nishi H, Furukawa Y, Matsushita YI. Comparison of Green’s functions for transition metal atoms using self-energy functional theory and coupled-cluster singles and doubles (CCSD). J Chem Phys 2018; 148:224103. [DOI: 10.1063/1.5029535] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Taichi Kosugi
- Department of Applied Physics, The University of Tokyo, Tokyo 113-8656, Japan
| | - Hirofumi Nishi
- Department of Applied Physics, The University of Tokyo, Tokyo 113-8656, Japan
| | - Yoritaka Furukawa
- Department of Applied Physics, The University of Tokyo, Tokyo 113-8656, Japan
| | | |
Collapse
|
22
|
Hoy EP, Mazziotti DA, Seideman T. Development and application of a 2-electron reduced density matrix approach to electron transport via molecular junctions. J Chem Phys 2018; 147:184110. [PMID: 29141419 DOI: 10.1063/1.4986804] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Can an electronic device be constructed using only a single molecule? Since this question was first asked by Aviram and Ratner in the 1970s [Chem. Phys. Lett. 29, 277 (1974)], the field of molecular electronics has exploded with significant experimental advancements in the understanding of the charge transport properties of single molecule devices. Efforts to explain the results of these experiments and identify promising new candidate molecules for molecular devices have led to the development of numerous new theoretical methods including the current standard theoretical approach for studying single molecule charge transport, i.e., the non-equilibrium Green's function formalism (NEGF). By pairing this formalism with density functional theory (DFT), a wide variety of transport problems in molecular junctions have been successfully treated. For some systems though, the conductance and current-voltage curves predicted by common DFT functionals can be several orders of magnitude above experimental results. In addition, since density functional theory relies on approximations to the exact exchange-correlation functional, the predicted transport properties can show significant variation depending on the functional chosen. As a first step to addressing this issue, the authors have replaced density functional theory in the NEGF formalism with a 2-electron reduced density matrix (2-RDM) method, creating a new approach known as the NEGF-RDM method. 2-RDM methods provide a more accurate description of electron correlation compared to density functional theory, and they have lower computational scaling compared to wavefunction based methods of similar accuracy. Additionally, 2-RDM methods are capable of capturing static electron correlation which is untreatable by existing NEGF-DFT methods. When studying dithiol alkane chains and dithiol benzene in model junctions, the authors found that the NEGF-RDM predicts conductances and currents that are 1-2 orders of magnitude below those of B3LYP and M06 DFT functionals. This suggests that the NEGF-RDM method could be a viable alternative to NEGF-DFT for molecular junction calculations.
Collapse
Affiliation(s)
- Erik P Hoy
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, USA
| | - David A Mazziotti
- Department of Chemistry and The James Frank Institute, University of Chicago, Chicago, Illinois 60637, USA
| | - Tamar Seideman
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, USA
| |
Collapse
|
23
|
Furukawa Y, Kosugi T, Nishi H, Matsushita YI. Band structures in coupled-cluster singles-and-doubles Green’s function (GFCCSD). J Chem Phys 2018; 148:204109. [DOI: 10.1063/1.5029537] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Yoritaka Furukawa
- Department of Applied Physics, The University of Tokyo, Tokyo 113-8656, Japan
| | - Taichi Kosugi
- Department of Applied Physics, The University of Tokyo, Tokyo 113-8656, Japan
| | - Hirofumi Nishi
- Department of Applied Physics, The University of Tokyo, Tokyo 113-8656, Japan
| | | |
Collapse
|
24
|
Warneke J, Hou GL, Aprà E, Jenne C, Yang Z, Qin Z, Kowalski K, Wang XB, Xantheas SS. Electronic Structure and Stability of [B12X12]2– (X = F–At): A Combined Photoelectron Spectroscopic and Theoretical Study. J Am Chem Soc 2017; 139:14749-14756. [DOI: 10.1021/jacs.7b08598] [Citation(s) in RCA: 46] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jonas Warneke
- Physical
Sciences Division, Pacific Northwest National Laboratory, 902 Battelle
Boulevard, P.O. Box 999, MS K8-88, Richland, Washington 99352, United States
| | - Gao-Lei Hou
- Physical
Sciences Division, Pacific Northwest National Laboratory, 902 Battelle
Boulevard, P.O. Box 999, MS K8-88, Richland, Washington 99352, United States
| | - Edoardo Aprà
- Environmental
Molecular Sciences Laboratory, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352 United States
| | - Carsten Jenne
- Anorganische
Chemie, Fakultät für Mathematik und Naturwissenschaften, Bergische Universität Wuppertal, Gaußstraße 20, 42119 Wuppertal, Germany
| | - Zheng Yang
- Physical
Sciences Division, Pacific Northwest National Laboratory, 902 Battelle
Boulevard, P.O. Box 999, MS K8-88, Richland, Washington 99352, United States
| | - Zhengbo Qin
- Physical
Sciences Division, Pacific Northwest National Laboratory, 902 Battelle
Boulevard, P.O. Box 999, MS K8-88, Richland, Washington 99352, United States
| | - Karol Kowalski
- Environmental
Molecular Sciences Laboratory, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352 United States
| | - Xue-Bin Wang
- Physical
Sciences Division, Pacific Northwest National Laboratory, 902 Battelle
Boulevard, P.O. Box 999, MS K8-88, Richland, Washington 99352, United States
| | - Sotiris S. Xantheas
- Advanced
Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box
999, MS K1-83, Richland, Washington 99352, United States
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
25
|
Maggio E, Kresse G. GW Vertex Corrected Calculations for Molecular Systems. J Chem Theory Comput 2017; 13:4765-4778. [DOI: 10.1021/acs.jctc.7b00586] [Citation(s) in RCA: 36] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Emanuele Maggio
- Faculty of Physics
and Center
for Computational Materials Science, University of Vienna, Sensengasse
8/12, A-1090 Vienna, Austria
| | - Georg Kresse
- Faculty of Physics
and Center
for Computational Materials Science, University of Vienna, Sensengasse
8/12, A-1090 Vienna, Austria
| |
Collapse
|
26
|
Hirata S, Doran AE, Knowles PJ, Ortiz JV. One-particle many-body Green’s function theory: Algebraic recursive definitions, linked-diagram theorem, irreducible-diagram theorem, and general-order algorithms. J Chem Phys 2017; 147:044108. [DOI: 10.1063/1.4994837] [Citation(s) in RCA: 44] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Alexander E. Doran
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Peter J. Knowles
- School of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, United Kingdom
| | - J. V. Ortiz
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849-5312, USA
| |
Collapse
|
27
|
Affiliation(s)
- Bo Peng
- William R. Wiley Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, WA, USA
| | - Karol Kowalski
- William R. Wiley Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, WA, USA
| |
Collapse
|
28
|
Berardo E, Kaplan F, Bhaskaran-Nair K, Shelton WA, van Setten MJ, Kowalski K, Zwijnenburg MA. Benchmarking the Fundamental Electronic Properties of small TiO2 Nanoclusters by GW and Coupled Cluster Theory Calculations. J Chem Theory Comput 2017. [DOI: 10.1021/acs.jctc.7b00538] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Enrico Berardo
- Department
of Chemistry, Imperial College London, South Kensington, London SW7 2AZ, United Kingdom
- Department
of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Ferdinand Kaplan
- Institute
of Nanotechnology, Karlsruhe Institute of Technology, Campus North, D-76344 Karlsruhe, Germany
| | - Kiran Bhaskaran-Nair
- Cain
Department of Chemical Engineering, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - William A. Shelton
- Cain
Department of Chemical Engineering, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - Michiel J. van Setten
- Nanoscopic
Physics, Institute of Condensed Matter and Nanosciences, Université Catholique de Louvain, 1348 Louvain-la-Neuve, Belgium
| | - Karol Kowalski
- William R.
Wiley Environmental Molecular Science Laboratory, Battelle, Pacific Northwest National Laboratory, K8-91, P. O. Box 999, Richland, Washington 99352, United States
| | - Martijn A. Zwijnenburg
- Department
of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
29
|
Peng B, Govind N, Aprà E, Klemm M, Hammond JR, Kowalski K. Coupled Cluster Studies of Ionization Potentials and Electron Affinities of Single-Walled Carbon Nanotubes. J Phys Chem A 2017; 121:1328-1335. [DOI: 10.1021/acs.jpca.6b10874] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Bo Peng
- William
R. Wiley Environmental Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory, K8-91, P.O. Box
999, Richland, Washington 99352, United States
| | - Niranjan Govind
- William
R. Wiley Environmental Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory, K8-91, P.O. Box
999, Richland, Washington 99352, United States
| | - Edoardo Aprà
- William
R. Wiley Environmental Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory, K8-91, P.O. Box
999, Richland, Washington 99352, United States
| | - Michael Klemm
- Intel Deutschland GmbH, Feldkirchen, Munich 08016, Germany
| | | | - Karol Kowalski
- William
R. Wiley Environmental Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory, K8-91, P.O. Box
999, Richland, Washington 99352, United States
| |
Collapse
|