The search for stationary points on a quantum mechanical/molecular mechanical potential-energy surface

 We propose a methodology to locate stationary points on a quantum mechanical/molecular mechanical potential-energy surface. This algorithm is based on a suitable approximation of an initial full Hessian matrix, either a modified Broyden–Fletcher–Goldfarg–Shanno or a Powell update formula for the location of, respectively, a minimum or a transition state, and the so-called rational function optimization. The latter avoids the Hessian matrix inversion required by a quasi-Newton–Raphson method. Some examples are presented and analyzed. Received: 16 July 2001 / Accepted: 9 October 2001 / Published online: 9 January 2002     

An overlap criterion for selection of core orbitals

 An overlap criterion is defined that connects the identification of core orbitals in a molecular system, which can be problematic, to that in isolated atoms, which is well defined. This approach has been tested on a variety of troublesome systems that have been identified in the literature, including molecules containing third-row main-group elements, and is shown to remove errors of up to 100 kcal/mol arising from an inconsistent treatment of core orbitals at different locations on a potential-energy surface. For some systems and choices of core orbitals, errors as large as 19 kcal/mol can be introduced even when consistent sets of orbitals are frozen, and the new method is shown to identify these cases of substantial core–valence mixing. Finally, even when there is limited core–valence mixing, the frozen-core approximation can introduce errors of more than 5 kcal/mol, which is much larger than the presumed accuracy of models such as G2 and CBS-QB3. The source of these errors includes interatomic core–core and core–valence dispersion forces. Received: 31 August 2001 / Accepted: 11 October 2001 / Published online: 9 January 2002     

Dual-level direct dynamics calculations of deuterium kinetic isotope effects for the Cl(2P)+C2H6 abstraction reaction

 Kinetic isotope effects, KIEs, for hydrogen abstraction from C₂H₆ and C₂D₆ by chlorine atom have been studied by the dual-level direct dynamics approach. A low-level potential energy surface is obtained with the MNDO-SRP method. High-level structural properties of the reactants, transition state, and products were obtained at the MP2 level with the cc-pVDZ, aug-cc-pVDZ, and the cc-pVTZ basis sets. Using the variational transition state theory with microcanonical optimized multidimensional tunneling, the values of deuterium KIE, at 300 K, range from 2.28 to 3.27, in good agreement with the experimental values (2.69–5.88). Received: 6 June 2001 / Accepted: 12 July 2001 / Published online: 19 November 2001     

Modeling the H5+ potential-energy surface: a first attempt

 Ab initio molecular electronic structure calculations are performed for H₅ ⁺ at the QCISD(T) level of theory, using a correlation-consistent quadruple-zeta basis set. Structures, vibrational frequencies and thermochemical properties are evaluated for ten stationary points of the H₅ ⁺ hypersurface and are compared with previous calculations. The features of the H₃ ⁺…H₂ interaction at intermediate and large intermolecular distances are also investigated. Furthermore, an analytical functional form for the potential-energy surface of H₅ ⁺ is derived using a first-order diatomics-in-molecule perturbation theory approach. Its topology is found to be qualitatively correct for the short-range interaction region. Received: 15 March 2001 / Accepted: 5 July 2001 / Published online: 11 October 2001     

Reaction path following by quadratic steepest descent

An efficient steepest descent algorithm for the integration of minimum energy paths, based on local quadratic approximations of the potential energy surface, is presented. The algorithm incorporates a selection procedure for the points at which the second derivatives of the energy are calculated fully or partially, thus minimizing the computational effort while maintaining high accuracy. This makes the method especially well suited for application in variational transition state theory calculations with tunnelling corrections, which have very high accuracy requirements. The performance of the algorithm is illustrated by ab initio calculations for four chemical reactions of differing complexity. The overall computational cost is less than for, or comparable to that of, first- or second-order algorithms published previously. Received: 28 July 1998 / Accepted: 10 August 1998 / Published online: 28 October 1998     

A theoretical study on the ionization of OCS with an analysis of vibrational structures of the photoelectron spectrum

 Ab initio calculations have been performed to study the molecular structures and vibrational levels of the four low-lying ionic states (1, 2²Π, and 1, 2²Σ⁺) of carbonyl sulfide. The global regions of the potential-energy surfaces have been obtained by multireference single and double excitation configuration interaction calculations. Vibrational calculations using explicit vibrational Hamiltonians have been used for vibrational analysis. The equilibrium molecular structures and a vibrational analysis of the four ionic states are presented. The theoretical ionization intensity curves including the vibrational structures of the ionic states are also presented and are compared with the photoelectron spectrum. Received: 20 January 2001 / Accepted: 22 August 2001 / Published online: 30 October 2001     

In situ Fourier transform IR spectroscopy and variable-temperature wide-angle X-ray diffraction studies on the crystalline transformation of melt-crystallized nylon 12 12

 The crystalline transformation in nylon 12 12 was monitored by variable-temperature wide-angle X-ray diffraction during heating, isothermal and cooling processes. It could be found that the α phase of the sample transforms into a γ phase at about 130 °C if the sample is heated from room temperature to a high temperature, which is the so-called Brill transition temperature. In addition, nylon 12 12 was found to crystallize into the γ phase when isothermally crystallized at high temperature from the melt. Upon further cooling from the crystallization temperature to room temperature, the reverse transition from the γ phase to the α phase can be detected. Furthermore, in situ Fourier transform IR spectroscopy was also used to study the Brill transition of nylon 12 12 samples on both heating and cooling. It is interesting that the hydrogen-bond strength in nylons decreases dramatically and that some Brill bands of nylon 12 12 disappear abruptly during the Brill transition on heating. Received: 3 April 2001 Accepted: 28 July 2001     

Energetics, structures, bonding, and kinetics of the HSCl–HClS system

The [H,S,Cl] potential-energy surface has been investigated at the self-consistent field (SCF), complete active space self-consistent field (CASSCF), second-order M?ller–Plesset, coupled-cluster single-double and perturbative triple excitation, [CCSD(T)]/6-31G(d,p), 6-31G(2df,2pd), and correlation-consistent polarized valence triple zeta (cc-pVTZ) levels of theory. CCSD(T)/ cc-pVTZ results predict a very stable HSCl species, an isomer HClS, 51.84 kcal/mol higher in energy, and a transition state 57.68 kcal/mol above HSCl. Independent of the level of theory, results with the smaller 6-31G(d,p) basis set turned out to be poor, especially for HClS. Vibrational analysis indicates that both species can be easily differentiated if isolated. Bonding differences between these molecules are illustrated by contour plots of valence orbitals. Viewed classically, bonding in HClS involves a dative bond. Transition-state rate constants, and equilibrium constants for the HSCl ↔ HClS isomerization have been estimated for various temperatures (200–1000 K). At 298.15 K, the forward rate is predicted to be 7.95 × 10⁻²⁹ s⁻¹, and the equilibrium constant to be 2.31 × 10⁻³⁸. Tunneling corrections vary from 1.57 at 298.15 K to 1.05 at 1000 K. Activation energies have been obtained by a two-points linear fit to the Arrhenius equation. Received: 7 May 1999 / Accepted: 22 July 1999 / Published online: 4 October 1999     

The surface tension of aqueous polyvinylamine and copolymers with <Emphasis Type="Italic">N</Emphasis>-vinylformamide

 The surface tensions of aqueous poly(N-vinylformamide) (PNVF), polyvinylamine (PVAM), and PNVF–PVAM copolymers were measured as functions of pH. The nonionic PNVF gave a pH-independent surface tension of 68 mN/m. The surface tension of PVAM was pH dependent; at pH 10 it was 56 mN/m, whereas it was 71.5 mN/m at pH 3.5. The transition from higher to lower surface tension values occurred most dramatically between pH 8 and 9, reflecting the dissociation behavior of the amine groups. The copolymers showed intermediate behavior. Received: 20 August 2001 Accepted: 26 September 2001     

Direct ab initio dynamics calculations on the rate constants for the hydrogen-abstraction reaction of C2H5F with O (3P)

The hydrogen-abstraction reaction C₂H₅F+O → C₂H₄F+OH has been studied by a dual-level direct dynamics method. For the reaction, three reaction channels, one for α-abstraction and two for β-abstraction, have been identified. The potential-energy surface information is obtained at the MP2(full)/6-311G(d,p) and PMP2(full)/6-311G(3df,3pd) (single-point) levels. By canonical variational transition-state theory, rate constants for each reaction channel are calculated with a small-curvature tunneling correction. The total rate constant is calculated from the sum of the individual rate constants and the temperature dependence of the branching ratios is obtained over a wide range of temperatures from 300 to 5,000 K. The agreement of the rate constants with experiment is good in the experimental temperature range from 1,000 to 1,250 K. The calculated results indicate that at low temperatures α-abstraction is most likely to be the major reaction channel, while β-abstraction channels will significantly contribute to the whole reaction rate as the temperature increases. Received: 23 January 2002 / Accepted: 23 June 2002 / Published online: 20 September 2002     

Transition state for intramolecular C–H bond cleavage in [(LCu)2(μ-O)2]2+ (L = 1,4,7-tribenzyl-1,4,7-triazacyclononane)

Hybrid quantum mechanical/molecular mechanical electronic structure calculations reveal the transition state for C–H bond cleavage in [(LCu)₂ (μ-O)₂]²⁺ (L=1,4,7-tribenzyl-1,4,7-triazacyclononane) to be consistent with a hydrogen-atom-transfer mechanism from carbon to oxygen. At the MPW1K/double-zeta effective core potential(+)|univeral force field level, 0 K activation enthalpies for the parent, p-CF3, and p-OH substituted benzyl systems are predicted to be 8.8, 9.5, and 7.8 kcal/mol. Using a one-dimensional Eckart potential to estimate quantum effects on the reaction coordinate, reaction in the unsubstituted system is predicted to proceed with a primary kinetic isotope effect of 22 at 233 K. Structural parameters associated with the hydrogen-atom transfer are consistent with the Hammond postulate. Received: 10 October 2000 / Accepted: 3 November 2000 / Published online: 3 April 2001     

Recombination of silicon ions by electron capture from atomic hydrogen and helium

 Ab initio potential-energy curves and coupling matrix elements of the Σ and Π molecular states involved in the collision of the Si²⁺, Si³⁺ and Si⁴⁺ multicharged ions on atomic hydrogen and helium have been determined by means of configuration interaction methods. The total and partial electron capture cross sections have been determined using a semiclassical or a quantal approach in the 0.002–0.1 au velocity range. A detailed comparison with very recent theoretical and experimental rate coefficient results is made. Received: 14 September 1999 / Accepted: 3 February 2000 / Published online: 2 May 2000     

RISM-SCF study of the free-energy profile of the Menshutkin-type reaction NH3+CH3Cl→NH3CH3++Cl− in aqueous solution

The free-energy profile for the Menshutkin-type reaction NH₃ + CH₃Cl → NH₃CH₃ ⁺ + Cl⁻ in aqueous solution is studied using the RISM-SCF method. The effect of electron correlation on the free-energy profile is estimated by the RISM-MP2 method at the HF optimized geometries along the reaction coordinate. Solvation was found to have a large influence on the vibrational frequencies at the reactant, transition state and product; these vibrational frequencies are utilized to calculate the zero-point energy correction of the free-energy profile. The computed barrier height and reaction exothermicity are in reasonable agreement with those of experiment and previous calculations. The change of solvation structure along the reaction path is represented by radial distribution functions between solute-solvent atomic sites. The mechanisms of the reaction are discussed from the view points of solute electronic and solvation structures. Received: 26 June 1998/Accepted: 28 August 1998 / Published online: 2 November 1998     

On the automatic restricted-step rational-function-optimization method

The rational function optimization algorithm is one of the widely used methods to search stationary points on surfaces. However, one of the drawbacks of this method is the step reduction procedure to deal with the overstepping problem. We present and comment on a method such that the step obtained from the solution of the rational function equations possesses the desired correct length. The analysis and discussion of the method is mainly centered on the location and optimization of transition states. Received: 18 June 1998 / Accepted: 17 September 1998 / Published online: 23 November 1998     

A fast algorithm to compute atomic charges based on the topology of the electron density

This work proposes a novel algorithm to compute atomic charges as defined by the theory of “atoms in molecules” (AIM). Using the divergence theorem it is possible to express the 3D volume integral over an atomic basin purely in terms of 2D surface integrals. Hence, it can be proven that an atomic charge is equal to the flux of the electric field of the whole molecule through the atom's complete boundary. This boundary consists of the interatomic surfaces and the so-called outeratomic surface, which is the open side of the atom. When fine-tuned the algorithm can generate atomic charges in the order of minutes without introducing any approximations. Moreover, the problem of the geometrical cusp occurring in atomic basins and that of multiple intersections is also eliminated. The computational overhead of computing the electric field (which is analytical) is compensated by the gain in computing time by eliminating one dimension of quadrature. The proposed algorithm opens an avenue to invalidate the oft-quoted drawback that AIM charges are computationally expensive. We explain the details of the implementation in MORPHY01 and illustrate the novel algorithm with a few examples. Received: 1 June 2000 / Accepted: 4 October 2000 / Published online: 23 January 2001     

Theoretical study of the potential-energy surface of C2NP

 The B3LYP/6-311G(d) and CCSD(T)/6-311G(2df) (single-point) methods have been used to investigate the singlet potential energy surface of C₂NP, in which seven stationary isomers and seventeen interconversion transition states are involved. At the final CCSD(T)/6-311G(2df)//B3LYP6-311G(d) level with zero-point vibrational energy correction the lowest-lying isomer is a linear NCCP followed by two linear CNCP isomers at 23.9  and CCNP at 65.8 kcal mol⁻¹, respectively. The three isomers are kinetically very stable towards both isomerization and dissociation, and CCNP is even more kinetically stable than CNCP – by 14.3 kcal mol⁻¹ despite its high energy. Further comparative calculations were performed at the QCISD and QCISD(T) levels with the 6-311G(d) and 6-311G(2d) basis sets to obtain more reliable structures and spectroscopy for the three isomers. The calculated bond lengths, rotational constant, and dipole moment for NCCP were in close agreement with the experimentally determined values. Finally, similarities and discrepancies between the potential energy surface of C₂NP and those of the analogous species C₂N₂ and C₂P₂ were compared. The results presented in this paper might be helpful for future identification of the two still unknown yet kinetically very stable isomers CNCP and CCNP, both in the laboratory and in interstellar space. Received: 3 January 2001 / Accepted: 6 June 2001 / Published online: 30 October 2001     

The occurrence of electron transfer in aromatic nitration: dynamical aspects

Electron transfer (ET) from toluene to the nitronium ion in the region of van der Waals intermolecular distances has been investigated by a quantum dynamical analysis performed on potential-energy surfaces computed at the ab initio multireference configuration interaction level. The results show that ET is very fast, occurring on a timescale of a few picoseconds. This has important implications for the mechanism of aromatic nitration: the ET path can compete efficiently with the direct attack of the nitronium ion to the aromatic substrate to yield the Wheland intermediate and that could explain some unsettled points in the mechanism of aromatic nitration. Received: 16 September 1999 / Accepted: 3 February 2000 / Published online: 12 May 2000     

Extracting parameters for base-pair level models of DNA from molecular dynamics simulations

 A method is described to extract a complete set of sequence-dependent energy parameters for a rigid base-pair model of DNA from molecular dynamics (MD) simulations. The method is properly consistent with equilibrium statistical mechanics and leads to effective inertia parameters for the base-pair units as well as stacking and stiffness parameters for the base-pair junctions. We give explicit formulas that yield a complete set of base-pair model parameters in terms of equilibrium averages that can be estimated from a time series generated in an MD simulation. The expressions to be averaged depend strongly both on the choice of coordinates used to describe rigid-body orientations and on the choice of strain measures at each junction. Received: 12 July 2000 / Accepted: 5 January 2001 / Published online: 3 May 2001     

Ab-initio study on low-lying states of the TiSi molecule

 The electronic structure of the TiSi molecule was examined using two types of multireference single and double excitation configuration interactions with highly extended basis sets, one including valence correlation and the other including valence and core–valence correlation. A multireference coupled-pair approximation (MRCPA) was further applied to the latter. The calculations suggest a ⁵Δ ground state, and the lowest excited state is ³Π and is only slightly (0.12 eV as estimated by MRCPA) above the ground state. The spectroscopic constants of the low-lying ¹Δ, ³Δ, ¹Π, ⁵Π, and ⁷Σ⁺ states as well as the ⁵Δ ground state and the ³Π excited states were evaluated, and we found that the molecule has only a weak σ bond and that six of the eight valence electrons essentially do not contribute to the bonding. The bonding nature of TiSi in these states is discussed in comparison with the TiC molecule. Received: 7 October 2000 / Accepted: 8 January 2001 / Published online: 3 May 2001