Mesoscale Clusters in the Crystallisation of Organic Molecules.

The early stages of the molecular self-assembly pathway leading to crystal nucleation have a significant influence on the properties and purity of organic materials. This mini review collates the work on organic mesoscale clusters and discusses their importance in nucleation processes, with a particular focus on their critical properties and susceptibility to sample treatment parameters. This is accomplished by a review of detection methods, including dynamic light scattering, nanoparticle tracking analysis, small angle X-ray scattering, and transmission electron microscopy. Considering the challenges associated with crystallisation of flexible and large-molecule active pharmaceutical ingredients, the dynamic nature of mesoscale clusters has the potential to expand the discovery of novel crystal forms. By collating literature on mesoscale clusters for organic molecules, a more comprehensive understanding of their role in nucleation will evolve and can guide further research efforts.

Crystal Nucleation of Tolbutamide in Solution: Relationship to Solvent, Solute Conformation, and Solution Structure.

The influence of the solvent in nucleation of tolbutamide, a medium-sized, flexible and polymorphic organic molecule, has been explored by measuring nucleation induction times, estimating solvent-solute interaction enthalpies using molecular modelling and calorimetric data, probing interactions and clustering with spectroscopy, and modelling solvent-dependence of molecular conformation in solution. The nucleation driving force required to reach the same induction time is strongly solvent-dependent, increasing in the order: acetonitrile

Solute clustering in undersaturated solutions - systematic dependence on time, temperature and concentration.

Molecular clustering and solvent-solute interactions in isopropanol solutions of fenoxycarb have been thoroughly and systematically investigated by dynamic light scattering, small-angle X-ray scattering, and nanoparticle tracking, supported by infrared spectroscopy and molecular dynamics simulations. The existence of molecular aggregates, clusters, ranging in size up to almost a micrometre is clearly recorded at undersaturated as well as supersaturated conditions by all three analysis techniques. The results systematically reveal that the cluster size increases with solute concentration and time at stagnant conditions. For most concentrations the time scale of cluster growth is of the order of days. In undersaturated solutions the size appears to eventually reach a maximum value, higher the higher the concentration. Below a certain concentration threshold clusters are significantly smaller. Clusters are found to be smaller in solutions pre-heated at a higher temperature, which offers a possible explanation for the so-called "history of solution" effect. The cluster distribution is influenced by filtration through membranes with a pore size of 0.1 µm, offering an alternative explanation for the "foreign particle-catalysed nucleation" effect. At moderate concentrations larger clusters appear to be sheared into smaller ones, but the original size distribution is rapidly re-established. At higher concentrations, although still well below solubility, the cluster size as well as solute concentration are strongly affected, suggesting that larger clusters contain at least a core of more organized molecules not able to pass through the filter.

Influence of solvent on crystal nucleation of risperidone.

Over 2100 induction time experiments were carried out for the medium-sized, antipsychotic drug molecule, risperidone in seven different organic solvents. To reach the same induction time the required driving force increases in the order: cumene, toluene, acetone, ethyl acetate, methanol, propanol, and butanol, which reasonably well correlates to the interfacial energies as determined within classical nucleation theory. FTIR spectroscopy has been used to investigate any shifts in the spectra and to estimate the interaction of solute and solvent at the corresponding site. The solution condition has also been investigated by Density Functional Theory (DFT) calculations over (1 : 1) solvent-solute binding interactions at 8 different sites on the risperidone molecule. The DFT computational results agree with the spectroscopic data suggesting that these methods do capture the binding strength of solvent molecules to the risperidone molecule. The difficulty of nucleation correlates reasonably to the DFT computations and the spectroscopic measurements. The results of the different measurements suggest that the stronger the solvent binds to the risperidone molecule in solution, the slower the nucleation becomes.

Investigating the role of solvent-solute interaction in crystal nucleation of salicylic acid from organic solvents.

In previous work, it has been shown that the crystal nucleation of salicylic acid (SA) in different solvents becomes increasingly more difficult in the order: chloroform, ethyl acetate acetonitrile, acetone, methanol, and acetic acid. In the present work, vibration spectroscopy, calorimetric measurements, and density functional theory (DFT) calculations are used to reveal the underlying molecular mechanisms. Raman and infrared spectra suggest that SA exists predominately as dimers in chloroform, but in the other five solvents there is no clear evidence of dimerization. In all solvents, the shift in the SA carbonyl peak reflecting the strength in the solvent-solute interaction is quite well correlated to the nucleation ranking. This shift is corroborated by DFT calculated energies of binding one solvent molecule to the carboxyl group of SA. An even better correlation of the influence of the solvent on the nucleation is provided by DFT calculated energy of binding the complete first solvation shell to the SA molecule. These solvation shell binding energies are corroborated by the enthalpy of solvent-solute interaction as estimated from experimentally determined enthalpy of solution and calculated enthalpy of cavity formation using the scaled particle theory. The different methods reveal a consistent picture and suggest that the stronger the solvent binds to the SA molecule in solution, the slower the nucleation becomes.

Investigation into the Nucleation of the p-Hydroxybenzoic Acid:Glutaric Acid 1:1 Cocrystal from Stoichiometric and Non-Stoichiometric Solutions.

The nucleation in the p-hydroxybenzoic acid:glutaric acid 1:1 cocrystal (PHBA:GLU) system has been investigated in stoichiometric and non-stoichiometric acetonitrile solutions by induction time experiments. Utilizing the ternary phase diagram, the supersaturated non-stoichiometric solutions were created with compositions along the invariant point boundary lines. In all cases, the PHBA:GLU cocrystal was the nucleating phase, even though the non-stoichiometric solutions were also supersaturated with respect to the pure solid phases. The nucleation of the cocrystal from the mixed solutions is found to be more difficult than the nucleation of the pure compounds from the respective pure solutions, as captured by lower pre-exponential factors (A). However, if the driving force is defined per reactant molecule instead of per heterodimer, the cocrystal nucleation difficulty is close to that of the more difficult-to-nucleate pure compound. The difference in nucleation difficulty of the cocrystal from stoichiometric and non-stoichiometric solutions was captured by differences in the interfacial energy, while the pre-exponential factor remained unchanged. Apart from the pure GLU system, the relation between the experimentally determined pre-exponential factors for the different systems correlates with calculated values using theoretical expressions for volume-diffusion and surface-integration control.

Characterization and Crystal Nucleation Kinetics of a New Metastable Polymorph of Piracetam in Alcoholic Solvents.

A new polymorph of the drug active pharmaceutical ingredient piracetam (Form VI) has been discovered and characterized by X-ray powder diffraction (PXRD), solid-state Raman, attenuated total reflectance infrared spectroscopy, and differential scanning calorimetry. The PXRD diffractogram of Form VI shows a distinct peak at 24.2° (2θ) that distinguishes it from the previously known polymorphs and solvates. Form VI is metastable with respect to the previously known polymorphs Form II and Form III; in ethanol solution at 288 K, Form VI transforms into Form II within 15 min, while in isopropanol solution Form VI is kinetically stable for at least 6 h. A total of 1200 crystal nucleation induction time experiments of piracetam in ethanol and isopropanol solutions have been conducted, in sets of 40-80 repeat experiments carried out at different temperatures and solute concentrations. Each solution nucleated as a single polymorph, and each set of repeat experiments resulted in different proportions of Form II, Form III, and Form VI, with Form VI dominating at low nucleation temperatures and Form II at higher nucleation temperatures. The induction time data for Form VI at 288 K have been evaluated within the framework of the classical nucleation theory. At equal driving force, nucleation of Form VI is less obstructed in ethanol than in isopropanol, as captured by a lower interfacial energy and higher pre-exponential factor in ethanol. The proportion of Form VI obtained at a comparable driving force increases in the order ethanol < isopropanol.

Nucleation of the Theophylline:Salicylic Acid 1:1 Cocrystal.

The nucleation behavior of the theophylline-salicylic acid 1:1 (THP:SA) cocrystal in chloroform has been investigated and compared with the corresponding behavior of the pure compounds. Induction times have been determined at different supersaturations at 10 °C under each condition in approximately 40-80 repetition experiments in 20 mL vials. Nucleation times, extracted from the median induction times by accounting for a nucleus growth time, have been used to determine the interfacial energy and the pre-exponential factor within the classical nucleation theory. Results show that the cocrystal at equal driving force has a longer nucleation time, or to reach equal nucleation time, the cocrystal requires a higher driving force. Pure theophylline is easier to nucleate than pure salicylic acid, despite the latter having a smaller molecular size, higher solubility, and is expected to form dimers already in the solution. The cocrystal is found to have an interfacial energy in between the respective values for the pure compounds. However, the higher molecular volume of the cocrystal, taken as the volume of the 1:1 theophylline-salicylic acid assembly, leads to the highest nucleation work, which, together with a low pre-exponential factor, explains why the cocrystal is the most difficult to nucleate. The experimentally extracted pre-exponential factor of the cocrystal is very similar to that of THP, and similar trends are observed from theoretical expressions of volume-diffusion- and surface-integration-controlled nucleation, respectively.

Pure Curcumin Spherulites from Impure Solutions via Nonclassical Crystallization.

Crystallization experiments performed with highly supercooled solutions produced highly pure (>99 wt %) and highly crystalline mesocrystals of curcumin from impure solutions (∼22% of two structurally similar impurities) in one step. These mesocrystals exhibited a crystallographic hierarchy and were composed of perfectly or imperfectly aligned nanometer-thick crystallites. X-ray diffraction and spectroscopic analysis confirmed that the spherulites are a new solid form of curcumin. A theoretical hypothesis based on particle aggregation, double nucleation, and repeated secondary nucleation is proposed to explain the spherulite formation mechanism. The experimental results provide, for the first time, evidence for an organic molecule to naturally form spherulites without the presence of any stabilizing agents. Control experiments performed with highly supercooled pure solutions produced spherulites, confirming that the formation of spherulites is attributed to the high degree of supercooling and not due to the presence of impurities. Likewise, control experiments performed with a lower degree of supercooling produced impure crystals of curcumin via classical molecular addition mechanisms. Collectively, these experimental observations provide, for the first time, evidence for particle-mediated crystallization as an alternate and efficient method to purify organic compounds.

Thermodynamics of the Enantiotropic Pharmaceutical Compound Benzocaine and Solubility in Pure Organic Solvents.

The thermodynamic relationship between FI and FII of ethyl 4-aminobenzoate (benzocaine) has been investigated. Slurry conversion experiments show that the transition temperature below which FI is stable is located between 302 K-303 K (29 °C-30 °C). The polymorphs FI and FII have been characterised by infrared spectroscopy (IR), Raman spectroscopy, transmission powder X-ray diffraction (XRPD) and differential scanning calorimetry (DSC). The isobaric solid state heat capacities have been measured by DSC. The quantitative thermodynamic stability relationship has been determined in a comprehensive thermodynamic analysis of the calorimetric data. The solubility of both polymorphs has been determined in eight pure organic solvents over the temperature range 278 K-323 K by a gravimetric method. The mole fraction solubility of benzocaine decreases in the order: 1,4-dioxane, acetone, ethyl acetate, chloroform, acetonitrile, methanol, n-butanol and toluene. Comparison with the determined activity of solid benzocaine forms shows that negative deviation from Raoult's law ideality is found in dioxane, acetone and ethyl acetate solutions, and positive deviation in solutions of the other investigated solvents.

Solubility of Two Polymorphs of Tolbutamide in n-Propanol: Comparison of Methods.

The solid-liquid solubility of two polymorphs of the title compound has been measured in n-propanol over the temperature range (278 K-303 K) by an isothermal, gravimetric method and a low heating rate polythermal method. Due to marked differences in the settling behavior of crystals of the two polymorphs in the investigated solvent, it is found that the low heating rate polythermal method gives the overall best performance for this particular system. Systematic slurry conversion experiments show that FII is the stable polymorph over the investigated temperature range (268 K-308 K). Solubility data for both polymorphs is well correlated, and has been extrapolated to the melting point, by a previously proposed semi-empirical regression model based on solid-phase calorimetric data. The system exhibits a marked positive deviation from Raoult's law, with solute activity coefficients at equilibrium decreasing with increasing temperature.

Solubility and thermodynamic analysis of ketoprofen in organic solvents.

The solubility of the racemic solid phase of ketoprofen (KTP) in methanol, ethanol, isopropanol, butanol, acetonitrile, ethyl acetate, 1,4-dioxane and toluene has been determined between 273 and 303 K by a gravimetric method. FTIR and Raman spectroscopy, SEM and PXRD, have been used to characterise the solid phase. The melting data and heat capacity of solid and melt have been determined by DSC, and used to estimate fusion thermodynamics and the activity of the solid phase as functions of temperature. Empirical and semi-empirical models have been fitted to experimental solubility data. The solution activity coefficients reveal positive deviation from ideality in all solvents except for in dioxane, and very close to ideality in methanol. The solubility is fairly high in the alcohols but decrease with increasing hydrocarbon chain. Generally and due to the presence of the carboxylic acid group, KTP is more readily dissolved in polar protic solvents, followed in order by polar aprotic and non-polar solvents. However, the highest solubility is found in dioxane, classified as a non-polar solvent, but notably though the molecule having two strong hydrogen bond accepting functionalities, and no hydrogen bond donation capability.

Solid and Solution State Thermodynamics of Polymorphs of Butamben (Butyl 4-Aminobenzoate) in Pure Organic Solvents.

The solubility of butamben has been measured gravimetrically in pure methanol, 1-propanol, 2-propanol, 1-butanol, and toluene over the temperature range 268-298 K. Polymorph transition and melting temperatures, associated enthalpy changes, and the heat capacity of the solid forms and the supercooled melt have been measured by differential scanning calorimetry. Based on extrapolated calorimetric data, the Gibbs energy, enthalpy and entropy of fusion, and the activity of solid butamben (the ideal solubility) have been calculated from below ambient temperature up to the melting point. Activity coefficients of butamben at equilibrium in the different solvents have been estimated from solubility data and the activity of the solid, revealing that all investigated systems exhibit positive deviation from Raoult's law. Solubility data are well correlated by a semiempirical regression model. On a mass basis, the solubility is clearly higher in methanol than in the other solvents, but mole fraction solubilities are very similar across all 5 solvents. The 2 known polymorphs are enantiotropically related, and the transition point is located at 283 K. Polymorph interconversions occur within 0.3 K of the transition point even in the solid state, and the 2 forms exhibit strong similarities in investigated properties.

Spherical crystallization of benzoic acid.

This paper deals with the development of a method for spherical crystallization of benzoic acid. Benzoic acid is dissolved in ethanol, water is used as anti-solvent and chloroform is used as bridging liquid. After an introductory screening of different methods, the influence of the amount of the bridging liquid, the solute concentration and the stirring rate is investigated. The product particle characterization includes the particle size distribution, morphology and strength. The mechanical strength of single agglomerates has been determined by compression in a materials testing machine, using a 10N load cell. It is found that favourable properties are obtained if the bridging liquid is added during the crystallization. Larger and stronger well-shaped agglomerates are formed. The stress-strain curves are J-shaped with no clear fracturing of the particles, and are well correlated by an exponential-polynomial equation.

Oiling out or molten hydrate-liquid-liquid phase separation in the system vanillin-water.

Vanillin crystals in a saturated aqueous solution disappear and a second liquid phase emerges when the temperature is raised above 51 degrees C. The phenomenon has been investigated with crystallization and equilibration experiments, using DSC, TGA, XRD and hot-stage microscopy for analysis. The new liquid solidifies on cooling, appears to melt at 51 degrees C, and has a composition corresponding to a dihydrate. However, no solid hydrate can be detected by XRD, and it is shown that the true explanation is that a liquid-liquid phase separation occurs above 51 degrees C where the vanillin-rich phase has a composition close to a dihydrate. To our knowledge, liquid-liquid phase separation has not previously been reported for the system vanillin-water, even though thousands of tonnes of vanillin are produced globally every year.

Effect of preparation method on compactability of paracetamol granules and agglomerates.

The objective of this study was to investigate the effect of fracture strength of paracetamol particles on their compactability. For this purpose two series of paracetamol particles were prepared by crystal agglomeration and by granulation using different solvents. A free flowing particle size fraction of all types of particles was characterized with respect to their shape, degree of agglomeration and single fracture strength. The powders were compressed to tablets and the compression mechanism of the particles and the evolution in tablet micro-structure were assessed by compression parameters derived from the Heckel and Kawakita equations and by a tablet permeabililty coefficient. Tablet tensile strength and porosity were determined. The degree of deformation and fragmentation during compression varied between agglomerates and granules and was dependent on their failure strength. The granules varied in compactability with particle failure strength while the agglomerates showed limited variation. It is proposed that, the dominant mechanism of compression for the granules was permanent deformation while for the agglomerates it was fragmentation. It was thus found that the compression mechanism of the particles was dependent on both the degree of agglomeration and the particle failure strength.

Modelling and understanding powder flow properties and compactability of selected active pharmaceutical ingredients, excipients and physical mixtures from critical material properties.

The development of solid dosage forms and manufacturing processes are governed by complex physical properties of the powder and the type of pharmaceutical unit operation the manufacturing processes employs. Suitable powder flow properties and compactability are crucial bulk level properties for tablet manufacturing by direct compression. It is also generally agreed that small scale powder flow measurements can be useful to predict large scale production failure. In this study, predictive multilinear regression models were effectively developed from critical material properties to estimate static powder flow parameters from particle size distribution data for a single component and for binary systems. A multilinear regression model, which was successfully developed for ibuprofen, also efficiently predicted the powder flow properties for a range of batches of two other active pharmaceutical ingredients processed by the same manufacturing route. The particle size distribution also affected the compactability of ibuprofen, and the scope of this work will be extended to the development of predictive multivariate models for compactability, in a similar manner to the approach successfully applied to flow properties.

Phase equilibria and thermodynamics of p-hydroxybenzoic acid.

The prevalence of phases and associated solubilities of p-hydroxybenzoic acid have been investigated in methanol, acetonitrile, acetic acid, acetone, water, and ethyl acetate at temperatures from 10 to 50 degrees C. Thermodynamic data was acquired through determination of van't Hoff enthalpy of solution, enthalpy of fusion, and melting temperature. Indications of polymorphic enantiotropy were found primarily through solubility analysis and FTIR-ATR. A comprehensive thermodynamic investigation disclosed correlation between the van't Hoff enthalpy of solution and the solubility in different solvents. A higher solubility is linked to a lower van't Hoff enthalpy of solution. A thermodynamic analysis to discriminate between different solid phases is presented.

Polymorphism and thermodynamics of m-hydroxybenzoic acid.

Solution and solid-state properties of m-hydroxybenzoic acid have been investigated. Two polymorphs were found where the monoclinic modification exhibits a higher stability than the orthorhombic form. The solubility of the monoclinic polymorph was determined between 10 and 50 degrees C in methanol, acetonitrile, acetic acid, acetone, water and ethyl acetate. The solubility of the orthorhombic polymorph was determined between 10 and 50 degrees C in acetonitrile, acetic acid, acetone and ethyl acetate. A thermodynamic analysis revealed a marked correlation between the molar solubility and the van't Hoff enthalpy of solution at constant temperature. In addition, in each solvent increased temperature resulted in increased van't Hoff enthalpy of solution. It is shown that the solubility data can be used to estimate melting properties for both polymorphs. The solubility ratio of the two forms and the DSC thermogram of the orthorhombic form strongly suggest that the system is monotropic. However, according to the polymorph rules of Burger and Ramberger, the estimated higher melting enthalpy and lower melting temperature of the orthorhombic form points towards an enantiotropic system. Hence, this system appears to be an exception to the Burger and Ramberger melting enthalpy rule, and the probable reason for this is found in the difference in the heat capacity of the two solid forms.

Thermodynamic Stability Analysis of Tolbutamide Polymorphs and Solubility in Organic Solvents.

Melting temperatures and enthalpies of fusion have been determined by differential scanning calorimetry (DSC) for 2 polymorphs of the drug tolbutamide: FI(H) and FV. Heat capacities have been determined by temperature-modulated DSC for 4 polymorphs: FI(L), FI(H), FII, FV, and for the supercooled melt. The enthalpy of fusion of FII at its melting point has been estimated from the enthalpy of transition of FII into FI(H) through a thermodynamic cycle. Calorimetric data have been used to derive a quantitative polymorphic stability relationship between these 4 polymorphs, showing that FII is the stable polymorph below approximately 333 K, above which temperature FI(H) is the stable form up to its melting point. The relative stability of FV is well below the other polymorphs. The previously reported kinetic reversibility of the transformation between FI(L) and FI(H) has been verified using in situ Raman spectroscopy. The solid-liquid solubility of FII has been gravimetrically determined in 5 pure organic solvents (methanol, 1-propanol, ethyl acetate, acetonitrile, and toluene) over the temperature range 278 to 323 K. The ideal solubility has been estimated from calorimetric data, and solution activity coefficients at saturation in the 5 solvents determined. All solutions show positive deviation from Raoult's law, and all van't Hoff plots of solubility data are nonlinear. The solubility in toluene is well below that observed in the other investigated solvents. Solubility data have been correlated and extrapolated to the melting point using a semiempirical regression model.