1
|
Yuan K, Zhou S, Li N, Li T, Ding B, Guo D, Ma Y. Fault-tolerant quantum chemical calculations with improved machine-learning models. J Comput Chem 2024; 45:2640-2658. [PMID: 39072777 DOI: 10.1002/jcc.27459] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2023] [Revised: 05/30/2024] [Accepted: 06/18/2024] [Indexed: 07/30/2024]
Abstract
Easy and effective usage of computational resources is crucial for scientific calculations. Following our recent work of machine-learning (ML) assisted scheduling optimization [J. Comput. Chem. 2023, 44, 1174], we further propose (1) the improved ML models for the better predictions of computational loads, and as such, more elaborate load-balancing calculations can be expected; (2) the idea of coded computation, that is, the integration of gradient coding, in order to introduce fault tolerance during the distributed calculations; and (3) their applications together with re-normalized exciton model with time-dependent density functional theory (REM-TDDFT) for calculating the excited states. Illustrated benchmark calculations include P38 protein, and solvent model with one or several excitable centers. The results show that the improved ML-assisted coded calculations can further improve the load-balancing and cluster utilization, owing primarily profit in fault tolerance that aims at the automated quantum chemical calculations for both ground and excited states.
Collapse
Affiliation(s)
- Kai Yuan
- College of Information Science and Technology, Beijing University of Chemical Technology, Beijing, China
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - Shuai Zhou
- Institute of Computing Technology, Chinese Academy of Sciences, Beijing, China
- School of Computer Science and Technology, University of Chinese Academy of Sciences, Beijing, China
| | - Ning Li
- College of Chemistry and Materials Engineering, Wenzhou University, Wenzhou, China
| | - Tianyan Li
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - Bowen Ding
- Institute of Chemistry, Chinese Academy of Sciences, Beijing, China
| | - Danhuai Guo
- College of Information Science and Technology, Beijing University of Chemical Technology, Beijing, China
| | - Yingjin Ma
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| |
Collapse
|
2
|
Wang K, Xu Y, Xie X, Ma H. Theoretical investigation of distal charge separation in a perylenediimide trimer. J Chem Phys 2024; 160:164303. [PMID: 38647303 DOI: 10.1063/5.0205671] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2024] [Accepted: 04/03/2024] [Indexed: 04/25/2024] Open
Abstract
An exciton-phonon (ex-ph) model based on our recently developed block interaction product basis framework is introduced to simulate the distal charge separation (CS) process in aggregated perylenediimide (PDI) trimer incorporating the quantum dynamic method, i.e., the time-dependent density matrix renormalization group. The electronic Hamiltonian in the ex-ph model is represented by nine constructed diabatic states, which include three local excited (LE) states and six charge transfer (CT) states from both the neighboring and distal chromophores. These diabatic states are automatically generated from the direct products of the leading localized neutral or ionic states of each chromophore's reduced density matrix, which are obtained from ab initio quantum chemical calculation of the subsystem consisting of the targeted chromophore and its nearest neighbors, thus considering the interaction of the adjacent environment. In order to quantum-dynamically simulate the distal CS process with massive coupled vibrational modes in molecular aggregates, we used our recently proposed hierarchical mapping approach to renormalize these modes and truncate those vibrational modes that are not effectively coupled with electronic states accordingly. The simulation result demonstrates that the formation of the distal CS process undergoes an intermediate state of adjacent CT, i.e., starts from the LE states, passes through an adjacent CT state to generate the intermediates (∼200 fs), and then formalizes the targeted distal CS via further charge transference (∼1 ps). This finding agrees well with the results observed in the experiment, indicating that our scheme is capable of quantitatively investigating the CS process in a realistic aggregated PDI trimer and can also be potentially applied to exploring CS and other photoinduced processes in larger systems.
Collapse
Affiliation(s)
- Ke Wang
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Yihe Xu
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Xiaoyu Xie
- Qingdao Institute for Theoretical and Computational Sciences, School of Chemistry and Chemical Engineering, Shandong University, Qingdao, Shandong 266237, China
| | - Haibo Ma
- Qingdao Institute for Theoretical and Computational Sciences, School of Chemistry and Chemical Engineering, Shandong University, Qingdao, Shandong 266237, China
| |
Collapse
|
3
|
Fedorov DG. Analysis of Site Energies and Excitonic Couplings: The Role of Symmetry and Polarization. J Phys Chem A 2024; 128:1154-1162. [PMID: 38302431 DOI: 10.1021/acs.jpca.3c06293] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2024]
Abstract
An excitonic coupling model is developed based on an equation-of-motion coupled cluster combined with the fragment molecular orbital method. The effects of polarization and excitonic coupling on the splitting of quasi-degenerate levels in systems containing multiple chromophores are elucidated on dimers of formaldehyde, water, formic acid, hydrogen fluoride, and carbon monoxide. It is shown that the level structure is mainly determined by the mutual polarization of chromophores and to a lesser extent by the excitonic coupling. The role of symmetry in excitonic coupling in dimers is discussed. The excitonic coupling between all residues in the photoactive yellow protein (PDB: 2PHY) is analyzed.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
4
|
Fedorov DG. Site-Specific Ionization Potentials and Electron Affinities in Large Molecular Systems at Coupled Cluster Level. J Phys Chem A 2023; 127:9357-9364. [PMID: 37782030 DOI: 10.1021/acs.jpca.3c04847] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/03/2023]
Abstract
A many-body expansion of ionization potentials and electron affinities is developed based on a combination of the fragment molecular orbital method and equation-of-motion coupled-cluster (EOM-CC). In addition to site-specific values, obtained as one-body properties, pair and triple corrections are added to account for nonlocal EOM-CC contributions of the molecular environment of a chromophore. The developed method is applied to carboxylic acids, alkyl cations, a protein ubiquitin (Protein Data Bank ID 1UBQ), and a nano ribbon of white graphene elucidating the effect of environment on ionization potential and electron affinity.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
5
|
Braunscheidel NM, Abraham V, Mayhall NJ. Generalization of the Tensor Product Selected CI Method for Molecular Excited States. J Phys Chem A 2023; 127:8179-8193. [PMID: 37733948 DOI: 10.1021/acs.jpca.3c03161] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/23/2023]
Abstract
In a recent paper [JCTC, 2020, 16, 6098], we introduced a new approach for accurately approximating full CI ground states in large electronic active-spaces called Tensor Product Selected CI (TPSCI). In TPSCI, a large orbital active space is first partitioned into disjoint sets (clusters) for which the exact, local many-body eigenstates are obtained. Tensor products of these locally correlated many-body states are taken as the basis for the full, global Hilbert space. By folding correlation into the basis states themselves, the low-energy eigenstates become increasingly sparse, creating a more compact selected CI expansion. While we demonstrated that this approach can improve accuracy for a variety of systems, there is even greater potential for applications to excited states, particularly those which have some excited-state character. In this paper, we report on the accuracy of TPSCI for excited states, including a far more efficient implementation in the Julia programming language. In traditional SCI methods that use a Slater determinant basis, accurate excitation energies are obtained only after a linear extrapolation and at a large computational cost. We find that TPSCI with perturbative corrections provides accurate excitation energies for several excited states of various polycyclic aromatic hydrocarbons with respect to the extrapolated result (i.e., near exact result). Further, we use TPSCI to report highly accurate estimates of the lowest 31 eigenstates for a tetracene tetramer system with an active space of 40 electrons in 40 orbitals, giving direct access to the initial bright states and the resulting 18 doubly excited (biexcitonic) states.
Collapse
Affiliation(s)
| | - Vibin Abraham
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Nicholas J Mayhall
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24060, United States
| |
Collapse
|
6
|
Wang YC, Feng S, Kong Y, Huang X, Liang W, Zhao Y. Electronic Couplings for Singlet Fission Processes Based on the Fragment Particle-Hole Densities. J Chem Theory Comput 2023. [PMID: 37296507 DOI: 10.1021/acs.jctc.3c00243] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
A new diabatization scheme is proposed to calculate the electronic couplings for the singlet fission process in multichromophoric systems. In this approach, a robust descriptor that treats single and multiple excitations on an equal footing is adopted to quantify the localization degree of the particle and hole densities of the electronic states. By maximally localizing the particles and holes in terms of predefined molecular fragments, quasi-diabatic states with well-defined characters (locally excited, charge transfer, correlated triplet pair, etc.) can be automatically constructed as the linear combinations of the adiabatic ones, and the electronic couplings can be directly obtained. This approach is very general in that it applies to electronic states with various spin multiplicities and can be combined with various kinds of preliminary electronic structure calculations. Due to the high numerical efficiency, it is able to manipulate more than 100 electronic states in diabatization. The applications to the tetracene dimer and trimer reveal that high-lying multiply excited charge transfer states have significant influences on both the formation and separation of the correlated triplet pair and can even enlarge the coupling for the latter process by 1 order of magnitude.
Collapse
Affiliation(s)
- Yu-Chen Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Shishi Feng
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Yi Kong
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Xunkun Huang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - WanZhen Liang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Yi Zhao
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| |
Collapse
|
7
|
Ma Y, Li Z, Chen X, Ding B, Li N, Lu T, Zhang B, Suo B, Jin Z. Machine-learning assisted scheduling optimization and its application in quantum chemical calculations. J Comput Chem 2023; 44:1174-1188. [PMID: 36648254 DOI: 10.1002/jcc.27075] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Revised: 12/16/2022] [Accepted: 12/23/2022] [Indexed: 01/18/2023]
Abstract
Easy and effective usage of computational resources is crucial for scientific calculations, both from the perspectives of timeliness and economic efficiency. This work proposes a bi-level optimization framework to optimize the computational sequences. Machine-learning (ML) assisted static load-balancing, and different dynamic load-balancing algorithms can be integrated. Consequently, the computational and scheduling engine of the ParaEngine is developed to invoke optimized quantum chemical (QC) calculations. Illustrated benchmark calculations include high-throughput drug suit, solvent model, P38 protein, and SARS-CoV-2 systems. The results show that the usage rate of given computational resources for high throughput and large-scale fragmentation QC calculations can primarily profit, and faster accomplishing computational tasks can be expected when employing high-performance computing (HPC) clusters.
Collapse
Affiliation(s)
- Yingjin Ma
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - ZhiYing Li
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - Xin Chen
- ShenZhen Bay Laboratory, Shenzhen, China
| | - Bowen Ding
- Institute of Chemistry, Chinese Academy of Sciences, Beijing, China
| | - Ning Li
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
- College of Chemistry and Materials Engineering, Wenzhou University, Wen Zhou, China
| | - Teng Lu
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - Baohua Zhang
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - BingBing Suo
- Department of Physics, Northwest University, Xi'an, China
| | - Zhong Jin
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| |
Collapse
|
8
|
Nakata H, Fedorov DG. Analytic Gradient for Time-Dependent Density Functional Theory Combined with the Fragment Molecular Orbital Method. J Chem Theory Comput 2023; 19:1276-1285. [PMID: 36753486 DOI: 10.1021/acs.jctc.2c01177] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023]
Abstract
The analytic energy gradient of energy with respect to nuclear coordinates is derived for the fragment molecular orbital (FMO) method combined with time-dependent density functional theory (TDDFT). The response terms arising from the use of a polarizable embedding are derived. The obtained analytic FMO-TDDFT gradient is shown to be accurate in comparison to both numerical FMO-TDDFT and unfragmented TDDFT gradients, at the level of two- and three-body expansions. The gradients are used for geometry optimizations, molecular dynamics, vibrational calculations, and simulations of IR and Raman spectra of excited states. The developed method is used to optimize the geometry of the ground and excited electronic states of the photoactive yellow protein (PDB: 2PHY).
Collapse
Affiliation(s)
- Hiroya Nakata
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
9
|
Du J, Liao K, Ma J, Li W, Li S. Generalized Energy-Based Fragmentation Approach for the Electronic Emission Spectra of Large Systems. J Chem Theory Comput 2022; 18:7630-7638. [PMID: 36399522 DOI: 10.1021/acs.jctc.2c00911] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Abstract
The excited-state (ES) geometry optimization and electronic emission (fluorescence and phosphorescence) spectra and the ES vibrational spectra of large systems are great challenges in quantum chemistry. In this work, we develop a generalized energy-based fragmentation (GEBF) approach to compute the localized ES structures and vibrational frequencies of large systems. In this approach, the ES energy derivatives (gradients or Hessians) for a localized ES of a large system can be obtained by combining the ES energy derivatives of the corresponding active subsystems (including local excitation center) and the ground-state energy derivatives of inactive subsystems. Two strategies are adopted to overcome two difficulties from state-classification and state-tracking for treating specific ESs. First, for state-classification, we develop an improved density-based spatial clustering applied with noise algorithm with a modified transition orbital projection (TOP) algorithm, which allow a certain ES energy and energy derivatives of the whole system to be calculated with different ES energies and energy derivatives of active subsystems. Furthermore, we also employ the TOP algorithm for tracking the ESs in their geometry optimizations at the time-dependent density functional theory (TDDFT) level. Then, the GEBF approach is applied to investigate the optimized ES geometries or ES vibrational frequencies for two typical systems. Our results show that the cost-effective GEBF approach can accurately reproduce the TDDFT fluorescence spectra of the cytosine derivative and the experimental phosphorescence spectra of the β-cyclodextrin derivative. The GEBF approach is expected to be routinely applied to investigate the electronic emission spectra of very large systems with local chromophores.
Collapse
Affiliation(s)
- Jiahui Du
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P. R. China
| | - Kang Liao
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P. R. China
| | - Jing Ma
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P. R. China
| | - Wei Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P. R. China
| | - Shuhua Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, P. R. China
| |
Collapse
|
10
|
Wang K, Ma J, Ma H. Characterizing the excited states of large photoactive systems by exciton models. J CHIN CHEM SOC-TAIP 2022. [DOI: 10.1002/jccs.202200381] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Ke Wang
- School of Chemistry and Chemical Engineering Nanjing University Nanjing China
| | - Jing Ma
- School of Chemistry and Chemical Engineering Nanjing University Nanjing China
| | - Haibo Ma
- School of Chemistry and Chemical Engineering Nanjing University Nanjing China
| |
Collapse
|