Interactions in ionic molecular crystals
Creators: name: Benedek, Nicole
Publisher: name: RMIT University
Year: 2024

description: We have used ab initio computational simulation techniques to investigate both intra- and intermolecular interactions in a novel family of ionic organophosphonate molecular crystals. We have examined the influence of various numerical approximations on the computed geometry and binding energies of a selection of well-characterised hydrogen bonded systems. It was found that numerical basis sets provided the efficiency required to study the large hydrogen bonded dimer anions present in the organophosphonate system, while also producing accurate geometries and binding energies. We then calculated the relaxed structures and binding energies of phenylphosphonic acid dimer in the two arrangements in which it is present in the bulk crystal. The computed geometries were in excellent agreement with the experimental structures and the binding energies were consistent with those found for other ionic hydrogen bonded systems. Electron density maps were used to gain insight into the nature of the hydrogen bonding interaction between phenylphosphonic acid dimers. We also examined the effect of aromatic ring substituents on the geometry and energetics of the hydrogen bonding interaction. The nitro-substituted dimer was predicted to have a stronger binding energy than its unsubstituted parent while the methyl-substituted dimer was predicted to have a similar binding energy to its unsubstituted parent. An analysis of crystal field effects showed that the structure of the phenylphosphonic acid dimers in the organophosphonates is a complex product of competing intra- and intermolecular forces and crystal field effects. Cooperative effects in the organophosphonate system were also investigated and it was found that the interactions were mostly one-body (local) in nature. We have examined the intramolecular charge-transfer interaction between copper-halogen cations in the organophosphonate materials. The origin of geometric differences between the Cu(I) starting material and Cu(II) product cations was attributed to the electronic configuration of the Cu ion, not crystal field effects. To gain further insight into the difference in electronic structure between the starting material and product, we attempted to simulate the step-by-step dissociation of the [CuI]+ system. Although this investigation was not successful, we were able to expose some of the pitfalls of simulating dissociating odd-electron systems. We also analysed and compared the charge-transfer interaction in the chloro-, bromo- and iodo-forms of the organophosphonate family. The charge-transfer interaction was predicted to increase on going from the chloro- to the iodo-form, consistent with solid-state UV-visible data. Finally, we used the highly accurate Quantum Monte Carlo (QMC) method to investigate the hydrogen bonding interaction in water dimer and to calculate the dissociation energy. The accuracy of the experimental estimate for the dissociation energy has recently been questioned and an alternative value has been put forward. Our results lend support to the validity of the alternative value and are also in excellent agreement with those from other high-level calculations. Our results also indicate that QMC techniques are a promising alternative to traditional wavefunction techniques in situations where both high accuracy and efficiency are important.