1
|
Sereda M, Allen T, Bradbury NC, Ibrahim KZ, Neuhauser D. Sparse-Stochastic Fragmented Exchange for Large-Scale Hybrid Time-Dependent Density Functional Theory Calculations. J Chem Theory Comput 2024; 20:4196-4204. [PMID: 38713513 DOI: 10.1021/acs.jctc.4c00260] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/09/2024]
Abstract
We extend our recently developed sparse-stochastic fragmented exchange formalism for ground-state near-gap hybrid DFT to calculate absorption spectra within linear-response time-dependent generalized Kohn-Sham DFT (LR-GKS-TDDFT) for systems consisting of thousands of valence electrons within a grid-based/plane-wave representation. A mixed deterministic/fragmented-stochastic compression of the exchange kernel, here using long-range explicit exchange functionals, provides an efficient method for accurate optical spectra. Both real-time propagation as well as frequency-resolved Casida-equation-type approaches for spectra are presented, and the method is applied to large molecular dyes.
Collapse
Affiliation(s)
- Mykola Sereda
- Department of Chemistry and Biochemistry, and California Nanoscience Institute, UCLA, Los Angeles, California 90095-1569, United States
| | - Tucker Allen
- Department of Chemistry and Biochemistry, and California Nanoscience Institute, UCLA, Los Angeles, California 90095-1569, United States
| | - Nadine C Bradbury
- Department of Chemistry and Biochemistry, and California Nanoscience Institute, UCLA, Los Angeles, California 90095-1569, United States
| | - Khaled Z Ibrahim
- Computer Science Department, Lawrence Berkeley National Laboratory, One Cyclotron Road, Berkeley, California 94720, United States
| | - Daniel Neuhauser
- Department of Chemistry and Biochemistry, and California Nanoscience Institute, UCLA, Los Angeles, California 90095-1569, United States
| |
Collapse
|
2
|
Zhang Z, Yin X, Hu W, Yang J. Machine Learning K-Means Clustering of Interpolative Separable Density Fitting Algorithm for Accurate and Efficient Cubic-Scaling Exact Exchange Plus Random Phase Approximation within Plane Waves. J Chem Theory Comput 2024; 20:1944-1961. [PMID: 38361423 DOI: 10.1021/acs.jctc.3c01157] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2024]
Abstract
The exact-exchange plus random-phase approximation (EXX+RPA) method has emerged as a crucial tool for precisely characterizing electronic structures in molecular and solid systems. We present an accurate and efficient implementation of EXX+RPA calculations that scale cubically and are conducted within plane waves. Our approach incorporates the interpolative separable density fitting (ISDF) algorithm, effectively mitigating the computational challenges associated with the plane wave basis set. To overcome the constraints of the conventional ISDF algorithm, characterized by the exceptionally high prefactor in QR factorization for interpolation point selection, we introduce an enhanced machine learning K-means method. This method incorporates a novel empirical weight function called "SSM+" for more precise interpolation point selection, capturing physical information more accurately across diverse systems. Our machine learning approach offers a quasiquadratic scaling alternative, effectively replacing the computationally demanding cubic-scaling QRCP algorithm in plane-wave-based EXX+RPA calculations. Furthermore, we enhance the method's capabilities by optimizing GPU acceleration using MATLAB's integrated GPU toolkit. In particular, our approach reduces the computational scaling of χ0 from 3.80 to 2.13 and the overall computational scaling of EXX from 2.74 to 2.10. We achieve a remarkable GPU acceleration speedup of up to 35×. Regarding CPU computation time, the standard quartic-scaling method requires 22 h to compute Si128, while QRCP completes the calculation in only around 1 h, achieving a speedup up to 20×. However, the utilization of the K-means algorithm reduces the time to 800 s, a substantial improvement of 100× compared to the standard algorithm. By employing the K-means algorithm, the computational time for interpolative point calculation using QRCP decreases from 1 h to 1 min, resulting in a 55× speed increase. With this improved algorithm, we successfully computed the dissociation curve of H2 and the equilibrium polyynic geometry of C18 molecules.
Collapse
Affiliation(s)
- Zhenlin Zhang
- Key Laboratory of Precision and Intelligent Chemistry, Department of Chemical Physics, and Anhui Center for Applied Mathematics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Xilin Yin
- Key Laboratory of Precision and Intelligent Chemistry, Department of Chemical Physics, and Anhui Center for Applied Mathematics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Wei Hu
- Key Laboratory of Precision and Intelligent Chemistry, Department of Chemical Physics, and Anhui Center for Applied Mathematics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Jinlong Yang
- Key Laboratory of Precision and Intelligent Chemistry, Department of Chemical Physics, and Anhui Center for Applied Mathematics, University of Science and Technology of China, Hefei, Anhui 230026, China
| |
Collapse
|
3
|
Wang L, Liu X, Wan L, Gao Y, Wang X, Liu J, Tan S, Guo Q, Zhao W, Hu W, Li Q, Yang J. Excitonic Effects of the Excited-State Photocatalytic Reaction at the Molecule/Metal Oxide Interface. J Phys Chem Lett 2024; 15:2096-2104. [PMID: 38358755 DOI: 10.1021/acs.jpclett.3c03075] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/16/2024]
Abstract
Excitonic effects caused by the Coulomb interaction between electrons and holes play a crucial role in photocatalysis at the molecule/metal oxide interface. As an ideal model for investigating the excitonic effect, coadsorption and photodissociation of water and methanol molecules on titanium dioxide involve complex ground-state thermalcatalytic and excited-state photocatalytic reaction processes. Herein, we systemically investigate the excited-state electronic structures of the coadsorption of H2O and CH3OH molecules on a rutile TiO2(110) surface by linear-response time-dependent density functional theory calculations and probe the reaction path for generating HCOOH or CO2, from ground-state and excited-state perspectives. The reaction barriers in excited-state calculations are significantly different from those in ground-state calculations during three processes, with the largest decrease being 0.94 eV for the Ti5c-O-CH2-O-Ti5c formation process.
Collapse
Affiliation(s)
- Lei Wang
- Department of Chemical Physics, and Key Laboratory of Precision and Intelligent Chemistry, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Xiaofeng Liu
- School of Physics, Hefei University of Technology, Hefei, Anhui 230009, China
| | - Lingyun Wan
- Hefei National Research Center for Physical Sciences at the Microscale and Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Yunzhi Gao
- Hefei National Research Center for Physical Sciences at the Microscale and Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Xiaoning Wang
- Department of Chemical Physics, and Key Laboratory of Precision and Intelligent Chemistry, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Jie Liu
- Hefei National Research Center for Physical Sciences at the Microscale and Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Shijing Tan
- Hefei National Research Center for Physical Sciences at the Microscale and Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Qing Guo
- Department of Chemistry, Southern University of Science and Technology, Shenzhen, Guangdong 518055, China
| | - Wenhui Zhao
- Department of Physics, Ningbo University, Ningbo, Zhejiang 315211, China
| | - Wei Hu
- Hefei National Research Center for Physical Sciences at the Microscale and Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Qunxiang Li
- Department of Chemical Physics, and Key Laboratory of Precision and Intelligent Chemistry, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Jinlong Yang
- Department of Chemical Physics, and Key Laboratory of Precision and Intelligent Chemistry, University of Science and Technology of China, Hefei, Anhui 230026, China
| |
Collapse
|
4
|
Liu J, Lu G, Zhang X. Exciton dispersion and exciton-phonon interaction in solids by time-dependent density functional theory. J Chem Phys 2023; 158:044116. [PMID: 36725491 DOI: 10.1063/5.0137326] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023] Open
Abstract
Understanding, predicting, and ultimately controlling exciton band structure and exciton dynamics are central to diverse chemical and materials problems. Here, we have developed a first-principles method to determine exciton dispersion and exciton-phonon interaction in semiconducting and insulating solids based on time-dependent density functional theory. The first-principles method is formulated in planewave bases and pseudopotentials and can be used to compute exciton band structures, exciton charge density, ionic forces, the non-adiabatic coupling matrix between excitonic states, and the exciton-phonon coupling matrix. Based on the spinor formulation, the method enables self-consistent noncollinear calculations to capture spin-orbital coupling. Hybrid exchange-correlation functionals are incorporated to deal with long-range electron-hole interactions in solids. A sub-Hilbert space approximation is introduced to reduce the computational cost without loss of accuracy. For validations, we have applied the method to compute the exciton band structure and exciton-phonon coupling strength in transition metal dichalcogenide monolayers; both agree very well with the previous GW-Bethe-Salpeter equation and experimental results. This development paves the way for accurate determinations of exciton dynamics in a wide range of solid-state materials.
Collapse
Affiliation(s)
- Junyi Liu
- Department of Physics and Astronomy, California State University Northridge, Northridge, California 91330, USA
| | - Gang Lu
- Department of Physics and Astronomy, California State University Northridge, Northridge, California 91330, USA
| | - Xu Zhang
- Department of Physics and Astronomy, California State University Northridge, Northridge, California 91330, USA
| |
Collapse
|