Experimental Data and Modeling of Solid–Liquid Equilibria of Binary Systems Containing Dibenzofuran and Long Chain n-alkanes

Solid–liquid equilibria (SLE) of binary mixtures of several n-alkanes (n-octadecane, n-eicosane, n-tetracosane, n-pentacosane, n-triacontane) and dibenzofuran covering the whole composition range were measured by differential scanning calorimetry (DSC) in a temperature range of 301–356 K. The dibenzofuran and the n-alkanes are completely miscible in the liquid state but non-miscible in the solid state, and in the context of this work, they seem to exhibit eutectic behavior. A linear trend is obtained for the eutectic temperature and eutectic composition versus the number of carbon atoms of n-alkane. The experimental data were compared to predictions made by using the ideal solution model, the DISQUAC model and several versions of the UNIFAC model, including the classical UNIFAC, the modified versions of Lyngby and Dortmund, and the recently proposed modified UNIFAC (NIST) model, to account for non-ideality in the liquid phase. Moreover, the experimental data were also compared with COSMO-RS predictions that utilize quantum chemical calculations and statistical thermodynamics to interpret the behavior of molecular interactions in the binary mixture. A good agreement was obtained between the predicted and the experimental temperatures.


Introduction
Solid-liquid equilibria (SLE) are important in chemical and petrochemical development, especially in avoiding problems with flow safety due to solid precipitation during extraction, transportation, and storage of complicated mixtures in liquid reservoirs and heavy oils [1][2][3][4][5][6].Polycyclic aromatic hydrocarbons (PAHs) and their derivatives containing oxygen, nitrogen, or sulfur atoms are present in crude oil, heavy oil products from coal, lubricants, or coal tar, and during the chemical conversion of crude oil or the incomplete combustion of hydrocarbon fuels.These chemicals are mutagenic and carcinogenic environmental Extended author information available on the last page of the article pollutants.The assessment of toxicity of PAHs and their derivatives necessitates the consideration of their solubility in various environmental media such as the atmosphere, soil, and water phases.PAHs are a class of organic compounds that exhibit volatility and have the ability to dissolve in atmospheric moisture.These compounds can also form aerosols and vapors, which can be inhaled or deposited on various surfaces.Furthermore, PAHs have the ability to adhere to the surfaces of particles or become physically integrated into particle structures, such as dust, or sediment.The solubility of PAHs in various environmental media plays a crucial role in determining their transport, fate, and bioavailability.This factor can have a considerable impact on the potential of PAHs to cause harm to both human health and the environment [7].
Systems containing polycyclic aromatic and long-chain normal alkanes are considered precursors to the complex systems found in fluid reservoirs.The representation of the thermodynamic properties of these fluids requires the consideration of the binary interactions between the different molecules involved.However, previous studies have been performed on binary mixtures related to the size of these molecules, either from two light components or from an asymmetric combination between heavy and light molecules.Although binary mixtures consisting of straight long-chain alkanes and PAH compounds have been extensively studied in the literature, studies on their thermodynamic properties have been limited.Djordjevic [8] measured the solubilities of five PAHs (naphthalene, acenaphthene, anthracene, phenanthrene, and pyrene) in n-octadecane.He concluded that these mixtures are slightly non-ideal, probably due to the compensation between enthalpic and entropic factors.
In addition, the excess free energy curves can be highly asymmetric, resulting in extremely high values for the limiting activity coefficients.For this reason, it is crucial from a practical point of view to have a thorough knowledge of the thermodynamic properties of these systems.Aoulmi et al. [9] investigated the thermodynamics of mixtures consisting of long-chain alkanes and PAHs.The authors investigated the mixing behavior of these compounds in terms of excess entropy of mixing, phase behavior, and heat of mixing while Mahmoud et al. [10] examined the SLE and excess enthalpy of the mixtures formed from long-chain n-alkanes and polycyclic heteropolyaromatics.Hafsaoui and Mahmoud [11] measured the SLE for naphthalene + n-tetracosane and dibenzofurane + n-tetracosane mixtures.A eutectic behavior was observed in these systems.The activity coefficients were obtained by readjusting the interaction parameters of the UNIFAC model.
PAHs are also significant to the refinery sector during hydrocracking because their deposition reduces heat exchange in these units.Aromatic-aromatic interactions stabilize the tertiary structure of proteins, for example and a deeper understanding requires the study of fundamental heteroatomic compounds such as indole, dibenzofuran, and dibenzothiophene [12].The study of mixtures containing a PAH and a long-chain alkane is essential because they play an important role in the flocculation and deposition of asphaltenes, a very critical issue during the exploitation, transport, and storage of crude oil.The study of these solutions is particularly useful as a starting point for the examination of mixtures containing heteroaromatic molecules, which are abundant in coal tar and coal liquids and are precursors to the complex systems found in fluid reservoirs.However, there are few studies on such binary mixtures in the literature [10,11,[13][14][15][16][17][18][19].
In this work, dibenzofuran was chosen as an example of a polyaromatic compound with a heteroatom and long-chain n-alkanes with more than 18 carbon atoms from n-octadecane to n-triacontane.Understanding the phase behavior of binary systems with dibenzofuran and long-chain n-alkanes is important for designing and optimizing separation techniques, which are needed to remove these contaminants from different industrial products and waste streams.Therefore, the aim of this work is to determine the SLE temperatures by means of differential scanning calorimetry (DSC) for dibenzofuran + n-octadecane, dibenzofuran + n-eicosane, dibenzofuran + n-tetracosane, dibenzofuran + n-pentacosane, and dibenzofuran + n-triacontane systems.In addition, the experimental results were correlated and compared with predicted values obtained with the ideal solution, UNIFAC [20], and several modifications of the UNIFAC [21][22][23][24][25], DISQUAC [26], and COSMO-RS [27] models to account for the non-ideality of the liquid phase.

Experimental Procedure
The dibenzofuran and the n-alkanes are completely miscible in the liquid state.Experimental methodology was followed as described in the literature [10,11,13]; for each binary system, the mixture was heated under constant agitation until complete melting, the sample was heated very slowly inside a glass cell, and the cell was immersed rapidly in a liquid nitrogen bath to obtain a homogeneous solid mixture.A small amount of solid (5 to 10 mg) was sealed in the aluminum pan of the DSC 7 (Perkin-Elmer).
The instrument was calibrated using Indium (purity = 99.99%,T m = 429.78K, ΔH f = 28.45J•g −1 ).The solubility measurements were conducted under a constant nitrogen flow (20 mL•min −1 ).The heating rate was fixed at 0.8 °C•min −1 .Data acquisition (enthalpy of fusion of the pure compounds and onsets corresponding to the solid-liquid equilibria temperatures) and processing was made with Perkin-Elmer's Pyris software.The measurement uncertainties are estimated to be ± 0.1 K for the temperatures and ± 0.35 kJ•mol −1 for the heat of fusion and solid-solid transition.i.e., the calculated reproducibility of the measurement without disturbance to the system.The calculated reproducibility (repeatability) of the measurement considering external factors such as sample mass, composition, and the fact of removing and replacing the sample in between each run was ± 0.4 K.The uncertainty of the mole fraction is 0.0005.The binary systems (dibenzofuran + n-eicosane), (dibenzofuran + n-tetracosane), and (dibenzofuran + n-pentacosane) have already been studied in the literature [10,11].Figures S6 to S8 (Supplementary material) show our values in comparison with those in the literature.The liquidus are very comparable, although our results are slightly higher.This is consistent with the higher melting temperatures of pure substances used during our experiments.It appears that the mean relative deviations on the liquidus temperatures are small.Values of the eutectic heat and of the heat of melting needed to determine the final composition of the eutectic points on the basis of the Tammann plots [14][15][16][17] (Figs.S1-S5, Supplementary material) are listed in Table S2, S3, S4, S5, and S6.

Determination of Solid-Liquid Phase Diagrams (SLPD)
For the five systems studied in this work, the components are not miscible in the solid phase and exhibit a simple eutectic diagram as the first peak appears at a constant temperature.The equation for the solid-liquid equilibrium is given by where ΔH m,i , T m,i , ΔC p,i , ΔH tr,i , and T tr,i are respectively the molar enthalpy of melting, the melting temperature, the molar heat capacity change (assumed to be independent of T) during the melting process, the enthalpy change corresponding to the transition, and the transition temperature of the pure component i .To a fair approximation, the term ∆C p,i , in Eq. 1 may be neglected [28].The last term corresponds to the enthalpic contribution of the first-order transition observed in the solid n-tetracosane, n-pentacosane, and n-triacontane.
According to Eq. 1, SLE data can be used to calculate activity coefficients, γ i , of the i-component in the mixture from experimental melting temperatures and pure component properties.Otherwise, when activity coefficients are calculated from a predictive model for the liquid phase, Eq. 1 can be used to predict the melting temperature for a given mixture.

UNIFAC Model
According to the UNIFAC model [20], Modified UNIFAC (Lyngby) by Larsen et al. [23,24], Modified UNIFAC (Dortmund) by Gmehling et al. [21][22][23], and Modified UNIFAC (NIST) by Kang et al. [25], activity coefficient γ i of component i in the liquid phase is given by For components i, r i, and q i are the volume and surface structural parameters, respectively, w i and w ′ i are volume fractions, and θ i is the surface fraction. (1) , where v (i) k is the number of groups k in the molecule i, Γ k is the group activity coefficient in the mixture, Γ (i)  k is the group activity coefficient of group k in pure i, z is the coordination number equal to 10, Q k and Θ k are the surface structural parameter and the surface fraction of group k, Ψ nm contains the interaction parameters: In the original UNIFAC version, α = 1 and the UNIFAC group interaction parameters between groups n and m, a nm , are assumed to be temperature independent.

Modified UNIFAC Model (Lyngby Version)
In Modified UNIFAC (Lyngby) by Larsen et al. [23,24], combinatorial contribution given by Eq. 3 is replaced by where α = 2/3 in w ′ i and Ψ nm is given by where T 0 is a reference temperature set equal to 298 K and the interaction between groups n and m is given by the three parameters a nm,1 , a nm,2 , and a nm,3 to improve the temperature effect.

Modified UNIFAC Model (Dortmund Version)
In the Modified UNIFAC (Dortmund) by Gmehling et al. version [21][22][23], α = 3/4 in w ′ i and Ψ nm is given by where interaction between groups n and m is given by parameters a nm,1 , a nm,2 , and a nm,3 .

New Modified UNIFAC Model (NIST Version)
In the New Modified UNIFAC by Kang et al. version [25], α = 3/4 in w ′ i and Ψ nm is identical to Eq 7″.This version included the specific group "furan" for describing cyclic ethers as dibenzofuran.
Structural parameters, R and Q, and interaction parameters for the several UNIFAC model versions for the ether group are listed in Table 1.

DISQUAC Model
Some important features of the model are the following.(i) DISQUAC is a group contribution model based on the rigid lattice theory developed by Guggenheim.[29] (ii) The total molecular volumes, r i , surfaces, q i , and the molecular surface fractions, α i , of the compounds present in the mixture, are calculated additively based on the group volumes R G and surfaces Q G recommended by Bondi [30] (see below).(iii) The partition function is factorized into two terms, in such a way that the excess functions are calculated as the sum of two contributions: a dispersive (DIS) term, which represents the contribution from the dispersive forces; and a quasichemical (QUAC) term which arises from the anisotropy of the field forces created by the solution molecules.
In the case of the Gibbs energy, G E m a combinatorial term G E,COMB m represented by the Flory-Huggins equation [31,32] must be considered.Thus, (iv) The interaction parameters depend on the molecular structure; (v) the value z = 4 for the coordination number is used for all the polar contacts.This represents one of the more important shortcomings of the model and is partially removed via the hypothesis of considering structure-dependent interaction parameters.(vi) The equations used to calculate the DIS and QUAC contributions to G E m and H E m in the framework of DISQUAC are given elsewhere [33].The temperature dependence of the interaction parameters, g where s ≠ t are two contact surfaces present in the mixture and l = 1 (Gibbs energy); l = 2 (enthalpy), l = 3 (heat capacity).The needed equations are [33] (vii) Dibenzofuran + n-alkane mixtures are built by three surfaces (groups): type a, aliphatic (CH 3 , CH 2 ), in n-alkanes; type b, aromatic in dibenzofuran, and type e, oxygen (-O-in dibenzofuran).The geometrical parameters of the groups referred to in this work are available in the literature [14].For dibenzofuran, r = 5.2766, q = 3.4140, b = 0.9399, e = 0.0606.(viii) These surfaces generate three contacts: (a,b); (a,e); and (b,e).The (a,b) contact is characterized by purely dispersive interaction parameters equal to those of naphthalene + n-alkane mixtures [34].Interactions parameters for the (a,e) and (b,e) contacts are represented by DIS and QUAC interchange coefficients.The fitting procedure has been explained in detail elsewhere [33].We remark that, due to the necessary experimental data, the QUAC interaction parameters used for the (a,e) and (b,e) contacts are those reported for aromatic ether + alkane systems [35].Table 2 shows our values of the interaction parameters for the mentioned contacts.The geometrical parameters for the groups referred to in this work are listed in Table S1 (Supplementary Material).( 8) .

COSMO-RS Model
The Conductor-like Screening Model for Real Solvents (COSMO-RS) was used to predict and compare the solid-liquid phase equilibria of dibenzofuran + n-alkanes [27,36].
The geometry optimization of the organic compound (dibenzofuran) and the n-alkanes (n-octadecane, n-eicosane, n-tetracosane, n-pentacosane, and n-triacosane) was performed using the Turbomole program package.The chemical structure of the target molecule was firstly drawn in three dimensional.Afterward, geometry optimization was performed at the Hartree-Fock level and 6-31G* basis set.In this work, only the lowest energy conformer was considered for the COSMO-RS calculation.During the optimization, the positions of atoms in the molecule were adjusted to find the lowest energy conformation of the molecule.This involved changing the three-dimensional coordinates of the atoms by varying the structural factors such as bond lengths, bond angles, bond rotations, and dihedral (torsional) angles.The generation of the.cosmo file was conducted through a single-point calculation using DFT with Becke-Perdew and the Triple-ζ Zeta Valence Potential (TZVP) basis set.Finally, the.cosmo files were exported to the COSMOthermX program of version 19.0.5, with parameterization BP_TZVP_19.ctd.COSMO-RS is a computational model which combines two methods, i.e., (i) quantum chemical considerations in the form of a conductor-like-screening tool and (ii) statistical thermodynamics.This enables the determination and prediction of thermodynamic properties without requiring experimental data.In this work, COSMO-RS was used to compute a range of binary concentrations for dibenzofuran + n-alkanes at which the solidification occurs.The comparison between COSMO-RS prediction data can be compared visually in phase equilibria involved in this study.
The following equation calculates the chemical potential of the surface segments: This set of COSMO-RS equations gives the chemical potential of all components in a mixture, which enables the estimation of several thermodynamic properties, such as activity coefficients, selectivities, distribution ratios, and phase equilibria.The detailed derivations have been described by the developers elsewhere [37][38][39], where the main resulting equations applied in COSMO-RS and their constants are summarized in Table 3.
In this work, when applying the COSMO-RS model, the second and the third terms in Eq. 1 were neglected (both the difference between the molar heat capacity of the solid and the subcooled liquid form and the enthalpic contribution of the first-order transition, respectively) yielding to the simplified equation: . So, the only input data for the COSMO-RS model are the enthalpy of fusion (∆H m ) and the melting point temperature (T m ) for each pure species manually assigned according to the data reported in Table 4.

Results and Discussion
The DSC heating thermograms corresponding to binary mixtures show two peaks.The onset temperature of the first peak (the one observed at the lowest temperature) is independent of the concentration and corresponds to the eutectic transition, while the second one is concentration dependent, and its onset temperature corresponds to the liquidus temperature.At the eutectic composition, the thermogram only presents a peak whose shape is similar to that of the pure compounds.
Phase change properties of the pure compounds are reported in Table 4, together with the values from the literature.These properties are globally in fair agreement with literature data [10,11,13,[40][41][42][43][44].The experimental data obtained for the SLE determination of dibenzofuran with n-alkane systems are listed in Table 5.
Table 5 lists experimental data obtained for the SLE determination of dibenzofuran with n-alkane systems.n-alkane range from n-octadecane up to n-triacontane and covered temperature range was 301-356 K. Eutectic points are determined, and they are (14)   , corresponding to the eutectic points depend on the chain length of n-alkane and increase with the increase of the number of carbon atoms.
Figure 1 represents the eutectic coordinates, composition, and temperature versus the number of carbon atoms for n-alkane, and a straight forward linear trend can be found for both magnitudes.
The objective of this study was a better understanding of the thermodynamic properties of mixtures formed by dibenzofuran and one long-chain n-alkane hydrocarbon.Calculated values were performed with ideal solution, classical UNIFAC [18], and modified UNIFAC versions of Lyngby [23,24], Dortmund [21][22][23], and NIST [25] models, with the parameters from Table 1.
Predicted melting temperatures were compared with those obtained experimentally by using DSC in terms of the mean percent relative deviation, σ T (%), as defined in Eq. 15: where n is the number of experimental data points, and subscripts exp and cal stand for experimental and calculated values.Values obtained for each studied system and the whole experimental data collection are reported in Table 6.  Figure 2 plots melting temperature versus mole fraction, and Fig. 3 plots the dibenzofuran activity coefficient for the system Dibenzofuran (1) + n-C18 (2); experimental and calculated values are plotted together.As expected, the ideal solution results in a clear underestimation of melting temperatures and composition for the eutectic point higher than the experimental value.Most UNIFAC model versions change such descriptions   and overestimate the deviations from ideal behavior, as shown in Fig. 3. Thus, they yield melting temperatures higher than experimental ones and shift the calculated eutectic composition to lower values.Similar plots and analyses could be obtained for the other binary systems studied.Calculated deviations from Table 6 show that the overall value for ideal solution, 2.58%, is similar to those obtained for UNIFAC model versions except for the Lyngby model, which yields the best results with a global deviation of around 1%, as it could also be seen in Fig. 2.
A similar comparison can be carried out for the activity coefficient values.Equation 16 gives the expression for the average standard deviation (%) between activity coefficients for each component calculated by the thermodynamic model and those obtained experimentally and listed in Table 5: Table 7 lists obtained values for the several models considered.As can be checked, the deviations in the activity coefficients are usually very high, mainly for component 1, despite the reasonably good deviations obtained for the melting temperatures.The obtained deviations showed that the application of this model with original interaction parameters could not reproduce a good behavior for the studied systems.
It was possible to represent better all systems studied by a readjustment of some of the UNIFAC interaction parameter models.Binary interaction parameter values between the aliphatic group with both aromatic and ether groups were assumed to have the main effect in these systems, according to S. L. Hafsaoui and R. Mahmoud [11].Alternative parameter values for both versions of the modified UNIFAC model were taken from [11], and values for the UNIFAC model were obtained in the present work.Parameter optimization was carried out by minimizing the sum of the two objective functions defined as follows: F 01 and F 02 indicate the objective function representing the part rich in dibenzofuran and that rich in n-alkane.Alternative values for interaction parameters are listed in Table 8.
Predicted melting temperatures and predicted activity coefficient values with these alternative parameters were compared with those obtained experimentally in terms of the mean percent relative deviation, σ(%), as defined in Eqs. 15 and 16, respectively.Accepted values are listed in Tables 9 and 10.The modified UNIFAC (Lyngby) does not improve the description and leads to an overall deviation of σ T = 1.21%, similar to the original model versions.However, the modified UNIFAC (Dortmund) version, even when given the least accurate description, σ T = 1.56%, represents a significant improvement from the original parameters.A better description is obtained with the original UNIFAC and the modified UNIFAC (Lyngby) with the furan group, with deviations ranging from 0.6 to 0.7%.Similar conclusions could be obtained when the comparison takes into account the deviations in activity coefficients (Table 10).Compared with (16 .

Table 7
Mean percent deviations in activity coefficients resulting from calculations using three UNIFAC versions with original parameters from  results provided by the other models, DISQUAC results are slightly better than those obtained from the UNIFAC models using original interaction parameters.In contrast, calculations from Dortmund UNIFAC with new interaction parameters slightly improve DISQUAC results.
Figures 4 and 5 plot the experimental and predicted values for the eutectic melting temperature and eutectic composition.The linear trend previously described in Fig. 1 is correctly described; however, deviations in eutectic composition are more significant than deviations in eutectic melting temperature, as expected from the low percent deviations obtained for melting temperatures.
The interactions between dibenzofuran and n-alkanes molecules can explain the solid-liquid phase equilibria based on the σ-profile and σ-potential.Its developers have demonstrated the validity of this approach in detail and reported in several recent studies involving molecular interactions analysis [37,38].In the σ-profile, when the screening charge density exceeds ± 0.0084 eÅ −2 , the molecule is considered adequately polar to induce hydrogen bonding.A higher absolute value of σ leads to a compound being a stronger hydrogen bond donor or hydrogen bond acceptor.For the horizontal axis, increasing negative and positive values for the hydrogen bonding threshold (± 0.0084 eÅ −2 )

Table 10
Mean percent deviations resulting from calculations using three UNIFAC versions with modified parameters (Table 8) and from DISQUAC application using the parameters listed in indicates the region of a molecule where interactions between hydrogen bond donors and hydrogen bond acceptors occur, respectively.On the other hand, the σ-potential indicates the affinity of a component in a mixture towards another.In the σ-potential plot, a higher negative value of µ indicates an increasing interaction between molecules, whereas a more elevated positive value indicates an increase in repulsive behavior.For the horizontal axis, increasing negative and positive values for the hydrogen bonding threshold (± 0.84 e/nm 2 ) means the region of a molecule where interactions between hydrogen bond donors and hydrogen bond acceptors occur, respectively.The σ-profile and σ-potential of all species involved in this study are depicted in Figs. 6 and 7.
The σ-profile of dibenzofuran has a curve that resolves mostly in the non-polar region.The three significant peaks from the left observed in the dibenzofuran curve correspond to the C-H bonds (− 0.006 eÅ −2 ), C atoms with π-electron of the benzene ring (− 0.001 eÅ −2 ), and O atom (0.004 eÅ −2 ).The molecular drawing also represents these peaks by the Fig. 4 Experimental eutectic composition of binary dibenzofuran + n-alkane and studied models versus the alkane number of carbon atoms Fig. 5 Experimental eutectic melting temperatures of binary mixtures dibenzofuran + n-alkane and studied models versus the alkane number of carbon atoms colors green-cyan, yellow, and orange, respectively.Only a small portion of the dibenzofuran curve falls on the polar regions.On the other hand, the sigma profiles of n-alkanes are totally in the non-polar region, with the area under the peak (A) increased proportionally to the carbon number, i.e., A n-C18 < A n-C20 < A n-C24 < A n-C25 < A n-C30 .This can be expected since the addition of carbon number is associated with the increased non-polar charge density of the -CH 2 group.
The polarity and interaction behaviors of dibenzofuran and n-alkanes are consistent with the σ-potential depicted in Fig. 7, where the attractive behavior between occurs only in the non-polar region.Although all species showed a similar shape of the curve, dibenzofuran (Fig. 7a) gave a much higher magnitude of repulsive (positive µ) and attractive (negative µ) behavior than those in n-alkanes (Fig. 7b).Dibenzofuran has a tiny portion of attractive forces in the polar region.Since dibenzofuran has a relatively higher magnitude of σ-potential than n-alkanes, it can be suggested that the increasing molar binary composition of n-alkanes (x 2 ) will have a low impact on reducing the eutectic point.This is in agreement with the curves observed in Fig. 2, where the melting point (T m ) of the binary mixture decreased only slightly at 0.5 ≤ x 1 < 1.0 and reduced rapidly when x 1 < 0.5.
COSMO-RS was also able to predict solid-liquid phase equilibria data comparable with, and in some plots, more accurate than other models.For dibenzofuran and n-C18 binary mixtures, it can be seen in Fig. 2 that the SLE data points predicted by COSMO-RS gave higher accuracy than those calculated by the UNIFAC and Mod UNIFAC (Dortmund), similar to those of Mod UNIFAC (NIST), and lower than those of Mod UNIFAC (Lyngby).Similarly, in Figs. 3, 4, and 5, COSMO-RS proved its ability to comparably predict the activity coefficients, eutectic compositions, and the eutectic melting temperature of binary mixtures of dibenzofuran and n-alkanes.These findings demonstrated the advantage of using COSMO-RS as a predictive tool for calculating thermodynamic properties of mixtures based on molecular information, which requires minimal experimental input.While being an 'a priori' method, its inclusion can improve the accuracy of the results by providing additional information about the system.It can also be concluded that COSMO-RS is a valuable tool for predicting the properties of complex mixtures where experimental data may be limited or unavailable, thus offering significant advantages for researchers and practitioners in various fields.

Conclusions
DSC technique is accurate enough to determine the solid-liquid equilibrium data for systems dibenzofuran + n-alkanes.Agreement with previous experimental data for pure components was achieved.The method allows the study of the composition effect and the accurate determination of the coordinates for the eutectic point.A linear trend is obtained for the eutectic melting temperature and the eutectic composition versus the alkane number of carbon atoms.
Experimental data show clear non-ideal behavior, verified by predicting the SLE using the ideal solution model.Non-idealities were computed using the group contribution method UNIFAC.Original UNIFAC and several modifications, Mod UNIFAC (Lyngby), Mod UNIFAC (Dortmund), and Mod UNIFAC (NIST), were considered.In general, nonidealities are overestimated, and predicted melting temperatures are higher than experimental values.The best description is obtained with Mod UNIFAC (Lyngby) method, which leads to the lowest deviations in melting temperature and activity coefficients.A better description is obtained when specific parameters are used.Previously available parameters obtained from data for similar systems lead to lower deviations for melting temperatures and activity coefficients.Data reported here were used to obtain specific parameters for the classical UNIFAC model, and the model's accuracy was improved.
The most recent version of the UNIFAC model, Mod UNIFAC (NIST), includes a specific group for cyclic ethers to define furan compounds.However, such a model version does not improve the results obtained from previous versions.Combining the Mod UNI-FAC (Lyngby) with the new group definition (furan group) requires the determination of four binary parameters.Still, the final version leads to an accuracy similar to classical UNI-FAC with specific parameters.DISQUAC results are slightly poorer than those obtained from Dortmund UNIFAC with modified interaction parameters.
COSMO-RS was also used in this work to generate the SLE data predictions based on quantum chemical calculations and statistical thermodynamics.The σ-profile and σ-potential developed in COSMO-RS were analyzed to view the behavior of molecular interactions in the binary mixture.This qualitative analysis supported the SLE data points obtained in all calculation models.Prediction from COSMO-RS are similar to UNIFAC (NIST) in spite of COSMO-RS not being a group contribution model.

Fig. 1
Fig. 1 Experimental eutectic composition (Filled circle) and eutectic melting temperature (Filled Square) of binary mixtures dibenzofuran + n-alkane versus the alkane number of carbon atoms

Table 1
Structural and interaction parameters for ether group involved in the dibenzofuran + n-alkane systems

Table 3
Equations and constants in COSMO-RS modelHydrogen Bond interaction energy E hb = eff c hb min 0, � + 2

Table 5 (
*Corresponding to the eutectic point

Table 6
Mean percent deviations in precipitation temperature resulting from calculations using UNIFAC versions, DISQUAC, and COSMO-RS model

Table 8
Alternative interaction parameters for ether group involved in the dibenzofuran + n-alkane systems

Table 9
Mean percent deviations in precipitation temperature resulting from calculations using three UNI-FAC versions with modified parameters (Table8) and the DISQUAC model