The experimental thermophysical properties (densities (ρ), speed of sound (u) and refractive (nD) indices) of binary mixtures containing 1:3 mol ratio of [1-butyl-3-dimethylimidazolium chloride to ethylene glycol] [BMIM]Cl: EG], deep eutectic solvent (DES), with selected alcohols (methanol and ethanol) were measured. The measurements of the binary mixtures and the pure solvents were conducted at temperatures ranging from (293.15 to 313.15) K in the interval of 5 K and at ambient pressure using an Anton Paar densitometer. Thermodynamic properties such as excess molar volumes (VmE), intermolecular free length (Lf), isentropic compressibilities (ks), deviation in isentropic compressibilities (Δks) and deviation in refractive indices (ΔnD) were computed from the measured data of physical properties. These excess thermodynamic properties are utilized to investigate the intermolecular interactions occurring between imidazolium DES and methanol or ethanol. Geometry optimization calculations were undertaken using Density Functional Theory (DFT) to model the structure of DES with a 1:1 mol ratio employing the B3LYP-D3 functional and a basis set of 6–31G++**. The DFT results revealed that the interaction energies within the DES in the presence of solvent favoured the formation of the complex. Furthermore, interrogation of the optimized structures at molecular level showed that the interactions were mediated by the electrostatic ionic bonding with the co-solvent playing a significant role in solubility. The excess molar volumes of binary mixtures were correlated using the Lorentz-Lorenz equation, and the densities and refractive indices were predicted using the Lorentz-Lorenz equation. The interaction energies and electrostatic interactions between the different solvents were predicted using quantum mechanics. The experimental data was found to be in good agreement with the theoretical data for the measured properties as well as the computational predictions of the stability of the systems.