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.

Classical density functional theory in the canonical ensemble.

Classical density functional theory for finite temperatures is usually formulated in the grand-canonical ensemble where arbitrary variations of the local density are possible. However, in many cases the systems of interest are closed with respect to mass, e.g., canonical systems with fixed temperature and particle number. Although the tools of standard, grand-canonical density functional theory are often used in an ad hoc manner to study closed systems, their formulation directly in the canonical ensemble has so far not been known. In this work, the fundamental theorems underlying classical DFT are revisited and carefully compared in the two ensembles showing that there are only trivial formal differences. The practicality of DFT in the canonical ensemble is then illustrated by deriving the exact Helmholtz functional for several systems: the ideal gas, certain restricted geometries in arbitrary numbers of dimensions, and, finally, a system of two hard spheres in one dimension (hard rods) in a small cavity. Some remarkable similarities between the ensembles are apparent even for small systems with the latter showing strong echoes of the famous exact of result of Percus in the grand-canonical ensemble.

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.

Reconsidering power functional theory.

The original derivation of power functional theory [M. Schmidt and J. M. Brader, J. Chem. Phys. 138, 214101 (2013)] is reworked in some detail with a view to clarifying and simplifying the logic and making explicit the various functional dependencies. We note various issues with the original development and suggest a modification that allows us to avoid them. In the process, we also suggest an alternative interpretation of our results, which bears surprising similarities to classical density functional theory.

Long-wavelength density fluctuations as nucleation precursors.

Recent theories of nucleation that go beyond classical nucleation theory predict that diffusion-limited nucleation of both liquid droplets and of crystals from a low-density vapor (or weak solution) begins with long-wavelength density fluctuations. This means that in the early stages of nucleation, "clusters" can have low density but large spatial extent, which is at odds with the classical picture of arbitrarily small clusters of the condensed phase. We present the results of kinetic Monte Carlo simulations using forward flux sampling to show that these predictions are confirmed, namely, that on average, nucleation begins in the presence of low-amplitude, but spatially extended density fluctuations thus confirming a significant prediction of the nonclassical theory.

Antagonistic cooperativity between crystal growth modifiers.

Ubiquitous processes in nature and the industry exploit crystallization from multicomponent environments1-5; however, laboratory efforts have focused on the crystallization of pure solutes6,7 and the effects of single growth modifiers8,9. Here we examine the molecular mechanisms employed by pairs of inhibitors in blocking the crystallization of haematin, which is a model organic compound with relevance to the physiology of malaria parasites10,11. We use a combination of scanning probe microscopy and molecular modelling to demonstrate that inhibitor pairs, whose constituents adopt distinct mechanisms of haematin growth inhibition, kink blocking and step pinning12,13, exhibit both synergistic and antagonistic cooperativity depending on the inhibitor combination and applied concentrations. Synergism between two crystal growth modifiers is expected, but the antagonistic cooperativity of haematin inhibitors is not reflected in current crystal growth models. We demonstrate that kink blockers reduce the line tension of step edges, which facilitates both the nucleation of crystal layers and step propagation through the gates created by step pinners. The molecular viewpoint on cooperativity between crystallization modifiers provides guidance on the pairing of modifiers in the synthesis of crystalline materials. The proposed mechanisms indicate strategies to understand and control crystallization in both natural and engineered systems, which occurs in complex multicomponent media1-3,8,9. In a broader context, our results highlight the complexity of crystal-modifier interactions mediated by the structure and dynamics of the crystal interface.

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.

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.

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.

Classical density functional theory, unconstrained crystallization, and polymorphic behavior.

While in principle, classical density functional theory (cDFT) should be a powerful tool for the study of crystallization, in practice this has not so far been the case. Progress has been hampered by technical problems which have plagued the study of the crystalline systems using the most sophisticated fundamental measure theory models. In this paper, the reasons for the difficulties are examined and it is proposed that the tensor functionals currently favored are in fact numerically unstable. By reverting to an older, more heuristic model it is shown that all of the technical difficulties are eliminated. Application to a Lennard-Jones fluid results in a demonstration of power of cDFT to describe crystallization in a highly inhomogeneous system. First, we show that droplets attached to a slightly hydrophobic wall crystallize spontaneously upon being quenched. The resulting crystallites are clearly faceted structures and are predominantly HCP structures. In contrast, droplets in a fully periodic calculational cell remain stable to lower temperatures and eventually show the same spontaneous localization of the density into "atoms" but in an amorphous structure having many of the structural characteristics of a glass. A small change of the protocol leads, at the same temperature, to the formation of crystals, this time with the fcc structure typical of bulk Lennard-Jones solids. The fcc crystals have lower free energy than the amorphous structures which in turn are more stable than the liquid droplets. It is demonstrated that as the temperature is raised, the free energy differences between the structures decrease until the solid clusters become less stable than the liquid droplets and spontaneously melt. The presence of energy barriers separating the various structures is therefore clearly demonstrated.

Lattice induced crystallization of nanodroplets: the role of finite-size effects and substrate properties in controlling polymorphism.

Targeting specific technological applications requires the control of nanoparticle properties, especially the crystalline polymorph. Freezing a nanodroplet deposited on a solid substrate leads to the formation of crystalline structures. We study the inherent mechanisms underlying this general phenomenon by means of molecular dynamics simulations. Our work shows that different crystal structures can be selected by finely tuning the solid substrate lattice parameter. Indeed, while for our system, face-centered cubic is usually the most preponderant structure, the growth of two distinct polymorphs, hexagonal centered packing and body-centered cubic, was also observed even when the solid substrate was face-centered cubic. Finally, we also demonstrated that the growth of hexagonal centered packing is conditioned by the appearance of large enough body-centered cubic clusters thus suggesting the presence of a cross-nucleation pathway. Our results provide insights into the impact of nanoscale effects and solid substrate properties towards the growth of polymorphic nanomaterials.

Solute particle near a nanopore: influence of size and surface properties on the solvent-mediated forces.

Nanoscopic pores are used in various systems to attract nanoparticles. In general the behaviour is a result of two types of interactions: the material specific affinity and the solvent-mediated influence also called the depletion force. The latter is more universal but also much more complex to understand since it requires modeling both the nanoparticle and the solvent. Here, we employed classical density functional theory to determine the forces acting on a nanoparticle near a nanoscopic pore as a function of its hydrophobicity and its size. A simple capillary model is constructed to predict those depletion forces for various surface properties. For a nanoscopic pore, complexity arises from both the specific geometry and the fact that hydrophobic pores are not necessarily filled with liquid. Taking all of these effects into account and including electrostatic effects, we establish a phase diagram describing the entrance and the rejection of the nanoparticle from the pore.

Step Crowding Effects Dampen the Stochasticity of Crystal Growth Kinetics.

Crystals grow by laying down new layers of material which can either correspond in size to the height of one unit cell (elementary steps) or multiple unit cells (macrosteps). Surprisingly, experiments have shown that macrosteps can grow under conditions of low supersaturation and high impurity density such that elementary step growth is completely arrested. We use atomistic simulations to show that this is due to two effects: the fact that the additional layers bias fluctuations in the position of the bottom layer towards growth and by a transition, as step height increases, from a 2D to a 3D nucleation mechanism.

Mechanism for the stabilization of protein clusters above the solubility curve.

Pan, Vekilov and Lubchenko [J. Phys. Chem. B, 2010, 114, 7620] have proposed that dense stable protein clusters appearing in weak protein solutions above the solubility curve are composed of protein oligomers. The hypothesis is that a weak solution of oligomer species is unstable with respect to condensation causing the formation of dense, oligomer-rich droplets which are stabilized against growth by the monomer-oligomer reaction. Here, we show that such a combination of processes can be understood using a simple capillary model yielding analytic expressions for the cluster properties which can be used to interpret experimental data. We also construct a microscopic Dynamic Density Functional Theory model and show that it is consistent with the predictions of the capillary model. The viability of the mechanism is thus confirmed and it is shown how the radius of the stable clusters is related to physically interesting quantities such as the monomer-oligomer rate constants.

Unification of classical nucleation theories via a unified Itô-Stratonovich stochastic equation.

Classical nucleation theory (CNT) is the most widely used framework to describe the early stage of first-order phase transitions. Unfortunately, the different points of view adopted to derive it yield different kinetic equations for the probability density function, e.g., Zeldovich-Frenkel or Becker-Döring-Tunitskii equations. Starting from a phenomenological stochastic differential equation, a unified equation is obtained in this work. In other words, CNT expressions are recovered by selecting one or another stochastic calculus. Moreover, it is shown that the unified CNT thus obtained produces the same Fokker-Planck equation as that from a recent update of CNT [J. F. Lutsko and M. A. Durán-Olivencia, J. Chem. Phys. 138, 244908 (2013)10.1063/1.4811490] when mass transport is governed by diffusion. Finally, we derive a general induction-time expression along with specific approximations of it to be used under different scenarios, in particular, when the mass-transport mechanism is governed by direct impingement, volume diffusion, surface diffusion, or interface transfer.

Mesoscopic Impurities Expose a Nucleation-Limited Regime of Crystal Growth.

Nanoscale self-assembly is naturally subject to impediments at the nanoscale. The recently developed ability to follow processes at the molecular level forces us to resolve older, coarse-grained concepts in terms of their molecular mechanisms. In this Letter, we highlight one such example. We present evidence based on experimental and simulation data that one of the cornerstones of crystal growth theory, the Cabrera-Vermilyea model of step advancement in the presence of impurities, is based on incomplete physics. We demonstrate that the piercing of an impurity fence by elementary steps is not solely determined by the Gibbs-Thomson effect, as assumed by Cabrera-Vermilyea. Our data show that for conditions leading up to growth cessation, step retardation is dominated by the formation of critically sized fluctuations. The growth recovery of steps is counter to what is typically assumed, not instantaneous. Our observations on mesoscopic impurities for lysozyme expose a nucleation-dominated regime of growth that has not been hitherto considered, where the system alternates between zero and near-pure velocity. The time spent by the system in arrest is the nucleation induction time required for the step to amass a supercritical fluctuation that pierces the impurity fence.

A two-parameter extension of classical nucleation theory.

A two-variable stochastic model for diffusion-limited nucleation is developed using a formalism derived from fluctuating hydrodynamics. The model is a direct generalization of the standard classical nucleation theory (CNT). The nucleation rate and pathway are calculated in the weak-noise approximation and are shown to be in good agreement with direct numerical simulations for the weak-solution/strong-solution transition in globular proteins. We find that CNT underestimates the time needed for the formation of a critical cluster by two orders of magnitude and that this discrepancy is due to the more complex dynamics of the two variable model and not, as often is assumed, a result of errors in the estimation of the free energy barrier.

Mesoscopic nucleation theory for confined systems: a one-parameter model.

Classical nucleation theory has been recently reformulated based on fluctuating hydrodynamics [J. F. Lutsko and M. A. Durán-Olivencia, Classical nucleation theory from a dynamical approach to nucleation, J. Chem. Phys. 138, 244908 (2013). The present work extends this effort to the case of nucleation in confined systems such as small pores and vesicles. The finite available mass imposes a maximal supercritical cluster size and prohibits nucleation altogether if the system is too small. We quantity the effect of system size on the nucleation rate. We also discuss the effect of relaxing the capillary-model assumption of zero interfacial width resulting in significant changes in the nucleation barrier and nucleation rate.

Microscopic theory of anomalous diffusion based on particle interactions.

We present a master equation formulation based on a Markovian random walk model that exhibits subdiffusion, classical diffusion, and superdiffusion as a function of a single parameter. The nonclassical diffusive behavior is generated by allowing for interactions between a population of walkers. At the macroscopic level, this gives rise to a nonlinear Fokker-Planck equation. The diffusive behavior is reflected not only in the mean squared displacement [~t(γ) with 0<γ≤1.5] but also in the existence of self-similar scaling solutions of the Fokker-Planck equation. We give a physical interpretation of sub- and superdiffusion in terms of the attractive and repulsive interactions between the diffusing particles and we discuss analytically the limiting values of the exponent γ. Simulations based on the master equation are shown to be in agreement with the analytical solutions of the nonlinear Fokker-Planck equation in all three diffusion regimes.