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]

CNDOL is an a priori, approximate Fockian for molecular wave functions. In this study, we employ several modes of singly excited configuration interaction (CIS) to model molecular excitation properties by using four combinations of the one electron operator terms. Those options are compared to the experimental and theoretical data for a carefully selected set of molecules. The resulting excitons are represented by CIS wave functions that encompass all valence electrons in the system for each excited state energy. The Coulomb–exchange term associated to the calculated excitation energies is rationalized to evaluate theoretical exciton binding energies. This property is shown to be useful for discriminating the charge donation ability of molecular and supermolecular systems. Multielectronic 3D maps of exciton formal charges are showcased, demonstrating the applicability of these approximate wave functions for modeling properties of large molecules and clusters at nanoscales. This modeling proves useful in designing molecular photovoltaic devices. Our methodology holds potential applications in systematic evaluations of such systems and the development of fundamental artificial intelligence databases for predicting related properties. [ABSTRACT FROM AUTHOR]

In this work, we introduce variational umbrella seeding, a novel technique for computing nucleation barriers. This new method, a refinement of the original seeding approach, is far less sensitive to the choice of order parameter for measuring the size of a nucleus. Consequently, it surpasses seeding in accuracy and umbrella sampling in computational speed. We test the method extensively and demonstrate excellent accuracy for crystal nucleation of nearly hard spheres and two distinct models of water: mW and TIP4P/ICE. This method can easily be extended to calculate nucleation barriers for homogeneous melting, condensation, and cavitation. [ABSTRACT FROM AUTHOR]

Clathrate hydrates are vital in energy research and environmental applications. Understanding their stability is crucial for harnessing their potential. In this work, we employ direct coexistence simulations to study finite-size effects in the determination of the three-phase equilibrium temperature (T3) for methane hydrates. Two popular water models, TIP4P/Ice and TIP4P/2005, are employed, exploring various system sizes by varying the number of molecules in the hydrate, liquid, and gas phases. The results reveal that finite-size effects play a crucial role in determining T3. The study includes nine configurations with varying system sizes, demonstrating that smaller systems, particularly those leading to stoichiometric conditions and bubble formation, may yield inaccurate T3 values. The emergence of methane bubbles within the liquid phase, observed in smaller configurations, significantly influences the behavior of the system and can lead to erroneous temperature estimations. Our findings reveal finite-size effects on the calculation of T3 by direct coexistence simulations and clarify the system size convergence for both models, shedding light on discrepancies found in the literature. The results contribute to a deeper understanding of the phase equilibrium of gas hydrates and offer valuable information for future research in this field. [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]

The importance of nitrate and ammonium salts both in the environment and in biological processes cannot be questioned. In this work, using the TIP4P/2005 water model, aqueous solutions of nitrate and ammonium electrolytes are parametrized using scaled charges while keeping a rigid structure and nonpolarizable charge distributions. The models are optimized by systematically testing a set of properties for twelve electrolytes—eight nitrate and four ammonium salts—thus, enlarging the number of potential chemical species encompassed within the Madrid-2019 force field for ions. The capacity of the force field for predicting densities, ion–ion and ion–solvent structures, and transport properties of the solutions comprised by the trial batch of salts was tested and discussed. Both the dependence of the densities with the salt concentration and the solution structure were nicely reproduced by the models in the whole concentration range without any trace of precipitating events and with improved accuracy in comparison with recently reported models, while the agreement of the simulated transport properties with experimental data ranges from good to reasonable, depending on the ion/counterion pair. These scaled charge models might be considered as force fields embodying a reasonable compromise between exactness and general applicability and also as an important step in the development of accurate models for polyatomic ions. [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]

A new approach for generating Gaussian basis sets is reported and tested for atoms from H to Ne. The basis sets thus calculated, named SIGMA basis sets, range from DZ to QZ sizes and have the same composition per shell as Dunning basis sets but with different treatment of the contractions. The standard SIGMA basis sets and their augmented versions have proven to be very suitable for providing good results in atomic and molecular calculations. The performance of the new basis sets is analyzed in terms of total, correlation, and atomization energies, equilibrium distances, and vibrational frequencies in several molecules, and the results are compared at several computational levels with those obtained with the corresponding Dunning and other basis sets. [ABSTRACT FROM AUTHOR]

One of the most accepted hypothesis to explain the anomalous behavior of water is the presence of a critical point between two liquids, the liquid–liquid critical point (LLCP), buried within the deep supercooled regime. Unfortunately, such hypothesis is hard to be experimentally confirmed due to fast freezing. Here, we show that the TIP4P/Ice water potential shifted by 400 bar can reproduce with unprecedented accuracy the experimental isothermal compressibility of water and its liquid equation of state for a wide pressure and temperature range. We find, both by extrapolation of response function maxima and by a Maxwell construction, that the location of the model LLCP is consistent with previous calculations. According to the pressure shift needed to recover the experimental behavior of supercooled water, we estimate the experimental LLCP to be located around 1250 bar and 195 K. We use the model to estimate the ice nucleation rate (J) in the vicinity of the hypothesized LLCP experimental location and obtain J = 1024 m−3 s−1. Thereby, experiments where the ratio between the cooling rate and the sample volume is equal or larger than the estimated nucleation rate could probe liquid–liquid equilibrium before freezing. Such conditions are not accessible in common experiments with microdroplets cooled at a few kelvin per second, but they could be, for instance, using nanodroplets of around 50 nm radius observed in a millisecond timescale. [ABSTRACT FROM AUTHOR]

In this paper, the solubility of carbon dioxide (CO2) in water along the isobar of 400 bar is determined by computer simulations using the well-known TIP4P/Ice force field for water and the TraPPE model for CO2. In particular, the solubility of CO2 in water when in contact with the CO2 liquid phase and the solubility of CO2 in water when in contact with the hydrate have been determined. The solubility of CO2 in a liquid–liquid system decreases as the temperature increases. The solubility of CO2 in a hydrate–liquid system increases with temperature. The two curves intersect at a certain temperature that determines the dissociation temperature of the hydrate at 400 bar (T3). We compare the predictions with T3 obtained using the direct coexistence technique in a previous work. The results of both methods agree, and we suggest 290(2) K as the value of T3 for this system using the same cutoff distance for dispersive interactions. We also propose a novel and alternative route to evaluate the change in chemical potential for the formation of hydrates along the isobar. The new approach is based on the use of the solubility curve of CO2 when the aqueous solution is in contact with the hydrate phase. It considers rigorously the non-ideality of the aqueous solution of CO2, providing reliable values for the driving force for nucleation of hydrates in good agreement with other thermodynamic routes used. It is shown that the driving force for hydrate nucleation at 400 bar is larger for the methane hydrate than for the carbon dioxide hydrate when compared at the same supercooling. We have also analyzed and discussed the effect of the cutoff distance of dispersive interactions and the occupancy of CO2 on the driving force for nucleation of the hydrate. [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]

Freezing of water is the most common liquid-to-crystal phase transition on Earth; however, despite its critical implications on climate change and cryopreservation among other disciplines, its characterization through experimental and computational techniques remains elusive. In this work, we make use of computer simulations to measure the nucleation rate (J) of water at normal pressure under different supercooling conditions, ranging from 215 to 240 K. We employ two different water models: mW, a coarse-grained potential for water, and TIP4P/ICE, an atomistic nonpolarizable water model that provides one of the most accurate representations of the different ice phases. To evaluate J, we apply the Lattice Mold technique, a computational method based on the use of molds to induce the nucleus formation from the metastable liquid under conditions at which observing spontaneous nucleation would be unfeasible. With this method, we obtain estimates of the nucleation rate for ice Ih and Ic and a stacking mixture of ice Ih/Ic, reaching consensus with most of the previously reported rates, although differing with some others. Furthermore, we confirm that the predicted nucleation rates obtained by the TIP4P/ICE model are in better agreement with experimental data than those obtained through the mW potential. Taken together, our study provides a reliable methodology to measure nucleation rates in a simple and computationally efficient manner that contributes to benchmarking the freezing behavior of two popular water models. [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]

Salt aqueous solutions are relevant in many fields, ranging from biological systems to seawater. Thus, the availability of a force-field that is able to reproduce the thermodynamic and dynamic behavior of salt aqueous solutions would be of great interest. Unfortunately, this has been proven challenging, and most of the existing force-fields fail to reproduce much of their behavior. In particular, the diffusion of water or the salt solubility are often not well reproduced by most of the existing force-fields. Recently, the Madrid-2019 model was proposed, and it was shown that this force-field, which uses the TIP4P/2005 model for water and non-integer charges for the ions, provides a good description of a large number of properties, including the solution densities, viscosities, and the diffusion of water. In this work, we assess the performance of this force-field on the evaluation of the freezing point depression. Although the freezing point depression is a colligative property that at low salt concentrations depends solely on properties of pure water, a good model for the electrolytes is needed to accurately predict the freezing point depression at moderate and high salt concentrations. The coexistence line between ice and several salt aqueous solutions (NaCl, KCl, LiCl, MgCl2, and Li2SO4) up to the eutectic point is estimated from direct coexistence molecular dynamics simulations. Our results show that this force-field reproduces fairly well the experimentally measured freezing point depression with respect to pure water freezing for all the salts and at all the compositions considered. [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]

We determine, for hard spheres, the Helmholtz free energy of a liquid that contains a solid cluster as a function of the size of the solid cluster by means of the formalism of the thermodynamics of curved interfaces. This is done at the constant total number of particles, volume, and temperature. We show that under certain conditions, one may have several local minima in the free energy profile, one for the homogeneous liquid and others for the spherical, cylindrical, and planar solid clusters surrounded by liquid. The variation of the interfacial free energy with the radius of the solid cluster and the distance between equimolar and tension surfaces are inputs from simulation results of nucleation studies. This is possible because stable solid clusters in the canonical ensemble become critical in the isothermal–isobaric ensemble. At each local minimum, we find no difference in chemical potential between the phases. At local maxima, we also find equal chemical potential, albeit in this case the nucleus is unstable. Moreover, the theory allows us to describe the stable solid clusters found in simulations. Therefore, we can use it for any combination of the total number of particles, volume, and global density as long as a minimum in the Helmholtz free energy occurs. We also study under which conditions the absolute minimum in the free energy corresponds to a homogeneous liquid or to a heterogeneous system having either spherical, cylindrical, or planar geometry. This work shows that the thermodynamics of curved interfaces at equilibrium can be used to describe nucleation. [ABSTRACT FROM AUTHOR]

Jiang and co-workers use high resolution STM to investigate the reaction and self-assembly of (3,6-dibromo-9,10-phenanthrenequinone, or DBPQ) molecules on Ag (100) and Ag (110) surfaces in order to understand the mechanism of bottom-up assembly on surfaces.[31] They show that, through the inclusion of multiple functional groups within a precursor molecule, it becomes possible to fabricate new low-dimensional materials with unique chemical, physical, and electronic properties. Herbst and Fransson consider the core-valence separation approximation that is often used in the calculation of core-level spectra.[5] They show how to quantify the errors in this approximation, thereby opening the door to error-quantified predictions relevant to x-ray spectroscopy. 153(16), 164108 (2020).10.1063/5.0019557 5 M. F. Herbst and T. Fransson, "Quantifying the error of the core-valence separation approximation", J. Chem. Phys. Zhu and co-workers tackle this problem for a model system containing a 2D semiconductor heterojunction and show convincingly the efficient hot electron transfer from photoexcited MoTe SB 2 sb to WS SB 2 sb .[30] This finding provides important insight into the competition between hot electron cooling and transfer at 2D semiconductor interfaces and suggests an intriguing possibility for the exploration of hot electron devices. [Extracted from the article]

In this work, polar soft-Statistical Associating Fluid Theory (SAFT) was used in a systematic manner to quantify the influence of polar interactions on the phase equilibria, interfacial, and excess properties of binary mixtures. The theory was first validated with available molecular simulation data and then used to isolate the effect of polar interactions on the thermodynamic behavior of the mixtures by fixing the polar moment of one component while changing the polar moment of the second component from non-polar to either highly dipolar or quadrupolar, examining 15 different binary mixtures. It was determined that the type and magnitude of polar interactions have direct implications on the vapor–liquid equilibria (VLE), resulting in azeotropy for systems of either dipolar or quadrupolar fluids when mixed with non-polar or low polar strength fluids, while increasing the polar strength of one component shifts the VLE to be more ideal. Additionally, excess properties and interfacial properties such as interfacial tension, density profiles, and relative adsorption at the interface were also examined, establishing distinct enrichment in the case of mixtures with highly quadrupolar fluids. Finally, polar soft-SAFT was applied to describe the thermodynamic behavior of binary mixtures of experimental systems exhibiting various intermolecular interactions (non-polar and polar), not only demonstrating high accuracy and robustness through agreement with experimental data but also providing insights into the effect of polarity on the interfacial properties of the studied mixtures. This work proves the value of having an accurate theory for isolating the effect of polarity, especially for the design of ad hoc polar solvents. [ABSTRACT FROM AUTHOR]

The application of the Young–Laplace equation to a solid–liquid interface is considered. Computer simulations show that the pressure inside a solid cluster of hard spheres is smaller than the external pressure of the liquid (both for small and large clusters). This would suggest a negative value for the interfacial free energy. We show that in a Gibbsian description of the thermodynamics of a curved solid–liquid interface in equilibrium, the choice of the thermodynamic (rather than mechanical) pressure is required, as suggested by Tolman for the liquid–gas scenario. With this definition, the interfacial free energy is positive, and the values obtained are in excellent agreement with previous results from nucleation studies. Although, for a curved fluid–fluid interface, there is no distinction between mechanical and thermal pressures (for a sufficiently large inner phase), in the solid–liquid interface, they do not coincide, as hypothesized by Gibbs. [ABSTRACT FROM AUTHOR]

Antifreeze proteins (AFPs) are biopolymers capable of interfering with ice growth. Their antifreeze action is commonly understood considering that the AFPs, by pinning the ice surface, force the crystal–liquid interface to bend forming an ice meniscus, causing an increase in the surface free energy and resulting in a decrease in the freezing point ΔTmax. Here, we present an extensive computational study for a model protein adsorbed on a TIP4P/Ice crystal, computing ΔTmax as a function of the average distance d between AFPs, with simulations spanning over 1 µs. First, we show that the lower the d, the larger the ΔTmax. Then, we find that the water–ice–protein contact angle along the line ΔTmax(d) is always larger than 0°, and we provide a theoretical interpretation. We compute the curvature radius of the stable solid–liquid interface at a given supercooling ΔT ≤ ΔTmax, connecting it with the critical ice nucleus at ΔT. Finally, we discuss the antifreeze capability of AFPs in terms of the protein–water and protein–ice interactions. Our findings establish a unified description of the AFPs in the contest of homogeneous ice nucleation, elucidating key aspects of the antifreeze mechanisms and paving the way for the design of novel ice-controlling materials. [ABSTRACT FROM AUTHOR]

We analyze the changes in the interfacial free energy between a spherical solid cluster and a fluid due to the change of the radius of the solid. Interfacial free energies from nucleation studies using the seeding technique for four different systems, being hard spheres, Lennard-Jones, and two models of water (mW and TIP4P/ICE), were plotted as a function of the inverse of the radius of the solid cluster. In all cases, the interfacial free energy was a linear function of the inverse of the radius of the solid cluster and this is consistent with Tolman's equation. This linear behavior is shown not only in isotherms but also along isobars. The effect of curvature on the interfacial free energy is more pronounced in water, followed by hard spheres, and smaller for Lennard-Jones particles. We show that it is possible to estimate nucleation rates of Lennard-Jones particles at different pressures by using information from simple NpT simulations and taking into account the variation of the interfacial free energy with the radius of the solid cluster. Neglecting the effects of the radius on the interfacial free energy (capillarity approximation) leads to incorrect values of the nucleation rate. For the Lennard-Jones system, the homogeneous nucleation curve is not parallel to the melting curve as was found for water in previous work. This is due to the increase in the interfacial free energy along the coexistence curve as the pressure increases. This work presents a simple and relatively straightforward way to approximately estimate nucleation rates. [ABSTRACT FROM AUTHOR]

In this work, a force field for several ions in water is proposed. In particular, we consider the cations Li+, Na+, K+, Mg2+, and Ca2+ and the anions Cl− and S O 4 2 − . These ions were selected as they appear in the composition of seawater, and they are also found in biological systems. The force field proposed (denoted as Madrid-2019) is nonpolarizable, and both water molecules and sulfate anions are rigid. For water, we use the TIP4P/2005 model. The main idea behind this work is to further explore the possibility of using scaled charges for describing ionic solutions. Monovalent and divalent ions are modeled using charges of 0.85 and 1.7, respectively (in electron units). The model allows a very accurate description of the densities of the solutions up to high concentrations. It also gives good predictions of viscosities up to 3 m concentrations. Calculated structural properties are also in reasonable agreement with the experiment. We have checked that no crystallization occurred in the simulations at concentrations similar to the solubility limit. A test for ternary mixtures shows that the force field provides excellent performance at an affordable computer cost. In summary, the use of scaled charges, which could be regarded as an effective and simple way of accounting for polarization (at least to a certain extend), improves the overall description of ionic systems in water. However, for purely ionic systems, scaled charges will not adequately describe neither the solid nor the melt. [ABSTRACT FROM AUTHOR]

The transformation of liquid water into solid ice is arguably the most important phase transition on Earth. A key aspect of such transformation is the speed with which ice grows once it is nucleated. There are contradictory experimental results as to whether the ice growth rate shows a maximum on cooling. Previous simulation results point to the existence of such a maximum. However, simulations were performed at constant temperature with the aid of a thermostat that dissipates the heat released at the ice-water interface unrealistically fast. Here, we perform simulations of ice growth without any thermostat. Large systems are required to perform these simulations at constant overall thermodynamic conditions (pressure and temperature). We obtain the same growth rate as in previous thermostatted simulations. This implies that the dynamics of ice growth is not affected by heat dissipation. Our results strongly support the experiments predicting the existence of a maximum in the ice growth rate. By using the Wilson-Frenkel kinetic theory, we argue that such maximum is due to a competition between an increasing crystallization thermodynamic driving force and a decreasing molecular mobility on cooling. [ABSTRACT FROM AUTHOR]

The phase diagram of molecular or colloidal systems depends strongly on the range and angular dependence of the interactions between the constituent particles. For instance, it is well known that the critical density of particles with "patchy" interactions shifts to lower values as the number of patches is decreased [see Bianchi et al. Phys. Rev. Lett. 97, 168301 (2006)]. Here, we present simulations that show that the phase behavior of patchy particles is even more interesting than had been appreciated. In particular, we find that, upon cooling below the critical point, the width of the liquid-vapor coexistence region of a system of particles with tetrahedrally arranged patches first increases, then decreases, and finally increases again. In other words, this system exhibits a doubly re-entrant liquid-vapor transition. As a consequence, the system exhibits a very large deviation from the law of rectilinear diameter, which assumes that the critical density can be obtained by linear extrapolation of the averages of the densities of the coexisting liquid and vapor phases. We argue that the unusual behavior of this system has the same origin as the density maximum in liquid water and is not captured by the Wertheim theory. The Wertheim theory also cannot account for our observation that the phase diagram of particles with three patches depends strongly on the geometrical distribution of the patches and on the degree to which their position on the particle surface is rigidly constrained. However, the phase diagram is less sensitive to small angular spreads in the patch locations. We argue that the phase behavior reported in this paper should be observable in experiments on patchy colloids and may be relevant for the liquid-liquid equilibrium in solutions of properly functionalized dendrimers. [ABSTRACT FROM AUTHOR]

We present a numerical study of the relative solubility of cholesterol in octanol and water. Our calculations allow us to compare the accuracy of the computed values of the excess chemical potential of cholesterol for several widely used water models (SPC, TIP3P, and TIP4P). We compute the excess solvation free energies by means of a cavity-based method [L. Li et al., J. Chem. Phys. 146(21), 214110 (2017)] which allows for the calculation of the excess chemical potential of a large molecule in a dense solvent phase. For the calculation of the relative solubility ("partition coefficient," log10 P o / w ) of cholesterol between octanol and water, we use the OPLS/AA force field in combination with the SPC, TIP3P, and TIP4P water models. For all water models studied, our results reproduce the experimental observation that cholesterol is less soluble in water than in octanol. While the experimental value for the partition coefficient is log10 P o / w = 3.7, SPC, TIP3P, and TIP4P give us a value of log10 P o / w = 4.5, 4.6, and 2.9, respectively. Therefore, although the results for the studied water models in combination with the OPLS/AA force field are acceptable, further work to improve the accuracy of current force fields is needed. [ABSTRACT FROM AUTHOR]

This work reexamines seeded simulation results for NaCl nucleation from a supersaturated aqueous solution at 298.15 K and 1 bar pressure. We present a linear regression approach for analyzing seeded simulation data that provides both nucleation rates and uncertainty estimates. Our results show that rates obtained from seeded simulations rely critically on a precise driving force for the model system. The driving force vs. solute concentration curve need not exactly reproduce that of the real system, but it should accurately describe the thermodynamic properties of the model system. We also show that rate estimates depend strongly on the nucleus size metric. We show that the rate estimates systematically increase as more stringent local order parameters are used to count members of a cluster and provide tentative suggestions for appropriate clustering criteria. [ABSTRACT FROM AUTHOR]

We use computer simulations to investigate the effect of salt on homogeneous ice nucleation. The melting point of the employed solution model was obtained both by direct coexistence simulations and by thermodynamic integration from previous calculations of the water chemical potential. Using a seeding approach, in which we simulate ice seeds embedded in a supercooled aqueous solution, we compute the nucleation rate as a function of temperature for a 1.85 NaCl mol per water kilogram solution at 1 bar. To improve the accuracy and reliability of our calculations, we combine seeding with the direct computation of the ice-solution interfacial free energy at coexistence using the Mold Integration method. We compare the results with previous simulation work on pure water to understand the effect caused by the solute. The model captures the experimental trend that the nucleation rate at a given supercooling decreases when adding salt. Despite the fact that the thermodynamic driving force for ice nucleation is higher for salty water for a given supercooling, the nucleation rate slows down with salt due to a significant increase of the ice-fluid interfacial free energy. The salty water model predicts an ice nucleation rate that is in good agreement with experimental measurements, bringing confidence in the predictive ability of the model. We expect that the combination of state-of-the-art simulation methods here employed to study ice nucleation from solution will be of much use in forthcoming numerical investigations of crystallization in mixtures. [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]

Despite considerable efforts over more than two decades, our knowledge of the interactions in electrolyte solutions is not yet satisfactory. Not even one of the most simple and important aqueous solutions, NaCl(aq), escapes this assertion. A requisite for the development of a force field for any water solution is the availability of a good model for water. Despite the fact that TIP4P/2005 seems to fulfill the requirement, little work has been devoted to build a force field based on TIP4P/2005. In this work, we try to fill this gap for NaCl(aq). After unsuccessful attempts to produce accurate predictions for a wide range of properties using unity ionic charges, we decided to follow recent suggestions indicating that the charges should be scaled in the ionic solution. In this way, we have been able to develop a satisfactory non-polarizable force field for NaCl(aq). We evaluate a number of thermodynamic properties of the solution (equation of state, maximum in density, enthalpies of solution, activity coefficients, radial distribution functions, solubility, surface tension, diffusion coefficients, and viscosity). Overall the results for the solution are very good. An important achievement of our model is that it also accounts for the dynamical properties of the solution, a test for which the force fields so far proposed failed. The same is true for the solubility and for the maximum in density where the model describes the experimental results almost quantitatively. The price to pay is that the model is not so good at describing NaCl in the solid phase, although the results for several properties (density and melting temperature) are still acceptable. We conclude that the scaling of the charges improves the overall description of NaCl aqueous solutions when the polarization is not included. [ABSTRACT FROM AUTHOR]

Symmetry plays an important role in the retention or annihilation of a desired interaction Hamiltonian in NMR experiments. Here, we explore the role of symmetry in the radio-frequency interaction frame Hamiltonian of the refocused-continuous-wave (rCW) pulse scheme that leads to efficient 1H heteronuclear decoupling in solid-state NMR. It is demonstrated that anti-periodic symmetry of single-spin operators (Ix, Iy, Iz) in the interaction frame can lead to complete annihilation of the ¹H-¹H homonuclear dipolar coupling effects that induce line broadening in solid-state NMR experiments. This symmetry also plays a critical role in cancelling or minimizing the effect of ¹H chemical-shift anisotropy in the effective Hamiltonian. An analytical description based on Floquet theory is presented here along with experimental evidences to understand the decoupling efficiency of supercycled (concatenated) rCW scheme. [ABSTRACT FROM AUTHOR]

By using the seeding technique the nucleation rate for the formation of ice at room pressure will be estimated for the TIP4P/ICE model using longer runs and a smaller grid of temperatures than in the previous work. The growth rate of ice will be determined for TIP4P/ICE and for the mW model of water. Although TIP4P/ICE and mW have a similar melting point and melting enthalpy, they differ significantly in the dynamics of freezing. The nucleation rate of mW is lower than that of TIP4P/ICE due to its higher interfacial free energy. Experimental results for the nucleation rate of ice are between the predictions of these two models when obtained from the seeding technique, although closer to the predictions of TIP4P/ICE. The growth rate of ice for the mW model is four orders of magnitude larger than for TIP4P/ICE. Avrami's expression is used to estimate the crystallization time from the values of the nucleation and growth rates. For mW the minimum in the crystallization time is found at approximately 85 K below the melting point and its value is of about a few ns, in agreement with the results obtained from brute force simulations by Moore and Molinero. For the TIP4P/ICE the minimum is found at about 55 K below the melting point, but its value is about ten microseconds. This value is compatible with the minimum cooling rate required to avoid the formation of ice and obtaining a glass phase. The crossover from the nucleation controlled crystallization to the growth controlled crystallization will be discussed for systems of finite size. This crossover could explain the apparent discrepancy between the values of J obtained by different experimental groups for temperatures below 230 K and should be considered as an alternative hypothesis to the two previously suggested: internal pressure and/or surface freezing effects. A maximum in the compressibility was found for the TIP4P/ICE model in supercooled water. The relaxation time is much smaller than the crystallization time at the temperature at which this maximum occurs, so this maximum is a real thermodynamic feature of the model. At the temperature of minimum crystallization time, the crystallization time is larger than the relaxation time by just two orders of magnitude. [ABSTRACT FROM AUTHOR]

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 molal at T = 298.15Kand 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. Our results have broad implications for the determination of solubilities of molecular models of ionic systems. [ABSTRACT FROM AUTHOR]

The effect of a fraction of Hartree-Fock exchange on the calculated spin-spin coupling constants involving fluorine through a hydrogen bond is analyzed in detail. Coupling constants calculated using wavefunction methods are revisited in order to get high-level calculations using the same basis set. Accurate MCSCF results are obtained using an additive approach. These constants and their contributions are used as a reference for density functional calculations. Within the density functional theory, the Hartree-Fock exchange functional is split in short- and long-range using a modified version of the Coulomb-attenuating method with the SLYP functional as well as with the original B3LYP. Results support the difficulties for calculating hydrogen bond coupling constants using density functional methods when fluorine nuclei are involved. Coupling constants are very sensitive to the Hartree-Fock exchange and it seems that, contrary to other properties, it is important to include this exchange for short-range interactions. Best functionals are tested in two different groups of complexes: those related with anionic clusters of type [F(HF)n]- and those formed by difluoroacetylene and either one or two hydrogen fluoride molecules. [ABSTRACT FROM AUTHOR]

In this work, a new rigid-nonpolarizable model of methanol is proposed. The model has three sites, located at the same positions as those used in the OPLS model previously proposed by Jorgensen [J. Phys. Chem. 90, 1276 (1986)]. However, partial charges and the values of the Lennard-Jones parameters were modified by fitting to an adequately selected set of target properties including solid-fluid experimental data. The new model was denoted as OPLS/2016. The overall performance of this model was evaluated and compared to that obtained with other popular models of methanol using a similar test to that recently proposed for water models. In the test, a certain numerical score is given to each model. It was found that the OPLS/2016 obtained the highest score (7.4 of a maximum of 10) followed by L1 (6.6), L2 (6.4), OPLS (5.8), and H1 (3.5) models. The improvement of OPLS/2016 with respect to L1 and L2 is mainly due to an improvement in the description of fluid-solid equilibria (the melting point is only 14 K higher than the experimental value). In addition, it was found that no methanol model was able to reproduce the static dielectric constant and the isobaric heat capacity, whereas the better global performance was found for models that reproduce the vaporization enthalpy once the so-called polarization term is included. Similar conclusions were suggested previously in the analysis of water models and are confirmed here for methanol. [ABSTRACT FROM AUTHOR]

We present, based on state of the art density functional theoretic calculations, a new putative ground state (GS) for the cluster (TiO2)10, which results more than 1 eV lower in energy than all those previously reported in the literature. The geometric and electronic properties of this new cluster are discussed in detail and in comparison with the rest. We analyze the implications of the new GS in the context of recent experiments of reactivity regarding oxygen exchange with gaseous CO2 in TiO2 nanostructures, and also in connection with a recent interpretation of photoelectron spectroscopic measurements of the band gap of gas phase TiO2- clusters. [ABSTRACT FROM AUTHOR]

We present, based on state of the art density functional theoretic calculations, a new putative ground state (GS) for the cluster (TiO2)10, which results more than 1 eV lower in energy than all those previously reported in the literature. The geometric and electronic properties of this new cluster are discussed in detail and in comparison with the rest. We analyze the implications of the new GS in the context of recent experiments of reactivity regarding oxygen exchange with gaseous CO2 in TiO2 nanostructures, and also in connection with a recent interpretation of photoelectron spectroscopic measurements of the band gap of gas phase TiO2- clusters. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We present a study of homogeneous crystal nucleation from metastable fluids via the seeding technique for four different systems: mW water, Tosi-Fumi NaCl, Lennard-Jones, and Hard Spheres. Combining simulations of spherical crystal seeds embedded in the metastable fluid with classical nucleation theory, we are able to successfully describe the nucleation rate for all systems in a wide range of metastability. The crystal-fluid interfacial free energy extrapolated to coexistence conditions is also in good agreement with direct calculations of such parameter. Our results show that seeding is a powerful technique to investigate crystal nucleation. [ABSTRACT FROM AUTHOR]

The role of cubic ice, ice Ic, in the nucleation of ice from supercooled water has been widely debated in the past decade. Computer simulations can provide insightful information about the mechanism of ice nucleation at a molecular scale. In this work, we use molecular dynamics to study the competition between ice Ic and hexagonal ice, ice Ih, in the process of ice nucleation. Using a seeding approach, in which classical nucleation theory is combined with simulations of ice clusters embedded in supercooled water, we estimate the nucleation rate of ice for a pathway in which the critical nucleus has an Ic structure. Comparing our results with those previously obtained for ice Ih [Sanz et al., J. Am. Chem. Soc. 135, 15008 (2013)], we conclude that within the accuracy of our calculations both nucleation pathways have the same rate for the studied water models (TIP4P/Ice and TIP4P/2005). We examine in detail the factors that contribute to the nucleation rate and find that the chemical potential difference with the fluid, the attachment rate of particles to the cluster, and the ice-water interfacial free energy are the same within the estimated margin of error for both ice polymorphs. Furthermore, we study the morphology of the ice clusters and conclude that they have a spherical shape. [ABSTRACT FROM AUTHOR]

In this work, we calculate the crystal-fluid interfacial free energy, ycf, for the Tosi-Fumi model of NaCl using three different simulation techniques: seeding, umbrella sampling, and mold integration. The three techniques give an orientationaly averaged ycf of about 100 mJ/m². Moreover, we observe that the shape of crystalline clusters embedded in the supercooled fluid is spherical. Using the mold integration technique, we compute ycf for four different crystal orientations. The obtained interfacial free energies range from 100 to 114 mJ/m², being (100) and (111) the crystal planes with the lowest and highest ycf, respectively. Within the accuracy of our calculations, the interfacial free energy either does not depend on temperature or changes very smoothly with it. Combining the seeding technique with classical nucleation theory, we also estimate nucleation free energy barriers and nucleation rates for a wide temperature range (800-1040 K). The obtained results compare quite well with brute force calculations and with previous results obtained with umbrella sampling [Valeriani et al., J. Chem. Phys, 122, 194501 (2005)] [ABSTRACT FROM AUTHOR]

The structure and spectroscopic parameters of the most relevant [C,H,N,Zn] isomers have been studied employing high-level quantum chemical methods. For each isomer, we provide predictions for their molecular structure, thermodynamic stabilities as well as vibrational and rotational spectroscopic parameters which could eventually help in their experimental detection. In addition, we have carried out a detailed study of the bonding situations by means of a topological analysis of the electron density in the framework of the Bader's quantum theory of atoms in molecules. The analysis of the relative stabilities and spectroscopic parameters suggests two linear isomers of the neutral [C,H,N,Zn] composition, namely, cyanidehydridezinc HZnCN (¹Σ) and hydrideisocyanidezinc HZnNC (¹Σ), as possible candidates for experimental detections. For the cationic [C,H,N,Zn]+ composition, the most stable isomers are the ion-molecule complexes arising from the direct interaction of the zinc cation with either the nitrogen or carbon atom of either hydrogen cyanide or hydrogen isocyanide, namely, HCNZn+ (2Σ) and HCNZn+ (²Σ). [ABSTRACT FROM AUTHOR]

The three phase equilibrium line (hydrate-liquid water-liquid carbon dioxide) has been estimated for the water + carbon dioxide binary mixture using molecular dynamics simulation and the direct coexistence technique. Both molecules have been represented using rigid nonpolarizable models. TIP4P/2005 and TIP4P/Ice were used for the case of water, while carbon dioxide was considered as a three center linear molecule with the parameterizations of MSM, EPM2, TraPPE, and ZD. The influence of the initial guest occupancy fraction on the hydrate stability has been analyzed first in order to determine the optimal starting configuration for the simulations, paying attention to the influence of the two different cells existing in the sI hydrate structure. The three phase coexistence temperature was then determined for a pressure range from 2 to 500 MPa. The qualitative shape of the equilibrium curve estimated is correct, including the high pressure temperature maximum that determines the hydrate re-entrant behaviour. However, in order to obtain quantitative agreement with experimental results, a positive deviation from the classical Lorentz-Berthelot combining rules must be considered. [ABSTRACT FROM AUTHOR]

In this work, we evaluate by means of computer simulations the rate for ice homogeneous nucleation for several water models such as TIP4P, TIP4P/2005,TIP4P/ICE, and mW (following the same procedure as in Sanz et al. [J. Am. Chem. Soc.135, 15008 (2013)]) in a broad temperature range. We estimate the ice-liquid interfacial free-energy, and conclude that for all water models γ decreases as the temperature decreases. Extrapolating our results to the melting temperature, we obtain a value of the interfacial free-energy between 25 and 32 mN/m in reasonable agreement with the reported experimental values. Moreover, we observe that the values of γ depend on the chosen water model and this is a key factor when numerically evaluating nucleation rates, given that the kinetic prefactor is quite similar for all water models with the exception of the mW (due to the absence of hydrogens). Somewhat surprisingly the estimates of the nucleation rates found in this work for TIP4P/2005 are slightly higher than those of the mW model, even though the former has explicit hydrogens. Our results suggest that it may be possible to observe in computer simulations spontaneous crystallization of TIP4P/2005 at about 60 K below the melting point. [ABSTRACT FROM AUTHOR]

In this work, we study the breakage process of an unknotted three-arm star-shaped polymer when it is pulled from its free ends by a constant force. The star polymer configuration is described through an array of monomers coupled by anharmonic bonds, while the rupture process is tracked in threedimensional space by means of Langevin Molecular Dynamics simulations. The interaction between monomers is described by a Morse potential, while a Weeks-Chandler-Anderson energetic contribution accounts for the excluded volume interaction. We explore the effect of the molecular architecture on the distributions of rupture times over a broad interval of pulling forces and star configurations. It was found that the rupture time distribution of the individual star arms is strongly affected by the star configuration imposed by the pulling forces and the length of the arms. We also observed that for large pulling forces the rupture time distributions resemble the dominant features observed for linear polymer chains. The model introduced here provides the basic ingredients to describe the effects of tensile forces on stress-induced degradation of branched macromolecules and polymer networks. [ABSTRACT FROM AUTHOR]