P.W. Ayers

A procedure for determining force constants that is independent of the internal redundant coordinate choice is presented. The procedure is based on solving each bond and angle term separately, using the Wilson B matrix. The method only requires a single ab initio frequency calculation at the minimum energy structure and is made available in the software "parafreq". The methodology is validated with a set of small molecules, by showing it can reproduce ab initio frequencies better than other methods such as taking the diagonal terms of the Hessian in internal coordinates or by using standard AMBER force fields. Finally, the utility of the method is demonstrated by parametrizing the dizinc scaffold of bis-dipicolylamine (BDPA) bound to phosphotyrosine, which is then functionalized into promising antitumor drug proteomimetics.

The Fukui function is considered as the diagonal element of the Fukui matrix in position space, where the Fukui matrix is the derivative of the one particle density matrix (1DM) with respect to the number of electrons. Diagonalization of the Fukui matrix, expressed in an orthogonal orbital basis, explains why regions in space with negative Fukui functions exist. Using a test set of molecules, electron correlation is found to have a remarkable effect on the eigenvalues of the Fukui matrix. The Fukui matrices at the independent electron model level are mathematically proven to always have an eigenvalue equal to exactly unity while the rest of the eigenvalues possibly differ from zero but sum to zero. The loss of idempotency of the 1DM at correlated levels of theory causes the loss of these properties. The influence of electron correlation is examined in detail and the frontier molecular orbital concept is extended to correlated levels of theory by defining it as the eigenvector of the Fukui matrix with the largest eigenvalue. The effect of degeneracy on the Fukui matrix is examined in detail, revealing that this is another way by which the unity eigenvalue and perfect pairing of eigenvalues can disappear.

For the Hirshfeld-I atom in the molecule (AIM) model, associated single-atom energies and interaction energies at the Hartree–Fock level are efficiently determined in one-electron Hilbert space. In contrast to most other approaches, the energy terms are fully consistent with the partitioning of the underlying one-electron density matrix (1DM). Starting from the Hirshfeld-I AIM model for the electron density, the molecular 1DM is partitioned with a previously introduced double-atom scheme (Vanfleteren et al., J Chem Phys 2010, 132, 164111). Single-atom density matrices are constructed from the atomic and bond contributions of the double-atom scheme. As the Hartree–Fock energy can be expressed solely in terms of the 1DM, the partitioning of the latter over the AIM naturally leads to a corresponding partitioning of the Hartree–Fock energy. When the size of the molecule or the molecular basis set does not grow too large, the method shows considerable computational advantages compared with other approaches that require cumbersome numerical integration of the molecular energy integrals weighted by atomic weight functions. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011

A novel force-field development strategy is proposed that tackles the well-known difficulty of parameter correlations arising in a conventional least-squares optimization. In the first step of the new gradient curves method (GCM), continuity criteria are imposed to transform the raw multidimensional ab initio training data to distinct sets of one-dimensional data, each associated with an individual energy term. In the second step, the transformed data suggest suitable analytical expressions, and the parameters in these expressions are fitted to the transformed data; that is, one does not have to postulate a priori analytical expressions for the force-field energy terms. This approach facilitates the derivation of valence terms. Benchmarks have been performed on a set of small molecules. The results show that the new method yields physically acceptable energy terms exactly when a conventional parametrization would suffer from parameter correlations, that is, when an increasing number of redundant internal coordinates is used in the force-field model. The generic treatment of parameter correlations in the proposed method facilitates an intuitive physical interpretation of the individual terms in the force-field expression, which is a prerequisite for the transferability of force-field models.