Dielectric Properties of Binary Solvent Mixtures of Dimethyl Sulfoxide with Water

In this paper, the dielectric properties of water-dimethylsulfoxide (DMSO) mixtures with different mole ratios have been investigated in the range of 1 GHz to 40 GHz at 298 K by using a molecular dynamics (MD) simulation. Only one dielectric loss peak was observed in the frequency range and the relaxation in these mixtures can be described by a single relaxation time of the Davidson-Cole. It was observed that within experimental error the dielectric relaxation can be described by the Debye-like model (β ≈ 1, S.M. Puranik, et al. J. Chem. Soc. Faraday Trans. 1992, 88, 433 – 435). In general, the results are very consistent with the experimental measurements.


Introduction
Over the past few years, the use of microwave heating for promoting organic chemical transformations has been widely accepted by scientists [1,2]. The heating characteristics of a particular material via microwave irradiation are dependent on its dielectric properties. Therefore, the study of materials' dielectric properties seems a crucial need [3].
Dimethylsulfoxide (DMSO) and its mixtures with other solvents (particularly water) have aroused much interest among scientists in the last decades [4][5][6][7][8][9][10][11][12]. The dielectric behavior of supercooled aqueous solutions of DMSO was investigated by Murthy. It was found that a 2:1 complex that was present (in DMSO solutions) in the liquid state was thermolabile and existed in an undissociated state [13]. The melting temperatures were measured by using the dielectric technique in combination with a differential scanning calorimeter. The equilibrium phase diagram of DMSO was found to be eutectic, with two compounds formed by water and DMSO at the ratios of 3:1 or 2:1 [14]. Rasmussen et al. studied the phase equilibrium and non-equilibrium behaviour of solutions dissolved in water and indicated that the water-DMSO system can, in all proportions, be crystallized completely, and a stable hydrate (DMSO·3H 2 O) was formed under certain conditions [15]. The energetic and structural properties of DMSO-water clusters in gas-phase by DFT and the polarizable force field have been developed by Zhang et al. The hydrogen bond lengths, binding energies and many-body energies of a series of DMSO-water clusters were calculated using high level ab initio methods [16,17].
In fact, the long-standing interest in all mixed solvents is largely due to their importance as tunable reaction media. The dielectric constant and the relaxation time, as well as the refractive indexes and other transport properties of mixed solvents, can be conveniently tuned by changing the composition [18][19][20]. From a molecular perspective, solvent mixtures are also interesting because they often exhibit complex structural and dynamical features, especially when the components are capable of specific interactions with one another, which may result in strong interspecies molecular associations.
A dielectric study of DMSO-water using molecular dynamics simulation (MD) is introduced in this paper. The static dielectric constant ε 0 , dielectric constant at high frequency ε ∞ , the relaxation time τ, the Cole-Cole curve and the complex permittivity spectrum have been obtained.

Interaction Potentials and Simulation Details
The mixtures with DMSO mole fractions x DMSO (x D ) = 0.1, 0.21, 0.35, 0.48 and 0.91, the intermolecular potentials are described by site-site interactions given by a sum of Lennard-Jones plus Coulomb terms: ( ) ( ) 12 6 where r is the distance between sites i and j in different molecules, i q and j q are the effective charges of sites i and j, ij ε and ij σ are the energy and distance parameters in the Lennard-Jones potential. They are determined from the Lennard-Jones parameters of the different sites according to the ''mixing'' rule: where i q , i σ and i ε for the oxygen, sulfur, and methyl sites respectively in DMSO molecule and oxygen, hydrogen sites in water are shown in Table 1 [21]. Simulations reported in this paper were run in the NVT ensemble on systems consisting of a thousand molecules placed in a cubic box with periodic boundary conditions at an average temperature of 298 K. The box length was chosen to match the experimental density of DMSO-water mixtures solvent at 298 K. The Lennard-Jones forces were cut off at half of the box length and long-range Coulomb interactions were treated using the Ewald sum method [22]. The equations of motion were solved using combined SHAKE [23] and leap-frog [24] algorithms with a time step of 1 fs, during uninterrupted production runs lasted about 50 ps for each mixture. Then simulation were run in the NVE ensemble, all the simulations were extended up to 1.2 ns, where the first 50 ps were considered as equilibration. The starting solvent configurations for the five compositions of the mixtures have been chosen according to the specifications given in Table 2. x is mole fraction of solvent, i n is number of molecules of solvent i in the cubic simulation cell of edge length L ; ρ is the density at 298K, D n is the refractive index.

The static dielectric constant
To calculate the effective dielectric constant of a complex solute from computer simulation, it is necessary to map the properties of the solute observed in the simulation onto a simpler geometry, one amenable to analytical treatment. Specifically, the solute is approximated as a spherical cavity of volume V and dielectric constant ε embedded in a uniform dielectric continuum with dielectric constant ε RF . The charge distribution of the solute is represented as a point charge and point dipole placed at the center of the spherical cavity. From the fluctuations of the solute dipole moment, M, observed during a computer simulation, the absolute temperature, T, the volume, V, of the solute cavity, and the external dielectric, ε RF , it is then possible to determine the dielectric constant ε inside the solute cavity: where M is the dipole moment of the computational box, < > denotes Boltzmann ensemble averaging,   Figure 1 describes the dielectric constant ε as a function of the composition. The error bars were estimated using the blocking method, considering that the intermolecular potentials used here lack of any polarizability effects beyond the usual dipole moment enhancement and the fact that the SPC/E underestimates water's dielectric constant by about 10%, the MD values of ε for the mixtures are small than those from Ref. [2] and [27].

The relaxation time
Thermodynamic and statistical perturbation theory originally developed by Zwanzig [28] can be used to calculate relative and absolute free energies for molecularly detailed systems [29]. The change of the Gibbs free energy of a system with N particles at a given temperature, T, and pressure, P, is given by [30]: Where ΔV is a perturbation of the Hamiltonian. The composition dependence of the calculated values of G Δ is given in Figure 2. From this Figure, it can be seen that the free energy of activation increases with the concentration of the nonelectrolyte in water, reaching a maximum value of 3.5 at the concentration fraction of 0.35, whereafter it decreases. This difference can be attributed to the fact that the components are capable of specific interactions with one another. Eyring et al. have applied the theory of rate processes to molecular dipole relaxation with the following equation [31]: where G Δ is the free energy change of activation for the dipole relaxation, B k is Boltzmann's constant, h, R are Planck's constant and the molar gas constant, respectively.
The mixture's overall dielectric relaxation time τ is in good agreement with the experimental, except for compositions between 0.3 and 0.42, where the simulations predict significantly higher relaxation times seen from Figure 3. In Figure 3, the open symbols show mole fraction dependence of the relaxation time is calculated by Equation 5.

The permittivities vs. Frequency
The complex permittivities are fitted by the nonlinear least-squares fit method to the Havriliak-Negami expression to obtain various dielectric parameters: Where * ( ) ε ω is the complex permittivity at an angular frequencyω , ε ∞  is the permittivity at high frequency, ε ∞ is the static dielectric constant, α is the shape parameter representing symmetrical distribution of relaxation time and β  is the shape parameter of an asymmetric relaxation curve.
where the complex permittivity * ( ) For the Debye Equation, the real ε ′ and imaginary ε ′′ contributions to the dielectric behavior can be expressed by: Permittivity and losses for DMSO-water mixtures at 298 K are shown in Figure 4, where a single relaxation peak is observed for the entire concentration range of all the DMSO-water mixtures in the range of 1~40 GHz. Figure 4(a) obviously indicates that the higher the frequency the lower of the real ε'. Nevertheless, the real ε' decreases with the concentration of DMSO increasing, however, a completely opposite result were obtained when D x is more than 0.35. It can be addressed that the real ε' of the mixtures strongly depend on the DMSO concentration. It is expected that this phenomenon is due to the cooperative motion of DMSO-water molecules through hydrogen bonds. From Figure 4(b), the imaginary ε ′′ of dielectric constant increases with the frequency till it reaches the summit, after which the imaginary part reduce gradually with the frequency. Generally, the primary process observed for various associated liquids and those water mixtures exhibit an asymmetric loss peak. In a theoretical study, it has been suggested that the motional units in the correlated domains cooperatively move and the heterogeneity of the size distribution of the domain results in the asymmetric shape of the loss peak. The concentration dependence of DMSO-water mixture solution had been calculated by Equation 9, the Cloe-Cloe diagrams can be drawn as shown in Figure 5. It can be concluded that the relation between ε ′ and ε ′′ are circle, but those centre aren't the cross axle ε ′ .

Conclusions
This paper presents the picosecond dynamical profile of water-DMSO solutions over the entire range of composition. The static dielectric constant, dielectric relaxation in the frequency of microwave, as well as the microwave spectra of DMSO-water mixtures have been investigated through molecular dynamics simulation at room conditions, and this simulation has been proved to be an efficient tool for the study of molecular processes in solutions under the sufficiently wide frequency range.