Normal modes for large molecules with arbitrary link constraints in the mobile block Hessian approach

In a previous paper [Ghysels et al., J. Chem. Phys. 126, 224102 (2007)] the mobile block Hessian (MBH) approach was presented. The method was designed to accurately compute vibrational modes of partially optimized molecular structures. The key concept was the introduction of several blocks of atoms, which can move as rigid bodies with respect to a local, fully optimized subsystem. The choice of the blocks was restricted in the sense that none of them could be connected, and also linear blocks were not taken into consideration. In this paper an extended version of the MBH method is presented that is generally applicable and allows blocks to be adjoined by one or two common atoms. This extension to all possible block partitions of the molecule provides a structural flexibility varying from very rigid to extremely relaxed. The general MBH method is very well suited to study selected normal modes of large macromolecules (such as proteins and polymers) because the number of degrees of freedom can be greatly reduced while still keeping the essential motions of the molecular system. The reduction in the number of degrees of freedom due to the block linkages is imposed here directly using a constraint method, in contrast to restraint methods where stiff harmonic couplings are introduced to restrain the relative motion of the blocks. The computational cost of this constraint method is less than that of an implementation using a restraint method. This is illustrated for the $α$-helix conformation of an alanine-20-polypeptide.