An official website of the United States government
Solution crystallization is a part of the synthesis of materials ranging from geological and biological minerals to pharmaceuticals, fine chemicals, and advanced electronic components. Attempts to predict the structure, growth rates and properties of emerging crystals have been frustrated, in part, by the poor understanding of the correlations between the oligomeric state of the solute, the growth unit, and the crystal symmetry. To explore how a solute monomer or oligomer is selected as the unit that incorporates into kinks and how crystal symmetry impacts this selection, we combine scanning probe microscopy, optical spectroscopy, and all-atom molecular simulations using as examples two organic materials, olanzapine (OZPN) and etioporphyrin I (EtpI). The dominance of dimeric structures in OZPN crystals has spurred speculation that the dimers preform in the solution, where they capture the majority of the solute, and then assemble into crystals. By contrast, EtpI in crystals aligns in parallel stacks of flat EtpI monomers unrelated by point symmetry. Raman and absorption spectroscopies show that solute monomers are the majority solute species in solutions of both compounds. Surprisingly, the kinetics of incorporation of OZPN into kinks is bimolecular, indicating that the growth unit is a solute dimer, a minority solution component. The disconnection between the dominant solute species, the growth unit, and the crystal symmetry is even stronger with EtpI, for which the (010) face grows by incorporating monomers, whereas the growth unit of the (001) face is a dimer. Collectively, the crystallization kinetics results with OZPN and EtpI establish that the structures of the dominant solute species and of the incorporating solute complex do not correlate with the symmetry of the crystal lattice. In a broader context, these findings illuminate the immense complexity of crystallization scenarios that need to be explored on the road to the understanding and control of crystallization.
more »
« less
How to Identify the Crystal Growth Unit
Abstract The structure and composition of the crystal growth unit are of huge fundamental and practical consequence. We propose a method to identify the solute species that incorporates into the growth site on crystal surfaces, the kinks, which rests on the kinetics of the elementary reaction at the kinks. We use as model crystals olanzapine, an antipsychotic medication, and etioporphyrin I, a field‐effect transistor. We combine time‐resolvedin situatomic force microscopy with Raman and absorption spectroscopies, complemented by density functional theory and all‐atom molecular dynamics modeling of the solutions. We show that the structure of the growth unit cannot be deduced neither from the solute oligomers nor from the crystal structure. Chemical kinetics analyses reveal that if the dominant solute species is the one that incorporates into the crystal growth sites, then the kinetics of layer growth complies with a monomolecular rate law. By contrast, if the crystal growth unit assembles from two units of the dominant solute form, a bimolecular rate law ensues. Solutions of both olanzapine and etioporphyrin I are dominated by solute monomers, which exist in equilibrium with a minority of dimers. Whereas numerous olanzapine crystal structures incorporate dimer motifs, etioporphyrin I crystals organize as stacks of monomers. Olanzapine crystal grow by incorporation of dimers. One of the studied face of etioporphyrin I grows by incorporation of the majority monomers, whereas the other one selects the minority dimers as a growth unit. The results highlight the power of the crystallization kinetics analyses to identify the growth unit and illuminate one of the most challenging issues of crystal growth.
more »
« less
- Award ID(s):
- 2128121
- PAR ID:
- 10308314
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Israel Journal of Chemistry
- Volume:
- 61
- Issue:
- 11-12
- ISSN:
- 0021-2148
- Page Range / eLocation ID:
- p. 818-827
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Pharmaceutical solid solutions are gaining increased interest as alternatives to salts and co-crystals for the enhancement of drug solubility and dissolution kinetics. Industrially, they are also responsible for the entrapment of potentially toxic impurities in drug substances. The accidental incorporation of process impurities into the lattice of a growing crystal, or the intentional incorporation of an additive, can vastly alter the product's properties. Reported effects include solubility enhancements, changes in melting point, shifting polymorph stabilities, growth inhibition, and change in crystal habit, among others. This work combines the fields of impurity rejection, solid solutions, and dyeing crystals, to provide visual evidence of those effects, and to further demonstrate how impure regions in a single crystal can present vastly different behaviors to the purified regions of the same crystal. The work revolves around four model host–guest pairs, two of them previously unreported. These include mixed crystals of acetaminophen with curcumin, sulforhodamine B, and acid fuchsin, as well as potassium sulfate dyed with acid fuchsin. Results challenge common assumptions in the study of multicomponent crystals, demonstrating how neglecting composition anisotropy may lead to misdiagnosing solid solutions as surface adsorbed impurities in impurity retention diagnostics, and how neglecting the habit-modifying effects of dissolved impurities may lead to the use of erroneous models for growth inhibition. At the same time, we present opportunities for the development of novel impurity rejection and crystal engineering strategies, aiding the growth of anisotropic crystals with properties that can be fine-tuned in continuum.more » « less
-
Co-crystal engineering is a promising method to create new classes of advanced materials. Co-crystal structure prediction is more challenging when one or more of the lattice constituents (tectons) are flexible molecules. This study reports four co-crystals that were prepared by mixing HAuCl 4 or HAuBr 4 with C 3 -symmetric tectons based on a 1,3,5-(methylacetamide)benzene scaffold. X-ray analysis of the co-crystals revealed the presence of three dominant supramolecular interactions; (a) hydrogen bonding between tecton amide NH residues and the AuX 4 − anion, (b) electrostatic stacking of the Au center against the tecton's π-electrons, (c) very short hydrogen bonds within a proton-bridged-carbonyls motif. Within all four co-crystals, the sterically-geared tecton was trapped in a high energy molecular conformation, which increased the number of favorable intermolecular interactions in the lattice. We infer from the results that the likelihood of high energy molecular conformations within a co-crystal increases if there are multiple dominant intermolecular interactions. Application of this generalizable rule should lead to improved crystal structure prediction.more » « less
-
null (Ed.)Abstract Numerical cloud models require estimates of the vapor growth rate for ice crystals. Current bulk and bin microphysical parameterizations generally assume that vapor growth is diffusion limited, though some parameterizations include the influence of surface attachment kinetics through a constant deposition coefficient. A parameterization for variable deposition coefficients is provided herein. The parameterization is an explicit function of the ambient ice supersaturation and temperature, and an implicit function of crystal dimensions and pressure. The parameterization is valid for variable surface types including growth by dislocations and growth by step nucleation. Deposition coefficients are predicted for the two primary growth directions of crystals, allowing for the evolution of the primary habits. Comparisons with benchmark calculations of instantaneous mass growth indicate that the parameterization is accurate to within a relative error of 1%. Parcel model simulations using Lagrangian microphysics as a benchmark indicate that the bulk parameterization captures the evolution of mass mixing ratio and fall speed with typical relative errors of less than 10%, whereas the average axis lengths can have errors of up to 20%. The bin model produces greater accuracy with relative errors often less that 10%. The deposition coefficient parameterization can be used in any bulk and bin scheme, with low error, if an equivalent volume spherical radius is provided.more » « less
-
High-level quantum chemical computations have provided significant insight into the fundamental physical nature of non-covalent interactions. These studies have focused primarily on gas-phase computations of small van der Waals dimers; however, these interactions frequently take place in complex chemical environments, such as proteins, solutions, or solids. To better understand how the chemical environment affects non-covalent interactions, we have undertaken a quantum chemical study of π– π interactions in an aqueous solution, as exemplified by T-shaped benzene dimers surrounded by 28 or 50 explicit water molecules. We report interaction energies (IEs) using second-order Møller–Plesset perturbation theory, and we apply the intramolecular and functional-group partitioning extensions of symmetry-adapted perturbation theory (ISAPT and F-SAPT, respectively) to analyze how the solvent molecules tune the π– π interactions of the solute. For complexes containing neutral monomers, even 50 explicit waters (constituting a first and partial second solvation shell) change total SAPT IEs between the two solute molecules by only tenths of a kcal mol −1 , while significant changes of up to 3 kcal mol −1 of the electrostatic component are seen for the cationic pyridinium–benzene dimer. This difference between charged and neutral solutes is attributed to large non-additive three-body interactions within solvated ion-containing complexes. Overall, except for charged solutes, our quantum computations indicate that nearby solvent molecules cause very little “tuning” of the direct solute–solute interactions. This indicates that differences in binding energies between the gas phase and solution phase are primarily indirect effects of the competition between solute–solute and solute–solvent interactions.more » « less
