Structure and energetics of ammonia clusters (NH3)n (n=3-20) investigated using a rigid-polarizable model derived from ab initio calculations

An analytical model has been developed to describe the interaction between rigid ammonia molecules including the explicit description of induction. The parameters of the model potential were chosen by fitting high quality ab initio data obtained using second-order Moller-Plesset (MP2) perturbation theory and extended basis sets. The description of polarization effects is introduced by using a noniterative form of the "charge on spring model", the latter accounting for more than 95% of the dipole induction energy and of the increased molecular dipole. Putative global minima for (NH3)(n) (n = 3-20) have been optimized using this new model, the structure and energetics of the clusters with n = 3-5 being found in good agreement with previous ab initio results including electronic correlation. Results for larger species have been compared with previous structural studies where only nonpolarizable models were employed. Our model predicts larger binding energies for any cluster size than previous analytical surfaces, the results often suggesting a reorganization of the relative energy ranking and a different structure for the global minimum

Study of the structure, energetics, and vibrational properties of small ammonia clusters (NH3)n (n=2-5) using correlated ab initio methods

Equilibrium geometries, interaction energies, and harmonic frequencies of (NH3)(n) isomers (n = 2-5) have been computed using correlated calculations (MP2) in conjunction with Dunning's aug-cc-pVXZ (X = D, T, Q) basis sets and the Counterpoise procedure. Whenever available, literature values for the binding energy and geometry of dimers and trimers agree well with our data. Low lying isomers for (NH3)(4) and (NH3)(5) have been found to have similar binding energies (roughly 16 and 20 kcal/mol for the tetramer and pentamer, respectively), perhaps suggesting the presence of a very smooth energy landscape. Using BSSE corrected forces or freezing the monomer structure to its gas phase geometry have been found to have only a weak impact on the energetic and structural properties of the clusters. The effect of zero-point energy (ZPE) on the relative stability of the clusters has been estimated using harmonic frequencies. The latter also highlighted the presence of vibrational fingerprints for the presence of double acceptor ammonia molecules. Many-body effects for (NH3)(n) isomers (n = 2-4) have been investigated to explore the possibility of building a pairwise interaction model for ammonia. In the frame of the work presented, we have found the 3-body effect to account for 10-15% of the total interaction energy, whereas the 4-body effects may be neglected as first approximation