D. Van Neck

An efficient protocol is presented to identify signals in vibrational spectra of silica oligomers based on theoretical molecular dynamics (MD) simulations. The method is based on the projection of the atomic velocity vectors on the tangential directions of the trajectories belonging to a predefined set of internal coordinates. In this way only contributions of atomic motions along these internal coordinates are taken into consideration. The new methodology is applied to the spectra of oligomers and rings, which play an important role in zeolite synthesis. A suitable selection of the relevant internal coordinates makes the protocol very efficient but relies on intuition and theoretical insight. The simulation data necessary to compute vibrational spectra of relevant silica species are obtained through MD using proper force fields. The new methodology—the so-called velocity projection method—makes a detailed analysis of vibrational spectra possible by establishing a one-to-one correspondence between a spectral signal and a proper internal coordinate. It offers valuable perspectives in understanding the elementary steps in silica organization during zeolite nanogrowth. The so-called velocity projection method is generally applicable on data obtained from all types of MD and is a highly valuable alternative to normal-mode analysis which has its limitations due to the presence of many local minima on the potential energy surface. In this work the method is exclusively applied to inelastic neutron scattering, but extension to the infrared power spectrum is apparent.

Based on the so-called Hirshfeld atom in the molecule scheme, a new AIM method is presented. The method is similar to the Hirshfeld-I scheme, with the AIM weight function being constructed by minimizing the information loss upon formation of the molecule, but now requiring explicitly that the promolecular densities integrate to the same number of electrons as the AIM densities. This new weight function leads to a new iterative AIM scheme, and the resulting operative scheme is examined and discussed. The final results indicate that the newly proposed method does not perform as well as the Hirshfeld-I method.