Kamil Polok

Free Energy of Mixing of Acetone and Methanol: A Computer Simulation Investigation

The change of the Helmholtz free energy, internal energy, and entropy accompanying the mixing of acetone and methanol is calculated in the entire composition range by the method of thermodynamic integration using three different potential model combinations of the two compounds. In the first system, both molecules are described by the OPLS, and in the second system, both molecules are described by the original TraPPE force field, whereas in the third system a modified version of the TraPPE potential is used for acetone in combination with the original TraPPE model of methanol. The results reveal that, in contrast with the acetone-water system, all of these three model combinations are able to reproduce the full miscibility of acetone and methanol, although the thermodynamic driving force of this mixing is very small. It is also seen, in accordance with the finding of former structural analyses, that the mixing of the two components is driven by the entropy term corresponding to the ideal mixing, which is large enough to overcompensate the effect of the energy increase and entropy loss due to the interaction of the unlike components in the mixtures. Among the three model combinations, the use of the original TraPPE model of methanol and modified TraPPE model of acetone turns out to be clearly the best in this respect, as it is able to reproduce the experimental free energy, internal energy, and entropy of mixing values within 0.15 kJ/mol, 0.2 kJ/mol, and 1 J/(mol K), respectively, in the entire composition range. The success of this model combination originates from the fact that the use of the modified TraPPE model of acetone instead of the original one in these mixtures improves the reproduction of the entropy of mixing, while it retains the ability of the original model of excellently reproducing the internal energy of mixing.

A comparison of force fields for ethanol–water mixtures

Aqueous ethanol mixtures are studied through molecular dynamics simulations with the focus on exploring how various force field models reproduce the association and its influence on selected thermo-physical properties of these mixtures. The most important conclusion seems to be the inadequacy of all classical force fields to reproduce the very peculiar shape of the excess enthalpy of these mixtures, as a function of the ethanol concentration, neither quantitatively nor qualitatively. The Kirkwood–Buff (KB) integrals calculated using the simulation data follow the same trends as the experimental ones. This suggests complicated correlation of the excess enthalpy with the concentration fluctuation and clustering in these mixtures. The KB force field shows better overall agreement with experimental results than the other studied models.

Detailed insight into the hydrogen bonding interactions in acetone–methanol mixtures. A molecular dynamics simulation and Voronoi polyhedra analysis study

Voronoi polyhedra (VP) analysis of mixtures of acetone and methanol is reported on the basis of molecular dynamics computer simulations, performed at 300 K and 1 bar. The composition of the systems investigated covers the entire range from neat acetone to neat methanol. Distribution of the volume, reciprocal volume and asphericity parameter of the VP as well as that of the area of the individual VP faces and of the radius of the empty voids located between the molecules are calculated. To investigate the tendency of the like molecules to self-associate the analyses are repeated by disregarding one of the two components. The self-aggregates of the disregarded component thus turn into large empty voids, which are easily detectable in VP analysis. The obtained results reveal that both molecules show self-association, but this behavior is considerably stronger among the acetone than among the methanol molecules. The strongest self-association of the acetone and methanol molecules is found in their mole fraction ranges of 02-0.5 and 0.5-0.6, respectively. The caging effect around the methanol molecules is found to be stronger than around acetones. Finally, the local environment of the acetone molecules turns out to be more spherical than that of the methanols, not only in the respective neat liquids, but also in their mixtures.

Inhomogeneous Distribution in Methanol/Acetone Mixture: Vibrational and NMR Spectroscopy Analysis

The main aim of this paper is to quantify the inhomogeneous distribution of the components of acetone/methanol mixture and to give detailed insight into the interplay between the dipole-dipole and hydrogen bonding interactions inducing this inhomogeneity. To this end, we used the concept of infrared excess molar absorption of a given vibrational mode as an observable which contains all the information on the collective interactions in the mixture. Indeed, the changes in the infrared excess molar absorption may be associated with the inhomogeneous distribution (clustering, self-association, or high-density domains) of the components and consequently with the interaction between the two components of the mixture. The results show that acetone molecules are not homogeneously distributed in the mixture, particularly in the mole fraction range of acetone between 0.05 and 0.55. The spectral signature of this inhomogeneity is associated with the appearance of a shoulder in the C═O and C-C stretching vibrational profiles of acetone. This inhomogeneity is driven by the prevalence of the dipole-dipole interactions over those of hydrogen bonding between acetone and methanol molecules. The inhomogeneous distribution of methanol molecules is found to occur in the mole fraction range of acetone between 0.55 and 1. In this case, the hydrogen bond interactions between methanol molecules prevail over those between methanol and acetone. However, the extent of this inhomogeneity is small compared with that of acetone in the low mole fraction range. The spectral signature of this inhomogeneity is not visible in the O-H stretching vibrational mode; however, a second peak appears as a shoulder of the C-O stretching vibrational mode in this range of acetone mole fraction.

Simulations of the OKE Response in Simple Liquids Using a Polarizable and a Nonpolarizable Force Field

Recently polarizable force fields are becoming increasingly popular for molecular dynamics simulations. As the signal obtained in the optical Kerr effect (OKE) experiment is due to the polarizability dynamics of the investigated system, a study is conducted in order to compare the experimental results with those obtained with the polarizable AMOEBA force field. The comparison is made in the frequency domain; however, time domain data are also included. The selected molecular systems are the isotropic carbon tetrachloride molecule, the anisotropic chloroform, carbon disulfide and acetone molecules, and the hydrogen-bonded water and methanol molecules. Different dipole-induced-dipole (DID) method variants are used for calculation of the OKE response, showing the importance of use of the all-atom approach with preoptimized atomic polarizabilities. In order to obtain a good intermolecular to intramolecular components amplitude ratio, the isotropic polarizability in the Thole correction needs to be updated between iterations. The convergence of the spectra calculated with different DID variants is also considered, and the approach that appears to be the best gives a very good approximation after three iterations. The comparison of the experimental and simulated spectra shows a rather good agreement for the non-hydrogen-bonded molecules, although the contribution of the reorientation of anisotropic molecules is overestimated. In the case of the hydrogen-bonded molecules, the theoretical spectra are far from the experimental ones. The highly overestimated librational bands indicate excessive polarizability anisotropy introduced by the potential model. Finally, in order to verify the significance of different components of the AMOEBA model, it is gradually simplified and compared with a simple reference potential model. Removal of polarizability shows a tremendous change in the case of hydrogen-bonded liquids, whereas for the other molecules it is of minor importance. The non-hydrogen-bonded liquids are, however, more sensitive to the presence of atomic multipoles in the model.

The influence of interactions between isotopoloques on coherent, ultrafast vibrational dynamics of liquid C2Cl4

The dynamics of intramolecular and intermolecular vibrations in liquid tetrachloroethylene are studied for the first time by use of femtosecond time-resolved techniques, such as transient transmission spectroscopy and optical Kerr effect spectroscopy. Fourier transforms of time signals are compared with spontaneous Raman spectra for both isotropic and anisotropic components. The isotopic effect resulting from natural abundance of chlorine isotopes manifests itself as splitting of the isotropic spectra of intramolecular symmetric vibrations. Application of windowed Fourier transform enables us to study the dynamics of both spectral responses in real time and to analyze the role of intermolecular interactions on the coherence in the system. In order to describe the dynamics of molecules in a liquid and to explain the experimental results, we use a simple theoretical model taking into account intermolecular interactions, which allowed us to find vibrational and rotational life times.