A molecular mechanics implementation of the cyclic cluster model.

The implementation of the cyclic cluster model (CCM) for molecular mechanics is presented in the framework of the computational chemistry program deMon2k. Because the CCM is particularly well-suited for the description of periodic systems with defects, it can be used for periodic QM/MM approaches where the non-periodic QM part is treated as a defect in a periodic MM surrounding. To this end, we present here the explicit formulae for the evaluation of the Ewald sum and its first- and second-order derivatives as implemented in deMon2k. The outlined implementation was tested in molecular dynamics (MD) simulations and periodic structure optimization calculations. MD simulations of an argon system were carried out using the Nosé-Hoover chain (NHC) thermostat and the Martyna-Tobias-Klein (MTK) barostat to control the temperature and pressure of the system, respectively. For the validation of CCM structure optimization a set of molecular crystals were optimized using the Ewald method for the evaluation of the electrostatic interactions. Two optimization procedures for the determination of the atomic positions and CCM cell parameters were tested. Our results show that the simultaneous optimization of the atomic positions and cell parameters is most efficient. [ABSTRACT FROM AUTHOR]

Copyright of Zeitschrift für Naturforschung B: A Journal of Chemical Sciences is the property of De Gruyter and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)