Interpolating moving least-squares methods for fitting potential energy surfaces: using classical trajectories to explore configuration space.

We develop two approaches for growing a fitted potential energy surface (PES) by the interpolating moving least-squares (IMLS) technique using classical trajectories. We illustrate both approaches by calculating nitrous acid (HONO) cis-->trans isomerization trajectories under the control of ab initio forces from low-level HF/cc-pVDZ electronic structure calculations. In this illustrative example, as few as 300 ab initio energy/gradient calculations are required to converge the isomerization rate constant at a fixed energy to approximately 10%. Neither approach requires any preliminary electronic structure calculations or initial approximate representation of the PES (beyond information required for trajectory initial conditions). Hessians are not required. Both approaches rely on the fitting error estimation properties of IMLS fits. The first approach, called IMLS-accelerated direct dynamics, propagates individual trajectories directly with no preliminary exploratory trajectories. The PES is grown "on the fly" with the computation of new ab initio data only when a fitting error estimate exceeds a prescribed tight tolerance. The second approach, called dynamics-driven IMLS fitting, uses relatively inexpensive exploratory trajectories to both determine and fit the dynamically accessible configuration space. Once exploratory trajectories no longer find configurations with fitting error estimates higher than the designated accuracy, the IMLS fit is considered to be complete and usable in classical trajectory calculations or other applications.

Electrochemical and theoretical investigation of corannulene reduction processes.

The voltammetric generation of corannulene anions was investigated over a large range of experimental conditions comprising either "traditional" electrochemical solvents, such as dimethylformamide, acetonitrile, and tetrahydrofuran, or "unconventional" solvents, such as liquid ammonia, liquid methylamine, or liquid dimethylamine, and several different supporting electrolytes. Strong ion pairing effects were found to dominate the electrochemical generation of corannulene higher anions, and through the suitable choice of the solvent/electrolyte system, we observed, for the first time, the reversible electrochemical generation of up to the triply reduced corannulene. The standard potentials obtained experimentally compared rather well with the theoretical values calculated by ab initio and density functional methods, in which solvation and ion pairing effect were explicitly taken into account. In particular, the calculations considered the effect of the electrolyte cation size on ion pairing in order to rationalize the occurrence of the third reduction within the experimental potential window.

Parity violation energy of biomolecules--I: polypeptides.

A computational algorithm is developed to compute the energy of parity-violation (E(pv)) due to weak-nuclear forces acting within chain polymers. The method is applied to estimate the magnitude of E(pv) associated with the folding of alpha-helices in polypeptide chains and, in turn, of typical proteins. Implications to the development of biological homochirality in nature are discussed in terms of a simple evolutionary model to assess the likelihood for energies of the order of the computed E(pv) to have determined the direction of natural homochirality.