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‐resolved
Precrystallization solute assemblies and crystal symmetry
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
- Award ID(s):
- 2128121
- NSF-PAR ID:
- 10347462
- Date Published:
- Journal Name:
- Faraday Discussions
- Volume:
- 235
- ISSN:
- 1359-6640
- Page Range / eLocation ID:
- 307 to 321
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract in situ atomic 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
Abstract Modifiers are commonly used in natural, biological, and synthetic crystallization to tailor the growth of diverse materials. Here, we identify tautomers as a new class of modifiers where the dynamic interconversion between solute and its corresponding tautomer(s) produces native crystal growth inhibitors. The macroscopic and microscopic effects imposed by inhibitor-crystal interactions reveal dual mechanisms of inhibition where tautomer occlusion within crystals that leads to natural bending, tunes elastic modulus, and selectively alters the rate of crystal dissolution. Our study focuses on ammonium urate crystallization and shows that the keto-enol form of urate, which exists as a minor tautomer, is a potent inhibitor that nearly suppresses crystal growth at select solution alkalinity and supersaturation. The generalizability of this phenomenon is demonstrated for two additional tautomers with relevance to biological systems and pharmaceuticals. These findings offer potential routes in crystal engineering to strategically control the mechanical or physicochemical properties of tautomeric materials.
-
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
-
Sequence control in synthetic copolymers remains a tantalizing objective in polymer science due to the influence of sequence on material properties and self-organization. A greater understanding of sequence development throughout the polymerization process will aid the design of simple, generalizable methods to control sequence and tune supramolecular assembly. In previous simulations of solution-based step-growth copolymerizations, we have shown that weak, non-bonding attractions between monomers of the same type can produce a microphase separation among the lengthening nascent oligomers and thereby alter sequence. This work explores the phenomenon further, examining how effective attractive interactions, mediated by a solvent selective for one of the reacting species, impact the development of sequence and the supramolecular assembly in a simple A–B copolymerization. We find that as the effective attractions between monomers increase, an emergent self-organization of the reactants causes a shift in reaction kinetics and sequence development. When the solvent-mediated interactions are selective enough, the simple mixture of A and B monomers oligomerize and self-assemble into structures characteristic of amphiphilic copolymers. The composition and morphology of these structures and the sequences of their chains are sensitive to the relative balance of affinities between the comonomer species. Our results demonstrate the impact of differing A–B monomer–solvent affinities on sequence development in solution-based copolymerizations and are of consequence to the informed design of synthetic methods for sequence controlled amphiphilic copolymers and their aggregates.more » « less
-
DNA-coated colloids can self-assemble into an incredible diversity of crystal structures, but their applications have been limited by poor understanding and control over the crystallization dynamics. To address this challenge, we use microfluidics to quantify the kinetics of DNA-programmed self-assembly along the entire crystallization pathway, from thermally activated nucleation through reaction-limited and diffusion-limited phases of crystal growth. Our detailed measurements of the temperature and concentration dependence of the kinetics at all stages of crystallization provide a stringent test of classical theories of nucleation and growth. After accounting for the finite rolling and sliding rates of micrometer-sized DNA-coated colloids, we show that modified versions of these classical theories predict the absolute nucleation and growth rates with quantitative accuracy. We conclude by applying our model to design and demonstrate protocols for assembling large single crystals with pronounced structural coloration, an essential step in creating next-generation optical metamaterials from colloids.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/D1FD00080B
