Computational Study of Solvent Effects on the Molecular Self-Assembly of Tetrolic Acid in Solution and Implications for the Polymorph Formed from Crystallization

Biased Bowl-Direction of Monofluorosumanene in the Solid State

A new curved π-conjugated molecule 1-fluorosumanene (1) was designed and synthesized that possesses one fluorine atom on the benzylic carbon of sumanene. This compound can exhibit bowl inversion in solution, leading to the formation of two diastereomers, 1endo and 1exo, with different dipole moments. Experimental and theoretical investigation revealed an energetical relationship among 1exo, 1endo, and solvent to realize the various endo:exo ratios in the single crystals of 1 depending on the crystallization solvent. Significantly, the molecular dynamics (MD) simulations revealed that 1exo positively worked for the elongation of the stacking structure and the final endo:exo ratio was affected by the relative stability difference between 1endo and 1exo derived by solvation. Such an arrangeable endo:exo ratio of 1 realized the preparation of unique materials showing a different dielectric response from the same molecule 1 just by changing the crystallization solvent.

Molecular Mechanism of Organic Crystal Nucleation: A Perspective of Solution Chemistry and Polymorphism

Crystal nucleation determining the formation and assembly pathway of first organic materials is the central science of various scientific disciplines such as chemical, geochemical, biological, and synthetic materials. However, our current understanding of the molecular mechanisms of nucleation remains limited. Over the past decades, the advancements of new experimental and computational techniques have renewed numerous interests in detailed molecular mechanisms of crystal nucleation, especially structure evolution and solution chemistry. These efforts bifurcate into two categories: (modified) classical nucleation theory (CNT) and non-classical nucleation mechanisms. In this review, we briefly introduce the two nucleation mechanisms and summarize current molecular understandings of crystal nucleation that are specifically applied in polymorphic crystallization systems of small organic molecules. Many important aspects of crystal nucleation including molecular association, solvation, aromatic interactions, and hierarchy in intermolecular interactions were examined and discussed for a series of organic molecular systems. The new understandings relating to molecular self-assembly in nucleating systems have suggested more complex multiple nucleation pathways that are associated with the formation and evolution of molecular aggregates in solution.

Fourier Transform Infrared Spectroscopy and Vibrational Circular Dichroism Assisted Elucidation of the Solution-State Supramolecular Speciation in Racemic and Enantiopure Ketoprofen

The molecular structure and solution-state molecular interactions in the popular non-steroidal anti-inflammatory drug, ketoprofen, are extensively studied with the aim of gaining a better understanding of the chemical behavior of its solution state and its connection to its nucleation pathway and crystallization outcome. Using as reference solid-state X-ray structures of enantiomeric and racemic forms of ketoprofen, a set of self-assembly models underpinned by density functional theory calculations has been considered for the analysis of spectroscopic data, infrared (IR) and vibrational circular dichroism (VCD), obtained for solutions of the samples as a function of composition and solvent. From our results it can be concluded that, contrary to the general belief for generic carboxylic acids, there are no cyclic dimeric structures of ketoprofen present in solution, but rather linear arrays made up of two (in high polar or diluted media) or more units (in low polar or low dilution media). This observation is in line with the idea that the weak contacts (other than H-bonding) would hold the key to molecular self-assembly, in agreement with recent studies on other aromatic carboxylic acids.

Pre-Nucleation Clusters Predict Crystal Structures in Models of Chiral Molecules

Kinetics can play an important role in the crystallization of molecules and can give rise to polymorphism, the tendency of molecules to form more than one crystal structure. Current computational methods of crystal structure prediction, however, focus almost exclusively on identifying the thermodynamically stable polymorph. Kinetic factors of nucleation and growth are often neglected because the underlying microscopic processes can be complex and accurate rate calculations are numerically cumbersome. In this work, we use molecular dynamics computer simulations to study simple molecular models that reproduce the crystallization behavior of real chiral molecules, including the formation of enantiopure and racemic crystals, as well as polymorphism. A significant fraction of these molecules forms crystals that do not have the lowest free energy. We demonstrate that at high supersaturation crystal formation can be accurately predicted by considering the similarities between oligomeric species in solution and molecular motifs in the crystal structure. For the case of racemic mixtures, we even find that knowledge of crystal free energies is not necessary and kinetic considerations are sufficient to determine if the system will undergo spontaneous chiral separation. Our results suggest conceptually simple ways of improving current crystal structure prediction methods.

Unveiling the self-association and desolvation in crystal nucleation

As the first step in the crystallization process, nucleation has been studied by many researchers. In this work, phenacetin (PHEN) was selected as a model compound to investigate the relationship between the solvent and nucleation kinetics. Induction times at different supersaturation in six solvents were measured. FTIR and NMR spectroscopy were employed to explore the solvent-solute interactions and the self-association properties in solution. Density functional theory (DFT) was adopted to evaluate the strength of solute-solvent interactions and the molecular conformations in different solvents. Based on these spectroscopy data, molecular simulation and nucleation kinetic results, a comprehensive understanding of the relationship between molecular structure, crystal structure, solution chemistry and nucleation dynamics is discussed. Both the solute-solvent interaction strength and the supramolecular structure formed by the self-association of solute molecules affect the nucleation rate. The findings reported here shed new light on the molecular mechanism of nucleation in solution.

Observation of crystallisation dynamics by crystal-structure-sensitive room-temperature phosphorescence from Au(I) complexes

The aggregation behaviour of Au(I) complexes in condensed phases can affect their emission properties. Herein, aggregation-induced room-temperature phosphorescence (RTP) is observed from the crystals of trinuclear Au(I) complexes. The RTP is highly sensitive to the crystal structure, with a slight difference in the alkyl side chains causing not only a change in the crystal structure but also a shift in the RTP maximum. Furthermore, in nanocrystals, reversible RTP colour changes are induced by phase transitions between crystal polymorphs during crystal growth from solution or the pulverisation of bulk crystals. The colour change mechanism is discussed in terms of intermolecular interactions in the crystal structure of the luminescent aggregates. The results suggest that the behaviour in nanocrystals may differ from that in bulk crystals. These insights will advance the fundamental understanding of crystallisation mechanisms and may aid in the discovery of new materials properties for solids with nano- to micrometre sizes.

Insight into the role of pre-assembly and desolvation in crystal nucleation: a case of p-nitrobenzoic acid

As one of the most important phenomena in crystallization, the crystal nucleation process has always been the focus of research. In this work, influences of pre-assembly species and the desolvation process on the crystal nucleation process were studied. p-Nitrobenzoic acid (PNBA) was taken as a model compound to investigate the relationship between solution chemistry and nucleation kinetics in seven different solvents. One unsolvated form and four solvates of PNBA were obtained and one of the solvates was newly discovered. The nucleation behaviours and nucleation kinetics of PNBA in the seven solvents were studied and analyzed. Density functional theory (DFT) and solvation energy calculation were adopted to evaluate the strength of solute-solvent interactions. Vibrational spectroscopy combined with molecular simulation was applied to reveal the pre-assembly species in the solution. Based on these results, a comprehensive understanding of the relationship between molecular structure, crystal structure, solution chemistry and nucleation dynamics was proposed and discussed. It was found that the structural similarity between solution chemistry and crystal structure, the interaction between specific sites and the overall strength of solvation will jointly affect the nucleation process.

Face indexing and shape analysis of salicylamide crystals grown in different solvents

The effect of solvent on salicylamide's crystal habit was investigated. It is deduced that ethyl acetate is adsorbed more strongly on the faces, the increased size of which, can explain the shape change.

Tautomer selection through solvate formation: the case of 5-hydroxynicotinic acid

Different 5-hydroxynicotinic acid tautomers were selectively captured through solvate formation. The selectivity is lost once the memory of solvation is erased by removing the solvent from the crystal lattice.

Molecular interaction transfer among solvents and solutes modulates the formation of linezolid crystals

Molecular interaction transfer from solute–solvent to solute–solute pairs is the vital step for the nucleation of linezolid crystals.

Effects of Hydrogen Bond Acceptor Ability of Solvents on Molecular Self-Assembly of Sulfadiazine Solvates

The solvate formation of sulfadiazine (SDZ) was systematically studied in the 4 selected solvents with the aids of experiment and simulation methods. The intermolecular interactions between solute and solvent molecules in different solid states were analyzed and compared through their single crystal structures, and the solution behavior of SDZ was discussed using molecular dynamics simulations. The results indicated that SDZ was easy to form solvates with the solvents having strong hydrogen bond acceptor ability, which determined the formation of hydrogen bonding synthon. Furthermore, the SDZ molecules conformation and packing were compared in various crystal structures. In addition, the desolvation processes of SDZ solvates were studied to investigate the role of solvent in different solvate structures.

Aromatic stacking - a key step in nucleation

Crystal nucleation from solution is of central importance in the chemical and biological sciences. Linking nucleation kinetics to the properties of solutes and solvents remains a grand-challenge in physical chemistry. Through a unique dataset of compounds able to self-assemble via both hydrogen-bonds and aromatic stacking, we are able to compare the importance of these two types of interaction in driving the nucleation process. Contrary to previous reports in which solution chemistry and hydrogen bonding have been seen as controlling factors, we are now able to show that cluster growth via aromatic stacking holds the key.

Role of Solvent Selection on Crystal Habit of 5-Aminosalicylic Acid-Combined Experimental and Computational Approach

Many active pharmaceutical ingredients exhibit a needle-like (acicular) crystal habit, which can significantly complicate their downstream processing. In this study, the acicular crystal habit of a model active pharmaceutical ingredient, 5-aminosalicylic acid (5-ASA), was modified by addition of selected organic solvents to the typical aqueous crystallization process. 5-ASA was crystallized by a pH shift from 7.5-8 to 4 in the presence of methanol, acetonitrile, acetone, tetramethylurea, tetrahydrofuran or dimethyl sulfoxide at 25% v/v, or butanol at 9% v/v. Changes to the experimentally observed crystal habit are rationalized by considering adsorption energy calculations for the solvent molecules onto the morphologically important crystal faces. The crystal habit was influenced most significantly by organic solvents containing a good H-bond acceptor atom, particularly oxygen in acetone, tetramethylurea, tetrahydrofuran, and dimethyl sulfoxide. Such solvents have strongly stabilizing adsorption energies onto the fast-growing crystal faces, and their presence in solution thereby serves to modify the acicular habit of 5-ASA. The developed knowledge base on crystal interface-solvent interactions can form a basis for further engineering of an optimal crystal habit for 5-ASA.

Density functional theory study of the oligomerization of carboxylic acids

We present a density functional theory [M06-2X/6-31+G(d,p)] study of the structures and free energies of formation of oligomers of four carboxylic acids (formic acid, acetic acid, tetrolic acid, and benzoic acid) in water, chloroform, and carbon tetrachloride. Solvation effects were treated using the SMD continuum solvation model. The low-lying energy structures of molecular complexes were located by adopting an efficient search procedure to probe the potential energy surfaces of the oligomers of carboxylic acids (CA)n (n = 2-6). The free energies of the isomers of (CA)n in solution were determined as the sum of the electronic energy, vibrational-rotational-translational gas-phase contribution, and solvation free energy. The assessment of the computational protocol adopted in this study with respect to the dimerization of acetic acid, (AA)2, and formic acid, (FA)2, located new isomers of (AA)2 and (FA)2 and gave dimerization constants in good agreement with the experimental values. The calculation of the self-association of acetic acid, tetrolic acid, and benzoic acid shows the following: (i) Classic carboxylic dimers are the most stable isomer of (CA)2 in both the gas phase and solution. (ii) Trimers of carboxylic acid are stable in apolar aprotic solvents. (iii) Molecular clusters consisting of two interacting classic carboxylic dimers (CA)4,(D+D) are the most stable type of tetramers, but their formation from the self-association of classic carboxylic dimers is highly unfavorable. (iv) For acetic acid and tetrolic acid the reactions (CA)2 + 2CA → (CA)4,(D+D) and (CA)3 + CA → (CA)4,(D+D) are exoergonic, but these aggregation pathways go through unstable clusters that could hinder the formation of tetrameric species. (v) For tetrolic acid the prenucleation species that are more likely to form in solution are dimeric and trimeric structures that have encoded structural motifs resembling the α and β solid forms of tetrolic acid. (vi) Stable tetramers of benzoic acid could form in carbon tetrachloride from the aggregation of trimers and monomers. (vii) Higher order clusters such as acetic acid pentamers and tetrolic acid hexamers are highly unstable in all solvation environments.

Mechanistic insight into the selective crystallization of the metastable polymorph of tolbutamide in ethanol–water solution

A combined experimental and molecular dynamics simulation study discloses the effects of solvents and supersaturation on the tolbutamide polymorphs outcome.

Nucleation of organic crystals--a molecular perspective

The outcome of synthetic procedures for crystalline organic materials strongly depends on the first steps along the molecular self-assembly pathway, a process we know as crystal nucleation. New experimental techniques and computational methodologies have spurred significant interest in understanding the detailed molecular mechanisms by which nuclei form and develop into macroscopic crystals. Although classical nucleation theory (CNT) has served well in describing the kinetics of the processes involved, new proposed nucleation mechanisms are additionally concerned with the evolution of structure and the competing nature of crystallization in polymorphic systems. In this Review, we explore the extent to which CNT and nucleation rate measurements can yield molecular-scale information on this process and summarize current knowledge relating to molecular self-assembly in nucleating systems.

Size-dependent phase stability of a molecular nanocrystal: a proxy for investigating the early stages of crystallization

We make the link between the size-dependent phase stability of a nanocrystal and the phase-transition behavior of emerging crystallites during the earliest stages of crystallization, by using the former as a proxy for the latter. We outline an extension of the classical nucleation theory to describe crystal nucleation and subsequent transformations of competing polymorphic phases that characterize Ostwald's rule of stages. The theoretical framework reveals that the relative stability of the competing phases is a function of cluster size, which in turn varies with time, and therefore explains the complex transformation behavior observed for some systems. We investigated the stability of a nanocrystal of dl-norleucine by means of molecular simulation as a proxy for post-nucleation phase-transformation behavior in emerging crystallites. The simulations reveal that, for nanocrystals, the surface energy of the transition state of a transformation can dominate the barrier to phase change, thus causing metastable phases to be stabilized, not because they are thermodynamically stable, but rather due to kinetic hindering. Therefore, in the context of the earliest stages of crystal growth, not only does phase stability vary as a function of cluster size, and hence time, but thermodynamically feasible transformations are also prone to kinetic hindering.

A computational study of the mechanism of the selective crystallization of α- and β-glycine from water and methanol-water mixture

Understanding the control of polymorphism in organic crystals is of paramount importance to the pharmaceutical, chemical, and food industries. In this work, we investigated two mechanisms described in the literature about the selective crystallization of α- and β-glycine from water and from mixtures of water and methanol using molecular simulations. The link hypothesis (J. Phys. Chem. B 2008, 112, 7794; Cryst. Growth Des. 2006, 6, 1788; J. Inclusion Phenom. Mol. Recognit. Chem. 1990, 8, 395; J. Am. Chem. Soc. 1986, 108, 5871.), which tries to relate the structure of the polymorph obtained from crystallization to the structure of the prenucleation aggregates in the solutions, says the abundance of glycine cyclic dimers in aqueous solutions leads to the crystallization of α-glycine, the polymorph using cyclic dimers as the packing units. This hypothesis was studied first. We revisited the self-assembly of glycine molecules in solution using molecular dynamics to address the debate (Phys. Rev. Lett. 2007, 99, 115702; J. Phys. Chem. B 2008, 112, 7280; J. Am. Chem. Soc. 2008, 130, 13973.) about which is the dominating species in the glycine aqueous solutions and whether there is a link between the solution chemistry and the polymorphic outcome of crystallization. The structures of the glycine clusters were characterized using a structural parameter called cyclic dimer fraction. The glycine clusters in methanol-water mixtures have higher cyclic dimer compositions than those in the pure aqueous solutions. Moreover, the glycine open-chain dimer is more stable than the cyclic dimer regardless of the presence of methanol. All these suggest that the link hypothesis does not work for the polymorphic system of glycine, and the selective crystallization of α- and β-glycine from water and methanol-water mixture, respectively, is not due to the abundance of glycine aggregates in the solution phase with a similar structure to the crystallizing solid form. The hypothesis of the methanol inhibition on the growth of α-glycine {010} and {010} faces, proposed by Weissbuch (Angew. Chem., Int. Ed. 2005, 44, 3226.), was also studied. The interfaces between the {010} and {010} faces of both crystal forms (α and β) and both solvents (water and methanol-water 3:7 mixture) were studied using molecular simulation. No strong binding of methanol onto the {010} and {010} faces of both crystal forms was observed, and the addition of methanol dilutes the crystal-solvent interactions on all faces. Therefore, the selective crystallization of β and α-glycine with and without methanol does not follow either of the two mechanisms in the literature.

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.