Thermodynamic properties of ammonia clusters (NH3)n n = 2-11: comparing classical and quantum simulation results for hydrogen bonded species.

Classical and quantum simulations of ammonia clusters in the dimer through the hendecamer range are performed using the stereographic projection path integral. Employing the most recent polarizable potential to describe intermolecular interactions, energetic and structural data obtained with our simulations provide support for a more fluxional or flexible nature at low temperature of the ammonia dimer, pentamer, and hexamer than in the other investigated species. The octamer and the hendecamer display a relatively strong melting peak in the classical heat capacity and a less intense but significant melting peak in the quantum heat capacity. The latter are shifted to lower temperature (roughly 15 and 40 K lower, respectively) by the quantum effects. The features present in both classical and quantum constant volume heat capacity are interpreted as an indication of melting even in the octamer case, where a large energy gap is present between its global minimum and second most stable species. We develop a first order finite difference algorithm to integrate the geodesic equations in the inertia ellipsoid generated by n rigid nonlinear bodies mapped with stereographic projections. We use the technique to optimize configurations and to explore the potential surface of the hendecamer.

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 Møller-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.