Molecular dynamics simulations of p-tert-butylcalix[4]arene with small guest molecules

Classical molecular dynamics simulations were used to study p-tert-butylcalix[4]arene inclusion compounds with xenon, nitrogen, hydrogen, methane, and sulfur dioxide guest molecules. The calixarene units were taken to be rigid and the intermolecular molecular interactions were modeled as a sum of the van der Waals interactions with parameters from the AMBER force field and electrostatic interactions. Simulations of the high-density alpha phase and low-density beta0 phase of p-tert-butylcalix[4]arene were used to test the force field. The predicted densities of the two phases were found to agree with experimental measurements at 173 K to within 5 %. Simulations were performed with guests placed inside the calixarene cages of the beta0 phase. Guest-host ratios of 1:1 to 1:4 were considered. Changes in the unit-cell volume and density of the phases with the addition of guest molecules and the inclusion energies for the guests were determined. Finally, the dynamics of the guest motion inside the cages were characterized by determining the root-mean-square displacements and velocity autocorrelation functions of the xenon and nitrogen guests.