On the mechanism of solid-state phase transitions in molecular crystals - the role of cooperative motion in (quasi)racemic linear amino acids.

During single-crystal-to-single-crystal (SCSC) phase transitions, a polymorph of a compound can transform to a more stable form while remaining in the solid state. By understanding the mechanism of these transitions, strategies can be developed to control this phenomenon. This is particularly important in the pharmaceutical industry, but also relevant for other industries such as the food and agrochemical industries. Although extensive literature exists on SCSC phase transitions in inorganic crystals, it is unclear whether their classications and mechanisms translate to molecular crystals, with weaker interactions and more steric hindrance. A comparitive study of SCSC phase transitions in aliphatic linear-chain amino acid crystals, both racemates and quasi-racemates, is presented. A total of 34 transitions are considered and most are classified according to their structural change during the transition. Transitions without torsional changes show very different characteristics, such as transition temperature, enthalpy and free energy, compared with transitions that involve torsional changes. These differences can be rationalized using classical nucleation theory and in terms of a difference in mechanism; torsional changes occur in a molecule-by-molecule fashion, whereas transitions without torsional changes involve cooperative motion with multiple molecules at the same time.

Inhibition of the Vapor-Mediated Phase Transition of the High Temperature Form of Pyrazinamide.

Tailor-made additives can prove an effective method to prolong the lifetime of metastable forms of pharmaceutical compounds by surface stabilization. Pyrazinamide (PZA) is a pharmaceutical compound with four polymorphic forms. The high temperature γ form, which can be produced by spray drying or sublimation growth, is metastable at room temperature and transforms within days when produced by spray drying, and within several months up to years for single crystals produced by sublimation. However, when PZA is cospray dried with 1,3-dimethylurea (DMU), it has been reported to remain in its γ form for several years. Scanning electron microscopy (SEM) images reveal that the phase transition from γ-PZA to the low temperature forms involves a vapor-mediated recrystallization, while the reverse phase transition upon heating is a nucleation-and-growth solid-solid phase transition. The lifetime-extending effect of DMU on spray-dried PZA has been investigated in more detail and compared with high-energy ball milling of sublimation-grown γ-PZA crystals. Co-ball milling of PZA and DMU is found to extend the lifetime of the high temperature form of PZA to a few months, while separate ball milling leads to an extension of merely a few weeks. DMU acts as an additive that most likely stabilizes the surface of γ-PZA, which would reduce the vapor pressure of PZA, thereby reducing the transition rate. Alternatively, DMU could prevent nucleation of low temperature forms of PZA.

Searching for stereoisomerism in crystallographic databases: algorithm, analysis and chiral curiosities.

The automated identification of chiral centres in molecular residues is a non-trivial task. Current tools that allow the user to analyze crystallographic data entries do not identify chiral centres in some of the more complex ring structures, or lack the possibility to determine and compare the chirality of multiple structures. This article presents an approach to identify asymmetric C atoms, which is based on the atomic walk count algorithm presented by Rücker & Rücker [(1993), J. Chem. Inf. Comput. Sci. 33, 683-695]. The algorithm, which we implemented in a computer program named ChiChi, is able to compare isomeric residues based on the chiral centres that were identified. This allows for discrimination between enantiomers, diastereomers and constitutional isomers that are present in crystallographic databases. ChiChi was used to process 254 354 organic entries from the Cambridge Structural Database (CSD). A thorough analysis of stereoisomerism in the CSD is presented accompanied by a collection of chiral curiosities that illustrate the strength and versatility of this approach.

The solubility behaviour and thermodynamic relations of the three forms of Venlafaxine free base.

The polymorphic and solubility behaviour of the active pharmaceutical ingredient Venlafaxine free base, which is used as an antidepressant, is studied. Using differential scanning calorimetry and slurry experiments, an enantiotropic relation between the three forms was found. Transition temperatures were determined using solubility data and compared with calculated transition temperatures based on the melting enthalpies and temperatures of the different forms. The solubility of Venlafaxine in heptane, toluene and methanol shows a large deviation from ideal behaviour. The deviations are to a large extent determined by the temperature dependence of the difference in fusion enthalpy of the undercooled melt and the solid.

Experimental and computational morphology of three polymorphs of the free base of Venlafaxine: a comparison of morphology prediction methods.

In this paper the experimental and the computational studies of the morphology of three polymorphs of the free base of Venlafaxine ((N,N-dimethyl)-2-(1-hydroxy cyclohex-1-yl)-2-(4-methoxyphenyl) ethylamine) are reported. The morphology of all polymorphs has been predicted using the Bravais-Friedel-Donnay-Harker method, the attachment energy method and kinetic Monte Carlo growth simulations and these predictions have been compared with experimental observations. The Monte Carlo simulations allow for a detailed simulation of the growth process, including driving force and growth mechanism, which leads to a semi-quantitative prediction of the growth morphologies of all three phases. For phase I two distinct growth habits are found experimentally under the same conditions. This is explained by the occurrence of a spiral growth mechanism in one of the two, which was observed using AFM and which is also supported by the Monte Carlo simulations. The habit of phase II could only be explained from simulations when a spiral growth mechanism is assumed; the shape of phase III could not be modeled accurately from the Monte Carlo simulations. Although the shape of the crystal is reproduced accurately, some of the indices of the faces predicted are not in agreement with the indices measured. The deviations are interpreted to be due to the presence of domains in the crystals as a result of the layered structure.

The hydration/dehydration behavior of aspartame revisited.

Aspartame, l-aspartyl-l-phenylalanine methyl ester, has two hydrates (IA and IB), a hemi-hydrate (IIA) and an anhydrate (IIB). The hydration/dehydration behavior of aspartame was investigated using hot-humidity stage X-ray powder diffraction (XRPD) and molecular mechanics modeling in combination with differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA). The results of this study are compared to earlier studies on aspartame as described in literature. It is shown that earlier transition studies were hampered by incomplete conversions and wrong assignment of the forms. The combination of the techniques applied in this study now shows consistent results for aspartame and yields a clear conversion scheme for the hydration/dehydration behavior of the four forms.

Solubility of molecular crystals: polymorphism in the light of solubility theory.

The thermodynamic theory of solubility of molecular crystals in solvents is reviewed with an emphasis on solutes showing polymorphism as in case of many pharmaceuticals. The relation between solubility and binary phase diagrams of the solute solvent system is treated. The astonishing variety of possible solubility curves as a function of temperature is explained using simple models for the solution thermodynamics assuming no mixing between the solvent and solute in the solid phase, though including the case of solvates or pseudo polymorphs. In passing a new method is introduced that allows to estimate the transition temperature of enantiotropically related polymorphs from melting temperatures and enthalpies of the polymorphs.

On the definition of a Monte Carlo model for binary crystal growth.

We show that consistency of the transition probabilities in a lattice Monte Carlo (MC) model for binary crystal growth with the thermodynamic properties of a system does not guarantee the MC simulations near equilibrium to be in agreement with the thermodynamic equilibrium phase diagram for that system. The deviations remain small for systems with small bond energies, but they can increase significantly for systems with large melting entropy, typical for molecular systems. These deviations are attributed to the surface kinetics, which is responsible for a metastable zone below the liquidus line where no growth occurs, even in the absence of a 2D nucleation barrier. Here we propose an extension of the MC model that introduces a freedom of choice in the transition probabilities while staying within the thermodynamic constraints. This freedom can be used to eliminate the discrepancy between the MC simulations and the thermodynamic equilibrium phase diagram. Agreement is achieved for that choice of the transition probabilities yielding the fastest decrease of the free energy (i.e., largest growth rate) of the system at a temperature slightly below the equilibrium temperature. An analytical model is developed, which reproduces quite well the MC results, enabling a straightforward determination of the optimal set of transition probabilities. Application of both the MC and analytical model to conditions well away from equilibrium, giving rise to kinetic phase diagrams, shows that the effect of kinetics on segregation is even stronger than that predicted by previous models.

Crystal structure prediction of organic pigments: quinacridone as an example.

The structures of the alpha, beta and gamma polymorphs of quinacridone (Pigment Violet 19) were predicted using Polymorph Predictor software in combination with X-ray powder diffraction patterns of limited quality. After generation and energy minimization of the possible structures, their powder patterns were compared with the experimental ones. On this basis, candidate structures for the polymorphs were chosen from the list of all structures. Rietveld refinement was used to validate the choice of structures. The predicted structure of the gamma polymorph is in accordance with the experimental structure published previously. Three possible structures for the beta polymorph are proposed on the basis of X-ray powder patterns comparison. It is shown that the alpha structure in the Cambridge Structural Database is likely to be in error, and a new alpha structure is proposed. The present work demonstrates a method to obtain crystal structures of industrially important pigments when only a low-quality X-ray powder diffraction pattern is available.

Epitaxial 2D nucleation of the stable polymorphic form of the steroid 7alphaMNa on the metastable form: implications for Ostwald's rule of stages.

This paper presents in situ observations of the epitaxial nucleation and growth of the stable polymorph of a steroid, 7alphaMna, on a specific face of the metastable form at low supersaturation, using optical microscopy and in situ Raman spectroscopy. The presence of the metastable polymorph is essential for the nucleation and growth of the stable one. The order of the metastable zones of the stable and metastable polymorphs is reversed for the epitaxial growth process as compared to the case of 3D nucleation. The rate of transformation of the metastable polymorph to the stable one can be controlled by the supersaturation.

Experimental and computational growth morphology of two polymorphs of a yellow isoxazolone dye.

We report on the crystal structures and the experimentally found and the computationally predicted growth morphologies of two polymorphs of a yellow isoxazolone dye. The stable polymorph I has a blocklike habit, and the metastable polymorph II grows as fine needles, nucleating only by heterogeneous or contact nucleation. The habits of both polymorphs depend on the supersaturation during growth. The experimental observations are compared with predictions of the attachment energy model and kinetic Monte Carlo lattice simulations in which the growth is modeled as an "atomistic process", governed by surface energetics. These Monte Carlo simulations correctly predict the shape and the dependence on supersaturation of the crystal morphology for both polymorphs: for polymorph I, a strong dependence on supersaturation is found from the simulations. For polymorph II, the order of morphological importance is reproduced correctly, as well as the needlelike morphology.

Thermodynamics of epitaxial calcite nucleation on self-assembled monolayers.

Classical heterogeneous nucleation theory is used to describe the epitaxial nucleation of calcite on self-assembled monolayers (SAMs). Both spherical and faceted clusters are considered. The use of faceted clusters reveals a useful relation between the shape of very small crystals and the ratio of the heterogeneous and homogeneous nucleation barriers. The experimental approach of this paper concerns the measurement of the threshold driving forces for both homogeneous and heterogeneous nucleation of calcite. This is accomplished by preparing solutions with well-defined driving forces and by measuring the resulting types of nucleation that are observed after a fixed experimental time. The results of the experiments and the theoretical shape analysis are compared, and it is shown that in the experiments no homogeneous nucleation of calcite occurs for driving forces up to at least Deltamu/k(B)T approximately equal to 6.0. A calculation of the critical cluster size for heterogeneous nucleation results in a range of 2-28 growth units and faceted critical clusters from 3-28 growth units, depending on the value of the surface free energy of calcite. These sizes are 50-100 times smaller than the crystalline domain sizes of SAMs and therefore small enough to explain the promoting effect of the substrate.

On the prediction of crystal morphology. III. Equilibrium and growth behaviour of crystal faces containing multiple connected nets.

In this paper, the equilibrium and growth behaviour of faces (hkl) with more than one connected net is studied. It is shown that for these types of orientation different surface phases exist under equilibrium conditions as a function of temperature. Depending on the exact bonding topology at the surface, flat, rough or disordered flat phases are found. Moreover, the growth rate R(hkl) of such faces can differ significantly from the usually calculated relative growth rates based on the attachment energy. Monte Carlo simulations confirm the results from the Hartman-Perdok analyses and offer a tool for the prediction of the crystal habit as a function of supersaturation.

Microtubules in root hairs.

The microtubules of root hairs of Raphanus sativus, Lepidium sativum, Equisetum hyemale, Limnobium stoloniferum, Ceratopteris thalictroides, Allium sativum and Urtica dioica were investigated using immunofluorescence and electron microscopy. Arrays of cortical microtubules were observed in all hairs. The microtubules in the hairs show net axial orientations, but in Allium and Urtica helical microtubule patterns are also present. Numerical parameters of microtubules in Raphanus, Equisetum and Limnobium were determined from dry-cleave preparations. The results are discussed with respect to cell wall deposition and cell morphogenesis.