The fundamental nature and role of the electrostatic potential in atoms and molecules

 A variety of atomic and molecular properties can be expressed in terms of the electrostatic potential. These include energies, covalent and anionic radii, electronegativities (chemical potentials) and a variety of properties that depend upon noncovalent interactons. We present a survey of such relationships, which may be exact or approximate; they may involve the potential in three-dimensional space, along the axes between bonded atoms, at nuclei or on molecular surfaces. Thus, the electrostatic potential, which is rigorously related to the electronic density by Poisson's equation, can be regarded as, effectively, another fundamental determinant of atomic and molecular properties. Received: 6 March 2002 / Accepted: 15 May 2002 / Published online: 29 July 2002     

The Antiepileptic Lamotrigine and Its Analogues: Comparative Theoretical Electronic Properties

Electronic properties of lamotrigine (LTG) and two analogues (A1 and A2) are compared through MOPAC-AM1 calculations. Two stable conformers of LTG are calculated to exist in agreement with X-ray crystallography. In the three compounds and the two conformers for each of them, the more favorable protonation sites are N₂ and N₄; these should then be the sites appropriate for interaction with a receptor, and group valence reinforces the supposition. The molecular electrostatic potentials show that a region between the two chlorine atoms in LTG could be the site for an electrostatic interaction with a corresponding site in the receptor. The fluorine atom in A1 would play an equivalent role. A simple model for LTG-receptor interaction is proposed.     

Nuclear Coordinate dependence of electronic transition moments in orbitally forbidden transitions: A new interpretation of deuterium isotope effects

The nuclear coordinate dependence of electronic transtion moments has been investigated for the purpose of finding new interpretations of deuterium isotope effects on spectral intensities and radiative decay rates in orbitally forbidden electronic transitions. By using “AO following nuclei” wavefunctions as the building block for the electronic wavefunction in the adiabatic BO vibronic wavefunction, the spin-free hamiltonian is diagonalized to generate eigenfunctions and eigen-energies. It is found that the electronic transtion moments based on these eigenfunctions show dependences upon the vibrational modes which are not directly involved in vibronic coupling. This leads to interpretations of the deuterium isotope effects in T₁ → S₀ radiative transitions of aromatic hydrocarbons and S₀ → S₁ absorption in pyrazine which are not based on the conventional Herzberg—Teller or non-BO coupling.     

Tuning the polarization along linear polyaromatic strands for rationally inducing mesomorphism in lanthanide nitrate complexes

The opposite orientation of the ester spacers in the rodlike ligands L 4C12 (benzimidazole-OOC-phenyl) and L 5C12 (benzimidazole-COO-phenyl) drastically changes the electronic structure of the aromatic systems, without affecting their meridional tricoordination to trivalent lanthanides, Ln(III), and their thermotropic liquid crystalline (i.e., mesomorphic) behaviors. However, the rich mesomorphism exhibited by the complexes [Ln(L 4C12)(NO3)3] (Ln=La-Lu) vanishes in [Ln(L 5C12)(NO3)3], despite superimposable molecular structures and comparable photophysical properties. Density functional theory (DFT) and time-dependant DFT calculations performed in the gas phase show that the inversion of the ester spacers has considerable effects on the electronic structure and polarization of the aromatic groups along the strands, which control residual intermolecular interactions responsible for the formation of thermotropic liquid-crystalline phases. As a rule of thumb, an alternation of electron-poor and electron-rich aromatic rings favors intermolecular interactions between the rigid cores and consequently mesomorphism, a situation encountered for L 4C12, L 5C12, [Ln(L 4C12)(NO3)3], but not for [Ln(L 5C12)(NO3)3]. The intercalation of an additional electron-rich diphenol ring on going from [Ln(L 5C12)(NO3)3] to [Ln(L 6C12)(NO3)3] restores mesomorphism despite an unfavorable orientation of the ester spacers, in agreement with our simple predictive model.     

Corannulenes with Electron-Withdrawing Substituents: Synthetic Approaches and Resulting Structural and Electronic Properties

Corannulene is a multifaceted polyaromatic compound. It has many interesting properties; for example, it has a bowl-shaped molecular structure that, in addition, undergoes a dynamic inversion process. It has attracted much attention within the last decades. This is not only due to its structural properties but also its electronic properties and its various potential applications to materials chemistry. Here, synthetic approaches towards corannulene derivatives with electron-withdrawing substituents are summarized. This includes both selective and unselective methods. Further, the electrochemical properties, that is, the reduction potentials, are analyzed and compared. As a main conclusion, one can state that the electron affinity depends roughly linearly on the number of substituents. Finally, the structural behavior of the substituted buckybowls in the solid state is highlighted. This also allows a general statement about the influence of the electronic and steric nature of substituents on the molecular structures and the solid-state packing of the corannulene derivatives.     

Molecular structure calculations without clamping the nuclei

A number of recent papers have considered ways in which molecular structure may be calculated when both the electrons and the nuclei are treated from the outset as quantum particles. This is in contrast to the conventional approach in which the nuclei initially have their positions fixed and so merely provide a potential for electronic motion. The usual approach is generally assumed to be justified by the 1927 work of Born and Oppenheimer. In this paper we discuss what precisely might be anticipated in the way of molecular structure from a mathematical consideration of the spectral properties of the full Coulomb Hamiltonian, to what extent the very idea of molecular structure might be dependent upon treating the nuclei simply as providing a potential and the extent to which the work of Born and Oppenheimer can be used to support this position.     

Relevance between the Chemical Structure Modifications and Physicochemical Descriptors of Chemical Reactivity for Series of Nitroxide Radicals

This work presents an experimental and theoretical study to address the chemical reactivity of series of nitroxide radicals. For that purpose two physicochemical properties: the half-wave potential and the hyperfine coupling constants of the nitrogen nuclei, were analyzed. Experimental values are compared with electronic structure calculations at the BHandHLYP/6-311++G(2d,2p) level. E_1/2 values were in good agreement with the adiabatic ionization potential when including the solvent effects by the Cramer and Truhlar Solvation Model. Preeliminar experimental electron spin deslocalization studies suggest that structural hindrance plays an important role in their deslocatization mechanism.     

Building photoactive molecular-scale wires

Polyalkynylene groups are known to function as excellent electronic conductors at the molecular level. Such moieties have now been used to interconnect redox and photoactive transition metal oligopyridine complexes so that the efficiency of light-induced energy or electron transfer along the molecular axis can be monitored. The important issues that control the effectiveness of electronic coupling through the alkyne are discussed. In particular, attention is given to separating the effects of electron delocalization within the triplet manifold from the more general decoupling of metal-centered and charge-transfer excited states that occurs upon lowering the triplet energy. The role of the auxiliary ligands is considered, as is the effect of nuclearity. Similarly, the size of the nuclear reorganization energy has to be taken into account in a proper discussion of the photophysical properties of such systems. A second issue of importance to the design of photoelectronic devices concerns the use of interspersed groups to modulate the electronic coupling properties of the alkyne spacer. Such electron relays may be aryl hydrocarbons or platinum bis-acetylides, both groups being able to curtail electron flow along the molecular axis.     

Can Counter-Intuitive Halogen Bonding Be Coulombic?

We use the term “counter-intuitive” to describe an intermolecular interaction in which the electrostatic potentials of the interacting regions of the ground-state molecules have the same sign, both positive or both negative. In the present work, we consider counter-intuitive halogen bonding with nitrogen bases, in which both the halogen σ-hole and the nitrogen lone pair have negative potentials on their molecular surfaces. We show that these interactions can be treated as Coulombic despite the apparent repulsion between the ground-state molecules, provided that both electrostatics and polarization are explicitly taken into account. We demonstrate first that the energies of 20 counter-intuitive interactions with four nitrogen bases can be expressed very well in terms of just two molecular properties: the electrostatic potential of the halogen σ-hole and the average polarizability of the nitrogen base. Then we show that the same two properties can also represent the energies of an expanded data base that includes the 20 counter-intuitive plus an additional 20 weak and moderately-strong intuitive halogen bonding interactions (in which the σ-hole potentials are now positive).     

What mathematicians know about the solutions of Schrodinger Coulomb Hamiltonian. Should chemists care?

An attempt is made to determine the relationship between the full Schrödinger Coulomb Hamiltonian and the clamped nuclei form that is usually the basis of electronic structure calculations, without treating identical nuclei as distinguishable. It is concluded that it is not at present possible to establish such a relationship in a mathematically secure way.     

Quantum chemical study on the interaction of some bisphosphonates and Ca2+: The role of molecular electrostatic potentials in the prediction of binding geometry

Summary Molecular electrostatic potentials have been used to model the calcium binding properties of some bisphosphonate drugs, which are used to treat various bone diseases. The mechanism of action involves the binding of bisphosphonates to the bone surface, where calcium plays an important role. Electrostatic potential maps derived from ab initio partial charges have been compared with both the crystal structure and the fully optimized ab initio structure of (dichloro)methylenebisphosphonate-calcium ion complex. Molecular electrostatic potentials can correctly predict the calcium binding geometry of bisphosphonate-type compounds and this type of information can be used in the practical drug design work.     

Guest Binding with Sulfated Cyclodextrins: Does the Size of Cavity Matter?

Sulfated cyclodextrins have recently emerged as potential candidates for producing host–induced guest aggregation with properties better than p-sulfonatocalixarenes that have previously shown numerous applications involving the phenomena of host-induced guest aggregation. In the class of sulfated cyclodextrins (SCD), sulfated β-cyclodextrin (β-SCD) remains the most extensively investigated host molecule. Although it is assumed that the host-induced guest aggregation is predominantly an outcome of interaction of the guest molecule with the charges on the exterior of SCD cavity, it has not been deciphered whether the variation in the cavity size will make a difference in the efficiency of host-induced guest-aggregation process. In this investigation, we present a systematic study of host–induced guest aggregation of a cationic molecular rotor dye, Thioflavin T (ThT) with three different sulfated cyclodextrin molecules, α-SCD, β-SCD and γ-SCD, which differ in their cavity size, using steady-state emission, ground-state absorption and time-resolved emission measurements. The obtained photophysical properties of ThT, upon interaction with different SCD molecules, indicate that the binding strength of ThT with different SCD molecules correlate with the cavity size of the host molecule, giving rise to the strongest complexation of ThT with the largest host molecule (γ-SCD). The binding affinity of ThT towards different host molecules has been supported by molecular docking calculations. The results obtained are further supported with the temperature and ionic strength dependent studies performed on the host-guest complex. Our results indicate that for host–induced guest aggregation, involving oppositely charged molecules, the size of the cavity also plays a crucial role beside the charge density on the exterior of host cavity.     

Computational investigation of the role of counterions and reorganization energy in a switchable bistable [2]rotaxane

Switchable bistable [2]rotaxanes, such as those of the Stoddart-Heath-type, show promise for the development of molecular electronic devices and functional prototypes have been demonstrated. Herein, one such switchable rotaxane system is studied computationally at the AM1-FS1 and DFT levels of theory. The results show that the computationally efficient AM1-FS1 method, (efficient relative to DFT) is capable of reliably predicting properties such as binding site preference and coconformational relative stabilities as well as the barrier to isomerization between the different coconformational states. These properties play important roles in the functionality of rotaxane-based molecular electronic devices. In addition, the role of the counterions is assessed from a computational standpoint. The results reveal that inclusion of counterions is not as significant as has been previously suggested. Finally, the reorganization energy associated with oxidation/reduction of the complex is studied. This provides a possible link to the origin of the observed conductivity difference between the two coconformational states, the property upon which device functionality is based.     

Tunneling currents that increase with molecular elongation

We present a model molecular system with an unintuitive transport-extension behavior in which the tunneling current increases with forced molecular elongation. The molecule consists of two complementary aromatic units (1,4-anthracenedione and 1,4-anthracenediol) hinged via two ether chains and attached to gold electrodes through thiol-terminated alkenes. The transport properties of the molecule as it is mechanically elongated in a single-molecule pulling setting are computationally investigated using a combination of equilibrium molecular dynamics simulations of the pulling with gDFTB computations of the transport properties in the Landauer limit. Contrary to the usual exponential decay of tunneling currents with increasing molecular length, the simulations indicate that upon elongation electronic transport along the molecule increases 10-fold. The structural origin of this inverted trend in the transport is elucidated via a local current analysis that reveals the dual role played by H-bonds in both stabilizing π-stacking for selected extensions and introducing additional electronic couplings between the complementary aromatic rings that also enhance tunneling currents across the molecule. The simulations illustrate an inverted electromechanical single-molecule switch that is based on a novel class of transport-extension behavior that can be achieved via mechanical manipulation and highlight the remarkable sensitivity of conductance measurements to the molecular conformation.     

Stereoelectronic properties of spiroquinazolinones in differential PDE7 inhibitory activity

A detailed computational study on a series of spiroquinazolinones showing phosphodiesterase 7 (PDE7) inhibitory activity was performed to understand the binding mode and the role of stereoelectronic properties in binding. Our docking studies reproduced the essential hydrogen bonding and hydrophobic interactions for inhibitors of this class of enzymes. The N1 proton of the quinazolinone scaffold was involved in H-bonding to an amide side chain of the conserved glutamine residue in the active site. The central bicyclic ring of the molecules showed hydrophobic and pi-stacking interactions with hydrophobic and aromatic amino acid residues, respectively, present in the PDE7 active site. The docked conformations were optimized with density functional theory (DFT) and DFT electronic properties were calculated. Comparison of molecular electrostatic potential (MEP) plots of inhibitors with the active site of PDE7 suggested that the electronic distribution in the molecules is as important as steric factors for binding of the molecules to the receptor. The hydrogen bonding ability and nucleophilic nature of N1 appeared to be important for governing the interaction with PDE7. For less active inhibitors (pIC(50) < 6.5), the MEP maximum at N1 of the spiroquinazolinone ring was high or low based on the electronic properties of the substituents. All the more active molecules (pIC(50) > 6.5) had MEP highest at N3, not N1. Efficient binding of these inhibitors may need some rearrangement of side chains of active-site residues, especially Asn365. This computational modeling study should aid in design of new molecules in this class with improved PDE7 inhibition.     

Nanoscale structural and electronic properties of ultrathin blends of two polyaromatic molecules: a Kelvin probe force microscopy investigation.

We describe a Kelvin Probe Force Microscopy (KPFM) study on the morphological and electronic properties of complex mono and bi-molecular ultrathin films self-assembled on mica. These architectures are made up from an electron-donor (D), a synthetic all-benzenoid polycyclic aromatic hydrocarbon, and an electron-acceptor (A), perylene-bis-dicarboximide. The former molecule self-assembles into fibers in single component films, while the latter molecule forms discontinuous layers. Taking advantage of the different solubility and self-organizing properties of the A and D molecules, multicomponent ultrathin films characterized by nanoscale phase segregated fibers of D embedded in a discontinuous layer of A are formed. The direct estimation of the surface potential, and consequently the local workfunction from KPFM images allow a comparison of the local electronic properties of the blend with those of the monocomponent films. A change in the average workfunction values of the A and D nanostructures in the blend occurs which is primarily caused by the intimate contact between the two components and the molecular order within the nanostructure self-assembled at the surface. Additional roles can be ascribed to the molecular packing density, to the presence of defects in the film, to the different conformation of the aliphatic peripheral chains that might cover the conjugated core and to the long-range nature of the electrostatic interactions employed to map the surface by KPFM limiting the spatial and potential resolution. The local workfunction studies of heterojunctions can be of help to tune the electronic properties of active multicomponent films, which is crucial for the fabrication of efficient organic electronic devices as solar cells.     

Topography of approximate molecular electrostatic potentials

A new method is described for the approximation of the molecular electrostatic potential (MESP). This method is used for the study of the topography of small molecules. The critical points of the approximate and the exact MESP are compared. It is found that most of the critical points of the exact MESP are retained, but in regions where the exact MESP changes slowly near critical points the number of critical points of the approximate MESP can be reduced.