Predicting reduction potentials of 1,3,6‐triphenyl fulvenes using molecular electrostatic potential analysis of substituent effects

The influence of mono‐ and multiple substituent effect on the reduction potential (E⁰) of 1,3,6‐triphenyl fulvenes is investigated using B3LYP‐SMD/6‐311+G(d,p) level density functional theory. The molecular electrostatic potential (MESP) minimum at the fulvene π‐system (V_min) and the change in MESP at any of the fulvene carbon atoms (ΔV_C) for both neutral and reduced forms are used as excellent measures of substituent effect from the para and meta positions of the 1,3 and 6‐phenyl moieties. Substitution at 6‐phenyl para position has led to significant change in E⁰ than any other positions. By applying the additivity rule of substituent effects, an equation in ΔV_C is derived to predict E⁰ for multiply substituted fulvenes. Further, E⁰ is predicted for a set of 2000 hexa‐substituted fulvene derivatives where the substituents and their positions in the system are chosen in a random way. The calculated E⁰ agreed very well with the experimental E⁰ reported by Godman et al. Predicting E⁰ solely by substituent effect offers a simple and powerful way to select suitable combinations of substituents on fulvene system for light harvesting applications. © 2018 Wiley Periodicals, Inc.     

密度泛函理论和从头算方法对富勒烯分子静电势的比较研究

分别在Hartree-Fock和密度泛函B3LYP理论下,用6-31G*基组研究了C60和C70分子的静电势,比较了这方法计算得到上述分子静电势值的大小,静电势图形和静电势差值曲线,分析了富勒烯的电子相关效应.     

A new method to determine electrostatic potential around a macromolecule in solution from molecular wave functions

The three-dimensional reference interaction site model integral equation theory (3D-RISM) combined with the ab initio molecular orbital method (3D-RISM-SCF) is applied to a solvated macromolecular system. The solvation structure around a solute molecule is obtained from the 3D-RISM integral equation under the electrostatic potential of the solute molecule, calculated by the ab initio molecular orbital theory. The electrostatic potential should be calculated on each grid point in the three-dimensional real space. Therefore, the calculation of the electrostatic potential is the most time consuming part in this method. In this article, we propose a new procedure to save the computational cost for calculating the electrostatic potential and the solvated fock matrix. The strategy of this procedure is to evaluate the electrostatic potential and the solvated fock matrix in different ways, depending on the distance between solute and solvent. Inside the repulsive cores of solute atoms, it is possible to avoid the calculation of electrostatic potential and solvated Fock matrix by assuming the potential to be infinity. In the region sufficiently far from solute, they are evaluated classically by putting the effective point charge on each atom. In the intermediate region, the electrostatic potential is evaluated directly by integrating the molecular orbitals of the solute molecule. The electronic structure and the energy gradient of Methionine-Enkephalin and solvation structure are estimated by using this procedure in aqueous solution, and are compared with the results from other procedures. The results are compared also with those from the continuum model.     

The electrostatical molecular potential — A tool for the prediction of electrostatic molecular interaction properties

Starting from the molecular potential we get, by using elementary electrostatics, information about energetically favoured regions for interaction with ions and dipoles around H₂O and H₂CO. The molecule-dipole interaction is represented by the electric field patterns.     

Atomic and molecular energies as functionals of the electrostatic potential

Starting from the Hohenberg–Kohn theorem, atomic and molecular energies have been expressed rigorously as functionals of the electronic electrostatic potential, V_elec(r). Explicit formulations have been derived for the functionals representing the kinetic energy and electronic interaction contributions to the total energies.Acknowledgements. The assistance of Dr. Jane S. Murray is greatly appreciated.Contribution to the Jacopo Tomasi Honorary Issue     

Prediction of densities for solid energetic molecules with molecular surface electrostatic potentials

The densities of high energetic molecules in the solid state were calculated with a simplified scheme based on molecular surface electrostatic potentials (MSEP). The MSEP scheme for density estimation, originally developed by Politzer et al., was further modified to calculate electrostatic potential on a simpler van der Waals surface. Forty-one energetic molecules containing at least one nitro group were selected from among a variety of molecular types and density values, and were used to test the suitability of the MSEP scheme for predicting the densities of solid energetic molecules. For comparison purposes, we utilized the group additivity method (GAM) incorporating the parameter sets developed by Stine (Stine-81) and by Ammon (Ammon-98 and -00). The absolute average error in densities from our MSEP scheme was 0.039 g/cc. The results based on our MSEP scheme were slightly better than the GAM results. In addition, the errors in densities generated by the MSEP scheme were almost the same for various molecule types, while those predicted by GAM were somewhat dependent upon the molecule types.     

富勒烯结构B_(80)分子静电势的理论研究

在混合密度泛函B3LYP理论下,用6-31G*基函数对富勒烯结构B80分子的3个异构体(1个具有Ih对称性,2个具有Th对称性)构型进行优化和分子静电势计算.结果表明:3个异构体球内全部为正电势,球外五元环中心所对应的区域都为负电势,B80Ih,Th(A)和Th(B)球外静电势的最大负值分别对应于20个六元环中心的B原子,五元环中心和12个六元环中心的B原子周围,它们组成了化学反应中最可能的活性点.     

Computed electrostatic potentials and average local ionization energies on the molecular surfaces of some tetracyclines

The electrostatic potentials and average local ionization energies computed on the molecular surfaces of four tetracyclines have been investigated with the objective of identifying common features as well as subtle differences that may be related to their biological activities. The four are the parent molecule tetracycline, chlortetracycline, oxytetracycline, and doxycyline. The calculations were carried out at the HF/STO‐3G*//STO‐3G* level. Our electrostatic potential results show that each molecule has a large negative region that extends along its lower portion, consistent with its ability to complex Ca²⁺ and Mg²⁺ ions. Although the surface electrostatic potentials of the four tetracyclines show many similarities, our statistical measure of local polarity allows us to label doxycycline as the one with the lowest degree of local polarity, consistent with its longer half‐life in vivo. The regions in the tetracyclines with the most reactive electrons are the amide nitrogen lone pairs and certain carbons of the outermost rings. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 160–169, 2001     

Near ab initio quality molecular electrostatic potential maps using hybridization displacement charges at PM3 level and effects of geometrical changes in amino groups

Charge distributions, dipole moments, and molecular electrostatic potentials (MEP) around several molecules consisting of carbon, nitrogen, oxygen, fluorine, sulfur, and chlorine atoms were studied using the PM3 semiempirical method and the results compared with those obtained using ab initio calculations at the RHF/6‐31G** level. Thus it is shown that relative MEP values near different atoms can be obtained using hybridization displacement charges (HDC) obtained by employing PM3 density matrices that usually agree quite satisfactorily with the ab initio ones. Further, positions of ab initio MEP minima are correctly located and the corresponding relative MEP values usually correctly predicted using the PM3(HDC) charges distributed continuously in three dimensions according to the forms of squares of valence s atomic orbitals. The necessary parameters for HDC calculations using the PM3 method were optimized. It is shown how within the frameworks of both PM3 and AM1 methods the π electrons or lone pairs associated with amino group nitrogen atoms and ring atoms can be satisfactorily treated in different situations. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 299–312, 2001     

The molecular electrostatic potential and steric accessibility of C-DNA

The molecular electrostatic potentials and steric accessibilities associated with reactive sites of C-DNA are calculated for the sequences poly(dG · dC) and poly(dA · dT). The distribution of potential on the surface envelopes of the double helices are also presented. The results are compared with those obtained for B-DNA.