A six-site intermolecular potential scheme for the azabenzene molecules, derived by crystal structure analysis

Computational Procedure of Lattice Energy Using the Ab Initio MO Method

Recently, we have proposed a computational procedure for calculations of lattice energies of molecular crystals using the ab initio MO method. This procedure does not use potential functions and is applicable to a variety of molecular crystals. The procedure has been successfully applied to calculation of packing structure of electron donor-acceptor complex, H3N–BF3, and hydrogen bonding crystal, CH3OH. In this work, we present a full account of the computational procedure. This method is applied to the packing structure calculations of hydrocarbon crystals, C2H2, C2H4 and C6H6. The lattice parameters optimized at the MP2/6-311++G** level are in good agreement with the experimental values. The basis set dependence of the lattice constants is also discussed for several crystals.

Anisotropic, Polarizable Molecular Mechanics Studies of Inter- and Intramolecular Interactions and Ligand-Macromolecule Complexes. A Bottom-Up Strategy

We present an overview of the SIBFA polarizable molecular mechanics procedure, which is formulated and calibrated on the basis of quantum chemistry (QC). It embodies nonclassical effects such as electrostatic penetration, exchange-polarization, and charge transfer. We address the issues of anisotropy, nonadditivity, and transferability by performing parallel QC computations on multimolecular complexes. These encompass multiply H-bonded complexes and polycoordinated complexes of divalent cations. Recent applications to the docking of inhibitors to Zn-metalloproteins are presented next, namely metallo-beta-lactamase, phosphomannoisomerase, and the nucleocapsid of the HIV-1 retrovirus. Finally, toward third-generation intermolecular potentials based on density fitting, we present the development of a novel methodology, the Gaussian electrostatic model (GEM), which relies on ab initio-derived fragment electron densities to compute the components of the total interaction energy. As GEM offers the possibility of a continuous electrostatic model going from distributed multipoles to densities, it allows an inclusion of short-range quantum effects in the molecular mechanics energies. The perspectives of an integrated SIBFA/GEM/QM procedure are discussed.

Vapor−Liquid and Vapor−Solid Phase Equilibria for United-Atom Benzene Models near Their Triple Points: The Importance of Quadrupolar Interactions

Gibbs ensemble Monte Carlo simulations were used to calculate the vapor-liquid and vapor-solid coexistence curves for benzene using two simple united-atom models. An extension of the Gibbs ensemble method that makes use of an elongated box containing a slab of the condensed phase with a vapor phase along one axis was employed for the simulations of the vapor-solid equilibria and the vapor-liquid equilibria at very low reduced temperatures. Configurational-bias and aggregation-volume-bias Monte Carlo techniques were applied to improve the sampling of particle transfers between the two simulation boxes and between the vapor and condensed-phase regions of the elongated box. An isotropic united-atom representation with six Lennard-Jones sites at the positions of the carbon atoms was used for both force fields, but one model contained three additional out-of-plane partial charge sites to explicitly represent benzene's quadrupolar interactions. Both models were fitted to reproduce the critical temperature and density of benzene and yield a fair representation of the vapor-liquid coexistence curve. In contrast, differences between the models are very large for the vapor-solid coexistence curve. In particular, the lack of explicit quadrupolar interactions for the 6-site model greatly reduces the energetic differences between liquid and solid phases, and this model yields a triple point temperature that is about a factor of 2 too low. In contrast, the 9-site model predicts a triple point of benzene at T = 253 +/- 6 K and p = 2.3 +/- 0.8 kPa in satisfactory agreement with the experimental data (T = 278.7 K and p = 4.785 kPa).