David Carriere*,
Jade Raimbault,
Mark A. Levenstein
and
Corinne Chevallard
Université Paris-Saclay, CEA, CNRS, NIMBE, LIONS, 91191 Gif-sur-Yvette, France. E-mail: david.carriere@cea.fr
First published on 18th September 2025
Liquid–liquid phase separation (LLPS) is being increasingly recognized as a critical intermediate step in non-classical crystallization pathways, representing a paradigm shift from early descriptions invoking single-step nucleation. While extensively documented in organic compounds, LLPS in mineral systems presents unique experimental challenges due to accelerated crystallization kinetics, thus complicating physical characterizations (viscosity, chemical composition, etc.). This review synthesizes current evidence for LLPS across diverse mineral systems, examining cases by decreasing order of experimental confidence—from the well-known calcium carbonate system to emerging discussions in oxalates, oxides, metallic nanoparticles, sulfur sols, sulfates, calcium phosphates, metal organic frameworks and sodium chloride. A fundamental challenge lies in definitively establishing liquid character, as cryogenic Transmission Electron Microscopy (TEM) and X-ray scattering methods cannot distinguish between liquid and solid amorphous structures, while liquid-phase TEM observations are prone to interfere with the real crystallization process. Understanding when and why LLPS occurs remains challenging, complicated by inconsistent reporting practices and the predominant use of thermodynamic interpretations where kinetic factors may actually govern the process. Operating far from equilibrium, these systems may require alternative mechanisms beyond classical thermodynamic treatments. Key research frontiers include rigorous demonstration of true liquid character, systematic exploration of structure and dynamics across mineral systems down to the atom and sub-millisecond scales, and integrated experimental–theoretical approaches capturing both thermodynamic and kinetic factors—essential for the rational design of materials and controlled nanoparticle morphologies through LLPS-mediated pathways.
Among the transient states identified in non-classical crystallization pathways (solid-state amorphous particles, clusters, nanocrystals, etc.), dense reactant-rich liquid precursors formed via liquid–liquid phase separation (LLPS) have been particularly discussed. LLPS holds special interest for the crystallization community as it potentially enables the rational design of synthetic crystalline (nano)materials and the understanding of natural phenomena (biomineralization, geological processes), while simultaneously challenging established theoretical frameworks.
While LLPS has been extensively reviewed and documented in organic systems, particularly proteins and molecular assemblies,4,5 where droplets formed by LLPS can develop to macroscopic dimensions, the phenomenon in inorganic or mineral systems has received comparatively less systematic attention. The calcium carbonate system represents the seminal example of mineral LLPS, establishing the foundational framework for understanding liquid–liquid phase separation in inorganic crystallization.6 Beyond this well-documented case, various mineral systems—some discovered even prior to CaCO3—have been reported to exhibit LLPS behavior, expanding the scope of non-classical crystallization pathways across diverse inorganic compositions.
However, direct observation and characterization of LLPS in mineral systems presents significantly greater experimental challenges than compared to those of their organic counterparts. These challenges arise primarily from accelerated crystallization kinetics that drastically limit the temporal window for detection and analysis, often reducing observable lifetimes to milliseconds or seconds rather than the minutes or hours typical in organic systems. Additionally, distinguishing between amorphous solid precursors and true liquid droplets based on dynamic information remains technically demanding, requiring sophisticated experimental approaches compatible with the short lifetime of these transient states. As a consequence, the experimental conditions and mechanisms governing mineral LLPS remain debated, creating a significant knowledge gap in our understanding of non-classical crystallization pathways.
This review aims to address these knowledge gaps by synthesizing existing evidence and identifying common patterns in LLPS behavior across various mineral systems, while critically evaluating the degree of confidence with which LLPS has been demonstrated in different minerals. We describe the systems by decreasing degree of confidence—as we perceive and justify it—examining the proposed mechanisms of LLPS, the transition to more stable mineral phases, and highlighting promising research perspectives for future investigation (Fig. 1).
System | Supporting techniques | Confidence for LLPS |
---|---|---|
Calcium carbonate (no additives) | Liquid-like morphology7–10 (SEM, cryo-TEM), diffusion dynamics11,12 (NMR, MD) | Very high. Debate on condensation of ion pairs vs. pre-nucleation clusters |
Growth dynamics13 (SEM) | ||
Droplet coalescence14 (LP TEM) | ||
Calcium carbonate (with additives) | Liquid-like morphology in bulk, on surfaces and in porous matrices6,15–21 (SEM, cryo-TEM, AFM) | Very high. Polymer-stabilized colloidal liquid |
Diffusion dynamics (NMR11) | ||
Droplet coalescence (LP TEM22) | ||
Cerium oxalate | Liquid-like morphologies in bulk and porous matrices (SEM, cryo-TEM),23–25 droplet coalescence25 (LP TEM) | Very high |
Metal nanoparticles | Liquid-like morphologies in bulk and on surfaces26,27 (cryo-TEM), soft droplets deposited on substrates27 (AFM), liquid-like dynamics28–30 (LP TEM) | Very high. Possible colloidal liquid |
Apatite | Liquid-like morphology in bulk, and in porous matrices31–33 (SEM, cryo-TEM), dense liquid (LP TEM)34,35 | Supportive. Granular structure not systematically assigned to colloidal liquid |
Barium sulfate | Liquid-like morphologies36,37 (TEM) | Suggestive, static images after ethanol quenching and drying under vacuum |
Sulfur hydrosols | Macroscopic behavior of an emulsion38,39 | Suggestive (thermal behavior), lack of sub-micrometer dynamical characterizations |
Oxides | Liquid-like morphologies in bulk40,41 (cryo-TEM) | Tentatively supportive. Colloidal liquid postulated |
Metal organic frameworks | Amorphous phase (electron and X-ray diffraction),42–45 morphology (cryo-TEM)45 | Tentatively supportive |
Either PNCs or LLPS have been reported for crystallization from all the common preparation methods, either in the presence or absence of polymer additives: 1) in situ production of CO2 in a CaCl2 aqueous solution, by reaction of dimethyl carbonate with NaOH,46 2) the Kitano method, where a saturated solution of calcium bicarbonate precipitates CaCO3 upon slow evaporation of water or decrease in CO2 partial pressure,7,8,47,48 3) the ammonia diffusion technique, where thermal decomposition of ammonium carbonate initiates diffusion of carbon dioxide and ammonia into an aqueous CaCl2 solution,6,15,17,20,21,49–55 and 4) direct mixing,9–11,13,14,18,19,22,56–62 where an aqueous solution of calcium chloride is mixed with a (bi)carbonate solution. This final method allows controlling the precursor concentrations, the initial pH of the (bi)carbonate solution, the pH during the whole mineralization process (by titration and simultaneous addition of NaOH/HCl), injection speed, and is easily implemented with various characterization techniques.
Similarly, for simulations,10,12,22,47,57,63–65 it is essential to keep track of the conditions under which the various results were obtained—not only the total simulation time (often on the order of tens of nanoseconds), but also the strategies used to explore the configurational landscape within available computational time. Two strategies include 1) initializing the system at artificially high ion concentrations to enhance sampling efficiency, then transferring to lower concentrations or reconstructing free energy profiles to account for finite-size and concentration effects57,65 and 2) simulating directly dilute conditions with enhanced sampling techniques10,12,63 and/or computationally less demanding methods.47,64
For clarity of presentation, we begin by analyzing reports of LLPS occurring in the absence of organic additives, before turning to the reports in presence of additives, including the so-called PILP (polymer-induced liquid precursor), although PILPs were historically first described.
![]() | ||
Fig. 2 Experimental investigation of LLPS in CaCO3. a and b) Cryo-TEM of a carbonate buffer solution titrated with a CaCl2 solution. Adapted with permission from Smeets et al.10 Copyright 2017 National Academy of Sciences. c and d) Cryo-TEM of carbonate and CaCl2 solutions mixed rapidly, frozen after 100 ms and 2 min after reaction, respectively. Adapted from Rieger et al.9 with permission from The Royal Society of Chemistry. e) TEM of a calcium carbonate solution collected 400 s into the evaporation process (Kitano method). Adapted with permission from Wolf et al.8 Copyright 2011 American Chemical Society. f) LP TEM of carbonate and CaCl2 solutions ca. 500 s after mixing in the observation cell. Adapted with permission from Ramnarain et al.14 Copyright 2022 American Chemical Society. g) SEM of carbonate and CaCl2 solutions mixed gently, 1 min after mixing. Adapted from Rieger et al.9 with permission from The Royal Society of Chemistry. |
To unambiguously assess the liquid character of transient structures, viscosity measurements represent the most direct approach. Bewernitz et al. reported two diffusion coefficients measured by nuclear magnetic resonance (NMR) in calcium chloride/NaOH-titrated carbonate solutions, assigning the slower coefficient (5 × 10−6 cm2 s−1) to liquid structures—only a factor of 2 lower than coexisting (bi)carbonate species in water.11 Molecular dynamics simulations provide detailed insights into liquid properties: computed ion diffusivity ranges from 10−7 to 10−8 cm2 s−1, representing two orders of magnitude higher mobility than in amorphous CaCO3·H2O, yet two to three orders of magnitude lower than ions in solution.10,12,22
More recently, Ramnarain et al. (2022) directly visualized coalescence of these structures using liquid-phase TEM (LP TEM), providing the first unequivocal evidence of their liquid nature (Fig. 2f).14 While the existence of liquid structures is now well-established, direct experimental characterization of their rheological properties remains limited. Nevertheless, the available evidence—spanning diffusion measurements, molecular dynamics predictions, and direct visualization—conclusively validates the initial hypothesis of liquid precursors originally based on indirect morphological observations.
The two competing interpretations to explain the reduction in calcium activity are:
1) sequestration in calcium (bi)carbonate ion pairs. This was reported as compatible with both regular solution theory and classical nucleation theory,10,47 the latter expecting that the concentration of ion pairs exceeds that of larger aggregates by orders of magnitude.
2) sequestration in pre-nucleation clusters (PNCs).11,48,56–60,65,67 Their specific properties have been refined over time, to reach the current status:67 PNCs are defined as disordered nanometer-sized oligomers of several tens of calcium (bi)carbonate moieties; they are reported to form in both under- and supersaturated solutions with respect to the first nucleating species (amorphous calcium carbonate); at fixed supersaturation, their size distribution remains time-independent; as supersaturation increases, the distribution shifts to larger sizes, and the dynamics (mobility inside the PNCs and rate of exchange with coexisting ion pairs) slows down.
The debate between both interpretations is still lively as illustrated by two recent reports on the matter.10,67 In the initial report of Gebauer et al.56 and subsequent work,11,48 oligomerization has been supported by analytical ultracentrifugation (AUC), with sedimentation coefficients compatible with monodisperse 2 nm-sized clusters (≃70 calcium atoms). Authors have, however, acknowledged that the sedimentation was only detected above saturation with respect to the nucleating phase, and that “the clusters cannot be detected in the undersaturated stage of the experiments”, despite the reduction in activity already being evident in this stage. Overlooking the metrological details of size determination, one can also reasonably object, like Smeets et al.,10 that nucleation could happen during centrifugation (>8 h in Gebauer et al.,56 unreported in Bewernitz et al.11). To our knowledge, evidence for or against nucleation 8 h after arrested titration has not been discussed in depth – presumably, owing to the difficulty in mitigating CO2 degassing and maintaining constant supersaturation over such extended periods of time.
Later work has reported direct observation of PNCs by cryo-TEM, showing 0.6 nm sized structures trapped in vitreous ice.48 However, a subsequent analysis of the operation conditions has shown that the achieved resolution was significantly above 1 nm.10 An attempt to observe PNCs with better resolution failed, which is compatible with either the ion pair interpretation, or the loss of PNCs upon grid preparation. In any case, if PNCs exist for CaCO3, their size is significantly smaller than the 2 nm estimated by AUC.
Support to either confirm or disprove PNC formation has been actively pursued through computational approaches. On one hand, Demichelis et al. reported PNCs in molecular dynamics with calcium concentrations in the 60 mM to 500 mM range during 50 to 70 ns, which could be observed within 1 ns after transfer to a 0.4 mM concentration.65 They, however, also reported that the PNCs do not form spontaneously in the 1 ns time window when simulated at 0.4 mM. On the other hand, Smeets et al.10 demonstrated that all initial configurations (random ion distributions and clusters) relax towards ion pairs in the 20–40 mM concentration range and 20–60 ns simulation time, with no significant clustering. Prominent ion pairing was later confirmed by Monte Carlo simulations using mean force potentials and validated experimentally by X-ray absorption spectroscopy.47 However, Gebauer has highlighted a fundamental challenge in bridging density functional theory calculations to molecular mechanics approaches, questioning the seamless integration of different computational scales and their collective ability to accurately represent PNCs behavior in solution.67
Additionally, PNCs formation has been rationalized with a chemical equilibrium model where clusters are produced by attachment /detachment of calcium carbonate ion pairs, each of the equilibria having an identical chemical constant independent of the cluster size,56 compatible with both measured enthalpies and activity reduction.11,56 However, this model alone is not discriminant as it remains compatible with all clustering sizes – including simple ion pairing. On the contrary, Smeets et al. have calculated that standard thermodynamic databases of calcium carbonate are sufficient to account for the reduced calcium activity without invoking clustering other than ion pairing, provided the ion–ion electrostatic interactions are correctly taken into account in the speciation model. They therefore produced a state-of-the-art speciation calculation without invoking additional thermodynamic species like PNCs. It is however not perfectly clear if this rules out the existence of PNCs, as part of the thermodynamic constants used have been constructed by fitting solubility values to experimental data, assuming there are no PNCs.
In summary, the debate surrounding pre-nucleation clusters (PNCs) remains active and will likely require continued experimental developments at synchrotron sources, reassessment of thermodynamic databases, and new theoretical considerations to be fully resolved. Notably, the term “pre-nucleation cluster” has gained widespread acceptance regardless of the preferred mechanistic hypothesis, reflecting the consensus that calcium activity is indeed reduced in these systems. The fundamental disagreement centers not on the sequestration of ions into larger assemblies (potentially reduced to ion pairs), but rather on their thermodynamic interpretation: whether PNCs represent entities with special thermodynamic status distinct from classical theory, or whether they simply describe the sub-critical clusters already encompassed by both classical nucleation theory and regular solution theory.
The critical temperature of the phase separation remains debated, with conflicting evidence supporting both low and high temperature scenarios. A low critical temperature has been proposed based on ACC size evolution with temperature in polymer-containing systems46 and ACC solubility–temperature relationships.59 In contrast, a high-temperature critical point appears more compatible with size variations observed under different supersaturation conditions13 and has been postulated from the retrograde solubility behavior of calcium carbonate.12
![]() | ||
Fig. 3 Experimental investigation of LLPS in CaCO3 in presence of polymers. a) LP TEM of a coalescence event. Adapted with permission from Jin et al.22 Copyright 2024 Springer Nature. b) Cryo-TEM of a solution aliquot after 200–250 min infusion. Adapted with permission from Xu et al.21 Copyright 2018 Springer Nature. c) SEM image of a glass substrate in contact with CaCl2 + polyaspartate solution exposed to ammonium carbonate vapor (ammonia diffusion method, reaction time not recorded). Adapted from Gower et al.,15 with permission from Elsevier. d) Collagen I matrix exposed to a calcium carbonate solution, without polymer, unrecorded reaction time (>3 days) and e) with poly(acrylic acid), after 15 days. Adapted with permission from Olzsta et al.16 Copyright 2003 Springer. |
Polymer addition suppresses direct crystal nucleation, enabling non-crystalline morphologies to develop through highly hydrated intermediates. This polyelectrolyte-induced liquid transient was termed PILP (polymer-induced liquid precursor).15 Various research teams have since reported liquid transients during CaCO3 synthesis using diverse additives: sodium polyacrylate,17,21,55 polycarboxylates,9 double-stranded DNA,21,54 ovalbumin,8,49,53 oligopeptides,52,63 poly(allylamine hydrochloride).21,51 The latter involves a positively charged polyelectrolyte, contrasting with other studies using negatively charged polyelectrolytes that mimic carboxylate-rich biopolymers in biomineralization.
Owing to the difficulty of directly observing transient liquid structures during syntheses, their liquid nature has most often been inferred from characteristic non-crystallographic morphologies (such as films, helices, and cylindrical nanorods),6,15 their ability to infiltrate nanometric pores,16,18 or their selective wetting behavior on various substrates.17,20,55 Berg et al. notably characterized PILP phase wetting on electrospun fibers functionalized with different chemical groups:55 OH-functionalized fibers exhibited full wetting (zero contact angle), while amine-functionalized surfaces induced non-wetting behavior (contact angle >150°). As with additive-free CaCO3 systems, direct confirmation of liquid character emerged only recently, with coalescence events confirmed by Jin et al. in 2024 by LP TEM (Fig. 3a).22
The dense liquid phase was initially interpreted as a coacervate from polymer phase micro-separation with counter-ion calcium (for negatively charged polymers)52,70 or carbonate (for positively charged polymers).51 This was considered a direct consequence of polyelectrolyte addition. However, NMR detection of intermediate liquid phases both with and without polymer additives challenged this view.11 Current understanding shows that polymers primarily stabilize developing liquid phases against amorphous conversion rather than altering underlying thermodynamics or shifting binodal lines.8,58,63 This supports Faatz et al.'s original hypothesis that LLPS can occur in CaCO3 systems without polymers, despite their experiments using additives.46
Several studies have attempted to characterize the composition and viscoelastic properties of the liquid phase.11,19,50 However, characterization remains limited due to time-dependent composition changes from progressive water and polymer ejection.50 The bicarbonate–carbonate composition remains debated: Bewernitz et al. reported bicarbonate enrichment in dense liquid structures relative to reactant-poor solutions,11 while Jin et al. recently detected only bicarbonate at high pH, with composition Ca(HCO3)2·nH2O (n = 7.5 ± 1.7).22 This contrasts with recent NMR reports indicating carbonate-rich phases.61 Similar to cerium oxalate (see below), the viscosity estimated from coalescence events directly recorded by LP TEM is reported to be nine orders of magnitude larger than water.22
Recent observations have questioned the “true” liquid nature of the PILP phase: the absence of intermediate contact angles, gel-like elasticity in AFM measurements,19 and nanogranularity in all PILP-formed structures (Fig. 3b).21,68 Cryo-TEM and NMR studies by Sommerdijk and coworkers on various polymer additives revealed no liquid phase but the systematic presence of ca. 50 nm amorphous clusters composed of 2 nm nanoparticles. They concluded that “PILP is actually a polymer-driven assembly of ACC clusters, and [that] its liquid-like behavior at the macroscopic level is due to the small size and surface properties of the assemblies”. The authors suggested redefining PILP as “polymer-induced liquid-like precursor”. Similarly, Gower coined the process “colloid assembly and transformation (CAT)” in 2017.68 Recent LP TEM studies with non-polymer additives confirm structural differences between polymer-free and polymer-stabilized LLPS:14,21 the former generates truly homogeneous liquid droplets of ion complexes, while the latter produces polymer-stabilized nanoparticle assemblies that change dynamically in shape.
In summary, like for additive-free calcium carbonate crystallization, LLPS propositions based on indirect observations have been validated by recent LP TEM characterization. However, the process involves colloidal assembly rather than a dense molecular liquid, though the nature of this shift in the presence of polymer remains unclear and certainly calls for active research.
![]() | ||
Fig. 4 Experimental investigation of LLPS in cerium oxalate. a and b) SEM of deposits collected on top of a 0.2 μm pore filter membrane. Repoduced from Rodríguez-Ruiz et al.23 with permission from The Royal Society of Chemistry. c) Cryo-TEM of a freeze-quenched solution 7 s after reaction. Adapted with permission from Raimbault et al.25 Copyright 2024 American Chemical Society. Identical image sizes. d) LP TEM of a coalescence event between two droplets wetting the observation window. Adapted with permission from Raimbault et al.25 copyright 2024 American Chemical Society. The two images are separated by ca. 1.2 s. |
![]() | ||
Fig. 5 Experimental investigation of LLPS in syntheses of metal nanoparticles. a) LP TEM of a HAuCl4 aqueous solution. Blue arrow: gold-poor liquid. Orange arrows: gold-rich structures. Adapted with permission from Loh et al.28 Copyright 2017 Springer Nature. b and c) LP TEM of a HAuCl4 aqueous solution. The arrow indicates the structure assigned to a dense liquid corona made of nanometer-sized clusters. Adapted with permission from Jin et al.29 Copyright 2019 Wiley-VCH. d) Liquid-cell AFM of mixture of aqueous HAuCl4 and sodium tricitrate in contact with a graphite substrate, imaged 150 min after mixing. Adapted with permission from Mikhlin et al.27 Copyright 2011 Elsevier. e) Same by TEM f and g) LP TEM of a Na2PdCl4 aqueous solution. The red arrow indicates a darker structure assigned to a crystal nucleus. The inset shows nanometer-scaled structures assigned to ultrasmall clusters constituting the liquidlike ‘cloud’. Adapted with permission from Jin et al.29 Copyright 2019 Wiley-VCH. h and i) LP TEM of a solution of palladium(II) 2,4-pentanedionate in methyl pyrrolidone. Red dashed lines outline the dense liquid structure around the nanocrystal. Adapted with permission from Li et al.30 Copyright 2021 American Chemical Society. |
Transmission electron microscopy cannot investigate additive-induced gold reduction syntheses because the electron beam dominates the reduction process. However, liquid-phase electron microscopy has become an extremely convenient in situ characterization tool using the beam as the electron source. Loh et al. reduced tetrachloroaurate by a TEM electron beam, and observed the formation of liquid structures of tens of-nanometers, rich in hydrated gold complexes with an estimated 4 M concentration (Fig. 5a).28 Jin et al. later observed similar structures by liquid-phase TEM, and additionally reported nanocrystals embedded in the parent Au-rich liquid phase (Fig. 5b and c).29
Liquid–liquid phase separation has been reported for other metal nanoparticles including silver28 and platinum.30 In palladium (Fig. 5f–i), additional substructure was resolved within the liquid: the liquid phase was characterized as a “cluster cloud”—a disordered, dynamic assembly of nanometer-sized amorphous cluster (Fig. 5f and g).29 For nickel nanoparticles, Yang et al. describe amorphous particles with “randomly packed granular features” resembling LLPS with cluster clouds, although they make no explicit LLPS connection.76
The generality of LLPS in metal ion solutions upon reduction both induced chemically and by liquid-phase TEM makes little doubt. However, challenges remain regarding the appropriate mechanistic description. Loh et al. notoriously proposed that the gold-rich liquid forms by spinodal decomposition. However, the lack of quantitative data prevented validation through, for example, the expected structure factors. As they acknowledged, confinement in LP TEM cells slows diffusion processes by orders of magnitude, potentially changing the rate-limiting step compared to bulk syntheses. Regarding the presence of the ‘cluster cloud’ outside a TEM, a positive indication was given by Ramamoorthy et al.77 They used Small-Angle X-ray Scattering (SAXS) to evidence aggregates containing 15–70 Au atoms in organic gold precursors before reduction (tetrachloroaurate stabilized by oleylamine in hexane). They described these nanometer-scale structures as pre-nucleation clusters, referencing well-known CaCO3 observations. Although they did not report possible aggregation into liquid phases, this may indicate a common driving force between PNCs at equilibrium prior to reaction and the amorphous primary units of the “cluster cloud”, but more direct evidence with simultaneous validation by TEM would be beneficial.
![]() | ||
Fig. 6 Snapshots from a molecular dynamics simulation of a NaCl solution in the spinodal domain (18.5 mol kg−1). a) after 10 ns, and b) after 70 ns. Ions with less than 8 neighbors are omitted. Orange particles show crystalline domains. Reprinted with permission from Jiang et al.78 Copyright 2019 American Chemical Society. |
![]() | ||
Fig. 7 Experimental investigation of LLPS in phosphate phases. a) Cryo-TEM of a crystallizing calcium phosphate solution after 10 min of reaction. Adapted with permission from Nudelman et al.32 Copyright 2010 Springer Nature. The limited resolution stems from a magnification of the original image, to ensure consistency with the following images. b) SEM of a collagen matrix impregnated for 4 days with a calcium phosphate solution. Adapted with permission from Jee et al.31 Copyright 2010 Elsevier. The calcium phosphate particles are interpreted as an inefficient impregnation in the absence of LLPS. c) Cryo-TEM of a crystallizing calcium phosphate solution in presence of polyaspartate after 10 min of reaction. Adapted with permission from Nudelman et al.32 Copyright 2010 Springer Nature. d) SEM of a collagen matrix impregnated 3 days with a calcium phosphate solution in the presence of polyaspartate. The absence of particles is interpreted as a better wetting of a PILP of calcium phosphate. Adapted with permission from Jee et al.31 Copyright 2010 Elsevier. e) LP TEM of a CaP solution, unspecified reaction time, evidencing a less electron-dense corona around particles. Adapted with permission from Dalmônico et al.35 Copyright 2022 American Chemical Society. f) Cryo-TEM and diffraction of a crystallizing struvite solution after 5 s of reaction, from Karafiludis et al.33 Copyright 2023 American Chemical Society. |
In polymer-free systems, early studies reported similar nanocluster assemblies (Fig. 7a), although neither Dey et al.85 nor Habraken et al.86 explicitly invoked liquid–liquid phase separation. Instead, these authors debated whether the observed nanoclusters represented ‘pre-nucleation clusters’ following the Betts and Posner model,87 or ion association complexes. Because of their corrugated morphologies, the shapes observed in these systems are not typically reminiscent of liquid phases – unless observed at coarser resolution, as is typically the case in LP TEM. Similarly, collagen matrices impregnated with CaP without polymers show surface deposits of mineral nodules, presented as an evidence for the lack of liquid–liquid phase separation (Fig. 7b). Two recent LP TEM observations finally claimed evidence for LLPS: He et al.34 evidenced a reactant-rich liquid as dark-contrast zones from which the formation of amorphous and crystalline apatite occurs, and Dalmônico et al.35 additionally reported a liquid phase forming mainly onto the observation windows. However, in our opinion, the out-of-equilibrium morphologies and the lack of shape fluctuation are compatible with solid-state nanocluster assemblies (Fig. 7e) – or at least, extremely low capillary velocities. Critically, no dynamic behaviors characteristic of liquids—such as coalescence or relaxation to equilibrium shapes—have been documented, leaving the actual viscosity of these liquid states uncharacterized.
CaP formation with and without additives involves aggregated nanometer-sized structures. At least with polymers, sufficient indirect evidence indicates that separating structures are ‘liquid-like’. However, whether such hybrid assemblies can be properly described using liquid–liquid phase separation diagrams typically applied to molecular liquid systems remains undetermined.
Recent advances have extended beyond calcium-based systems to other phosphate-based phases. Karafiludis et al.33 reported the formation of struvite (MgNH4PO4·6H2O) through a liquid precursor pathway, observing the formation of dense liquid particles in cryo-TEM (Fig. 7f). The liquid character of these precursor phases was demonstrated through the observation of coalescence morphologies. Notably, the SAXS patterns presented in their study appear to rule out the possibility that these liquid phases consist of discrete nanometer-sized clusters, but instead support nanometer-scale concentration fluctuations (ca. 1 nm) that preceded the formation of the dense liquid.
![]() | ||
Fig. 8 Experimental investigation of LLPS in barium sulfate. a) TEM of a solution quenched by ethanol addition. Reproduced from Ruiz-Agudo.36 Copyright 2015 American Chemical Society. b) STEM-HAADF of a solution quenched by ethanol addition. Reproduced from Ruiz-Agudo et al.37 with permission from The Royal Society of Chemistry. |
A consensus exists regarding the liquid nature of the transient sulfur particles and their subsequent transition to crystalline phases. However, experimental characterization remains limited: literature reports focus primarily on the bulk suspension viscosity (reaching “honey-like viscosity”39 at concentrations of 600 g L−1) rather than the intrinsic properties of the liquid forming the individual particles. The liquid character is primarily assessed through thermal behavior analysis of concentrated sols at sulfur's melting point, where the absence of abrupt viscosity changes is interpreted as evidence that sulfur exists in a liquid state. This interpretation requires careful consideration, as concentrated colloidal suspensions of solid-state nanoparticles (such as SiO2) can exhibit macroscopically liquid-like behavior.89
The mechanism governing the transition from liquid droplets to crystals remains a subject of ongoing debate. LaMer and Dinegar originally postulated what would now be termed classical (direct) nucleation from dissolved monomers, accompanied by secondary nucleation through contact between crystal nuclei and metastable droplets. In contrast, Steudel proposes a multistep nucleation mechanism, reasoning that “the unit cell consists of not less than 16 molecules which have to be oriented in a special manner. This ordering process is hindered first by the impurities present in the dirty cluster and second by the high viscosity liquid sulfur must exhibit at 20 °C”.90 Regarding the formation of the liquid droplets themselves, Lamer has established this system as a paradigm for nucleation, supported by the critical (in the physical sense) nature of the transition: the white milk of sulfur forms abruptly above a given concentration. However, a modern analysis by Whitehead et al. has since objected that the large colloids observed by LaMer are far from critical sizes typically expected, and that the quantitative evolution of turbidity deviates from a classical description of nucleation.88
In summary, sulfur represents a widely accepted paradigm of liquid-mediated nucleation in mineral systems, with the liquid nature of the intermediate phase being generally acknowledged within the crystallization community. However, despite this widespread acceptance, to our knowledge, no modern in situ characterization techniques have been systematically applied to definitively confirm the liquid properties of sulfur droplets or to elucidate the detailed transition mechanisms at the molecular level.
H2SO4 + 2Na[Al(OH)4] → 2AlO(OH) + Na2SO4 + 4H2O |
![]() | ||
Fig. 9 Experimental investigation of LLPS in oxides a) cryo-TEM prior to boehmite crystallization, freezed ca. 10 ms after reaction. Adapted with permission from Haberkorn et al.40 Copyright 2003 Elsevier. Bubbling is induced by the electron beam to demonstrate the difference between the dense liquid precursor and solid-state amorphous particles. b) Cryo-TEM prior to magnetite crystallization, aliquot of a dilute HCl solution titrated with a FeCl3 solution. Adapted with permission from Scheck et al.41 Copyright 2016 American Chemical Society. The granular structure is assigned to pre-nucleation clusters. |
Regarding the formation of the liquid structures themselves, Haberkorn et al. posit that mixing heterogeneity leads to the formation of reagent droplets at the output of the mixer nozzle. Indeed, both reagents are mixed thought a simple T-junction, and the authors raised concerns about the solution reaching molecular-scaled homogeneity before cryogenic quenching. Back in 2003, the transient liquid structures were therefore not assigned to an underlying thermodynamic phase diagram, but rather to a kinetic hydrodynamic origin – a position overlooked nowadays.
For magnetite formation via iron salt hydrolysis, Kashyap et al. reported a “liquid-like, disordered pre-nucleation phase” prior to crystallization in the presence of Mms6 protein.92 While they provided evidence for an amorphous layer forming at vesicle surfaces, the actual liquid character remains unconfirmed. Subsequently, Scheck et al. observed aggregated nanometer-sized structures interpreted as pre-nucleation clusters in additive-free syntheses (Fig. 9b).93 Drawing on controversial calcium carbonate developments, they cautiously proposed that interface formation “may be conceived as passing a liquid–liquid binodal”. Building on this hypothesis, Kuhrts et al. rationalized organic additive effects on magnetite nanoparticle size by invoking the wetting properties of transient liquid-like ferrihydrite:94 “approximating the ferrihydrite precursor as a liquid, [Scheck 2016]93 we treat its interaction with magnetite as a wetting process”. Although magnetite formation via non-crystalline particles is well-established,95,96 the liquid-state nature of these transient intermediates remains hypothetical and requires direct experimental validation.
There is no doubt that ZIF-8 formation follows a non-classical pathway via amorphous intermediates. Cravillon et al. first reported non-crystalline transient phases in ZIF-8 synthesis using time-resolved Small- and Wide Angle X-ray Scattering (SAXS/WAXS), observing the immediate formation of ca. 1 nm clusters followed by the growth of larger nanoscale particles that crystallized within 15 s, providing strong evidence for the initial precipitation of an amorphous precursor phase.42 Subsequent studies by Terban et al.43 and Carraro et al.44 confirmed these results, unequivocally demonstrating a long-lived amorphous (implicitly, solid) phase prior to ZIF-8 crystallization, with increased stability of this transient phase in water compared to methanol.
In this context, Liu et al.45 first suggested a liquid precursor phase in ZIF-8 synthesis based on liquid-cell and cryo-TEM observations of aqueous synthesis (Fig. 10). However, while their cryo-TEM images show “liquid-like” structures, their actual liquid character remains undemonstrated: liquid-cell TEM videos lack temporal and spatial resolution and show no clear structural progression. Additionally, they provided no quantification of density, chemistry, or putative fluid properties, and few subsequent reports corroborate their findings. Jongert et al. reported LLPS based on in situ electrochemical measurements of particles with positive z-potential during the first two minutes of methanolic syntheses,97 yet these results do not unambiguously confirm a transient liquid phase. Another cryo-TEM study by Ogata et al.100 investigated aqueous ZIF-8 synthesis at varying imidazole:
Zn ratios with and without bovine serum albumin, observing various amorphous phases leading to crystals through dissolution–reprecipitation and solid-state transformations, but reported no liquid precursor phase.
![]() | ||
Fig. 10 Experimental investigation of LLPS in ZIF-8. Cryo-TEM of the solution frozen a) 10 s and b) 60 s after mixing, respectively. Adapted with permission from Liu et al.45 Copyright 2021 National Academy of Sciences. The scale bar applies to both images. |
Dok et al. used in situ harmonic light scattering and NMR to follow the nucleation of ZIF-8 crystals in methanol,101 demonstrating formation of positively charged clusters ‘instantaneously’ upon mixing. They interpreted these results as PNCs that release protons and condense into amorphous particles before crystallizing. This supports earlier claims of amorphous transient phases; however, more work is required to reconcile competing theories regarding condensed phase formation through PNCs versus LLPS. Recent MOF and ZIF-8 nucleation modeling has examined non-classical mechanisms but has not evidenced or considered LLPS, instead focusing on cluster formation, including disordered “glassy clusters”, dependent on solvent and kinetic factors.102,103 Therefore, evidence for non-classical nucleation pathways in ZIF-8 crystallization is strong, but more data are needed to confirm the involvement of liquid–liquid phase separation.
First, characterizing not only the presence but the emergence itself of transient structures (liquid or not) remains a huge experimental challenge (Fig. 1a), motivated by questions regarding LLPS mechanisms via ion pairs, PNCs, nucleation, spinodal decomposition, etc. Although nanometer-scale resolution is technically accessible today, time resolution remains limited. Reaction times in LP-TEM are ill-defined because of overly simple mixers in state-of-the-art liquid sample holders, creating a fundamental bottleneck due to supersaturation drops when the parent solution starts evolving. Although millisecond resolution has been reached in SAXS, some reports show this is still insufficient.24,25,74,96,104 A first perspective in the field is to push the boundaries of time resolution, presumably with advanced microfluidic tools.105
Second, definitively establishing the ‘liquid’ character of observed precursor phases, in terms of physical dynamics (e.g., viscosity), remains a challenge (Fig. 1b). While state-of-the-art cryo-TEM and SAXS/WAXS analyses ensure detection down to a few milliseconds after reaction, these static techniques are blind to visco-elastic properties and therefore unsuitable for distinguishing between liquid and solid amorphous structures. Past intuitions and postulates based on macroscopic and indirect evidence have proven correct in several emblematic cases, suggesting that simpler analytical approaches should not be overlooked in the pursuit of more sophisticated characterization methods. These simple approaches should nowadays be systematically validated with AFM and LP TEM, which enable testing directly true liquid behavior, including droplet coalescence. However, caution should be taken because these techniques are only able to access viscosities starting around five orders of magnitude higher than water. Thus, the typical ‘liquid-like’ displacements accessible by electron microscopy (nanometer scale) over experimental time resolutions (seconds, up to several days in impregnation experiments) correspond to capillary velocities similar to those of tar (nm s−1), and would be classified as solid or at least “gel-like” at typical laboratory scales. The “liquid-like” behavior even observed directly may actually reflect extremely viscous flow over nanometer distances and minute timescales rather than true liquid dynamics. This interpretation challenge is compounded by the already discussed issues of representativeness in LP TEM observations, including electron beam effects and confinement artifacts. A second perspective in the field is to use systematically methods (AFM, LP-TEM) that distinguish solid-state and true liquid amorphous phases, with the even more challenging prospect of systematically measuring viscosity. Development of alternative techniques less prone to electron beam damage or interference from SiN membrane surfaces is strongly needed.
Further, while the term “liquid” usually refers to molecular or ionic liquids, many have reported LLPS involving aggregated primary nanoparticles (such as in PILP systems and cluster clouds, and in some of the debated PNC pathways). In these colloidal assemblies, the bulk phase may exhibit flow behavior, while individual primary particles might remain in solid-state configurations. This distinction between macroscopic flow properties and microscopic structural states represents a critical aspect of accurate phase characterization. A third perspective in the field is therefore to systematically assess internal structure combined with local dynamics (Fig. 1c), for instance coupling SAXS/WAXS and NMR.
Finally, validating an appropriate theoretical framework represents the last paramount challenge (Fig. 1d). The complexity of this understanding is exacerbated by inconsistent reporting practices, which raises questions about whether LLPS is genuinely absent in unreported cases or simply falls beyond the scope of individual studies. More systematic exploration across different mineral systems is needed as a prerequisite to developing a better phenomenological understanding of LLPS occurrence and, ultimately, comprehensive theoretical frameworks. The dominant consensus describes LLPS through thermodynamic concepts (phase diagrams, equilibrium constants), despite the parent solution being extremely far from equilibrium. In electrolyte systems, this thermodynamic approach aligns with theoretical predictions, particularly for primitive model systems where LLPS is predicted in +z/−z electrolytes with z > 3.8 in water using Gibbs ensemble methods.106 However, alternative mechanisms merit consideration. In the field of thermodynamic models, Stawski et al. summarized how the emergence of (solid-state) amorphous transient structures can be rationalized, relying on a balance between either the solubility or interfacial tension of the amorphous phases and the crystal,107 which are applicable to liquid structures. However, some evidence suggests that kinetic rather than purely thermodynamic factors may govern LLPS processes. Reinvestigation of alternative pathways similar to mixing heterogeneity, irreversible condensation or out-of-equilibrium fluctuations could prove fruitful. Since these systems operate far from equilibrium conditions, the traditional Gibbs phase rule may not apply, potentially allowing for multiple competing pathways and more complex behavior than classical thermodynamic treatments would predict. A final perspective in the field is therefore to integrate experimental and theoretical approaches that can capture both out-of-equilibrium thermodynamics and kinetic factors governing LLPS.
This journal is © The Royal Society of Chemistry 2025 |