Search

The shift in the temperature of maximum in density (TMD) at room pressure of aqueous solutions of a set of five salts containing NO 3 − and/or NH 4 + groups is studied both through experiments and through molecular dynamics simulations using the Madrid-2019 force field for ions and the TIP4P/2005 model for water. The experiments demonstrate the potential transferability and limitations of the Madrid-2019 force field for nitrate and ammonium ions recently developed by our group at different temperatures and add updated information to the reported datasets of TMDs for strong electrolytes. By using the Despretz law, individual ion contributions are extracted for predictive purposes from the experimental values of the shift in the TMD. Interesting findings for the behavior of the shift in the TMD in nitrate salts expose that this property might be particularly challenging for modelization approaches when dealing with polyatomic species. [ABSTRACT FROM AUTHOR]

The last generation of four center non-polarizable models of water can be divided into two groups: those reproducing the dielectric constant of water, as OPC, and those significantly underestimating its value, as TIP4P/2005. To evaluate the global performance of OPC and TIP4P/2005, we shall follow the test proposed by Vega and Abascal in 2011 evaluating about 40 properties to fairly address this comparison. The liquid–vapor and liquid–solid equilibria are computed, as well as the heat capacities, isothermal compressibilities, surface tensions, densities of different ice polymorphs, the density maximum, equations of state at high pressures, and transport properties. General aspects of the phase diagram are considered by comparing the ratios of different temperatures (namely, the temperature of maximum density, the melting temperature of hexagonal ice, and the critical temperature). The final scores are 7.2 for TIP4P/2005 and 6.3 for OPC. The results of this work strongly suggest that we have reached the limit of what can be achieved with non-polarizable models of water and that the attempt to reproduce the experimental dielectric constant deteriorates the global performance of the water force field. The reason is that the dielectric constant depends on two surfaces (potential energy and dipole moment surfaces), whereas in the absence of an electric field, all properties can be determined simply from just one surface (the potential energy surface). The consequences of the choice of the water model in the modeling of electrolytes in water are also discussed. [ABSTRACT FROM AUTHOR]

In this work, we present a parameterization of Sr2+ and Ba2+ cations, which expands the alkali earth set of cations of the Madrid-2019 force field. We have tested the model against the experimental densities of eight different salts, namely, SrCl2, SrBr2, SrI2, Sr(NO3)2, BaCl2, BaBr2, BaI2, and Ba(NO3)2. The force field is able to reproduce the experimental densities of all these salts up to their solubility limit. Furthermore, we have computed the viscosities for two selected salts, finding that the experimental values are overestimated, but the predictions are still reasonable. Finally, the structural properties for all the salts have been calculated with this model and align remarkably well with experimental observations. [ABSTRACT FROM AUTHOR]

In this work, we present a parameterization of Sr2+ and Ba2+ cations, which expands the alkali earth set of cations of the Madrid-2019 force field. We have tested the model against the experimental densities of eight different salts, namely, SrCl2, SrBr2, SrI2, Sr(NO3)2, BaCl2, BaBr2, BaI2, and Ba(NO3)2. The force field is able to reproduce the experimental densities of all these salts up to their solubility limit. Furthermore, we have computed the viscosities for two selected salts, finding that the experimental values are overestimated, but the predictions are still reasonable. Finally, the structural properties for all the salts have been calculated with this model and align remarkably well with experimental observations. [ABSTRACT FROM AUTHOR]

CO2 and CH4 hydrates are of great importance both from an energetic and from an environmental point of view. It is therefore highly relevant to quantify and understand the rate with which they grow. We use molecular dynamics simulations to shed light on the growth rate of these hydrates. We put the solid hydrate phase in contact with a guest aqueous solution in equilibrium with the pure guest phase and study the growth of both hydrates at 400 bars with temperature. We compare our results with previous calculations of the ice growth rate. We find a growth rate maximum as a function of the supercooling in all cases. The incorporation of guest molecules into the solid structure strongly decelerates hydrate growth. Consistently, ice grows faster than either hydrate and the CO2 hydrate grows faster than the CH4 one because of the higher solubility of CO2. We also quantify the molecular motion required to build the solids under study and find that the distance traveled by liquid molecules exceeds by orders of magnitude that advanced by any solid. Less molecular motion is needed in order for ice to grow as compared to the hydrates. Moreover, when temperature increases, more motion is needed for solid growth. Finally, we find a good agreement between our growth rate calculations and experiments of hydrate growth along the guest–solution interface. However, more work is needed to reconcile experiments of hydrate growth toward the solution among each other and with simulations. [ABSTRACT FROM AUTHOR]

In this work, we performed Molecular Dynamics simulations to investigate the solubility of $ {\rm CO}_2 $ CO 2 in both pure water and NaCl aqueous solutions, employing various force fields for water and NaCl. We first focussed on the solubility of $ {\rm CO}_2 $ CO 2 in TIP4P/2005 and TIP4P/Ice water models. We found that it was necessary to apply positive deviations to the energetic Lorentz-Berthelot combining rules for the cross interactions of water and $ {\rm CO}_2 $ CO 2 in both models to accurately replicate the experimental solubility of the gas. Then we found a negative adsorption of ions at the NaCl solution– $ {\rm CO}_2 $ CO 2 interface. We also observed an increase of the interfacial tensions of the system when adding salt which is related with the observed negative adsorption of the ions at the interface. Furthermore, we found that unit charge models tend to highly overestimate the change in the interfacial tension compared to scaled charge models. Finally, we also explored the salting out effect of $ {\rm CO}_2 $ CO 2 using different force fields for NaCl. Our findings indicate that unit charge force fields significantly overestimate the salting out effect of $ {\rm CO}_2 $ CO 2 while the scaled charge model Madrid-2019 accurately reproduces the experimental salting out effect of $ {\rm CO}_2 $ CO 2 in NaCl solutions. [ABSTRACT FROM AUTHOR]

In this work, we shall estimate via computer simulations the homogeneous nucleation rate for the methane hydrate at 400 bars for a supercooling of about 35 K. The TIP4P/ICE model and a Lennard-Jones center were used for water and methane, respectively. To estimate the nucleation rate, the seeding technique was employed. Clusters of the methane hydrate of different sizes were inserted into the aqueous phase of a two-phase gas–liquid equilibrium system at 260 K and 400 bars. Using these systems, we determined the size at which the cluster of the hydrate is critical (i.e., it has 50% probability of either growing or melting). Since nucleation rates estimated from the seeding technique are sensitive to the choice of the order parameter used to determine the size of the cluster of the solid, we considered several possibilities. We performed brute force simulations of an aqueous solution of methane in water in which the concentration of methane was several times higher than the equilibrium concentration (i.e., the solution was supersaturated). From brute force runs, we infer the value of the nucleation rate for this system rigorously. Subsequently, seeding runs were carried out for this system, and it was found that only two of the considered order parameters were able to reproduce the value of the nucleation rate obtained from brute force simulations. By using these two order parameters, we estimated the nucleation rate under experimental conditions (400 bars and 260 K) to be of the order of log10 (J/(m3 s)) = −7(5). [ABSTRACT FROM AUTHOR]

In this work, we discuss the use of scaled charges when developing force fields for NaCl in water. We shall develop force fields for Na+ and Cl− using the following values for the scaled charge (in electron units): ±0.75, ±0.80, ±0.85, and ±0.92 along with the TIP4P/2005 model of water (for which previous force fields were proposed for q = ±0.85 and q = ±1). The properties considered in this work are densities, structural properties, transport properties, surface tension, freezing point depression, and maximum in density. All the developed models were able to describe quite well the experimental values of the densities. Structural properties were well described by models with charges equal to or larger than ±0.85, surface tension by the charge ±0.92, maximum in density by the charge ±0.85, and transport properties by the charge ±0.75. The use of a scaled charge of ±0.75 is able to reproduce with high accuracy the viscosities and diffusion coefficients of NaCl solutions for the first time. We have also considered the case of KCl in water, and the results obtained were fully consistent with those of NaCl. There is no value of the scaled charge able to reproduce all the properties considered in this work. Although certainly scaled charges are not the final word in the development of force fields for electrolytes in water, its use may have some practical advantages. Certain values of the scaled charge could be the best option when the interest is to describe certain experimental properties. [ABSTRACT FROM AUTHOR]

By using the direct coexistence method, we have calculated the melting points of ice Ih at normal pressure for three recently proposed water models, namely, TIP3P-FB, TIP4P-FB, and TIP4P-D. We obtained Tm = 216 K for TIP3P-FB, Tm = 242 K for TIP4P-FB, and Tm = 247 K for TIP4P-D. We revisited the melting point of TIP4P/2005 and TIP5P obtaining Tm = 250 and 274 K, respectively. We summarize the current situation of the melting point of ice Ih for a number of water models and conclude that no model is yet able to simultaneously reproduce the melting temperature of ice Ih and the temperature of the maximum in density at room pressure. This probably points toward our both still incomplete knowledge of the potential energy surface of water and the necessity of incorporating nuclear quantum effects to describe both properties simultaneously. [ABSTRACT FROM AUTHOR]

In this work, we studied the effect of Li+, Na+, K+, Mg2+, and Ca2+ chlorides and sulfates on the temperature of maximum density (TMD) of aqueous solutions at room pressure. Experiments at 1 molal salt concentration were carried out to determine the TMD of these solutions. We also performed molecular dynamics simulations to estimate the TMD at 1 and 2 m with the Madrid-2019 force field, which uses the TIP4P/2005 water model and scaled charges for the ions, finding an excellent agreement between experiment and simulation. All the salts studied in this work shift the TMD of the solution to lower temperatures and flatten the density vs temperature curves (when compared to pure water) with increasing salt concentration. The shift in the TMD depends strongly on the nature of the electrolyte. In order to explore this dependence, we have evaluated the contribution of each ion to the shift in the TMD concluding that Na+, Ca2+, and S O 4 2 − seem to induce the largest changes among the studied ions. The volume of the system has been analyzed for salts with the same anion and different cations. These curves provide insight into the effect of different ions upon the structure of water. We claim that the TMD of electrolyte solutions entails interesting physics regarding ion–water and water–water interactions and should, therefore, be considered as a test property when developing force fields for electrolytes. This matter has been rather unnoticed for almost a century now and we believe it is time to revisit it. [ABSTRACT FROM AUTHOR]

In this work, an extension of the Madrid-2019 force field is presented. We have added the cations Rb+ and Cs+ and the anions F−, Br−, and I−. These ions were the remaining alkaline and halogen ions, not previously considered in the Madrid-2019 force field. The force field, denoted as Madrid-2019-Extended, does not include polarizability and uses the TIP4P/2005 model of water and scaled charges for the ions. A charge of ±0.85e is assigned to monovalent ions. The force field developed provides an accurate description of aqueous solution densities over a wide range of concentrations up to the solubility limit of each salt studied. Good predictions of viscosity and diffusion coefficients are obtained for concentrations below 2 m. Structural properties obtained with this force field are also in reasonable agreement with the experiment. The number of contact ion pairs has been controlled to be low so as to avoid precipitation of the system at concentrations close to the experimental solubility limit. A comprehensive comparison of the performance for aqueous solutions of alkaline halides of force fields of electrolytes using scaled and integer charges is now possible. This comparison will help in the future to learn about the benefits and limitations of the use of scaled charges to describe electrolyte solutions. [ABSTRACT FROM AUTHOR]

In this work, an extension of the Madrid-2019 force field is presented. We have added the cations Rb+ and Cs+ and the anions F−, Br−, and I−. These ions were the remaining alkaline and halogen ions, not previously considered in the Madrid-2019 force field. The force field, denoted as Madrid-2019-Extended, does not include polarizability and uses the TIP4P/2005 model of water and scaled charges for the ions. A charge of ±0.85e is assigned to monovalent ions. The force field developed provides an accurate description of aqueous solution densities over a wide range of concentrations up to the solubility limit of each salt studied. Good predictions of viscosity and diffusion coefficients are obtained for concentrations below 2 m. Structural properties obtained with this force field are also in reasonable agreement with the experiment. The number of contact ion pairs has been controlled to be low so as to avoid precipitation of the system at concentrations close to the experimental solubility limit. A comprehensive comparison of the performance for aqueous solutions of alkaline halides of force fields of electrolytes using scaled and integer charges is now possible. This comparison will help in the future to learn about the benefits and limitations of the use of scaled charges to describe electrolyte solutions. [ABSTRACT FROM AUTHOR]

The Madrid-2019 force field was recently developed to perform simulations of electrolytes in water. The model was specifically parameterized for TIP4P/2005 water and uses scaled charges for the ions. In this work, we test the compatibility of the Madrid-2019 force field with another water model: TIP4P/Ice. We shall denote this combination as Madrid-2019(TIP4P/Ice) force field. The key idea of this combination is to keep the ion-ion (Madrid-2019) and water-water (TIP4P/Ice) interactions unaltered with respect to the original models and taking the Lennard-Jones parameters for the ion-water interactions from the Madrid-2019 force field. By implementing this approach, we have maintained a reasonably good performance of the model regarding the densities and structural features of aqueous solutions, albeit yielding a moderately higher viscosity than the original model. However, the standout achievement of this new combination lies in its effective reproduction of the absolute values of the freezing temperatures of a number of ionic aqueous solutions, which could also be useful when studying hydrate formation from a two-phase system containing an aqueous solution in contact with a gas., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)