skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Selective desolvation in two-step nucleation mechanism steers crystal structure formation
The two-step nucleation (TSN) theory and crystal structure prediction (CSP) techniques are two disjointed yet popular methods to predict nucleation rate and crystal structure, respectively. The TSN theory is a well-established mechanism to describe the nucleation of a wide range of crystalline materials in different solvents. However, it has never been expanded to predict the crystal structure or polymorphism. On the contrary, the existing CSP techniques only empirically account for the solvent effects. As a result, the TSN theory and CSP techniques continue to evolve as separate methods to predict two essential attributes of nucleation – rate and structure. Here we bridge this gap and show for the first time how a crystal structure is formed within the framework of TSN theory. A sequential desolvation mechanism is proposed in TSN, where the first step involves partial desolvation to form dense clusters followed by selective desolvation of functional groups directing the formation of crystal structure. We investigate the effect of the specific interaction on the degree of solvation around different functional groups of glutamic acid molecules using molecular simulations. The simulated energy landscape and activation barriers at increasing supersaturations suggest sequential and selective desolvation. We validate computationally and experimentally that the crystal structure formation and polymorph selection are due to a previously unrecognized consequence of supersaturation-driven asymmetric desolvation of molecules.  more » « less
Award ID(s):
2132022
PAR ID:
10351271
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Nanoscale
Volume:
14
Issue:
5
ISSN:
2040-3364
Page Range / eLocation ID:
1723 to 1732
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; (, Chemical Engineering Research and Design)
    Heng, Jerry (Ed.)
    The morphological evolution of organic crystals during crystallization depends on the face-specific growth rates. Classical growth rate models relate the face-specific growth rates to the crystal lattice, energy of stable facets, growth mechanism, and supersaturation. The complexities of these models have increased over time to account accurately for solution conditions, the structure of growth units, and their attachment rates. Such advanced growth rate models require several layers of computations to obtain attachment energies of facets, nucleation rates, kink density, and attachment rates. Among these, the most intensive and time-consuming computation is for attachment rates, which require molecular dynamic simulations. This substantially increases the overall computation time to predict the absolute growth rate for even one crystallization condition. Since it is nearly impossible to iterate such a growth rate model, optimization schemes cannot be implemented to identify solution conditions that favor specific crystal growth. To reduce the computational time for attachment rate calculations, we implement a group contribution method (GCM) that relates the properties of functional groups in a molecule to their attachment rates to the crystal lattice, thereby rapidly estimating the growth rates of organic crystals. The process of molecular attachment involves partial desolvation of a solvated molecule, referred to as a transition state, followed by total desolvation via spontaneous attachment to a crystal facet. The first step in GCM is to identify the equilibrium states of fully solvated and partially desolvated solute molecules. The degree of supersaturation dictates the extent of this equilibrium and, thereby, the activation barrier for the growth of crystals, according to transition state theory. Identifying this equilibrium phenomenon allows for capturing the functional-group-specific interactions that depend on molecular motion, which could be related to operating conditions such as temperature and pressure. The stochastic optimization technique with Monte-Carlo sampling allows an efficient optimization problem solution to obtain the group interaction parameters. The GCM approach is first validated for the estimation of growth rates of glutamic acid and L-histidine, and then extended to predict growth rates of alanine and glycine rapidly. The optimized parameters and GCM scheme can be used to estimate growth rates in other crystallization systems. 
    more » « less
  2. ; (, Nature Communications)
    Abstract An inexpensive and reliable method for molecular crystal structure predictions (CSPs) has been developed. The new CSP protocol starts from a two-dimensional graph of crystal’s monomer(s) and utilizes no experimental information. Using results of quantum mechanical calculations for molecular dimers, an accurate two-body, rigid-monomer ab initio-based force field (aiFF) for the crystal is developed. Since CSPs with aiFFs are essentially as expensive as with empirical FFs, tens of thousands of plausible polymorphs generated by the crystal packing procedures can be optimized. Here we show the robustness of this protocol which found the experimental crystal within the 20 most stable predicted polymorphs for each of the 15 investigated molecules. The ranking was further refined by performing periodic density-functional theory (DFT) plus dispersion correction (pDFT+D) calculations for these 20 top-ranked polymorphs, resulting in the experimental crystal ranked as number one for all the systems studied (and the second polymorph, if known, ranked in the top few). Alternatively, the polymorphs generated can be used to improve aiFFs, which also leads to rank one predictions. The proposed CSP protocol should result in aiFFs replacing empirical FFs in CSP research. 
    more » « less
  3. ; ; ; ; ; ; ; ; ; ; et al (, Acta Crystallographica Section B Structural Science, Crystal Engineering and Materials)
    A seventh blind test of crystal structure prediction was organized by the Cambridge Crystallographic Data Centre featuring seven target systems of varying complexity: a silicon and iodine-containing molecule, a copper coordination complex, a near-rigid molecule, a cocrystal, a polymorphic small agrochemical, a highly flexible polymorphic drug candidate, and a polymorphic morpholine salt. In this first of two parts focusing on structure generation methods, many crystal structure prediction (CSP) methods performed well for the small but flexible agrochemical compound, successfully reproducing the experimentally observed crystal structures, while few groups were successful for the systems of higher complexity. A powder X-ray diffraction (PXRD) assisted exercise demonstrated the use of CSP in successfully determining a crystal structure from a low-quality PXRD pattern. The use of CSP in the prediction of likely cocrystal stoichiometry was also explored, demonstrating multiple possible approaches. Crystallographic disorder emerged as an important theme throughout the test as both a challenge for analysis and a major achievement where two groups blindly predicted the existence of disorder for the first time. Additionally, large-scale comparisons of the sets of predicted crystal structures also showed that some methods yield sets that largely contain the same crystal structures. 
    more » « less
  4. ; ; ; ; (, ChemRxiv)
    Polymorphism in molecular crystals influences their properties and performance. Crystal structure prediction (CSP) can help explore the crystal structure landscape and discover potentially stable polymorphs computationally. We present a new version of the Genarris open-source code, which generates random molecular crystal structures in all space groups and applies physical constraints on intermolecular distances. The main new feature in Genarris 3.0 is the ``Rigid Press algorithm, which uses a regularized hard-sphere potential to compress the unit cell and achieve a maximally close-packed structure based on purely geometric considerations without performing any energy evaluations. In addition, Genarris 3.0 is interfaced with machine-learned interatomic potentials (MLIPs) to accelerate the exploration of the potential energy landscape. We present a new clustering and down-selection workflow that employs the MACE-OFF23(L) MLIPs to perform geometry optimization and energy ranking in the early stages. We use Genarris 3.0 to successfully predict the structure of six targets: aspirin, Target I and Target XXII from previous CSP blind tests, and the energetic materials HMX, CL-20, and DNI. We further analyze the performance of MACE-OFF23(L) compared to dispersion-inclusive density functional theory (DFT) for geometry relaxation and energy ranking. We find significant variability in the performance of MACE-OFF23(L) across chemically diverse targets with particularly poor performance for energetic materials, which is mitigated by our clustering and down-selection procedure. Genarris 3.0 can thus be used effectively to perform CSP and to generate molecular crystal datasets for training ML models. 
    more » « less
  5. ; ; ; ; (, Proceedings of the National Academy of Sciences)
    The Earth's inner core started forming when molten iron cooled below the melting point. However, the nucleation mechanism, which is a necessary step of crystallization, has not been well understood. Recent studies have found that it requires an unrealistic degree of undercooling to nucleate the stable, hexagonal, close-packed (hcp) phase of iron that is unlikely to be reached under core conditions and age. This contradiction is referred to as the inner core nucleation paradox. Using a persistent embryo method and molecular dynamics simulations, we demonstrate that the metastable, body-centered, cubic (bcc) phase of iron has a much higher nucleation rate than does the hcp phase under inner core conditions. Thus, the bcc nucleation is likely to be the first step of inner core formation, instead of direct nucleation of the hcp phase. This mechanism reduces the required undercooling of iron nucleation, which provides a key factor in solving the inner core nucleation paradox. The two-step nucleation scenario of the inner core also opens an avenue for understanding the structure and anisotropy of the present inner core. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email