Atomistic simulations of co-crystal formation via mechanochemistry

点击翻译按钮获取中文摘要

Most chemical synthesis is performed in solution because in this phase it is easy to ensure that there are a large number of reactive collisions between reactant molecules. In addition, solution chemistry is well understood and we thus have a high degree of control over the reactions that can be performed and the products that can be synthesised. The problem with this approach is twofold. Firstly, the solvents many solvents are environmentally unfriendly and secondly separating the product from the solution at the end of the reaction often requires distillation, which requires a large input of energy and which introduces an extra step to the whole process. It would thus be enormously beneficial if this step could be avoided and if the solvent could be eliminated. Mechanochemical reactions allow for just this possibility. In these processes the reactants are powdered crystals. These powders are mixed together and mechanical work is done on the mixture in, for example, a mortar and pestle, a ball mill or an extruder. Experiments have demonstrated that it is possible to do a wide range of reactions in this way i.e., "mechanochemically". Furthermore, these mechanochemical processes are seen in some quarters to be the best way to synthesise systems known as co-crystals in which one or more chemical components are packed together into an ordered, crystalline structure. However, wider use of these processes and commercialization of these technologies is prevented because of the relative lack of understanding of the fundamental mechanisms that are in play in these reactions. The aim of this project is to examine what happens in a mechanochemical reaction by performing molecular dynamics simulations using a computer. Such simulations are useful because it is possible to keep track of the positions of all the atoms at all times. This, however, is also the difficulty as specialized tools are required to make sense of large volume of high dimensional data that emerges from such simulations. One of our intentions is, therefore, to develop computational tools for studying these highly complex processes. Throughout the work a reaction between two pharmaceutically active molecules, aspirin and meloxicam, will be studied. We will construct models for nanoparticles composed of each of these molecule types and will use non-equilibrium molecular dynamics simulations to force collisions between these particles to occur. Collisions will be performed for a range of collision velocities and for a number of different collision geometries. We will investigate head on collisions between the particles and glancing collisions as well as collisions in which we will change the relative orientations of the two crystal structures. For all these various kinds of collisions we will investigate the degree to which the two chemical components mix and the degree to which the crystallinity of the structure is disrupted by the collision. This work will give us one of the first visualizations of the zone of reaction in a mechanochemical process. More importantly, however, it will provide us with a way of rationalising what is being observed in the reactive zone. This work will thus provide new fundamental insights into how and why these reactions proceed and will serve as a basis for future work on the comercial exploitation of these reactions.