Toward a semiempirical density functional theory of chemical binding

The new ideas ofbond electronegativity andbond hardness are introduced, and a semiempirical density functional approach to the theory of molecular electronic structure and chemical binding is outlined. There result effective electronegativity equalization procedures that permit calculation of binding energies as well as partial charges. By a modelling of the bond electronegativity and bond hardness, a density functional interpretation of earlier bond charge models is established. Some numerical results are given for diatomic molecules.Dedicated to Professor J. Koutecký on the occasion of his 65th birthday     

Solvation effect on conformations of 1,2:dimethoxyethane: charge-dependent nonlinear response in implicit solvent models

The physical content of and, in particular, the nonlinear contributions from the Langevin-Debye model are illustrated using two applications. First, we provide an improvement in the Langevin-Debye model currently used in some implicit solvent models for computer simulations of solvation free energies of small organic molecules, as well as of biomolecular folding and binding. The analysis is based on the implementation of a charge-dependent Langevin-Debye (qLD) model that is modified by subsequent corrections due to Onsager and Kirkwood. Second, the physical content of the model is elucidated by discussing the general treatment within the LD model of the self-energy of a charge submerged in a dielectric medium for three different limiting conditions and by considering the nonlinear response of the medium. The modified qLD model is used to refine an implicit solvent model (previously applied to protein dynamics). The predictions of the modified implicit solvent model are compared with those from explicit solvent molecular dynamics simulations for the equilibrium conformational populations of 1,2-dimethoxyethane (DME), which is the shortest ether molecule to reproduce the local conformational properties of polyethylene oxide, a polymer with tremendous technological importance and a wide variety of applications. Because the conformational population preferences of DME change dramatically upon solvation, DME is a good test case to validate our modified qLD model. The present analysis of the modified qLD model provides the motivation and tools for studying a wide variety of other interesting systems with heterogeneous dielectric properties and spatial anisotropy.     

On the relation between electronic and nuclear vibrations via the response matrix derived from semiempirical density functional-based sensitivity analysis

Energetic and physical aspects of the intuitive (mapping) relations between independent (normal) charge polarization patterns and nuclear vibration modes are derived within the atoms-in-molecules (AIM ) electronegativity equalization method (EEM ) formulation of the charge sensitivity analysis (CSA ). The concepts are demonstrated with calculations on homologous series of small organic molecules. © 1995 John Wiley & Sons, Inc.     

Spin density calculations in semiempirical restricted and unrestricted SCF MO methods

In the present paper spin densities are calculated by ten methods on chosen radical systems. - interaction parameters for the individual methods are determined by McConnel's relation a _i =Q _i and by that of Colpa-Bolton a _i =(Q+K _i ) _i and the correlation between the results of the individual methods is studied.

---

Zusammenfassung In der vorliegenden Arbeit werden Spindichten anhand von zehn Methoden an ausgewählten Radikalsystemen ermittelt. Die --Kopplungsparameter für die einzelnen Methoden werden nach den McConnell a _i =Q _i - und Colpa-Bolton a _i =(Q+K _i ) _i -Beziehungen bestimmt. Die zwischen den Ergebnissen der einzelnen Methoden bestehende Korrelation wird untersucht.

---

Résumé Dans ce travail les densitées de spin sont calculées par dix méthodes sur un système des radicaux choisis. On détermine des parametres - d'interaction pour les méthodes individuelles dans la relation de McConnell a _i =Q _i et celle de Colpa-Bolton a _i =(Q+K _i ) _i et étudie la correlation entre les résultats des méthodes individuelles.

     

Magnetic effects of disulfide bridges: a density functional and semiempirical study

Density functional chemical shielding calculations are reported for methane and hydrogen disulfide dimers. The calculations show that the contributions of disulfide bridges to the chemical shielding of neighboring protons is sizable at distances that are frequently sampled in protein structures. A semiempirical model of the quantum chemical data is developed. It is shown that magnetic anisotropy effects of disulfide are poorly described by the McConnell equation, both qualitatively and quantitatively. In particular, the ratio of magnetic anisotropy contributions to shielding along and perpendicular to the magnetic anisotropy principal axis do not conform to the predictions of the McConnell equation, and magnetic anisotropy effects are not null along the magic angle axis. A sulfur-based model of the magnetic anisotropy of the disulfide is developed and shown to give much better agreement with the quantum chemical data.     

Looking at self-consistent-charge density functional tight binding from a semiempirical perspective

The self-consistent-charge density functional tight binding (SCC-DFTB) method is compared with other semiempirical methods (MNDO, AM1, PM3, OM1, OM2, OM3). Despite the differences in the underlying philosophy and derivation, these methods share many common features. Systematic evaluations of their performance are reported for standard test sets that are in common use. The overall accuracy of SCC-DFTB and the other methods is in the same range, with the overall tendency AM1相似文献   

Application of semiempirical long-range dispersion corrections to periodic systems in density functional theory

Ewald summation is used to apply semiempirical long-range dispersion corrections (Grimme, J Comput Chem 2006, 27, 1787; 2004, 25, 1463) to periodic systems in density functional theory. Using the parameters determined before for molecules and the Perdew-Burke-Ernzerhof functional, structure parameters and binding energies for solid methane, graphite, and vanadium pentoxide are determined in close agreement with observed values. For methane, a lattice constant a of 580 pm and a sublimation energy of 11 kJ mol(-1) are calculated. For the layered solids graphite and vanadia, the interlayer distances are 320 pm and 450 pm, respectively, whereas the graphite interlayer energy is -5.5 kJ mol(-1) per carbon atom and layer. Only when adding the semiempirical dispersion corrections, realistic values are obtained for the energies of adsorption of C(4) alkenes in microporous silica (-66 to -73 kJ mol(-1)) and the adsorption and chemisorption (alkoxide formation) of isobutene on acidic sites in the micropores of zeolite ferrierite (-78 to -94 kJ mol(-1)). As expected, errors due to missing self-interaction correction as in the energy for the proton transfer from the acidic site to the alkene forming a carbenium ion are not affected by the dispersion term. The adsorption and reaction energies are compared with the results from M?ller-Plesset second-order perturbation theory with basis set extrapolation.     

Theory of relaxed density matrices: Application to second-order response properties

The evaluation of second derivatives of the electronic energy for nonvariational wave functions using an energy functional is discussed. It is shown that, in certain cases, the formation of the first-order relaxed density matrix leads to an efficient algorithm for the calculation of second-order response properties. Detailed formulas are given for second-order Møller–Plesset perturbation theory. © 1994 John Wiley & Sons, Inc.     

Nanomechanical properties of mechanical double-layers: a novel semiempirical analysis

Force-displacement curves have been acquired with a commercial atomic force microscope on a thin film of poly(n-butyl methacrylate) on glass substrates. The film thickness is nonuniform, ranging in the measured area from 0 to 30 nm, and gives the possibility to survey the so-called "mechanical double-layer" topic, i.e., the influence of the substrate on the mechanical properties of the film in dependence of the film thickness. The stiffness and the deformation for each force-distance curve were determined and related to the film thickness. We were able to estimate the resolution of the film thickness that can be achieved by means of force-distance curves. By exploiting the data acquired in the present and in a previous experiment, a novel semiempirical approach to describe the mechanical properties of a mechanical double-layer is introduced. The mathematical model, with which deformation-force curves can be described, permits to calculate the Young's moduli of film and substrate in agreement with literature values and to determine the film thickness in agreement with the topography.     

Frequency-dependent response properties and excitation energies from one-electron density matrix functionals

The recent formulation of the time-dependent density matrix functional theory (TD-DMFT) has opened an avenue to calculations of frequency-dependent response properties and excitation energies of atoms and molecules. In practice, the accuracy of the computed data is limited by both the errors inherent to the adiabatic approximation or its modifications and the quality of the energy functionals. The relative importance of these two factors is carefully assessed with test calculations on diatomic molecules with few electrons. The test results clearly demonstrate the superiority of an ad hoc approach that corrects the improper behavior of the adiabatic approximation at the low-frequency limit. Even more importantly, TD-DMFT convincingly removes the ambiguity in the choice of the two-electron integrals that is present in the stationary-state case. On the other hand, paralleling the previously reached conclusions pertinent to ionization potentials, the presently available BBC-type functionals are found to be insufficiently accurate to provide reliable quantitative predictions of excitation energies.     

A new scaling procedure to correct semiempirical MEP and MEP-derived properties

Summary The MNDO, AM1, PM3, and ab initio 6–31G^* and 6–31+G^* MEPs for 21 neutral and 12 charged molecules were computed in layers ranging from 1.2 to 2.0 times the van der Waals radii of atoms. Semiempirical and ab initio MEPs for each layer and two groups of layers were compared to gain insight into the relationships between semiempirical and ab initio MEPs. A detailed statistical study allowed us to obtain a new set of scaling coefficients able to correct the semiempirical MEPs to provide better representations of the ab initio values. The corrected semiempirical MEPs were used to obtain electrostatic charges, whose quality was tested by the comparison between semiempirical Coulombic MEPs and ab initio quantum mechanical MEPs.     

Mechanical properties of blends of low density with linear low density polyethylene

Torque values and mechanical properties are reported for blends of low density with linear low density polyethylene (PE). The torque values show that incompatible mixtures in the molten state are formed when the torques of the pure polymers are almost equal. Semicompatible behaviour is shown in the other cases. The mechanical properties indicate that semicompatible mixtures are formed in the solid state. The mechanical properties of these blends are strongly influenced by the linear low density polyethylene only for contents exceeding 25%.     

Pair density analysis of semiempirical wave functions: Study of water molecule including environmental effects

The pair-density-analysis method is applied to the study of CNDO /2 wave functions of the water molecule in various configurations. Particular attention is devoted to the modifications introduced by the consideration of environmental effects. The validity of the separated pair approximation is discussed and a detailed analysis of group functions is given.     

多氰基立方烷生成热的DFT-B3LYP和半经验MO研究

运用密度泛函理论(DFT)B3LYP方法和半经验MO(MINDO/3,MNDO,AM1和PM3)方法系统计算了全部21种多氰基立方烷的生成热,首先,在DFT-B3LYP/6-31G^*水平下通过不破裂立方烷笼状骨架(亦即选择立方烷为参考物)的等键反应设计,精确计算了9种多氰基立方烷的生成热;发现B3LYP/6-31G^*结果分别地均与上述四种半经验MO方法求得的生成热之间存在良好的线性关系(相关系数均在0.9994以上),且以AM1生成热与B3LYP/6-31G^*计算值最为接近,其次,其它12种多氰基立方烷的精确生成热借助上述线性关系通过校正对应的AM1结果而获得,多氰基立方烷的生成热很高,且随-CN基数目的增加而线性地增大,表明它们属于极具潜力的“新一低高能炸药”而具开发价值。     

Periodic calculations of excited state properties for solids using a semiempirical approach

The semiempirical SCF MO method MSINDO (modified symmetrically orthogonalized intermediate neglect of differential overlap) [T. Bredow and K. Jug, Electronic Encyclopedia of Computational Chemistry, 2004] is extended to the calculation of excited state properties through implementation of the configuration interaction singles (CIS) approach. MSINDO allows the calculation of periodic systems via the cyclic cluster model (CCM) [T. Bredow et al., J. Comput. Chem., 2001, 22, 89] which is a direct-space approach and therefore can be in principle combined with all molecular quantum-chemical techniques. The CIS equations are solved for a cluster with periodic boundary conditions using the Davidson-Liu iterative block diagonalization approach. As a proof-of-principle, MSINDO-CCM-CIS is applied for the calculation of optical spectra of ZnO and TiO(2), oxygen-defective rutile, and F-centers in NaCl. The calculated spectra are compared to available experimental and theoretical literature data. After re-adjustment of the empirical parameters the quantitative agreement with experiment is satisfactory. The present approximate approach is one of the first examples of a quantum-chemical methodology for solids where excited states are correctly described as n-electron state functions. After careful benchmark testing it will allow calculation of photophysical and photochemical processes relevant to materials science and catalysis.     

Photophysical properties and semiempirical calculations of perylene-3,4,9,10-tetracarboxylic tetramethylester (PTME)

The spectral behavior and fluorescence quantum yield of perylene-3,4,9,10-tetracarboxylic tetramethylester (PTME) have been measured in different solvents. Both electronic absorption and fluorescence spectra are not sensitive to medium polarity. The dye exhibits high fluorescence quantum yield and high photostable. Crystalline solid of PTME gives excimer-like emission at 530 nm. The laser activity of PTME has been investigated. The dye solution in N,N-dimethylformamide (DMF) gives laser emission around 480 nm upon excitation by 337.1 nm nitrogen laser pulse. The excitation energy transfer from 7-dimethylamino-4-methylcoumarine (DMC) to PTME has also has been studied and the value of energy transfer rate constant, k(ET), and critical transfer distance, R(0) indicate a F?rster-type mechanism. The photodecomposition of PTME in chloromethane solvents has been also studied. We applied semiempirical MO calculations using (PM3 and ZINDO-CI) calculations to explain the geometric and electronic behaviors of the PTME molecule in both ground and excited states and make a correlation with the experimental observations.