Surface Tension and Viscosity of Binary Mixtures of the Fluorinated and Non-fluorinated Ionic Liquids [PFBMIm][PF6] and [C4C1Im][PF6] by the Pendant Drop Method and Surface Light Scattering

Mixtures of fluorinated and non-fluorinated ionic liquids (ILs) show a distinct structural organization in the bulk and at the surface. To understand how such microscopic effects influence the macroscopic bulk and surface properties of IL mixtures, knowledge of corresponding thermophysical properties including viscosity and surface tension is required yet lacking. With the intention of investigating surface enrichment effects of the fluorinated IL [PFBMIm][PF6] (3-methyl-1-(3,3,4,4,4-pentafluorobutyl)imidazolium hexafluorophosphate) in mixtures with the structurally similar, non-fluorinated IL [C4C1Im][PF6] (1-butyl-3-methylimidazolium hexafluorophosphate) observed with angle-resolved X-ray photoelectron spectroscopy (ARXPS), the pendant drop method and surface light scattering (SLS) were applied in the present study to determine surface tension and dynamic viscosity between (293 and 368) K. By adding small amounts of [PFBMIm][PF6] up to 9 mol %, a distinct increase in the viscosity and decrease in the surface tension of the mixtures relative to the properties of pure [C4C1Im][PF6] was found. This behavior reflects the nanosegregated structure in the bulk and at the surface of the binary IL mixtures. Using the results about the pronounced surface enrichment of the fluorinated chain of [PFBMIm][PF6] quantified by ARXPS, a linear mixing rule for the surface tension of the IL mixtures based on the surface tensions of the pure ILs and the surface concentration of their most surface-active groups is suggested.


Introduction
Ionic liquids (ILs) are interesting working fluids in different fields of chemistry and engineering such as electrolysis [1,2], catalysis [3,4], or separation technology [5,6]. Often, ILs consist of an inorganic anion and an organic cation with at least one alkyl chain of varying length attached to the charged head group of the cation [7]. Within the last years, ILs featuring fluorinated alkyl chains in the anion or cation, so-called fluorinated ILs, received considerable attention due to their, for example, high chemical stabilities and large solubilities for gases [8,9]. By using, e.g., mixtures of fluorinated and non-fluorinated ILs, their properties can be tailored for specific applications [10][11][12]. Therefore, it is necessary to know the thermophysical properties associated with the bulk and the surface of such IL mixtures. In this study, the viscosity and especially the surface tension as representative bulk and surface properties are of main interest. The viscosity as one key transport property describes momentum transfer [13] and the surface tension is relevant in connection with, for example, mass transport across phase boundaries and wetting [14].
To represent thermophysical properties of IL mixtures, experiments, simulations, and theoretical approaches have been used; see, e.g., Refs. [10,11,15]. In this connection, reliable data obtained from measurement techniques of proven validity are inevitable to validate and develop molecular simulations [16,17] and to establish prediction schemes [12,16,18]. Experimental thermophysical property research on IL mixtures is strongly focused on their phase equilibria including melting points and densities [15,[19][20][21][22][23][24], while viscosity [15,19,21,22,24] and in particular surface tension [12,[25][26][27] have been studied to much lesser extent. In most cases related to fluorinated ILs, only the pure fluids featuring fluorinated chains, in the anion rather than in the cation, were investigated with respect to viscosity and surface tension [8,9,28,29]. For these two properties, only a single study on binary mixtures of fluorinated and non-fluorinated ILs with a common anion is given by Merrigan et al. [25]. By adding 0.3 mass % of four different fluorinated ILs to a non-fluorinated IL, a decrease in the surface tensions of the corresponding binary IL mixtures ranging between (10 and 15) % relative to the value of the non-fluorinated IL was reported [25].
For binary mixtures of non-fluorinated ILs, the viscosity [10,11,15] and surface tension [12,26] often show a distinct non-ideal behavior as a function of composition. This is reflected in deviations of the experimental data from linear mixing rules. Deviations are found especially for the surface tension as the surface-near region can show a composition which is significantly different from that in the bulk [12,30]. Though very speculative, it is common practice that thermophysical property data of ILs or IL mixtures are used to draw conclusions about the fluid structure. Vice versa, the macroscopic thermophysical properties are governed by the molecular interactions, resulting from a specific fluid structure on a microscopic level [10,12,31].
Owing to the different types of interactions present [7,10], ILs exhibit a highly organized "nanostructure" in the bulk and at the surface. For non-fluorinated ILs, the liquid bulk features polar nanodomains consisting of the anion and the cationic head group as well as nonpolar ones made up of sufficiently long alkyl chains [10,31]. The latter are preferentially enriched at the IL surface given their relatively low cohesive energy [10,16,32]. By including fluorinated chains in the anions or cations of pure ILs, experiments [30,33] and molecular dynamics (MD) simulations [9,17,28] showed that the bulk and surface characteristics are altered. In the liquid bulk, the polar and nonpolar colloid-like domains are accompanied by nonpolar "fluorous" domains containing the fluorinated chains [9,28,33]. At the surface, a competition between the nonpolar alkyl chains in the cation and the fluorinated chains in the anion was observed [17]. For relatively short alkyl chains, the outer surface is preferentially occupied by the fluorinated chains, while the latter are gradually replaced by alkyl chains with increasing chain length [17].
To investigate whether such surface enrichment effects are also present in mixtures of fluorinated and non-fluorinated ILs, Heller et al. [30] have recently applied angle-resolved X-ray photoelectron spectroscopy (ARXPS). With this method, the composition and molecular arrangement at the surface-near region with an information depth of about (1.0 to 1.5) nm can be analyzed under ultra-high vacuum conditions; see, e.g., Refs. [30,[34][35][36]. For binary mixtures of the non-fluorinated IL [C 4 C 1 Im][PF 6 ] (1-butyl-3-methylimidazolium hexafluorophosphate) and the structurally similar fluorinated IL [PFBMIm][PF 6 ] (3-methyl-1-(3,3,4,4,4-pentafluorobutyl) imidazolium hexafluorophosphate) with a partial fluorination in the cationic alkyl chain, a pronounced enrichment of the fluorinated chain at the surface was observed [30], which is in accordance with the aforementioned findings for pure ILs [17] and with those reported for another binary mixture of fluorinated and nonfluorinated ILs [12]. The predominant surface presence of the fluorinated chain, being more distinct with decreasing [PFBMIm][PF 6 ] content, was associated with the lower surface tension of the fluorinated IL than of the non-fluorinated IL [30], an assumption that has not yet been proven experimentally.
Motivated by the aforementioned structural effects observed for IL systems containing fluorinated and non-fluorinated alkyl chains, our interest within the present work was directed to the question to which extent such effects influence the surface tension and viscosity of binary mixtures of fluorinated and non-fluorinated ILs. For answering this question, binary mixtures of the ILs [PFBMIm] [PF 6 ] and [C 4 C 1 Im][PF 6 ] were selected as model system for the experimental investigation by surface light scattering (SLS) and the pendant drop (PD) method under consistent conditions at macroscopic thermodynamic equilibrium. Here, the surface tension data obtained from the analysis of the contour of pendant drops were combined with information about the dynamics of surface fluctuations probed by SLS in order to determine the viscosity of the bulk phase. For evaluating the measurement results from both techniques relying on rigorous working equations, also density measurements were performed. To probe whether the observed surface enrichment effects are reflected in the surface tension data and in the dynamics of the surface fluctuations, we focused on the influence of small amounts of [PFBMIm][PF 6 ] up to 9 mol % on the surface tension and viscosity of the IL mixtures between (293 and 368) K. Based on a discussion of the thermophysical properties in connection with the detailed structural information available by ARXPS, a new strategy for relating the surface tension with the surface composition of the binary IL mixtures is proposed.

Materials and Sample Preparation
The non-fluorinated IL [C 4 C 1 Im][PF 6 ] (molecular weight M = 284.19 g mol −1 , melting temperature T m ≈ 265 K) was purchased from IoLiTec with a specified mole fraction purity of 0.995. For the fluorinated IL [PFBMIm][PF 6 ] (M = 374.14 g mol −1 , T m ≈ 339 K), the synthesis of the sample with an estimated mole fraction purity of about 0.990 can be found in Ref. [37]. To remove particulate contaminations from the pure IL samples, syringe filters consisting of poly(tetrafluoroethylene) and having a pore size of 450 nm were applied at temperatures of (298 and 353) K for [C 4  For the smaller of the two compositions, two different samples with similar contents of x F,B = 0.028 and 0.029 were prepared for investigations with the different techniques requiring different sample amounts, as it is discussed in the Supplementary Material. For the mixture preparation, the corresponding IL masses were added to flasks and weighed with a balance (Sartorius Entris 224i-1S) with an expanded uncertainty on a 95 % confidence level (coverage factor k = 2) of 1 mg. Considering the weighing procedure, the absolute uncertainty (k = 2) in the mole fractions x F,B of the binary mixtures can be estimated to be 0.001.
For reducing the amount of water or other volatile compounds, all samples were then degassed by an oil-sealed vacuum pump over a period of 3 h at temperatures above the melting point of the samples and at a pressure of 10 −5 MPa. For sample transfer and handling and as inert gas in the sample cells, argon purchased from Linde AG with a mole fraction purity of 0.99999 was used. Except for the pure [PFBMIm][PF 6 ], the concentration of water by mass, w H2O , was determined by Karl Fischer coulometric titration (Metrohm, 756 KF coulometer) before and after the experiments and averaged. For all water contents with an estimated relative expanded uncertainty (k = 2) of 10 %, average values ranging between 9.8·10 −4 and 1.8·10 −3 were obtained. This can be considered to be sufficiently small enough to have no impact on the properties of the IL samples [38]. Detailed information about the sources, purities, and compositions of the samples investigated is given in Table 1.
Example IL samples investigated in this work have been analyzed in ARXPS measurements concerning their purity. Within the error margins, the results matched the previous measurements reported by Heller et al. [30]. Thus, the structural results detailed therein can reliably be discussed in connection with the present investigations.

Vibrating-Tube Densimeter
The liquid density ρ required for the evaluation of the PD and SLS experiments was measured at atmospheric pressure of 0.  6 ] with a relatively large melting point of 339 K, the DMA 5000 instrument was applied between (340.65 and 363.15) K in steps of 2.5 K. For further details on the adjustment and measurement procedures, the measurement results, and the uncertainties of the density data, the reader is referred to the Supplementary Material.

Pendant Drop (PD) Method
The surface tension σ of the pure ILs and their mixtures in an argon atmosphere at about 0.11 MPa was measured by the PD method. Although the term surface tension is only connected to pure fluids consisting of a liquid and a vapor phase under saturation conditions, it is also employed for all investigated systems in the following for the sake of convenience. Besides being the key property for analyzing the effects of an enrichment of [PFBMIm][PF 6 ] at the liquid surface in mixtures with [C 4 C 1 Im] [PF 6 ], the surface tension is also required to evaluate the SLS measurements with respect to viscosity in the present case of overdamped surface fluctuations. The PD method relies on the evaluation of the contour of a pendant drop based on the Young-Laplace equation. This equation describes the profile of axisymmetric fluid bodies, such as bubbles or droplets, that are immersed in another fluid, are exposed to a gravitational field, and form a stable phase boundary [14]. In connection with the measurements for pure [C 4 6 ] as well as their mixtures at mole fractions of 0.029 and 0.090 of the fluorinated IL for temperatures between (298 and 368) K, the corresponding details on the measurement procedures and conditions are reported in the Supplementary Material. In the following, the evaluation of the surface tension data of the IL samples is described.
The calibration of the pictures with respect to the transition from the pixel to the metric scale was realized by the outer diameter of the capillary. For this, the diameter of the capillary tip visible in the pictures was measured several times with a caliper, resulting in an average value of 1.57 mm. Here, the expanded (k = 2) uncertainty is 0.03 mm, which corresponds to a relative expanded uncertainty of 1.91 %. In the analysis of the recorded pictures of the pendant drops, the Canny edge detector [39] was used to identify the droplet contour. To find suitable start values for the numerical solution of the Young-Laplace equation, the radius of a circle representing the curvature of the drop shape at the apex was estimated with a circle fit function [40]. This radius was combined with the approximation of the Bond number as given by Hansen et al. [41] and the Young-Laplace equation to estimate the surface tension of the system studied. Here, the middle of all pixels representing the lower edge of the droplet picture is assumed to be the horizontal position of the apex. For applying the Bond number approximation, the maximum width of the droplet and the droplet width at the vertical distance from the apex being equal to the maximum droplet width were determined from the detected droplet contour. In this and in the further evaluation steps, the acceleration of gravity g = 9.81 m·s −2 and the densities of the liquid and the surrounding gas phase are needed.
To determine the final surface tension value for a given droplet picture, theoretical drop profiles were compared to the experimental one. For generating suitable theoretical drop profiles, the Young-Laplace equation was solved numerically by using the 4 th order Runge-Kutta algorithm. Here, the radius representing the curvature at the apex and the surface tension value are varied independently within limited intervals around the approximation results obtained by the procedure described above. The corresponding grid featured total ranges and step sizes of (10 and 0.13) pixels for the radius representing the curvature at the apex and (8 and 0.1) mN·m −1 for σ. For each combination, the deviation of the theoretical drop profile from the experimental one was calculated considering all contour points up to the droplet height being equal to the maximum droplet diameter. Here, the theoretical drop profiles were also shifted horizontally and vertically within ± 0.5 pixels in steps of 0.25 pixels relative to the experimental profile to account for any errors of the apex location determined from the droplet contour. The surface tension in the theoretical parameter set, which represented the experimental drop profile best, was considered to be the correct one. The σ values listed in the Supplementary Material are the averages of the 3 to 5 droplets analyzed for each state. Repetition measurements performed for the mixture with x F,B = 0.090 at (338 and 299) K agreed with the initial results clearly within the expanded (k = 2) uncertainties of the surface tension data of 4 %. Furthermore, the maximum and average absolute relative deviations of the σ data obtained for the individual droplets from the reported mean values were (0.96 and 0.31) %.
Before the experiments for the IL systems, validation measurements were performed for deionized and degassed water as well as for n-dodecane in argon atmosphere at 303 K and compared with reference data. For the data evaluation, the densities of water and n-dodecane based on the publications of Wagner and Pruss [42] as well as Lemmon and Huber [43] were extracted from the REFPROP database [44]. The same holds for the reference values for the surface tension. The measured values for water and n-dodecane of (69.2 ± 0.2 and 23.5 ± 0.2) mN·m −1 showing larger and smaller σ values than the studied ILs deviated by (-2.7 and -3.8) % from the reference values [45,46]. Thus, the expanded (k = 2) uncertainties of the surface tension data for the IL systems are estimated to be 4 %. This value is considered to reflect reasonably the two main uncertainty contributions attributable to the calibration procedure and the evaluation of the droplet profile.
For the evaluation of the pendant drops of the IL systems, the correlation of our experimental density data according to Eq. 2 with the parameters given in Table 2 were used assuming negligible solubility of argon in the ILs. As the density data were not measured for exactly the same mixture compositions studied with the PD method, ideal mixing of the pure ILs was assumed to calculate the corresponding mixture densities. This includes an extrapolation of the liquid density of [PFBMIm] [PF 6 ] to temperatures below its melting point. Given that densities based on the ideal mixing rule do not deviate by more than 0.07 % from those measured for the mixture with a [PFBMIm][PF 6 ] mole fraction of 0.028, the impact of this simplification on the resulting σ data is assumed to be negligible. For the density of the argon atmosphere, the equation of state of Tegeler et al. [47] implemented in the REFPROP database [44] was employed at 0.11 MPa assuming that the vapor pressure of the ILs can be neglected. The uncertainties related to the assumptions made in the determination of density data for the evaluation of the PD experiments should be covered by the conservative estimation of the uncertainty of the obtained surface tension values.

Surface Light Scattering
With surface light scattering (SLS), the dynamics of thermal surface fluctuations present at the phase boundary of fluid systems is analyzed. In detail, the dynamics in the form of the frequency ω q and damping Γ probed at a defined modulus of scattering vector q, i.e., at a defined wavelength of the surface fluctuations is reflected by the temporal behavior of the scattered light emerging from the interaction of the incident laser light and the fluctuating surface [48,49]. For fluids with relatively small viscosities and/or large surface or interfacial tensions, an oscillatory behavior of surface fluctuations is found by examining typically investigated wavelengths of Table 2 Coefficients of Eq. 2 for the liquid density ρ calc (T) of the investigated IL systems at 0.1 MPa a Average absolute relative deviation between ρ and the correlation based on Eq. 2 fluctuations below about 10 µm [48][49][50]. This enables a simultaneous determination of viscosity and surface tension by SLS [48,[50][51][52][53][54]. Fluids with relatively large viscosities and/or small surface or interfacial tensions show an overdamped behavior of surface fluctuations (ω q = 0) [48,49], which holds also for all present measurements. In this case, the accessible mean lifetime of surface fluctuations, τ C = 1/Γ ≈ 2η/(σ q), is connected to the dynamic viscosity η and surface tension σ at a first approximation [48,55,56] and is of interest in this work. Here, combining the information about τ C with the surface tension data obtained by the PD method also enabled to evaluate the viscosity of the true bulk phase. The latter property is accessed since the penetration depth of the surface fluctuations being on the order of µm is distinctly larger than the dimension of the surface-near region, where surface enrichment is found and which exhibits extensions on the order of 1 nm or less. For both the overdamped and oscillatory cases, the SLS technique represents an absolute method requiring no calibration procedure [48] and is based on rigorous working equations for the description of the dynamics of hydrodynamic surface fluctuations [57,58]. For details about the fundamentals and methodological principles of the SLS technique, specialized literature [48,49,57,58] is recommended. The SLS setup used in this work is the same as in our former investigations of systems consisting of a vapor and liquid phase including ILs [48,51,52,55,59]. Experimental procedures and conditions for the present measurements for pure [C 4  Within the present study, heterodyne conditions, where stronger reference light is superimposed to the light scattered at the phase boundary, and a sufficient suppression of line broadening effects are fulfilled. In this case, the normalized intensity correlation function for the analysis of overdamped surface fluctuations of a defined q value takes the form [48]: In Eq. 1, the parameters a and b depend on the experimental boundary conditions. Besides the exponentially decaying mode with the characteristic decay time τ C , in theory a second such mode with a shorter decay time and negative contribution to g (2) (τ) is present [48]. Due to the typically much smaller signal strength of this second mode compared to the one shown in Eq. 1, this mode can usually not be resolved, which is also the case for the present measurements. Figure 1 shows three example normalized pseudo-cross correlation functions of the scattered light intensity obtained for the three studied IL samples with [PFBMIm][PF 6 ] mole fractions of x F,B = 0, 0.028, and 0.090 at T = 308.22 K and at a defined incident angle Θ E = 3.0°, i.e., a defined wave number q = 6.18 × 10 5 m −1 . From a nonlinear regression, the values and expanded uncertainties (k = 2) for τ C are calculated and are also given in Fig. 1. Here and in all further measurements, no systematic deviations of the measured correlation function from Eq. 1 could be observed; see the lower part of Fig. 1.
An increase in the concentration of the fluorinated IL is accompanied by an increase in the decay times of the surface fluctuations, which are clearly outside combined expanded uncertainties. At x F = 0.028, for example, τ C is about 11 % larger than the value associated with the pure non-fluorinated IL [C 4 C 1 Im][PF 6 ]. According to the first-order approximation, τ C ≈ 2η/(σ q), which is a sound representation for the present fluids due to the relatively small reduced capillary numbers Y ≪ 0.145 [48], increasing τ C values with increasing [PFBMIm][PF 6 ] mole fraction are most probably related to an increase in the viscosity and/or a decrease in the surface tension. The latter behavior can be expected taking into account that the fluorinated IL is enriched at the surface of the IL mixtures relative to the nonfluorinated one [30]. To resolve the separate influences of the liquid viscosity and surface tension on the dynamics of the surface fluctuations, it is necessary to combine the results from SLS with the surface tension data from the PD method in order to determine the viscosity.
For a reliable determination of viscosity, the dispersion relation for hydrodynamic surface fluctuations for a free liquid surface, i.e., neglecting the presence of a second fluid phase, was solved in an exact manner [57,58]. The assumption of a free liquid surface is reliably given for a system consisting of IL and argon at atmospheric pressure, where the density and viscosity of argon are distinctly smaller than the respective liquid quantities [44,55]. Furthermore, the presence of viscoelastic effects and mono or double layers at the interface can be neglected for the studied systems. Besides the directly measured data for τ C and q, reference values for the liquid density ρ and the surface tension σ are needed as input to solve the dispersion equation  [57,58]. For ρ, the same correlations according to Eq. 2 or related predictions based on the measured densities were used as it has been performed in connection with the evaluation of the PD measurements. For σ, the T-dependent correlations, Eq. 3, developed in this work were applied for the temperature states investigated by SLS. Here, extrapolations of at most 6.0 K to the lowest temperature of 293 K studied by SLS were performed. Furthermore, the correlations for σ based on the data measured at x F,B = 0.029 were used to evaluate the SLS results obtained at x F,B = 0.028 due to the negligible difference in the compositions within combined uncertainties. To represent the final η values, the averages over six independent measurements at different wave numbers q of the surface fluctuations probed were calculated. The uncertainties for η were obtained from an error propagation based on the first-order approximation for τ C . For this, the relative expanded uncertainties (k = 2) associated with the primary measured variables U r (τ C ) between (0.5 and 3.0) % and U r (q) = 0.2 % and with the reference data for the surface tension U r (σ) = 4 % were taken into account. Due to the good representation of the experimental data by the fits, the correlated σ data are assumed to have the same uncertainties as the measured ones. For the liquid dynamic viscosity η of the IL systems, a relative expanded uncertainty (k = 2) of U r (η) = 4.4 % averaged over all studied states was obtained, which mainly originates from the uncertainty of the surface tension.

Results and Discussion
The measured results for the density, viscosity, and surface tension listed in the Supplementary Material (Tables S1 to S3)

Density
The experimental data for the liquid mass densities ρ of the non-fluorinated and fluorinated ILs as well as their mixture at x F,B = 0.028 shown in the upper part of Fig. 2 were fitted as a function of temperature T by a second-order polynomial, In the fitting, all data points were considered with the same statistical weight. Table 2 gives the fit parameters ρ 0 , ρ 1 , and ρ 2 and the average absolute deviations (AADs) of the experimental ρ data from the correlations based on Eq. 2. . Here, the light hydrogen atoms in the alkyl chain of the cation are substituted by heavier but still relatively small fluorine atoms, without having a strong impact on the molecular volume of the cation which is only slightly increased [29]. Increasing densities of pure ILs  6 ]. For the density and viscosity, error bars representing the expanded uncertainties (k = 2) are within the symbols. For the surface tension, error bars of 4 % (k = 2) are only shown exemplarily for one data set related to x F,B = 0 due to visibility purposes by fluorination of the alkyl chains in the cations has also been reported by other authors [8,29,60]. A mole-based linear mixing rule for the molar volume of the mixtures, which is based on the densities of the pure substances correlated by Eq. 2 and assumes no excess molar volume, describes the correlated mixture densities for x F,B = 0.028 within 0.07 % between temperatures of (293.15 and 368.15) K. This finding that binary IL mixtures show negligible excess volume, in particular for ILs with different cations or anions of relatively similar sizes, is in agreement with literature [10,15,19,20]. Thus, the linear mixing rule can also be applied to represent the mixture densities at x F,B = 0.029 studied with the PD method and at x F,B = 0.090 where no measurements could be performed in this work.
For pure [C 4 C 1 Im][PF 6 ], the experimental density data reported in literature and mainly determined via vibrating-tube densimetry agree in most cases within combined uncertainties with our present measurement results based on Eq. 2, as can be seen in Fig. S1 of the Supplementary Material. The few large deviations between the different data sets of up to about 3 % seem not to be caused by differences in the water contents of the samples, but rather by different and/or undefined impurities, which can affect the density of ILs to this extent.

Viscosity
For the dynamic viscosity of the bulk phase of the IL systems studied, the temperature dependency of the experimental data obtained by SLS is represented by a Vogel-type equation, For the correlation, it was considered that the statistical weight of each experimental data point is the same. The fit parameters η 0 , B, and C and the AAD values of the measured η data from the correlations according to Eq. 3 are given in Table 3. In all cases, the absolute percentage residuals of the measured viscosity data from the fit, which are 0.53 % on average, are clearly smaller than the experimental expanded uncertainties.
The middle part of Fig. 2 shows that the viscosity of the three IL systems decreases strongly with increasing temperature, which is a common behavior of ILs [7]. With increasing mole fraction of the fluorinated IL up to 0.090, significantly increasing liquid viscosities are found. Comparing the results at x F,B = 0.090 with those at x F,B = 0, an increase in temperature results in smaller   6 ] − anion as well as non-fluorinated and fluorinated imidazolium-based cations of longer chain length than investigated here. This behavior also corresponds to larger viscosities found for ILs with fluorinated cations compared to their non-fluorinated analogues [8,29]. In general, the viscosity of ILs is strongly governed by intramolecular entanglement as well as by intermolecular interactions including hydrogen bonding, electrostatic Coulomb interactions, dipole-dipole interactions, and van der Waals (vdW) interactions [7,10]. Due to the complex interplay of such interactions over different spatial extensions, the bulk phase of ILs shows a distinct heterogeneity on a nanoscale in form of segregated colloid-like domains [7,10,31]. The polar domains mainly feature the anions and the charged head group of the cations as a result of attractive electrostatic interactions, while the nonpolar domains are formed by sufficiently long alkyl chains interacting predominantly via dispersive vdW forces [31].
By the introduction of fluorinated alkyl chains in ILs, a third segregated "fluorous" domain can be formed, as was observed for pure ILs in experiments [33] and MD simulations [28]. This so-called "fluorophobic" effect is caused by the exclusion of the fluorinated chains from other domains and especially from the nonpolar alkyl chains which show stronger vdW interactions between each other than to the fluorinated chains [61,62]. Among the latter, rather weak or even repulsive interactions have been reported [62]. This behavior is apparently overcompensated by stronger attractive electrostatic cation-cation interactions induced by fluorinated cations than by non-fluorinated ones in the bulk of IL systems. The strongly negatively charged fluorine atoms in the fluorinated chains allow for a favorable electrostatic interaction with positively charged atoms in the imidazolium ring [29], which may also imply weak hydrogen bonding [61]. This results in stronger cation-cation interactions in [PFBMIm][PF 6 ] than in [C 4 C 1 Im] [PF 6 ] as an additional contribution to the predominant anion-cation interactions.
In conclusion, the more pronounced electrostatic interactions between the fluorinated cations than between the non-fluorinated ones seem to result in increased viscosities for pure [ 6 ] scatter around our data correlated by Eq. 3 with relative deviations between (+79 and −32) %. These deviations and those among the different data sets are often clearly outside combined measurement uncertainties, where our data from SLS seem to form the median of the literature data. Such large discrepancies, which have also been found for other ILs [55,63,64], may be attributed to the samples, especially the presence of water [38], and to the measurement procedure, where often conventional techniques requiring a calibration procedure and/or the application of a macroscopic shear gradient are used.

Summary and Discussion of Experimental Results
For the surface tension as characteristic macroscopic property of IL surfaces, the experimental data of the different IL systems determined by the PD method are described by a linear equation as a function of temperature, Besides the fit parameters σ 0 and σ 1 , the AAD values of the measured σ data from the correlations based on Eq. 4 are given in Table 4.
As can be seen from the lower part of Fig. 2, the surface tensions of the studied IL systems show the typical trend of decreasing values with increasing temperature. For all temperatures studied, the surface tension decreases with increasing mole fraction of the fluorinated IL. This behavior is in agreement with the significantly lower surface tension of the fluorinated IL [PFBMIm][PF 6 ] (σ = 29.8 mN·m −1 at 368 K) compared to that of the non-fluorinated IL [C 4 C 1 Im][PF 6 ] (σ = 38.3 mN·m −1 at 368 K). A significant reduction in the surface tension of pure ILs by fluorination of alkyl chains in cations [12] or anions [17] has also been found for other IL types. For the same binary IL mixtures as mentioned before in connection with the viscosity, Merrigan et al. [25] observed a distinct decrease in the surface tension of the mixtures relative to that of the non-fluorinated IL by adding small amounts of the fluorinated IL.
For the surface tension of [C 4 C 1 Im][PF 6 ] shown in Fig. S3 of the Supplementary Material, most of the literature data show positive deviations from the present data represented by Eq. 4. The deviations ranging between (+28 and −2) % are, however, often still within combined expanded uncertainties of the data sets, in particular considering the majority of the data sets which feature deviations of about 5 % and less. Differences in published surface tension data of ILs present in this work and in previous studies [55,65] appear to be mainly caused by the underlying method and by impurities, out of which water as common impurity has no significant impact up to mole fractions of about 0.4 [38]. The tendency of somewhat lower σ values of the IL obtained in this study could also be seen for the reference fluids water and n-dodecane and is predominantly related to the calibration procedure and the evaluation of the droplet profile. For ILs, the surface tension is often discussed in the context of the Langmuir's principle [66]. A comprehensive analysis of this principle can be found in the corresponding literature [10,16,67,68]. In the following, only the information relevant for this study is given. Langmuir's principle [66] states that the surface tension is affected by the outer region of the surface. Since any fluid aims to reduce its surface energy and, thus, the surface tension, the segments of the IL with the lowest cohesive energy are preferentially enriched at the surface [16]. Besides this enthalpic effect, also entropic effects associated with the molecular orientation at the surface need to be considered [16]. For pure non-fluorinated ILs consisting of 1-n-alkyl-3-methylimidazolium ([C n C 1 Im] + ) cations and various types of anions, the alkyl side chains with their relatively weak vdW interactions align themselves preferentially towards the gas phase. This outer nonpolar layer is above an ionic sublayer mainly containing the charged head groups of the imidazolium rings and the anions, as could be evidenced in experiments [16,30,32] and MD simulations [16,69]. In the study of Heller et al. [30], such behavior has also been observed by ARXPS measurements for the IL [C 4 C 1 Im][PF 6 ] relevant in this work. For this IL, a slight enrichment of the relatively short butyl chain (C 4 ) and a slight depletion of the [PF 6 ] − anion is present at the outer surface region.
For pure imidazolium-based ILs with fluorinated chains in the cation [12,30] or anion [17], these chains were also found to be strongly enriched at the surface. This enrichment is even larger than the enrichment of the alkyl chains, as it could be revealed by ARXPS for [PFBMIm][PF 6 ] compared to [C 4 C 1 Im][PF 6 ] [30]. The same behavior was also observed with MD simulations for pure ILs with fluorinated chains in the anions and alkyl chains in the cation, if the latter ones are relatively short [17]. The preferential enrichment of the fluorinated chains can mainly be related to the weaker intermolecular interactions between fluorinated chains compared to those between alkyl chains. In view of Langmuir's principle [66], the surface tension of [PFBMIm] [PF 6 ] is smaller than that of [C 4 C 1 Im] [PF 6 ], which is given by the PD measurement results at 368 K.
For binary mixtures of ILs with the same anion, but with non-fluorinated and fluorinated cations, experiments evidenced that the surface is preferentially occupied by the fluorinated chains [12,30]. This includes the present mixture of [C 4 C 1 Im][PF 6 ] and [PFBMIm][PF 6 ], where ARXPS showed a clear enrichment of the fluorinated chains at the surface relative to the bulk composition, while the alkyl chains are depleted [30]. The relative enrichment of the fluorinated chains is more pronounced with decreasing mole fraction of [PFBMIm][PF 6 ] and with decreasing temperature [30]. To analyze these structural effects in connection with the behavior of the surface tension, the measured σ data are shown in Fig. 3 as a function of the bulk mole fraction x F,B of [PFBMIm][PF 6 ] for comparable temperatures. In all cases, the surface tension data show a significant decrease with increasing content of the fluorinated IL, which corresponds with the surface enrichment effects found by ARXPS qualitatively well. The somewhat larger enrichment at lower temperatures is not reflected in a more distinct decrease of the surface tension data, which is similar for the different temperatures studied.

Prediction Method for IL Mixtures
To quantify the effect of the surface enrichment on the surface tension, it is of interest to probe correlations between these microscopic and macroscopic properties of IL surfaces. In a recent approach [16], experimental ARXPS results were combined with MD simulations to represent the surface tension values of the homologous series of [C n C 1 Im] + -based ILs by contributions from different segments terminating the surface, considering their volume fractions. For binary IL mixtures being of interest here, different theoretical and empirical approaches have been discussed in the literature, as it has been reviewed by, e.g., Smoll et al. [12] and Nakajima et al. [70]. A comparison with relevant approaches is not possible for the present IL mixture due to the lack of required input information. Instead, a new, phenomenological approach based on the surface tensions of the  6 ] at comparable temperatures. For 368 K, a linear mole fraction-based mixing rule using the surface tension data of the pure ILs is also shown as dashed line. For the two IL molecules given in the graph, the same nomenclature has been applied as in Ref. [30]. Errors bars for the surface tension data are not shown due to visibility purposes pure ILs and their "effective" surface concentrations is suggested in the present work, which is shown in the following.
It is obvious that a relation between the experimental σ values and the bulk mole fractions in the IL mixtures fails. The discrepancy can be seen from the trend of the two measured surface tension results of the IL mixtures at 368 K which are systematically lower than the values predicted by a linear mixing rule based on the bulk mole fractions and the surface tensions of the pure ILs (see dashed line in Fig. 3). This non-ideal behavior has also been found by, e.g., Smoll et al. [12] for binary IL mixtures with a common anion and two non-fluorinated [C n C 1 Im] + cations of different alkyl chain lengths (n = 2 and 12). All this indicates that instead of bulk compositions, surface compositions need to be considered in appropriate prediction schemes [12,70].
With the objective of developing a relatively simple relationship between surface tension and surface composition, we make use of the detailed structural information about the surface gained by ARXPS. This technique provides the quantitative determination of the mole fractions of all elements typically found in ILs except for hydrogen. Using an emission angle of 80° with respect to the surface normal, the information depth is (1.0 to 1.5) nm, i.e., the outer surface region is probed [34,71]. By normalizing the measured content of each element in the outer surface region to the corresponding nominal content which corresponds to the bulk composition, so-called normalized contents can be calculated, which show relative enrichment (value > 1) or depletion (value < 1) in the outer surface region. For the mixtures of [C 4 C 1 Im][PF 6 ] and [PFBMIm] [PF 6 ] studied at 368.15 K, the upper part of Fig. 4 shows the normalized contents of the fluorine atoms F CFx attached to carbon atoms C CFx in the fluorinated cation (shown in green color in Fig. 3) and of the carbon atoms C alkyl in the alkyl chain of the non-fluorinated cation (shown in gray color in Fig. 3) as a function of the bulk mole fraction x F,B , as they are reported in Ref. [30]. In this work, three significant digits were used for the corresponding measured contents contributing to the normalized contents. In addition, the uncertainty of the measured contents of 10 % [30] is assigned to the normalized contents.
In line with the previous discussion, a clear surface enrichment of the fluorinated chains is found, while the alkyl chains are depleted at the surface; see the upper part of Fig. 4a. Only for the pure non-fluorinated IL, i.e., at x F,B = 0, the C alkyl group exhibits a slight enrichment. Furthermore, the normalized contents for the atoms in the [PF 6 ] − anion and in the imidazolium ring (not shown in the upper part of Fig. 4) are similar and close to 1 over the entire composition range, which indicates no distinct surface enrichment or depletion of both groups [30]. Thus, among the four relevant molecular segments in the IL mixtures, i.e., the anion, the imidazolium ring, the alkyl chain, and the fluorinated chain, it is the idea of considering only the two competing, most surface-active moieties in form of the F CFx atoms in the fluorinated group and the C alkyl group associated with the fluorinated and non-fluorinated IL, respectively. The mole fractions of these groups in such a pseudo-binary mixture account for their molar composition in the surface-near region and are considered to reflect the "effective" contributions of each associated IL to the surface tension of the IL mixture.
Based on the normalized F CFx and C alkyl surface contents measured by ARXPS under an emission angle of 80° at 368.15 K for varying x F,B values, the surface mole fraction x F,S of the fluorinated chain associated with the fluorinated IL, can be derived analytically. Given the larger normalized content of F CFx than that of C alkyl , Eq. 5 provides always larger values for x F,S than for x F,B . The relative surface exposure of the fluorinated chain can be specified by an "enrichment factor" E F = x F,S /x F,B , which shows a similar concentration-dependent trend as the normalized F CFx content and ranges between 2.9 at x F,B = 0.10 and 1.0 at x F,B = 1.0. Based on the E F values at x F,B = 0.10, 0.25, 0.50, 0.75, and 1.0 calculated from Eq. 5 with the aid of the ARXPS results from the upper part of Fig. 4 [30], a simple empirical fit according to E F = a F + b F ·x F,B −0.5 with the parameters a F = 0.10 and b F = 0.90 was performed. With this fit and Eq. 5, the surface mole fraction x F,S can be calculated as a function of the bulk mole fraction x F,B , which is shown by the solid line in the lower part of Fig. 4. As can be seen, the fit represents the underlying data points indicated by the triangles well and within the uncertainty of 10 % for the x F,S data, which is also associated with the normalized contents in Eq. 5. For relatively low x F,B values which are of main interest here, an initial strong deviation of x F,S from the nominal stoichiometric contents (see dashed diagonal line) is found. For the mixtures investigated at 368 K, the values for x F,B of 0.029 and 0.090 correspond to larger values of x F,S of 0.16 and 0.28, respectively. The compositions for the pure ILs, i.e., at x F,B = 0 and 1, are not affected because the corresponding x F,S values are the same. As a result, a surface mole fraction-based, linear mixing rule for the representation of the surface tension of the binary IL mixtures according to is suggested. This rule is based on the surface tension of the pure fluorinated (σ F ) and non-fluorinated (σ non-F ) IL and their corresponding "effective" surface mole fractions x F,S and (1 − x F,S ), and can, thus, be understood within the concept of Langmuir [66]. Figure 5 shows that the experimental surface tension data for the mixtures at 368 K represented as a function of x F,S follow the linear mole fractionbased mixing rule Eq. 6 given by the solid line well. By converting the x F,S values to The results for the surface tension and surface composition can also be discussed in the concept of the Gibbs adsorption equation [72]. This equation describes the excess of the surface-active compound, i.e., the fluorinated IL [PFBMIm][PF 6 ], per unit area of the surface relative to the amount which would be present in case the corresponding bulk concentration is extended to the surface. Assuming an infinite dilution of the fluorinated IL in the binary mixtures where its activity coefficient approaches unity, the relative excess Γ F at constant temperature T and pressure p is given by [72].
where R is the ideal gas constant. To calculate the partial derivative in Eq. 7, the model for the surface tension in Eq. 6 was combined with the fit for E F which connects x F,S with x F,B . For the two studied mole fractions x F,B of 0.029 and 0.09 at 368 K, the corresponding Γ F values are about (0.22 and 0.40) µmol·m −2 . Positive Γ F values indicate that the fluorinated IL with the lower surface tension is the surfaceactive substance and reduces the surface tension of the mixture. This is in accordance with the behavior of the experimental surface tensions and the predominant presence of the fluorinated IL observed by the ARXPS measurements [30].
To test the predictive character of Eq. 6 over a wider composition range, the surface tensions and surface compositions for binary mixtures of the non-fluorinated ILs [C 2 C 1 Im][Tf 2 N] and [C 12 C 1 Im][Tf 2 N] with the common [Tf 2 N] − (bis(trifluoromethylsulfonyl)imide) anion measured by Smoll et al. [12] at about 296 K were used. For this, the reported values for the relative surface exposure of the surface-active dodecyl chains in [C 12 C 1 Im] + , which is larger compared to the corresponding shorter ethyl chains in [C 2 C 1 Im] + , correspond to the surface mole fractions required in Eq. 6. For the 13 different mixture compositions ranging from pure [C 2 C 1 Im][Tf 2 N] to pure [C 12 C 1 Im][Tf 2 N], agreement between the experimental surface tensions [12] and the σ values calculated by Eq. 6 with an average absolute deviation of 0.8 % was found. All this indicates that the use of the surface mole fraction associated with the characteristic surface-active groups of each IL seems to be a reasonable parameter to describe the surface tension in binary mixtures of ILs containing a common anion via Eq. 6. Since this molar surface concentration is directly accessible by, e.g., ARXPS, it may be of advantage in comparison to commonly used surface volume or area concentrations of IL molecules (see in Refs. [12,70]) which can hardly be determined.

Conclusions
The present study contributes to an improved understanding of the surface tension and viscosity of binary mixtures of fluorinated and non-fluorinated ILs. Our work was not only motivated by the lack of detailed information about the thermophysical properties of such mixtures [25], but also by the question how the macroscopic properties viscosity and surface tension are influenced by the distinct structural organization in the bulk and at the surface of ILs containing fluorinated and nonfluorinated alkyl chains [17,28,30] 6 ]. By using the detailed microscopic information accessible by ARXPS, a new, phenomenological approach for the surface tension of the studied IL mixtures was suggested. Based on a linear mixing rule considering the surface tensions of the pure ILs and the surface mole fractions of the corresponding most surface-active group of each IL, i.e., the fluorinated and non-fluorinated chain, good agreement between the correlated and experimental surface tensions was found. Despite its reasonable representation for the present binary IL mixture and for a further one consisting of two non-fluorinated ILs [12], the relatively simple approach suggested in this work, Eq. 6, requires further validation by application to other IL mixtures over broad ranges of composition and temperature, which is intended in our future activities. Furthermore, the present contribution may stimulate further research, where the usage of IL mixtures including, e.g., fluorinated ILs serves as a promising way to tailor the surface and bulk properties for specific applications related to surface and colloid chemistry.