Can existing models qualitatively describe the mixing behavior of acetone with water?

Surface Properties of N,N-Dimethylformamide-Water Mixtures, As Seen from Computer Simulations

The liquid-vapor interface of N,N-dimethylformamide (DMF)-water mixtures, spanning the entire composition range, is investigated in detail at 298 K by molecular dynamics simulation and intrinsic surface analysis. DMF molecules are found to adsorb strongly at the liquid surface, but this adsorption extends only to the first molecular layer. Water and DMF molecules mix with each other on the molecular scale even in the surface layer; thus, no marked self-association of any of the components is seen at the liquid surface. The major surface component prefers such orientation in which the molecular dipole vector lays parallel with the macroscopic plane of the surface. On the other hand, the preferred orientation of the minor component is determined, at both ends of the composition range, by the possibility of H-bond formation with the major component. The lack of H-donating ability of DMF leads to a rapid breakup of the percolating H-bond network at the surface; due to the strong adsorption of DMF, this breakup occurs below the bulk phase DMF mole fraction of 0.03. The disruption of the surface H-bond network also accelerates the exchange of both species between the liquid surface and bulk liquid phase, although, for water, this effect becomes apparent only above a bulk phase DMF mole fraction of 0.4. H-bonds formed by a DMF and a water molecule live, on average, 25-60% longer than those formed by two water molecules at the liquid surface. A similar, but smaller (i.e., about 10-20%) difference is seen in the bulk liquid phase. The enhanced surface mobility of the molecules results in 2-6 times larger diffusion coefficient and 2-5 times shorter H-bond lifetime values at the liquid surface than in the bulk liquid phase. The diffusion of both molecules is slowed down in the presence of the other species; in the case of DMF, this effect is caused by the formation of water-DMF H-bonds, whereas for water, steric hindrances imposed by the bulky DMF neighbors are responsible for this slowing down.

Computer Simulation of the Surface of Aqueous Ionic and Surfactant Solutions

The surface of aqueous solutions of simple salts was not the main focus of scientific attention for a long while. Considerable interest in studying such systems has only emerged in the past two decades, following the pioneering finding that large halide ions, such as I^-, exhibit considerable surface affinity. Since then, a number of issues have been clarified; however, there are still several unresolved points (e.g., the effect of various salts on lateral water diffusion at the surface) in this respect. Computer simulation studies of the field have largely benefited from the appearance of intrinsic surface analysis methods, by which the particles staying right at the boundary of the two phases can be unambiguously identified. Considering complex ions instead of simple ones opens a number of interesting questions, both from the theoretical point of view and from that of the applications. Besides reviewing the state-of-the-art of intrinsic surface analysis methods as well as the most important advances and open questions concerning the surface of simple ionic solutions, we focus on two such systems in this Perspective, namely, the surface of aqueous mixtures of room temperature ionic liquids and that of ionic surfactants. In the case of the former systems, for which computer simulation studies have still scarcely been reported, we summarize the theoretical advances that could trigger such investigations, which might well be of importance also from the point of view of industrial applications. Computer simulation methods are, on the other hand, widely used in studies of the surface of surfactant solutions. Here we review the most important theoretical advances and issues to be addressed and discuss two areas of applications, namely, the inclusion of information gathered from such simulations in large scale atmospheric models and the better understanding of the airborne transmission of viruses, such as SARS-CoV-2.

Calculation of the Free Energy of Mixing as a Tool for Assessing and Improving Potential Models: The Case of the N,N-Dimethylformamide-Water System

The Helmholtz free energy, energy, and entropy of mixing of N,N-dimethylformamide (DMF) and water are calculated in the entire composition range by means of Monte Carlo computer simulations and thermodynamic integration using all possible combinations of five DMF and three widely used water models. Our results reveal that the mixing of DMF and water is highly non-ideal. Thus, in their dilute solutions, both molecules induce structural ordering of the major component, as evidenced by the concomitant decrease in the entropy. Among the 15 model combinations considered, only 4 reproduce the well-known full miscibility of DMF and water, 3 of which strongly exaggerate the thermodynamic driving force of the miscibility. Thus, the combination of the CS2 model of DMF and the TIP4P/2005 water model reproduces the properties of the DMF-water mixtures far better than the other combinations tested. Our results also reveal that moving a fractional negative charge from the N atom to the O atom of the DMF molecule, leading to the increase in its dipole moment, improves the miscibility of the model with water. Starting from the CS2 model and optimizing the charge to be moved, we propose a new model of DMF that reproduces very accurately both the Helmholtz free energy of mixing of aqueous DMF solutions in the entire composition range (when used in combination with the TIP4P/2005 water model) and also the internal energy of neat DMF.

Thermodynamics of mixing methanol with supercritical CO2 as seen from computer simulations and thermodynamic integration

The changes in extensive thermodynamic quantities, such as volume, energy, Helmholtz free energy and entropy, occurring upon mixing liquid methanol with supercritical CO₂, are calculated using Monte Carlo simulations and thermodynamic integration for all eight combinations of four methanol and two CO₂ potential models in the entire composition range at 313 K. The obtained results are also compared with experimental data whenever possible. The transition of the system from liquid to a supercritical state is found to occur at this temperature around a CO₂ mole fraction value of 0.95 with all model combinations considered. This liquid to supercritical transition is always accompanied by positive Helmholtz free energy of mixing values and, consequently, by the non-miscibility of the two components. Furthermore, both this non-miscibility around the liquid to supercritical transition and also the miscibility of the two components below this transition, in the liquid regime, are found to be primarily of the energetic rather than entropic origin; the entropy of mixing turns out to be very close to zero, and around the liquid to supercritical transition even its qualitative behaviour is strongly model dependent. Finally, it is found that the methanol expansion coefficient is not sensitive to the details of the potential models, and it is always in excellent agreement with the experimental data. On the other hand, both the volume and the energy of mixing depend strongly on the molar volume of neat CO₂ in the model being used, and in this respect the TraPPE model of CO₂ [J. J. Potoff and J. I. Siepmann, AIChE J., 2001, 47, 1676] performs considerably better than that of Zhang and Duan [Z. Zhang and Z. Duan, J. Chem. Phys., 2005, 122, 214507].

Nanoformulation-by-design: an experimental and molecular dynamics study for polymer coated drug nanoparticles

The formulation of drug compounds into nanoparticles has many potential advantages in enhancing bioavailability and improving therapeutic efficacy. However, few drug molecules will assemble into stable, well-defined nanoparticulate structures. Amphiphilic polymer coatings are able to stabilise nanoparticles, imparting defined surface properties for many possible drug delivery applications. In the present article we explore, both experimentally and in silico, a potential methodology to coat drug nanoparticles with an amphiphilic co-polymer. Monomethoxy polyethylene glycol-polycaprolactone (mPEG-b-PCL) diblock copolymers with different mPEG lengths (M w 350, 550, 750 and 2000), designed to give different levels of colloidal stability, were used to coat the surface of indomethacin nanoparticles. Polymer coating was achieved by a flow nanoprecipitation method that demonstrated excellent batch-to-batch reproducibility and resulted in nanoparticles with high drug loadings (up to 78%). At the same time, in order to understand this modified nanoprecipitation method at an atomistic level, large-scale all-atom molecular dynamics simulations were performed in parallel using the GROMOS53a6 forcefield parameters. It was observed that the mPEG-b-PCL chains act synergistically with the acetone molecules to dissolve the indomethacin nanoparticle while after the removal of the acetone molecules (mimicking the evaporation of the organic solvent) a polymer-drug nanoparticle was formed (yield 99%). This work could facilitate the development of more efficient methodologies for producing nanoparticles of hydrophobic drugs coated with amphiphilic polymers. The atomistic insight from the MD simulations in tandem with the data from the drug encapsulation experiments thus leads the way to a nanoformulation-by-design approach for therapeutic nanoparticles.

All-atom molecular dynamics study on the non-solvent induced phase separation: Thermodynamics of adding water to poly(vinylidene fluoride)/N-methyl-2-pyrrolidone solution

The change in the thermodynamics when adding water in poly(vinylidene fluoride) (PVDF)/N-methyl-2-pyrrolidone (NMP) solution is studied from all atom molecular dynamics (MD) simulations. This is done by estimating the free energy of mixing of PVDF/NMP solution with increasing volume fraction of water (ϕ_w) using an appropriately chosen thermodynamic cycle and the Bennett acceptance ratio method. The MD calculations predict the thermodynamic phase separation point of water/NMP/PVDF to be at ϕ_w = 0.08, in close agreement with the experimental cloud point measurement (ϕ_w = 0.05). Examining the enthalpic and entropic components of the free energy of mixing reveals that at low concentrations of water, the enthalpy term has the most significant contribution to the miscibility of the ternary system, whereas at higher concentrations of water, the entropy term dominates. Finally, the free energy of mixing was compared with the Flory-Huggins (FH) free energy of mixing by computing the concentration-dependent interaction parameters from MD simulations. The FH model inadequately predicted the miscibility of the PVDF solution, mainly due to its negligence of the excess entropy of mixing.

P-NIPAM in water–acetone mixtures: experiments and simulations

The lower critical solution temperature (LCST) of poly-N-isopropylacrylamide (p-NIPAM) diminishes when a small volume of acetone is added to the aqueous polymer solution, and then increases for further additions, producing a minimum at a certain acetone concentration.

Thermodynamics of Mixing Primary Alkanolamines with Water

The volume, energy, entropy, and Helmholtz free energy of mixing of the seven simplest primary alkanolamine molecules, i.e., monoethanolamine, monoisopropanolamine, 2-amino-propan-1-ol, 2-amino-butan-1-ol, 2-amino-2-methyl-propan-1-ol, 1-amino-2-methyl-propan-2-ol, and 1-amino-butan-2-ol, with water is investigated by extensive computer simulations and thermodynamic integration. To check the force field dependence of the results, all calculations are repeated with two commonly used water models, namely, SPC/E and TIP4P. The obtained results show that the thermodynamics of mixing of alkanolamines and water is largely independent from the type of the alkanolamine molecule. The Helmholtz free energy of mixing is found to be negative for all alkanolamines at every composition, in accordance with the experimentally known full miscibility of these molecules and water. This free energy decrease occurring upon mixing is found to be clearly of energetic origin, as the energy of mixing always turns out to be negative in the entire composition range, while the entropy of mixing is also negative up to high alkanolamine mole fractions. The obtained results suggest that alkanolamines form, on average, stronger hydrogen bonds with water than what is formed by two water molecules, and they induce some ordering of the hydrating water molecules both through the hydrophobic hydration of their side chains and through the strong hydrogen bonding.

Miscibility and Thermodynamics of Mixing of Different Models of Formamide and Water in Computer Simulation

The thermodynamic changes that occur upon mixing five models of formamide and three models of water, including the miscibility of these model combinations itself, is studied by performing Monte Carlo computer simulations using an appropriately chosen thermodynamic cycle and the method of thermodynamic integration. The results show that the mixing of these two components is close to the ideal mixing, as both the energy and entropy of mixing turn out to be rather close to the ideal term in the entire composition range. Concerning the energy of mixing, the OPLS/AA_mod model of formamide behaves in a qualitatively different way than the other models considered. Thus, this model results in negative, while the other ones in positive energy of mixing values in combination with all three water models considered. Experimental data supports this latter behavior. Although the Helmholtz free energy of mixing always turns out to be negative in the entire composition range, the majority of the model combinations tested either show limited miscibility, or, at least, approach the miscibility limit very closely in certain compositions. Concerning both the miscibility and the energy of mixing of these model combinations, we recommend the use of the combination of the CHARMM formamide and TIP4P water models in simulations of water-formamide mixtures.

Thermodynamics of mixing water with dimethyl sulfoxide, as seen from computer simulations

The Helmholtz free energy, energy, and entropy of mixing of eight different models of dimethyl sulfoxide (DMSO) with four widely used water models are calculated at 298 K over the entire composition range by means of thermodynamic integration along a suitably chosen thermodynamic path, and compared with experimental data. All 32 model combinations considered are able to reproduce the experimental values rather well, within RT (free energy and energy) and R (entropy) at any composition, and quite often the deviation from the experimental data is even smaller, being in the order of the uncertainty of the calculated free energy or energy, and entropy values of 0.1 kJ/mol and 0.1 J/(mol K), respectively. On the other hand, none of the model combinations considered can accurately reproduce all three experimental functions simultaneously. Furthermore, the fact that the entropy of mixing changes sign with increasing DMSO mole fraction is only reproduced by a handful of model pairs. Model combinations that (i) give the best reproduction of the experimental free energy, while still reasonably well reproducing the experimental energy and entropy of mixing, and (ii) that give the best reproduction of the experimental energy and entropy, while still reasonably well reproducing the experimental free energy of mixing, are identified.

Modeling of mixing acetone and water: how can their full miscibility be reproduced in computer simulations?

The free energy of mixing of acetone and water is calculated at 298 K by means of thermodynamic integration considering combinations of three acetone and six water potentials. The Anisotropic United Atom 4 (AUA4) and Transferable Potential for phase Equilibria (TraPPE) models of acetone are found not to be miscible with any of the six water models considered, although the free energy cost of the mixing of any of these model pairs is very small, being below the mean kinetic energy of the molecules along one degree of freedom of 0.5RT. On the other hand, the combination of the Pereyra, Asar, and Carignano (PAC) acetone and TIP5P-E water models turns out to be indeed fully miscible, and it is able to reproduce the change of the energy, entropy, and Helmholtz free energy of mixing of the two neat components very accurately (i.e., within 0.8 kJ/mol, 2.5 J/(mol K), and 0.3 kJ/mol, respectively) in the entire composition range. The obtained results also suggest that the PAC model of acetone is likely to be fully miscible with other water models, at least with SPC and TIP4P, as well.

Comment on "Can existing models quantitatively describe the mixing behavior of acetone with water" [J. Chem. Phys. 130, 124516 (2009)]

A recent publication indicated that simulations of acetone-water mixtures using the KBFF model for acetone indicate demixing at mole fractions less than 0.28 of acetone, in disagreement with experiment and two previously published studies. Here, we indicate some inconsistancies in the current study which could help to explain these differences.