Ab initio studies of the ground-state potential energy surface of formamide

Ab initio calculations employing a standard double-zeta basis set augmented with various polarization functions have been used to investigate the lowest energy region of the ground-state potential energy surface of the formamide molecule. Hartree-Fock calculations with only d polarization functions on the nonhydrogen atoms located two stable minima, that with geometry distorted from planarity having slightly lower energy; only one stable minimum with planar structure is found when p polarization functions on the hydrogens are included. In contrast optimizations, which account approximately for the correlation energy using second-order Møller-Plesset perturbation theory consistently favor a single slightly nonplanar minimum energy geometry.     

Long-range interaction between 1s and 2s or 2p hydrogen atoms

Long-range interaction energy between two hydrogen atoms has been computed in the second order of the perturbation theory. All states of the system arising when one of the atoms is in the 1s and the other in the 2s or 2p state have been considered. The energy represented by a series expansion in inverse powers of the internuclear distance, R, has been computed up to the terms in R^?8. The results are believed to give reliable interaction energies for R > 15 a.u. Accurate interaction energy for two ground-state hydrogen atoms has also been obtained up to the terms in R^?10. Results for the B′ ¹∑ state are employed to discuss the experimental ground-state dissociation energy of H₂, D₂, and HD. For H₂ all values of the dissociation energy obtained from various experimental absorption limits, by using the computed potential energy curve to separate off the effect of rotation, are shown to be satisfactorily consistent. The resulting total energy of H₂ is, however, higher than the most accurate theoretical value.     

Molecular structures of acetylene derivatives of tin

The geometry and force fields of the bis(trimethylstannyl)acetylene molecule (a conformer withD _3d symmetry corresponding to a minimum of the total energy of the molecule) were calculated by the RHF and MP2(fc) methods. The effective core potential in SBK form with the optimized 31G^* valence basis set was employed in the case of Sn atoms. The 6–31G^** and 6–311G^** basis sets were used for carbon and hydrogen atoms. Vibrational spectra of the light and perdeuterated isotopomers of bis(trimethylstannyl)acetylene were interpreted using the procedure of scaling the quantum-chemical force fields. For Part 5, see Ref. 1. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 4, pp. 616–626, April, 2000.     

Intramolecular nucleophilicS N 2 substitution at the tetrahedral carbon atom: Anab initio study

Pentacoordination of carbon atom in bicyclic organic compounds of the pentalene type was studied by theab initio RHF/6-31G^** and MP2(full)/6-31G^** methods. It was shown that intramolecularS _N 2 reactions with energy barriers within the energy scale of NMR spectroscopy can occur in systems in which a linear orientation of the attacking and leaving groups is realized. The barrier to the intramolecular nucleophilic substitution reaction in 2,3-dihydro-3-formylmethylenefuran is 36.9 (RHF) and 27.7 kcal mol⁻¹ (MP2) and decreases to 16.4 and 19.4 kcal mol⁻¹, respectively, in the case of diprotonation at the O atoms in this system. For model pentalene type compounds containing electron-deficient B atoms in the ring, theab initio calculations predict a further decrease in the barrier height (down to less than 10 kcal mol⁻¹). Translated fromIzvestiya Akademii Nauk, Seriya Khimicheskaya, No. 7, pp. 1246–1256, July, 1999.     

Theoretical calculation of the low laying electronic states of the molecular ion KH

Using an ab initio method the potential energy has been calculated for the 25 lowest molecular states of symmetries ²Σ⁺, ²Π, ²Δ for the molecular ion KH⁺. The calculation is based on nonempirical pseudopotentials and parameterized ℓ-dependent polarization potentials. Gaussian basis sets have been used for both atoms. The spectroscopic constants for 18 electronic sates have been calculated by fitting the calculated energy values to a polynomial in terms of the internuclear distance R. Through the canonical functions approach the eigenvalue E_v, the abscissas of the corresponding turning points (R_min and R_max) and the rotational constants B_v have been calculated up to 24 vibrational levels for the considered bound states. The comparison of the present results with those available in literature shows a very good agreement.     

Crash-Conformational Kinematics of the Covalent Lattice of Diamond on Melting

Using the model potential approach we have obtained a generalized potential of interatomic interactions for two condensed states providing an abrupt change in thermodynamic functions at a point of a solid to liquid transition. A kinematic model is considered for the transition of the diamond lattice with pure covalent bonds to the melted state. The model is based on the concept of crash destruction of the lattice during the conformational transformation of the initial symmetry of the chosen atomic group. The kinetic energy of the thermal motions of atoms does not exceed the potential energy of lattice atomization. The kinematic patterns of the positions of 150 atoms at T T _m and T T _m are given for the two-dimensional hexagonal net.     

Statistical Adiabatic Channel Model Calculations of the Reaction of H Atoms and FO Radicals on Ab Initio Potential Energy Surfaces

The statistical adiabatic channel model has been used to calculate rate coefficients for the reaction between H atoms and FO radicals on potential energy surfaces based on density functional theory and ab initio calculations between 200 and 1000 K. The rate coefficient calculated at 300 K with the ab initio potential is in good agreement with a recently reported experimental value.     

Wavelength-decomposition-based embedded cluster density approximation for systems with nonlocal electron correlation

Local correlation methods rely on the assumption that electron correlation is nearsighted. In this work, we develop a method to alleviate this assumption. This new method is demonstrated by calculating the random phase approximation (RPA) correlation energies in several one-dimensional model systems. In this new method, the first step is to approximately decompose the RPA correlation energy to the nearsighted and farsighted components based on the wavelength decomposition of electron correlation developed by Langreth and Perdew. The short-wavelength (SW) component of the RPA correlation energy is then considered to be nearsighted, and the long-wavelength (LW) component of the RPA correlation energy is considered to be farsighted. The SW RPA correlation energy is calculated using a recently developed local correlation method: the embedded cluster density approximation (ECDA). The LW RPA correlation energy is calculated globally based on the system's Kohn-Sham orbitals. This new method is termed λ-ECDA, where λ indicates the wavelength decomposition. The performance of λ-ECDA is examined on a one-dimensional model system: a H₂₄ chain, in which the RPA correlation energy is highly nonlocal. In this model system, a softened Coulomb interaction is used to describe the electron-electron and electron-ion interactions, and slightly stronger nuclear charges (1.2e ) are assigned to the pseudo-H atoms. Bond stretching energies, RPA correlation potentials, and Kohn-Sham eigenvalues predicted by λ-ECDA are in good agreement with the benchmarks when the clusters are made reasonably large. We find that the LW RPA correlation energy is critical for obtaining accurate prediction of the RPA correlation potential, even though the LW RPA correlation energy contributes to only a few percent of the total RPA correlation energy.     

The homogeneous gas phase H2?D2 metathesis at room temperature: Reaction induced by specific vibrational excitation

We have demonstrated that reactions for which substantial activation energies are needed can be induced to occur at room temperature via specific vibrational excitation. Indeed, the indications are that the atom-switching reactions for which E_a > 25 kcal take place with high probability only when the activation energy is localized in the vibrational mode. In this preliminary report on the utilization of the stimulated Raman effect to generate substantial populations in the critical vibrational states required for the homogeneous atom exchange between H₂ and D₂, we first summarized the historical development of the concept. The experimental arrangement is then described and the analytical results tabulated; the observed dependence on relative concentrations is semiquantitatively rationalized on the basis of a model proposed in 1964. Independent shock tube and molecular beam investigations were similarly accounted for. Attention is called to the discrepancy between the generally concordant experimental results and the ab initio quantum mechanical calculations of the potential energy surface for 4H atoms.     

Quantum chemical study of the reaction mechanism of ozone and methane with fluorine and chlorine atoms

Ab initio calculations of the potential energy surface for the F + O₃ and Cl + O₃ reactions have been performed using the G3 and G3MP2 methods, which optimize the geometry configuration of reactants, products, intermediates, and transition states. The results show that fluorine atoms react with ozone as violently as chlorine atoms. At the same time, we have studied the reaction mechanisms of F atoms and Cl atoms with methane. It is found that fluorine atoms prefer to react with methane and chlorine atoms with ozone when there is competition between ozone and methane. Therefore, we can reasonably explain why chlorine atoms play the main role of reactants depleting ozone, while the more active fluorine atoms deplete less ozone. © 2002 Wiley Periodicals, Inc.; DOI 10.1002/qua.10119     

Beiträge zur Chemie der Pyrrolpigmente, 42. Mitt. Kraftfeldrechnungen an Gallenfarbstoffen: Die Energiehyperfläche rubinoider Pigmente

The energy hypersurface of rubinoid bile pigments is calculated using a force field model previously described. Two cases are observed. I: The pigment is substituted symmetrically or unsymmetrically by apolar groups. This results in a very shallow energy valley containing several enantiomeric conformers of approximately equal energies. II: Substitution by polar groups, especially in position 8 and 12 of the rubin skeleton, (e.g.–CH₂–CH₂–COOH) is followed by a lock in bonding between such groups and the lactam ring functions. Thereby only two enantiomeric conformers, which are energetically stabilized, are possible. The barrier between these species amounts to about 40 kJ/mol. These results are compared with the experimental facts available so far. An analysis of the corresponding energy hypersurfaces of the diastereomeric forms of (Z,Z)-, (E,Z)-, (Z,E)- and (E,E)-configurations is given as well.

41. Mitt.:Falk H., Müller N., Mh. Chem.112, 791 (1981).     

Halogen and Chalcogen Bond Energies Evaluated Using Electron Density Properties

Halogen (X-bond) and chalcogen bond (Ch−bond) energies for 36 complexes have been obtained at the RI-MP2/def2−TZVP level of theory, involving the heavier halogen and chalcogen atoms (Br, I, Se, Te). We have explored the existence of linear relationships between the interaction energies and the local kinetic energy densities at the bond critical points that characterize the σ-hole interactions (both electronic G(r) and potential V(r) energy densities). Interestingly, we have found strong relationships for halogen and chalcogen bonding energies, especially for the V(r) energy density, thus allowing to estimate the interaction energy without computing the separate monomers. This is also useful to estimate the interaction in monomeric systems (intramolecular X/Ch-bonds), as illustrated using several examples. Remarkably, we have also found a good relationship when in the same representation both halogen and chalcogen atoms are included, thus allowing to use the same empirical correlation for both interactions.     

The chemisorption of hydrogen on Cu(111): A dynamical study

Ab initio band-structure calculations within a density functional formalism were performed to compute the binding energy curves of atomic hydrogen with the high-symmetry adsorption sites of the (111) surface of copper. For a two-layer slab of Cu atoms and H coverage equal to 0.25, the binding energies are 2.25, 3.12, and 3.24 eV, for on-top, bridge, and threefold sites, so that the chemisorption of H₂ on Cu(111) is exothermic for threefold and bridge sites, but endothermic for on-top sites. Starting from these results, an LEPS potential for the interaction of H₂ with the Cu(111) surface was built. In this model potential, the most favored approaches correspond to a H₂ molecule parallel to the Cu surface, and for them, the activation barrier is located at the corner between the entrance and the exit channels of the reaction, and its lowest value is 0.6 eV. The LEPS potential was used in quasi-classical trajectories calculations to simulate the adsorption of a beam of H₂ molecules on Cu(111). The results show that (a) when H₂ is in the ground vibrational state the dissociative adsorption probability P_a increases from 0 to .90 along a roughly sigmoidal curve by increasing the collision kinetic energy from 0.4 to 1.3 eV, and (b) the vibrational energy can be as effective as the translational one in promoting dissociative chemisorption, in agreement with the experimental results. © 1994 John Wiley & Sons, Inc.     

Evaluation of the origin of rotational barrier in NH2X (X = NO, NS)

The origin of rotational barrier in N-thionitrosamine NH₂NS (TNA) has been examined with the aid of topological theory of atoms in molecules (AIM) and natural bond orbital (NBO) analyses and has been compared with N-nitrosamine NH₂NO (NA). Frequency calculations show that the rotational barrier for TNA is greater than NA. This can be attributed to the more charge transfer from nitrogen of amino group to sulfur in TNA than from nitrogen to oxygen in NA. NBO analysis reveals pyramidalization of nitrogen atom of NH₂ group leads to decrease of delocalization energy contribution and increase of Lewis energy contribution on total energy.     

Quantum-chemical calculation of a spillover model on a graphite support

The simplest quantum-chemical models of hydrogen spillover over a graphite-like surface as a proton or radical have been considered. The condensed planar C₂₄H₁₂ molecule was used as a model surface. Theab initio calculations of the interaction of hydrogen with the model surface were carried out by the restricted Hartree-Fock (HF) method in the STO-3G and 6-31 G^* basis sets. The radical hydrogen can not bind to such a surface, whereas the proton binds to it with an energy release of 186 kcal mol⁻¹. The activation energy of the transfer of the proton between two neighboring carbon atoms (10 kcal mol⁻¹) has been determined. The simplest model of the hydrogen migration as a proton over the model surface can be used for describing the spillover of hydrogen over the graphite surface. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 3, pp. 428–430, March, 1997.     

An accurate analytical formula for the van der Waals potentials of homonuclear rare-gas dimers with one adjustable parameter

In the present article, the Tang–Toennies–Yiu (TTY) potential model is modified by introducing one adjustable parameter. Then, the van der Waals potentials of He₂, Ne₂, Ar₂, Kr₂, and Xe₂ are calculated by this model with the adjustable parameter being determined by the well determined well depth D_e of these systems. Based on the derived potentials, the vibrational energy spacings of these systems are also calculated. It is shown that the present derived potentials and vibrational energy spacings agree well with experiment and other theoretical calculations. Finally, the normalization constant A in the asymptotic wave function of rare-gas atoms is estimated. The present derived normalization constant A is very close to the one by calculating the ratio between the Hartree–Fock function and the asymptotic wave function. The results confirm that absorbing the first-order polarization energy into the exchange energy expression is a well approximation for the present systems.     

An efficient newton-like method for molecular mechanics energy minimization of large molecules

Techniques from numerical analysis and crystallographic refinement have been combined to produce a variant of the Truncated Newton nonlinear optimization procedure. The new algorithm shows particular promise for potential energy minimization of large molecular systems. Usual implementations of Newton's method require storage space proportional to the number of atoms squared (i.e., O(N²)) and computer time of O(N³). Our suggested implementation of the Truncated Newton technique requires storage of less than O(N^1.5) and CPU time of less than O(N²) for structures containing several hundred to a few thousand atoms. The algorithm exhibits quadratic convergence near the minimum and is also very tolerant of poor initial structures. A comparison with existing optimization procedures is detailed for cyclohexane, arachidonic acid, and the small protein crambin. In particular, a structure for crambin (662 atoms) has been refined to an RMS gradient of 3.6 × 10^?6 kcal/mol/Å per atom on the MM2 potential energy surface. Several suggestions are made which may lead to further improvement of the new method.     

13-atom Ni-Al alloy clusters: Structures and dynamics

Structural and dynamical properties of model 13-atom Ni_nAl_m alloy clusters derived from a many-body potential are presented and discussed. Characterization of the structures corresponding to a given stoichiometric composition (i.e., chosen number of Ni and Al atoms) is carried out in terms of isomeric (geometric) forms and different distributions of the two types of atoms between the sites of a chosen isomer. We use the term homotops (“the same topography or geometry”) to label the structural forms that differ only by these distributions. The number and the energy spectra of the homotops are sensitive functions of the stoichiometric composition and isomeric form. Similarly to homogeneous clusters, alloy clusters undergo a solid-to-liquidlike transition as their energy is increased. Individual stages in the transition, such as isomerizations involving only surface atoms, isomerizations involving all atoms, surface melting (in a system as small as 13 atoms), and complete melting are identified and characterized. The actual occurrence of some or all of these stages in the meltinglike transition of a given cluster depends on the character of the energy spectra of its homotops, i.e., ultimately, on its stoichiometric composition. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 62: 185–197, 1997     

Computational studies of stable hexanuclear CulAgmAun (l + m + n = 6; l,m, n > 0) clusters

A DFT study was carried out on the ground state structures of ternary Cu_lAg_mAu_n (l + m + n = 6) clusters, with the aim of investigating changes of thermal and kinetic stabilities as an effect of composition, as well as the composition dependence of the electrostatic potential, of stable planar structures. DFT optimizations were performed using the PBE functional and the SDD basis set. All the optimized structures adopt planar geometries with bent triangular structures. Calculated binding energy values are in the range 1.5–1.9 eV/atom, which shows their thermal stability. The predicted HOMO‐LUMO energy gap values are in the semiconductor region, providing a qualitative indication of a moderate kinetic stability. NBO analyses indicate the existence of two mechanisms promoting planar structural stability, one due to bonding‐antibonding orbital interaction, and the other one due to the well‐known spd hybridization. Wiberg indices were obtained showing interatomic bonding. Electrostatic potential calculations show the existence of nucleophilic attack regions preferentially around silver and copper atoms located at the vertices while electrophilic attack regions are found in the vicinity of gold atoms over the cluster plane. Apparently, charge transfer occurs toward gold from silver and copper atoms when the concentration is favorable in the proximity of gold atoms. In particular, if the small ternary clusters discussed here contain only one gold atom, then a high electron density is observed at the site of this gold atom. © 2016 Wiley Periodicals, Inc.     

Extension of a simplified method for molecular correlation energy calculations to molecules containing third row atoms. I. Methodological developments

An extension of a simplified method for molecular correlation energy calculations to molecules containing third row atoms is presented. In addition to the use of pseudo-potentials in the calculations, the consequences of this extension on the different components of the energy partition which is the basic idea of the method, is analysed. Particular emphasis is placed on the specific role played by the 3d orbitals in each of the energy components. First, at the zeroth order, the energy is found to be very sensitive to the optimization of the 3d polarization functions. Secondly, the internal correlation energy, calculated by CI, requires the optimization of distinct 3d correlation orbitals to describe adequately the strong near-degeneracy effects that occur within the valence space. Finally it is shown that the 3d orbitals contribute partially to the non-internal correlation energy and that, the atoms-in-molecule structures corresponding typically to all-external contributions are negligible. The concept of error energy is introduced in place of the non-internal correlation energy: it includes the relativistic contributions within the semi-empirical tables. Such tables are presented for second row atoms and for the chlorine atom. From these tables, predicted values for some atomic term energies, experimentally undetermined, are derived. The methodological tests are limited here to the chlorine atom which is chosen for further applications in the next paper of this series. The conclusions concerning the applicability of the method to third row atoms are however quite general.Boursier I.R.S.I.A