1
|
Mashkovtsev D, Orimoto Y, Aoki Y. Fast and Accurate Calculation of the UV-Vis Spectrum with the Modified Local Excitation Approximation. J Chem Theory Comput 2023; 19:5548-5562. [PMID: 37471461 DOI: 10.1021/acs.jctc.3c00266] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/22/2023]
Abstract
The local excitation approximation (LEA), a method for the calculation of electronic excitations localized in a specific region of a molecule, has been modified with new approaches to enhance the accuracy of the original method. The primary concept behind LEA involves isolating the region of interest as a submolecule from the full molecule using a localization method, followed by calculating electronic excitations solely within this submolecule. In this study, we examined approaches that improve the accuracy in describing the region of interest, particularly its molecular orbital energies. Additionally, the localization method was extended with a new projection technique to accelerate calculations. These approaches were studied in time-dependent density functional theory (TDDFT) calculations applied to four testing systems with a chromophore as the region of interest: two basic linear molecules, acrolein surrounded by 24 water molecules, and a model of a green fluorescent protein. For all studied systems, the results of TDDFT calculations combined with LEA exhibited near-zero error when groups of atoms adjacent to the chromophore were explicitly included in the submolecule. This was achieved with at least a quadratic speedup of the calculation time as a function of the submolecule size.
Collapse
Affiliation(s)
- Denis Mashkovtsev
- Department of Molecular and Material Sciences, Interdisciplinary Graduate School of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
| | - Yuuichi Orimoto
- Department of Material Sciences, Faculty of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
| | - Yuriko Aoki
- Department of Material Sciences, Faculty of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
| |
Collapse
|
2
|
Chen WK, Fang WH, Cui G. Extending multi-layer energy-based fragment method for excited-state calculations of large covalently bonded fragment systems. J Chem Phys 2023; 158:044110. [PMID: 36725521 DOI: 10.1063/5.0129458] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
Recently, we developed a low-scaling Multi-Layer Energy-Based Fragment (MLEBF) method for accurate excited-state calculations and nonadiabatic dynamics simulations of nonbonded fragment systems. In this work, we extend the MLEBF method to treat covalently bonded fragment ones. The main idea is cutting a target system into many fragments according to chemical properties. Fragments with dangling bonds are first saturated by chemical groups; then, saturated fragments, together with the original fragments without dangling bonds, are grouped into different layers. The accurate total energy expression is formulated with the many-body energy expansion theory, in combination with the inclusion-exclusion principle that is used to delete the contribution of chemical groups introduced to saturate dangling bonds. Specifically, in a two-layer MLEBF model, the photochemically active and inert layers are calculated with high-level and efficient electronic structure methods, respectively. Intralayer and interlayer energies can be truncated at the two- or three-body interaction level. Subsequently, through several systems, including neutral and charged covalently bonded fragment systems, we demonstrate that MLEBF can provide accurate ground- and excited-state energies and gradients. Finally, we realize the structure, conical intersection, and path optimizations by combining our MLEBF program with commercial and free packages, e.g., ASE and SciPy. These developments make MLEBF a practical and reliable tool for studying complex photochemical and photophysical processes of large nonbonded and bonded fragment systems.
Collapse
Affiliation(s)
- Wen-Kai Chen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
3
|
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
|
4
|
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
|
5
|
Du J, Ma Y, Ma J, Li S, Li W. Transition orbital projection approach for excited state tracking. J Chem Phys 2022; 156:214104. [DOI: 10.1063/5.0081207] [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
Quantitively comparing the features between different electronic excited states (ESs) is a crucial task in both potential energy surface (PES) studies and excited-state fragmentation approaches. However, it is still a challenging problem in regard to the comparison of complex and highly degenerate systems. Herein, we present a transition orbital projection (TOP) method to calculate the similarity of different ESs based on the configuration vectors of two types of transition densities. It fully considers four significant problems, including phase, hole-particle bijectivity, orbital permutation, and sign of configuration coefficients. TOP state-tracking-based excited-state optimization shows high robustness in several high-symmetric systems, which are difficult to describe with traditional state-tracking approaches. The TOP state-tracking method is expected to be widely applied to the PES of photochemical reactions, ES molecular dynamics to track the diabatic states, and fragmentation approaches for local excitation of large systems.
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, People’s Republic of China
| | - Yixuan 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, People’s Republic of 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, People’s Republic of 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, People’s Republic of 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, People’s Republic of China
| |
Collapse
|
6
|
Wang K, Xie Z, Luo Z, Ma H. Low-Scaling Excited State Calculation Using the Block Interaction Product State. J Phys Chem Lett 2022; 13:462-470. [PMID: 35015548 DOI: 10.1021/acs.jpclett.1c03445] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
We develop an automatic and efficient scheme for the accurate construction of the bases for excitonic models, which can enable "black-box" excited state structure calculations for large molecular systems. These new and optimized bases, which are named the block interaction product state (BIPS), can be expressed as the direct products of the local states for each chromophore. Each chromophore's local states are selected by diagonalization of its reduced density matrix, which is obtained by quantum chemical calculation of the small subsystem composed of the chromophore and its nearest neighbors. We implemented the BIPS framework with fragment-based calculations considering two- and three-body interactions. Test calculations for eight different molecular aggregates demonstrate that this framework provides an accurate description of not only the excitation energies but also the first-order wave function properties (dipole moment and transition dipole moment) of the low-lying excited states at a low-scaling computational cost.
Collapse
Affiliation(s)
- Ke Wang
- School of Chemistry and Chemical Engineering, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University, Nanjing 210023, China
| | - Zhaoxuan Xie
- School of Chemistry and Chemical Engineering, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University, Nanjing 210023, China
| | - Zhen Luo
- School of Chemistry and Chemical Engineering, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University, Nanjing 210023, China
| | - Haibo Ma
- School of Chemistry and Chemical Engineering, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University, Nanjing 210023, China
| |
Collapse
|
7
|
Shen C, Wang X, He X. Fragment-Based Quantum Mechanical Calculation of Excited-State Properties of Fluorescent RNAs. Front Chem 2022; 9:801062. [PMID: 35004616 PMCID: PMC8727457 DOI: 10.3389/fchem.2021.801062] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2021] [Accepted: 11/24/2021] [Indexed: 11/13/2022] Open
Abstract
Fluorescent RNA aptamers have been successfully applied to track and tag RNA in a biological system. However, it is still challenging to predict the excited-state properties of the RNA aptamer–fluorophore complex with the traditional electronic structure methods due to expensive computational costs. In this study, an accurate and efficient fragmentation quantum mechanical (QM) approach of the electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) scheme was applied for calculations of excited-state properties of the RNA aptamer–fluorophore complex. In this method, the excited-state properties were first calculated with one-body fragment quantum mechanics/molecular mechanics (QM/MM) calculation (the excited-state properties of the fluorophore) and then corrected with a series of two-body fragment QM calculations for accounting for the QM effects from the RNA on the excited-state properties of the fluorophore. The performance of the EE-GMFCC on prediction of the absolute excitation energies, the corresponding transition electric dipole moment (TEDM), and atomic forces at both the TD-HF and TD-DFT levels was tested using the Mango-II RNA aptamer system as a model system. The results demonstrate that the calculated excited-state properties by EE-GMFCC are in excellent agreement with the traditional full-system time-dependent ab initio calculations. Moreover, the EE-GMFCC method is capable of providing an accurate prediction of the relative conformational excited-state energies for different configurations of the Mango-II RNA aptamer system extracted from the molecular dynamics (MD) simulations. The fragmentation method further provides a straightforward approach to decompose the excitation energy contribution per ribonucleotide around the fluorophore and then reveals the influence of the local chemical environment on the fluorophore. The applications of EE-GMFCC in calculations of excitation energies for other RNA aptamer–fluorophore complexes demonstrate that the EE-GMFCC method is a general approach for accurate and efficient calculations of excited-state properties of fluorescent RNAs.
Collapse
Affiliation(s)
- Chenfei Shen
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, China
| | - Xianwei Wang
- College of Science, Zhejiang University of Technology, Hangzhou, China
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, China.,New York University-East China Normal University Center for Computational Chemistry at New York University Shanghai, Shanghai, China
| |
Collapse
|
8
|
Generalized Energy-Based Fragmentation Approach for the Accurate Binding Energies and Raman Spectra of Methane Hydrate Clusters. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2111256] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
9
|
Li W, Ma H, Li S, Ma J. Computational and data driven molecular material design assisted by low scaling quantum mechanics calculations and machine learning. Chem Sci 2021; 12:14987-15006. [PMID: 34909141 PMCID: PMC8612375 DOI: 10.1039/d1sc02574k] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2021] [Accepted: 10/12/2021] [Indexed: 12/11/2022] Open
Abstract
Electronic structure methods based on quantum mechanics (QM) are widely employed in the computational predictions of the molecular properties and optoelectronic properties of molecular materials. The computational costs of these QM methods, ranging from density functional theory (DFT) or time-dependent DFT (TDDFT) to wave-function theory (WFT), usually increase sharply with the system size, causing the curse of dimensionality and hindering the QM calculations for large sized systems such as long polymer oligomers and complex molecular aggregates. In such cases, in recent years low scaling QM methods and machine learning (ML) techniques have been adopted to reduce the computational costs and thus assist computational and data driven molecular material design. In this review, we illustrated low scaling ground-state and excited-state QM approaches and their applications to long oligomers, self-assembled supramolecular complexes, stimuli-responsive materials, mechanically interlocked molecules, and excited state processes in molecular aggregates. Variable electrostatic parameters were also introduced in the modified force fields with the polarization model. On the basis of QM computational or experimental datasets, several ML algorithms, including explainable models, deep learning, and on-line learning methods, have been employed to predict the molecular energies, forces, electronic structure properties, and optical or electrical properties of materials. It can be conceived that low scaling algorithms with periodic boundary conditions are expected to be further applicable to functional materials, perhaps in combination with machine learning to fast predict the lattice energy, crystal structures, and spectroscopic properties of periodic functional materials.
Collapse
Affiliation(s)
- 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 China
| | - Haibo 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 China
- Jiangsu Key Laboratory of Advanced Organic Materials, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University Nanjing 210023 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 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 China
- Jiangsu Key Laboratory of Advanced Organic Materials, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University Nanjing 210023 China
| |
Collapse
|
10
|
Wang Z, Liu W. iOI: An Iterative Orbital Interaction Approach for Solving the Self-Consistent Field Problem. J Chem Theory Comput 2021; 17:4831-4845. [PMID: 34240856 DOI: 10.1021/acs.jctc.1c00445] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
An iterative orbital interaction (iOI) approach is proposed to solve, in a bottom-up fashion, the self-consistent field problem in quantum chemistry. While it belongs grossly to the family of fragment-based quantum chemical methods, iOI is distinctive in that (1) it divides and conquers not only the energy but also the wave function and that (2) the subsystem sizes are automatically determined by successively merging neighboring small subsystems until they are just enough for converging the wave function to a given accuracy. Orthonormal occupied and virtual localized molecular orbitals are obtained in a natural manner, which can be used for all post-SCF purposes.
Collapse
Affiliation(s)
- Zikuan Wang
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, P. R. China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, P. R. China
| |
Collapse
|
11
|
Li W, Dong H, Ma J, Li S. Structures and Spectroscopic Properties of Large Molecules and Condensed-Phase Systems Predicted by Generalized Energy-Based Fragmentation Approach. Acc Chem Res 2021; 54:169-181. [PMID: 33350806 DOI: 10.1021/acs.accounts.0c00580] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
ConspectusThe structures and spectroscopic properties of molecules and condensed-phase systems are usually experimentally characterized by X-ray, infrared (IR), Raman, nuclear magnetic resonance (NMR), and electronic absorption/emission spectra. Quantum mechanics (QM) calculations are critical in quantitatively understanding the relationship between the structure and physicochemical properties of various chemical systems. However, it is very challenging to apply traditional QM methods to large molecules and condensed-phase systems with large unit cells due to their steep computational scaling with the system size. To overcome this difficulty, theoretical chemists have developed various linear (or low) scaling QM methods, among which energy-based fragmentation methods have achieved great success for large molecules or clusters. One of the most popular energy-based fragmentation methods is the generalized energy-based fragmentation (GEBF) approach developed by us.In this approach, the ground-state energy of a large molecule can be evaluated from the ground-state energies of a series of embedded subsystems. In this Account, we focus on the recent developments and applicability of the GEBF approach for the structures and spectroscopic properties of complicated large molecules and condensed-phase systems. With new fragmentation schemes, the GEBF approach can now describe ionic liquid clusters and metal-containing supramolecular systems accurately and can provide accurate binding energies for host-guest complexes. In addition, the GEBF approach is now available for describing the localized excited states of large systems including a chromophore. More importantly, the GEBF approach under periodic boundary conditions (PBC-GEBF) has been developed to deal with periodic molecular crystals and liquids. Then, the ground-state energy (or property) per unit cell of a periodic condensed phase system can be predicted with QM calculations on nonperiodic embedded subsystems. This feature enables accurate electron correlation calculations on molecular crystals and liquids to be feasible on ordinary workstations. The PBC-GEBF approach has been applied to predict the crystal structures, lattice energies, and spectroscopic properties of some typical molecular crystals and solutions. By combining the GEBF method and machine learning (ML) method, a GEBF-ML force field has been developed for long normal alkanes, and the IR spectra of long alkanes can be obtained from the GEBF-ML molecular dynamics (MD) simulations. The GEBF and its periodic variant are expected to play increasingly important roles in investigating real-life chemical systems of broad interests at the ab initio levels.
Collapse
Affiliation(s)
- Wei Li
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Hao Dong
- Kuang Yaming Honors School and Institute for Brain Sciences, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Jing Ma
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Shuhua Li
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| |
Collapse
|
12
|
Kim Y, Bui Y, Tazhigulov RN, Bravaya KB, Slipchenko LV. Effective Fragment Potentials for Flexible Molecules: Transferability of Parameters and Amino Acid Database. J Chem Theory Comput 2020; 16:7735-7747. [PMID: 33236635 DOI: 10.1021/acs.jctc.0c00758] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
An accurate but efficient description of noncovalent interactions is a key to predictive modeling of biological and materials systems. The effective fragment potential (EFP) is an ab initio-based force field that provides a physically meaningful decomposition of noncovalent interactions of a molecular system into Coulomb, polarization, dispersion, and exchange-repulsion components. An EFP simulation protocol consists of two steps, preparing parameters for molecular fragments by a series of ab initio calculations on each individual fragment, and calculation of interaction energy and properties of a total molecular system based on the prepared parameters. As the fragment parameters (distributed multipoles, polarizabilities, localized wave function, etc.) depend on a fragment geometry, straightforward application of the EFP method requires recomputing parameters of each fragment if its geometry changes, for example, during thermal fluctuations of a molecular system. Thus, recomputing fragment parameters can easily become both computational and human bottlenecks and lead to a loss of efficiency of a simulation protocol. An alternative approach, in which fragment parameters are adjusted to different fragment geometries, referred to as "flexible EFP", is explored here. The parameter adjustment is based on translations and rotations of local coordinate frames associated with fragment atoms. The protocol is validated on extensive benchmark of amino acid dimers extracted from molecular dynamics snapshots of a cryptochrome protein. A parameter database for standard amino acids is developed to automate flexible EFP simulations in proteins. To demonstrate applicability of flexible EFP in large-scale protein simulations, binding energies and vertical electron ionization and electron attachment energies of a lumiflavin chromophore of the cryptochrome protein are computed. The results obtained with flexible EFP are in a close agreement with the standard EFP procedure but provide a significant reduction in computational cost.
Collapse
Affiliation(s)
- Yongbin Kim
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Yen Bui
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Ruslan N Tazhigulov
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Ksenia B Bravaya
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| | - Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| |
Collapse
|
13
|
Jin X, Glover WJ, He X. Fragment Quantum Mechanical Method for Excited States of Proteins: Development and Application to the Green Fluorescent Protein. J Chem Theory Comput 2020; 16:5174-5188. [PMID: 32551640 DOI: 10.1021/acs.jctc.9b00980] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Understanding the excited-state properties of luminescent biomolecules is of central importance to their biophysical applications. In this study, we develop the Electrostatically Embedded Generalized Molecular Fractionation with Conjugate Caps (EE-GMFCC) method for quantitatively characterizing properties of covalently bonded systems with localized excitations (i.e., involving a single chromophore), such as fluorescent proteins. The excitation energy, transition dipole moment, and oscillator strength of wild-type Green Fluorescent Protein (wt-GFP) calculated by EE-GMFCC are found to be in excellent agreement with full system time-dependent density functional theory results. We also applied the Polarized Protein-Specific Charge model to wt-GFP, and found that electronic polarization of the protein is critical in stabilizing hydrogen bonding interactions in wt-GFP, which influences its absorption spectrum. The predicted absorption spectra of wt-GFP in the A and B states qualitatively agree with experiment. The fragmentation approach further allows a straightforward per residue decomposition of the excitation which reveals the influence of the protein environment on the absorption spectra of wt-GFP A and B states. Our results demonstrate that the EE-GMFCC method is both accurate and efficient for excited-state property calculations on proteins.
Collapse
Affiliation(s)
- Xinsheng Jin
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - William J Glover
- NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China.,Department of Chemistry, New York University, New York, New York 10003, United States
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
14
|
Chen WK, Zhang Y, Jiang B, Fang WH, Cui G. Efficient Construction of Excited-State Hessian Matrices with Machine Learning Accelerated Multilayer Energy-Based Fragment Method. J Phys Chem A 2020; 124:5684-5695. [DOI: 10.1021/acs.jpca.0c04117] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Wen-Kai Chen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Yaolong Zhang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
15
|
Zhang W, Liu J, Jin X, Gu X, Zeng XC, He X, Li H. Quantitative Prediction of Aggregation‐Induced Emission: A Full Quantum Mechanical Approach to the Optical Spectra. Angew Chem Int Ed Engl 2020; 59:11550-11555. [DOI: 10.1002/anie.202003326] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2020] [Indexed: 12/23/2022]
Affiliation(s)
- Wei Zhang
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering Beijing University of Chemical Technology Beijing 100029 China
- Department of Chemistry University of Nebraska Lincoln NE 68588 USA
| | - Jinfeng Liu
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development School of Chemistry and Molecular Engineering East China Normal University Shanghai 200062 China
- Department of Basic Medicine and Clinical Pharmacy China Pharmaceutical University Nanjing 210009 China
| | - Xinsheng Jin
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development School of Chemistry and Molecular Engineering East China Normal University Shanghai 200062 China
| | - Xinggui Gu
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering Beijing University of Chemical Technology Beijing 100029 China
| | - Xiao Cheng Zeng
- Department of Chemistry University of Nebraska Lincoln NE 68588 USA
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development School of Chemistry and Molecular Engineering East China Normal University Shanghai 200062 China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai Shanghai 200062 China
| | - Hui Li
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering Beijing University of Chemical Technology Beijing 100029 China
| |
Collapse
|
16
|
Zhang W, Liu J, Jin X, Gu X, Zeng XC, He X, Li H. Quantitative Prediction of Aggregation‐Induced Emission: A Full Quantum Mechanical Approach to the Optical Spectra. Angew Chem Int Ed Engl 2020. [DOI: 10.1002/ange.202003326] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Wei Zhang
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering Beijing University of Chemical Technology Beijing 100029 China
- Department of Chemistry University of Nebraska Lincoln NE 68588 USA
| | - Jinfeng Liu
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development School of Chemistry and Molecular Engineering East China Normal University Shanghai 200062 China
- Department of Basic Medicine and Clinical Pharmacy China Pharmaceutical University Nanjing 210009 China
| | - Xinsheng Jin
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development School of Chemistry and Molecular Engineering East China Normal University Shanghai 200062 China
| | - Xinggui Gu
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering Beijing University of Chemical Technology Beijing 100029 China
| | - Xiao Cheng Zeng
- Department of Chemistry University of Nebraska Lincoln NE 68588 USA
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development School of Chemistry and Molecular Engineering East China Normal University Shanghai 200062 China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai Shanghai 200062 China
| | - Hui Li
- Beijing Advanced Innovation Center for Soft Matter Science and Engineering Beijing University of Chemical Technology Beijing 100029 China
| |
Collapse
|
17
|
Chen WK, Fang WH, Cui G. Integrating Machine Learning with the Multilayer Energy-Based Fragment Method for Excited States of Large Systems. J Phys Chem Lett 2019; 10:7836-7841. [PMID: 31786927 DOI: 10.1021/acs.jpclett.9b03113] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
In this work we have combined machine learning techniques with our recently developed multilayer energy-based fragment method for studying excited states of large systems. The photochemically active and inert regions are separately treated with the complete active space self-consistent field method and the trained models. This method is demonstrated to provide accurate energies and gradients leading to essentially the same excited-state potential energy surfaces and nonadiabatic dynamics compared with full ab initio results. Furthermore, in conjunction with the use of machine learning models, this method is highly parallel and exhibits low-scaling computational cost. Finally, the present work could encourage the marriage of machine learning with fragment-based electronic structure methods to explore photochemistry of large systems.
Collapse
Affiliation(s)
- Wen-Kai Chen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry , Beijing Normal University , Beijing 100875 , People's Republic of China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry , Beijing Normal University , Beijing 100875 , People's Republic of China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry , Beijing Normal University , Beijing 100875 , People's Republic of China
| |
Collapse
|
18
|
Chen WK, Fang WH, Cui G. A multi-layer energy-based fragment method for excited states and nonadiabatic dynamics. Phys Chem Chem Phys 2019; 21:22695-22699. [PMID: 31595910 DOI: 10.1039/c9cp04842a] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We developed a multi-layer energy-based fragment (MLEBF) method within the many-body energy expansion framework. It supplies accurate energies and gradients, and accurately reproduces excited-state topological structures. Moreover, MLEBF-based nonadiabatic dynamics simulations give nearly the same results compared with full ab initio ones. The present work could stimulate developing energy-based fragment methods for photochemistry of large systems.
Collapse
Affiliation(s)
- Wen-Kai Chen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China.
| | | | | |
Collapse
|
19
|
Izsák R. Single‐reference coupled cluster methods for computing excitation energies in large molecules: The efficiency and accuracy of approximations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1445] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Róbert Izsák
- Max‐Planck‐Institut für Kohlenforschung Mülheim an der Ruhr Germany
| |
Collapse
|
20
|
Liu J, Sun H, Glover WJ, He X. Prediction of Excited-State Properties of Oligoacene Crystals Using Fragment-Based Quantum Mechanical Method. J Phys Chem A 2019; 123:5407-5417. [PMID: 31187994 DOI: 10.1021/acs.jpca.8b12552] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
A fundamental understanding of the excited-state properties of molecular crystals is of central importance for their optoelectronics applications. In this study, we developed the electrostatically embedded generalized molecular fractionation (EE-GMF) method for the quantitative characterization of the excited-state properties of locally excited molecular clusters. The accuracy of the EE-GMF method is systematically assessed for oligoacene crystals. Our result demonstrates that the EE-GMF method is capable of providing the lowest vertical singlet (S1) and triplet excitation energies (T1), in excellent agreement with the full-system quantum mechanical calculations. Using this method, we also investigated the performance of different density functionals in predicting the excited-state properties of the oligoacene crystals. Among the 13 tested functionals, B3LYP and MN15 give the two lowest overall mean unsigned errors with reference to the experimental S1 and T1 excitation energies. The EE-GMF approach can be readily utilized for studying the excited-state properties of large-scale organic solids at diverse ab initio levels.
Collapse
Affiliation(s)
- Jinfeng Liu
- Department of Basic Medicine and Clinical Pharmacy , China Pharmaceutical University , Nanjing 210009 , China
| | | | - William J Glover
- NYU Shanghai , Shanghai 200122 , China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai 200062 , China.,Department of Chemistry , New York University , New York , New York 10003 , United States
| | - Xiao He
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai 200062 , China
| |
Collapse
|
21
|
Fu F, Liao K, Ma J, Cheng Z, Zheng D, Gao L, Liu C, Li S, Li W. How intermolecular interactions influence electronic absorption spectra: insights from the molecular packing of uracil in condensed phases. Phys Chem Chem Phys 2019; 21:4072-4081. [DOI: 10.1039/c8cp06152a] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Intermolecular interactions in terms of molecular packing are crucial for the investigation of the absorption spectra of uracil in different environments.
Collapse
Affiliation(s)
- Fangjia Fu
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Kang Liao
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Jing Ma
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Zheng Cheng
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Dong Zheng
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Liuzhou Gao
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Chungen Liu
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Shuhua Li
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Wei Li
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| |
Collapse
|
22
|
Dutta AK, Saitow M, Riplinger C, Neese F, Izsák R. A near-linear scaling equation of motion coupled cluster method for ionized states. J Chem Phys 2018; 148:244101. [DOI: 10.1063/1.5029470] [Citation(s) in RCA: 40] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Achintya Kumar Dutta
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstr. 34-36, D-45470 Mülheim an der Ruhr, Germany
| | - Masaaki Saitow
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstr. 34-36, D-45470 Mülheim an der Ruhr, Germany
| | | | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Róbert Izsák
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
23
|
Li Y, Yuan D, Wang Q, Li W, Li S. Accurate prediction of the structure and vibrational spectra of ionic liquid clusters with the generalized energy-based fragmentation approach: critical role of ion-pair-based fragmentation. Phys Chem Chem Phys 2018; 20:13547-13557. [DOI: 10.1039/c8cp00513c] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The GEBF method with the ion-pair-based fragmentation has been developed to facilitate ab initio calculations of general ionic liquid clusters.
Collapse
Affiliation(s)
- Yunzhi Li
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Dandan Yuan
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Qingchun Wang
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Wei Li
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Shuhua Li
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| |
Collapse
|
24
|
Zhao D, Song R, Li W, Ma J, Dong H, Li S. Accurate Prediction of NMR Chemical Shifts in Macromolecular and Condensed-Phase Systems with the Generalized Energy-Based Fragmentation Method. J Chem Theory Comput 2017; 13:5231-5239. [DOI: 10.1021/acs.jctc.7b00380] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Dongbo Zhao
- 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, People’s Republic of China
- Kuang
Yaming Honors School, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Ruiheng Song
- Kuang
Yaming Honors School, Nanjing University, Nanjing 210023, People’s Republic of 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, People’s Republic of 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, People’s Republic of China
| | - Hao Dong
- Kuang
Yaming Honors School, Nanjing University, Nanjing 210023, People’s Republic of 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, People’s Republic of China
| |
Collapse
|
25
|
Zhang L, Li W, Fang T, Li S. Accurate Relative Energies and Binding Energies of Large Ice–Liquid Water Clusters and Periodic Structures. J Phys Chem A 2017; 121:4030-4038. [DOI: 10.1021/acs.jpca.7b03376] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Lei Zhang
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Wei Li
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Tao Fang
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Shuhua Li
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| |
Collapse
|