Experimental and Computational Study of the Properties of Imidazole Compounds with Branched and Cycloalkyl Substituents

N-functionalized imidazole compounds with linear alkyl groups have been widely utilized precursors for imidazolium ionic liquids (ILs) while the effects of branched and cycloalkyl substituents on properties of imidazole compounds have not been studied; however, such compounds are just as synthetically accessible as those with linear alkyl groups. In this work, two fundamental properties, density and viscosity, of selected N-functionalized imidazoles bearing iso-propyl, iso-butyl, sec-butyl methylcyclopropyl, cyclopentyl, and methylcyclohexyl groups have been measured in the temperature range of 293.15–353.15 K for the guidance of molecular design for future applications. A linear and parabolic model were used for temperature-density correlation while temperature dependence of viscosity was summarized using the Andrade Equation and the Vogel-Fulcher-Tammann equation. In addition to experimental data, density, viscosity, vapor pressure and vaporization of enthalpies of target imidazole compounds were predicted using COSMOtherm calculations and compared with experimental data. It was found that the calculated densities were quite close to the experimental data, while viscosity data, obtained from COSMOtherm, underestimated experimental measurements and a scaling factor provided agreement with experiments. Predictions of vapor pressure were relatively reliable at low temperature, although the difference between experiment and prediction tended to expand with increasing temperature. Variances of vaporization enthalpies were small upon temperature change and a maximum error of ~12.3% was observed for all compounds studied.


Introduction
Ionic liquids (ILs), a well-known class of molten salts, have attracted great research interest and have been the subject of much study during the last 20+ years owing to their unique combinations of properties. Among the cations that can be used for IL design, the 1-n-alkyl-3-methylimidazolium ([C n mim]) motif is a common choice and the properties of these ILs relating to the length of the n-alkyl substituent have been thoroughly measured and modeled. Imidazolium ILs have been used as substitutes for conventional organic solvents and have been successfully used in applications including polymer solvation [1], gas separation [2][3][4], metal extraction [5,6], biomaterial preparation [7][8][9] and energy storage including fuel cells [10,11], Li+-ion batteries [11,12] and supercapacitors [13].
While imidazolium ILs are tunable through modification of the cation and/or anion moieties, the customizability of the imidazolium cation provides many more possibilities, compared to a relatively limited reservoir of candidate anions. Imidazole is an excellent building block which allows convenient grafting of functional groups to the N atoms of the five-membered ring [14], allowing for great versatility in the design of N-functionalized imidazoles and corresponding imidazolium ILs. Moreover, N-functionalized imidazoles have a plethora of other uses such as metal extraction [15,16], gas separation [17], temperature swing solvent extraction (TSSE) [18] and can provide correlations for prediction on densities of corresponding ILs [19].
In the family of N-functionalized imidazoles, most research efforts have focused on understanding the properties of N-alkylimidazole compounds with linear/straight chains including density [20,21], viscosity [22], vapor pressure [23][24][25][26], heat capacity [27], and pK a [28]. However, the properties of N-alkylimidazoles with branched or cycloalkyl functional groups have been rarely studied, although the synthesis of such compounds is no more difficult or expensive than the linear alkyl counterparts. The study of Nalkylimidazoles with branched or cycloalkyl groups will further understanding of imidazoles and corresponding ILs.
In this work, three branched and three cycloalkyl imidazole compounds were studied as shown in Table 1. Two fundamental properties, density and viscosity, were measured in temperature range of 293.15-353.15 K and equations for correlation between temperature and corresponding properties were derived from experimental data for extrapolation and guidance for molecular design of neutral alkylimidazole analogues. In addition to experimental data, COSMOtherm calculations were conducted and data from different sources were compared. Liquids 2022, 2, FOR PEER REVIEW 2 the five-membered ring [14], allowing for great versatility in the design of N-functionalized imidazoles and corresponding imidazolium ILs. Moreover, N-functionalized imidazoles have a plethora of other uses such as metal extraction [15,16], gas separation [17], temperature swing solvent extraction (TSSE) [18] and can provide correlations for prediction on densities of corresponding ILs [19]. In the family of N-functionalized imidazoles, most research efforts have focused on understanding the properties of N-alkylimidazole compounds with linear/straight chains including density [20,21], viscosity [22], vapor pressure [23][24][25][26], heat capacity [27], and pKa [28]. However, the properties of N-alkylimidazoles with branched or cycloalkyl functional groups have been rarely studied, although the synthesis of such compounds is no more difficult or expensive than the linear alkyl counterparts. The study of N-alkylimidazoles with branched or cycloalkyl groups will further understanding of imidazoles and corresponding ILs.
In this work, three branched and three cycloalkyl imidazole compounds were studied as shown in Table 1. Two fundamental properties, density and viscosity, were measured in temperature range of 293.15-353.15 K and equations for correlation between temperature and corresponding properties were derived from experimental data for extrapolation and guidance for molecular design of neutral alkylimidazole analogues. In addition to experimental data, COSMOtherm calculations were conducted and data from different sources were compared.

Materials
All the N-functionalized imidazole compounds studied in this work were synthesized as shown in Scheme 1, in which imidazole was deprotonated by NaOH, creating sodium imidazolate as an intermediate towards the formation of modified imidazoles via reaction with corresponding halide compounds. Detailed synthetic procedures for these compounds can be found in our prior works [29,30]. Compounds 1-6 used for characterizations were determined to contain <5 × 10 −4 mass fraction of impurities by Hewlett-Packard 5890 Series II gas chromatography with a flame ionization detector. 1 H NMR data obtained from a Bruker AVANCE 500 MHz NMR spectrometer in the Chemistry Department at the University of Alabama showed good consistency with previously published spectra and can be found in our previous publication [30] as they were from the same 2 N-isobutylimidazole 16245-89-9 C 7 H 12 N 2 Liquids 2022, 2, FOR PEER REVIEW 2 the five-membered ring [14], allowing for great versatility in the design of N-functionalized imidazoles and corresponding imidazolium ILs. Moreover, N-functionalized imidazoles have a plethora of other uses such as metal extraction [15,16], gas separation [17], temperature swing solvent extraction (TSSE) [18] and can provide correlations for prediction on densities of corresponding ILs [19]. In the family of N-functionalized imidazoles, most research efforts have focused on understanding the properties of N-alkylimidazole compounds with linear/straight chains including density [20,21], viscosity [22], vapor pressure [23][24][25][26], heat capacity [27], and pKa [28]. However, the properties of N-alkylimidazoles with branched or cycloalkyl functional groups have been rarely studied, although the synthesis of such compounds is no more difficult or expensive than the linear alkyl counterparts. The study of N-alkylimidazoles with branched or cycloalkyl groups will further understanding of imidazoles and corresponding ILs.
In this work, three branched and three cycloalkyl imidazole compounds were studied as shown in Table 1. Two fundamental properties, density and viscosity, were measured in temperature range of 293.15-353.15 K and equations for correlation between temperature and corresponding properties were derived from experimental data for extrapolation and guidance for molecular design of neutral alkylimidazole analogues. In addition to experimental data, COSMOtherm calculations were conducted and data from different sources were compared.

Materials
All the N-functionalized imidazole compounds studied in this work were synthesized as shown in Scheme 1, in which imidazole was deprotonated by NaOH, creating sodium imidazolate as an intermediate towards the formation of modified imidazoles via reaction with corresponding halide compounds. Detailed synthetic procedures for these compounds can be found in our prior works [29,30]. Compounds 1-6 used for characterizations were determined to contain <5 × 10 −4 mass fraction of impurities by Hewlett-Packard 5890 Series II gas chromatography with a flame ionization detector. 1 H NMR data obtained from a Bruker AVANCE 500 MHz NMR spectrometer in the Chemistry Department at the University of Alabama showed good consistency with previously published spectra and can be found in our previous publication [30] as they were from the same the five-membered ring [14], allowing for great versatility in the design of N-functionalized imidazoles and corresponding imidazolium ILs. Moreover, N-functionalized imidazoles have a plethora of other uses such as metal extraction [15,16], gas separation [17], temperature swing solvent extraction (TSSE) [18] and can provide correlations for prediction on densities of corresponding ILs [19]. In the family of N-functionalized imidazoles, most research efforts have focused on understanding the properties of N-alkylimidazole compounds with linear/straight chains including density [20,21], viscosity [22], vapor pressure [23][24][25][26], heat capacity [27], and pKa [28]. However, the properties of N-alkylimidazoles with branched or cycloalkyl functional groups have been rarely studied, although the synthesis of such compounds is no more difficult or expensive than the linear alkyl counterparts. The study of N-alkylimidazoles with branched or cycloalkyl groups will further understanding of imidazoles and corresponding ILs.
In this work, three branched and three cycloalkyl imidazole compounds were studied as shown in Table 1. Two fundamental properties, density and viscosity, were measured in temperature range of 293.15-353.15 K and equations for correlation between temperature and corresponding properties were derived from experimental data for extrapolation and guidance for molecular design of neutral alkylimidazole analogues. In addition to experimental data, COSMOtherm calculations were conducted and data from different sources were compared.

Materials
All the N-functionalized imidazole compounds studied in this work were synthesized as shown in Scheme 1, in which imidazole was deprotonated by NaOH, creating sodium imidazolate as an intermediate towards the formation of modified imidazoles via reaction with corresponding halide compounds. Detailed synthetic procedures for these compounds can be found in our prior works [29,30]. Compounds 1-6 used for characterizations were determined to contain <5 × 10 −4 mass fraction of impurities by Hewlett-Packard 5890 Series II gas chromatography with a flame ionization detector. 1 H NMR data obtained from a Bruker AVANCE 500 MHz NMR spectrometer in the Chemistry Department at the University of Alabama showed good consistency with previously published spectra and can be found in our previous publication [30] as they were from the same 4 N-cyclopropylmethylimidazole 717908-74-2 C 7 H 10 N 2 Liquids 2022, 2, FOR PEER REVIEW 2 the five-membered ring [14], allowing for great versatility in the design of N-functionalized imidazoles and corresponding imidazolium ILs. Moreover, N-functionalized imidazoles have a plethora of other uses such as metal extraction [15,16], gas separation [17], temperature swing solvent extraction (TSSE) [18] and can provide correlations for prediction on densities of corresponding ILs [19]. In the family of N-functionalized imidazoles, most research efforts have focused on understanding the properties of N-alkylimidazole compounds with linear/straight chains including density [20,21], viscosity [22], vapor pressure [23][24][25][26], heat capacity [27], and pKa [28]. However, the properties of N-alkylimidazoles with branched or cycloalkyl functional groups have been rarely studied, although the synthesis of such compounds is no more difficult or expensive than the linear alkyl counterparts. The study of N-alkylimidazoles with branched or cycloalkyl groups will further understanding of imidazoles and corresponding ILs.
In this work, three branched and three cycloalkyl imidazole compounds were studied as shown in Table 1. Two fundamental properties, density and viscosity, were measured in temperature range of 293.15-353.15 K and equations for correlation between temperature and corresponding properties were derived from experimental data for extrapolation and guidance for molecular design of neutral alkylimidazole analogues. In addition to experimental data, COSMOtherm calculations were conducted and data from different sources were compared.

Materials
All the N-functionalized imidazole compounds studied in this work were synthesized as shown in Scheme 1, in which imidazole was deprotonated by NaOH, creating sodium imidazolate as an intermediate towards the formation of modified imidazoles via reaction with corresponding halide compounds. Detailed synthetic procedures for these compounds can be found in our prior works [29,30]. Compounds 1-6 used for characterizations were determined to contain <5 × 10 −4 mass fraction of impurities by Hewlett-Packard 5890 Series II gas chromatography with a flame ionization detector. 1 H NMR data obtained from a Bruker AVANCE 500 MHz NMR spectrometer in the Chemistry Department at the University of Alabama showed good consistency with previously published spectra and can be found in our previous publication [30] as they were from the same 5 N-cyclopentylimidazole 71614-58-9 C 8 H 12 N 2 Liquids 2022, 2, FOR PEER REVIEW 2 the five-membered ring [14], allowing for great versatility in the design of N-functionalized imidazoles and corresponding imidazolium ILs. Moreover, N-functionalized imidazoles have a plethora of other uses such as metal extraction [15,16], gas separation [17], temperature swing solvent extraction (TSSE) [18] and can provide correlations for prediction on densities of corresponding ILs [19]. In the family of N-functionalized imidazoles, most research efforts have focused on understanding the properties of N-alkylimidazole compounds with linear/straight chains including density [20,21], viscosity [22], vapor pressure [23][24][25][26], heat capacity [27], and pKa [28]. However, the properties of N-alkylimidazoles with branched or cycloalkyl functional groups have been rarely studied, although the synthesis of such compounds is no more difficult or expensive than the linear alkyl counterparts. The study of N-alkylimidazoles with branched or cycloalkyl groups will further understanding of imidazoles and corresponding ILs.
In this work, three branched and three cycloalkyl imidazole compounds were studied as shown in Table 1. Two fundamental properties, density and viscosity, were measured in temperature range of 293.15-353.15 K and equations for correlation between temperature and corresponding properties were derived from experimental data for extrapolation and guidance for molecular design of neutral alkylimidazole analogues. In addition to experimental data, COSMOtherm calculations were conducted and data from different sources were compared.

Materials
All the N-functionalized imidazole compounds studied in this work were synthesized as shown in Scheme 1, in which imidazole was deprotonated by NaOH, creating sodium imidazolate as an intermediate towards the formation of modified imidazoles via reaction with corresponding halide compounds. Detailed synthetic procedures for these compounds can be found in our prior works [29,30]. Compounds 1-6 used for characterizations were determined to contain <5 × 10 −4 mass fraction of impurities by Hewlett-Packard 5890 Series II gas chromatography with a flame ionization detector. 1 H NMR data obtained from a Bruker AVANCE 500 MHz NMR spectrometer in the Chemistry Department at the University of Alabama showed good consistency with previously published spectra and can be found in our previous publication [30] as they were from the same 6 N-cyclohexylmethylimidazole 71621-00-6 C 10 H 16 N 2 Liquids 2022, 2, FOR PEER REVIEW 2 the five-membered ring [14], allowing for great versatility in the design of N-functionalized imidazoles and corresponding imidazolium ILs. Moreover, N-functionalized imidazoles have a plethora of other uses such as metal extraction [15,16], gas separation [17], temperature swing solvent extraction (TSSE) [18] and can provide correlations for prediction on densities of corresponding ILs [19]. In the family of N-functionalized imidazoles, most research efforts have focused on understanding the properties of N-alkylimidazole compounds with linear/straight chains including density [20,21], viscosity [22], vapor pressure [23][24][25][26], heat capacity [27], and pKa [28]. However, the properties of N-alkylimidazoles with branched or cycloalkyl functional groups have been rarely studied, although the synthesis of such compounds is no more difficult or expensive than the linear alkyl counterparts. The study of N-alkylimidazoles with branched or cycloalkyl groups will further understanding of imidazoles and corresponding ILs.
In this work, three branched and three cycloalkyl imidazole compounds were studied as shown in Table 1. Two fundamental properties, density and viscosity, were measured in temperature range of 293.15-353.15 K and equations for correlation between temperature and corresponding properties were derived from experimental data for extrapolation and guidance for molecular design of neutral alkylimidazole analogues. In addition to experimental data, COSMOtherm calculations were conducted and data from different sources were compared.

Materials
All the N-functionalized imidazole compounds studied in this work were synthesized as shown in Scheme 1, in which imidazole was deprotonated by NaOH, creating sodium imidazolate as an intermediate towards the formation of modified imidazoles via reaction with corresponding halide compounds. Detailed synthetic procedures for these compounds can be found in our prior works [29,30]. Compounds 1-6 used for characterizations were determined to contain <5 × 10 −4 mass fraction of impurities by Hewlett-Packard 5890 Series II gas chromatography with a flame ionization detector. 1 H NMR data obtained from a Bruker AVANCE 500 MHz NMR spectrometer in the Chemistry Department at the University of Alabama showed good consistency with previously published spectra and can be found in our previous publication [30] as they were from the same

Materials
All the N-functionalized imidazole compounds studied in this work were synthesized as shown in Scheme 1, in which imidazole was deprotonated by NaOH, creating sodium imidazolate as an intermediate towards the formation of modified imidazoles via reaction with corresponding halide compounds. Detailed synthetic procedures for these compounds can be found in our prior works [29,30]. Compounds 1-6 used for characterizations were determined to contain <5 × 10 −4 mass fraction of impurities by Hewlett-Packard 5890 Series II gas chromatography with a flame ionization detector. 1 H NMR data obtained from a Bruker AVANCE 500 MHz NMR spectrometer in the Chemistry Department at the University of Alabama showed good consistency with previously published spectra and can be found in our previous publication [30] as they were from the same batch of reactions. All compounds were further purified by distillation under reduced pressure prior to density and viscosity measurements.
batch of reactions. All compounds were further purified by distillation under reduced pressure prior to density and viscosity measurements. Scheme 1. Synthesis of N-functionalized imidazoles using imidazolate method.

Density Measurements
Densities of imidazole compounds were measured using a Mettler Toledo DM45 Del-taRange density meter through oscillation of a U-tube induced by electromagnetic effect. The measurements were conducted under atmosphere pressure with variations compensated automatically. The glass U-tube was rinsed with acetone before and after each run, followed by air flow to guarantee the dryness of the tube. A minimum volume of 1.2 mL was required for the injection of samples to fill the U-tube. Samples were loaded at ambient temperature (~293 K) and the density meter started to collect data starting at 293.15 K, and more data were collected at different temperatures with an increment of 10 K controlled by the preset program until it reached 353.15 K. As a result, seven data points were collected for each compound.

Viscosity Measurements
Viscosities of imidazole compounds were measured using a Brookfield DV-II + Pro viscometer. Viscosity was calculated automatically based on corresponding torque value and shear rate of the "ULA" spindle. For each measurement, around 20 mL liquid was used to fill the sample chamber such that the spindle was fully submerged. Temperature was controlled by an external Brookfield TC-602P circulation bath, allowing viscosity data to be collected between 273.15 and 353.15 K as programmed. The cylinder and spindle were washed with acetone after the measurement of each sample was completed. The cylinder and spindle were then fully dried with air flow prior to starting the next run.

Simulation Method
Thermophysical properties of the studied set of N-functionalized imidazole compounds were calculated using the COSMOtherm software package (BIOVIA COS-MOtherm, Release 2020) [31], which had been shown to be a practical method of predicting the physical properties of imidazoles [20]. Due to the absence of the molecular structures of the selected N-functionalized imidazole compounds in the COSMObase library provided by the software package, the necessary COSMO files were generated based on Density Functional Theory (DFT) calculations carried out with Gaussian 16, revision A.03 [32]. The optimized geometries of the structures were obtained using the BP86 functional [33,34] and the TZVP basis set [35]. Single-point energy calculations were performed using the BP86 functional and TZVP basis set with the COSMO-RS solvation model [36,37] to obtain the necessary COSMO files. This work followed the procedure of previous studies that utilized COSMOtherm to predict the properties of imidazole-based ILs [20,[38][39][40]. Finally, we performed COSMO calculations in COSMOtherm (BIOVIA COSMOtherm, Release 2020) at the TZVP level.
Densities ( ) and dynamic viscosities ( ) were calculated for each of the N-functionalized imidazole compounds over a temperature range of 293.15-353.15 K. Enthalpy of vaporization (∆ ) and vapor pressures (P) were also calculated for each of the N-functionalized imidazole compounds over a temperature range of 283.4-323.0 K for N-iso-Scheme 1. Synthesis of N-functionalized imidazoles using imidazolate method.

Density Measurements
Densities of imidazole compounds were measured using a Mettler Toledo DM45 DeltaRange density meter through oscillation of a U-tube induced by electromagnetic effect. The measurements were conducted under atmosphere pressure with variations compensated automatically. The glass U-tube was rinsed with acetone before and after each run, followed by air flow to guarantee the dryness of the tube. A minimum volume of 1.2 mL was required for the injection of samples to fill the U-tube. Samples were loaded at ambient temperature (~293 K) and the density meter started to collect data starting at 293.15 K, and more data were collected at different temperatures with an increment of 10 K controlled by the preset program until it reached 353.15 K. As a result, seven data points were collected for each compound.

Viscosity Measurements
Viscosities of imidazole compounds were measured using a Brookfield DV-II + Pro viscometer. Viscosity was calculated automatically based on corresponding torque value and shear rate of the "ULA" spindle. For each measurement, around 20 mL liquid was used to fill the sample chamber such that the spindle was fully submerged. Temperature was controlled by an external Brookfield TC-602P circulation bath, allowing viscosity data to be collected between 273.15 and 353.15 K as programmed. The cylinder and spindle were washed with acetone after the measurement of each sample was completed. The cylinder and spindle were then fully dried with air flow prior to starting the next run.

Simulation Method
Thermophysical properties of the studied set of N-functionalized imidazole compounds were calculated using the COSMOtherm software package (BIOVIA COSMOtherm, Release 2020) [31], which had been shown to be a practical method of predicting the physical properties of imidazoles [20]. Due to the absence of the molecular structures of the selected N-functionalized imidazole compounds in the COSMObase library provided by the software package, the necessary COSMO files were generated based on Density Functional Theory (DFT) calculations carried out with Gaussian 16, revision A.03 [32]. The optimized geometries of the structures were obtained using the BP86 functional [33,34] and the TZVP basis set [35]. Single-point energy calculations were performed using the BP86 functional and TZVP basis set with the COSMO-RS solvation model [36,37] to obtain the necessary COSMO files. This work followed the procedure of previous studies that utilized COSMOtherm to predict the properties of imidazole-based ILs [20,[38][39][40]. Finally, we performed COSMO calculations in COSMOtherm (BIOVIA COSMOtherm, Release 2020) at the TZVP level.
Densities (ρ) and dynamic viscosities (η) were calculated for each of the N-functionalized imidazole compounds over a temperature range of 293. 15 where η c c is the corrected viscosity (in mPa·s), and η c is the viscosity (in mPa·s) calculated in COSMOtherm.

Densities of Branched and Cycloalkyl Imidazole Compounds
Density values of branched and cycloalkyl imidazole compounds from measurement and COSMOtherm calculation results over the temperature range of 293.15-353.15 K, along with the error between computation and measurement, are shown in Table 2. Density data for N-cyclohexylmethylimidazole were not collected as this compound is solid at room temperature and could not be injected into the densitometer. As expected, density decreases with increasing temperature, which inversely correlates to volume expansion caused by increased molecular motion. For branched imidazole compounds 1, 2 and 3, the isopropyl group is denser than the iso/sec butyl group. To explain this phenomenon, the addition of carbon from the isopropyl to iso/sec-butyl group is regarded as a volume dilution of the imidazole ring which exhibits high density. However, this trend does not apply to cycloalkyl imidazole compounds in that density increases from compound 4 to compound 5 as one more carbon is incorporated to the alkyl group moiety. It is also noted that cyclization contributes to density increment as compound 4 is denser than compounds 2 and 3, although they share the same number of carbon atoms. To illustrate the effect of branching and cyclization on compound density, data for corresponding linear alkyl imidazoles published in our previous work [22] are compiled in Table S1. When comparing alkyl compounds with the same number of carbon atoms, density of imidazole compounds follows the order of branched < linear < cycloalkyl. The decrease in density from linear to branched imidazole might be attributed to the decreased packing efficiency, in which alkyl chain is reorganized and occupies more volume. This phenomenon agrees well with the trend of density change in isomers with different configurations (i.e., butanol vs. isobutanol/sec-butanol) and even in polymers (i.e., HDPE vs. LDPE). However, the density increments from linear to cycloalkyl imidazole indicate a reduction of molar volume, which is opposite to the branching effect when carbocyles are regarded as a special form of branching, where branches share both starting and end point. Thus, the effects of carbocyclic substituents should be considered. Here the increase in density of cycloalkyl imidazole can be explained as the result of reduction in degree of freedom brought by cyclization, causing the reduction of possible conformations associated with enhanced approachability between molecules and reduced intermolecular space.
Apart from experimental density data, results from COSMOtherm calculations in Table 2 showed the same trend associated with changes in molecule configuration and temperature and are close to measured values with maximum absolute percentage error as 2.64%, indicating the suitability of using COSMOtherm for density simulation.
Measured density data of imidazole compounds are presented in Figure 1, from which a clear sequence in order of magnitude can be determined. Moreover, a linear fit model was first applied to each compound for its simplicity, and corresponding fit parameters are shown in Table 3 was then utilized and corresponding parameters are shown in Table 4 Table 4. Parameters of parabolic density-temperature equation (ρ = A + B·(T − T 0 ) + C·(T − T 0 ) 2 ) derived from measurement data for imidazole compounds studied in this work within temperature range of 293.15-353.15 K.

Viscosity of Branched and Cycloalkyl Imidazole Compounds
While COSMOtherm shows good reliability in predictions on compound density, the agreement of viscosity calculations with experiments are not as close. As a result, a correction factor of 1.3 was used to better approximate experimental data from computational estimates. Viscosity values of branched and cycloalkyl imodazole compounds from measurements and corrected COSMOtherm calculations over the temperature range of 293.15-353.15 K, along with the percentage errors, are shown in Table 5 with an exception of N-cyclohexylmethylimidazole starting from 323.15 K to 353.15 K as it melts when it's heated over 318.15 K. Unscaled data from the COSMOtherm calculations are provided in Table S2.  As shown in Table 4, almost all measured viscosities were < 10 cP with the only exception of compound 6 at 323.15 K, the first viscosity data point obtained above the compound's melting point. Similar to the trend of change in density, viscosity for each compound decreases when heated due to enhanced molecular diffusion. For all six compounds shown in Table 5, viscosities of either branched or cyclized imidazoles increase upon addition of carbon atoms and this phenomenon has also been observed for linear alkyl imidazole analogues [22]. However, when comparing viscosity of linear (Table S3) and branched imidazoles, there was not a certain trend of behavior as found in density analysis. Viscosity of compound 1 is lower than its linear isomer, N-propylimidazole, across the temperature range of measurement, while compound 2 and compound 3 are of higher viscosity at low temperature and become less viscous than N-butylimidazole when heated over 343.15 K and 318.15 K, respectively. Unlike branched imidazole compounds exhibiting viscosities close to corresponding linear analogues, notable increases in viscosity are observed for the transformation from linear to cycloalkyl imidazoles. For instance, at 323.15 K, viscosity increases by 56.8% from N-butylimidazole to compound 4, 35.1% from N-pentylimidazole to compound 5 and 203.7% from N-octylimidazole to compound 6. As indicated by density, N-cycloalkylimidazole molecules share less intermolecular volume than linear alkylimidazoles, which means they are packed more closely, and thus increased van der Waals interactions result in the increment in viscosity of cycloalkylimidazoles. Apart from comparisons on viscosity values, it is found that viscosity of linear alkyl imidazoles change less with temperature. For example, viscosity of N-butylimidazole decreases by 64.1% from 293.15 to 353.15 K, while the reductions for compounds 2, 3, and 4 are 74.6%, 76.5 and 73.5%, respectively.
To better understand the tendency of compound viscosity with respect to temperature, experimental data are plotted in Figure 2. Unlike density, viscosity followed a non-linear trend. Among several empirical equations that can be used to model viscosity-temperature relationship, the Andrade Equation [41], in the form of η = a·exp·(b·T −1 ), was chosen first, for its simplicity and effectiveness in our prior work [42][43][44]. Parameters for the Andrade equation are presented in Table 6, coefficients of determination (R 2 ) for all equations derived for branched and cycloalkyl imidazoles considered in this work are >0.99 while the broad residual deviations indicated by the maximium absolute residual and SSR values call for model improvement. Consequently, the Vogel-Fulcher-Tammann equation [45] in the form of was used correlating viscosity and temperature. Fit parameters can be found in Table 7.
For each imidazole compound, the maximium absolute residual value is smaller than corresponding data in Table 6

Simulations on Vapor Pressure and Vaporization Enthalpy
Vapor pressures of branched and cycloalkyl imidazoles considered have been measured previously and enthalpies of vaporization were calculated accordingly [30]. COS-MOtherm calculations were applied in the current study with the complete results shown in Table S4 with compilation of measured values for comparison purpose. Data obtained at relatively low temperature (283-326 K) in the whole temperature range of consideration (280-350 K) were selected and presented in Table 8. It is noted that that COSMOtherm results on vapor pressure are close to measured values at lower temperature (<300 K) for branched imidazoles while variances in percentage error are much more scattered between cycloalkyl imidazoles. Moreover, as can be seen in Table S4, the error associated with vapor pressure predictions for all compounds tended to expand along with increased temperature. Unlike vapor pressure, variances in the simulation results for vaporization enthalpies were relatively consistent and acceptable, with absolute percentage error ranging from 0.84-12.31% for all compounds considered. Among these, compound 4 exhibited the best accuracy on vaporization enthalpy prediction (percentage error: 0.84-3.42%), despite a considerable gap found in vapor pressure estimations.

Conclusions
In this work, two fundamental properties, density and viscosity, of selected branched and cycloalkyl imidazole compounds were characterized in temperature range of 293.15-353.15 K. Densities of alkyl imidazole compounds bearing the same number of carbon atoms are in the order of branched < linear < cycloalkyl. It was also found that the density of branched and linear alkyl imidazoles decrease upon addition of hydrocarbon moiety while cycloalkyl imidazoles showed an opposite tendency. No distinct variance in viscosity was observed for linear and branched imidazole, but cyclic groups were found to make imidazole compounds considerably more viscous. Densities and viscosities of imidazole compounds considered were in reasonable range. Data for density were fitted into a linear and parabolic model. Viscosity data were summarized using the Andrade Equation and the Vogel-Fulcher-Tammann equation for understanding the temperature dependence of imidazole compounds on corresponding properties. Data obtained from experimental characterizations can be used for guidance of molecular design on branched and cycloalkyl imidazole analogues. Apart from property characterizations, COSMOtherm calculations applied for comparison showed excellent reliability for density simulations, while there might be more space for improvement in predictions of viscosity and vapor pressure.