Further, we ensure equilibration at isothermal-isobaric (NPT) ensemble by attaining a steady pressure of ~1 atm using the Parrinello-Rahman barostat , considering a pressure relaxation time of 1ps. We found that this period of equilibration is sufficient for the convergence of thermodynamic parameters like pressure, temperature. Short-range Lennard-Jones interactions are accounted for using the minimum image convention. To estimate the non-bonding interactions including electrostatic as well as van der Waals interactions, a spherical cut-off distance of 1 nm is taken into account. In all three directions, periodic boundary conditions have been used to remove the edge effects. SHAKE algorithm is incorporated to constrain the bonds involving the hydrogen atom of the water molecules. Long-range electrostatic interactions are calculated using the particle mesh Ewald (PME) method. The production run time considered for all ternary mixture systems is 100ns (for bulk water it is 50ns) with a time step of 2fs and the trajectory is saved at every 2ps for analysis purpose. To extract a wide range of microscopic structural properties of the solvent, in-built modules of GROMACS 5.0.5 and some in-house scripts are used. For visualization purpose, we have used VMD 1.9.3. 88