1
|
Lu T. A comprehensive electron wavefunction analysis toolbox for chemists, Multiwfn. J Chem Phys 2024; 161:082503. [PMID: 39189657 DOI: 10.1063/5.0216272] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2024] [Accepted: 08/07/2024] [Indexed: 08/28/2024] Open
Abstract
Analysis of electron wavefunction is a key component of quantum chemistry investigations and is indispensable for the practical research of many chemical problems. After more than ten years of active development, the wavefunction analysis program Multiwfn has accumulated very rich functions, and its application scope has covered numerous aspects of theoretical chemical research, including charge distribution, chemical bond, electron localization and delocalization, aromaticity, intramolecular and intermolecular interactions, electronic excitation, and response property. This article systematically introduces the features and functions of the latest version of Multiwfn and provides many representative examples. Through this article, readers will be able to fully understand the characteristics and recognize the unique value of Multiwfn. The source code and precompiled executable files of Multiwfn, as well as the manual containing a detailed introduction to theoretical backgrounds and very rich tutorials, can all be downloaded for free from the Multiwfn website (http://sobereva.com/multiwfn).
Collapse
Affiliation(s)
- Tian Lu
- Beijing Kein Research Center for Natural Sciences, Beijing 100024, People's Republic of China
| |
Collapse
|
2
|
Kenouche S, Bachir N, Bouchal W, Martínez-Araya JI. Aromaticity of six-membered nitro energetic compounds through molecular electrostatic potential, magnetic, electronic delocalization and reactivity-based indices. J Mol Graph Model 2024; 129:108728. [PMID: 38412811 DOI: 10.1016/j.jmgm.2024.108728] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2023] [Revised: 02/06/2024] [Accepted: 02/13/2024] [Indexed: 02/29/2024]
Abstract
The electron density depletion associated with π-hole at the ring center typical of energetic compounds was clearly revealed by the molecular electrostatic potential (ESP). In addition, the spatial arrangement of NO2 groups appears to affect the ESP value in the ring center, and therefore sensitivity to detonation. Indeed, for monocyclic nitrobenzene compounds with the same number of NO2 groups, the ESP value in the ring center decreases as the NO2 groups are more closely spaced. As expected, the central rings become less aromatic as NO2 groups are added. The MCI, PDI, PLR, NICSzz(1), FLU indices are all strongly correlated with the ESP values observed in the ring center of the set of nitrobenzenes. Aromaticity indices based on electron delocalization criteria appear to be very sensitive to small variations in aromaticity. Among magnetic-based indices, only NICSzz(1) is capable to predict small changes in aromaticity. The PLR index derived from conceptual DFT is quite relevant for predicting small variations in aromaticity. According to our results, the most suitable aromaticity index is not based on a single criterion, and that selecting it is more subtle. Therefore, it is important to combine information from several criteria to obtain a more complete description of the aromaticity of the studied compounds.
Collapse
Affiliation(s)
- Samir Kenouche
- Group of Modeling of Chemical Systems using Quantum Calculations, Applied Chemistry Laboratory (LCA). University M. Khider of Biskra, 07000 Biskra, Algeria
| | - Nassima Bachir
- Group of Modeling of Chemical Systems using Quantum Calculations, Applied Chemistry Laboratory (LCA). University M. Khider of Biskra, 07000 Biskra, Algeria
| | - Wissam Bouchal
- Molecular Chemistry and Environment Laboratory, University of Mohammed Khider of Biskra, BP 145 RP, Biskra 07000, Algeria
| | - Jorge I Martínez-Araya
- Departamento de Ciencias Químicas, Facultad de Ciencias Exactas, Universidad Andres Bello (UNAB), Av. República 275, 8370146 Santiago, Chile; Centro de Química Teórica y Computacional (CQT&C). Facultad de Ciencias Exactas, Santiago, Chile.
| |
Collapse
|
3
|
Wang B, Geerlings P, Liu S, De Proft F. Extending the Scope of Conceptual Density Functional Theory with Second Order Analytical Methodologies. J Chem Theory Comput 2024. [PMID: 38310523 DOI: 10.1021/acs.jctc.3c01184] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2024]
Abstract
In the context of the growing impact of conceptual density functional theory (DFT) as one of the most successful chemical reactivity theories, response functions up to second order have now been widely applied; in recent years, among others, particular attention has been focused on the linear response function and also extensions to higher order have been put forward. As the larger part of these studies have been carried using a finite difference approach to compute these concepts, we now embarked on (an extension of) an analytical approach to conceptual DFT. With the ultimate aim of providing a complete set of analytically computable second order properties, including the softness and hardness kernels, the hardness as the simplest second order response function is scrutinized again with numerical results highlighting the difference in nature between the analytical hardness (referred to as hardness condition) and the Parr-Pearson absolute chemical hardness. The hardness condition is investigated for its capability to gauge the (de)localization error of density functional approximations (DFAs). The analytical Fukui function, besides overcoming the difficulties in the finite difference approach in treating negatively charged systems, also showcases the errors of deviating from the straight-line behavior using fractional occupation number calculations. Subsequently, the softness kernel and its atom-condensed inverse, the hardness matrix, are accessed through the Berkowitz-Parr relation. Revisiting the softness kernel confirms and extends previous discussions on how Kohn's Nearsightedness of Electronic Matter principle can be retrieved and identified as the physicist's version of the chemist's "transferability of functional groups" concept. The accurate, analytical hardness matrix evaluation on the other hand provides further support for the basics of Nalewajski's charge sensitivity analysis. Based on Parr and Liu's functional expansion of the energy functional, a new energy decomposition is introduced with an order of magnitude analysis of the different terms for a series of simple molecules both at their equilibrium geometry and upon variation in bond length and dihedral angle. Finally, for the first time, the perturbation expansion of the energy functional is studied numerically up to second order now that all response functions and integration techniques are at hand. The perturbation expansion energies are in excellent agreement with those obtained directly from DFA calculations giving confidence in the convergence of the perturbation series and its use in judging the importance of the different terms in reactivity investigations.
Collapse
Affiliation(s)
- Bin Wang
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, B-1050 Brussels, Belgium
| | - Paul Geerlings
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, B-1050 Brussels, Belgium
| | - 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
| | - Frank De Proft
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, B-1050 Brussels, Belgium
| |
Collapse
|
4
|
Kozáková S, Alharzali N, Černušák I. Cyclo[ n]carbons and catenanes from different perspectives: disentangling the molecular thread. Phys Chem Chem Phys 2023; 25:29386-29403. [PMID: 37901943 DOI: 10.1039/d3cp03887d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/31/2023]
Abstract
All-carbon atomic rings, cyclo[n]carbons, have recently attracted vivid attention of experimentalists and theoreticians. Among them, cyclo[18]carbon is the most studied system. In this paper, we summarize and review various properties of cyclo[n]carbons, emphasising the aspects of their aromaticity/antiaromaticity. In the first part, the trends in bonding patterns and selected aromaticity indices with the increasing size of the rings are discussed. In the second part we explore the properties of catenane models based on interlocked cyclo[18]carbon rings from different perspectives and investigate their behaviour under the action of external force using computational experiments.
Collapse
Affiliation(s)
- Silvia Kozáková
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava, Ilkovičova 6, 84215 Bratislava, Slovakia.
| | - Nissrin Alharzali
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava, Ilkovičova 6, 84215 Bratislava, Slovakia.
| | - Ivan Černušák
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava, Ilkovičova 6, 84215 Bratislava, Slovakia.
| |
Collapse
|
5
|
Martínez JP, Trzaskowski B. Electrophilicity of Hoveyda-Grubbs Olefin Metathesis Catalysts as the Driving Force that Controls Initiation Rates. Chemphyschem 2022; 23:e202200580. [PMID: 36062870 DOI: 10.1002/cphc.202200580] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2022] [Revised: 09/01/2022] [Indexed: 01/05/2023]
Abstract
The dissociative mechanism of initiation for a series of Hoveyda-Grubbs type metathesis catalysts modified at the para and meta positions in the isopropoxybenzylidene ligand is investigated by means of DFT calculations. The electron donating/withdrawing capacity of the ligand was screened through the incorporation of various substituents such as halogens, nitro, alkoxides, ketones, esters, amines, and amides. Variations in structural parameters, energy barriers for the Ru-O bond dissociation, and Ru-O bond strength were examined as a function of the Hammett constant. It was found that electronic properties of the catalysts such as chemical potential, hardness, and electrophilicity correlate linearly with the dissociative energy barriers. These findings enable a systematic rationalization and prediction of rate of precatalyst initiation through the calculation of only the HOMO-LUMO gap of catalysts, as the faster the initiation, the more electrophilic the catalyst.
Collapse
|
6
|
Apebende CG, idante PS, Magu TO, Asogwa FC, Onyebuenyi IB, Unimuke TO, Gber TE. Density Functional Theory Study of the Influence of Activating and Deactivating Groups on Naphthalene. RESULTS IN CHEMISTRY 2022. [DOI: 10.1016/j.rechem.2022.100669] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
|
7
|
Báez-Grez R, Pino-Rios R. Evaluation of Slight Changes in Aromaticity through Electronic and Density Functional Reactivity Theory-Based Descriptors. ACS OMEGA 2022; 7:21939-21945. [PMID: 35785290 PMCID: PMC9245093 DOI: 10.1021/acsomega.2c02291] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/12/2022] [Accepted: 05/13/2022] [Indexed: 06/15/2023]
Abstract
Aromaticity is a useful tool to rationalize the structure, stability, and reactivity in several compounds. Although aromaticity is not directly an observable, it is well accepted that electronic delocalization around the molecular ring is a key stabilizing feature of aromatic compounds. This contribution presents a systematic evaluation of the capability of delocalization and reactivity criteria to describe aromaticity in a set of fluorinated benzenes. The aromaticity indices are compared with quantities obtained from the magnetic criteria of aromaticity, i.e., the strength of the ring current induced by an external magnetic field and the popular NICS zz (1) index. In this evaluation, the indices based on delocalization criteria used are aromatic fluctuation index (FLU), para-delocalization index (PDI), PDIπ, and the multicenter delocalization index (MCI). In addition, indices based on the bifurcation values of scalar functions are derived from electron density such as electron localization function (the π contribution, ELFπ) and the π contribution of the localized orbital locator (LOLπ). Furthermore, reactivity indices based on chemical reactivity and the information-theoretic (reactivity) approach are para-linear response (PLR), Shannon entropy, Fisher information, and Ghosh-Berkowitz-Parr (GBP) entropy. The results obtained show that the delocalization-based indicators present a high sensitivity to slight changes in aromaticity and that the reactivity criterion can be considered as a complementary tool for the study of this phenomenon, even when these changes are minimal. These results encourage the use of multiple indicators for a complete understanding of aromaticity in various chemical compounds.
Collapse
Affiliation(s)
- Rodrigo Báez-Grez
- Computational
and Theoretical Chemistry Group, Departamento de Ciencias Químicas,
Facultad de Ciencias Exactas, Universidad
Andres Bello, República 275, Santiago 8370146, Chile
| | - Ricardo Pino-Rios
- Instituto
de Ciencias Químicas Aplicadas, Facultad de Ingeniería, Universidad Autónoma de Chile, 7500912 Santiago, Chile
| |
Collapse
|
8
|
Eno EA, Louis H, Unimuke TO, Egemonye TC, Adalikwu SA, Agwupuye JA, Odey DO, Abu AS, Eko IJ, Ifeatu CE, Ntui TN. Synthesis, characterization, and theoretical investigation of 4-chloro-6(phenylamino)-1,3,5-triazin-2-yl)asmino-4-(2,4-dichlorophenyl)thiazol-5-yl-diazenyl)phenyl as potential SARS-CoV-2 agent. PHYSICAL SCIENCES REVIEWS 2022. [DOI: 10.1515/psr-2021-0161] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Abstract
The synthesis of 4-chloro-6(phenylamino)-1,3,5-triazin-2-yl)amino-4-(2,4 dichlorophenyl)thiazol-5-yl-diazenyl)phenyl is reported in this work with a detailed structural and molecular docking study on two SARS-COV-2 proteins: 3TNT and 6LU7. The studied compound has been synthesized by the condensation of cyanuric chloride with aniline and characterized with various spectroscopic techniques. The experimentally obtained spectroscopic data has been compared with theoretical calculated results achieved using high-level density functional theory (DFT) method. Stability, nature of bonding, and reactivity of the studied compound was evaluated at DFT/B3LYP/6-31 + (d) level of theory. Hyper-conjugative interaction persisting within the molecules which accounts for the bio-activity of the compound was evaluated from natural bond orbital (NBO) analysis. Adsorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) properties of the experimentally synthesized compound was studied to evaluate the pharmacological as well as in silico molecular docking against SARS-CoV-2 receptors. The molecular docking result revealed that the investigated compound exhibited binding affinity of −9.3 and −8.8 for protein 3TNT and 6LU7 respectively. In conclusion, protein 3TNT with the best binding affinity for the ligand is the most suitable for treatment of SARS-CoV-2.
Collapse
Affiliation(s)
- Ededet A. Eno
- Computational and Bio-Simulation Research Group , University of Calabar, Calabar , Nigeria
- Department of Pure and Applied Chemistry, Faculty of Physical Sciences , University of Calabar , Calabar , Nigeria
| | - Hitler Louis
- Computational and Bio-Simulation Research Group , University of Calabar, Calabar , Nigeria
- Department of Pure and Applied Chemistry, Faculty of Physical Sciences , University of Calabar , Calabar , Nigeria
| | - Tomsmith O. Unimuke
- Computational and Bio-Simulation Research Group , University of Calabar, Calabar , Nigeria
- Department of Pure and Applied Chemistry, Faculty of Physical Sciences , University of Calabar , Calabar , Nigeria
| | - ThankGod C. Egemonye
- Computational and Bio-Simulation Research Group , University of Calabar, Calabar , Nigeria
- Department of Pure and Applied Chemistry, Faculty of Physical Sciences , University of Calabar , Calabar , Nigeria
| | - Stephen A. Adalikwu
- Computational and Bio-Simulation Research Group , University of Calabar, Calabar , Nigeria
| | - John A. Agwupuye
- Computational and Bio-Simulation Research Group , University of Calabar, Calabar , Nigeria
- Department of Pure and Applied Chemistry, Faculty of Physical Sciences , University of Calabar , Calabar , Nigeria
| | - Diana O. Odey
- Computational and Bio-Simulation Research Group , University of Calabar, Calabar , Nigeria
- Department of Biochemistry, Faculty of Physical Sciences , Cross River University of Technology , Calabar , Nigeria
| | - Abu Solomon Abu
- Computational and Bio-Simulation Research Group , University of Calabar, Calabar , Nigeria
- Department of Marine Biology, Faculty of Biology Sciences , University of Calabar , Calabar , Nigeria
| | - Ishegbe J. Eko
- Department of Polymer and Textile Engineering , Ahmadu Bello University Zaria , Kaduna , Nigeria
| | - Chukwudubem E. Ifeatu
- Computational and Bio-Simulation Research Group , University of Calabar, Calabar , Nigeria
| | - Tabe N. Ntui
- Computational and Bio-Simulation Research Group , University of Calabar, Calabar , Nigeria
- Department of Chemistry, Faculty of Physical Sciences , Cross River University of Technology , Calabar , Nigeria
| |
Collapse
|
9
|
Martínez JP, Trzaskowski B. Olefin Metathesis Catalyzed by a Hoveyda-Grubbs-like Complex Chelated to Bis(2-mercaptoimidazolyl) Methane: A Predictive DFT Study. J Phys Chem A 2022; 126:720-732. [PMID: 35080885 PMCID: PMC8842278 DOI: 10.1021/acs.jpca.1c09336] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
Although highly selective
complexes for the cross-metathesis of
olefins, particularly oriented toward the productive metathesis of Z-olefins, have been reported in recent years, there is
a constant need to design and prepare new and improved catalysts for
this challenging reaction. In this work, guided by density functional
theory (DFT) calculations, the performance of a Ru-based catalyst
chelated to a sulfurated pincer in the olefin metathesis was computationally
assessed. The catalyst was designed based on the Hoveyda–Grubbs
catalyst (SIMes)Cl2Ru(=CH–o–OiPrC6H4) through the substitution
of chlorides with the chelator bis(2-mercaptoimidazolyl)methane. The
obtained thermodynamic and kinetic data of the initiation phase through
side- and bottom-bound mechanisms suggest that this system is a versatile
catalyst for olefin metathesis, as DFT predicts the highest energy
barrier of the catalytic cycle of ca. 20 kcal/mol, which is comparable
to those corresponding to the Hoveyda–Grubbs-type catalysts.
Moreover, in terms of the stereoselectivity evaluated through the
propagation phase in the metathesis of propene–propene to 2-butene,
our study reveals that the Z isomer can be formed
under a kinetic control. We believe that this is an interesting outcome
in the context of future exploration of Ru-based catalysts with sulfurated
chelates in the search for high stereoselectivity in selected reactions.
Collapse
Affiliation(s)
- J Pablo Martínez
- Centre of New Technologies, University of Warsaw, 02-097 Warszawa, Poland
| | | |
Collapse
|
10
|
Bradley D, Branley C, Peeks MD. A straightforward method to quantify the electron-delocalizing ability of pi-conjugated molecules. Phys Chem Chem Phys 2022; 24:11486-11490. [DOI: 10.1039/d2cp01497a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Electronic delocalization is essential to the properties of pi-conjugated molecules. We introduce the inter-fragment delocalization index (IFDI) as an easy-to-use computational method for quantifying the electronic delocalization in pi-conjugated oligomers...
Collapse
|
11
|
Menant S, Guégan F, Tognetti V, Merzoud L, Joubert L, Chermette H, Morell C. Polarisation of Electron Density and Electronic Effects: Revisiting the Carbon-Halogen Bonds. Molecules 2021; 26:6218. [PMID: 34684795 PMCID: PMC8538414 DOI: 10.3390/molecules26206218] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2021] [Revised: 10/02/2021] [Accepted: 10/03/2021] [Indexed: 11/25/2022] Open
Abstract
Electronic effects (inductive and mesomeric) are of fundamental importance to understand the reactivity and selectivity of a molecule. In this article, polarisation temperature is used as a principal index to describe how electronic effects propagate in halogeno-alkanes and halogeno-alkenes. It is found that as chain length increases, polarisation temperature decreases. As expected, polarisation is much larger for alkenes than for alkanes. Finally, the polarisation mode of the carbon-fluorine bond is found to be quite different and might explain the unusual reactivity of fluoride compounds.
Collapse
Affiliation(s)
- Sébastien Menant
- Institut des Sciences Analytiques, UMR 5280, CNRS, Université de Lyon, 5 rue de la Doua, F-69100 Villeurbanne, France; (S.M.); (L.M.); (H.C.)
| | - Frédéric Guégan
- IC2MP UMR 7285, Université de Poitiers-CNRS, 4 rue Michel Brunet TSA, CEDEX 9, 86073 Poitiers, France
| | - Vincent Tognetti
- COBRA UMR 6014-FR 3038, Université de Rouen, INSA Rouen, CNRS, 1 rue Tesnière, CEDEX, 76821 Mont St Aignan, France;
| | - Lynda Merzoud
- Institut des Sciences Analytiques, UMR 5280, CNRS, Université de Lyon, 5 rue de la Doua, F-69100 Villeurbanne, France; (S.M.); (L.M.); (H.C.)
| | - Laurent Joubert
- COBRA UMR 6014-FR 3038, Université de Rouen, INSA Rouen, CNRS, 1 rue Tesnière, CEDEX, 76821 Mont St Aignan, France;
| | - Henry Chermette
- Institut des Sciences Analytiques, UMR 5280, CNRS, Université de Lyon, 5 rue de la Doua, F-69100 Villeurbanne, France; (S.M.); (L.M.); (H.C.)
| | - Christophe Morell
- Institut des Sciences Analytiques, UMR 5280, CNRS, Université de Lyon, 5 rue de la Doua, F-69100 Villeurbanne, France; (S.M.); (L.M.); (H.C.)
| |
Collapse
|
12
|
Tsuyuki M, Kugaya Y, Kanamori H, Yabushita S. Size-Dependent π g + π u Combination Band Intensities of Polyynes C 2nH 2 ( n = 1-9) Analyzed by the Local CCH Bending and the Linear Response Functions. J Phys Chem A 2021; 125:6404-6419. [PMID: 34275276 DOI: 10.1021/acs.jpca.1c04168] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Polyynes (C2nH2) show the unusually strong πg + πu combination bands in the infrared absorption spectra. We calculated them as the first overtone of the local CCH bending; the strong intensities are interpreted as a consequence of the large-amplitude bending vibration of the acidic acetylenic hydrogen combined with the size-dependent π electron conjugation. Our theoretical calculations show that the absorption intensity increases steadily and their increase rate is gradually slowed down by increasing the number of acetylene units up to n = 9. However, the calculated vibrational wavenumber converges quickly in agreement with the experimental observation. The second-order electron density deformation caused by the local CCH bending was analyzed using the linear response functions, including the linear and nonlinear contributions, to explain the n dependence. The easily polarizable π electron density caused two kinds of deformation-dominant but dark δxx-yy type and minor but bright σ type. Both of them exhibit interesting zigzag sign alternations, consistent with the law of alternating polarity of Coulson and Longuet-Higgins. The electron density polarization in these intra- and interacetylene units induces a large axial component molecular dipole moment, contributing to the intensity that increases with n. The difference between the curvilinear and rectilinear bending coordinates is interpreted within the present theoretical scheme.
Collapse
Affiliation(s)
- Masafumi Tsuyuki
- Department of Chemistry, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan
| | - Yuto Kugaya
- Department of Chemistry, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan
| | - Hideto Kanamori
- Department of Physics, Tokyo Institute of Technology, Ohokayama 2-12-1, Tokyo 152-8551, Japan
| | - Satoshi Yabushita
- Department of Chemistry, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan
| |
Collapse
|
13
|
Witek HA. Zhang-Zhang Polynomials of Multiple Zigzag Chains Revisited: A Connection with the John-Sachs Theorem. Molecules 2021; 26:molecules26092524. [PMID: 33925975 PMCID: PMC8123625 DOI: 10.3390/molecules26092524] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2021] [Revised: 04/09/2021] [Accepted: 04/22/2021] [Indexed: 12/05/2022] Open
Abstract
Multiple zigzag chains Zm,n of length n and width m constitute an important class of regular graphene flakes of rectangular shape. The physical and chemical properties of these basic pericondensed benzenoids can be related to their various topological invariants, conveniently encoded as the coefficients of a combinatorial polynomial, usually referred to as the ZZ polynomial of multiple zigzag chains Zm,n. The current study reports a novel method for determination of these ZZ polynomials based on a hypothesized extension to John–Sachs theorem, used previously to enumerate Kekulé structures of various benzenoid hydrocarbons. We show that the ZZ polynomial of the Zm,n multiple zigzag chain can be conveniently expressed as a determinant of a Toeplitz (or almost Toeplitz) matrix of size m2×m2 consisting of simple hypergeometric polynomials. The presented analysis can be extended to generalized multiple zigzag chains Zkm,n, i.e., derivatives of Zm,n with a single attached polyacene chain of length k. All presented formulas are accompanied by formal proofs. The developed theoretical machinery is applied for predicting aromaticity distribution patterns in large and infinite multiple zigzag chains Zm,n and for computing the distribution of spin densities in biradical states of finite multiple zigzag chains Zm,n.
Collapse
Affiliation(s)
- Henryk A. Witek
- Department of Applied Chemistry, Institute of Molecular Science, National Yang Ming Chiao Tung University, 1001 Ta-Hsueh Road, Hsinchu 300092, Taiwan;
- Center for Emergent Functional Matter Science, National Yang Ming Chiao Tung University, 1001 Ta-Hsueh Road, Hsinchu 300092, Taiwan
| |
Collapse
|
14
|
Tsuyuki M, Furudate S, Kugaya Y, Yabushita S. Graphical Transition Moment Decomposition and Conceptual Density Functional Theory Approaches to Study the Fundamental and Lower-Level Overtone Absorption Intensities of Some OH Stretching Vibrations. J Phys Chem A 2021; 125:2101-2113. [PMID: 33663218 DOI: 10.1021/acs.jpca.0c11619] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The investigation of electron density migrations caused by molecular structure changes is of central importance in various fields of chemistry. To address this topic in general and to study absorption intensities of vibrations, we analyze sensitive dipole moment functions (DMFs) of a molecule by combining the linear response function of conceptual DFT and bond dipoles separated by the quantum theory of atoms in molecules with a graphical transition moment decomposition scheme. The fundamental intensities of OH stretching vibrations depend strongly on the substituents but only weakly on the molecular conformations. Interestingly, in some alcohols, completely opposite trends have been observed for the lower-level overtone intensities: a weak substituent dependence but a stronger conformation dependence. It is well known that the formation of a hydrogen-bonded complex increases the OH stretching fundamental intensity, but less well known is the decrease in their overtone intensities. To investigate these characteristics comprehensively, we calculated their intensities (Δv = 1, 2, and 3) for conformers of ethanol and trifluoroethanol (TFE) and hydrogen-bonded phenol (PhOH) systems via the DFT method in the local mode model for the OH stretching coordinate ΔR. Their first and second derivatives of the electron density with respect to ΔR were calculated and interpreted using their bond moments. For ethanol and TFE, the OH, CC, and CH bond moments were found to make an important contribution to the molecular DMF derivatives parallel to the OH bond. The OH bond contributes only to the first derivative of DMF, and its conformational dependence is determined by the magnitude of the charge polarization of each structure. The electron density derivatives in the CC bond region were largely maintained during the internal rotation; thus, their conformation-dependent contributions were expressed by a geometrical factor of the CC bond direction. The CH bond at the antiperiplanar position of the OH bond was found to make a remarkably large contribution to the second derivative of DMF in the gauche conformer. The importance of electron density migration on substituents was also identified in the hydrogen-bonded phenol, in which the π-electron density change on the aromatic ring was clearly shown. This migration creates the DMF derivatives both perpendicular and parallel to the OH bond and strongly affects the absorption intensities. In all the cases, some bond moments on the substituents contribute to the first and second DMF derivatives in a structure-dependent manner, thus explaining their stereoelectronic effects.
Collapse
Affiliation(s)
- Masafumi Tsuyuki
- Department of Chemistry, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoku-Ku, Yokohama 223-8522, Japan
| | - Shunki Furudate
- Department of Chemistry, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoku-Ku, Yokohama 223-8522, Japan
| | - Yuto Kugaya
- Department of Chemistry, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoku-Ku, Yokohama 223-8522, Japan
| | - Satoshi Yabushita
- Department of Chemistry, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoku-Ku, Yokohama 223-8522, Japan
| |
Collapse
|
15
|
Vu KB, Le Phuc Nhi T, Vu VV, Tung Ngo S. How do magnetic, structural, and electronic criteria of aromaticity relate to HOMO – LUMO gap? An evaluation for graphene quantum dot and its derivatives. Chem Phys 2020. [DOI: 10.1016/j.chemphys.2020.110951] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
|
16
|
Stuyver T, Danovich D, Shaik S. Captodative Substitution Enhances the Diradical Character of Compounds, Reduces Aromaticity, and Controls Single-Molecule Conductivity Patterns: A Valence Bond Study. J Phys Chem A 2019; 123:7133-7141. [PMID: 31318209 DOI: 10.1021/acs.jpca.9b06096] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The present contribution uses a valence bond (VB) perspective to consider the captodative substitution strategy, a method to enhance the diradical character of (potentially aromatic) compounds. We confirm the qualitative reasoning that has generally been used to rationalize the diradical-character-enhancing effect of captodative substitution: this type of substitution scheme disproportionally stabilizes specific Dewar/diradical(oid) VB structures, thus increasing their weight in the full ground-state wave function. Furthermore, we assess the effect of captodative substitution on the aromaticity of the considered compound. We observe a clear trade-off between diradical character and aromaticity for our model systems: as one of these properties increases, the other decreases. This finding is especially significant within the field of single-molecule electronics because it enables unification of the previously observed inverse proportionality between the aromaticity of a compound and the magnitude of conductance through that molecule, with the observed proportionality between diradical character and the magnitude of conductance associated with a compound. To some extent, both properties, i.e., aromaticity and diradical character, appear to be the flip-sides of the same coin.
Collapse
Affiliation(s)
- Thijs Stuyver
- Department of Organic Chemistry and the Lise Meitner-Minerva Centre for Computational Quantum Chemistry , The Hebrew University , Jerusalem 91904 , Israel.,Algemene Chemie , Vrije Universiteit Brussel , Pleinlaan 2 , 1050 Brussels , Belgium
| | - David Danovich
- Department of Organic Chemistry and the Lise Meitner-Minerva Centre for Computational Quantum Chemistry , The Hebrew University , Jerusalem 91904 , Israel
| | - Sason Shaik
- Department of Organic Chemistry and the Lise Meitner-Minerva Centre for Computational Quantum Chemistry , The Hebrew University , Jerusalem 91904 , Israel
| |
Collapse
|
17
|
Yu D, Stuyver T, Rong C, Alonso M, Lu T, De Proft F, Geerlings P, Liu S. Global and local aromaticity of acenes from the information-theoretic approach in density functional reactivity theory. Phys Chem Chem Phys 2019; 21:18195-18210. [PMID: 31389933 DOI: 10.1039/c9cp01623f] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
In this work, we report a systematic study on the global and local aromaticity of acenes using a series of model structures from 2-acene to 11-acene. A recently developed ansatz, an information-theoretic approach coached into density functional reactivity theory has been employed, which essentially provides different density functionals characterizing the molecular electron density distribution. Based on the correlation analysis of six conventional aromaticity indices with eight information-theoretic quantities, we examined the aromaticity of acenes from both global and local perspectives. From the global aromaticity viewpoint, our results suggest that different descriptors based on various physicochemical properties are intrinsically dependent. A novel laminated feature ruling local aromaticity of acenes has been unveiled, from which we found that the distance from the terminal rings plays the critical role. Based on the shape of the correlation plots between the conventional aromaticity indices and information-theoretic quantities, the latter could be separated into three subgroups. The seemingly contradictory results from global and local aromaticity perspectives not only present us the uniqueness of the acene systems but all demonstrate the effectiveness of the information-theoretic approach from density functional reactivity theory. Besides strengthening the validity of a series of new aromaticity descriptors, our results should lead to more clear insights into the chemical significance of the information-theoretic quantities.
Collapse
Affiliation(s)
- Donghai Yu
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), College of Chemistry and Chemical Engineering, Hunan Normal University, Changsha, 410081, China.
| | | | | | | | | | | | | | | |
Collapse
|
18
|
Firouzi R, Shafie H, Tohidnia H. Characterization of Local Aromaticity in Polycyclic Conjugated Hydrocarbons Based on Anisotropy of π-Electron Density. ChemistrySelect 2017. [DOI: 10.1002/slct.201702407] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Rohoullah Firouzi
- Department of Physical Chemistry; Chemistry and Chemical Engineering Research Center of Iran; P.O. Box 14968-13151 Tehran Iran
| | - Hoda Shafie
- Department of Physical Chemistry; Chemistry and Chemical Engineering Research Center of Iran; P.O. Box 14968-13151 Tehran Iran
| | - Hassan Tohidnia
- Department of Physical Chemistry; Chemistry and Chemical Engineering Research Center of Iran; P.O. Box 14968-13151 Tehran Iran
| |
Collapse
|
19
|
Tognetti V, Guégan F, Luneau D, Chermette H, Morell C, Joubert L. Structural effects in octahedral carbonyl complexes: an atoms-in-molecules study. Theor Chem Acc 2017. [DOI: 10.1007/s00214-017-2116-9] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
20
|
Setiawan D, Kraka E, Cremer D. Quantitative Assessment of Aromaticity and Antiaromaticity Utilizing Vibrational Spectroscopy. J Org Chem 2016; 81:9669-9686. [DOI: 10.1021/acs.joc.6b01761] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Dani Setiawan
- Computational and Theoretical
Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave., Dallas, Texas 75275-0314, United States
| | - Elfi Kraka
- Computational and Theoretical
Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave., Dallas, Texas 75275-0314, United States
| | - Dieter Cremer
- Computational and Theoretical
Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave., Dallas, Texas 75275-0314, United States
| |
Collapse
|
21
|
|
22
|
The nucleophilicity equalization principle and new algorithms for the evaluation of molecular nucleophilicity. COMPUT THEOR CHEM 2016. [DOI: 10.1016/j.comptc.2016.02.006] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
23
|
Feixas F, Matito E, Poater J, Solà M. Quantifying aromaticity with electron delocalisation measures. Chem Soc Rev 2015; 44:6434-51. [PMID: 25858673 DOI: 10.1039/c5cs00066a] [Citation(s) in RCA: 291] [Impact Index Per Article: 32.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Aromaticity cannot be measured directly by any physical or chemical experiment because it is not a well-defined magnitude. Its quantification is done indirectly from the measure of different properties that are usually found in aromatic compounds such as bond length equalisation, energetic stabilisation, and particular magnetic behaviour associated with induced ring currents. These properties have been used to set up the myriad of structural-, energetic-, and magnetic-based indices of aromaticity known to date. The cyclic delocalisation of mobile electrons in two or three dimensions is probably one of the key aspects that characterise aromatic compounds. However, it has not been until the last decade that electron delocalisation measures have been widely employed to quantify aromaticity. Some of these new indicators of aromaticity such as the PDI, FLU, ING, and INB were defined in our group. In this paper, we review the different existing descriptors of aromaticity that are based on electron delocalisation properties, we compare their performance with indices based on other properties, and we summarise a number of applications of electronic-based indices for the analysis of aromaticity in interesting chemical problems.
Collapse
Affiliation(s)
- Ferran Feixas
- Institut de Química Computacional i Catàlisi (IQCC) and Departament de Química, Universitat de Girona, Campus Montilivi, 17071 Girona, Spain.
| | | | | | | |
Collapse
|
24
|
Boisdenghien Z, Fias S, Da Pieve F, De Proft F, Geerlings P. The polarisability of atoms and molecules: a comparison between a conceptual density functional theory approach and time-dependent density functional theory. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1021110] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Zino Boisdenghien
- General Chemistry (ALGC), Vrije Universiteit Brussel (Free University Brussels – VUB) Brussel, Belgium
| | - Stijn Fias
- General Chemistry (ALGC), Vrije Universiteit Brussel (Free University Brussels – VUB) Brussel, Belgium
| | - Fabiana Da Pieve
- European Theoretical Spectroscopy Facility (ETSF), Ecole Polytechnique, Palaiseau, France
- Laboratoire des Solides Irradiés, Ecole Polytechnique, Palaiseau, France
| | - Frank De Proft
- General Chemistry (ALGC), Vrije Universiteit Brussel (Free University Brussels – VUB) Brussel, Belgium
| | - Paul Geerlings
- General Chemistry (ALGC), Vrije Universiteit Brussel (Free University Brussels – VUB) Brussel, Belgium
| |
Collapse
|
25
|
|
26
|
Stuyver T, Fias S, De Proft F, Fowler PW, Geerlings P. Conduction of molecular electronic devices: Qualitative insights through atom-atom polarizabilities. J Chem Phys 2015; 142:094103. [DOI: 10.1063/1.4913415] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
27
|
Fias S, Boisdenghien Z, De Proft F, Geerlings P. The spin polarized linear response from density functional theory: Theory and application to atoms. J Chem Phys 2014; 141:184107. [DOI: 10.1063/1.4900513] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
28
|
Geerlings P, Fias S, Boisdenghien Z, De Proft F. Conceptual DFT: chemistry from the linear response function. Chem Soc Rev 2014; 43:4989-5008. [PMID: 24531142 DOI: 10.1039/c3cs60456j] [Citation(s) in RCA: 112] [Impact Index Per Article: 11.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
Within the context of reactivity descriptors known in conceptual DFT, the linear response function (χ(r,r')) remained nearly unexploited. Although well known, in its time dependent form, in the solid state physics and time-dependent DFT communities the study of the "chemistry" present in the kernel was, until recently, relatively unexplored. The evaluation of the linear response function as such and its study in the time independent form are highlighted in the present review. On the fundamental side, the focus is on the approaches of increasing complexity to compute and represent χ(r,r'), its visualisation going from plots of the unintegrated χ(r,r') to an atom condensed matrix. The study on atoms reveals its physical significance, retrieving atomic shell structure, while the results on molecules illustrate that a variety of chemical concepts are retrieved: inductive and mesomeric effects, electron delocalisation, aromaticity and anti-aromaticity, σ and π aromaticity,…. The applications show that the chemistry of aliphatic (saturated and unsaturated) chains, saturated and aromatic/anti-aromatic rings, organic, inorganic or metallic in nature, can be retrieved via the linear response function, including the variation of the electronic structure of the reagents along a reaction path. The connection of the linear response function with the concept of nearsightedness and the alchemical derivatives is also highlighted.
Collapse
Affiliation(s)
- Paul Geerlings
- General Chemistry (ALGC), Vrije Universiteit Brussel (Free University Brussels-VUB), Pleinlaan 2, 1050 Brussel, Belgium.
| | | | | | | |
Collapse
|
29
|
Boisdenghien Z, Fias S, Van Alsenoy C, De Proft F, Geerlings P. Evaluating and interpreting the chemical relevance of the linear response kernel for atoms II: open shell. Phys Chem Chem Phys 2014; 16:14614-24. [DOI: 10.1039/c4cp01331j] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
30
|
Firouzi R, Sharifi Ardani S. Description of heteroaromaticity on the basis of π-electron density anisotropy. Phys Chem Chem Phys 2014; 16:11538-48. [DOI: 10.1039/c4cp01125b] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
31
|
Balawender R, Welearegay MA, Lesiuk M, De Proft F, Geerlings P. Exploring Chemical Space with the Alchemical Derivatives. J Chem Theory Comput 2013; 9:5327-40. [PMID: 26592270 DOI: 10.1021/ct400706g] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
In this paper, we verify the usefulness of the alchemical derivatives in the prediction of chemical properties. We concentrate on the stability of the transmutation products, where the term "transmutation" means the change of the nuclear charge at an atomic site at constant number of electrons. As illustrative transmutations showing the potential of the method in exploring chemical space, we present some examples of increasing complexity starting with the deprotonation, continuing with the transmutation of the nitrogen molecule, and ending with the substitution of isoelectronic B-N units for C-C units and N units for C-H units in carbocyclic systems. The basis set influence on the qualitative and quantitative accuracies of the alchemical predictions was investigated. The alchemical deprotonation energy (from the second order Taylor expansion) correlates well with the vertical deprotonation energy and can be used as a preliminary indicator for the experimental deprotonation energy. The results of calculations for the BN derivatives of benzene and pyrene show that this method has great potential for efficient and accurate scanning of chemical space.
Collapse
Affiliation(s)
- Robert Balawender
- Institute of Physical Chemistry, Polish Academy of Sciences , Kasprzaka 44/52, PL-01-224 Warsaw, Poland
| | - Meressa A Welearegay
- Institute of Physical Chemistry, Polish Academy of Sciences , Kasprzaka 44/52, PL-01-224 Warsaw, Poland
| | - Michał Lesiuk
- Faculty of Chemistry, University of Warsaw , Pasteura 1, PL-02-093 Warsaw, Poland
| | - Frank De Proft
- Eenheid Algemene Chemie, Vrije Universiteit Brussel, Faculteit Wetenschappen , Pleinlaan 2, B-1050 Brussels, Belgium
| | - Paul Geerlings
- Eenheid Algemene Chemie, Vrije Universiteit Brussel, Faculteit Wetenschappen , Pleinlaan 2, B-1050 Brussels, Belgium
| |
Collapse
|
32
|
|
33
|
Fias S, Boisdenghien Z, Stuyver T, Audiffred M, Merino G, Geerlings P, de Proft F. Analysis of Aromaticity in Planar Metal Systems using the Linear Response Kernel. J Phys Chem A 2013; 117:3556-60. [DOI: 10.1021/jp401760j] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Stijn Fias
- General Chemistry, Free University of Brussels (VUB), Pleinlaan 2, B-1050 Brussels, Belgium
| | - Zino Boisdenghien
- General Chemistry, Free University of Brussels (VUB), Pleinlaan 2, B-1050 Brussels, Belgium
| | - Thijs Stuyver
- General Chemistry, Free University of Brussels (VUB), Pleinlaan 2, B-1050 Brussels, Belgium
| | - Martha Audiffred
- Departamento de Física Aplicada, Centro de Investigación y de Estudios Avanzados, Unidad Mérida, Km. 6 Antigua Carretera a Progreso, Apdo.
Postal 73, Cordemex, 97310, Mérida, Yuc., México
| | - Gabriel Merino
- Departamento de Física Aplicada, Centro de Investigación y de Estudios Avanzados, Unidad Mérida, Km. 6 Antigua Carretera a Progreso, Apdo.
Postal 73, Cordemex, 97310, Mérida, Yuc., México
| | - Paul Geerlings
- General Chemistry, Free University of Brussels (VUB), Pleinlaan 2, B-1050 Brussels, Belgium
| | - Frank de Proft
- General Chemistry, Free University of Brussels (VUB), Pleinlaan 2, B-1050 Brussels, Belgium
| |
Collapse
|
34
|
Stachowicz A, Rogalski M, Korchowiec J. Charge sensitivity approach to mutual polarization of reactants: molecular mechanics perspective. J Mol Model 2013; 19:4163-72. [PMID: 23392762 PMCID: PMC4271185 DOI: 10.1007/s00894-013-1757-4] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2012] [Accepted: 01/03/2013] [Indexed: 11/26/2022]
Abstract
Charge sensitivity analysis (CSA) in force-field atoms resolution was applied to describe the mutual polarization of reactants as well as charge-transfer (CT) effects. An inclusion complex of β-cyclodextrin with salicylic acid was used as a model system. Three CSA models were taken into account and verified on a Born–Oppenheimer molecular dynamics (BOMD) trajectory. The models differed in terms of the equilibrium conditions imposed on the system. It was demonstrated that mutual polarization is an important source of stabilization, in contrast to the results obtained from static charge calculations. The energy lowering induced by CT was small and comparable to the CT stabilization that occurs in hydrogen-bonded systems. All models correctly described the main topological features of the BOMD energy surface. CSA in force-field atoms resolution qualitatively reproduced the charge reorganization accompanying hydrogen-bond formation. It was shown that CSA parameters are very sensitive to the bond formation process, which suggests that they could be applied in reactive force fields as detectors of newly formed chemical bonds. Fukui function detector of bond formation during molecular dynamics simulations of inclusion complex of b-cyclodextrin with salicylic acid. ![]()
Collapse
Affiliation(s)
- Anna Stachowicz
- K. Gumiński Department of Theoretical Chemistry, Faculty of Chemistry, Jagiellonian University, R. Ingardena 3, Kraków, Poland
| | - Marek Rogalski
- Laboratoire de Chimie et Physique - Approche Multi-Echelle des Milieux Complexes, EA 4164, Université de Lorraine, 1 Boulevard Arago, 57070 Metz, France
| | - Jacek Korchowiec
- K. Gumiński Department of Theoretical Chemistry, Faculty of Chemistry, Jagiellonian University, R. Ingardena 3, Kraków, Poland
| |
Collapse
|
35
|
Fias S, Geerlings P, Ayers P, De Proft F. σ, π aromaticity and anti-aromaticity as retrieved by the linear response kernel. Phys Chem Chem Phys 2013; 15:2882-9. [PMID: 23337925 DOI: 10.1039/c2cp43612d] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The chemical importance of the linear response kernel from conceptual Density Functional Theory (DFT) is investigated for some σ and π aromatic and anti-aromatic systems. The effect of the ring size is studied by looking at some well known aromatic and anti-aromatic molecules of different sizes, showing that the linear response is capable of correctly classifying and quantifying the aromaticity for five- to eight-membered aromatic and anti-aromatic molecules. The splitting of the linear response in σ and π contributions is introduced and its significance is illustrated using some σ-aromatic molecules. The linear response also correctly predicts the aromatic transition states of the Diels-Alder reaction and the acetylene trimerisation and shows the expected behavior along the reaction coordinate, proving that the method is accurate not only at the minimum of the potential energy surface, but also in non-equilibrium states. Finally, the reason for the close correlation between the linear response and the Para Delocalisation Index (PDI), found in previous and the present study, is proven mathematically. These results show the linear response to be a reliable DFT-index to probe the σ and π aromaticity or anti-aromaticity of a broad range of molecules.
Collapse
Affiliation(s)
- Stijn Fias
- General Chemistry, Free University of Brussels, Pleinlaan 2, B-1050 Brussels, Belgium.
| | | | | | | |
Collapse
|
36
|
Boisdenghien Z, Van Alsenoy C, De Proft F, Geerlings P. Evaluating and Interpreting the Chemical Relevance of the Linear Response Kernel for Atoms. J Chem Theory Comput 2013; 9:1007-15. [DOI: 10.1021/ct300861r] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Zino Boisdenghien
- Eenheid Algemene Chemie, Vrije
Universiteit Brussel, Pleinlaan 2, 1050 Brussels, Belgium
| | - Christian Van Alsenoy
- Departement Chemie, Universiteit
Antwerpen, Campus Drie Eiken, Universiteitsplein 1, 2610 Wilrijk,
Belgium
| | - Frank De Proft
- Eenheid Algemene Chemie, Vrije
Universiteit Brussel, Pleinlaan 2, 1050 Brussels, Belgium
| | - Paul Geerlings
- Eenheid Algemene Chemie, Vrije
Universiteit Brussel, Pleinlaan 2, 1050 Brussels, Belgium
| |
Collapse
|
37
|
Tognetti V, Morell C, Ayers PW, Joubert L, Chermette H. A proposal for an extended dual descriptor: a possible solution when Frontier Molecular Orbital Theory fails. Phys Chem Chem Phys 2013; 15:14465-75. [DOI: 10.1039/c3cp51169c] [Citation(s) in RCA: 59] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
38
|
Yang W, Cohen AJ, De Proft F, Geerlings P. Analytical evaluation of Fukui functions and real-space linear response function. J Chem Phys 2012; 136:144110. [DOI: 10.1063/1.3701562] [Citation(s) in RCA: 56] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|