Structural transition on cooling of plastic adamantane: A molecular-dynamics study

Density Functional Computations and Molecular Dynamics Simulations of the Triethylammonium Triflate Protic Ionic Liquid

Systematic molecular dynamics simulations based on an empirical force field have been carried out for samples of triethylammonium trifluoromethanesulfonate (triethylammonium triflate, [TEA][Tf]), covering a wide temperature range 200 K ≤ T ≤ 400 K and analyzing a broad set of properties, from self-diffusion and electrical conductivity to rotational relaxation and hydrogen-bond dynamics. The study is motivated by recent quasi-elastic neutron scattering and differential scanning calorimetry measurements on the same system, revealing two successive first order transitions at T ≈ 230 and 310 K (on heating), as well as an intriguing and partly unexplained variety of subdiffusive motions of the acidic proton. Simulations show a weakly discontinuous transition at T = 310 K and highlight an anomaly at T = 260 K in the rotational relaxation of ions that we identify with the simulation analogue of the experimental transition at T = 230 K. Thus, simulations help identifying the nature of the experimental transitions, confirming that the highest temperature one corresponds to melting, while the one taking place at lower T is a transition from the crystal, stable at T ≤ 260 K, to a plastic phase (260 ≤ T ≤ 310 K), in which molecules are able to rotate without diffusing. Rotations, in particular, account for the subdiffusive motion seen at intermediate T both in the experiments and in the simulation. The structure, distribution, and strength of hydrogen bonds are investigated by molecular dynamics and by density functional computations. Clustering of ions of the same sign and the effect of contamination by water at 1% wgt concentration are discussed as well.

Self-assembly of diamondoid molecules and derivatives (MD simulations and DFT calculations)

We report self-assembly and phase transition behavior of lower diamondoid molecules and their primary derivatives using molecular dynamics (MD) simulation and density functional theory (DFT) calculations. Two lower diamondoids (adamantane and diamantane), three adamantane derivatives (amantadine, memantine and rimantadine) and two artificial molecules (ADM•Na and DIM•Na) are studied separately in 125-molecule simulation systems. We performed DFT calculations to optimize their molecular geometries and obtained atomic electronic charges for the corresponding MD simulation, by which we predicted self-assembly structures and simulation trajectories for the seven different diamondoids and derivatives. Our radial distribution function and structure factor studies showed clear phase transitions and self-assemblies for the seven diamondoids and derivatives.

Orientational Melting and Reorientational Motion in a Cubane Molecular Crystal: A Molecular Simulation Study

Detailed molecular simulations are carried out to investigate the effect of temperature on orientational order in cubane molecular crystal. We report a transition from an orientationally ordered to an orientationally disordered plastic crystalline phase in the temperature range 425-450 K. This is similar to the experimentally reported transition at 395 K. The nature of this transition is first order and is associated with a 4.8% increase in unit cell volume that is comparable to the experimentally reported unit cell volume change of 5.4% (Phys. Rev. Lett. 1997, 78, 4938). An orientational order parameter, eta(T), has been defined in terms of average angle of libration of a molecular 3-fold axis and the orientational melting has been characterized by using eta(T). The orientational melting is associated with an anomaly in specific heat at constant pressure (C(P)) and compressibility (kappa). The enthalpy of transition and entropy of transition associated with this orientational melting are 20.8 J mol(-1) and 0.046 J mol(-1) K(-1), respectively. The structure of crystalline as well as plastic crystalline phases is characterized by using various radial distribution functions and orientational distribution functions. The coefficient of thermal expansion of the plastic crystalline phase is more than twice that of the crystalline phase.

High-pressure study of adamantane: variable shape simulations up to 26 GPa

We report simulations of adamantane by carefully combining ab initio and empirical approaches to enable simulations with internal degrees of freedom on crystalline adamantane up to a pressure of 26 GPa. Two sets of simulations, assuming the adamantane molecule as a rigid (RB) and flexible body (FB), have been carried out as a function of pressure up to 26 GPa to understand changes in the crystal structure as well as molecular structure. The flexible body simulations have been performed by including 6 lowest frequency internal modes (obtained from DFT calculations performed with Gaussian98) out of the total of 72. The calculated variation in c/a and V/V(0) from the RB and FB calculations as a function of pressure have been compared with the experimental curve. Other relevant thermodynamic and structural properties reported are the radial distribution functions, deviation in the position of a given type of atom with respect to its position at standard pressure, delta(s), cell parameters, volume, and energy. With an increase in pressure, three additional peaks are seen to develop gradually at three different pressures in the center of mass (com)-com radial distribution function (rdf). We attribute these changes to structural transformations (probably second-order phase transitions) which is consistent with the three phase transitions reported by Vijayakumar et al. for adamantane in the pressure range of 1 atm-15 GPa. Our simulations also show that these additional peaks in the rdf's are associated with the differences between opposite and parallel spin neighbors of Greig and Pawley as well as the crystallographic directional dependence of intermolecular distances in the first three shells of the neighbors. Also, the structural quantities from the RB calculation show considerable deviation from the FB calculation for pressures greater than 5 GPa, which suggests that the rigid body assumption for molecules may not be valid above this pressure.

Pressure-Induced Ordering in Adamantane: A Monte Carlo Simulation Study

Isothermal-isobaric ensemble Monte Carlo simulation of adamantane has been carried out with a variable shape simulation cell. The low-temperature crystalline phase and the room-temperature plastic crystalline phases have been studied employing the modified Williams potential. We show that at room temperature, the plastic crystalline phase transforms to the crystalline phase on increase in pressure. Further, we show that this is the same phase as the low-temperature ordered tetragonal phase of adamantane. The high-pressure ordered phase appears to be characterized by a slightly larger shift of the first peak toward a lower value of r in C-C, C-H, and H-H radial distribution functions as compared to the low-temperature tetragonal phase. The coexistence curve between the crystalline and plastic crystalline phase has been obtained approximately up to a pressure of 4 GPa.