Classical density functional theory, unconstrained crystallization, and polymorphic behavior

Diffusive first-order phase transition: nucleation, growth and coarsening in solids

The phenomena of nucleation and growth, which fall into the category of first-order phase transitions, are of great importance. They are present everywhere in our daily lives. They enable us to understand and model a vast number of phenomena, from the formation of raindrops, to the gelling of polymers, the evolution of a virus population and the formation of galaxies. Surprisingly, this whole range of phenomena can be described by two seemingly antagonistic approaches: classical nucleation theory, which highlights the atomistic approach of the diffusion process, and the phase-field (PF) approach, which erases the discrete nature of the diffusion process. Although there is an huge quantity of articles and review papers dealing with the problem of first-order phase transition, the subject is so important and vast that it is very difficult to provide nowadays exhaustive syntheses on the subject. The revival over the past 20 years in the condensed matter world of PF approaches such as PF crystal, or the recent development of optimization methods such as gentle ascend dynamics, as well as the emergence of atom probe tomography, have enabled us to better understand the links between these antagonistic approaches, and above all to provide new experimental results to test the limits of both. This renewal has motivated the writing of this review, both to take stock of current knowledge on these two approaches. This review has two distinct objectives: summarizing generic previous models applies to discuss the nucleation, the growth and the coarsening processes. Despite some reviews already exist on these different subject, few of them present the different logical links between these models and their limitations, unifying them within the framework of the theory of macroscopic fluctuations, which has been developed over the last 20 years. In particular, we present the extension of the Cahn-Hilliard formalism to model the nucleation and growth process and we discuss the relevance of the notion of pseudo-spinodal and discuss. Such an extension allows interpreting experiments performed fat from the solubility limit and the spinodal line. Finally, this work proposes some clues to make this unified approach more predictive.

Statistical mechanics of crystal nuclei of hard spheres

In the study of crystal nucleation via computer simulations, hard spheres are arguably the most extensively explored model system. Nonetheless, even in this simple model system, the complex thermodynamics of crystal nuclei can sometimes give rise to counterintuitive results, such as the recent observation that the pressure inside a critical nucleus is lower than that of the surrounding fluid, seemingly clashing with the strictly positive Young-Laplace pressure we would expect in liquid droplets. Here, we re-derive many of the founding equations associated with crystal nucleation and use the hard-sphere model to demonstrate how they give rise to this negative pressure difference. We exploit the fact that, in the canonical ensemble, a nucleus can be in a (meta)stable equilibrium with the fluid and measure the surface stress for both flat and curved interfaces. Additionally, we explain the effect of defects on the chemical potential inside the crystal nucleus. Finally, we present a simple, fitted thermodynamic model to capture the properties of the nucleus, including the work required to form critical nuclei.

A microscopic approach to crystallization: Challenging the classical/non-classical dichotomy

We present a fundamental framework for the study of crystallization based on a combination of classical density functional theory and fluctuating hydrodynamics that is free of any assumptions regarding order parameters and that requires no input other than molecular interaction potentials. We use it to study the nucleation of both droplets and crystalline solids from a low-concentration solution of colloidal particles using two different interaction potentials. We find that the nucleation pathways of both droplets and crystals are remarkably similar at the early stages of nucleation until they diverge due to a rapid ordering along the solid pathways in line with the paradigm of "non-classical" crystallization. We compute the unstable modes at the critical clusters and find that despite the non-classical nature of solid nucleation, the size of the nucleating clusters remains the principle order parameter in all cases, supporting a "classical" description of the dynamics of crystallization. We show that nucleation rates can be extracted from our formalism in a systematic way. Our results suggest that in some cases, despite the non-classical nature of the nucleation pathways, classical nucleation theory can give reasonable results for solids but that there are circumstances where it may fail. This contributes a nuanced perspective to recent experimental and simulation work, suggesting that important aspects of crystal nucleation can be described within a classical framework.

Confinement Effects in Droplet Formation on a Solid Particle

Formation of a droplet around a spherical solid particle in supersaturated vapor is considered. The number and stability of equilibrium solutions in a closed small system are studied in the canonical ensemble in comparison to an open system in the grand canonical ensemble. Depending on the system's parameters, two modes exist in the canonical ensemble: the first one with only one solution and the second one with three solutions; the presence of the third solution is due to confinement. The analysis is conducted first on a macroscopic thermodynamic level of description, and then the results are supported by studies within two versions of classical density functional theory: the square-gradient approximation with the Carnahan-Starling equation of state for hard spheres on a completely wettable particle and the random-phase approximation with the fundamental measure theory on a poorly wettable particle. In the latter case, a solution breaking the spherical symmetry is observed at a small total number of molecules.

Effect of Laser-Exposed Volume and Irradiation Position on Nonphotochemical Laser-Induced Nucleation of Potassium Chloride Solutions

Herein, we study the influences of the laser-exposed volume and the irradiation position on the nonphotochemical laser-induced nucleation (NPLIN) of supersaturated potassium chloride solutions in water. The effect of the exposed volume on the NPLIN probability was studied by exposing distinct milliliter-scale volumes of aqueous potassium chloride solutions stored in vials at two different supersaturations (1.034 and 1.050) and laser intensities (10 and 23 MW/cm2). Higher NPLIN probabilities were observed with increasing laser-exposed volume as well as with increasing supersaturation and laser intensity. The measured NPLIN probabilities at different exposed volumes are questioned in the context of the dielectric polarization mechanism and classical nucleation theory. No significant change in the NPLIN probability was observed when samples were irradiated at the bottom, top, or middle of the vial. However, a significant increase in the nucleation probability was observed upon irradiation through the solution meniscus. We discuss these results in terms of mechanisms proposed for NPLIN.

Design and Validation of a Droplet-based Microfluidic System To Study Non-Photochemical Laser-Induced Nucleation of Potassium Chloride Solutions

Non-photochemical laser-induced nucleation (NPLIN) has emerged as a promising primary nucleation control technique offering spatiotemporal control over crystallization with potential for polymorph control. So far, NPLIN was mostly investigated in milliliter vials, through laborious manual counting of the crystallized vials by visual inspection. Microfluidics represents an alternative to acquiring automated and statistically reliable data. Thus we designed a droplet-based microfluidic platform capable of identifying the droplets with crystals emerging upon Nd:YAG laser irradiation using the deep learning method. In our experiments, we used supersaturated solutions of KCl in water, and the effect of laser intensity, wavelength (1064, 532, and 355 nm), solution supersaturation (S), solution filtration, and intentional doping with nanoparticles on the nucleation probability is quantified and compared to control cooling crystallization experiments. Ability of dielectric polarization and the nanoparticle heating mechanisms proposed for NPLIN to explain the acquired results is tested. Solutions with lower supersaturation (S = 1.05) exhibit significantly higher NPLIN probabilities than those in the control experiments for all laser wavelengths above a threshold intensity (50 MW/cm2). At higher supersaturation studied (S = 1.10), irradiation was already effective at lower laser intensities (10 MW/cm2). No significant wavelength effect was observed besides irradiation with 355 nm light at higher laser intensities (≥50 MW/cm2). Solution filtration and intentional doping experiments showed that nanoimpurities might play a significant role in explaining NPLIN phenomena.

Crystal Polymorphism Induced by Surface Tension

We use classical density functional theory (cDFT) to calculate fluid-solid surface tensions for fcc and bcc crystals formed by hard spheres and Lennard-Jones (LJ) particles. For hard spheres, our results show that the recently introduced "explicitly stable" functionals perform as well as the state of the art, and for both interaction potentials, our results compare well to simulation. We use the resulting bulk and interfacial energies for LJ to parametrize a capillary model for the free energy of small solid clusters and thereby determine the relative stability of bcc and fcc LJ clusters. We show a crossover from bcc to fcc stability as cluster size increases, thus providing insight into long-standing tension between simulation results and theoretical expectations. We also confirm that the bcc phase in contact with a vapor is unstable, thus extending earlier zero-temperature results. Our Letter demonstrates the potential of cDFT as an important tool in understanding crystallization and polymorphism.

Using classical density functional theory to determine crystal-fluid surface tensions

Classical density functional theory is used to determine the fluid-solid surface tensions for low-index faces of crystals of hard spheres and Lennard-Jones particles. The calculations make use of the recently introduced explicitly stable fundamental measure theory model for hard spheres, and we show that this gives state-of-the-art accuracy compared to simulation. For the Lennard-Jones system, results are presented for both solid-liquid and solid-vapor interfaces, and in both cases the FCC results compare favorably with existing results from the literature. We find that the BCC crystal has significantly lower solid-liquid surface tension than the FCC structure. For the solid-vapor interface, our results indicate that the BCC phase is unstable with respect to transition to the HCP structure, in agreement with various zero-temperature results in the literature.

On Capacitance and Energy Storage of Supercapacitor with Dielectric Constant Discontinuity

The classical density functional theory (CDFT) is applied to investigate influences of electrode dielectric constant on specific differential capacitance Cd and specific energy storage E of a cylindrical electrode pore electrical double layer. Throughout all calculations the electrode dielectric constant varies from 5, corresponding to a dielectric electrode, to εwr= 10⁸ corresponding to a metal electrode. Main findings are summarized as below. (i): By using a far smaller value of the solution relative dielectric constant εr=10, which matches with the reality of extremely narrow tube, one discloses that a rather high saturation voltage is needed to attain the saturation energy storage in the ultra-small pore. (ii): Use of a realistic low εr=10 value brings two obvious effects. First, influence of bulk electrolyte concentration on the Cd is rather small except when the electrode potential is around the zero charge potential; influence on the E curve is almost unobservable. Second, there remain the Cd and E enhancing effects caused by counter-ion valency rise, but strength of the effects reduces greatly with dropping of the εr value; in contrast, the Cd and E reducing effects coming from the counter-ion size enhancing remain significant enough for the low εr value. (iii) A large value of electrode relative dielectric constant εrw always reduces both the capacitance and energy storage; moreover, the effect of the εrw value gets eventually unobservable for small enough pore when the εrw value is beyond the scope corresponding to dielectric electrode. It is analyzed that the above effects take their rise in the repulsion and attraction on the counter-ions and co-ions caused by the electrode bound charges and a strengthened inter-counter-ion electrostatic repulsion originated in the low εr value.

Entropy Pair Functional Theory: Direct Entropy Evaluation Spanning Phase Transitions

We prove that, within the class of pair potential Hamiltonians, the excess entropy is a universal, temperature-independent functional of the density and pair correlation function. This result extends Henderson’s theorem, which states that the free energy is a temperature dependent functional of the density and pair correlation. The stationarity and concavity of the excess entropy functional are discussed and related to the Gibbs–Bugoliubov inequality and to the free energy. We apply the Kirkwood approximation, which is commonly used for fluids, to both fluids and solids. Approximate excess entropy functionals are developed and compared to results from thermodynamic integration. The pair functional approach gives the absolute entropy and free energy based on simulation output at a single temperature without thermodynamic integration. We argue that a functional of the type, which is strictly applicable to pair potentials, is also suitable for first principles calculation of free energies from Born–Oppenheimer molecular dynamics performed at a single temperature. This advancement has the potential to reduce the evaluation the free energy to a simple modification to any procedure that evaluates the energy and the pair correlation function.

Explicitly stable fundamental-measure-theory models for classical density functional theory

The derivation of the state of the art tensorial versions of Fundamental Measure Theory (a form of classical Density Functional Theory for hard spheres) is reexamined in the light of the recently introduced concept of global stability of the density functional based on its boundedness [Lutsko and Lam, Phys. Rev. E 98, 012604 (2018)2470-004510.1103/PhysRevE.98.012604]. It is shown that within the present paradigm, explicit stability of the functional can be achieved only at the cost of giving up accuracy at low densities. It is argued that this is an acceptable trade-off since the main value of DFT lies in the study of dense systems. Explicit calculations for a wide variety of systems show that a proposed explicitly stable functional is competitive in all ways with the popular White Bear models while sharing some of their weaknesses when applied to non-close-packed solids.

Equivalence between condensation and boiling in a Lennard-Jones fluid

Condensation and boiling are phase transitions highly relevant to industry, geology, and atmospheric science. These phase transitions are initiated by the nucleation of a drop in a supersaturated vapor and of a bubble in an overstretched liquid, respectively. The surface tension between both phases, liquid and vapor, is a key parameter in the development of such nucleation stage. Whereas the surface tension can be readily measured for a flat interface, there are technical and conceptual limitations to obtain it for the curved interface of the nucleus. On the technical side, it is quite difficult to observe a critical nucleus in experiments. From a conceptual point of view, the interfacial free energy depends on the choice of the dividing surface, being the surface of tension the one relevant for nucleation. We bypass the technical limitation by performing simulations of a Lennard-Jones fluid where we equilibrate critical nuclei (both drops and bubbles). Regarding the conceptual hurdle, we find the relevant cluster size by searching the radius that correctly predicts nucleation rates and nucleation free energy barriers when combined with Classical Nucleation Theory. With such definition of the cluster size we find the same value of the surface tension for drops and bubbles of a given radius. Thus, condensation and boiling can be viewed as two sides of the same coin. Finally, we combine the data coming from drops and bubbles to obtain, via two different routes, estimates of the Tolman length, a parameter that allows describing the curvature dependence of the surface tension in a theoretical framework.

Classical density-functional theory applied to the solid state

The standard model of classical density-functional theory (cDFT) for pair potentials consists of a hard-sphere functional plus a mean-field term accounting for long ranged attraction. However, most implementations using sophisticated fundamental measure hard-sphere functionals suffer from potential numerical instabilities either due to possible instabilities in the functionals themselves or due to implementations that mix real- and Fourier-space components inconsistently. Here we present a new implementation based on a demonstrably stable hard-sphere functional that is implemented in a completely consistent manner. The present work does not depend on approximate spherical integration schemes and so is much more robust than previous algorithms. The methods are illustrated by calculating phase diagrams for the solid state using the standard Lennard-Jones potential as well as a new class of potentials recently proposed by Wang et al. [Phys. Chem. Chem. Phys. 22, 10624 (2020)PPCPFQ1463-907610.1039/C9CP05445F]. The latter span the range from potentials for small molecules to those appropriate to colloidal systems simply by varying a parameter. We verify that cDFT is able to semiquantitatively reproduce the phase diagram in all cases. We also show that for these problems computationally cheap Gaussian approximations are nearly as good as full minimization based on finite differences.

Phase coexistence of active Brownian particles

We investigate motility-induced phase separation of active Brownian particles, which are modeled as purely repulsive spheres that move due to a constant swim force with freely diffusing orientation. We develop on the basis of power functional concepts an analytical theory for nonequilibrium phase coexistence and interfacial structure. Theoretical predictions are validated against Brownian dynamics computer simulations. We show that the internal one-body force field has four nonequilibrium contributions: (i) isotropic drag and (ii) interfacial drag forces against the forward motion, (iii) a superadiabatic spherical pressure gradient, and (iv) the quiet life gradient force. The intrinsic spherical pressure is balanced by the swim pressure, which arises from the polarization of the free interface. The quiet life force opposes the adiabatic force, which is due to the inhomogeneous density distribution. The balance of quiet life and adiabatic forces determines bulk coexistence via equality of two bulk state functions, which are independent of interfacial contributions. The internal force fields are kinematic functionals which depend on density and current but are independent of external and swim forces, consistent with power functional theory. The phase transition originates from nonequilibrium repulsion, with the agile gas being more repulsive than the quiet liquid.

How crystals form: A theory of nucleation pathways

Recent advances in classical density functional theory are combined with stochastic process theory and rare event techniques to formulate a theoretical description of nucleation, including crystallization, that can predict nonclassical nucleation pathways based on no input other than the interaction potential of the particles making up the system. The theory is formulated directly in terms of the density field, thus forgoing the need to define collective variables. It is illustrated by application to diffusion-limited nucleation of macromolecules in solution for both liquid-liquid separation and crystallization. Both involve nonclassical pathways with crystallization, in particular, proceeding by a two-step mechanism consisting of the formation of a dense-solution droplet followed by ordering originating at the core of the droplet. Furthermore, during the ordering, the free-energy surface shows shallow minima associated with the freezing of liquid into solid shells, which may shed light on the widely observed metastability of nanoscale clusters.

Solvent-mediated interactions between nanostructures: From water to Lennard-Jones liquid

Solvent-mediated interactions emerge from complex mechanisms that depend on the solute structure, its wetting properties, and the nature of the liquid. While numerous studies have focused on the first two influences, here, we compare the results from water and Lennard-Jones liquid in order to reveal to what extent solvent-mediated interactions are universal with respect to the nature of the liquid. Besides the influence of the liquid, the results were obtained with classical density functional theory and brute-force molecular dynamics simulations which allow us to contrast these two numerical techniques.