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.

Thermodynamic properties of ammonia clusters (NH(3))(n) n=2-11: Comparing classical and quantum simulation results for hydrogen bonded species

MELLA, MASSIMO
2009-01-01

Abstract

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.
2009
http://scitation.aip.org/content/aip/journal/jcp/131/3/10.1063/1.3159398
C., Lubombo; E., Curotto; P. E., Janeiro; Mella, Massimo
File in questo prodotto:
File Dimensione Formato  
AmmoC_thermoPIMC.pdf

accesso aperto

Descrizione: PDF editoriale
Tipologia: Altro materiale allegato
Licenza: Creative commons
Dimensione 592.93 kB
Formato Adobe PDF
592.93 kB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11383/1735840
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 19
  • ???jsp.display-item.citation.isi??? 19
social impact