1
|
Fu J, Li M, Rong C, Zhao D, Liu S. Information-theoretic quantities as effective descriptors of electrophilicity and nucleophilicity in density functional theory. J Mol Model 2024; 30:341. [PMID: 39289254 DOI: 10.1007/s00894-024-06116-7] [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: 05/30/2024] [Accepted: 08/19/2024] [Indexed: 09/19/2024]
Abstract
CONTEXT Electrophilicity and nucleophilicity are two vastly important chemical concepts gauging the capability of atoms in molecules to accept and donate the maximal number of electrons. In our earlier studies, we proposed to simultaneously quantify them using the Kullback-Leibler divergence from the information-theoretic approach in density functional theory. However, several issues with this scheme remain to be clarified such as its general validity, predictability, and relationship with other information-theoretic quantities. In this work, we revisit the matter with bigger datasets and deeper theoretical insights. Five information-theoretic quantities including Kullback-Leibler divergence, Hirshfeld charge, Ghost-Berkowitz-Parr entropy, and second and third orders of relative Onicescu information energy are found to be reliable and robust descriptors of electrophilicity and nucleophilicity propensities. Employing these five descriptors, we design a list of new compounds and predict their electrophilicity and nucleophilicity scales. This work should markedly improve our confidence and capability in applying information-theoretic quantities to evaluate electrophilicity and nucleophilicity propensities and henceforth pave the route for more applications of these quantities from information-theoretic approach in density functional theory in the future. METHODS All structures were fully optimized at the M06-2X/6-311 + G(d) level of DFT functional using the Gaussian 16 package (version C01) with integration grids and tight self-consistent-field convergence. The solvent effect was taken into account by using the implicit solvent model (CPCM) in the CH2Cl2 solvent, and all 3D contour surfaces of Fukui function, local temperature, and ITA (information-theoretic approach) quantities were generated by GaussView. The Multiwfn 3.8 program was used to calculate the ITA indexes and atomic charges.
Collapse
Affiliation(s)
- Jia Fu
- College of Chemistry and Chemical Engineering, Hunan Normal University, Changsha, Hunan, 410081, People's Republic of China
| | - Meng Li
- College of Chemistry and Chemical Engineering, Hunan Normal University, Changsha, Hunan, 410081, People's Republic of China
| | - Chunying Rong
- College of Chemistry and Chemical Engineering, Hunan Normal University, Changsha, Hunan, 410081, People's Republic of China.
| | - Dongbo Zhao
- Institute of Biomedical Research, Yunnan University, Kunming, 650500, People's Republic of China.
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, NC, 27599-3420, USA.
- Department of Chemistry, University of North Carolina, Chapel Hill, NC, 27599-3290, USA.
| |
Collapse
|
2
|
Góger S, Karimpour MR, Tkatchenko A. Four-Dimensional Scaling of Dipole Polarizability: From Single-Particle Models to Atoms and Molecules. J Chem Theory Comput 2024; 20:6621-6631. [PMID: 39015013 PMCID: PMC11325554 DOI: 10.1021/acs.jctc.4c00582] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/18/2024]
Abstract
Scaling laws enable the determination of physicochemical properties of molecules and materials as a function of their size, density, number of electrons or other easily accessible descriptors. Such relations can be counterintuitive and nonlinear, and ultimately yield much needed insight into quantum mechanics of many-particle systems. In this work, we show on the basis of single-particle models, multielectron atoms and molecules that the dipole polarizability of quantum systems is generally proportional to the fourth power of a characteristic length, computed from the ground-state wave function. This four-dimensional (4D) scaling is independent of the ratio of bound-to-bound and bound-to-continuum electronic transitions and applies to many-electron atoms when a correlated length metric is used. Finally, this scaling law is applied to predict the polarizability of molecules by electrostatically coupled atoms-in-molecules approach, obtaining approximately 8% absolute and relative accuracy with respect to hybrid density functional theory (DFT) on the QM7-X data set of organic molecules, providing an efficient and scalable model for the anisotropic polarizability tensors of extended (bio)molecules.
Collapse
Affiliation(s)
- Szabolcs Góger
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Mohammad Reza Karimpour
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| |
Collapse
|
3
|
Valatoon S, Alipour M. Reaction rate constant: a theoretical description from local temperature. Phys Chem Chem Phys 2024; 26:14839-14846. [PMID: 38726725 DOI: 10.1039/d4cp01251h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/23/2024]
Abstract
Application of various descriptors based on electron density and its associated quantities to quantify chemical reactivity within the conceptual density functional theory has recently come into spotlight. Among others and particularly relevant to our study, local temperature based on electron density as well as kinetic energy density, as a measure of the kinetic energy of an electron moving in the Kohn-Sham potential of systems, should be mentioned. In this work, we propose to use the local temperature for describing the reaction rate constant, where our main idea originates from the point that the smaller the local temperature at the reaction center, the easier the electron removal, leading to a larger rate constant. On the basis of theoretical considerations, it is proved that the rate constant variations caused by the substituent effects can well be proportional to the local temperature at the reaction center. In order to numerically validate our proposed approach, we have taken the phenol derivatives with the available experimental rate constants of their O-methylation reaction as working models. The reason for this choice is that one of the most versatile approaches for labeling biologically active compounds with the 11C nuclide for positron emission tomography (PET) is methylation by methyl iodide including 11C nuclide, [11C]MeI, where methylation of phenolic oxygen with [11C]MeI is utilized to label some important tracers for PET studies. Our results unveil that the local temperature changes at the reaction center of the aforementioned reaction are reasonably correlated with the rate constant variations. Hopefully, incorporating the proposed correlations between the local temperature and the kinetics data into a computer control algorithm not only provides a simple tool for predicting the rate constant of the O-methylation reaction for other substituted phenols, but also, as a part of the chemical artificial intelligence, the optimum [11C]MeI labeling conditions for a wide variety of phenol derivatives can be controlled.
Collapse
Affiliation(s)
- Saba Valatoon
- Department of Chemistry, School of Science, Shiraz University, Shiraz 71946-84795, Iran.
| | - Mojtaba Alipour
- Department of Chemistry, School of Science, Shiraz University, Shiraz 71946-84795, Iran.
| |
Collapse
|
4
|
Zhao D, Zhao Y, Xu E, Liu W, Ayers PW, Liu S, Chen D. Fragment-Based Deep Learning for Simultaneous Prediction of Polarizabilities and NMR Shieldings of Macromolecules and Their Aggregates. J Chem Theory Comput 2024; 20:2655-2665. [PMID: 38441881 DOI: 10.1021/acs.jctc.3c01415] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/08/2024]
Abstract
Simultaneous prediction of the molecular response properties, such as polarizability and the NMR shielding constant, at a low computational cost is an unresolved issue. We propose to combine a linear-scaling generalized energy-based fragmentation (GEBF) method and deep learning (DL) with both molecular and atomic information-theoretic approach (ITA) quantities as effective descriptors. In GEBF, the total molecular polarizability can be assembled as a linear combination of the corresponding quantities calculated from a set of small embedded subsystems in GEBF. In the new GEBF-DL(ITA) protocol, one can predict subsystem polarizabilities based on the corresponding molecular wave function (thus electron density and ITA quantities) and DL model rather than calculate them from the computationally intensive coupled-perturbed Hartree-Fock or Kohn-Sham equations and finally obtain the total molecular polarizability via a linear combination equation. As a proof-of-concept application, we predict the molecular polarizabilities of large proteins and protein aggregates. GEBF-DL(ITA) is shown to be as accurate enough as GEBF, with mean absolute percentage error <1%. For the largest protein aggregate (>4000 atoms), GEBF-DL(ITA) gains a speedup ratio of 3 compared with GEBF. It is anticipated that when more advanced electronic structure methods are used, this advantage will be more appealing. Moreover, one can also predict the NMR chemical shieldings of proteins with reasonably good accuracy. Overall, the cost-efficient GEBF-DL(ITA) protocol should be a robust theoretical tool for simultaneously predicting polarizabilities and NMR shieldings of large systems.
Collapse
Affiliation(s)
- Dongbo Zhao
- Institute of Biomedical Research, Yunnan University, Kunming, Yunnan 650500, P. R. China
| | - Yilin Zhao
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton ONL8S4M1, Canada
| | - Enhua Xu
- Graduate School of System Informatics, Kobe University, Nada-ku, Kobe, Hyogo 657-8501, Japan
| | - Wenqi Liu
- Institute of Biomedical Research, Yunnan University, Kunming, Yunnan 650500, P. R. China
| | - Paul W Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton ONL8S4M1, Canada
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420, United States
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290, United States
| | - Dahua Chen
- Institute of Biomedical Research, Yunnan University, Kunming, Yunnan 650500, P. R. China
| |
Collapse
|
5
|
He X, Li M, Rong C, Zhao D, Liu W, Ayers PW, Liu S. Some Recent Advances in Density-Based Reactivity Theory. J Phys Chem A 2024; 128:1183-1196. [PMID: 38329898 DOI: 10.1021/acs.jpca.3c07997] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2024]
Abstract
Establishing a chemical reactivity theory in density functional theory (DFT) language has been our intense research interest in the past two decades, exemplified by the determination of steric effect and stereoselectivity, evaluation of electrophilicity and nucleophilicity, identification of strong and weak interactions, and formulation of cooperativity, frustration, and principle of chirality hierarchy. In this Featured Article, we first overview the four density-based frameworks in DFT to appreciate chemical understanding, including conceptual DFT, use of density associated quantities, information-theoretic approach, and orbital-free DFT, and then present a few recent advances of these frameworks as well as new applications from our studies. To that end, we will introduce the relationship among these frameworks, determining the entire spectrum of interactions with Pauli energy derivatives, performing topological analyses with information-theoretic quantities, and extending the density-based frameworks to excited states. Applications to examine physiochemical properties in external electric fields and to evaluate polarizability for proteins and crystals are discussed. A few possible directions for future development are followed, with the special emphasis on its merger with machine learning.
Collapse
Affiliation(s)
- Xin He
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, China
| | - Meng Li
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), Hunan Normal University, Changsha, Hunan 410081, China
| | - Chunying Rong
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), Hunan Normal University, Changsha, Hunan 410081, China
| | - Dongbo Zhao
- Institute of Biomedical Research, Yunnan University, Kunming 650500, China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, China
| | - Paul W Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton ONL8S, Canada
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420, United States
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290, United States
| |
Collapse
|
6
|
Li M, Wan X, Rong C, Zhao D, Liu S. Directionality and additivity effects of molecular acidity and aromaticity for substituted benzoic acids under external electric fields. Phys Chem Chem Phys 2023; 25:27805-27816. [PMID: 37814823 DOI: 10.1039/d3cp02982d] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/11/2023]
Abstract
Our recent study [M. Li et al.Phys. Chem. Chem. Phys., 2023, 25, 2595-2605] unveiled that the impact of an external electric field on molecular acidity and aromaticity for benzoic acid is directional, which can be understood using changes in frontier orbitals and partial charges. However, it is unclear if the effect will disappear when substituting groups are present and whether new patterns of changes will show up. In this work, as a continuation of our efforts to appreciate the impact of external electric fields on physiochemical properties, we find that the directionality effect is still in place for substituted benzoic acid derivatives and that there exists the additivity effect with respect to the number of substituent groups, regardless of the direction of the applied field and the type of substituting groups. We confirm the findings using electron-donating and electron-accepting groups with the electric field applied either parallelly or perpendicularly to the carboxyl group along the benzene ring. The directionality and additivity effects uncovered from this work should enrich the body of our knowledge about the impact of external electric fields on physiochemical properties and could be applicable to other systems and properties as well.
Collapse
Affiliation(s)
- Meng Li
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), Hunan Normal University, Changsha, Hunan 410081, P. R. China.
| | - Xinjie Wan
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), Hunan Normal University, Changsha, Hunan 410081, P. R. China.
| | - Chunying Rong
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), Hunan Normal University, Changsha, Hunan 410081, P. R. China.
| | - Dongbo Zhao
- Institute of Biomedical Research, Yunnan University, Kunming 650500, P. R. China.
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, NC 27599-3420, USA.
- Department of Chemistry, University of North Carolina, Chapel Hill, NC 27599-3290, USA
| |
Collapse
|
7
|
Zhao D, Zhao Y, He X, Li Y, Ayers PW, Liu S. Accurate and Efficient Prediction of Post-Hartree-Fock Polarizabilities of Condensed-Phase Systems. J Chem Theory Comput 2023; 19:6461-6470. [PMID: 37676647 DOI: 10.1021/acs.jctc.3c00646] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/08/2023]
Abstract
To accurately and efficiently predict the molecular response properties (such as polarizability) at post-Hartree-Fock levels for condensed-phase systems under periodic boundary conditions (PBC) is still an unaccomplished and ongoing task. We demonstrate that static isotropic polarizabilities can be cost-effectively predicted at post-Hartree-Fock levels by combining the linear-scaling generalized energy-based fragmentation (GEBF) and information-theoretic approach (ITA) quantities. In PBC-GEBF, the total molecular polarizability of an extended system is obtained as a linear combination of the corresponding quantities of a series of small embedded subsystems of several monomers. Here, we show that in the PBC-GEBF-ITA framework, one can obtain the molecular polarizabilities and establish linear relations to ITA quantities. Once these relations are established for smaller subsystems, one can predict the polarizabilities of larger subsystems directly from the molecular wavefunction (or electron density) via ITA quantities. Alternatively, one can determine the total molecular polarizability via a linear combination equation in PBC-GEBF. We have corroborated that this newly proposed PBC-GEBF-ITA protocol is much more efficient than the original PBC-GEBF approach but is not much less accurate and that this conclusion holds for both many-body perturbation theory and the coupled cluster calculations. Good efficiency and transferability of the PBC-GEBF-ITA protocol are demonstrated for periodic systems with several hundred atoms in a unit cell.
Collapse
Affiliation(s)
- Dongbo Zhao
- Institute of Biomedical Research, Yunnan University, Kunming 650500, P. R. China
| | - Yilin Zhao
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON L8S 4M1, Canada
| | - Xin He
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao 266237, P. R. China
| | - Yunzhi Li
- School of Chemistry and Chemical Engineering, Linyi University, Linyi 276000, P. R. China
| | - Paul W Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON L8S 4M1, Canada
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420, United States
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290, United States
| |
Collapse
|
8
|
Zhao D, He X, Ayers PW, Liu S. Excited-State Polarizabilities: A Combined Density Functional Theory and Information-Theoretic Approach Study. Molecules 2023; 28:molecules28062576. [PMID: 36985548 PMCID: PMC10058485 DOI: 10.3390/molecules28062576] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2023] [Revised: 03/07/2023] [Accepted: 03/10/2023] [Indexed: 03/14/2023] Open
Abstract
Accurate and efficient determination of excited-state polarizabilities (α) is an open problem both experimentally and computationally. Following our previous work, (Phys. Chem. Chem. Phys. 2023, 25, 2131−2141), in which we employed simple ground-state (S0) density-related functions from the information-theoretic approach (ITA) to accurately and efficiently evaluate the macromolecular polarizabilities, in this work we aimed to predict the lowest excited-state (S1) polarizabilities. The philosophy is to use density-based functions to depict excited-state polarizabilities. As a proof-of-principle application, employing 2-(2′-hydroxyphenyl)benzimidazole (HBI), its substituents, and some other commonly used ESIPT (excited-state intramolecular proton transfer) fluorophores as model systems, we verified that either with S0 or S1 densities as an input, ITA quantities can be strongly correlated with the excited-state polarizabilities. When transition densities are considered, both S0 and S1 polarizabilities are in good relationships with some ITA quantities. The transferability of the linear regression model is further verified for a series of molecules with little or no similarity to those molecules in the training set. Furthermore, the excitation energies can be predicted based on multivariant linear regression equations of ITA quantities. This study also found that the nature of both the ground-state and excited-state polarizabilities of these species are due to the spatial delocalization of the electron density.
Collapse
Affiliation(s)
- Dongbo Zhao
- Institute of Biomedical Research, Yunnan University, Kunming 650500, China
| | - Xin He
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao 266237, China
| | - Paul W. Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON L8S 4M1, Canada
- Correspondence: (P.W.A.); (S.L.)
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, NC 27599-3420, USA
- Department of Chemistry, University of North Carolina, Chapel Hill, NC 27599-3290, USA
- Correspondence: (P.W.A.); (S.L.)
| |
Collapse
|