Lattice Dynamics of Solid Cubane within the Quasi-Harmonic Approximation

Solid cubane, which is composed of weakly interacting cubic molecules, exhibits many unusual and interesting properties, such as a very large thermal expansion and a first-order phase transition at T$_{p}$=394 K from an orientationally-ordered phase of R$\bar{3}$ symmetry to a {\it non-cubic} disordered phase of the same symmetry with a volume expansion of 5.4%, among the largest ever observed. We study the lattice dynamics of solid cubane within the quasi-harmonic and rigid-molecule approximation to explain some of these unusual dynamical properties. The calculated phonon density of states, dispersion curves and thermal expansion agree surprisingly well with available experimental data. We find that the amplitude of thermally excited orientational excitations (i.e. librons) increases rapidly with increasing temperature and reaches about 35$^{\rm o}$ just before the orientational phase transition. Hence, the transition is driven by large-amplitude collective motions of the cubane molecules. Similarly the amplitude of the translational excitations shows a strong temperature dependence and reaches one tenth of the lattice constant at T=440 K. This temperature is in fair agreement with the experimental melting temperature of 405 K, indicating that the Lindemann criterion works well even for this unusual molecular solid.Comment: 15 pages, 6 figures (devoted to Prof. Ciraci in honor of his sixtieth birthday

From dimers to the solid-state: Distributed intermolecular force-fields for pyridine

A.A. thanks A.W.E. financial support through the EngDoc studentship from M3S Centre for Doctoral Training (EPSRC Grant No. EP/G036675/1). General computational infrastructure used is developed under No. EPSRC EP/K039229/1

Quasi Harmonic Lattice Dynamics and Molecular Dynamics calculations for the Lennard-Jones solids

We present Molecular Dynamics (MD), Quasi Harmonic Lattice Dynamics (QHLD) and Energy Minimization (EM) calculations for the crystal structure of Ne, Ar, Kr and Xe as a function of pressure and temperature. New Lennard-Jones (LJ) parameters are obtained for Ne, Kr and Xe to reproduce the experimental pressure dependence of the density. We employ a simple method which combines results of QHLD and MD calculations to achieve densities in good agreement with experiment from 0 K to melting. Melting is discussed in connection with intrinsic instability of the solid as given by the QHLD approximation. (See http://www.fci.unibo.it/~valle for related papers)Comment: 7 pages, 5 figures, REVte

When a proton attacks cellobiose in the gas phase: ab initio molecular dynamics simulations

Investigations of reaction pathways between a proton and cellobiose (CB), a glucose disaccharide of importance, were carried out in cis and trans CB using Ab Initio Molecular Dynamics (AIMD) simulations starting from optimized configurations where the proton is initially placed near groups with affinity for it. Near and above 300 K, protonated CB (H(+)CB) undergoes several transient reactions including charge transfer to the sugar backbone, water formation and dehydration, ring breaking and glycosidic bond breaking events as well as mutarotation and ring puckering events, all on a 10 ps timescale. cis H(+)CB is energetically favoured over trans H(+)CB in vacuo, with an energy gap larger than for the neutral CB