8 Search Results

Here, we use large-scale molecular dynamics simulations coupled to free energy calculations to identify for the first time a limit of stability (spinodal) and a change in the nucleation mechanism in aqueous NaCl solutions. This is a system of considerable atmospheric, geological, and technical significance. We find that the supersaturated metastable NaCl solution reaches its limit of stability at sufficiently high salt concentrations, as indicated by the composition dependence of the salt chemical potential, indicating the transition to a phase separation by spinodal decomposition. However, the metastability limit of the NaCl solution does not correspond to spinodal decomposition with respectmore » to crystallization. We find that beyond this spinodal, a liquid/amorphous separation occurs in the aqueous solution, whereby the ions first form disordered clusters. We term these clusters as “amorphous salt.” We also identify a transition from one- to two-step crystallization mechanism driven by a spinodal. In particular, crystallization from aqueous NaCl solution beyond the spinodal is a two-step process, in which the ions first phase-separate into disordered amorphous salt clusters, followed by the crystallization of ions in the amorphous salt phase. By contrast, in the aqueous NaCl solution at concentrations lower than the spinodal, crystallization occurs via a one-step process as the ions aggregate directly into crystalline nuclei. The change of mechanism with increasing supersaturation underscores the importance of an accurate determination of the driving force for phase separation. The study has broader implications on the mechanism for nucleation of crystals from solutions at high supersaturations.« less

In this work, we use molecular dynamics simulations with a polarizable force field, namely, the modified AH/BK3 model [J. Kolafa, J. Chem. Phys. 145, 204509 (2016)], in combination with the forward flux sampling technique, to calculate the rates of homogeneous nucleation of NaCl from supersaturated aqueous solutions at 298 K and 1 bar. A non-polarizable model that reproduces the experimental equilibrium solubility {AH/TIP4P-2005 of Benavides et al. [J. Chem. Phys. 147, 104501 (2017)]} is also used for comparison. Nucleation rates calculated from the polarizable force field are found to be in good agreement with experimental measurements, while the non-polarizable modelmore » severely underestimates the nucleation rates. These results, in combination with our earlier study of a different non-polarizable force field [H. Jiang et al., J. Chem. Phys. 148, 044505 (2018)], lead to the conclusion that nucleation rates are sensitive to the details of force fields, and a good representation of nucleation rates may not be feasible using available non-polarizable force fields, even if these reproduce the equilibrium salt solubility. Finally, inclusion of polarization could be important for an accurate prediction of nucleation rates in salt solutions.« less

We used molecular dynamics simulations and the path sampling technique known as forward flux sampling to study homogeneous nucleation of NaCl crystals from supersaturated aqueous solutions at 298 K and 1 bar. Nucleation rates were obtained for a range of salt concentrations for the Joung-Cheatham NaCl force field combined with the Extended Simple Point Charge (SPC/E) water model. The calculated nucleation rates are significantly lower than the available experimental measurements. The estimates for the nucleation rates in this work do not rely on classical nucleation theory, but the pathways observed in the simulations suggest that the nucleation process is bettermore » described by classical nucleation theory than an alternative interpretation based on Ostwald’s step rule, in contrast to some prior simulations of related models. In addition to the size of NaCl nucleus, we find that the crystallinity of a nascent cluster plays an important role in the nucleation process. Nuclei with high crystallinity were found to have higher growth probability and longer lifetimes, possibly because they are less exposed to hydration water.« less

In this work, we propose a method to calculate the equilibrium contact angle of heterogeneous 3-phase solid/fluid/fluid systems using molecular dynamics simulations. The proposed method, which combines the phantom-wall method [F. Leroy and F. Müller-Plathe, J. Chem. Phys. 133, 044110 (2010)] and Bennett’s acceptance ratio approach [C. H. Bennett, J. Comput. Phys. 22, 245 (1976)], is able to calculate the solid/fluid surface tension relative to the solid surface energy. The calculated relative surface tensions can then be used in Young’s equation to estimate the equilibrium contact angle. A fluid droplet is not needed for the proposed method, in contrast tomore » the situation for direct simulations of contact angles. In addition, while prior free-energy based methods for contact angles mainly focused on the wetting of fluids in coexistence with their vapor on solid surfaces, the proposed approach was designed to study the contact angles of fluid mixtures on solid surfaces above the fluid saturation pressures. Using the proposed approach, the contact angles of binary Lennard-Jones fluid mixtures on a non-polar solid substrate were calculated at various interaction parameters and the contact angle of water in equilibrium with CO₂ on a hydrophilic polar silica surface was obtained. For both non-polar and polar systems, the calculated contact angles from the proposed method were in agreement with those obtained from the geometry of a cylindrical droplet. The computational cost of the proposed method was found to be comparable to that of simulations that use fluid droplets, but the new method provides a way to calculate the contact angle directly from Young’s equation without ambiguity.« less

Molecular simulation techniques using classical force-fields occupy the space between ab initio quantum mechanical methods and phenomenological correlations. In particular, Monte Carlo and molecular dynamics algorithms can be used to provide quantitative predictions of thermodynamic and transport properties of fluids relevant for geologic carbon sequestration at conditions for which experimental data are uncertain or not available. These methods can cover time and length scales far exceeding those of quantum chemical methods, while maintaining transferability and predictive power lacking from phenomenological correlations. The accuracy of predictions depends sensitively on the quality of the molecular models used. Many existing fixed-point-charge models formore » water and aqueous mixtures fail to represent accurately these fluid properties, especially when descriptions covering broad ranges of thermodynamic conditions are needed. Recent work on development of accurate models for water, CO₂, and dissolved salts, as well as their mixtures, is summarized in this Account. Polarizable models that can respond to the different dielectric environments in aqueous versus nonaqueous phases are necessary for predictions of properties over extended ranges of temperatures and pressures. Phase compositions and densities, activity coefficients of the dissolved salts, interfacial tensions, viscosities and diffusivities can be obtained in near-quantitative agreement to available experimental data, using relatively modest computational resources. In some cases, for example, for the composition of the CO₂-rich phase in coexistence with an aqueous phase, recent results from molecular simulations have helped discriminate among conflicting experimental data sets. The sensitivity of properties on the quality of the intermolecular interaction model varies significantly. Properties such as the phase compositions or electrolyte activity coefficients are much more sensitive than phase densities, viscosities, or component diffusivities. Strong confinement effects on physical properties in nanoscale media can also be directly obtained from molecular simulations. Future work on molecular modeling for CO₂ and aqueous brines is likely to be focused on more systematic generation of interaction models by utilizing quantum chemical as well as direct experimental measurements. New ion models need to be developed for use with the current generation of polarizable water models, including ion–ion interactions that will allow for accurate description of dense, mixed brines. Methods will need to be devised that go beyond the use of effective potentials for incorporation of quantum effects known to be important for water, and reactive force fields developed that can handle bond creation and breaking in systems with carbonate and silicate minerals. Lastly, another area of potential future work is the integration of molecular simulation methods in multiscale models for the chemical reactions leading to mineral dissolution and flow within the porous media in underground formations.« less

Not provided.

Direct coexistence molecular dynamics simulations of NaCl solutions and Lennard-Jones binary mixtures were performed to explore the origin of reported discrepancies between solubilities obtained by direct interfacial simulations and values obtained from the chemical potentials of the crystal and solution phases. We find that the key cause of these discrepancies is the use of crystal slabs of insufficient width to eliminate finite-size effects. We observe that for NaCl crystal slabs thicker than 4 nm (in the direction perpendicular to the interface), the same solubility values are obtained from the direct coexistence and chemical potential routes, namely, 3.7 ± 0.2 molalmore » at T = 298.15 K and p = 1 bar for the JC-SPC/E model. Such finite-size effects are absent in the Lennard-Jones system and are likely caused by surface dipoles present in the salt crystals. We confirmed that μs-long molecular dynamics runs are required to obtain reliable solubility values from direct coexistence calculations, provided that the initial solution conditions are near the equilibrium solubility values; even longer runs are needed for equilibration of significantly different concentrations. We do not observe any effects of the exposed crystal face on the solubility values or equilibration times. For both the NaCl and Lennard-Jones systems, the use of a spherical crystallite embedded in the solution leads to significantly higher apparent solubility values relative to the flat-interface direct coexistence calculations and the chemical potential values. Finally, our results have broad implications for the determination of solubilities of molecular models of ionic systems« less