1
|
Williams-Young DB, Yuwono SH, DePrince III AE, Yang C. Approximate Exponential Integrators for Time-Dependent Equation-of-Motion Coupled Cluster Theory. J Chem Theory Comput 2023; 19:9177-9186. [PMID: 38086060 PMCID: PMC10753770 DOI: 10.1021/acs.jctc.3c00911] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2023] [Revised: 10/09/2023] [Accepted: 10/16/2023] [Indexed: 12/27/2023]
Abstract
With a growing demand for time-domain simulations of correlated many-body systems, the development of efficient and stable integration schemes for the time-dependent Schrödinger equation is of keen interest in modern electronic structure theory. In this work, we present two approaches for the formation of the quantum propagator for time-dependent equation-of-motion coupled cluster theory based on the Chebyshev and Arnoldi expansions of the complex, nonhermitian matrix exponential, respectively. The proposed algorithms are compared with the short-iterative Lanczos method of Cooper et al. [J. Phys. Chem. A 2021 125, 5438-5447], the fourth-order Runge-Kutta method, and exact dynamics for a set of small but challenging test problems. For each of the cases studied, both of the proposed integration schemes demonstrate superior accuracy and efficiency relative to the reference simulations.
Collapse
Affiliation(s)
- David B. Williams-Young
- Applied
Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Stephen H. Yuwono
- Department
of Chemistry and Biochemistry, Florida State
University, Tallahassee, Florida 32306, United States
| | - A. Eugene DePrince III
- Department
of Chemistry and Biochemistry, Florida State
University, Tallahassee, Florida 32306, United States
| | - Chao Yang
- Applied
Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
2
|
Perrella F, Coppola F, Rega N, Petrone A. An Expedited Route to Optical and Electronic Properties at Finite Temperature via Unsupervised Learning. Molecules 2023; 28:molecules28083411. [PMID: 37110644 PMCID: PMC10144358 DOI: 10.3390/molecules28083411] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2023] [Revised: 04/06/2023] [Accepted: 04/07/2023] [Indexed: 04/29/2023] Open
Abstract
Electronic properties and absorption spectra are the grounds to investigate molecular electronic states and their interactions with the environment. Modeling and computations are required for the molecular understanding and design strategies of photo-active materials and sensors. However, the interpretation of such properties demands expensive computations and dealing with the interplay of electronic excited states with the conformational freedom of the chromophores in complex matrices (i.e., solvents, biomolecules, crystals) at finite temperature. Computational protocols combining time dependent density functional theory and ab initio molecular dynamics (MD) have become very powerful in this field, although they require still a large number of computations for a detailed reproduction of electronic properties, such as band shapes. Besides the ongoing research in more traditional computational chemistry fields, data analysis and machine learning methods have been increasingly employed as complementary approaches for efficient data exploration, prediction and model development, starting from the data resulting from MD simulations and electronic structure calculations. In this work, dataset reduction capabilities by unsupervised clustering techniques applied to MD trajectories are proposed and tested for the ab initio modeling of electronic absorption spectra of two challenging case studies: a non-covalent charge-transfer dimer and a ruthenium complex in solution at room temperature. The K-medoids clustering technique is applied and is proven to be able to reduce by ∼100 times the total cost of excited state calculations on an MD sampling with no loss in the accuracy and it also provides an easier understanding of the representative structures (medoids) to be analyzed on the molecular scale.
Collapse
Affiliation(s)
- Fulvio Perrella
- Scuola Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
| | - Federico Coppola
- Scuola Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
| | - Nadia Rega
- Scuola Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M.S. Angelo, via Cintia 21, I-80126 Napoli, Italy
- Istituto Nazionale di Fisica Nucleare, Sezione di Napoli, Complesso Universitario di M.S. Angelo ed. 6, via Cintia 21, I-80126 Napoli, Italy
| | - Alessio Petrone
- Scuola Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M.S. Angelo, via Cintia 21, I-80126 Napoli, Italy
- Istituto Nazionale di Fisica Nucleare, Sezione di Napoli, Complesso Universitario di M.S. Angelo ed. 6, via Cintia 21, I-80126 Napoli, Italy
| |
Collapse
|
3
|
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
|
4
|
Nascimento DR, Govind N. Computational approaches for XANES, VtC-XES, and RIXS using linear-response time-dependent density functional theory based methods. Phys Chem Chem Phys 2022; 24:14680-14691. [PMID: 35699090 DOI: 10.1039/d2cp01132h] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
The emergence of state-of-the-art X-ray light sources has paved the way for novel spectroscopies that take advantage of their atomic specificity to shed light on fundamental physical, chemical, and biological processes both in the static and time domains. The success of these experiments hinges on the ability to interpret and predict core-level spectra, which has opened avenues for theory to play a key role. Over the last two decades, linear-response time-dependent density functional theory (LR-TDDFT), despite various theoretical challenges, has become a computationally attractive and versatile framework to study excited-state spectra including X-ray spectroscopies. In this context, we focus our discussion on LR-TDDFT approaches for the computation of X-ray Near-Edge Structure (XANES), Valence-to-Core X-ray Emission (VtC-XES), and Resonant Inelastic X-ray Scattering (RIXS) spectroscopies in molecular systems with an emphasis on Gaussian basis set implementations. We illustrate these approaches with applications and provide a brief outlook of possible new directions.
Collapse
Affiliation(s)
- Daniel R Nascimento
- Department of Chemistry, The University of Memphis, Memphis, TN, 38152, USA.
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, USA.
| |
Collapse
|
5
|
Li X, Govind N, Isborn C, DePrince AE, Lopata K. Real-Time Time-Dependent Electronic Structure Theory. Chem Rev 2020; 120:9951-9993. [DOI: 10.1021/acs.chemrev.0c00223] [Citation(s) in RCA: 64] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195-1700, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Christine Isborn
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, United States
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - Kenneth Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| |
Collapse
|
6
|
Williams‐Young DB, Petrone A, Sun S, Stetina TF, Lestrange P, Hoyer CE, Nascimento DR, Koulias L, Wildman A, Kasper J, Goings JJ, Ding F, DePrince AE, Valeev EF, Li X. The Chronus Quantum software package. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1436] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- David B. Williams‐Young
- Computational Research Division Lawrence Berkeley National Laboratory Berkeley California
- Department of Chemistry University of Washington Seattle Washington
| | - Alessio Petrone
- Dipartimento di Scienze Chimiche Università di Napoli “Federico II”, Complesso Universitario di M.S. Angelo Naples Italy
| | - Shichao Sun
- Department of Chemistry University of Washington Seattle Washington
| | - Torin F. Stetina
- Department of Chemistry University of Washington Seattle Washington
| | | | - Chad E. Hoyer
- Department of Chemistry University of Washington Seattle Washington
| | - Daniel R. Nascimento
- Department of Chemistry and Biochemistry Florida State University Tallahassee Florida
| | - Lauren Koulias
- Department of Chemistry University of Washington Seattle Washington
| | - Andrew Wildman
- Department of Chemistry University of Washington Seattle Washington
| | - Joseph Kasper
- Department of Chemistry University of Washington Seattle Washington
| | | | - Feizhi Ding
- Division of Chemistry and Chemical Engineering California Institute of Technology Pasadena California
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry Florida State University Tallahassee Florida
| | | | - Xiaosong Li
- Department of Chemistry University of Washington Seattle Washington
| |
Collapse
|
7
|
Stetina TF, Kasper JM, Li X. Modeling L2,3-edge X-ray absorption spectroscopy with linear response exact two-component relativistic time-dependent density functional theory. J Chem Phys 2019; 150:234103. [DOI: 10.1063/1.5091807] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Torin F. Stetina
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Joseph M. Kasper
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| |
Collapse
|
8
|
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
|
9
|
Kasper JM, Lestrange PJ, Stetina TF, Li X. Modeling L2,3-Edge X-ray Absorption Spectroscopy with Real-Time Exact Two-Component Relativistic Time-Dependent Density Functional Theory. J Chem Theory Comput 2018; 14:1998-2006. [DOI: 10.1021/acs.jctc.7b01279] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Joseph M. Kasper
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Patrick J. Lestrange
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Torin F. Stetina
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
10
|
Kasper JM, Williams-Young DB, Vecharynski E, Yang C, Li X. A Well-Tempered Hybrid Method for Solving Challenging Time-Dependent Density Functional Theory (TDDFT) Systems. J Chem Theory Comput 2018; 14:2034-2041. [DOI: 10.1021/acs.jctc.8b00141] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Joseph M. Kasper
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | | | - Eugene Vecharynski
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Chao Yang
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|