Predictive Modeling of Thermo Physical Properties in Deep Eutectic Solvent Systems Using Jouyban–Acree and Mcallister Correlations

Liquids and liquid mixtures are central to many chemical, pharmaceutical, and industrial processes and understanding their transport properties especially viscosity is essential for design, optimization, and safety considerations. In particular, binary mixtures of alcohols often exhibit non-ideal behavior due to molecular interactions such as hydrogen bonding, dipole–dipole forces, and steric effects, making them a useful tested for studying mixing behavior and correlation modeling. Binary systems involving alcohols are of particular interest because of their strong tendency for hydrogen bonding and association–dissociation equilibria. 2-Propenol (allyl alcohol) and benzyl alcohol differ significantly in molecular size, polarity, and hydrogen-bonding capabilities. The interaction between these components often leads to complex thermodynamic behavior, reflected in deviations from ideal mixing. Both 2-propenol and benzyl alcohol possess wide industrial and pharmaceutical importance due to their physicochemical versatility. 2-propenol is commonly used as a solvent, antiseptic, and cleaning agent in pharmaceutical formulations and laboratory applications because of its volatility and bactericidal properties [1]. Benzyl alcohol, on the other hand, acts as a preservative, local anesthetic, and solvent in drug preparations, cosmetics, and perfumery products due to its low toxicity and aromatic characteristics [2]. In industrial contexts, 2-propenol serves as a raw material in paints, coatings, and fuel additives, while benzyl alcohol is utilized in resins, dye carriers, and photographic materials [3]. The binary mixture of 2-propenol and benzyl alcohol provides an excellent system for investigating non-ideal molecular interactions and thermo physical behavior, which are crucial for designing efficient solvent systems in chemical and pharmaceutical industries [4]. Such studies also aid in understanding molecular association, hydrogen bonding, and deviation from ideality, enabling optimized process control and enhanced formulation design. In continuation of our previous investigations [5-10] on binary liquid mixtures, the present study aims to analyze the thermo physical and viscous behavior of the 2-Propenol + Benzyl alcohol system at three different temperatures 298.15 K, 308.15 K, and 318.15 K and range of mole fractions. The experimental data were analyzed to evaluate the excess viscosity (η?) by six mixing predictive models, such as, Eyring [11-12], Kendall–Monroe [13], Frenkel [14], Dey [15] (Modified), Hind [16], and Wijk[17]. Excess viscosity (η?) and excess molar volume (V?), are very sensitive parameters for understanding molecular interactions and deviations from ideal behavior. The Redlich–Kister [18] polynomial was employed to correlate the excess properties with composition and to obtain polynomial coefficients that describe the strength and type of interactions between unlike molecules. Additionally, the McAllister [19] four-body and Jouyban–Acree [20-21] correlation models were applied to theoretically interpret the viscosity–composition relationship and assess their predictive accuracy through comparison with literature [22] values. The main objectives of the present work are to (i) analyze molecular interactions using excess viscosity and excess molar volume, (ii) correlate the data using the Redlich–Kister polynomial and suitable theoretical models, (iii) examine the temperature dependence of these interactions, and (iv)Identify the most accurate correlation model for predicting viscosity in hydrogen-bonded binary systems. This study provides valuable information regarding the molecular association and thermodynamic behavior of alcohol mixtures, contributing to their efficient use in industrial and pharmaceutical formulations.

Theoretical Modeling          

Kendall-Munroe:   

Hind et al: 

Frenkel relation:   

Dey et al (modified Frenkel relation):    

Wijk (WK) relation: 

Where  η1 and η2  are viscosity of pure liquids and η  is the viscosity of binary liquid mixture.

RESULTS AND DISCUSSION

The Redlich–Kister polynomial equation (7) is widely used to correlate the excess properties of binary liquid mixtures and to elucidate the nature of molecular interactions between the components. The obtained coefficients (A?–A?) for the binary mixture of 2-propenol + benzyl alcohol at three temperatures 298.15, 308.15, and 318.15 K along with the corresponding standard deviations (δ ) calculated by equation (8) are presented in Table 1.A close observation reveals that at 298.15 K, the first coefficient A? = –1.909 exhibits a relatively large negative value, indicating strong specific interactions such as hydrogen bonding and dipole–dipole forces between the component molecules. The magnitude of A? decreases progressively with increasing temperature from –1.909 to –0.736, suggesting a weakening of molecular association as thermal agitation becomes more significant at higher temperatures.

Table1. Standard deviation (δ) and Coefficient of Redlich-Kister polynomial

The subsequent coefficients (A?–A?) show both positive and negative deviations, revealing the non-ideal nature of the mixture and the presence of composition-dependent interactions. The sign changes of these coefficients imply a complex interplay between self-association of alcohol molecules and cross-interactions between 2-propenol and benzyl alcohol. The standard deviation (δ ) values remain very low (0.005–0.007) across all temperatures, confirming that the Redlich–Kister polynomial provides an excellent fit to the experimental data and accurately represents the variation of excess viscosity with mole fraction. Deviation in viscosity (?η ) were represented in form of polynomial by redlich kister for correlating measured data as;

Where A_i represent the polynomial coefficient computed by least square multiple regression method and compared in term of standard deviation (δ) by the equation given below;

Excess viscosity was computed by the following relation;

The experimental and theoretically computed values of excess viscosity (η?) for the binary mixture of 2-propenol + benzyl alcohol at temperatures 298.15, 308.15, and 318.15 K are presented in table2. The theoretical estimations were made using six predictive mixing models such as Eyring, Kendall–Monroe (KM), Frenkel (Fre), Dey (modified) relation, Hind–Ubbelhode (H–U), and Wijk (WK) equations were used to calculate the theoretical values of viscosity over the entire range of mole fractions at different temperature. At all studied temperatures, the experimental η? values are negative throughout the composition range, reaching a minimum near the equimolar region (X? ≈ 0.45–0.55). This negative deviation indicates structure-breaking interactions between unlike molecules, where the dipole–dipole and hydrogen-bonding forces between 2-propenol and benzyl alcohol molecules are weaker than the self-associative interactions present in the pure components. Such behavior is characteristic of systems where dispersion forces dominate over specific interactions, leading to a decrease in the flow resistance of the mixture compared to the ideal behavior. Clear temperature dependence is observed in the data. As the temperature increases from 298.15 K to 318.15 K, the magnitude of negative η? decreases from–0.48 to –0.20, signifying a weakening of associative forces with rising temperature. This trend can be attributed to thermal agitation disrupting intermolecular hydrogen bonds, resulting in less ordered molecular structures and a reduction in non-ideality. Among the six theoretical models, the Eyring and KM relations show closer agreement with the experimental values across the composition range, particularly at higher temperatures, suggesting that these models better represent the viscous behavior of the system under study.