Molecular structural formulas as one-electron density and hamiltonian operators: the VIF method extended

The valency interaction formula (VIF) method is given a broader and more general interpretation in which these simple molecular structural formulas implicitly include all overlaps between valence atomic orbitals even for interactions not drawn in the VIF picture. This applies for VIF pictures as one-electron Hamiltonian operators as well as VIF pictures as one-electron density operators that constitute a new implementation of the VIF method simpler in its application and more accurate in its results than previous approaches. A procedure for estimating elements of the effective charge density-bond order matrix, Pmunu, from electron configurations in atoms is presented, and it is shown how these lead to loop and line constants in the VIF picture. From these structural formulas, one finds the number of singly, doubly, and unoccupied molecular orbitals, as well as the number of molecular orbitals with energy lower, equal, and higher than -1/2Eh, the negative of the hydrogen atom's ionization energy. The VIF results for water are in qualitative agreement with MP2/6311++G3df3pd, MO energy levels where the simple VIF for water presented in the earlier literature does not agree with computed energy levels. The method presented here gives the simplest accurate VIF pictures for hydrocarbons. It is shown how VIF can be used to predict thermal barriers to chemical reactions. Insertion of singlet carbene into H2 is given as an example. VIF pictures as one-electron density operators describe the ground-state multiplicities of B2, N2, and O2 molecules and as one-electron Hamiltonian operators give the correct electronegativity trend across period two. Previous implementations of VIF do not indicate singly occupied molecular orbitals directly from the pictorial VIF rules for these examples. The direct comparison between structural formulas that represent electron density and those that represent energy is supported by comparison of a simple electronegativity scale, chiD=N/n2, with well-known electronegativity scales of Pauling, Mulliken, and Allen. This scale comes from the method used to calculate Pmumu for sp3 hybridized period-two elements and is comparable to electronegativity because it has the same form as <1/r> for hydrogenic orbitals. It therefore provides a physical basis for the representation of one electron density and Hamiltonian operators by the same VIF picture.     

Molecular structure, spectroscopic studies and first-order molecular hyperpolarizabilities of p-amino acetanilide

The infrared absorption, Raman spectra and SERS spectra of p-amino acetanilide have been analyzed with the aid of density functional theory calculations at B3LYP/6-311G(d,p) level. The electric dipole moment (mu) and the first hyperpolarizability (beta) values of the investigated molecule have been computed using ab initio quantum mechanical calculations. The calculation results also show that the synthesized molecule might have microscopic nonlinear optical (NLO) behavior with non-zero values. Computed geometries reveal that the PAA molecule is planar, while secondary amide group is twisted with respect to the phenyl ring is found, upon hydrogen bonding. The hyperconjugation of the C=O group with adjacent C-C bond and donor-acceptor interaction associated with the secondary amide have been investigated using computed geometry. The carbonyl stretching band position is found to be influenced by the tendency of phenyl ring to withdraw nitrogen lone pair, intermolecular hydrogen bonding, conjugation and hyperconjugation. The existence of intramolecular C=O...H hydrogen bonded have been investigated by means of the natural bonding orbital (NBO) analysis. The influence of the decrease of N-H and C=O bond orders and increase of C-N bond orders due to donor-acceptor interaction has been identified in the vibrational spectra. The SERS spectral analysis reveals that the large enhancement of in-plane bending, out of plane bending and ring breathing modes in the surface-enhanced Raman scattering spectrum indicates that the molecule is adsorbed on the silver surface in a 'atleast vertical' configuration, with the ring perpendicular to the silver surface.     

Are oxygen rings (On) and their negative ions (O) unstable?

Oxygen clusters of the equilateral ring types O_n and their negative ions O are studied by the recently developed pictorial-topological quantum chemistry method, VIF (valency interaction formulas), of this author. The species are found to be of high energy relative to separated oxygen atoms, the cause being evident from the VIF pictures. The odd n rings should be at local minima but are not likely to have observable negative ions. The even n rings have distortional instabilities, yet their negative ions are more likely to be observed as transients in beam experiments.     

Theoretical investigations on chalcogen-chalcogen interactions: what makes these nonbonded interactions bonding?

To understand the intermolecular interactions between chalcogen centers (O, S, Se, Te), quantum chemical calculations on pairs of model systems were carried out. For the oxygen derivatives, one of the components of the supermolecules consists of dimethyl ether, while the second component is either dimethyl ether (1) or ethynyl methyl ether (2) or methyl cyanate (3). The model calculations were also extended to the sulfur (4-6), selenium (7-9), and tellurium congeners (10-12). The MP2/SDB-cc-pVTZ, 6-311G level of theory was used to derive the geometrical parameters and the global energies of the model systems. A detailed analysis based on symmetry adapted perturbation theory (SAPT) reveals that induction and dispersion forces contribute to the bonding in each case. For 1-3 the electrostatic energy also contributes to the intermolecular bonding, but not for 4-12. The NBO analysis reveals that the interaction in the dimers 1-3 is mainly due to weak hydrogen bonding between methyl groups and chalcogen centers. Similar hydrogen bonding is also found in the case of 4 and to a lesser extent in 5 and 7. For the aggregates with heavier centers the chalcogen-chalcogen interaction dominates, and hydrogen bonding only plays a minor role. Electron-withdrawing groups on the chalcogen centers increase the interaction energy and reduce the intermolecular distance dramatically. The one-electron picture of an interaction between the lone pair of the donor and the chalcogen carbon sigma orbital allows a qualitatively correct reproduction of the observed trend.     

Chemical bonding in view of electron charge density and kinetic energy density descriptors

Stalke's dilemma, stating that different chemical interpretations are obtained when one and the same density is interpreted either by means of natural bond orbital (NBO) and subsequent natural resonance theory (NRT) application or by the quantum theory of atoms in molecules (QTAIM), is reinvestigated. It is shown that within the framework of QTAIM, the question as to whether for a given molecule two atoms are bonded or not is only meaningful in the context of a well‐defined reference geometry. The localized‐orbital‐locator (LOL) is applied to map out patterns in covalent bonding interaction, and produces results that are consistent for a variety of reference geometries. Furthermore, LOL interpretations are in accord with NBO/NRT, and assist in an interpretation in terms of covalent bonding. © 2008 Wiley Periodicals, Inc.J Comput Chem, 2009.     

Theoretical investigations on heteronuclear chalcogen-chalcogen interactions: on the nature of weak bonds between chalcogen centers

To understand the intermolecular interactions between chalcogen centers (O, S, Se, Te), quantum chemical calculations on model systems were carried out. These model systems were pairs of monomers of the composition (CH3)2X1 (X1 = O, S, Se, Te) as the donors and CH3X2Z (with X2 = O, S, Se, Te and Z = Me, CN) as the acceptors. The variation of X1, X2, and Z leads to 32 pairs with 8 homonuclear cases (X1 = X2 = O, S, Se, Te) and 24 heteronuclear cases (X1 not equal X2). The MP2/SDB-cc-pVTZ, 6-311G* level of theory was used to derive the geometrical parameters and the interaction energies of the model systems. The pairs with Z = CN (17-32) show a considerably higher interaction energy than the pairs with CH3 groups only (1-16). Natural bond orbital (NBO) analysis revealed that the interaction of the dimers 1, 2, 5, 6, 9, 10, 13, 14, 17, 21, 25, and 29 is mainly due to weak hydrogen bonding between methyl groups and chalcogen centers. These systems all contain hard chalcogen atoms as acceptors. For all other systems, the chalcogen-chalcogen interaction dominates. The one-electron picture of an interaction between the lone pair of the donor chalcogen atom and the chalcogen-carbon antibonding sigma* orbital serves as a model to qualitatively rationalize trends found in many of these systems. However, it has to be applied with some amount of skepticism. A detailed analysis based on symmetry-adapted perturbation theory (SAPT) reveals that induction and dispersion forces dominate and contribute to the bonding in each case. Hydrogen-bonded compounds involve bonding electrostatic contributions. Compounds dominated by chalcogen-chalcogen interactions exhibit bonding due to electrostatic interactions only if one of the chalcogen atoms involved is sulfur or oxygen.     

铂配合物与DNA碱基对间相互作用的理论研究

用量子化学方法研究一系列Pt(II)配合物作用于嘌呤碱基N7位点后对Watson-Crick碱基对AT、GC的影响. 计算结果显示铂配体与碱基对AT、GC的作用以静电作用为主，同时极化作用也是影响GC碱基对的重要因素. 静电作用极大地增强了铂化嘌呤碱基与嘧啶碱基间的相互作用，而嘌呤碱基与嘧啶碱基间作用与未铂化碱基对作用相近. Pd(II) 和 Ni(II)的相关研究得到类似的结果. 碱基对间氢键作用“二阶微扰能”分析结果与氢键强弱变化一致.     

Ab initio cluster calculations of silica surface sites

Silica surface sites, which can be formed in cleavage processes, and their hydrolyzed counterparts are investigated with ab initio cluster calculations. Natural Bond Orbital (NBO) theory is used to characterize bonding around silica surface sites. Higher energy lone pairs of electrons on oxygen atoms either hyperconjugate to vicinal silanol/siloxane antibonding orbitals or backdonate electron density via donor–acceptor π-type bonding with participation of pd or p hybrids on silicon atoms. Upon substitution of hydroxyl groups of orthosilicic acid with silica monomers the strength of siloxane and silanol Si–O bonding increases as energies of bonding orbitals and contributions from p-orbitals decrease. Silanone sites and a complementary pair of silyl/siloxy radical sites are found to be the most stable geminal and single non-hydrolyzed sites, respectively. Atomic charges based on natural wavefunctions and on fitting to electrostatic potential, and characteristic bands of IR spectra associated with siloxane and silanol stretching vibrations of silica surface sites are reported.     

Dyad algebra and multiplication of graphs. I. Directed graphs

The ket-bra algebra for quantum mechanics and for the quantum chemistry.in valence shells was made by this author fully covariant recently. The resulting principle of linear covariance allowed diverse approaches such as molecular orbital, valence bond, localized orbital theories to come out as special cases in particular basis frames not necessarily orthonormal. The principal also led to the pictorial VIF (valency interaction formula) methods for deducing qualitative quantum chemistry directly from the structural formulas of molecules. The present set of two papers (II on undirected graphs) develops graphs and graph rules for abstract linear vector spaces, bras, kets, and abstract operators as ket-bra dyads. Multiplications of such operators are carried out with graphs of two kinds of lines and two kinds of vertices. The theorems are demonstrated on some examples and are useful, e.g., with the recent method of moments and in deriving Lie algebras pertinent to quantum chemistry.     

Natural bond orbital analysis and vibrational spectroscopic studies of H-bonded N,N′-diphenylguanidinium nitrate

NIR-FT Raman and FT-IR spectra of the crystallized biologically active molecule N,N′-diphenylguanidinium nitrate (DGN) have been recorded and analyzed using quantum chemical computations based on density functional theory. The extraordinary basicity and strong stability of this novel bioactive compound has been discussed as the consequence of resonance stabilization leading to Y-aromaticity and hydrogen bonding. This peculiar Y-delocalization character of DGN is well reflected in the optimized geometry and bond order (BO) calculations. The observance of the equality of C–N bond lengths in the protonated species indicates delocalization of the π-electron system. The spectroscopic and natural bonds orbital (NBO) analysis confirms the occurrence of strong network of inter molecular hydrogen bonds. The changes in electron density in the global minimum and in the energy of hyperconjugative interactions of DGN calculated by second order perturbation theory have been studied extensively in comparison with the values of the neutral species. The observed characteristic ring vibrations are well fit with the theoretical values calculated at B3LYP/6-31G* level.     

A Theoretical Study of NBO,NICS, and 14N NQR Parameters of Adenine Tautomers in the Gas Phase via DFT

The natural bond orbital (NBO) analysis, nucleus independent chemical shift (NICS), and ¹⁴N NQR parameters of the most stable tautomers of adenine in the gas phase were predicted using density functional theory method. The NBO analysis revealed that the resonance interaction between lone pair of the nitrogen atom and empty non‐Lewis NBO increases with increasing the p character of the nitrogen lone pair. The present investigation indicated the π clouds in both the considered heterocyclic rings containing six electrons, and these tautomers has the aromatic character. The NICS study utilizing the gauge‐invariant atomic orbital method showed that there are diatropic currents in the heterocyclic rings of the tautomers, so we determined the order of overall aromaticity of these tautomers. The results of NQR parameter calculations showed three parameters are effective on nuclear quadrupole coupling constant; the p character value of lone pair electrons of nitrogens, and the related occupancies and whenever, the lone pair electrons of nitrogens participate in the formation of chemical bond and/or π system of the ring, the q_zz and consequently its χ decreases.     

Natural bond orbital analysis: A critical overview of relationships to alternative bonding perspectives

We sketch the basic principles of natural bond orbital (NBO) theory, including critical discussion of its relationship to alternative bonding concepts and selected illustrations of its application to a broad spectrum of chemical bonding motifs. Particular emphasis is placed on the close NBO connections to prequantal bonding, and electromerism concepts, as well as the deep roots in quantal eigenvalue, superposition, and Pauli exclusion concepts that are manifested in many aspects of NBO donor–acceptor analysis. With respect to leading alternative perspectives, we identify similarities and differences that distinguish NBO theory from the corresponding precepts of valence bond theory, molecular orbital theory, and Bader's quantum theory of atoms in molecules, with critical discussion of the assumptions underlying characteristic differences in each case. © 2012 Wiley Periodicals, Inc.     

FT-IR and theoretical study of 3,5-dimethyl-1H-pyrazole-1-carboxamidine (L) and the complexes CoL2(H2O)2(NO3)2, NiL2(H2O)2(NO3)2

In the paper a joint experimental and theoretical study of 3,5-dimethyl-1H-pyrazole-1-carboxamidine (L) as well as its complexes CoL2(H2O)2(NO3)2 and NiL2(H2O)2(NO3)2 is reported. On the basis of FT-IR experiments and a DFT-derived scaled quantum mechanical force field the normal coordinate analysis of L was carried out. The FT-IR spectra of the two complexes were interpreted using the present assignment of L and computed vibrational data of the complexes. The ionic and charge transfer interactions in the complexes were assessed by means of natural bond orbital (NBO) analysis.     

Hydrogen bonding and molecular vibrations of 3,5-diamino-1,2,4-triazole

This work deals with the analysis of hydrogen bonding and the vibrational spectroscopy of 3,5-diamino-1,2,4-triazole by means of quantum chemical calculations. The mid and far FTIR and FT-Raman spectra were measured in the condensed state. The fundamental vibrational frequencies were calculated under different possible symmetries by applying the density functional theory with the B3LYP functional and the 6-31G* basis set. The results of the calculations obtained under C(2) symmetry produces the global minimum on the potential energy surface. The vibrational spectra were interpreted with the aid of normal coordinate analysis based on scaled density functional force field. The infrared and Raman spectra were also predicted from the calculated intensities.     

核酸中碱基的氢键作用机理及电子特征理论研究

采用耦合簇量子化学方法 CCSD/aug-cc-pVDZ研究了嘧啶与嘌呤之间的相互作用,利用基函数叠加误差法(BSSE)消除相互作用能误差,并进行了几何结构优化;采用Gaussian 03程序包中的NBO程序分析了二阶稳定化能及自然键轨道.与此同时,应用约化密度函数(RDG)填色等值面图对体系进行了图形化分析,分析了氢键相互作用所在的空间位置和相对强度,以及氢键相互作用的性质,以进一步了解二者的相互作用.结果表明,嘧啶-嘌呤体系的相互作用属于闭合壳层静电相互作用.电子密度跃迁矩阵分析结果表明,激发区域主要集中在N原子和O原子处,涉及的空间广度很大,第一激发态主要涉及前线分子轨道,属于σ→π*或n→π*类型跃迁.     

On the quantum chemical origin for the nonvalidity of Koopmans' theorem in transitionmetal compounds

The quantum chemical origin for the nonvalidity of Koopmans' theorem in transitionmetal compounds of the 3d series is analyzed by means of the Green's function formalism applied in the framework of a semiempirical INDO Hamiltonian. In the case of ferrocene (1), cyclobutadiene iron tricarbonyl (2) and irontetracarbonyl dihydride (3) the self-energy part of a geometric approximation has been partitioned into relaxation and correlation (pair removal, pair relaxation) increments. The breakdown of Koopmans' theorem for strongly localized MOs with large Fe 3d amplitudes is predominantly the result of electronic relaxation lowering the calculated ionization potentials. On the other hand the variation of the pair correlation energy in the cationic hole-state is by no means negligible and acts into the opposite direction as the relaxation increment. These significant pair relaxation contributions explain the wellknown failtures of the ΔSCF approach in combination with large scaleab initio bases. The loss of ground state pair correlation in the outer valence region is small in comparison to relaxation and pair relaxation. The magnitude of the aforementioned reorganization increments has been studied as a function of the localization properties of the MOs and as a function of the one-electron energies of the available particle- and hole-states. The computational findings derived with the INDO model are compared with recentab initio studies.     

Toward a physical understanding of electron-sharing two-center bonds. I. General aspects

In 1916, Lewis and Kossel laid the empirical ground for the electronic theory of valence, whose quantum theoretical foundation was uncovered only slowly. We can now base the classification of the various traditional chemical bond types in a threefold manner on the one- and two-electron terms of the quantum-physical Hamiltonian (kinetic, atomic core attraction, electron repulsion). Bond formation is explained by splitting up the real process into two physical steps: (i) interaction of undeformed atoms and (ii) relaxation of this nonstationary system. We aim at a flexible bond energy partitioning scheme that can avoid cancellation of large terms of opposite sign. The driving force of covalent bonding is a lowering of the quantum kinetic energy density by sharing. The driving force of heteropolar bonding is a lowering of potential energy density by charge rearrangement in the valence shell. Although both mechanisms are quantum mechanical in nature, we can easily visualize them, since they are of one-electron type. They are however tempered by two-electron correlations. The richness of chemistry, owing to the diversity of atomic cores and valence shells, becomes intuitively understandable with the help of effective core pseudopotentials for the valence shells. Common conceptual difficulties in understanding chemical bonds arise from quantum kinematic aspects as well as from paradoxical though classical relaxation phenomena. On this conceptual basis, a dozen different bond types in diatomic molecules will be analyzed in the following article. We can therefore examine common features as well as specific differences of various bonding mechanisms.     

Properties of a new fluorescent cytosine analogue, pyrrolocytosine

Pyrrolocytosine is a novel, environment sensitive, fluorescent base that can be used in place of cytosine as a fluorescent marker in nucleic acids. In this work the results of a detailed computational investigation into the hybridization and photochemical properties of the base are reported. The interaction energy of the base pair formed between pyrrolocytosine and guanine, calculated at the MP2/6-31G(d,0.25)//HF/6-31G(d,p) level, was found to be -27.2 kcal mol(-1), comparing very favorably with the value calculated for the cytosine and guanine base pair, -25.8 kcal mol(-1). The wavelengths for the vertical transitions of pyrrolocytosine and cytosine were determined using both the configuration interaction technique, with singly excited configurations (CIS) and time-dependent density functional theory using the B3LYP functional (TDB3LYP). It was found that the spacing between the first pipi state and the first npi state was significantly larger in the case of pyrrolocytosine than cytosine, providing a rationale for the higher fluorescence quantum yield of the former. Hydrogen bonding of pyrrolocytosine to guanine did not affect the predicted fluorescence properties of pyrrolocytosine whereas stacking guanine above pyrrolocytosine, in a manner appropriate to B-form DNA, significantly reduced the predicted fluorescence. Calculations on the two base systems using the TDB3LYP method produced low-lying charge-transfer states which are not predicted when the CIS method is used and are not thought to be physically meaningful.     

On isomers and tautomers of Nitro-1-deazapurine: A DFT study

The 1-deazapurine molecules substituted by the NO₂ group at three different positions of the six-membered ring were subject of computational study performed at the B3LYP/cc-pVTZ level. For each substitution three tautomers were considered. The N3H tautomers are relatively instable in the gas and water phases which is due to significant decrease in aromaticity of the N3H forms. In the gas phase, the equilibria between the N7H and N9H tautomers are determined by a competition of the repulsive and attractive intermolecular interaction of different moieties of Nitro-1-deazapurine. The close neighborhood of the two tertiary N atoms and attractive close neighborhood of the tertiary N atom and the NH group result in preference of the N9H tautomer over the N7H one by ca. 3.5 kcal/mol. By comparing energy of different forms and proposed isodesmic reaction we showed that the Gibbs free energy of the attractive interaction between the NO₂ and HN groups is equal to ca. 1.0 kcal/mol, whereas the repulsive interaction between the NO₂ group and the tertiary N atom of the imidazole is equal to ca. 6.4 kcal/mol. It was shown also that the increase in dipole moment in the water media is the crucial effect influencing the N7H/N9H tautomeric equilibria of Nitro-1-deazapurines. For the three isomers dissolved in water, the two tautomers, N7H and N9H, are predicted to be observed and the former should dominate slightly for the 2-Nitro isomer whereas the latter for the 6-Nitro and 1-Nitro isomers. The NBO analysis showed that the NO₂ group withdraws 0.3 e of the σ-electrons from the pyridine ring of 1-deazapurine and has no influence on the σ-electrons of the imidazole ring. The NBO analysis shows also that the ratio between number of π-electrons which are withdrawed from imidazole vs. pyridine ring is characteristic for position of the substitution. The isodesmic reaction used revealed also that the NO₂ group destabilizes the N3H-1-deazapurine system, whereas it stabilizes the N7H- and N9H-1-deazapurine systems.     

Structural,electronic properties and intramolecular hydrogen bonding of substituted 2-[(E)-imino methyl] benzenethiol in ground and first excited state by quantum chemical methods

Quantum chemical calculations of geometric structure, the intramolecular hydrogen bond, harmonic vibrational frequencies, NMR spin–spin coupling constants, and physical properties such as chemical potential and chemical hardness of the 2-(E)-imino methyl benzenethiol and its nineteen derivatives were carried out using density functional theory (DFT/B3LYP/6-311++G**) method in the gas phase and the water solution. Furthermore, the topological properties of the electron density distributions for S–H···N intramolecular hydrogen bond have been analyzed in terms of the Bader’s theory of atoms in molecules (AIM). Natural bond orbital (NBO) analysis also performed for better understanding the nature of intramolecular interactions, the results of analysis by quantum theory of AIM and NBO method fairly supported the DFT results. Besides, MEP was performed by the DFT method. On the other hand, the aromaticity of the formed ring has been measured using several well-established indices of aromaticity such as nucleus-independent chemical shift, harmonic oscillator models of the aromaticity, para-delocalization index, average two-center indices, and aromatic fluctuation index. Also, the excited-state properties of intramolecular hydrogen bonding in these systems have been investigated theoretically using the time-dependent DFT method.