How Crystals Nucleate and Grow in Aqueous NaCl Solution

The rise and fall of adenine clusters in the gas phase: a glimpse into crystal growth and nucleation

The emergence of a crystal nucleus from disordered states is a critical and challenging aspect of the crystallization process, primarily due to the extremely short length and timescales involved. Methods such as liquid-cell or low-dose focal-series transmission electron microscopy (TEM) are often employed to probe these events. In this study, we demonstrate that ion mobility spectrometry-mass spectrometry (IMS-MS) offers a complementary and insightful perspective on the nucleation process by examining the sizes and shapes of small clusters, specifically those ranging from n = 2 to 40. Our findings reveal the significant role of sulfate ions in the growth of adeninediium sulfate clusters, which are the precursors to the formation of single crystals. Specifically, sulfate ions stabilize adenine clusters at the 1:1 ratio. In contrast, guanine sulfate forms smaller clusters with varied ratios, which become stable as they approach the 1:2 ratio. The nucleation size is predicted to be between n = 8 and 14, correlating well with the unit cell dimensions of adenine crystals. This correlation suggests that IMS-MS can identify critical nucleation sizes and provide valuable structural information consistent with established crystallographic data. We also discuss the strengths and limitations of IMS-MS in this context. IMS-MS offers rapid and robust experimental protocols, making it a valuable tool for studying the effects of various additives on the assembly of small molecules. Additionally, it aids in elucidating nucleation processes and the growth of different crystal polymorphs.

Is the Local Ion Density Sufficient to Drive NaCl Nucleation from the Melt and Aqueous Solution?

Even though nucleation is ubiquitous in different science and engineering problems, investigating nucleation is extremely difficult due to the complicated ranges of time and length scales involved. In this work, we simulate NaCl nucleation in both molten and aqueous environments using enhanced sampling of all-atom molecular dynamics with deep-learning-based estimation of reaction coordinates. By incorporating various structural order parameters and learning the reaction coordinate as a function thereof, we achieve significantly improved sampling relative to traditional ad hoc descriptions of what drives nucleation, particularly in an aqueous medium. Our results reveal a one-step nucleation mechanism in both environments, with reaction coordinate analysis highlighting the importance of local ion density in distinguishing solid and liquid states. However, although fluctuations in the local ion density are necessary to drive nucleation, they are not sufficient. Our analysis shows that near the transition states, descriptors such as enthalpy and local structure become crucial. Our protocol proposed here enables robust nucleation analysis and phase sampling and could offer insights into nucleation mechanisms for generic small molecules in different environments.

Variation of Surface Propensity of Halides with Droplet Size and Temperature: The Planar Interface Limit

The radial number density profiles of halide and alkali ions in aqueous clusters with equimolar radius ≲1.4 nm, which correspond to ≲255 H2O molecules, have been extensively studied by computations. However, the surface abundance of Cl-, Br-, and I- relative to the bulk interior in these smaller clusters may not be representative of the larger systems. Indeed, here we show that the larger the cluster is, the lower the relative surface abundance of chaotropic halides is. In droplets with an equimolar radius of ≈2.45 nm, which corresponds to ≈2000 H2O molecules, the polarizable halides show a clear number density maximum in the droplet's bulk-like interior. A similar pattern is observed in simulations of the aqueous planar interface with halide salts at room temperature. At elevated temperature the surface propensity of Cl- decreases gradually, while that of I- is partially preserved. The change in the chaotropic halide location at higher temperatures than the room temperature may considerably affect photochemical reactivity in atmospheric aerosols, vapor-liquid nucleation and growth mechanisms, and salt crystallization via solvent evaporation. We argue that the commonly used approach of nullifying parameters in a force field in order to find the factors that determine the ion location does not provide transferable insight into other force fields.

Resolving Temperature-Dependent Hydrate Nucleation Pathway: The Role of "Transition Layer"

Understanding the nucleation of natural gas hydrate (NGH) at different conditions has important implications to NGH recovery and other industrial applications, such as gas storage and separation. Herein, vast numbers of hydrate nucleation events are traced via molecular dynamics (MD) simulations at different degrees of supercooling (or driving forces). Specifically, to precisely characterize a hydrate nucleus from an aqueous system during the MD simulation, we develop an evolutionary order parameter (OP) to recognize the nucleus size and shape. Subsequently, the free energy landscapes of hydrate during nucleation are explored by using the newly developed OP. The results suggest that at 270 K (or 0.92 Tm supercooling, where Tm is the melting point), the near-rounded nucleus prevails during the nucleation, as described from the classical nucleation theory. In contrast, at relatively strong driving forces of 0.85 and 0.88 Tm, nonclassical nucleation events arise. Specifically, the pathway toward an elongated nucleus becomes as important as the pathway toward a near-rounded nucleus. To explain the distinct nucleation phenomena at different supercoolings, a notion of a "transition layer" (or liquid-blob-like layer) is proposed. Here, the transition layer is to describe the interfacial region between the nucleus and aqueous solution, and this layer entails two functionalities: (1) it tends to retain CH4 depending on the degrees of supercooling and (2) it facilitates collision among CH4, which thus promote the incorporation of CH4 into nucleus. Our simulation indicates that compared to the near-rounded nucleus, the transition layer surrounding the elongated nucleus is more evident with the higher collision rate among CH4 molecules. As such, the transition layer tends to promote the elongated nucleus pathway, while offsetting the cost of larger surface free energy associated with the elongated nucleus. At 0.92 Tm, however, the transition layer gradually disappears, and classical nucleation events dominate. Overall, the notion of "transition layer" offers deeper insight into the NGH nucleation at different degrees of supercooling and could be extended to describe other types of hydrate nucleation.

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.

Crystallization of FAPbI3: Polytypes and stacking faults

Molecular dynamics simulations are performed to study the crystallization of formamidinium lead iodide. From all-atom simulations of the crystal growth process and the δ-α-phase transitions, we try to reveal the formation of various stack-faulted intermediate defected structures and report various polytypes of formamidinium lead iodide that are observed from simulations.

Exploring the Barriers in the Aggregation of a Hexadecameric Human Prion Peptide through the Markov State Model

The prefibrillar aggregation kinetics of prion peptides are still an enigma. In this perspective, we employ atomistic molecular dynamics (MD) simulations of the shortest human prion peptide (HPP) (127GYMLGS132) at various temperatures and peptide concentrations and apply the Markov state model to determine the various intermediates and lag phases. Our results reveal that the natural mechanism of prion peptide self-assembly in the aqueous phase is impeded by two significant kinetic barriers with oligomer sizes of 6-9 and 12-13 peptides, respectively. The first one is the aggregation of unstructured lower-order oligomers, and the second is fibril nucleation, which impedes the further growth of prion aggregates. Among these two activation barriers, the second one is found to be dominant irrespective of the increase in temperature and peptide concentration. These lag phases are captured in all three different force-field parameters, namely, GROMOS-54a7, AMBER-99SB-ILDN, and CHARMMS 36m, at different concentrations. The GROMOS-54a7 and AMBER-99SB-ILDN force fields showed a comparatively higher percentage of β-sheet formation in the metastable aggregate that evolved during the aggregation process. In contrast, the CHARMM-36m force field showed mostly coil or turn conformations. The addition of a novel catecholamine derivative (naphthoquinone dopamine (NQDA)) arrests the aggregation process between the lag phases by increasing the activation barrier for the Lag1 and Lag2 phases in all of the force fields, which further validates the existence of these lag phases. The preferential binding of NQDA with the peptides increases the hydration of peptides and eventually disrupts the organized morphology of prefibrillar aggregates. It reduces the dimer dissociation energy by -24.34 kJ/mol.

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.

A Tale of Two Foulants: The Coupling of Organic Fouling and Mineral Scaling in Membrane Desalination

Membrane desalination that enables the harvesting of purified water from unconventional sources such as seawater, brackish groundwater, and wastewater has become indispensable to ensure sustainable freshwater supply in the context of a changing climate. However, the efficiency of membrane desalination is greatly constrained by organic fouling and mineral scaling. Although extensive studies have focused on understanding membrane fouling or scaling separately, organic foulants commonly coexist with inorganic scalants in the feedwaters of membrane desalination. Compared to individual fouling or scaling, combined fouling and scaling often exhibits different behaviors and is governed by foulant-scalant interactions, resembling more complex but practical scenarios than using feedwaters containing only organic foulants or inorganic scalants. In this critical review, we first summarize the performance of membrane desalination under combined fouling and scaling, involving mineral scales formed via both crystallization and polymerization. We then provide the state-of-the-art knowledge and characterization techniques pertaining to the molecular interactions between organic foulants and inorganic scalants, which alter the kinetics and thermodynamics of mineral nucleation as well as the deposition of mineral scales onto membrane surfaces. We further review the current efforts of mitigating combined fouling and scaling via membrane materials development and pretreatment. Finally, we provide prospects for future research needs that guide the design of more effective control strategies for combined fouling and scaling to improve the efficiency and resilience of membrane desalination for the treatment of feedwaters with complex compositions.

Species-selective nanoreactor molecular dynamics simulations based on linear-scaling tight-binding quantum chemical calculations

Here, extensions to quantum chemical nanoreactor molecular dynamics simulations for discovering complex reactive events are presented. The species-selective algorithm, where the nanoreactor effectively works for the selected desired reactants, was introduced to the original scheme. Moreover, for efficient simulations of large model systems with the modified approach, the divide-and-conquer linear-scaling density functional tight-binding method was exploited. Two illustrative applications of the polymerization of propylene and cyclopropane mixtures and the aggregation of sodium chloride from aqueous solutions indicate that species-selective quantum chemical nanoreactor molecular dynamics is a promising method to accelerate the sampling of multicomponent chemical processes proceeding under relatively mild conditions.

Nucleation Rate Theory for Coordination Number: Elucidating Water-Mediated Formation of a Zigzag Na2SO4 Morphology

Predicting and controlling nanostructure formation during nucleation can pave the way to synthesizing novel energy materials via crystallization. However, such control over nucleation and crystallization remains challenging due to an inadequate understanding of critical factors that govern evolving atomistic structures and dynamics. Herein, we utilize coordination number as a reaction coordinate and rate theory to investigate how sodium sulfate, commonly known as a phase-change energy material, nucleates in a supersaturated aqueous solution. In conjunction with ab initio and force field-based molecular dynamics simulation, the rate theoretical analysis reveals that sodium sulfate from an initially dissolved metastable state transits to a heterogeneous mixture of prenucleated clusters and finally to a large cylindrical zigzag morphology. Measurements of Raman spectra and their ab initio modeling confirm that this nucleated morphology contains a few waters for every sulfate. Rate processes such as solvent exchange and desolvation exhibit high sensitivity to the evolving prenucleation/nucleation structures, providing a means to distinguish between critical nucleation precursors. Desolvation and forming the first-shell interionic coordination structure via monomer-by-monomer addition around sulfates are found to explain the formation of large nuclei. Thus, a detailed understanding of the step-by-step structure formation across scales has been achieved. This can be leveraged to predict nucleation-related structures and dynamics and potentially control the synthesis of novel phase-change materials for energy applications.

Study of brine-halite phase separation through optical constringence and molecular dynamics

This work presents a study of the reciprocal dispersive power, also known as constringence or Abbe number of an aqueous solution of NaCl in a wide range of concentrations. The constringence exhibited a distinct behavior in the region close to the phase transition between a phase containing exclusively brine and a phase containing brine+halite. Molecular dynamics simulations of this system indicated the existence of halite formation below the known saturation curve, which agreed with the experimental measurements, indicating a crystal growth in the unsaturated region.

Salt Enrichment and Dynamics in the Interface of Supercooled Aqueous Droplets

The interconversion reaction of NaCl between the contact-ion pair (CIP) and the solvent-separated ion pair (SSIP) as well as the free-ion state in cold droplets has not yet been investigated. We report direct computational evidence that the lower is the temperature, the closer to the surface the ion interconversion reaction takes place. In supercooled droplets the enrichment of the subsurface in salt becomes more evident. The stability of the SSIP relative to the CIP increases as the ion-pairing is transferred toward the droplet's outer layers. In the free-ion state, where the ions diffuse independently in the solution, the number density of Cl^- shows a broad maximum in the interior in addition to the well-known maximum in the surface. In the study of the reaction dynamics, we find a weak coupling between the interionic NaCl distance reaction coordinate and the solvent degrees of freedom, which contrasts with the diffusive crossing of the free energy barrier found in bulk solution modeling. The H₂O self-diffusion coefficient is found to be at least an order of magnitude larger than that in the bulk solution. We propose to exploit the enhanced surface ion concentration at low temperature to eliminate salts from droplets in native mass spectrometry ionization methods.

An investigation of the kinetics and thermodynamics of NaCl nucleation through composite clusters

Having a good understanding of nucleation is critical for the control of many important processes, such as polymorph selection during crystallization. However, a complete picture of the molecular-level mechanisms of nucleation remains elusive. In this work, we take an in-depth look at the NaCl homogeneous nucleation mechanism through thermodynamics. Distinguished from the classical nucleation theory, we calculate the free energy of nucleation as a function of two nucleus size coordinates: crystalline and amorphous cluster sizes. The free energy surface reveals a thermodynamic preference for a nonclassical mechanism of nucleation through a composite cluster, where the crystalline nucleus is surrounded by an amorphous layer. The thickness of the amorphous layer increases with an increase in supersaturation. The computed free energy landscape agrees well with the composite cluster-free energy model, through which phase specific thermodynamic properties are evaluated. As the supersaturation increases, there is a change in stability of the amorphous phase relative to the solution phase, resulting in a change from one-step to two-step mechanism, seen clearly from the free energy profile along the minimum free energy path crossing the transition curve. By obtaining phase-specific diffusion coefficients, we construct the full mesoscopic model and present a clear roadmap for NaCl nucleation.

Fluorescence of KCl Aqueous Solution: A Possible Spectroscopic Signature of Nucleation

Ion pairing in water solutions alters both the water hydrogen-bond network and ion solvation, modifying the dynamics and properties of electrolyte water solutions. Here, we report an anomalous intrinsic fluorescence of KCl aqueous solution at room temperature and show that its intensity increases with the salt concentration. From the ab initio density functional theory (DFT) and time-dependent DFT modeling, we propose that the fluorescence emission could originate from the stiffening of the hydrogen bond network in the hydration shell of solvated ion-pairs that suppresses the fast nonradiative decay and allows the slower radiative channel to become a possible decay pathway. Because computations suggest that the fluorophores are the local ion-water structures present in the prenucleation phase, this band could be the signature of the incoming salt precipitation.

Effects of salinity and shear stress on clay deformation: A molecular dynamics study

The deformation of clay minerals is an important phenomenon that is relevant to many problems, particularly those that occur in subsurface geological formations. The salinity of the formations and external shear stress applied to them are two important factors that contribute to the deformation of such porous media. To gain a deeper understanding of such phenomena, we have carried out extensive molecular dynamics simulations using the Na-montmorillonite (Na-MMT) structure as the model of clay minerals and have studied the effect of salt concentration on its swelling. As the NaCl concentration increases, so also does the basal spacing. We demonstrate the effect of the coupling between the applied shear stress and NaCl salinity on the swelling behavior of Na-MMT, namely, deformation of the interlayer space that results in swelling. According to the results, the extent of Na-MMT deformation depends on both the brine salinity and the shear rate.

NaCl aggregation in water at elevated temperatures and pressures: Comparison of classical force fields

The properties of water vary dramatically with temperature and density. This can be exploited to control its effectiveness as a solvent. Thus, supercritical water is of keen interest as solvent in many extraction processes. The low solubility of salts in lower density supercritical water has even been suggested as a means of desalination. The high temperatures and pressures required to reach supercritical conditions can present experimental challenges during collection of required physical property and phase equilibria data, especially in salt-containing systems. Molecular simulations have the potential to be a valuable tool for examining the behavior of solvated ions at these high temperatures and pressures. However, the accuracy of classical force fields under these conditions is unclear. We have, therefore, undertaken a parametric study of NaCl in water, comparing several salt and water models at 200 bar-600 bar and 450 K-750 K for a range of salt concentrations. We report a comparison of structural properties including ion aggregation, hydrogen bonding, density, and static dielectric constants. All of the force fields qualitatively reproduce the trends in the liquid phase density. An increase in ion aggregation with decreasing density holds true for all of the force fields. The propensity to aggregate is primarily determined by the salt force field rather than the water force field. This coincides with a decrease in the water static dielectric constant and reduced charge screening. While a decrease in the static dielectric constant with increasing NaCl concentration is consistent across all model combinations, the salt force fields that exhibit more ionic aggregation yield a slightly smaller dielectric decrement.

Multiple Pathways in NaCl Homogeneous Crystal Nucleation

NaCl crystal nucleation from metastable solutions has long been considered to occur according to a single-step mechanism where the growth in the size and crystalline order of the emerging nuclei...

Hydration breaking and chemical ordering in a levitated NaCl solution droplet beyond the metastable zone width limit: evidence for the early stage of two-step nucleation

For over two decades, NaCl nucleation from a supersaturated aqueous solution has been predicted to occur via a two-step nucleation (TSN) mechanism, i.e., two sequential events, the formation of locally dense liquid regions followed by structural ordering. However, the formation of dense liquid regions in the very early stage of TSN has never been experimentally observed. By using a state-of-the-art technique, a combination of electrostatic levitation (ESL) and in situ synchrotron X-ray and Raman scatterings, we find experimental evidence that indicates the formation of dense liquid regions in NaCl bulk solution at an unprecedentedly high level of supersaturation (S = 2.31). As supersaturation increases, evolution of ion clusters leads to chemical ordering, but no topological ordering, which is a precursor for forming the dense disordered regions of ion clusters at the early stage of TSN. Moreover, as the ion clusters proceed to evolve under highly supersaturated conditions, we observe the breakage of the water hydration structure indicating the stability limit of the dense liquid regions, and thus leading to nucleation. The evolution of solute clusters and breakage of hydration in highly supersaturated NaCl bulk solution will provide new insights into the detailed mechanism of TSN for many other aqueous solutions.

This work provides evidence for two-step nucleation in highly supersaturated bulk NaCl solution, using electrostatic levitation combined with Raman/X-ray scatterings.

Effect of Choice of Solvent on Crystallization Pathway of Paracetamol: An Experimental and Theoretical Case Study

The choice of solvents influences crystalline solid formed during the crystallization of active pharmaceutical ingredients (API). The underlying effects are not always well understood because of the complexity of the systems. Theoretical models are often insufficient to describe this phenomenon. In this study, the crystallization behavior of the model drug paracetamol in different solvents was studied based on experimental and molecular dynamics data. The crystallization process was followed in situ using time-resolved Raman spectroscopy. Molecular dynamics with simulated annealing algorithm was used for an atomistic understanding of the underlying processes. The experimental and theoretical data indicate that paracetamol molecules adopt a particular geometry in a given solvent predefining the crystallization of certain polymorphs.

Lennard-Jones Parameters Determined to Reproduce the Solubility of NaCl and KCl in SPC/E, TIP3P, and TIP4P/2005 Water

Most classical nonpolarizable ion potential models underestimate the solubility values of NaCl and KCl in water significantly. We determine Lennard-Jones parameters of Na⁺, K⁺, and Cl^- that reproduce the solubility as well as the hydration free energy in dilute aqueous solutions for three water potential models, SPC/E, TIP3P, and TIP4P/2005. The ion-oxygen distance in the solution and the cation-anion distance in salt are also considered in the parametrization. In addition to the target properties, the hydration enthalpy, hydration entropy, self-diffusion coefficient, coordination number, lattice energy, enthalpy of solution, density, viscosity, and number of contact ion pairs are calculated for comparison with 17 frequently used or recently developed ion potential models. The overall performance of each ion model is represented by a global score using a scheme that was originally developed for comparison of water potential models. The global score is better for our models than for the other 17 models not only because of the quite good prediction for the solubility but also because of the relatively small deviation from the experimental value for many of the other properties.

Homogeneous Nucleation Mechanism of NaCl in Aqueous Solutions

In this study, molecular dynamic simulations are employed to investigate the homogeneous nucleation mechanism of NaCl crystal in solutions. According to the simulations, the dissolved behaviors of NaCl in water are dependent on ion concentrations. With increasing NaCl concentrations, the dissolved Na+ and Cl- ions tend to be aggregated in solutions. In combination with our recent studies, the aggregate of dissolved solutes is mainly ascribed to the hydrophobic interactions. Different from the two-step mechanism, no barrier is needed to overcome the formation of the aggregate. In comparison with the classical nucleation theory (CNT), because of the formation of solute aggregate, this lowers the barrier height of nucleation and affects the nucleation mechanism of NaCl crystal in water.

Molecular insights on NaCl crystal formation approaching PVDF membranes functionalized with graphene

Atomistic simulations of graphene–PVDF membranes speeding up NaCl crystal nucleation and growth in comparison to the pristine PVDF membranes.

Molecular dynamics simulation of homogeneous nucleation of supersaturated potassium chloride (KCl) in aqueous solutions

Molecular dynamics (MD) simulation is used to investigate the mechanism of crystal nucleation of potassium chloride (KCl) in a supersaturated aqueous solution at 293 K and 1 atm.

Hopper Growth of Salt Crystals

The growth of hopper crystals is observed for many substances, but the mechanism of their formation remains ill understood. Here we investigate their growth by performing evaporation experiments on small volumes of salt solutions. We show that sodium chloride crystals that grow very fast from a highly supersaturated solution form a peculiar form of hopper crystal consisting of a series of connected miniature versions of the original cubic crystal. The transition between cubic and such hopper growth happens at a well-defined supersaturation where the growth rate of the cubic crystal reaches a maximum (∼6.5 ± 1.8 μm/s). Above this threshold, the growth rate varies as the third power of supersaturation, showing that a new mechanism, controlled by the maximum speed of surface integration of new molecules, induces the hopper growth of cubic crystals in cascade.

Effects of nanoparticle heating on the structure of a concentrated aqueous salt solution

The effects of a rapidly heated nanoparticle on the structure of a concentrated aqueous salt solution are studied using molecular dynamics simulations. A diamond-like nanoparticle of radius 20 Å is immersed in a sodium-chloride solution at 20% above the experimental saturation concentration and equilibrated at T = 293 K and P = 1 atm. The nanoparticle is then rapidly heated to several thousand degrees Kelvin, and the system is held under isobaric-isoenthalpic conditions. It is observed that after 2-3 ns, the salt ions are depleted far more than water molecules from a proximal zone 15-25 Å from the nanoparticle surface. This leads to a transient reduction in molality in the proximal zone and an increase in ion clustering in the distal zone. At longer times, ions begin to diffuse back into the proximal zone. It is speculated that the formation of proximal and distal zones, and the increase in ion clustering, plays a role in the mechanism of nonphotochemical laser-induced nucleation.

Self-assembly of MoO3 needles in gas current for cubic formation pathway

Nucleation and subsequent rapid growth are enigmatic due to the unrevealed pathways. Despite the relatively simpler mechanism compared to nucleation and growth in solution, that in vapor has received little attention. The largest hindrance to unveiling this process may be observing the rapid and mesoscopic-scale phenomena. To overcome this hindrance, we combine an experimental approach with in situ spatial scanning Fourier-transform infrared spectroscopy, which reveals the nucleating and growing nanoparticles in vapor. The nanoparticles are then collected at different evolutionary stages and analyzed by ex situ transmission electron microscopy (TEM). Needle-shaped molybdenum oxide (MoO₃) nanoparticles were formed within ∼0.1 s after homogeneous nucleation from a highly supersaturated vapor. Over one second, the needle particles gradually evolved into a cubic shape by fusion in a crystallographically favored orientation in a free-flying state in vapor. The similar sizes of the elongated axes of the needle and cubic structures suggest an additional growth stage, in which the needle particles become the growth units of the cubic particles. The morphology of a final crystal should reflect the formation environment of the particle because growing crystals are sensitive to the formation conditions such as temperature, concentration, and impurities. Although nucleation under very high supersaturation induces the anisotropic growth of the needle particles, this information of the initial nucleation environment is lost in the final cubic crystal. These findings enrich our understanding of pathways in the nucleation and growth of nanoparticles and provide new insights into the growth stages driven by oriented attachment.

Consensus on the solubility of NaCl in water from computer simulations using the chemical potential route

The solubility of NaCl in water is evaluated by using three force field models: Joung-Cheatham for NaCl dissolved in two different water models (SPC/E and TIP4P/2005) and Smith Dang NaCl model in SPC/E water. The methodology based on free-energy calculations [E. Sanz and C. Vega, J. Chem. Phys. 126, 014507 (2007)] and [J. L. Aragones et al., J. Chem. Phys. 136, 244508 (2012)] has been used, except, that all calculations for the NaCl in solution were obtained by using molecular dynamics simulations with the GROMACS package instead of homemade MC programs. We have explored new lower molalities and made longer runs to improve the accuracy of the calculations. Exploring the low molality region allowed us to obtain an analytical expression for the chemical potential of the ions in solution as a function of molality valid for a wider range of molalities, including the infinite dilute case. These new results are in better agreement with recent estimations of the solubility obtained with other methodologies. Besides, two empirical simple rules have been obtained to have a rough estimate of the solubility of a certain model, by analyzing the ionic pairs formation as a function of molality and/or by calculating the difference between the NaCl solid chemical potential and the standard chemical potential of the salt in solution.

Birth of NaCl Crystals: Insights from Molecular Simulations

Molecular dynamics simulations are used to investigate the factors that influence the nucleation of NaCl crystals in a supersaturated aqueous solution. We describe a methodology for detecting solidlike NaCl clusters (potential nuclei) and following their evolution in time until they achieve nucleation (which is very rare) or dissolve back into solution. Through an analysis of cluster lifetimes and multiple nucleation events, we demonstrate that cluster size is not the only property that influences cluster stability and the probability of achieving nucleation. We introduce a parameter called cluster crystallinity, which is a measure of the solidlike order in a particular cluster. We show that cluster order (as measured by this parameter) has a strong influence on the lifetime and nucleation probability of clusters of equal sizes, with the lifetime and probability of nucleation increasing with increasing crystallinity. These observations remain true for clusters as small as six ions, showing that the structural factors are important even at the earliest stages of crystal birth.

Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations

The nucleation of crystals in liquids is one of nature's most ubiquitous phenomena, playing an important role in areas such as climate change and the production of drugs. As the early stages of nucleation involve exceedingly small time and length scales, atomistic computer simulations can provide unique insights into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that, in the past few decades, have unraveled crucial aspects of crystal nucleation in liquids. We put into context the theoretical framework of classical nucleation theory and the state-of-the-art computational methods by reviewing simulations of such processes as ice nucleation and the crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insights into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that, as a result, the general applicability of classical nucleation theory has been repeatedly called into question. We have attempted to identify the most pressing open questions in the field. We believe that, by improving (i) existing interatomic potentials and (ii) currently available enhanced sampling methods, the community can move toward accurate investigations of realistic systems of practical interest, thus bringing simulations a step closer to experiments.

Carbon-Encapsulated Co3O4 Nanoparticles as Anode Materials with Super Lithium Storage Performance

A high-performance anode material for lithium storage was successfully synthesized by glucose as carbon source and cobalt nitrate as Co3O4 precursor with the assistance of sodium chloride surface as a template to reduce the carbon sheet thickness. Ultrafine Co3O4 nanoparticles were homogeneously embedded in ultrathin porous graphitic carbon in this material. The carbon sheets, which have large specific surface area, high electronic conductivity, and outstanding mechanical flexibility, are very effective to keep the stability of Co3O4 nanoparticles which has a large capacity. As a consequence, a very high reversible capacity of up to 1413 mA h g(-1) at a current density of 0.1 A g(-1) after 100 cycles, a high rate capability (845, 560, 461 and 345 mA h g(-1) at 5, 10, 15 and 20 C, respectively, 1 C = 1 A g(-1)), and a superior cycling performance at an ultrahigh rate (760 mA h g(-1) at 5 C after 1000 cycles) are achieved by this lithium-ion-battery anode material.

Carbon-Encapsulated Co3O4 Nanoparticles as Anode Materials with Super Lithium Storage Performance

A high-performance anode material for lithium storage was successfully synthesized by glucose as carbon source and cobalt nitrate as Co₃O₄ precursor with the assistance of sodium chloride surface as a template to reduce the carbon sheet thickness. Ultrafine Co₃O₄ nanoparticles were homogeneously embedded in ultrathin porous graphitic carbon in this material. The carbon sheets, which have large specific surface area, high electronic conductivity, and outstanding mechanical flexibility, are very effective to keep the stability of Co₃O₄ nanoparticales which has a large capacity. As a consequence, a very high reversible capacity of up to 1413 mA h g⁻¹ at a current density of 0.1 A g⁻¹ after 100 cycles, a high rate capability (845, 560, 461 and 345 mA h g⁻¹ at 5, 10, 15 and 20 C, respectively, 1 C = 1 A g⁻¹), and a superior cycling performance at an ultrahigh rate (760 mA h g⁻¹ at 5 C after 1000 cycles) are achieved by this lithium-ion-battery anode material.