D. Van Neck

The isoelectronic series of Be, Ne and Si are investigated using a variational determination of the second-order density matrix. A semidefinite program was developed that exploits all rotational and spin symmetries in the atomic system. We find that the method is capable of describing the strong static electron correlations due to the incipient degeneracy in the hydrogenic spectrum for increasing central charge. Apart from the ground-state energy various other properties are extracted from the variationally determined second-order density matrix. The ionization energy is constructed using the extended Koopmans' theorem. The natural occupations are also studied, as well as the correlated Hartree-Fock-like single particle energies. The exploitation of symmetry allows to study the basis set dependence and results are presented for correlation-consistent polarized valence double, triple and quadruple zeta basis sets.

In an earlier work, the authors developed a new method, the mobile block Hessian (MBH) approach, to accurately calculate vibrational modes for partially optimized molecular structures [ J. Chem. Phys. 2007, 126 (22), 224102.]. It is based on the introduction of blocks, consisting of groups of atoms, that can move as rigid bodies. The internal geometry of the blocks need not correspond to an overall optimization state of the total molecular structure. The standard MBH approach considers free blocks with six degrees of freedom. In the extended MBH approach introduced herein, the blocks can be connected by one or two adjoining atoms, which further reduces the number of degrees of freedom. The new approach paves the way for the normal-mode analysis of biomolecules such as proteins. It rests on the hypothesis that low-frequency modes of proteins can be described as pure rigid-body motions of blocks of consecutive amino acid residues. The method is validated for a series of small molecules and further applied to alanine dipeptide as a prototype to describe vibrational interactions between two peptide units; to crambin, a small protein with 46 amino acid residues; and to ICE/caspase-1, which contains 518 amino acid residues.

In an earlier paper, the authors have developed a new method, the mobile block Hessian (MBH), to accurately calculate vibrational modes for partially optimized molecular structures [J. Chem. Phys. 2007, 126 (22), 224102]. The proposed procedure remedies the artifact of imaginary frequencies, occurring in standard frequency calculations, when parts of the molecular system are optimized at different levels of theory. Frequencies are an essential ingredient in predicting reaction rate coefficients due to their input in the vibrational partition functions. The question arises whether the MBH method is able to describe the chemical reaction kinetics in an accurate way in large molecular systems where a full quantum chemical treatment at a reasonably high level of theory is unfeasible due to computational constraints. In this work, such a validation is tested in depth. The MBH method opens a lot of perspectives in predicting accurate kinetic parameters in chemical reactions where the standard full Hessian procedure fails.