An approach to developing a force field for molecular simulation of martensitic phase transitions between phases with subtle differences in energy and structure

Uncovering molecular processes in crystal nucleation and growth by using molecular simulation

Exploring nucleation processes by molecular simulation provides a mechanistic understanding at the atomic level and also enables kinetic and thermodynamic quantities to be estimated. However, whilst the potential for modeling crystal nucleation and growth processes is immense, there are specific technical challenges to modeling. In general, rare events, such as nucleation cannot be simulated using a direct "brute force" molecular dynamics approach. The limited time and length scales that are accessible by conventional molecular dynamics simulations have inspired a number of advances to tackle problems that were considered outside the scope of molecular simulation. While general insights and features could be explored from efficient generic models, new methods paved the way to realistic crystal nucleation scenarios. The association of single ions in solvent environments, the mechanisms of motif formation, ripening reactions, and the self-organization of nanocrystals can now be investigated at the molecular level. The analysis of interactions with growth-controlling additives gives a new understanding of functionalized nanocrystals and the precipitation of composite materials.

Modelling organic crystal structures using distributed multipole and polarizability-based model intermolecular potentials

Crystal structure prediction for organic molecules requires both the fast assessment of thousands to millions of crystal structures and the greatest possible accuracy in their relative energies. We describe a crystal lattice simulation program, DMACRYS, emphasizing the features that make it suitable for use in crystal structure prediction for pharmaceutical molecules using accurate anisotropic atom-atom model intermolecular potentials based on the theory of intermolecular forces. DMACRYS can optimize the lattice energy of a crystal, calculate the second derivative properties, and reduce the symmetry of the spacegroup to move away from a transition state. The calculated terahertz frequency k = 0 rigid-body lattice modes and elastic tensor can be used to estimate free energies. The program uses a distributed multipole electrostatic model (Q, t = 00,...,44s) for the electrostatic fields, and can use anisotropic atom-atom repulsion models, damped isotropic dispersion up to R(-10), as well as a range of empirically fitted isotropic exp-6 atom-atom models with different definitions of atomic types. A new feature is that an accurate model for the induction energy contribution to the lattice energy has been implemented that uses atomic anisotropic dipole polarizability models (alpha, t = (10,10)...(11c,11s)) to evaluate the changes in the molecular charge density induced by the electrostatic field within the crystal. It is demonstrated, using the four polymorphs of the pharmaceutical carbamazepine C(15)H(12)N(2)O, that whilst reproducing crystal structures is relatively easy, calculating the polymorphic energy differences to the accuracy of a few kJ mol(-1) required for applications is very demanding of assumptions made in the modelling. Thus DMACRYS enables the comparison of both known and hypothetical crystal structures as an aid to the development of pharmaceuticals and other speciality organic materials, and provides a tool to develop the modelling of the intermolecular forces involved in molecular recognition processes.

Solid-solid phase transitions in DL-norvaline studied by differential scanning calorimetry and Raman spectroscopy

The structural modifications of the amino acid DL-Norvaline have been studied using differential scanning calorimetry (DSC) and Raman spectroscopy. DSC results showed that this amino acid undergoes two solid-solid phase transitions at -116.9 and -76.1 degrees C in the temperature range -130 to +300 degrees C. Raman spectroscopy was applied to complement DSC results. The combination of the two methodologies point out that the observed phase transitions correspond to an increment of disordering in the aliphatic side chain of amino acid, an augmentation of the rotational motion of the amino group and a decrease of the strength of the intramolecular hydrogen bonding of the initial dimers at low temperatures. The observed phase transitions of DL-norvaline are compared with those found in DL-norleucine.

In Silico Prediction of Drug Solubility. 3. Free Energy of Solvation in Pure Amorphous Matter

The solubility of drugs in water is investigated in a series of papers. In this work, we address the process of bringing a drug molecule from the vapor into a pure drug amorphous phase. This step enables us to actually calculate the solubility of amorphous drugs in water. In our general approach, we, on one hand, perform rigorous free energy simulations using a combination of the free energy perturbation and thermodynamic integration methods. On the other hand, we develop an approximate theory containing parameters that are easily accessible from conventional Monte Carlo simulations, thereby reducing the computation time significantly. In the theory for solvation, we assume that DeltaG* = DeltaGcav + ELJ + EC/2, where the free energy of cavity formation, DeltaGcav, in pure drug systems is obtained using a theory for hard-oblate spheroids, and ELJ and EC are the Lennard-Jones and Coulomb interaction energies between the chosen molecule and the others in the fluid. The theoretical predictions for the free energy of solvation in pure amorphous matter are in good agreement with free energy simulation data for 46 different drug molecules. These results together with our previous studies support our theoretical approach. By using our previous data for the free energy of hydration, we compute the total free energy change of bringing a molecule from the amorphous phase into water. We obtain good agreement between the theory and simulations. It should be noted that to obtain accurate results for the total process, high precision data are needed for the individual subprocesses. Finally, for eight different substances, we compare the experimental amorphous and crystalline solubility in water with the results obtained by the proposed theory with reasonable success.

Concerted molecular displacements in a thermally-induced solid-state transformation in crystals of DL-norleucine

Martensitic transformations are of considerable technological importance, a particularly promising application being the possibility of using martensitic materials, possibly proteins, as tiny machines. For organic crystals, however, a molecular level understanding of such transformations is lacking. We have studied a martensitic-type transformation in crystals of the amino acid DL-norleucine using molecular dynamics simulation. The crystal structures of DL-norleucine comprise stacks of bilayers (formed as a result of strong hydrogen bonding) that translate relative to each other on transformation. The simulations reveal that the transformation occurs by concerted molecular displacements involving entire bilayers rather than on a molecule-by-molecule basis. These observations can be rationalized on the basis that at sufficiently high excess temperatures, the free energy barriers to concerted molecular displacements can be overcome by the available thermal energy. Furthermore, in displacive transformations, the molecular displacements can occur by the propagation of a displacement wave (akin to a kink in a carpet), which requires the molecules to overcome only a local barrier. Concerted molecular displacements are therefore considered to be a significant feature of all displacive transformations. This finding is expected to be of value toward developing strategies for controlling or modulating martensitic-type transformations.

What is Syncrystallization? States of the pH Indicator Methyl Red in Crystals of Phthalic Acid

The concept of syncrystallization was reinvestigated by focusing on phthalic acid (PA) grown with methyl red (MR). Crystals are alternately red and yellow in adjacent growth sectors. X-ray structures of MR and its cocrystals, revealing MR in the neutral, zwitterionic, and protonated states, as well as measurements of linear birefringence and linear dichroism of mixed crystals, were used to investigate mechanisms of PA coloring. These experiments were complemented by force field calculations of the lowest energy stable surfaces of expressed facets and energies of MR on and in crystals, as well as molecular orbital calculations of MR. Two MR species were detected in PA having distinct energies, polarizations, and face selectivities. Assignments of structures to these MRs, previously thought to be neutral and protonated, required a nuanced analysis of hydrogen bonds. The essential difference between yellow and red species is whether the MR carboxylic acid proton is inter- or intramolecularly hydrogen bound. Inferences about mixed crystal structure drawn from an examination of cocrystals of PA and MR are inconsistent with polarization spectroscopy signaling caution when using stoichiometric compounds as models of dilute solid solutions. Upon heating mixed crystals, linear dichroism diminishes and oriented, elongated pools of MR separate and pass through the bulk in directions perpendicular to the direction of elongation. These bâtonnets subsequently crystallize leaving macroscopic oriented crystals of a MR-rich phase within PA. No evidence was found for the simultaneous crystallization of MR and PA; however, the MR reorientation on heating as well as the separation and recrystallization of a MR-rich phase are distinct processes that could be embraced by the literal meaning of syncrystallization.