Large Electrostatic Interaction Research Articles

Multifacets of Fullerene-Metal Clusters: From Fundamental to Application.

ConspectusBuckminsterfullerene, C60, was discovered through a prominent mass peak containing 60 atoms produced from laser vaporization of graphite, driven by Kroto's interest in understanding the formation mechanisms of carbon-containing molecules in space. Inspired by the geodesic dome-shaped architecture designed by Richard Buckminster Fuller, after whom the particle was named, C60 was found to have a football-shaped structure comprising 20 hexagons and 12 pentagons. It sparked worldwide interest in understanding this new carbon allotrope, resulting in the awarding of the Noble Prize in Chemistry to Smalley, Kroto, and Curl in 1996.Intrinsically, C60 is an exceptional species because of its high stability and electron-accepting ability and its structural tunability by decorating or substituting either on its exterior surface or interior hollow cavity. For example, metal-decorated fullerene complexes have found important applications ranging from superconductivity, nanoscale electronic devices, and organic photovoltaic cells to catalysis and biomedicine. Compared to the large body of studies on atoms and molecules encapsulated by C60, studies on the exteriorly modified fullerenes, i.e., exohedral fullerenes, are scarcer. Surprisingly, to date, uncertainty exists about a fundamental question: what is the preferable exterior binding site of different kinds of single atoms on the C60 surface?In recent years, we have developed an experimental protocol to synthesize the desired fullerene-metal clusters and to record their infrared spectra via messenger-tagged infrared multiple photon dissociation spectroscopy. With complementary quantum chemical calculations and molecular dynamics simulations, we determined that the most probable binding site of a metal, specifically a vanadium cation, on C60 is above a pentagonal center in an η5 fashion. We explored the bonding nature between C60 and V+ and revealed that the high thermal stability of this cluster originates from large orbital and electrostatic interactions. Through comparing the measured infrared spectra of [C60-Metal]+ with the observational Spitzer data of several fullerene-rich planetary nebulae, we proposed that the complexes formed by fullerene and cosmically abundant metals, for example, iron, are promising carriers of astronomical unidentified spectroscopic features. This opens the door for a real consideration of Kroto's 30-year-old hypothesis that complexes involving cosmically abundant elements and C60 exhibit strong charge-transfer bands, similar to those of certain unidentified astrophysical spectroscopic features. We compiled a VibFullerene database and extracted a set of vibrational frequencies and intensities for fullerene derivatives to facilitate their potential detection by the James Webb Space Telescope. In addition, we showed that upon infrared irradiation C60V+ can efficiently catalyze water splitting to generate H2. This finding is attributed to the novel geometric-electronic effects of C60, acting as "hydrogen shuttle" and "electron sponge", which illustrates the important role of carbon-based supports in single-atom catalysts. Our work not only unveils the basic structures and bonding nature of fullerene-metal clusters but also elucidates their potential importance in astrophysics, astrochemistry, and catalysis, showing the multifaceted character of this class of clusters. More exciting and interesting aspects of the fullerene-metal clusters, such as ultrafast charge-transfer dynamics between fullerene and metal and their relevance to designing hybrid fullerene-metal junctions for electronic devices, are awaiting exploration.

Pentagon, Hexagon, or Bridge? Identifying the Location of a Single Vanadium Cation on Buckminsterfullerene Surface.

Buckminsterfullerene C60 has received extensive research interest since its discovery. In addition to its interesting intrinsic properties of exceptional stability and electron-accepting ability, the broad chemical tunability by decoration or substitution on the C60-fullerene surface makes it a fascinating molecule. However, to date, there is uncertainty about the binding location of such decorations on the C60 surface, even for a single adsorbed metal atom. In this work, we report the gas-phase synthesis of the C60V+ complex and its in situ characterization by mass spectrometry and infrared spectroscopy with the help of quantum chemical calculations and molecular dynamics simulations. We identify the most probable binding position of a vanadium cation on C60 above a pentagon center in an η5-fashion, demonstrate a high thermal stability for this complex, and explore the bonding nature between C60 and the vanadium cation, revealing that large orbital and electrostatic interactions lie at the origin of the stability of the η5-C60V+ complex.

Effect of Thermal Treatment on the Phase Composition and Surface Properties of WO3‐TiO2 Nanocomposites Synthesized via Hydro‐Thermal Method

AbstractThere is a huge interest in the synthesis of nanocomposites with specific physicochemical properties. Calcination temperature affects the size, morphology and other surface characteristics. WO3‐TiO2 nanocomposites were synthesized by hydrothermal method using titanium (IV) butoxide and tungsten (VI) oxide (WO3). The effects of drying time and calcination temperature on the hydrodynamic particle diameter, crystallite size, shape and zeta potential of the nanocomposites was investigated. The samples were dried in oven at 80 °C for 18 hours, calcinated at 300 to 1000 °C and analyzed by X‐ray diffraction (XRD), field emission scanning electron microscope (FE‐SEM) and fourier‐transform infrared spectroscopy (FTIR). Hydrodynamic diameter and zeta potential were measured by particle analyzer (DLS). The mean crystallite size increased from 25.2 nm to 54.9 nm when calcination temperature increased from 300 to 1000 °C. Hydrodynamic size of the nanocomposite decreased with increasing drying time, became constant after 12 hours of drying and did not change significantly after that. Moreover, the hydrodynamic size increased with increase in the calcination temperature due to growth of the crystals and formation of inter particle necks. The large zeta values (−35.8 to −41.2 mV) suggest that the nanocomposites show large electrostatic repulsive interaction and enhanced stability in aqueous solutions.

Pulling a DNA molecule through a nanopore embedded in an anionic membrane: tension propagation coupled to electrostatics

We consider the influence of electrostatic forces on driven translocation dynamics of a flexible polyelectrolyte being pulled through a nanopore by an external force on the head monomer. To this end, we augment the iso-flux tension propagation theory with electrostatics for a negatively charged biopolymer pulled through a nanopore embedded in a similarly charged anionic membrane. We show that in the realistic case of a single-stranded DNA molecule, dilute salt conditions characterized by weak charge screening, and a negatively charged membrane, the translocation dynamics is unexpectedly accelerated despite the presence of large repulsive electrostatic interactions between the polymer coil on the cis side and the charged membrane. This is due to the rapid release of the electrostatic potential energy of the coil during translocation, leading to an effectively attractive force that assists end-driven translocation. The speedup results in non-monotonic polymer length and membrane charge dependence of the exponent α characterizing the translocation time of the polymer with length N0. In the regime of long polymers N0 ≳ 500, the translocation exponent exceeds its upper limit α = 2 previously observed for the same system without electrostatic interactions.

Densely Grafted Polyelectrolyte Brushes Trigger “Water-in-Salt”-like Scenarios and Ultraconfinement Effect

Summary This paper presents an all-atom molecular dynamics (MD) simulation of highly charged and densely grafted polyelectrolyte (PE) brushes on a flat surface. Simulation findings for the variation of brush height with polymer size and grafting density are explained by the scaling laws for brushes in the non-linear osmotic regime. More importantly, this study establishes the triggering of an ultraconfinement effect by the densely grafted brushes. This effect leads to significant changes in the distribution, structure, and properties as well as a massive mobility reduction of both the counterions and water. Furthermore, the interplay of the ultraconfinement effect and the large counterion-PE electrostatic interactions trigger a “water-in-salt”-like scenario (witnessed in highly concentrated aqueous electrolyte solutions), which is characterized by the counterion-PE-functional group serving as the “salt” with the “salt” overwhelming the water in both mass and volume and affecting the solvation structure of the counterions by replacing the solvation water with the PE-functional group.

Salting Up of Proteins at the Air/Water Interface.

Vibrational sum frequency generation (VSFG) spectroscopy and surface pressure measurements are used to investigate the adsorption of a globular protein, bovine serum albumin (BSA), at the air/water interface with and without the presence of salts. We find at low (2 to 5 ppm) protein concentrations, which is relevant to environmental conditions, both VSFG and surface pressure measurements of BSA behave drastically different from at higher concentrations. Instead of emerging to the surface immediately, as observed at 1000 ppm, protein adsorption kinetics is on the order of tens of minutes at lower concentrations. Most importantly, salts strongly enhance the presence of BSA at the interface. This "salting up" effect differs from the well-known "salting out" effect as it occurs at protein concentrations well-below where "salting out" occurs. The dependence on salt concentration suggests this effect relates to a large extent electrostatic interactions and volume exclusion. Additionally, results from other proteins and the pH dependence of the kinetics indicate that salting up depends on the flexibility of proteins. This initial report demonstrates "salting up" as a new type of salt-driven interfacial phenomenon, which is worthy of continued investigation given the importance of salts in biological and environmental aqueous systems.

Effects of Electronic Structure of Adjacent Carbon on the Strength of C─F⋯H─F Organofluorine Hydrogen Bonds.

We investigate the effects of the electronic structure of carbon atom on the organofluorine hydrogen bonds, C─F⋯H─F. Our results show that we can modulate the strength of organofluorine hydrogen bonds by adjusting the volume of fluorine atom in C─F via changing the electronic structure of adjacent carbon atoms. Different with the conventional hydrogen bonds, we found that instead of carbon rehybridization and hyperconjugative effects, the magnitude of fluorine atomic volume plays important roles in determining the strength of the C─F⋯H─F organofluorine hydrogen bonds. The lone pair electrons at both the proximal and the vicinal carbon dramatically reinforce the strength of C─F⋯H─F organofluorine hydrogen bond with its interaction energy in the range of about 15-25 kcal/mol, that is, the carbanion-mediated organofluorine hydrogen bond could be very strong. Due to the high electronegativity of fluorine atom, it easily attracts the excess electron from the proximal and vicinal carbon, which results in the increase of its volume and negative charge. The enhanced volume of fluorine atom gives rise to the large polarization energy, and its enhanced negative charge favors the large electrostatic interaction, both of which substantially contribute to making the organofluorine hydrogen bonds strong. © 2019 Wiley Periodicals, Inc.

Structural basis of fullerene derivatives as novel potent inhibitors of protein acetylcholinesterase without catalytic active site interaction: insight into the inhibitory mechanism through molecular modeling studies

Acetylcholinesterase (AChE) is an important kind of esterase that plays a key biological role in the central and peripheral nervous systems. Recent research has demonstrated that some fullerene derivatives serve as a new nanoscale class of potent inhibitors of AChE, but the specific inhibition mechanism remains unclear. In the present article, several molecular modeling methods, such as molecular docking, molecular dynamics simulations and molecular mechanics/generalized Born surface area calculations, were used for the investigation of the binding mode and inhibition mechanism of fullerene inhibitions for AChE. Results revealed that fullerene inhibitors block the entrance of substrates by binding with the peripheral anionic site (PAS) region. Thus, fullerene derivatives might mainly serve as competitive inhibitors. The interactions of a fullerene backbone with AChE are at the same level in different single side chain systems and seem to be related to the length or aromaticity of the side chain. The inhibitor with multihydroxyl side chains shows an obviously large electrostatic interaction as it forms additional hydrogen bonds with AChE. Moreover, fullerene derivatives might exhibit noncompetitive inhibition partly by affecting some secondary structures around them. Thus, the destructions of these secondary structures can lead to conformational changes in some important regions, such as the catalytic triad and acyl pocket. The investigation is of great importance to the discovery of good fullerene inhibitors.Communicated by Ramaswamy H. Sarma

Understanding and Quantifying London Dispersion Effects in Organometallic Complexes.

Quantum chemical methods are nowadays able to determine properties of larger chemical systems with high accuracy and Kohn-Sham density functional theory (DFT) in particular has proven to be robust and suitable for everyday applications of electronic structure theory. A clear disadvantage of many established standard density functional approximations like B3LYP is their inability to describe long-range electron correlation effects. The inclusion of such effects, also termed London dispersion, into DFT has been extensively researched in recent years, resulting in some efficient and routinely used correction schemes. The well-established D3 method has demonstrated its efficiency and accuracy in numerous applications since 2010. Recently, it was improved by developing the successor (termed D4) which additionally includes atomic partial charge information for the generation of pairwise dispersion coefficients. These coefficients determine the leading-order (two-body) and higher-order (three- or many-body) terms of the D4 dispersion energy which is simply added to a standard DFT energy. With its excellent accuracy-to-cost ratio, the DFT-D4 method is well suited for the determination of structures and chemical properties for molecules of most kinds. While dispersion effects in organic molecules are nowadays well studied, much less is known for organometallic complexes. For such systems, there has been a growing interest in designing dispersion-controlled reactions especially in the field of homogeneous catalysis. Here, efficient electronic structure methods are necessary for screening of promising model complexes and quantifying dispersion effects. In this Account, we describe the quality of calculated structural and thermodynamic properties in gas-phase obtained with DFT-D4 corrected methods, specifically for organometallic complexes. The physical effects leading to London dispersion interactions are briefly discussed in the picture of second-order perturbation theory. Subsequently, basic theoretical aspects of the D4 method are introduced followed by selected case studies. Several chemical examples are presented starting with the analysis of transition metal thermochemistry and noncovalent interactions for small, heavy element containing main group compounds. Computed reaction energies can only match highly accurate reference values when all energy contributions are included in the DFT treatment, thus highlighting the major role of dispersion interactions for the accurate description of thermochemistry in gas-phase. Furthermore, the correlation between structural and catalytic properties is emphasized where the accessibility of high quality structures is essential for reaction planning and catalyst design. We present calculations for aggregates of organometallic systems with intrinsically large repulsive electrostatic interactions which can be stabilized by London dispersion effects. The newly introduced inclusion of atomic charge information in the DFT-D4 model robustly leads to quantitatively improved dispersion energies in particular for metallic systems. By construction it yields results which are easily understandable due to a clear separation into hybridization and charge (oxidation) state and two- and many-body effects, respectively. Due to its high computational efficiency, the D4 dispersion model is even applicable to low-cost classical and semiempirical theoretical methods.

Doping of alkali metals in carbon frameworks for enhancing CO2 capture: A theoretical study

We use a multi-scale simulation method that combines the grand canonical Monte Carlo (GCMC) simulations and density function theory (DFT), a series of alkali metals doping in carbon surfaces, and the influence of alkali metal-doped PCs on CO2 uptake were comprehensive study. The results show that all alkali metals can significantly increase the CO2 uptake in carbon surfaces. The enhancement of CO2 adsorption mainly originates from the effects of the large electrostatic interaction and the strong adsorption energy. By the comparative studies above, it is found that the doping of alkali metals into carbon surface for CO2 capture reach 8.43–12.46 mmol g−1, which is about 2.6–3.8 times those in the nondoped one. In combination with its superior selectivity for CO2/N2 and CO2/CH4 at 25 °C, the doping of alkali metals on carbon surface has been shown to be one of the promising adsorbents for CO2 uptake and separation, and provides an effective and excellent alternative strategy for the design and screening of CO2 adsorbents.

Fragment Molecular Orbital Calculations with Implicit Solvent Based on the Poisson-Boltzmann Equation: Implementation and DNA Study.

In this study, an ab initio fragment molecular orbital (FMO) methodology was developed to evaluate the solvent effects on electrostatic interactions, which make a significant contribution to the physical and chemical processes occurring in biological systems. Here, a fully polarizable solute consisting of the FMO electron density was electrostatically coupled with an implicit solvent based on the Poisson-Boltzmann (PB) equation; in addition, the nonpolar contributions empirically obtained from the molecular surface area (SA) were added. Interaction analysis considering solvent-screening and dispersion effects is now available as a powerful tool to determine the local stabilities inside solvated biomolecules. This methodology is applied to a deoxyribonucleic acid (DNA) duplex known as the Dickerson dodecamer. We found that excessively large electrostatic interactions inside the duplex are effectively damped by the screening, and the frontier molecular orbital energies are also successfully lowered. These observations indicate the stability of highly charged DNA duplexes in solution. Moreover, the solvation free energies in the implicit model show fairly good agreement with those in the explicit model while avoiding the costly statistical sampling of the electrolyte distribution. Consequently, our FMO-PBSA approach could yield new insights into biological phenomena and pharmacological problems via this ab initio methodology.

Exploring surface landscapes with molecules: rotationally induced diffraction of H2 on LiF(001) under fast grazing incidence conditions.

Atomic diffraction by surfaces under fast grazing incidence conditions has been used for almost a decade to characterize surface properties with more accuracy than with more traditional atomic diffraction methods. From six-dimensional solutions of the time-dependent Schrödinger equation, we show that diffraction of H2 molecules under fast grazing incidence conditions could be even more informative for the characterization of ionic surfaces, due to the large anisotropic electrostatic interaction between the quadrupole moment of the molecule and the electric field created by the ionic crystal. Using the LiF(001) surface as a benchmark, we show that fast grazing incidence diffraction of H2 strongly depends on the initial rotational state of the molecule, while rotationally inelastic processes are irrelevant. We demonstrate that, as a result of the anisotropy of the impinging projectile, initial rotational excitation leads to an increase in intensity of high-order diffraction peaks at incidence directions that satisfy precise symmetry constraints, thus providing a more detailed information on the surface characteristics than that obtained from low-order atomic diffraction peaks under fast grazing incidence conditions. As quadrupole-ion surface potentials are expected to accurately represent the interaction between H2 and any surface with a marked ionic character, our analysis should be of general applicability to any of such surfaces. Finally, we show that a density functional theory description of the molecule-ion surface potential catches the main features observed experimentally.

Interaction between ions and substituted buckybowls: a comprehensive computational study.

Complexes formed by substituted buckybowls derived from corannulene and sumanene with sodium cation or chloride anion have been computationally studied by using a variety of methods. Best results have been obtained with the SCS-MP2 method extrapolated to basis set limit, which reproduces the highest-level values obtained with the MP2.X method. All bowls form stable complexes with chloride anion, with stabilities ranging from -6 kcal/mol in the methylated corannulene derivative to -45 kcal/mol in the CN-substituted sumanene. The opposite trend is observed in sodium complexes, going from deeply attractive complexes with the methylated derivatives (-36 kcal/mol with sumanene derivative) to slightly repulsive ones in the CN-substituted bowls (2 kcal/mol in the corannulene derivative). Anion complexes are stabilized by large electrostatic interactions combined with smaller though significant dispersion and induction contributions. Conversely, cation complexes are stabilized by large induction contributions capable of holding together the bowl and the cation even in cases where the electrostatic interaction is repulsive. The effect of substitution is mainly reflected on changes in the molecular electrostatic potential of the bowl and, thus, in the electrostatic contribution to the interaction. Therefore, the variations in the stability of the complexes on substitution could be roughly predicted just considering the changes in the electrostatic interaction. However, other contributions also register changes mainly as a consequence of displacements on the position of the ion at the minimum, so the accurate prediction of the stability of this kind of complexes requires going further than the electrostatic approach.

Self-consistent continuum solvation (SCCS): The case of charged systems

The recently developed self-consistent continuum solvation model (SCCS) [O. Andreussi, I. Dabo, and N. Marzari, J. Chem. Phys. 136, 064102 (2012)] is applied here to charged species in aqueous solutions. Describing ions in solution represents a great challenge because of the large electrostatic interactions between the solute and the solvent. The SCCS model is tested over 106 monocharged species, both cations and anions, and we demonstrate its flexibility, notwithstanding its much reduced set of parameters, to describe charged species in solution. Remarkably low mean absolute errors are obtained with values of 2.27 and 5.54 kcal/mol for cations and anions, respectively. These results are comparable or better than the state of the art to describe solvation of charged species in water. Finally, differences of behavior between cations and anions are discussed.

Origin of Attraction in Chalgogen–Nitrogen Interaction of 1,2,5-Chalcogenadiazole Dimers

Intermolecular interaction in the 1,2,5-chalcogenadiazole dimers was studied by ab initio molecular orbital calculations. Estimated CCSD(T) interaction energies for the thia-, selena- and tellura-diazole dimers are -3.14, -5.29, and -12.42 kcal/mol, respectively. The electrostatic and dispersion interactions are the major sources of the attraction in the dimers, although it was claimed that the orbital mixing (charge-transfer interaction) was the most prominent contribution to the stabilization. The induction (induced polarization) interaction also contributes largely to the attraction in the telluradiazole dimer. The large electrostatic and induction interactions are responsible for the strong attraction in the telluradiazole dimer. The short-range (orbital-orbital) interaction (sum of the exchange-repulsion and charge-transfer interactions) is repulsive. The directionality of the interactions increases in order of S < Se < Te. The electrostatic interaction is mainly responsible for the directionality. The strong directionality suggests that the chalcogen-nitrogen interaction plays important roles in controlling the orientation of molecules in those organic crystals. The nature of the chalcogen-nitrogen interaction in the chalcogenadiazole dimers is similar to that of the halogen bond, which is an electrostatically driven noncovalent interaction.

Systematic Trends in Non-resonant Inelastic X-ray Scattering of Actinides

ABSTRACTNonresonant inelastic x-ray scattering (NIXS) measurements probing the 5d → 5f electronic transitions at the uranium O4,5 edges in a localized electron system such as UO2 show a good agreement with many-electron atomic spectral calculations. The higher multipole spectra are split into two peaks that can be assigned as the 5d5/2 and 5d3/2 structures, despite the large electrostatic interactions. A new sum rule for the branching ratio of spin-orbit split core levels in electric multipole spectra, which generalizes the sum rule for dipole transitions in x-ray absorption spectroscopy, allows for a systematic analysis of the NIXS spectra. The branching ratio is linearly proportional to the expectation value of the angular part of the spin-orbit interaction in the initial state, where the rank of the multipole determines the coefficient of the linear dependence. This spin-orbit sum rule can be an important diagnostic tool for high-energy spectroscopies.

YVO4:Eu3+functionalized porous silica submicrospheres as delivery carriers of doxorubicin

Porous silica microspheres were fabricated by a facile surface-protected etching strategy. Polyvinylpyrrolidone (PVP) was used as a protecting polymer absorbed on the surface of silica microspheres and NaOH was employed as an etching agent. Owing to the protective action of PVP and inhomogeneous etching, mesopores were created in the silica microspheres. Then, based on the Pechini-type sol-gel and impregnating process, YVO(4):Eu(3+) nanocrystals were integrated into the channels to form highly luminescent YVO(4):Eu(3+)@SiO(2) composite microspheres. The biocompatibility tests on L929 fibroblast cells using MTT assay reveal low cytotoxicity of the system. Owing to the large interior space and electrostatic interaction, the porous microspheres show a relatively high loading capacity (438 mg DOX/YVO(4):Eu(3+)@SiO(2) g) and encapsulation efficiency (87.6%) for the anti-cancer drug doxorubicin hydrochloride (DOX). The drug release behavior and cytotoxic effect against human cervical carcinoma cells (HeLa cells) of the DOX-loaded YVO(4):Eu(3+)@SiO(2) carriers were investigated in vitro. It was found that the carriers present a highly pH-dependent drug release behavior due to electrostatic interaction between the silica surface and DOX molecules. The drug release rate became greater at low pH owing to the increased electrostatic repulsion. The DOX-loaded carriers demonstrate a similar or even greater anti-cancer activity with respect to the free DOX against HeLa cells. Furthermore, the PL intensity of the microspheres shows correlation with the cumulative release of DOX. These results suggest that the composite can potentially act as a multifunctional drug carrier system with luminescent tagging and pH-controlled release properties.

Dipole-Induced, Thermally Stable Lamellar Structure by Polar Aromatic Silane

Controlled self-assembly of polar aromatic silane leads to the formation of well-ordered lamellar structures. Graphite-like features are clearly visible with a scanning electron microscope (SEM). In addition, X-ray diffraction (XRD) patterns suggest a d spacing of 14.28 A along the z-axis and 4.42 A in the xy plane, which all agree with theoretical modeling. Constructing multistacks of silane molecules with a high degree of ordering is a daunting task. Amorphous monolayers are frequently reported. Aggravated van der Waals interaction, pi-pi electron overlapping, and solvophobic interactions can all lead to the formation of multistacks. The importance of a dipole to the ordered stacking is essentially unknown. This work suggests that a strong dipole-dipole interaction can be another important driving force in forming lamellar structures. The resulting large electrostatic interactions between the dipole and water provide an excellent thermal stability for these lamellas up to 350 degrees C. Organized, layered structures with a permanent dipole can be used in piezoelectric devices or as active surfaces to bind polar molecules, such as toxic gas, methanol, or DNA.

Structure and energy decomposition analysis of metal–metal bonding in [PhM–MPh] and [ClM–MCl] (M = Zn, Cd, Hg)

The metal–metal bonding in the title compounds have been investigated with the help of energy decomposition analysis at the DFT/TZ2P level. The optimized Zn–Zn and Cd–Cd bond distances in model compounds [Ph 2M 2] (M = Zn, Cd) are slightly longer than the experimental values in [ R 2 M 2 ] [ R = C 6 H 3 - 2 , 6 - ( C 6 H 3 - 2 , 6 - Pr 2 i ) 2 ] . The calculated data show that bond dissociation energies (−BDE) for Zn–Zn bonds −53.99 kcal/mol in [Ph 2Zn 2] and −59.41 kcal/mol in [Cl 2Zn 2] are greater amongst the compounds under study. In addition, [Cl 2M 2] is found to have a bonding energy slightly larger than those in [Ph 2M 2]. The values of interaction energy, Δ E int, electrostatic interactions, Δ E elstat, and orbital interactions, Δ E orb are arranged in the following order: Zn > Hg > Cd. There are strong electrostatic attractions between metals, with Δ E elstat (−69.23 to −103.55 kcal/mol). The electrostatic attractions, Δ E elstat are always greater than the orbital interactions, Δ E orb. The M–M bonding has more than half ionic character (53–62%). The reason for the large electrostatic interactions is the anisotropic charge distribution at the metal atom. The results demonstrate clearly that the atomic partial charges cannot be taken as a measure of the electrostatic interactions between the atoms. The orbital interactions are relatively weaker in the cadmium–cadmium bond than in the mercury–mercury bond, whereas the Cd–Cd bond distances are longer than the Hg–Hg bond distances. We want to point out that the values of Δ E orb is not simply a function of the interatomic distances.

Model transition states for methane diazonium ion methylation of guanine runs in oligomeric DNA

The DNA reaction pattern of the methane diazonium ion, which is the reactive intermediate formed from several carcinogenic methylating agents, was examined at N7 and O(6) sites in guanine runs occurring in oligonucleotides and model oligonucleotides. Density functional B3LYP/6-31G*, and SCF 3-21G and STO-3G energies of model transition states were calculated in the gas phase and in the CPCM reaction field. For nucleotides containing two, three, and four stacked guanines with counterions in the gas phase, O(6) reactivity is greater than N7 reactivity. In the reaction field, N7 reactivity is 9.0 to 9.8 times greater than O(6) reactivity. For a double-stranded oligonucleotide containing two stacked guanines with counterions in the reaction field, the N7 and O(6) reactivities of the 3'-guanine are 3.9 times greater than the corresponding sites in the 5'-guanine. For double-stranded oligonucleotides with three or four stacked guanines and counterions, the reactivities of the interior guanines are higher than corresponding reactivities of guanines at the ends. These reaction patterns agree with most of the available experimental data. Activation energy decomposition analysis for gas-phase reactions in a double-stranded dinucleotide containing two stacked guanines with counterions indicates that selectivity at O(6) is almost entirely due to electrostatic forces. Selectivity at N7 also has a large electrostatic interaction. However, the orbital interaction also contributes significantly to the gas-phase selectivity, accounting for 32% of the total interaction energy difference between the 3'- and 5'-guanine reactions. In aqueous solution, the relative orbital contribution to N7 selectivity is likely to be larger.