Tsui, Helen Hoi Yan; (2001) New models for intermolecular repulsion and their application to van der Waals complexes and crystals of organic molecules. Doctoral thesis (Ph.D), UCL (University College London).

Preview

Text New_models_for_intermolecular_.pdf Download (17MB) | Preview

Abstract

Model inter molecular potentials are required for simulations of molecules in the gas, liquid, or solid phase. The widely used isotropic atom-atom model potentials are empirically fitted and based on the assumptions of transferability, combining rules and that atoms in molecules are spherical. This thesis develops a non-empirical method of modelling repulsion by applying the overlap model, which we show as a general non-empirical method of deriving repulsion potentials for a specific molecule. In this thesis, the repulsion parameters for an exponential atom-atom model potential are obtained from the ab initio charge density of a small organic molecule by making the assumption that the repulsion is proportional to the overlap of a pair of molecules. The proportionality constant is fixed by a limited number of intermolecular perturbation theory (IMPT) calculations. To complete the model potential, the electrostatic interaction is represented by a distributed multipole analysis, and the Slater-Kirkwood formula is used for the dispersion. These non-empirical potentials can reproduce experimental crystal structure when applied to crystal structure prediction of an oxyboryl derivative. A detailed study on further improving the overlap model was carried out for phenol-water, by including other minor intermolecular contributions of charge-transfer and penetration. High quality ab initio calculations on the complex were performed for use in comparison. To compare with experimental data, diffusion Monte Carlo simulations were performed with the potential, so that the effects of anharmonic zero-point motion on structure and energy of the system are included. When the system is too large for an IMPT calculation, the proportionality constant can be determined empirically by fitting the cell volume as shown in our study of crystal structures of chlorothalonil. This is used with an anisotropic repulsion model that has been derived for Cl and N atoms in chlorothalonil. This model potential reproduces the crystal structure well, and is used to predict other polymorphs by a search for minima in the lattice energy.

Download activity - last month

Export as JSONXMLCSV

Download activity - last 12 months

Export as XMLJSONCSV

Downloads by country - last 12 months

Export as JSONXMLCSV

Archive Staff Only

View Item