1
|
Dakkouri M. A Theoretical Investigation of Novel Sila- and Germa-Spirocyclic Imines and Their Relevance for Electron-Transporting Materials and Drug Discovery. Molecules 2023; 28:6298. [PMID: 37687127 PMCID: PMC10489060 DOI: 10.3390/molecules28176298] [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: 07/23/2023] [Revised: 08/21/2023] [Accepted: 08/22/2023] [Indexed: 09/10/2023] Open
Abstract
A new class of spirocyclic imines (SCIs) has been theoretically investigated by applying a variety of quantum chemical methods and basis sets. The uniqueness of these compounds is depicted by various peculiarities, e.g., the incidence of planar six-membered rings each with two imine groups (two π bonds) and the incorporation of the isosteres carbon, silicon, or germanium spiro centers. Additional peculiarities of these novel SCIs are mirrored by their three-dimensionality, the simultaneous occurrence of nucleophilic and electrophilic centers, and the cross-hyperconjugative (spiro-conjugation) interactions, which provoke charge mobility along the spirocyclic scaffold. Substitution of SCIs with strong electron-withdrawing substituents, like the cyano group or fluorine, enhances their docking capability and impacts their reactivity and charge mobility. To gain thorough knowledge about the molecular properties of these SCIs, their structures have been optimized and various quantum chemical concepts and models were applied, e.g., full NBO analysis and the frontier molecular orbitals (FMOs) theory (HOMO-LUMO energy gap) and the chemical reactivity descriptors derived from them. For the assessment of the charge density distribution along the SCI framework, additional complementary quantum chemical methods were used, e.g., molecular electrostatic potential (MESP) and Bader's QTAIM. Additionally, using the aromaticity index NICS (nuclear independent chemical shift) and other criteria, it could be shown that the investigated cross-hyperconjugated sila and germa SCIs are spiro-aromatics of the Heilbronner Craig-type Möbius aromaticity.
Collapse
Affiliation(s)
- Marwan Dakkouri
- Department of Electrochemistry, University of Ulm, D-89069 Ulm, Germany
| |
Collapse
|
2
|
Shustov GV. Theoretical study of water-assisted reaction between methane and nitrosonium cation. CAN J CHEM 2022. [DOI: 10.1139/cjc-2021-0299] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
The water-assisted reaction between CH4 and NO+ has been studied employing post Hartree–Fock and density functional theory methods. Two reaction pathways were considered: (i) a frontside NO+ attack on the C–H bond of methane with retention of configuration and (ii) a backside attack with inversion of configuration. The first pathway leads to a thermodynamically more favorable hydrate of N-protonated nitrosomethane, and the second one leads to its O-protonated isomer. The catalytic effect of a single water molecule is expressed in decreasing the activation energy for the elementary step by about a factor of two for the frontside mechanism and by a factor of four to five for the backside one when compared with the calculated literature data on water-free methane nitrosation. The combination of the activation strain model of reactivity and the energy decomposition analysis reveals that the activation barriers are largely determined by the relative stability of the termolecular reactant complexes. The crucial factor that stabilizes these complexes comes from the electrostatic attraction. The catalytic effect of the key water molecule is decreased with introduction of additional one or two explicit water molecules, which form a coordinate O→N bonding with NO+. On the contrary, an additional water molecule hydrogen bonded to the key catalytic water molecule enhances the catalytic effect.
Collapse
Affiliation(s)
- Gennady V. Shustov
- Department of Chemistry, University of the Fraser Valley, Abbotsford, BC V2S 7M8, Canada
- Department of Chemistry, University of the Fraser Valley, Abbotsford, BC V2S 7M8, Canada
| |
Collapse
|
3
|
Han B, Isborn CM, Shi L. Determining Partial Atomic Charges for Liquid Water: Assessing Electronic Structure and Charge Models. J Chem Theory Comput 2021; 17:889-901. [PMID: 33405925 DOI: 10.1021/acs.jctc.0c01102] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Partial atomic charges provide an intuitive and efficient way to describe the charge distribution and the resulting intermolecular electrostatic interactions in liquid water. Many charge models exist and it is unclear which model provides the best assignment of partial atomic charges in response to the local molecular environment. In this work, we systematically scrutinize various electronic structure methods and charge models (Mulliken, natural population analysis, CHelpG, RESP, Hirshfeld, Iterative Hirshfeld, and Bader) by evaluating their performance in predicting the dipole moments of isolated water, water clusters, and liquid water as well as charge transfer in the water dimer and liquid water. Although none of the seven charge models is capable of fully capturing the dipole moment increase from isolated water (1.85 D) to liquid water (about 2.9 D), the Iterative Hirshfeld method performs best for liquid water, reproducing its experimental average molecular dipole moment, yielding a reasonable amount of intermolecular charge transfer, and showing modest sensitivity to the local water environment. The performance of the charge model is dependent on the choice of the density functional and the quantum treatment of the environment. The computed molecular dipole moment of water generally increases with the percentage of the exact Hartree-Fock exchange in the functional, whereas the amount of charge transfer between molecules decreases. For liquid water, including two full solvation shells of surrounding water molecules (within about 5.5 Å of the central water) in the quantum chemical calculation converges the charges of the central water molecule. Our final pragmatic quantum chemical charge-assigning protocol for liquid water is the Iterative Hirshfeld method with M06-HF/aug-cc-pVDZ and a quantum region cutoff radius of 5.5 Å.
Collapse
Affiliation(s)
- Bowen Han
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Liang Shi
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
4
|
Mason TG, Seeger ZL, Nguyen ALP, Fujita K, Izgorodina EI. Predicting Entropic Effects of Water Mixing with Ionic Liquids Containing Anions of Strong Hydrogen Bonding Ability: Role of the Cation. J Phys Chem B 2020; 124:9182-9194. [PMID: 33007160 DOI: 10.1021/acs.jpcb.0c07732] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Ionic liquids (ILs) such as choline dihydrogen phosphate exhibit an extraordinary solubilizing ability for proteins such as cytochrome C when mixed with 20 wt % water. Most widely used imidazolium-based ionic liquids coupled with dihydrogen phosphate do not exhibit the same solubilizing properties, suggesting that a multifunctional cation such as choline might play a key role in enhancing these properties of ionic liquid mixtures with water. In this theoretical work, we compare intermolecular interactions between the water molecule and ionic liquid ions in two ion-paired clusters of choline- and 1-butyl-3-methyl-imidazolium-based ionic liquids coupled with acetate, dihydrogen phosphate, and mesylate. Gibbs free energy (GFE) of solvation of water in these ionic liquids was calculated. Incorporation of a water molecule into ionic liquid clusters was accompanied by negative GFEs of solvation in both types of cations. These results were in good agreement with previously reported experimental GFEs of solvation of water in ILs. Compared to imidazolium-based clusters, strong interionic interactions of choline ionic liquids resulted in more negative GFEs due to their smaller deformation upon the addition of a water molecule, with dihydrogen phosphate and mesylate predicting the lowest GFEs of -30.1 and -43.5 kJ/mol-1, respectively. Lower GFEs of solvation of water in choline-based clusters were also accompanied with smaller entropic penalties, suggesting that water easily incorporates itself into the existing ionic network. Analysis of the intramolecular bonds within the water molecule showed that the choline hydroxyl group donates electron density to the neighboring water molecule, leading to additional polarization. The predicted infrared spectra of clusters of ionic liquids with water showed a pronounced red shift due to strongly polarized O-H bonds, in excellent agreement with the experimentally measured infrared spectra of ionic liquid mixtures with water. Increased polarization of water in choline-based ionic liquids undoubtedly creates more effective solvents for stabilizing biological molecules such as proteins.
Collapse
Affiliation(s)
- Thomas G Mason
- School of Chemistry, Monash University, 17 Rainforest Walk, Clayton, VIC 3800, Australia
| | - Zoe L Seeger
- School of Chemistry, Monash University, 17 Rainforest Walk, Clayton, VIC 3800, Australia
| | - Anh L P Nguyen
- School of Chemistry, Monash University, 17 Rainforest Walk, Clayton, VIC 3800, Australia
| | - Kyoko Fujita
- Department of Pathophysiology, School of Pharmacy, Tokyo University of Pharmacy and Life Sciences, 1432-1 Horinouchi, Hachioji, Tokyo 192-0392, Japan
| | - Ekaterina I Izgorodina
- School of Chemistry, Monash University, 17 Rainforest Walk, Clayton, VIC 3800, Australia
| |
Collapse
|
5
|
Sidler D, Meuwly M, Hamm P. An efficient water force field calibrated against intermolecular THz and Raman spectra. J Chem Phys 2018; 148:244504. [DOI: 10.1063/1.5037062] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- David Sidler
- Department of Chemistry, University of Zurich, Zurich, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Basel, Switzerland
| | - Peter Hamm
- Department of Chemistry, University of Zurich, Zurich, Switzerland
| |
Collapse
|
6
|
Heterogeneous nucleation of polymorphs on polymer surfaces: polymer-molecule interactions using a heterogeneous dielectric solvation model. J Mol Model 2018; 24:156. [PMID: 29886563 DOI: 10.1007/s00894-018-3657-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2017] [Accepted: 04/13/2018] [Indexed: 10/14/2022]
Abstract
We have investigated the mechanism of the nucleation of acetaminophen on poly(methyl-methacrylate) and poly(vinyl-acetate) utilizing a combination of quantum mechanical computations and electrostatic models. We have used a heterogeneous dielectric solvation model to determine the stability of different orientations of acetaminophen on polymer surfaces. We find that for the nucleation of acetaminophen on the polymer surfaces in vacuum, the most stable orientation is a flat orientation. For the nucleation process in solution where acetaminophen and the polymer surface are surrounded by a solvent, we find that the heterogeneous dielectric solvation model predicts that a sideways orientation is the most stable orientation.
Collapse
|
7
|
Molecular Properties of Sandwiched Molecules Between Electrodes and Nanoparticles. ADVANCES IN QUANTUM CHEMISTRY 2017. [DOI: 10.1016/bs.aiq.2017.03.002] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
8
|
Shahid HA, Jahangir S, Yousuf S, Hanif M, Sherwani SK. Synthesis, crystal structure, structural characterization and in vitro antimicrobial activities of 1-methyl-4-nitro-1H-imidazole. ARAB J CHEM 2016. [DOI: 10.1016/j.arabjc.2014.11.001] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022] Open
|
9
|
Phipps MJS, Fox T, Tautermann CS, Skylaris CK. Energy Decomposition Analysis Based on Absolutely Localized Molecular Orbitals for Large-Scale Density Functional Theory Calculations in Drug Design. J Chem Theory Comput 2016; 12:3135-48. [PMID: 27248370 DOI: 10.1021/acs.jctc.6b00272] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We report the development and implementation of an energy decomposition analysis (EDA) scheme in the ONETEP linear-scaling electronic structure package. Our approach is hybrid as it combines the localized molecular orbital EDA (Su, P.; Li, H. J. Chem. Phys., 2009, 131, 014102) and the absolutely localized molecular orbital EDA (Khaliullin, R. Z.; et al. J. Phys. Chem. A, 2007, 111, 8753-8765) to partition the intermolecular interaction energy into chemically distinct components (electrostatic, exchange, correlation, Pauli repulsion, polarization, and charge transfer). Limitations shared in EDA approaches such as the issue of basis set dependence in polarization and charge transfer are discussed, and a remedy to this problem is proposed that exploits the strictly localized property of the ONETEP orbitals. Our method is validated on a range of complexes with interactions relevant to drug design. We demonstrate the capabilities for large-scale calculations with our approach on complexes of thrombin with an inhibitor comprised of up to 4975 atoms. Given the capability of ONETEP for large-scale calculations, such as on entire proteins, we expect that our EDA scheme can be applied in a large range of biomolecular problems, especially in the context of drug design.
Collapse
Affiliation(s)
- M J S Phipps
- School of Chemistry, University of Southampton , Highfield, Southampton SO17 1BJ, United Kingdom
| | - T Fox
- Lead Identification and Optimization Support, Boehringer Ingelheim Pharma GmbH & Co. KG , 88397 Biberach, Germany
| | - C S Tautermann
- Lead Identification and Optimization Support, Boehringer Ingelheim Pharma GmbH & Co. KG , 88397 Biberach, Germany
| | - C-K Skylaris
- School of Chemistry, University of Southampton , Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
10
|
Rick SW. A polarizable, charge transfer model of water using the drude oscillator. J Comput Chem 2016; 37:2060-6. [DOI: 10.1002/jcc.24426] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2015] [Revised: 04/14/2016] [Accepted: 05/17/2016] [Indexed: 12/24/2022]
Affiliation(s)
- Steven W. Rick
- Department of ChemistryUniversity of New OrleansNew Orleans70148 Los Angeles
| |
Collapse
|
11
|
Giese TJ, York DM. Ambient-Potential Composite Ewald Method for ab Initio Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulation. J Chem Theory Comput 2016; 12:2611-32. [PMID: 27171914 DOI: 10.1021/acs.jctc.6b00198] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
A new approach for performing Particle Mesh Ewald in ab initio quantum mechanical/molecular mechanical (QM/MM) simulations with extended atomic orbital basis sets is presented. The new approach, the Ambient-Potential Composite Ewald (CEw) method, does not perform the QM/MM interaction with Mulliken charges nor electrostatically fit charges. Instead the nuclei and electron density interact directly with the MM environment, but in a manner that avoids the use of dense Fourier transform grids. By performing the electrostatics with the underlying QM density, the CEw method avoids self-consistent field instabilities that have been encountered with simple charge mapping procedures. Potential of mean force (PMF) profiles of the p-nitrophenyl phosphate dissociation reaction in explicit solvent are computed from PBE0/6-31G* QM/MM molecular dynamics simulations with various electrostatic protocols. The CEw profiles are shown to be stable with respect to real-space Ewald cutoff, whereas the PMFs computed from truncated and switched electrostatics produce artifacts. PBE0/6-311G**, AM1/d-PhoT, and DFTB2 QM/MM simulations are performed to generate two-dimensional PMF profiles of the phosphoryl transesterification reactions with ethoxide and phenoxide leaving groups. The semiempirical models incorrectly produce a concerted ethoxide mechanism, whereas PBE0 correctly produces a stepwise mechanism. The ab initio reaction barriers agree more closely to experiment than the semiempirical models. The failure of Mulliken-charge QM/MM-Ewald is analyzed.
Collapse
Affiliation(s)
- Timothy J Giese
- Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University , Piscataway, New Jersey 08854-8087, United States
| | - Darrin M York
- Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University , Piscataway, New Jersey 08854-8087, United States
| |
Collapse
|
12
|
Hanif M, Chen L, Zhu L, Zhao D, Xiong T, Hou H. Stille cross-coupling applied to get higher molecular weight polymers: Synthesis, optoelectronic, Vocproperties, and solar cell application. J Appl Polym Sci 2015. [DOI: 10.1002/app.42147] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Muddasir Hanif
- Department of Chemistry and Chemical Engineering; Jiangxi Normal University; Nanchang Jiangxi 330022 People's Republic of China
| | - Linlin Chen
- Department of Chemistry and Chemical Engineering; Jiangxi Normal University; Nanchang Jiangxi 330022 People's Republic of China
| | - Li Zhu
- Department of Chemistry and Chemical Engineering; Jiangxi Normal University; Nanchang Jiangxi 330022 People's Republic of China
| | - Dan Zhao
- Department of Chemistry and Chemical Engineering; Jiangxi Normal University; Nanchang Jiangxi 330022 People's Republic of China
| | - Tianrou Xiong
- Department of Chemistry and Chemical Engineering; Jiangxi Normal University; Nanchang Jiangxi 330022 People's Republic of China
| | - Haoqing Hou
- Department of Chemistry and Chemical Engineering; Jiangxi Normal University; Nanchang Jiangxi 330022 People's Republic of China
| |
Collapse
|
13
|
Cummings OT, Wick CD. Interfacial behavior of simple inorganic salts at the air-water interface investigated with a polarizable model with electrostatic damping. J Chem Phys 2013; 139:064708. [DOI: 10.1063/1.4817775] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
14
|
Teodoro TQ, Haiduke RLA. Atomic charge and atomic dipole fluxes during stretching displacements in small molecules. COMPUT THEOR CHEM 2013. [DOI: 10.1016/j.comptc.2012.11.017] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
15
|
Viana RB, Santos EDA, Valencia LJ, Cavalcante RM, Costa EB, Moreno-Fuquen R, da Silva ABF. 4-hydroxy-2,5-dimethylphenyl-benzophenone: conformational stability, FT-IR and Raman investigation. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2013; 102:386-392. [PMID: 23220682 DOI: 10.1016/j.saa.2012.09.094] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/23/2012] [Revised: 06/30/2012] [Accepted: 09/14/2012] [Indexed: 06/01/2023]
Abstract
In this paper we have studied the 4-hydroxyl-2,5 dimethylphenyl-benzophenone. Also, it was analyzed the vibrational assignment (FT-IR and Raman) in conjunction with computational results. The conformational analysis showed three barrier heights where two are due to the dihedral rotation and the other one is attributed to hydroxyl rotation. While the high dihedral rotational barrier (TS1) is 6.06-7.22 kcal mol(-1), the lower one (TS2) is almost three times smaller. The variations with the change of basis set is 5-8% to TS1, and 3-15% in the values predicted to TS2. In the case of OH rotational barrier, the values range from 3.70 to 4.86 kcal mol(-1), and it is also observed that this transition state is less sensitive to the change of basis set and to the method. Two isomers was detected due to the changes in the OH rotation with the gap energy lower than 0.7 kcal mol(-1), and at this point is seen that semi-empirical methods fail into describe the most stable conformation which may be due to the small energy gap. The enthalpy formation at 0 K and 298 K was 111.71 and 102.20 kcal mol(-1), respectively.
Collapse
Affiliation(s)
- Rommel B Viana
- Instituto de Química de São Carlos, Universidade de São Paulo, São Carlos, SP, Brazil.
| | | | | | | | | | | | | |
Collapse
|
16
|
Wick CD, Lee AJ, Rick SW. How intermolecular charge transfer influences the air-water interface. J Chem Phys 2012; 137:154701. [PMID: 23083178 DOI: 10.1063/1.4758457] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Affiliation(s)
- Collin D Wick
- Louisiana Tech University, Ruston, Louisina 71270, USA
| | | | | |
Collapse
|
17
|
Isegawa M, Gao J, Truhlar DG. Incorporation of charge transfer into the explicit polarization fragment method by grand canonical density functional theory. J Chem Phys 2011; 135:084107. [PMID: 21895159 DOI: 10.1063/1.3624890] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Molecular fragmentation algorithms provide a powerful approach to extending electronic structure methods to very large systems. Here we present a method for including charge transfer between molecular fragments in the explicit polarization (X-Pol) fragment method for calculating potential energy surfaces. In the conventional X-Pol method, the total charge of each fragment is preserved, and charge transfer between fragments is not allowed. The description of charge transfer is made possible by treating each fragment as an open system with respect to the number of electrons. To achieve this, we applied Mermin's finite temperature method to the X-Pol wave function. In the application of this method to X-Pol, the fragments are open systems that partially equilibrate their number of electrons through a quasithermodynamics electron reservoir. The number of electrons in a given fragment can take a fractional value, and the electrons of each fragment obey the Fermi-Dirac distribution. The equilibrium state for the electrons is determined by electronegativity equalization with conservation of the total number of electrons. The amount of charge transfer is controlled by re-interpreting the temperature parameter in the Fermi-Dirac distribution function as a coupling strength parameter. We determined this coupling parameter so as to reproduce the charge transfer energy obtained by block localized energy decomposition analysis. We apply the new method to ten systems, and we show that it can yield reasonable approximations to potential energy profiles, to charge transfer stabilization energies, and to the direction and amount of charge transferred.
Collapse
Affiliation(s)
- Miho Isegawa
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | | | | |
Collapse
|
18
|
Lee AJ, Rick SW. The effects of charge transfer on the properties of liquid water. J Chem Phys 2011; 134:184507. [PMID: 21568521 DOI: 10.1063/1.3589419] [Citation(s) in RCA: 88] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Affiliation(s)
- Alexis J Lee
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, USA
| | | |
Collapse
|
19
|
Rana MK, Pazzona FG, Suffritti GB, Demontis P, Masia M. Estimation of Partial Charges in Small Zeolite Imidazolate Frameworks from Density Functional Theory Calculations. J Chem Theory Comput 2011; 7:1575-82. [DOI: 10.1021/ct100685p] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Malay Kumar Rana
- Dipartimento di Chimica, Università degli Studi di Sassari, Via Vienna 2, 07100 Sassari, Italy
| | - Federico Giovanni Pazzona
- Dipartimento di Chimica, Università degli Studi di Sassari, and Consorzio Interuniversitario Nazionale per la Scienza e Tecnologia dei Materiali (INSTM), Unità di Ricerca di Sassari, Via Vienna 2, 07100 Sassari, Italy
| | - Giuseppe Baldovino Suffritti
- Dipartimento di Chimica, Università degli Studi di Sassari, and Consorzio Interuniversitario Nazionale per la Scienza e Tecnologia dei Materiali (INSTM), Unità di Ricerca di Sassari, Via Vienna 2, 07100 Sassari, Italy
| | - Pierfranco Demontis
- Dipartimento di Chimica, Università degli Studi di Sassari, and Consorzio Interuniversitario Nazionale per la Scienza e Tecnologia dei Materiali (INSTM), Unità di Ricerca di Sassari, Via Vienna 2, 07100 Sassari, Italy
| | - Marco Masia
- Dipartimento di Chimica, Università degli Studi di Sassari, Istituto Officina dei Materiali del CNR, UOS SLACS, Via Vienna 2, 07100 Sassari, Italy
| |
Collapse
|
20
|
Mukherjee G, Patra N, Barua P, Jayaram B. A fast empirical GAFF compatible partial atomic charge assignment scheme for modeling interactions of small molecules with biomolecular targets. J Comput Chem 2010; 32:893-907. [DOI: 10.1002/jcc.21671] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2009] [Revised: 07/17/2010] [Accepted: 08/12/2010] [Indexed: 11/07/2022]
|
21
|
From Molecules to Droplets. ADVANCES IN QUANTUM CHEMISTRY 2008. [DOI: 10.1016/s0065-3276(07)00217-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register]
|
22
|
Buttingsrud B, Alsberg BK, Astrand PO. Validation of critical points in the electron density as descriptors by building quantitative structure-property relationships for the atomic polar tensor. J Comput Chem 2007; 28:2130-9. [PMID: 17464968 DOI: 10.1002/jcc.20666] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
A crucial component of research in the field of quantitative structure-activity/property relationships is the identification of molecular descriptors relevant to the activity or property of interest. Descriptors based on the topology of the electron density as formulated in Bader's theory of atoms in molecules are investigated in detail in this work. In a model study, the authors investigate their ability to predict the atomic polar tensor (the gradient of the molecular dipole moment), which contains information on the vibrational intensities in infrared spectroscopy and constitutes a scheme for partitioning the total charge distribution into atomic charges. The atomic polar tensor may therefore be used to investigate whether the descriptors give adequate information on the local electronic structure in the molecule. Both the trace of the atomic polar tensor and for planar molecules its out-of-plane component may be interpreted as definitions of atomic charges suitable for prediction. Hydrogen and carbon atoms in a set of 60 aromatic compounds with various substituents have been studied. Excellent results for prediction of hydrogen and carbon charges have been achieved with cross-validated squared correlation coefficients between predicted and theoretical values varying from 0.92 and 0.977 for the most complex set of substituents when the value, Laplacian, and ellipticity of the electron density in the bond critical points are used as descriptors. The carbon charges defined from the trace of the atomic polar tensor are correlated with its out-of-plane component whereas such relationship is not observed for the hydrogen charges studied in this work.
Collapse
Affiliation(s)
- Bård Buttingsrud
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), 7491 Trondheim, Norway
| | | | | |
Collapse
|
23
|
Mrázek J, Burda JV. Can the pH value of water solutions be estimated by quantum chemical calculations of small water clusters? J Chem Phys 2006; 125:194518. [PMID: 17129134 DOI: 10.1063/1.2363383] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In the study, various water clusters were explored from the point of view of the proton transfer between H-bonded neighbors. A relatively modest approach--the MP26-31++G(d,p) level--was chosen as acceptable considering the fact that also larger systems have to be included. The tight ion-pair model (with usually three fixed O-O distances) was adopted for the autodissociation process. First, cluster-estimated pH values rapidly decrease as cluster size increases from 2 to 6. For larger clusters in gas phase, the topology of H bonds plays an important role, varying pH from 7 to 13 in hexamers and from 5 to 15 in octamer clusters. The relationship energy/distance was quantified, too. Enhancing our model with the conductorlike screening model (COSMO) approach brought significant improvement in description of the autodissociation reaction with a stable zwitterionic structure. However, when the O-O restrictions were released, the small barrier for backward reaction disappeared, reforming neutral cluster spontaneously. Also Boltzmann weighting procedure was applied, and for the explored clusters in vacuo, the series of pH 25-18-14-13-10 was obtained for cluster sizes n=2, 3, 4, 6, 8. Using the COSMO approach, the analogous series is 15-14-12-10-9. The limit of the series is still about two to three units above the experimentally known pH. In order to reach the size-independent (bulk) value, larger clusters are needed. However, the situation is far from hopeless since (as it was proven in the study) four-coordinated molecules are not involved in the proton transfer process directly; they can only be a part of the surrounding environment.
Collapse
Affiliation(s)
- Jirí Mrázek
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 121 16 Prague 2, Czech Republic
| | | |
Collapse
|
24
|
Mikkelsen KV. CORRELATED ELECTRONIC STRUCTURE NONLINEAR RESPONSE METHODS FOR STRUCTURED ENVIRONMENTS. Annu Rev Phys Chem 2006; 57:365-402. [PMID: 16599815 DOI: 10.1146/annurev.physchem.57.032905.104740] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
▪ Abstract This contribution concerns a brief outline of structural environment models where correlated electronic structure response methods are utilized for the determination of nonlinear optical properties of molecules. The presentation provides theory and applications of a heterogeneous dielectric media model and a quantum mechanical–classical mechanical model at the level of correlated electronic structure response methods. The correlated electronic structure response methods include the multiconfigurational self-consistent field method.
Collapse
Affiliation(s)
- Kurt V Mikkelsen
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen, Denmark.
| |
Collapse
|
25
|
|
26
|
Abstract
The charge distribution in the water molecule has been analyzed using a broad variety of basis sets, four different quantum mechanical methods (Hartree-Fock, Becke3LYP, MP2, and QCISD), and six population analysis methods (Mulliken, NPA, AIM, CHELPG, Merz-Kollman, and Resp). The influence of the molecular structure on the calculated atomic charges has been studied using small perturbations of the experimentally determined structure.
Collapse
Affiliation(s)
- F Martin
- Department Chemie, LMU München, Butenandstr. 13, D-81377 München, Germany
| | | |
Collapse
|
27
|
Solheim H, Ruud K, Astrand PO. Atomic dipole moments calculated using analytical molecular second-moment gradients. J Chem Phys 2004; 120:10368-78. [PMID: 15268064 DOI: 10.1063/1.1738631] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We have implemented analytical second-moment gradients for Hartree-Fock and multiconfigurational self-consistent-field wave functions. The code is used to calculate atomic dipole moments based on the generalized atomic polar tensor (GAPT) formalism [Phys. Rev. Lett. 62, 1469 (1989)], and the proposal of Dinur and Hagler (DH) for the calculation of atomic multipoles [J. Chem. Phys. 91, 2949 (1989)]. Both approaches display smooth basis-set convergence toward a well-defined basis-set limit and give reasonable electron correlation effects on the calculated atomic properties. However, the atomic charges and atomic dipole moments obtained from the GAPT partitioning scheme are unable to provide even qualitatively meaningful molecular quadrupole moments for some molecules, and thus the atomic multipole moments calculated in this scheme cannot be considered well suited for analyzing the electron density in molecules and for calculating intermolecular interaction energies. In contrast, the DH approach gives atomic charges and dipole moments that by definition exactly reproduce the molecular quadrupole moments. The approach of DH is, however, restricted to planar molecules and thus suffers from not being applicable to molecules of arbitrary shape. Both the GAPT and DH approaches give rather poor results for octupole and hexadecapole moments, indicating that at least atomic quadrupole moments are required for an accurate representation of the molecular charge distribution in terms of atomic electric moments.
Collapse
Affiliation(s)
- Harald Solheim
- Department of Chemistry, University of Tromsø, N-9037 Tromsø, Norway
| | | | | |
Collapse
|
28
|
Fonseca Guerra C, Handgraaf JW, Baerends EJ, Bickelhaupt FM. Voronoi deformation density (VDD) charges: Assessment of the Mulliken, Bader, Hirshfeld, Weinhold, and VDD methods for charge analysis. J Comput Chem 2004; 25:189-210. [PMID: 14648618 DOI: 10.1002/jcc.10351] [Citation(s) in RCA: 654] [Impact Index Per Article: 32.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
We present the Voronoi Deformation Density (VDD) method for computing atomic charges. The VDD method does not explicitly use the basis functions but calculates the amount of electronic density that flows to or from a certain atom due to bond formation by spatial integration of the deformation density over the atomic Voronoi cell. We compare our method to the well-known Mulliken, Hirshfeld, Bader, and Weinhold [Natural Population Analysis (NPA)] charges for a variety of biological, organic, and inorganic molecules. The Mulliken charges are (again) shown to be useless due to heavy basis set dependency, and the Bader charges (and often also the NPA charges) are not realistic, yielding too extreme values that suggest much ionic character even in the case of covalent bonds. The Hirshfeld and VDD charges, which prove to be numerically very similar, are to be recommended because they yield chemically meaningful charges. We stress the need to use spatial integration over an atomic domain to get rid of basis set dependency, and the need to integrate the deformation density in order to obtain a realistic picture of the charge rearrangement upon bonding. An asset of the VDD charges is the transparency of the approach owing to the simple geometric partitioning of space. The deformation density based charges prove to conform to chemical experience.
Collapse
Affiliation(s)
- Célia Fonseca Guerra
- Afdeling Theoretische Chemie, Scheikundig Laboratorium der Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| | | | | | | |
Collapse
|
29
|
Yang ZZ, Wu Y, Zhao DX. Atom-bond electronegativity equalization method fused into molecular mechanics. I. A seven-site fluctuating charge and flexible body water potential function for water clusters. J Chem Phys 2004; 120:2541-57. [PMID: 15268398 DOI: 10.1063/1.1640345] [Citation(s) in RCA: 83] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Recently, experimental and theoretical studies on the water system are very active and noticeable. A transferable intermolecular potential seven points approach including fluctuation charges and flexible body (ABEEM-7P) based on a combination of the atom-bond electronegativity equalization and molecular mechanics (ABEEM/MM), and its application to small water clusters are explored and tested in this paper. The consistent combination of ABEEM and molecular mechanics (MM) is to take the ABEEM charges of atoms, bonds, and lone-pair electrons into the intermolecular electrostatic interaction term in molecular mechanics. To examine the charge transfer we have used two models coming from the charge constraint types: one is a charge neutrality constraint on whole water system and the other is on each water molecule. Compared with previous water force fields, the ABEEM-7P model has two characters: (1) the ABEEM-7P model not only presents the electrostatic interaction of atoms, bonds and lone-pair electrons and their changing in respond to different ambient environment but also introduces "the hydrogen bond interaction region" in which a new parameter k(lp,H)(R(lp,H)) is used to describe the electrostatic interaction of the lone-pair electron and the hydrogen atom which can form the hydrogen bond; (2) nonrigid but flexible water body permitting the vibration of the bond length and angle is allowed due to the combination of ABEEM and molecular mechanics, and for van der Waals interaction the ABEEM-7P model takes an all atom-atom interaction, i.e., oxygen-oxygen, hydrogen-hydrogen, oxygen-hydrogen interaction into account. The ABEEM-7P model based on ABEEM/MM gives quite accurate predictions for gas-phase state properties of the small water clusters (H(2)O)(n) (n=2-6), such as optimized geometries, monomer dipole moments, vibrational frequencies, and cluster interaction energies. Due to its explicit description of charges and the hydrogen bond, the ABEEM-7P model will be applied to discuss properties of liquid water, ice, aqueous solutions, and biological systems.
Collapse
Affiliation(s)
- Zhong-Zhi Yang
- Department of Chemistry, Liaoning Normal University, Dalian 116029, People's Republic of China
| | | | | |
Collapse
|
30
|
Guillot B. A reappraisal of what we have learnt during three decades of computer simulations on water. J Mol Liq 2002. [DOI: 10.1016/s0167-7322(02)00094-6] [Citation(s) in RCA: 757] [Impact Index Per Article: 34.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
31
|
Jørgensen S, Ratner MA, Mikkelsen KV. Heterogeneous solvation: An ab initio approach. J Chem Phys 2001. [DOI: 10.1063/1.1387979] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
32
|
Kumar A, Mishra P. Study of hybridization displacement charge distribution using ab initio wave functions and evaluation of molecular electrostatic properties. ACTA ACUST UNITED AC 2001. [DOI: 10.1016/s0166-1280(00)00858-7] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
33
|
|
34
|
Rick SW. Simulations of ice and liquid water over a range of temperatures using the fluctuating charge model. J Chem Phys 2001. [DOI: 10.1063/1.1336805] [Citation(s) in RCA: 162] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
35
|
Kar T, Ángyán JG, Sannigrahi AB. Comparison of ab Initio Hartree−Fock and Kohn−Sham Orbitals in the Calculation of Atomic Charge, Bond Index, and Valence. J Phys Chem A 2000. [DOI: 10.1021/jp0008694] [Citation(s) in RCA: 80] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Tapas Kar
- Department of Chemistry, Southern Illinois University, Carbondale, Illinois 62901
| | - János G. Ángyán
- Laboratoire de Chimie Théorique, UMR CNRS No. 7565, Institut Nancéien de Chimie Moléculaire, Université Henri Poincaré, B.P. 239, 54506 Vandoeuvre-lés-Nancy, France
| | - A. B. Sannigrahi
- Department of Chemistry, Indian Institute of Technology, Kharagpur 721 302, India
| |
Collapse
|
36
|
Grigorenko BL, Nemukhin AV, Topol IA, Burt SK. Hydrogen bonding at the diatomics-in-molecules level: Water clusters. J Chem Phys 2000. [DOI: 10.1063/1.1303850] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
|
37
|
Maroulis G. Static hyperpolarizability of the water dimer and the interaction hyperpolarizability of two water molecules. J Chem Phys 2000. [DOI: 10.1063/1.481985] [Citation(s) in RCA: 129] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
|
38
|
de Oliveira AE, Haiduke RLA, Bruns RE. Atomic Mean Dipole Moment Derivatives and GAPT Charges. J Phys Chem A 2000. [DOI: 10.1021/jp994405l] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Anselmo E. de Oliveira
- Instituto de Química, Universidade Estadual de Campinas, CP 6154, 13083-970 Campinas, SP, Brazil
| | - Roberto L. A. Haiduke
- Instituto de Química, Universidade Estadual de Campinas, CP 6154, 13083-970 Campinas, SP, Brazil
| | - Roy E. Bruns
- Instituto de Química, Universidade Estadual de Campinas, CP 6154, 13083-970 Campinas, SP, Brazil
| |
Collapse
|
39
|
Hermida Ramón JM, Rı́os MA. A new intermolecular polarizable potential for cis-formic acid. Introduction of many-body interactions in condensed phases. Chem Phys 1999. [DOI: 10.1016/s0301-0104(99)00318-3] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
40
|
Ribeiro MCC, Almeida LCJ. Fluctuating charge model for polyatomic ionic systems: A test case with diatomic anions. J Chem Phys 1999. [DOI: 10.1063/1.479085] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
41
|
Bryce RA, Vincent MA, Hillier IH. Binding Energy of F(H2O)- and the Simulation of Fluoride Water Clusters Using a Hybrid QM/MM (Fluctuating Charge) Potential. J Phys Chem A 1999. [DOI: 10.1021/jp9900990] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Richard A. Bryce
- Department of Chemistry, University of Manchester, Manchester M13 9PL, U.K
| | - Mark A. Vincent
- Department of Chemistry, University of Manchester, Manchester M13 9PL, U.K
| | - Ian H. Hillier
- Department of Chemistry, University of Manchester, Manchester M13 9PL, U.K
| |
Collapse
|