Wenhao Sun

The thermodynamic scale of inorganic crystalline metastability

Data-mining the stability of 29,902 material phases reveals the thermodynamic landscape of inorganic crystalline metastability. The space of metastable materials offers promising new design opportunities for next-generation technological materials, such as complex oxides, semiconductors, pharmaceuticals, steels, and beyond. Although metastable phases are ubiquitous in both nature and technology, only a heuristic understanding of their underlying thermodynamics exists. We report a large-scale data-mining study of the Materials Project, a high-throughput database of density functional theory–calculated energetics of Inorganic Crystal Structure Database structures, to explicitly quantify the thermodynamic scale of metastability for 29,902 observed inorganic crystalline phases. We reveal the influence of chemistry and composition on the accessible thermodynamic range of crystalline metastability for polymorphic and phase-separating compounds, yielding new physical insights that can guide the design of novel metastable materials. We further assert that not all low-energy metastable compounds can necessarily be synthesized, and propose a principle of ‘remnant metastability’—that observable metastable crystalline phases are generally remnants of thermodynamic conditions where they were once the lowest free-energy phase.

Nucleation of metastable aragonite CaCO3 in seawater

Significance Crystallization from solution is a materials synthesis process common both in nature and in the laboratory. Unlike conventional high-temperature solid-state synthesis, solution-based syntheses often yield metastable phases, contrary to expectations from equilibrium thermodynamics. Using a recently developed ab initio scheme to calculate the surface energy of a critical nucleus in equilibrium with the aqueous environment, we present a framework to compare relative nucleation rates between competing polymorphs as a function of solution chemistry. We apply this approach to demonstrate how seawater chemistry can preferentially bias nucleation toward the metastable aragonite phase of calcium carbonate, rather than the stable phase calcite––which is of great relevance to biomineralization, carbon sequestration, paleogeochemistry, and the vulnerability of marine life to ocean acidification. Predicting the conditions in which a compound adopts a metastable structure when it crystallizes out of solution is an unsolved and fundamental problem in materials synthesis, and one which, if understood and harnessed, could enable the rational design of synthesis pathways toward or away from metastable structures. Crystallization of metastable phases is particularly accessible via low-temperature solution-based routes, such as chimie douce and hydrothermal synthesis, but although the chemistry of the solution plays a crucial role in governing which polymorph forms, how it does so is poorly understood. Here, we demonstrate an ab initio technique to quantify thermodynamic parameters of surfaces and bulks in equilibrium with an aqueous environment, enabling the calculation of nucleation barriers of competing polymorphs as a function of solution chemistry, thereby predicting the solution conditions governing polymorph selection. We apply this approach to resolve the long-standing “calcite–aragonite problem”––the observation that calcium carbonate precipitates as the metastable aragonite polymorph in marine environments, rather than the stable phase calcite––which is of tremendous relevance to biomineralization, carbon sequestration, paleogeochemistry, and the vulnerability of marine life to ocean acidification. We identify a direct relationship between the calcite surface energy and solution Mg–Ca ion concentrations, showing that the calcite nucleation barrier surpasses that of metastable aragonite in solutions with Mg:Ca ratios consistent with modern seawater, allowing aragonite to dominate the kinetics of nucleation. Our ability to quantify how solution parameters distinguish between polymorphs marks an important step toward the ab initio prediction of materials synthesis pathways in solution.

Band structure engineering through orbital interaction for enhanced thermoelectric power factor

Band structure engineering for specific electronic or optical properties is essential for the further development of many important technologies including thermoelectrics, optoelectronics, and microelectronics. In this work, we report orbital interaction as a powerful tool to finetune the band structure and the transport properties of charge carriers in bulk crystalline semiconductors. The proposed mechanism of orbital interaction on band structure is demonstrated for IV-VI thermoelectric semiconductors. For IV-VI materials, we find that the convergence of multiple carrier pockets not only displays a strong correlation with the s-p and spin-orbit coupling but also coincides with the enhancement of power factor. Our results suggest a useful path to engineer the band structure and an enticing solid-solution design principle to enhance thermoelectric performance.

Thermodynamics of Phase Selection in MnO2 Framework Structures through Alkali Intercalation and Hydration

While control over crystal structure is one of the primary objectives in crystal growth, the present lack of predictive understanding of the mechanisms driving structure selection precludes the predictive synthesis of polymorphic materials. We address the formation of off-stoichiometric intermediates as one such handle driving polymorph selection in the diverse class of MnO2-framework structures. Specifically, we build on the recent benchmark of the SCAN functional for the ab initio modeling of MnO2 to examine the effect of alkali-insertion, protonation, and hydration to derive the thermodynamic conditions favoring the formation of the most common MnO2 phases-β, γ, R, α, δ, and λ-from aqueous solution. We explain the phase selection trends through the geometric and chemical compatibility of the alkali cations and the available phases, the interaction of water with the system, and the critical role of protons. Our results offer both a quantitative synthesis roadmap for this important class of functional oxides, and a description of the various structural phase transformations that may occur in this system.

Induction time of a polymorphic transformation

When a solution is supersaturated with respect to multiple polymorphs, the polymorph with the lowest nucleation barrier will form first. If this is a metastable polymorph, it will persist until the induction of nuclei of a lower free-energy phase, marking the onset of a polymorphic transformation. Induction of a more stable polymorph can occur under two conditions: a.) during steady-state nucleation, or b.) after complete crystal growth of the metastable phase. Using the theory of competing stochastic processes, we derive the rare probability of forming a higher-barrier, more-stable polymorph during steady-state nucleation, and use this to explain variations in crystallization products between repeat experiments, which may underlie the phenomenon of “Disappearing Polymorphs”. We also derive the induction time of a stable phase from a solution equilibrated with a metastable phase, and show that once the stable phase nucleates, bulk metastable crystals will spontaneously ripen onto nanoscale nuclei or seeds of a more stable phase, thermodynamically driving dissolution–reprecipitation processes. Existing strategies to prolong or shorten the lifetimes of transient metastable phases are reviewed and interpreted within the context of polymorphic induction. The analyses in this work are conducted from classical nucleation and crystal growth theories, suggesting that paradigms shifts to ‘non-classical’ nucleation theories may not be necessary to rationalize multistage crystallization.