A comparative study of aqueous DMSO mixtures by computer simulations and integral equation theories

Aurélien Perera,Bernarda Lovrinčević

DOI: https://doi.org/10.1080/00268976.2018.1483040

2018-01-05

Abstract:Several computer simulation studies of aqueous dimethylsulfoxyde with different force field models, and conducted by different authors, point out to an anomalous depressing of second and third neighbour correlations of the water-water radial distribution functions. This seemingly universal feature can be interpreted as the formation of linear water clusters. We test here the ability of liquid state integral equation theories to reproduce this feature. It is found that the incorporation of the water bridge diagram function is required to reproduce this feature. These theories are generally unable to properly reproduce atom-atom distribution functions. However, the near-ideal Kirkwood-Buff integrals are relatively well reproduced. We compute the Xray scattering function and compare with available experimental results, with the particular focus to explain why this data does not reproduce the cluster pre-peak observed in the water-water structure factor.

Chemical Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is about the anomalously reduced second - and third - nearest - neighbor correlations among water molecules in water - dimethyl sulfoxide (DMSO) mixtures. This phenomenon has been observed in different force - field models and computer simulation studies carried out by different authors, and can be interpreted as the formation of linear water clusters. The main objective of the paper is to test whether the Integral Equation Theories (IET) in liquid can reproduce this feature and to explore why the existing X - ray scattering data fail to show the pre - peak associated with water clusters.

### Specific problem points:
1. **Anomalously reduced second - and third - nearest - neighbor correlations among water molecules**: This phenomenon is common in computer simulations of water - DMSO mixtures and is considered a sign of the formation of linear water clusters.
2. **Application of the Integral Equation Theories in liquid**: The paper tests whether the Integral Equation Theories in liquid can reproduce the above - mentioned features, especially by introducing the water bridge diagram function to improve theoretical predictions.
3. **Interpretation of X - ray scattering data**: The paper explores why the experimentally obtained X - ray scattering data do not show the pre - peak associated with water clusters, which is inconsistent with the computer simulation results.

### Solutions:
- **Computer simulation**: Use different force - field models for computer simulation to verify the behavior of water molecules in water - DMSO mixtures.
- **Integral Equation Theories in liquid**: Apply different closure equations (such as the Hypernetted Chain (HNC), Kovalenko - Hirata (KH), and Percus - Yevick (PY) closure equations) to calculate the atom - atom distribution function and the structure factor.
- **Calculation of X - ray scattering function**: Calculate the X - ray scattering intensity through the Debye formula and compare it with the experimental data to explain why no pre - peak is observed in the experimental data.

### Main findings:
- **Computer simulation results**: The simulation results show that water molecules tend to form linear clusters due to the hydrogen - bonding interactions among water molecules and the influence of DMSO molecules.
- **Results of the Integral Equation Theories in liquid**: The HNC and KH closure equations can partially reproduce the computer simulation results in some cases, but there are still large deviations overall. In particular, the HNC closure equation cannot be solved at low concentrations, while the KH closure equation can be solved, but predicts an unreasonable increase in the small - k region.
- **Interpretation of X - ray scattering data**: The paper proposes that the reason why no pre - peak is observed in the experimental data is that the negative pre - peak of the water - DMSO cross - term cancels out the positive pre - peak of the water - water term, resulting in no obvious pre - peak in the total scattering intensity.

### Conclusions:
Although the Integral Equation Theories in liquid have limitations in describing local order, they show relatively good accuracy in calculating average properties such as Kirkwood - Buff integrals and scattering intensities. These theories are applicable to weakly aggregated systems, such as water - DMSO mixtures, but are not effective for strongly aggregated systems. The paper emphasizes the importance of many - body correlations in describing complex liquid mixtures.