Bond critical points in the electronic structures of the main group diatomic hydrides of lithium through bromine

The bond critical points and associated electronic properties of the diatomic hydrides of the twenty-one main group elements from lithium to bromine have been calculated with large basis sets. As part of a systematic study of the polarity of chemical bonds, the position of the bond critical point, the charge density at the bond critical point, the Laplacian of the charge density at the bond critical point, and the molecular dipole moment of each molecule have been calculated. Particular attention has been paid to the effect of bond length elongation and contraction on the electronic properties. Variation of the bond length reveals that with atoms of low electronegativity, the bond critical point of AH tends to follow atom A, whereas with atoms of high electronegativity, the bond critical point tends to follow the hydrogen atom as the bond lengthens. Furthermore, it is shown that some properties of the diatomic hydrides vary monotonically within each row of the periodic table, while others effect a classification according to the character of the bond.     

Topology of charge density of flucytosine and related molecules and characteristics of their bond charge distributions

The molecular charge distribution of flucytosine (4-amino-5-fluoro-2-pyrimidone), uracil, 5-fluorouracil, and thymine was studied by means of density functional theory calculations (DFT). The resulting distributions were analyzed by means of the atoms in molecules (AIM) theory. Bonds were characterized through vectors formed with the charge density value, its Laplacian, and the bond ellipticity calculated at the bond critical point (BCP). Within each set of C=O, C-H, and N-H bonds, these vectors showed little dispersion. C-C bonds formed three different subsets, one with a significant degree of double bonding, a second corresponding to single bonds with a finite ellipticity produced by hyperconjugation, and a third one formed by a pure single bond. In N-C bonds, a decrease in bond length (an increase in double bond character) was not reflected as an increase in their ellipticity, as in all C-C bonds studied. It was also found that substitution influenced the N-C, C-O, and C-C bond ellipticity much more than density and its Laplacian at the BCP. The Laplacian of charge density pointed to the existence of both bonding and nonbonding maxima in the valence shell charge concentration of N, O, and F, while only bonding ones were found for the C atoms. The nonbonding maxima related to the sites for electrophilic attack and H bonding in O and N, while sites of nucleophilic attack were suggested by the holes in the valence shell of the C atoms of the carbonyl groups.     

Progress in the understanding of drug-receptor interactions, Part 1: experimental charge-density study of an angiotensin II receptor antagonist (C30H30N6O3S) at T = 17 K

An experimental study of the electron-density distribution rho(r) in an angiotensin II receptor antagonist 1 has been made on the basis of single-crystal X-ray diffraction data collected at a low temperature. The crystal structure of 1 consists of infinite ribbons in which molecules are connected by an N-H...N hydrogen bond and several interactions of the C-H...O, C-H...N, and C-H...S type. The molecular conformation, characterized by the syn orientation of a tetrazole and a pyrimidinone ring with respect to a phenyl spacer group, is stabilized by two short SO and SN intramolecular contacts between a substituted thiophene fragment and the other two heterocycles of 1. The electrostatic nature of these interactions is documented. Furthermore, the Laplacian of rho(r) in the plane defined by the sulfur, oxygen, and nitrogen atoms involved in these interactions shows their strongly directional character as the regions of charge concentration on the valence shell of the nitrogen and oxygen atoms directly face the regions of charge depletion on the valence shell of the sulfur atom. All the chemical bonds and the relevant intra- and intermolecular interactions of 1 have been quantitatively described by the topological analysis of rho(r). Simple relationships between the bond path lengths (R(b)) and the values of rho at the bond critical points (rho(bcp)) have been obtained for the 28 C-C bonds, the seven N-C bonds, and the four O-C bonds. For the first two classes of bonds the relationship is in the form of a straight line, whose parameters, for the C-C bonds, agree, within experimental uncertainty, with those previously derived in our laboratory from a 19 K X-ray diffraction study of crystals of a different compound. Maps of the molecular electrostatic potential phi(r) derived from the experimental charge density display features that are important for the drug-receptor recognition of 1.     

THEORETICAL STUDIES ON THE SUBSTITUTION EFFECT OF FLUORINE ON ETHANE

<正> The substitution effect of fluorine on ethane has been investigated by means of studying the properties of the charge distribution at the bond critical points with the theory of atoms in molecule.It is found that the major substitution effects of fluorine atom are positive a inductive and polarity effect.At the same time,fluorine atom partially provides π electrons to other chemical bonds by means of hy-perconjugation in molecules with two fluorine atoms and one or two carbon atoms in the same plane,and these effects are reflected in the quantity of bond ellipticity,Laplacian and the charge density of charge distribution at the bond critical points.The substitution of hydrogen by fluorine in ethane strengthens all the bonds in substituted ethanes.Other effects originating from the substitution of hydrogen by fluorine have also been discussed.     

Topological analysis of the electronic charge density in the ethene protonation reaction catalyzed by acidic zeolite

In the present work, the distribution of the electronic charge density in the ethene protonation reaction by a zeolite acid site is studied within the framework of the density functional theory and the atoms in molecules (AIM) theory. The key electronic effects such as topological distribution of the charge density involved in the reaction are presented and discussed. The results are obtained at B3LYP/6-31G(**) level theory. Attention is focused on topological parameters such as electron density, its Laplacian, kinetic energy density, potential energy density, and electronic energy density at the bond critical points (BCP) in all bonds involved in the interaction zone, in the reactants, pi-complex, transition state, and alkoxy product. In addition, the topological atomic properties are determined on the selected atoms in the course of the reaction (average electron population, N(Omega), atomic net charge, q(Omega), atomic energy, E(Omega), atomic volume, v(Omega), and first moment of the atomic charge distribution, M(Omega)) and their changes are analyzed exhaustively. The topological study clearly shows that the ethene interaction with the acid site of the zeolite cluster, T5-OH, in the ethene adsorbed, is dominated by a strong O-H...pi interaction with some degree of covalence. AIM analysis based on DFT calculation for the transition state (TS) shows that the hydrogen atom from the acid site in the zeolitic fragment is connected to the carbon atom by a covalent bond with some contribution of electrostatic interaction and to the oxygen atom by closed shell interaction with some contribution of covalent character. The C-O bond formed in the alkoxy product can be defined as a weaker shared interaction. Our results show that in the transition state, the dominant interactions are partially electrostatic and partially covalent in nature, in which the covalent contribution increases as the concentration and accumulation of the charge density along the bond path between the nuclei linked increases.     

DFT-B3LYP, NPA-, and QTAIM-based study of the physical properties of [M(II)(H2O)2(15-crown-5)] (M = Mn, Fe, Co, Ni, Cu, Zn) complexes

A density functional theory study of the structure of the title compounds with the divalent metal ions in their high-spin ground state, obtained using B3LYP/6-311++G(d,p) in vacuo and in aqueous solution simulated using a polarized continuum medium, is reported for the first time. The modeling reproduces the pseudo pentagonal bipyramidal crystallographic structures very well, including some asymmetry in the equatorial bonds lengths to the crown ether O donors. The very marked asymmetry in the Ni(2+) structure due to a Jahn-Teller distortion of a d(8) system in a D(5h) ligand field is also well reproduced. The gas phase binding energies of the complexes follow the order Mn(2+) < Fe(2+) < Co(2+) < Ni(2+) < Cu(2+) > Zn(2+), in precise agreement with the Irving-William series. Both the NPA and Bader charges show there is ligand-to-metal charge transfer; however, the values obtained from the NPA procedure, unlike those obtained from Bader's quantum theory of molecules approach, do not correlate with the electronegativity of the metal ions, the stabilization energies of the solvated complexes or the ionic radii of the metal ions, and so appear to be less reliable. The nature of the bonding between the ligands and the metal ions has been explored using the topological properties of the electron charge density. The metal-ligand bond distances were found to be exponentially correlated with the electron charge density, its Laplacian, and with its curvature in the direction of the bond path at M-O bond critical points. While the bonding with coordinated H(2)O is predominantly ionic, that to the crown ether donor atoms has some covalent character the extent of which increases across the first transition series. The delocalization indices of M-O bonds in these complexes correlate reasonably well with the electron density and its Laplacian at the bond critical points; this therefore provides a rapid and computationally very efficient way of determining these properties, from which insight into the nature of the bonding can be obtained, obviating the need for time-consuming integration over atomic basins.     

Strength and nature of hydrogen bonding interactions in mono- and di-hydrated formamide complexes

In this work, mono- and di-hydrated complexes of the formamide were studied. The calculations were performed at the MP2/6-311++G(d,p) level of approximation. The atoms in molecules theory (AIM), based on the topological properties of the electronic density distribution, was used to characterize the different types of bonds. The analysis of the hydrogen bonds (H-bonds) in the most stable mono- and di-hydrated formamide complexes shows a mutual reinforcement of the interactions, and some of these complexes can be considered as "bifunctional hydrogen bonding hydration complexes". In addition, we analyzed how the strength and the nature of the interactions, in mono-hydrated complexes, are modified by the presence of a second water molecule in di-hydrated formamide complexes. Structural changes, cooperativity, and electron density redistributions demonstrate that the H-bonds are stronger in the di-hydrated complexes than in the corresponding mono-hydrated complexes, wherein the σ- and π-electron delocalization were found. To explain the nature of such interactions, we carried out the atoms in molecules theory in conjunction with reduced variational space self-consistent field (RVS) decomposition analysis. On the basis of the local Virial theorem, the characteristics of the local electron energy density components at the bond critical points (BCPs) (the 1/4? (2)ρ(b) component of electron energy density and the kinetic energy density) were analyzed. These parameters were used in conjunction with the electron density and the Laplacian of the electron density to analyze the characteristics of the interactions. The analysis of the interaction energy components for the systems considered indicates that the strengthening of the hydrogen bonds is manifested by an increased contribution of the electrostatic energy component represented by the kinetic energy density at the BCP.     

Hydrogen bonding in substituted formic acid dimers

The hydrogen-bonded dimers of formic acid derivatives XCOOH (X = H, F, Cl, and CH3) have been investigated using density functional theory (B3LYP) and second-order M?ller-Plesset perturbation (MP2) methods, with the geometry optimization carried out using 6-311++G(2d,2p) basis set. The dimerization energies calculated using aug-cc-pVXZ (with X = D and T) basis have been extrapolated to infinite basis set limit using the standard methodology. The results indicate that the fluorine-substituted formic acid dimer is the most stable one in comparison to the others. Topological analysis carried out using Bader's atoms in molecules (AIM) theory shows good correlation of the values of electron density and its Laplacian at the bond critical points (BCP) with the hydrogen bond length in the dimers. Natural bond orbital (NBO) analysis carried out to study the charge transfer from the proton acceptor to the antibonding orbital of the X-H bond in the complexes reveals that most of the dimers are associated with conventional H-bonding except a few, where improper blue-shifting hydrogen bonds are found to be present.     

Electron density topological analysis of hydrogen bonding in glucopyranose and hydrated glucopyranose

Topological analysis of the electron density profiles and the atomic basin integration data for the most energetically favorable (4)C(1) and (1)C(4) conformers of beta-D-glucopyranose, calculated at the B3LYP/6-31+G(d), MPWlPW91/6-311+G(2d,p), and MP2/6-31+G(d) levels, demonstrates that intramolecular hydrogen bonding between adjacent ring OH groups does not occur in glucopyranose, given the need to demonstrate a bond critical point (BCP) of correct (3,-1) topology for such an interaction to be termed a hydrogen bond. On the other hand, pyranose ring OH groups separated by three, rather than just two, carbon atoms are able to form an intramolecular hydrogen bond similar in topological properties and geometry to that found for propane-1,3-diol. Vicinal, equatorial OH groups in the (4)C(1) conformer of glucopyranose are, however, able to form strong bidentate hydrogen bonds with water molecules in a cooperative manner, each water molecule acting simultaneously as both hydrogen bond donor and acceptor, and characterized by (3,-1) bond critical points with increased values for the electron density and the Laplacian of rho(r) compared to an isolated ethane-1,2-diol/water complex.     

Topological properties of the peptide bond in glycyl-L-threonine dihydrate based on a fast synchrotron/CCD-diffraction experiment at 100 K

The charge density of glycyl-L-threonine dihydrate is extracted from a synchrotron data set of 98405 reflections collected at 100 K with a Bruker CCD area detector up to a resolution of d=0.38 A (sintheta/lambda = 1.32 A 1). The data are interpreted in terms of the "rigid pseudoatom" model. The topology of the experimental density is analyzed and compared with the topology obtained experimentally for the constituting amino acids and to that derived from Hartree-Fock calculations for the isolated molecule. All critical points of the electron density at the covalent and hydrogen bonds, as well as those of the Laplacian, were located, thereby deriving quantitative topological data for the peptide and side chain bonds. Bond topological indices in the dipeptide compare well with those of the corresponding bonds in the building amino acids, thus suggesting transferability of electronic properties of atoms and functional groups when these are derived by Bader's partitioning. Discrepancies between theoretical and experimental results could be attributed to crystal field effects.     

羟基二酰亚胺二聚体分子间相互作用的量子化学研究

采用从头算MP2方法在6-311++G~**基组水平上讨论了CP梯度校正对两种羟基 二酰亚胺异构体所的相互作用能和几何结构的影响，并利用分子中的原子理论（ Atoms in molecules，AIM）计算了五个拓扑参数：键临荷密度、电荷密度的 Iaplacian值、氢键中氢原子的体积、氢原子集居数、氢原子能量来表征氢键的形 成．种构型氢键体系中还讨论了二聚体的相互作用能与氢键临界点的电荷密度、质 子供体X-H键长的线性相关性问题．表明这种线性相关性的存在有范围限制，复合 物和其中单体的构型能够影响这种关系的存在.     

Single‐electron halogen bond: Ab initio study

Ab initio calculations have been performed on single‐electron halogen bonds between methyl radical and bromine‐containing molecules to gain a deeper insight into the nature of such noncovalent interactions. Bader's atoms in molecules (AIM) theory have also been applied to the analysis of the linking of the single‐electron halogen bond. Various characteristics of the R? Br…CH₃ interaction, i.e., binding energies, geometrical parameters and topological properties of the electron density have been determined. The presence of the bond critical points (BCPs) between the bromine atom and methyl radical and the values of electron density and Laplacian of electron density at these BCPs indicate the closed‐shell interactions in the complexes. The single‐electron halogen bonds, which are significantly weaker than the normal halogen bonds, exhibit equally bond strength as compared to the single‐electron hydrogen bond. It has been also found that plotting of the binding energies versus topological properties of the electron density at the BCPs gives two straight lines. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Three-center-two-electron and four-center-four-electron bonds. A study by electron charge density over the structure of methonium cations

We study the electronic density charge topology of CH(5)(+) species 1 (C(s)()), 2 (C(s)()), and 3 (C(2)(v)) at ab initio level using the theory of atoms in molecules developed by Bader. Despite the reports of previous studies concerning carbocationic species, the methane molecule is protonated at the carbon atom, which clearly shows its pentacoordination. In addition to the fact that hydrogen atoms in the methonium molecule behave in a very fluxional fashion and that the energy difference among the species 1, 2, and 3 are very low, is important to point out that two different topological situations can be defined on the basis of our study of the topology of the electronic charge density. Then, the species 1 and 2 present a three-center-two-electron (3c-2e) bond of singular characteristics as compared with other carbocationic species, but in the species 3, the absence of a 3c-2e bond is noteworthy. This structure can be characterized through the three bond critical points found, corresponding to saddle points on the path bonds between the C-H(2,3,5) that lie in the same plane. These nuclei define a four-center interaction where the electronic delocalization produced among the sigma(C-H) bonds provide a stabilization of the three C-H bonds involved in this interaction (the remaining two C-H bonds are similar to those belonging to the nonprotonated species). Our results show that bonding situations with a higher number of atom arrays are possible in protonated hydrocarbons.     

Theoretical studies on the substitution effect of fluorine on ethylene

The substitution effect of fluorine on ethylene is investigated by means of studyingthe properties of the charge distribution at the bond critical points with the theory of atomsin molecules.It is found that fluorine atom acts not only as a σ electron acceptor,but also asa π electron donor,and these double effects are reflected in the quantity of ellipticity,Lap-lacian and the charge density of charge distribution at the bond critical points.For C—C,C—Fbonds,the major axis of elliptical contours is perpendicular to the molecular plane,but forC—H bond,it is parallel to the molecular plane.Other effects originating from the substi-tution have also been discussed.     

Experimental and theoretical charge density distribution of the colossal magnetoresistive transition metal sulfide FeCr2S4

The total charge density distribution rho(r) of the colossal magnetoresistive transition metal sulfide FeCr(2)S(4) was evaluated through a multipole formalism from a set of structure factors obtained both experimentally, by means of single crystal high-quality x-ray diffraction data collected at T=23 K, and theoretically, with an extended-basis unrestricted Hartree-Fock periodic calculation on the experimental geometry. A full topological analysis, followed by the calculation of local energy density values and net atomic charges, was performed using the quantum theory of atoms in molecules. The experimental and theoretical results were compared. Good agreement was found for the topological properties of the system, as well as for the atomic net charges and the nature of the chemical bonds. An analysis of the electron density rho(r), its Laplacian nabla(2)[rho(r)], and the total energy density H(r) at the bond critical points was employed to classify all the interactions that resulted as predominantly closed shell (ionic) in nature. The topological indicators of the bonded interactions for Fe are distinct from those for Cr. The Fe-S bond distances were found to be 0.145 A shorter than the ideal values computed on the basis of Shannon's crystal radii, much shorter than the Cr-S distances with respect to their ideal Shannon lengths. Concomitantly, rho(r) and |H(r)| at the bond critical points are greater for Fe-S interactions, indicating that the local concentration of charge density in the internuclear region is larger for the tetrahedrally coordinated iron than for the octahedrally coordinated chromium. The isosurface in the real space for nabla(2)[rho(r)]=0 was plotted for both iron and chromium, pointing out the local zones of valence shell charge concentration and relating them to the partial d-orbital occupancy of the two transition metal atoms.     

Electron density and energy decomposition analysis in hydrogen-bonded complexes of azabenzenes with water, acetamide, and thioacetamide

Ab initio and density functional theoretical studies on hydrogen-bonded complexes of azabenzenes with water, acetamide, and thioacetamide have been carried out to explore the controversy involved in the relative order of their stability in a systematic way. The interaction energies of these complexes have been analyzed using the Morokuma energy decomposition method, and the nature of the various hydrogen bonds formed has been investigated through topological aspects using Bader's atom in a molecule (AIM) theory. Morokuma energy decomposition analysis reveals that the major contributions to the energetics are from the polarization (PL) and charge transfer (CT) energies. From the calculated topological results, excellent linear correlation is shown to exist between the hydrogen-bond length, electron density [rho(r)], and its Laplacian [nabla(2)rho(r)] at the bond critical points for all the complexes considered.     

Continua of interactions between pairs of atoms in molecular crystals

The electron density distributions in crystals of five previously studied DMAN complexes and five Schiff bases (two new ones) have been analysed in terms of various properties of bond critical points (BCPs) found in the pair-wise interactions in their lattices. We analysed the continua of interactions including covalent/ionic bonds as well as hydrogen bonds and all other types of weak interactions for all pairs of interacting atoms. The charge density at BCPs and local kinetic and potential energy densities vary exponentially with internuclear distance (or other measures of separation). The parameters of the dependences appear to be characteristics of particular pairs of atom types. The Laplacian and the total (sum of kinetic and potential) energy density at BCPs show similar behaviour with the dependence being of the Morse type. The components lambda1, lambda2, lambda3 of the Laplacian at BCPs vary systematically with internuclear distance according to the type of atom pair. For lambda1 and lambda2 the distribution is of the exponential type, whereas lambda3 does not seem to follow any simple functional form, consistent with previous theoretical findings. Analytical nonlinear dependences of Laplacian on charge density have been found. They agree reasonably well with those obtained by least square fit of the Laplacian to charge density data. There are four distinct regions of the [symbol: see text]2rho(BCP)/rho(BCP) space, generated by E(BCP) = 0 and G(BCP)/rho(BCP) = 1 conditions. Two regions clearly correspond to the shared-shell and closed-shell interactions and the other two to some intermediate situation.     

A topological analysis of electron density and chemical bonding in cyclophosphazenes PnNnX2n (X = H, F, Cl; n = 2, 3, 4)

The restricted Hartree-Fock method was used to determine the cycle size effects on the geometric parameters of several inorganic templates, cyclophosphazenes P_nN_nX_2n (X = H, F, Cl; n = 2, 3, 4). A topological analysis of local electronic properties at the electron density critical points of bonds allowed us to quantitatively characterize the chemical bond in cyclophosphazenes and its dependence on the cycle size and substituents at phosphorus. The calculated distributions of the electron density Laplacian and electron pair localization functions revealed the special features of the electronic structure of the nitrogen and phosphorus atoms. These results explain the nature of noncovalent interactions between the P atoms of one cyclophosphazene molecule and the N atoms of the other.     

Electron density distribution of an oxamato bridged Mn(II)-Cu(II) bimetallic chain and correlation to magnetic properties

The electron density distribution of the ferrimagnetic MnCu(pba)(H2O)3.2H2O chain compound, where pba stands for 1,3-propylenebis(oxamato), has been derived from high resolution X-ray diffraction measurements at 114 K using a multipolar model. The analysis of the chemical bonding has been carried out through the "Atoms in Molecules" formalism and thoroughly interpreted with regards to the strong intrachain and weak interchain magnetic couplings. The topological properties of the electron density on the oxamato bridge indicate large electron delocalization and conjugation effects, in addition to high charge transfer from both metals to the bridge. The resulting positive charges on Mn (+1.45 e) and Cu (+1.56 e) induce charge polarization of the bridge, leading to a shift of electron density from the central C atoms to the metal coordinating O and N atoms. The Mn-bridge interactions are mainly closed-shell interactions with low electron density at the corresponding bond critical points, whereas the Cu-bridge interactions exhibit significant covalent character. The Cu-N bonds are moreover stronger than the Cu-O bonds. The 3d Cu and Mn orbital populations are consistent with pyramidal and regular octahedral environments, respectively, in agreement with the loss of degeneracy due to ligand field effects. Interchain interaction pathways are evidenced by the existence of four bond critical points in hydrogen bond regions. Finally, these intrachain and interchain bonding features are correlated to the results of experimental and theoretical spin density distributions, as well as magnetic measurements.     

Experimental electron density analysis of Mn2(CO)10: metal-metal and metal-ligand bond characterization

The experimental electron density rho(r) of Mn2(CO)10 was determined by a multipole analysis of accurate X-ray diffraction data at 120 K. The quantum theory of atoms in molecules (QTAM) was applied to rho(r) and its Laplacian [symbol: see text] 2 rho(r). The QTAM analysis of rho(r) showed the presence of a bond critical point (rc); its associated bond path connects the two Mn atoms, but no cross interaction line was found between one manganese and the equatorial carbonyls of the other. The distribution of [symbol: see text] 2 rho(r) indicated "closed-shell" interactions for the metallic Mn-Mn bond and the dative Mn-CO bonds. The values of the topological parameters of the density at rc, rho(rc), [symbol: see text] 2 rho(rc), G(rc) (kinetic energy density), and V(rc) (potential energy density), characterize the bonds and are intermediate to those corresponding to typical ionic and covalent bonds.